Файл: Алабин М.А. Корреляционно-регрессионный анализ статистических данных в двигателестроении.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.

(в)

Для нахождения значений коэффициетов 6124; 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.242,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