Файл: Алабин М.А. Корреляционно-регрессионный анализ статистических данных в двигателестроении.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.07.2024
Просмотров: 108
Скачиваний: 0
h + |
?3 Г23 + ?4 Г24 + ■••+ ?/п Гъ т ---- |
Г12 |
|
|||
?2Г32 + |
?4 Г34 + |
+ ? т Г3ш — Г13 |
(31") |
|||
|
|
|
|
|
|
|
(V m 2 + |
РзГ тЗ + |
?4Г ш4 + • • ■ + ?/» — |
Г1 т |
|
||
П ри ЭТОМ / m m — 1; |
Г(т—1ут— I т(т—1)- |
|
|
^ |
||
Перестановка |
членов позволяет оценить |
смысл параметров |
||||
|
. 2 ----'1 2 |
^3Г23 |
lJ4r 24 ................. |
|
|
|
|
®3= |
Г13 |
f^2r 32 |
i^4r 34 ................. |
|
|
|
Pm |
Пт |
^2^”m2 |
iJ3^m3 ........ |
|
|
Из этого следует, что каждый из параметров (32, (Зз, •• Рш, ха рактеризующий зависимость результирующего параметра от со ответствующего составляющего параметра, равен коэффициенту парной корреляции результирующего параметра с одним из сос тавляющих параметров за вычетом влияния остальных состав ляющих параметров, действующих косвенно через влияние на этот параметр. Абсолютная величина параметра позволяет су дить об относительной силе влияния каждого составляющего па раметра; знак минус указывает при этом на обратную зависи мость.
Параметры уравнения корреляции в натуральном масштабе вычисляются по, формулам
“X, |
их , |
' 12.34. ,.т ~ =f>2- |
'13.24...т " |
°Х , |
|
Свободный член уравнения корреляции после определения зна
чений &12.34 ... т \ |
&I3.24 ... т , • • • • НЭХОДИТСЯ |
ПО формуле |
|
- * 1 - |
' ^12.34. |
'^13.24. . . т ^ З |
'lm .2 3 ...m - 1 * Я |
1.7. Способы решения системы нормальных уравнений |
|||
Коэффициенты системы |
нормальных уравнений могут быть |
||
вычислены: |
|
|
|
—алгебраическим решением;
—механизированным путем с использованием ЭВМ. Наиболее достоверные результаты, особенно при решении
системы из четырех и более уравнений по большому количест ву наблюдений, могут быть получены только на ЭВМ.
Алгебраическим путем целесообразно решать системы с чис лом уравнений не более четырех. Первый путь алгебраического решения заключается в том, что дополнительными преобразова ниями обеспечивается решение уравнения с одним неизвестным,
23-
а идя в обратном , порядке, находят значения всех остальных неизвестных. Второй путь алгебраического решения основан на использовании определителей.
Ниже приведены типовые порядки алгебраического решения следующей системы уравнений (записанной в форме отклонений от среднего значения каждой случайной величины):
^12.34 ^ |
Х\"Т~ ^13.24 ^ Х2Х3+ +4.23 2 |
= ^ |
-ДЛ*2> |
+2.34 2 |
х 2х зТ“ +3.24 V x - -j- bx423 ^ |
+А:4= ^ |
XjA'g, |
+2.34 |
А2Х4+ ^13,24 ^ АГ3Х4-(- ^Х4.23 |
^4 = |
-И-П. |
Ввиду сравнительной громоздкости решения указанной си стемы уравнений в общем виде, произведем ее решение приме нительно к следующим уравнениям:
3 ,3^ 2.34+ 8 ,8^3 24 + 15614_;гз= 106,2; |
(а) |
|
8,8^12.34 |
25Z,2^x3.24 80^14.23:= 4-90,5; |
(б) |
15612.з4 |
5064324+ 10064423= 150. |
(в) |
Для нахождения значений коэффициетов 612.з4; 613.24; 614.2з над уравнениями последней системы производятся действия в сле дующем порядке.
1. Все члены уравнения (а) делятся на постоянный коэффи циент при 612.34, а знаки заменяются на обратные
—612.34—2,666613.24—4,5456i4 2з= —32,182. |
(а/) |
2. Все члены уравнения (а) умножаются на абсолютное зна чение коэффициента при втором члене уравнения (а')
8,86 i2.34 + 23,4608613.24+ 39,99614.23= 283,129. |
(&") |
3. Из уравнения (а'') вычитается уравнение (б)
280,666x3.244"89,99644 2з=773,629. |
(г) |
4. Делим все члены уравнения (г) на коэффициент при пер вом члене
613.24 + 0,32066i4.23= 2,756. |
(г7) |
5. Умножая все члены уравнения (а) на абсолютное значе ние коэффициента при третьем члене уравнения (а'), т. е. на
4,545:
15612.34+ 39,9966is.24 + 68,1756,4 2з = 482,679. |
(а7//) |
6. Вычитаем из (6) уравнение (а"')
—906,з.24+ 31,8256,4.2з= —332,679. ■ |
(г") |
24
7. Умножаем все члены уравнения (г) на коэффициент при втором члене уравнения (г'), т. е. на 0,3206:
90й13.24 + 28,851 £„.2з = 248,025. |
(г'") |
8. Складываем уравнения (г") и (г'")
60,676&14.2з= —84,654.
Тогда
^14.23= — 1,395.
9.Подставляем Ьц.23= — 1392 в (г')
*>13.24 + 0,32(— 1,392) =2,757.
Отсюда
£13.24= 3,204. |
|
|
'10. Подставляем значения 613.24 и £14.23 |
в уравнение (а') |
|
Ьх2.34 + 2,666-3,204 + 4,545 (— 1,392) = 32,182. |
|
|
Тогда |
|
|
£12.34= 30,0. |
|
|
Величины £12.34; £13.24; £i4.23> вычисленные таким образом, |
мо |
|
гут быть проверены подстановкой их значений _в уравнение |
(в). |
|
Подставляя соответствующие значения |
Х 2; Х 2; Хц и коэффи |
циентов £12.34, £13.24, £14.23 в уравнение для определения свобод ного члена, получим
6о=—213.
Тогда искомое уравнение регрессии, представляющее собой зависимость тяги на максимальном режиме от расхода воздуха через двигатель, диаметра реактивного сопла, температуры газов за турбиной, будет иметь вид
£>= — 213 + ЗО0В+ 3,2047*- 1,39а .
Система уравнений (а), (б), (в) может быть .решена и дру гим способом.
Разделим все члены каждого уравнения на коэффициенты при одной и той же неизвестной величине. Осуществляя это по отношению к £12.34, будем иметь
£12.34+ 2,666£13.24 + 4,545£i4.23 = 32,182;
£12.34— 29,227£13.24— 5,682^14.23= — 55,738; £12.34—3,333£i3.24+ 6,666614.23= 9,999.
Вычитая из второго и третьего уравнений первое, получим си стему двух уравнений с двумя неизвестными
31,893£i3.24+ 10,227£i4.23— 87,920;
5,999£ i3.24— 2,121 £14.23= 22, 183.
25
Тем же путем преобразуем эту систему к виду
3,118£*13.24 |
= 8,596; |
2,8286i3.24— 614.23 = 10,458. |
|
Откуда 613.24= 3,204. |
|
Дальнейший порядок определения |
коэффициентов 613.24 и 612.34 |
не вызывает затруднений. |
|
При алгебраическом решении системы уравнений вторым путем неизвестное значение каждого коэффициента уравнения регрессии находится как частное от деления двух определителей. Составление и решение определителей осуществляется в следую щем порядке.
1. Вначале составляется определитель знаменателя, предст ляющий собой все члены левой части системы уравнений, на писанные в табличной форме:
2 Х22 ХйХ3..... ■2 Xn-x 'i
2 Л*2Л'з 2 Х3 ‘ ’' ~
Дэнам
^ Х3Хт Х3Хт" - 2 -^
2. Определитель числителя получается из определителя зн менателя заменой всех членов столбца при определяемом зна чении коэффициента уравнения регрессии на соответствующие значения членов правой части системы уравнений. Тогда, напри мер, определитель числителя для нахождения коэффициента 612.34... т будет иметь вид
2 |
х & |
2 |
* |
» * |
з " |
|
2 |
2 |
х л |
2 |
* |
з |
........ • 2 |
% |
|
Д чи сл |
|
|
|
|
|
|
|
2 |
х 1х т 2 |
|
Х 3Х Щ- |
..........У х 1 |
|||
|
|
Jmi 171 |
|||||
|
|
|
|
|
|
|
3. Производится решение определителей высших порядков тем постепенного перехода к определителям более низших по рядков, используя следующее свойство: определитель порядка т равен сумме произведений элементов какого-либо столбца (или строки) на их алгебраическое дополнение.
Алгебраическое дополнение некоторого элемента есть минор этого элемента, взятый со своим знаком, если сумма номеров строки и столбца, на пересечении которых расположен элемент,
26
есть число четное и с обратным знаком, если это число — нечет
ное. |
Минором некоторого |
элемента определителя |
порядка |
т |
называется определитель |
порядка т— 1, получаемый из исход |
|||
ного |
определителя путем |
вычеркивания строки и |
столбца, |
на |
пересечении которых расположен этот элемент.
Таким образом, например, из определителя 3-го порядка по
лучается три определителя второго порядка |
|
|
||||||
|
У X2 |
У з д |
V a v |
|
|
У 4 |
У А3А4 |
|
. Дэнам |
V Г \ - |
V Л-2 |
V V X |
|
|
|||
_ ^ Л2Л3 |
_ £ Л3 |
^ А3А4 |
|
|
У -v t . У 4 |
|||
|
|
|
|
|
|
|
||
|
V* V* |
V * X X |
/ |
j r 2 |
|
|
|
|
|
7 j Л2Л4 |
7 j A 3A4 |
,Л4 |
|
|
|
|
|
|
У А2ХаУ Х2Х4 |
|
|
У х лх яу х2х4 |
||||
|
У Х„Х. У X2. |
+ У А2Х4 |
У х \ |
V х,х4 |
||||
|
1 |
J-Л |
" |
|||||
|
|
vw3 4 |
4 |
|
|
|
3 |
3 4 |
Получаемые в конечном итоге определители второго порядка ре шаются по схеме крест-накрест, т. е.
У4 У л'за4 =У4У 4-(У-
УА' з А'4 у 4
Покажем на примере предыдущей системы уравнений поря док определения коэффициентов уравнения регрессии.
Определитель знаменателя |
|
|
|
|
|
|
|||||||
|
|
3,3 |
|
|
8,8 |
|
15 |
|
— 257,2 |
— 00 |
|||
Дэнам = |
8,8 |
- |
257,2 |
- 5 0 |
= |
|
|||||||
3,3 |
|
100 |
|||||||||||
|
|
15 |
- 5 0 |
|
100 |
|
- 5 0 |
||||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
8,8 |
15 |
+ |
15 |
|
8,8 |
15 |
|
||
|
|
|
- 5 0 |
100 |
|
|
|
-2 5 7 ,2 |
- 5 0 |
|
|||
= 3,3 [ ( _ |
257,2) ЮО- |
( - |
50) ( - |
50)] - |
8,8 [8,8 •100 - |
15 ( - 50)] + |
|||||||
|
|
+ |
15 [8,8 ( - 5 0 ) |
— 15 ( — 257,2)]; |
|
||||||||
|
|
|
|
|
Дэнам= - 5 6 |
|
200. |
|
|
||||
Определитель числителя для 612.34 |
|
|
|
||||||||||
|
|
106,2 |
|
8,8 |
|
15 |
= |
106,2 |
-2 5 7 ,2 - 5 0 |
||||
Дчисл |
- |
490,5 |
- 257,2 - 5 0 |
||||||||||
- 5 0 |
100 |
||||||||||||
|
|
150 |
|
- 5 0 |
|
100 |
|
|
|
|
|||
|
- ( - 4 9 0 ,5 ) |
|
8,8 |
15 |
+ |
150 |
8,8 |
15 |
|||||
|
|
|
|
- 5 0 |
100 |
|
|
-2 5 7 ,2 - 5 0 |
27