вался удельный (отнесенный к изменению координаты) расход топлива G7…G10, вызванный изменением соответствующей управляемой координаты q7 … q10 соответственно.
Анализ полученных экспериментальных зависимостей G7…G10 показал, что они имеют существенно нелинейный характер по всем переменным-предикторам (q8, q9, mГР). Выбросы и ошибки в выборке отсутствовали, т.к. использовалась имитационная математическая модель и вычисления с двойной точностью.
Таблица 4.2. Значения коэффициентов bi (i [1;27]) уравнений регрессии
удельного расхода горючего G7 при произвольном направлении изменения координаты q7 и удельных расходов горючего G8 … G9 при изменении координат q8, q9 на подъем груза
Обозн. |
Коэффициенты уравнения регрессии для |
||
коэф. |
G7 |
G8 |
G9, (bi∙10–3) |
b1 |
0,001114309572183 |
0,010520260360517 |
0,635585822053724 |
b2 |
0,000000126996818 |
0,000005303960654 |
0,000067903845021 |
b3 |
0,000000000006045 |
0,000000000213632 |
–0,000000015008719 |
b4 |
0,000012740405459 |
0,000418003360941 |
0,027103231553251 |
b5 |
0,000000021451692 |
0,000000748974903 |
–0,000014045239655 |
b6 |
0,000000000001759 |
–0,000000000024556 |
–0,000000002076003 |
b7 |
0,000001633276816 |
0,000011242631491 |
–0,001770730237815 |
b8 |
0,000000002253500 |
–0,000000011202115 |
0,000000959041760 |
b9 |
0,000000000001078 |
0,000000000001278 |
–0,000000000103330 |
b10 |
–0,000122674836024 |
0,004423327921248 |
0,717348087006645 |
b11 |
–0,000000067638361 |
–0,000001195893002 |
0,000413822100503 |
b12 |
–0,000000000011689 |
–0,000000000274757 |
0,000000034964904 |
b13 |
–0,000013710606624 |
–0,000040605796028 |
–0,077632454807223 |
b14 |
–0,000000002776654 |
–0,000000645423255 |
0,000047480462805 |
b15 |
–0,000000000004185 |
0,000000000033372 |
0,000000004017220 |
b16 |
–0,000000487381117 |
–0,000020463350013 |
0,005287565497309 |
b17 |
–0,000000002656128 |
0,000000017580833 |
–0,000002952146167 |
b18 |
–0,000000000001937 |
0,000000000011043 |
0,000000000080907 |
b19 |
0,000012341644141 |
–0,006738098722693 |
–0,232493435483927 |
b20 |
–0,000000019596043 |
–0,000001720905606 |
–0,000102746683669 |
b21 |
0,000000000005445 |
0,000000000078589 |
–0,000000010053091 |
b22 |
0,000002314154312 |
–0,000187033946874 |
0,050500245691884 |
b23 |
–0,000000009561386 |
0,000000109880307 |
–0,000029322463849 |
b24 |
0,000000000002203 |
–0,000000000007180 |
–0,000000001419451 |
b25 |
–0,000000552283474 |
0,000009459885495 |
–0,003656292965043 |
b26 |
0,000000000807746 |
–0,000000004925676 |
0,000002028652411 |
b27 |
0,000000000000849 |
–0,000000000011943 |
–0,000000000044880 |
|
|
230 |
|
Таблица 4.3. Значения коэффициентов bi (i [1;27]) уравнений регрессии
удельных расходов горючего G8 … G9 при изменении координат q8, q9 на опускание груза
Обозн. |
Коэффициенты уравнения регрессии для |
|
коэф. |
G8 |
G9 |
b1 |
0,006150713211537 |
–0,000279520475304 |
b2 |
0,000007472775903 |
0,000000278130539 |
b3 |
–0,000000000463230 |
0,000000000049579 |
b4 |
0,001354704293717 |
0,000273872832935 |
b5 |
0,000002232490517 |
–0,000000005355766 |
b6 |
–0,000000000094972 |
0,000000000004368 |
b7 |
–0,000109038612316 |
–0,000012699802403 |
b8 |
–0,000000075334492 |
–0,000000000316236 |
b9 |
–0,000000000040923 |
–0,000000000002016 |
b10 |
–0,005870926219376 |
0,003174241418732 |
b11 |
–0,000004766275782 |
–0,000000153815563 |
b12 |
0,000000000517270 |
–0,000000000155266 |
b13 |
0,000245681370756 |
–0,000528490082249 |
b14 |
–0,000003632401732 |
–0,000000120689926 |
b15 |
0,000000000211808 |
0,000000000016244 |
b16 |
0,000190060145800 |
0,000020705601745 |
b17 |
0,000000226962927 |
0,000000014804293 |
b18 |
0,000000000032111 |
0,000000000001236 |
b19 |
0,006652575474157 |
–0,001877838166523 |
b20 |
–0,000000016848817 |
0,000000130321072 |
b21 |
–0,000000000142067 |
0,000000000094174 |
b22 |
–0,001212305714707 |
0,000445626257559 |
b23 |
0,000001412371695 |
0,000000084575233 |
b24 |
–0,000000000107578 |
–0,000000000017465 |
b25 |
–0,000064796689961 |
–0,000021680239735 |
b26 |
–0,000000119078526 |
–0,000000010864560 |
b27 |
–0,000000000000964 |
0,000000000000369 |
Предварительно были получены данные о неадекватности линейной модели множественной регрессии удельных расходов топлива G7 … G9 от предикторов, а также неадекватности нелинейных по переменным моделей, линеаризованных: полиномами сумм одночленов из одной независимой переменной до 2…4-й степеней включительно, логарифмическими и экспоненциальными функциями (погрешность δmax для всех перечисленных моделей составила свыше 40 %). Снижение размерности модели (уменьшение числа предикторов либо их преобразование в меньшее число переменных) без существенного снижения точности регрессии также оказалось невозможным.
231
Было принято решение об использовании нелинейной множественной регрессии формулой вида (4.144). Выражение представляет собой симметричный многочлен от трех переменных-предикторов в степенях [0; 1; 2] в 27-ми всевозможных сочетаниях. Степень многочлена (4.144) по совокупности всех переменных – 6 [171]. Использовалась реализация алгоритма Левенберга-Марквардта в программном продукте STATISTICA 8 [32, 194, 245].
Значения коэффициентов bi (i [1;27]) уравнений регрессии удельных расходов горючего G7 … G9 при изменении управляемых координат ГПК q8, q9, вызывающем подъем груза, приведены в табл. 4.2, а вызывающем опускание груза – в табл. 4.3.
Таблица 4.4. Значения показателей качества уравнений регрессии удельного
расхода горючего G7 … G10 при изменении управляемых координат ГПК q8 – q10, направленном на подъем груза
Показатель |
|
Координата |
|
||||
G7 |
G8 |
G9 |
G10 |
||||
|
|
|
|||||
Коэффициент детерминации R2 |
0,9968 |
0,9998 |
0,9999 |
0,9994 |
|||
Скорректированный коэффи- |
0,9959 |
0,9998 |
0,9999 |
0,9994 |
|||
циент детерминации |
|
2 |
|||||
R |
|||||||
Критерий Фишера F |
30305,6 |
123139,4 |
7351,1 |
413351,1 |
|||
Сумма квадратов остатков RSS |
3,0374∙10–8 |
1,4812∙10–6 |
4,0895∙10–8 |
1,1913∙10–6 |
|||
Стандартная ошибка уравнения |
1,7605∙10–5 |
1,2294∙10–4 |
2,0428∙10–5 |
9,8818∙10–5 |
|||
регрессии SEE |
|
|
|
|
|||
Максимальная относительная |
2,82 |
2,1 |
2,89 |
3,32 |
|||
погрешность аппроксимации |
|||||||
δmax, % |
|
|
|
|
|||
Таблица 4.5. Значения показателей качества уравнений регрессии удельного
расхода горючего G8 … G10 при изменении управляемых координат ГПК q8 – q10, направленном на опускание груза
Показатель |
|
|
Координата |
|
||
|
G8 |
G9 |
G10 |
|||
|
|
|
|
|||
Коэффициент детерминации R2 |
|
0,9976 |
0,9982 |
0,9914 |
||
Скорректированный коэффициент |
|
0,9971 |
0,9969 |
0,9914 |
||
детерминации |
|
2 |
|
|||
R |
||||||
Критерий Фишера F |
|
227949,8 |
16243,4 |
528362,3 |
||
Сумма квадратов остатков RSS |
|
2,4932∙10–6 |
6,3847∙10–7 |
3,9113∙10–7 |
||
Стандартная ошибка уравнения |
|
3,5274∙10–5 |
4,5627∙10–5 |
6,3428∙10–5 |
||
регрессии SEE |
|
|
|
|
||
Максимальная относительная по- |
|
3,3548 |
2,93 |
2,14 |
||
грешность аппроксимации δmax, % |
|
|||||
|
|
|
|
|||
|
|
|
232 |
|
|
|
Уравнение регрессии удельного расхода горючего G10 при подъеме груза (л/УЛЕ), отнесенного к изменению управляемой координаты ГПК q10, имеет более простой вид, т.к. зависит только от одного из трех учитываемых предикторов – от массы груза mГР:
G10=(0,0644832+0,00000453877∙mГР)2. |
(4.148) |
Удельный расход горючего G10 при опускании груза (л/УЛЕ), отнесенный к изменению управляемой координаты ГПК q10, аппроксимируется аналогичной зависимостью:
G10=(0,0454509+0,00000278131∙mГР)2. |
(4.149) |
Вид данных зависимостей – параболический.
Анализ показателей качества уравнений множественной нелинейной регрессии величин G7 … G9, аналогичных уравнению (4.144), и квадратичных уравнений вида (4.148) и (4.149) (табл. 4.4, 4.5) показал, что регрессия по уравнениям данного вида дает наилучшие результаты. Все коэффициенты уравнений регрессии величин G7 … G10, согласно t- статистике Стьюдента, значимы. Максимальная относительная погрешность аппроксимации δmax во всем рассматриваемом диапазоне изменения предикторов не превышает 3,4 % (см. табл. 4.4, 4.5, рис. 4.28).
|
G7, |
q9=10 УЛЕ |
|
|
Экспериментальные графики (∙∙∙∙∙) |
|
|
G7, л/рад |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
л/рад |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
mГР=750 кг |
||||||||||||
|
|
|
|
mГР=500 кг |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
mГР=250 кг |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q9=10 УЛЕ |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q8, рад |
q9, УЛЕ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
mГР=250 кг |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
Регрессионные графики (—) |
б) |
|||||||||
|
|
а) |
q8, рад |
||||||||||||||
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Рис. 4.28. Экспериментальные и регрессионные зависимости удельного расхода G7 (пример): а – поверхности, соответствующие mГР =250 кг; б – линии, соответствующие q9=10 УЛЕ и нескольким различным mГР
Полученные регрессионные уравнения позволяют в наглядной форме в виде графических зависимостей представить удельные расходы топлива при изменении управляемых координат как функции значимых с точки зрения энергозатрат параметров технологического процесса ГПК: q8, q9, mГР (рис. 4.29).
233