Моделирование и оптимизация объектов и процессов
..pdf
Глава 6. Модели I типа |
|
|
|
|
|
|
83 |
|
торе n, толщина b градиентных (метпллических) колец, толщина |
||||||||
диэлектрика в одной секции d. |
|
|
|
|
||||
Для определения оптимального числа секций в секционирован- |
||||||||
ном изоляторе найдем максимум значения U. Для этого продиф- |
||||||||
ференцируем выражение (6.17) и приравняем значение производной |
||||||||
к нулю: |
|
|
|
|
|
|
|
|
∂U |
= k(1 − α)n−α(H − bn)α − αkn1−αb(H − bn)α−1 = 0. |
|||||||
∂n |
||||||||
Преобразовав это выражение, получим |
|
|
|
|||||
∂U |
|
|
|
|
αbn |
(6.18) |
||
∂n = (1 − α) − αnb(H − bn)−1 = 1 |
− α − H |
− |
bn = 0. |
|||||
|
|
|
|
|
|
|
|
|
Из уравнения (6.18) после преобразования получим |
|
|||||||
|
|
|
|
H − αH = bn. |
|
|
(6.19) |
|
Из выражения (6.19) следует |
|
|
|
|
||||
|
|
N = nопт = H(1 − α) , |
|
|
(6.20) |
|||
|
|
|
|
|
b |
|
|
|
где N — оптимальное число секций в изоляторе. |
|
|
||||||
Коэффициенты k и α в выражении (6.15) для каждого конкрет- |
||||||||
ного случая можно рассчитать методом наименьших квадратов, ис- |
||||||||
пользуя экспериментальные значения, полученные при снятии за- |
||||||||
висимости пробивного напряжения U от толщины диэлектрика d. |
||||||||
Рассмотрим процесс определе- |
|
|
|
|
||||
ния указанных коэффициентов на |
|
|
|
|
||||
примере |
полиэтиленовых |
колец, |
|
|
|
|
||
помещенных в вакуум. |
Экспери- |
|
|
|
|
|||
ментальные значения |
|
пробивных |
|
|
|
|
||
напряжений U полиэтиленовых ко- |
|
|
|
|
||||
лец от их толщины d приведены в |
|
|
|
|
||||
табл. 6.3 и отображены на рис. 6.3 |
|
|
|
|
||||
кружками. |
|
|
|
|
|
|
|
|
Для |
упрощения |
нахождения |
|
|
|
|
||
коэффициентов k и α в формуле |
|
|
|
|
||||
(6.15) линеаризуем уравнение. Для |
Зависимость про- |
|||||||
этого прологарифмируем указанное |
бивного напряжения поверх- |
|||||||
выражение и получим |
|
|
|
ности полиэтиленового образ- |
||||
ln U = ln k + α ln d. |
|
(6.21) |
ца в вакууме от его толщины |
|||||
|
|
|
|
|
||||
Введем обозначения y = ln U, b0 |
= ln k, x = ln d. |
С учетом |
||||||
84
Таблица 6.3 Экспериментальные значения пробивных напряжений U
полиэтиленовых колец от их толщины d
n |
|
1 |
10 |
20 |
|
30 |
|
40 |
|
|
50 |
|
60 |
|
70 |
|
80 |
90 |
|
100 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U, кВ |
|
83 |
304 |
433 |
|
518 |
|
574 |
|
606 |
|
613 |
|
595 |
|
447 |
433 |
|
0 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d, мм |
|
297 |
27 |
12 |
|
7 |
|
4,5 |
|
3 |
|
|
2 |
|
1,29 |
|
0,75 |
0,33 |
|
0 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 6.4 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8 |
|
i |
|
1 |
|
2 |
|
|
3 |
|
|
|
|
4 |
|
5 |
|
6 |
|
7 |
8 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
d |
|
2 |
|
4 |
|
|
8 |
|
|
|
|
10 |
|
16 |
|
20 |
|
25 |
30 |
|
|
|||
x = ln d |
0,693 |
|
1,386 |
|
2,079 |
|
2,303 |
|
2,773 |
|
2,996 |
|
3,219 |
3,40 |
18,85 |
|||||||||
x2 = (ln d)2 |
0,480 |
|
1,922 |
|
4,324 |
|
5,302 |
|
7,687 |
|
8,974 |
|
10,36 |
11,6 |
50,62 |
|||||||||
U |
|
10 |
|
13 |
|
|
19 |
|
|
|
21 |
|
24 |
|
27 |
|
29 |
30 |
|
|
||||
y = ln U |
2,302 |
|
2,565 |
|
2,944 |
|
3,045 |
|
3,178 |
|
3,296 |
|
3,367 |
3,40 |
24,1 |
|||||||||
xy = ln d ln U |
1,596 |
|
3,555 |
|
6,121 |
|
7,012 |
|
8,813 |
|
9,874 |
|
10,84 |
11,6 |
59,38 |
|||||||||
Uтеор |
|
10,2 |
|
13,6 |
|
18,3 |
|
|
|
20,1 |
|
24,4 |
|
26,8 |
|
29,4 |
31,7 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
введенных обозначений уравнение (6.21) можно записать в виде |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
y = b0 + αx. |
|
|
|
|
(6.22) |
|||||||||||
Коэффициенты b0 и α в уравнении (6.22) с использованием метода наименьших квадратов могут быть определены по выражениям, взятых из [36]:
|
N |
N |
|
N |
N |
|
|
N |
yixi − |
N |
N |
|
|
|
b0 = ln k = |
yi |
xi2 − |
yixi |
2 |
xi |
N |
yi |
|
xi |
|||||
N |
|
N |
|
; α = |
N |
|
N |
2 . |
||||||
|
i=1 |
i=1 |
|
i=1 |
i=1 |
|
|
i=1 |
|
i=1 |
i=1 |
|
|
|
|
N i=1 xi2 |
− |
i=1 xi |
|
|
|
N i=1 xi2 |
− |
i=1 xi |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(6.23)
Для расчета численных значений коэффициентов b0 и α используем табл. 6.4.
Подставив численные значения в формулы (6.23), получим
|
|
N |
N |
|
|
N |
N |
|
|
|
b0 |
= ln k = |
yi |
xi2 − |
|
yixi |
2 |
xi |
|||
N |
|
|
N |
= |
||||||
|
|
i=1 |
i=1 |
|
|
i=1 |
i=1 |
|
|
|
|
|
N i=1 xi2 |
− |
i=1 xi |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
= |
24,098 · 50,618 − 59,378 |
· 18,85 |
= 2,0257. |
|||||||
|
|
49,621 |
|
|
|
|
|
|
|
|
Глава 6. Модели I типа |
|
|
|
|
|
|
|
|
|
|
|
|
|
85 |
||||
Проведя потенцирование, получим k = 7,648. Далее |
|
|
|
|||||||||||||||
α = |
N |
|
N |
N |
= |
|
· |
|
|
|
− |
|
|
|
· |
= |
||
|
|
2 |
|
|
|
|
|
|
|
|||||||||
|
N |
yixi − |
yi |
xi |
8 |
|
59,378 |
|
24,098 |
|
18,85 |
|
||||||
|
|
i=1 |
|
i=1 |
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
N |
|
|
|
|
8 |
· |
50,618 |
− |
355,323 |
|
||||
|
N i=1 xi2 |
− |
i=1 xi |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
= 0,4187. |
|
|
|
|
|
|
|
|
|
|
|||
На рис. 6.4 черными квадратиками обозначены расчетные значения пробивного напряжения поверхности полиэтиленовых колец, определяемого по выражению (6.15).
Окончательный вид уравнения (6.15) после подстановки в него численных значений k = 7,648 и α = 0,4187 примет вид
U = kdα = 7,648d0,4187. |
(6.24) |
Как следует из рис. 6.4, расчетные значения пробивного напряжения от толщины диэлектрика по формуле (6.24) практически полностью совпадают с экспериментальными значениями, т. е. формула (6.24) адекватно описывают эксперимент.
Задавшись высотой изолятора H = 300 мм и толщиной градиентного кольца b = 3 мм по формуле (6.17) был проведен расчет пробивных напряжений секционированных изоляторов при различной толщине d диэлектрических колец в секции. Результаты расчета приведены в табл. 6.4.
Как следует из табл. 6.4, для одной секции высотой H = 300 мм,
вкоторой было одно градиентное кольцо толщиной b = 3 мм и одно кольцо из полиэтилена толщиной d = 293 мм, расчетное пробивное напряжение составило 83 кВ. По мере увеличения числа секций, которому сопутствовало снижение толщины диэлектрика d в каждой секции, происходило повышение пробивного напряжения изолятора. Максимального значения пробивное напряжение достигало при числе секций 60. При дальнейшем увеличении числа секций происходило снижение пробивного напряжения. Это снижение было обусловлено тем, что при дроблении в изоляторе начиналось преобладание металлической составляющей, из-за возрастания в нем числа градиентных колец. Расчет оптимального число секций N
всекционированном изоляторе высотой H = 300 мм и толщиной градиентного кольца b = 3 мм по формуле (6.20)
N |
|
n |
|
H(1 − α) |
|
300(1 − 0,4187) |
, |
|
. |
|
|
= |
|
опт = |
b |
= |
3 |
|
= 58 13 |
≈ 58 |
|
86
Толщина одной секции d при оптимальном количестве секций в изоляторе будет равна
d = H − bnопт = 300 − 3 · 58 = 2,17 мм.
nопт 58
Расчетное пробивное напряжение такого секционированного изолятора определим по формуле (6.15):
H − bn α n
Рассчитанное предельное значение напряжения Uд = 613,6 кВ указывает на то, что данный изолятор может быть использован в установках, рабочее напряжение которых не превышает напряже-
ние Uд. |
В соответствии с |
проведенными |
|
|
|||
|
расчетами был выполнен секциониро- |
||
|
ванный изолятор из полиэтиленовых |
||
|
плоских колец, |
внутренний диаметр |
|
|
которых составлял 200 мм, а внешний |
||
|
300 мм, разделенных |
градиентными |
|
|
кольцами такого же размера, выпол- |
||
|
ненными из нержавеющей стали. Изо- |
||
|
лятор был испытан на электрическую |
||
|
прочность на вакуумной установке, |
||
|
приведенной на рис. 6.4 [63]. |
||
|
Проверка предлагаемой методики |
||
|
определения оптимального числа сек- |
||
|
ций в проходном изоляторе и сравне- |
||
|
ние её с аналогами осуществлялось |
||
|
следующим образом. |
Первоначально |
|
|
рассчитывалось |
оптимальное число |
|
|
секций в изоляторе по предлагаемой |
||
|
методике (6.24) и по методикам, пред- |
||
|
лагаемых аналогами, и рассчитыва- |
||
|
лось пробивное напряжения всей кон- |
||
Испытательная ва- |
струкции изолятора по формуле (6.24). |
||
куумная высоковольтная уста- |
После этого делались макеты изолято- |
||
новка |
ра, габариты и секции в которых были |
||
выполнены в соответствии с проведенными расчетами. После изготовления макетов они испытывались на пробивное напряжение на высоковольтной пробивной установке [38], приведенной на рис. 6.4. Установка включала в себя вакуумную камеру и проходной изоля-
Глава 6. Модели I типа |
87 |
тор. Вакуумная камера представляла собой цилиндр, выполненный из нержавеющей стали. Высота цилиндра равнялась 750 мм, а его диаметр 1200 мм. Высота изолятора была равна 2870 мм.
Внутренний диаметр верхнего и нижнего оснований были равны соответственно 400 и 600 мм. Верхняя часть изолятора заглушалась крышкой, имеющей связь с токоведущей трубой. Источник напряжения представлял собой генератор импульсных напряжений (ГИН) выполненный по схеме Аркадьева–Маркса. ГИН был рассчитан на напряжение 3 миллиона вольт.
При испытании на указанной установке изготовленного макета изолятора, выполненного в соответствии с формулой (6.17) пробой произошел при напряжении U = 615 кВ, что свидетельствовало об адекватности созданной нами методики расчета оптимального числа секций в изоляторе.
Для сравнения результатов, полученных выше с использованием предлагаемой нами методики, с результатами аналогов был произведен расчет оптимального числа секций в изоляторе по этим аналогам. Используя методику-аналог, предложенную в [35], определяли оптимальное число секций по выражению (6.15):
N = |
H ln En/Ei |
= |
|
300 · ln 1,116/0,278 |
|
= |
||
|
|
|
||||||
|
b[ln n + ln En/Ei] 3[ln 12 + ln 1,116/0,278] |
|
||||||
|
300 · 1,39 |
|
, |
|
. |
|
|
|
|
= 3[2,485 + 1,39] |
36 |
|
|
||||
|
= 35 87 ≈ |
|
|
|
||||
Толщина dспрот одной секции изолятора, выполненного в соот-
ветствии с аналогом, будет равна |
|
|
|
|
|
|||||||
|
|
dспрот = |
H − bN |
= |
300 − 3 |
· 36 |
= 5,3 мм. |
|
||||
|
|
|
N |
|
|
|
||||||
|
|
|
|
|
|
|
36 |
|
|
|
|
|
Пробивное напряжение секционированного изолятора, выпол- |
||||||||||||
ненного по аналогу [61] будет равно |
|
|
|
|
|
|||||||
U |
прот = |
k |
H − bN |
α N |
, |
, 0,4187 |
, |
. |
||||
|
|
N |
= 7 648 |
· 5 |
3 |
|
· 36 = 553 5 кВ |
|
||||
Как следует из приведенных цифр, предлагаемая в этом разделе методика, по сравнению методикой, изложенной в [35], позволяет повысить пробивное напряжение приведенного в примере изолятора высотой H = 300 мм и толщиной градиентной прокладки b = 3 мм на 11,7 %.
В [34] авторы в примере конкретного выполнения утверждают, что при высоте изолятора H = 385 мм и толщине прокладок из электропроводящего материала b = 5 мм оптимальное число секций
88
N = 10. Пробивное напряжение такого изолятора можно определить по формуле (6.17):
U |
|
k |
H − bN α |
N |
, |
|
385 − 5 · 10 0,4187 |
, |
. |
||
|
прот = |
|
|
|
|
= 7 648 |
|
|
|
· 10 = 332 7 кВ |
|
|
N |
|
10 |
|
|
||||||
В соответствии с предлагаемой нами методике по формуле (6.20) оптимальное число секций в изоляторе высотой H = 385 мм и толщине прокладок из электропроводящего материала b = 5 мм будет равно
N |
= |
n |
опт = |
H(1 − α) |
= |
385(1 − 0,4187) |
≈ 48 |
. |
|
|
|
b |
5 |
|
|
||||
Пробивное напряжение такого изолятора в соответствии с формулой (6.17) будет равно
U |
|
k |
H − bN α |
N |
, |
|
385 − 5 · 48 0,4187 |
, |
. |
||
|
прот = |
|
|
|
|
= 7 648 |
|
|
|
· 48 = 583 1 кВ |
|
|
N |
|
48 |
|
|
||||||
Таким образом, изолятор, оптимальное число секций в котором рассчитано по нашей методике, имеет в 1,75 раз более высокую пробивную прочность изолятора, чем изолятор, оптимальное число в котором определено по методике, предлагаемой в [34]. Рассмотренный подход к определению оптимального числа секций в изоляторе подробно отражён в патентах [36, 37] и в статье [38].
Таким образом, предложенный способ определения оптимального числа секций в проходном изоляторе по сравнению со способомпрототипом имеет более высокую точность определения оптимального числа секций в изоляторе, что позволяет при заданной высоте изолятора H и заданной толщине градиентной прокладки b получить максимально возможное пробивное напряжение для указанных габаритов изолятора. Кроме того, предлагаемый способ, по сравнению со способом-прототипом, достаточно прост и не требует для своего воплощения, особенно для изоляторов с большими габаритами, изготовления испытательного стенда со сверхвысоким источником напряжения, что не всегда возможно воплотить в жизнь. В частности, в рассмотренных нами примерах для реализации заявляемого способа достаточно того, чтобы экспериментальная высоковольтная установка позволяла получать напряжение не превышающего 30...35 кВ. Тогда как для реализации способа-прототипа необходимо, как следует из табл. 6.3, чтобы экспериментальная установка имела возможность экспериментировать при напряжениях 450...500 кВ, что более чем на порядок выше напряжения, необходимого для реализации заявляемого способа.
Глава 6. Модели I типа |
89 |
Предлагаемый способ позволяет, не меняя габаритов изолятора путем определения оптимального числа секций, повысить его электрическую прочность не менее чем на 11,7 %.
Кроме того, при заданном рабочем напряжении изолятора выбором оптимального числа его секций можно добиться сокращения габаритов и уменьшения стоимости изолятора.
КОРРЕЛЯЦИОННЫЙ АНАЛИЗ
7.1. Оценка тесноты нелинейной связи
Если считать, что уравнение регрессии найдено с достаточной степенью точности, то остаточная дисперсия обусловлена только наличием дисперсии воспроизводимости Sост2 ≈ Sвоспр2 . Чем меньше доля Sост2 в общей дисперсии Sy2, тем сильнее связь между y и x, ибо тем меньше случайности в этой связи. Поэтому силу связи можно характеризовать величиной
ξ = |
(N − l)Sост2 |
. |
|
|
|
|
(N − 1)Sy2 |
|
|
|
|||
|
√ |
|
|
, 0 |
θ 1, |
|
Связь тем сильнее, чем меньше ξ. Величина θ = |
1 |
− ξ |
||||
называется корреляционным отношением. Чем больше θ, тем сильнее связь.
Если θ = 1, существует функциональная зависимость между параметрами. Только в случае нормального распределения равенство θ = 0 однозначно свидетельствует об отсутствии связи между случайными величинами.
Корреляционное отношение, как и коэффициент корреляции в
линейной регрессии, характеризует тесноту связи между случайными величинами. Анализ силы связи по θ называется корреляцион-
ным анализом. В случае линейной регрессии корреляционное отно-
шение равно коэффициенту корреляции |
|
|||||
|
|
|
|
|
|
|
|
− N − 1 Sy2 |
|
||||
θ = 1 |
N − 2 |
|
Sост2 |
|
= r . |
|
|
|
|||||
7.2. Метод множественной корреляции
Если необходимо исследовать корреляционную связь между многими переменными, пользуются уравнением множественной регрессии
Y = b0 + b1x1 + b2x2 + b3x3 + . . . + bkxk. |
(7.1) |
Глава 7. Корреляционный анализ |
91 |
Таблица 7.1
№ опыта |
X1 |
X2 |
. . . |
Xk |
Y |
|
|
|
|
|
|
1 |
x11 |
x21 |
. . . |
xk1 |
y1 |
2 |
x12 |
x22 |
. . . |
xk2 |
y2 |
. . . |
. . . |
. . . |
. . . |
. . . |
. . . |
N |
x1N |
x2N |
. . . |
xkN |
yN |
Таблица 7.2
№ опыта |
X0 |
X0 |
. . . |
X0 |
Y 0 |
|
1 |
2 |
|
k |
|
1 |
x110 |
x210 |
. . . |
xk01 |
y |
2 |
x0 |
x0 |
. . . |
x0 |
y0 |
|
12 |
22 |
|
k2 |
2 |
. . . |
. . . |
. . . |
. . . |
. . . |
. . . |
N |
x0 |
x0 |
. . . |
x |
y0 |
|
1N |
2N |
|
|
N |
Здесь мы уже имеем дело не с линией регрессии, а с поверхностью
регрессии при k = 2 и с гиперповерхностью при k > 2. В общем случае эта поверхность называется поверхностью отклика.
При построении поверхности отклика на координатных осях факторного пространства откладываются численные значения параметров (факторов). Исходный статистический материал представляют в виде табл. 7.1.
Прежде всего, перейдём от натурального масштаба к новому, проводя нормировку всех значений случайных величин по форму-
лам: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yj0 = |
yj − |
|
|
; xji0 = |
xji − |
|
|
j |
, |
|
|
|
|
|
|
|
|
|||||||||||||
y |
x |
|
i = 1, 2, . . . , N; j = 1, 2, . . . , k, (7.2) |
|||||||||||||||||||||||||||
|
||||||||||||||||||||||||||||||
|
|
Sxj |
|
|
|
|
|
|||||||||||||||||||||||
|
|
Sy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
где y0 |
, x0 |
— нормированные значения соответствующих факторов; |
||||||||||||||||||||||||||||
j |
|
ji |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
ji0 — среднеквадратичные отклоне- |
||||||||||
средние значения факторов; |
|
|
||||||||||||||||||||||||||||
y |
x |
|||||||||||||||||||||||||||||
ния факторов от средних значений: |
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
; Sxj = |
|
|
|
|
|
|
|
|
|
|||||||||||
Sy = |
|
1 |
|
N (yi |
− |
|
|
)2 |
|
|
1 |
|
N (xji |
− |
|
j )2. |
||||||||||||||
|
|
|
y |
x |
||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
N |
− |
1 i=1 |
|
|
|
|
|
|
|
|
|
N |
− |
1 i=1 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
В табл. 7.2 приведен исходный статистический материал в но-
вом масштабе. В новом масштабе |
|
|
j0 = 0; |
|
|
0 = 0; S20 |
= 1; S20 = 1. |
|||||||||||||||||||||
x |
y |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
j |
|
|
y |
j |
|
Подставив в выражение |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
N |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
rхy = (N |
− |
1)SxSy |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
(xi − x)(yi − y) |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|||||||||
новые переменные, получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
1 |
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
N |
|
|
|
rxj0xm0 = |
|
|
|
|
|
|
|
|
(xj0 |
− |
|
|
j0)(xm0 − |
|
|
m0 ) = |
|
|
|
xji0 |
xmi0 ; |
|||||||
|
|
|
|
|
|
|
|
|
x |
x |
|
|
|
|||||||||||||||
(N |
− |
1)S |
0 S 0 |
|
N |
− |
1 |
|||||||||||||||||||||
|
|
|
x |
j |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
m i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(7.3) |
|
|
1 |
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
rxj0y0 = |
|
|
|
xij0 yj0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(7.4) |
|||||
N |
− |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вычисленный по формулам (7.3) и (7.4) выборочный коэффи-
92
циент корреляции равен коэффициенту корреляции между переменными, выраженными в натуральном масштабе. Уравнение регрессии в новых переменных не имеет свободного члена. Для того чтобы показать это, вернемся к выражению (7.1):
Y = b0 + b1x1 + b2x2 + b3x3 + . . . + bkxk.
Не записывая всю систему уравнений Гаусса для этого уравнения регрессии, выделим лишь первое из них, которое имеет вид
|
|
N |
|
|
|
|
N |
|
N |
N |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Nb0 + b1 |
x1i + b2 x2i + . . . + bk |
|
|
xki = yi. |
(7.5) |
|||||||||||||||
|
|
i=1 |
|
|
|
i=1 |
i=1 |
i=1 |
|
|||||||||||
Разделим обе части уравнения (7.5) на N, получим |
|
|||||||||||||||||||
|
|
|
b0 + b1 |
|
|
+ b2 |
|
+ . . . + bk |
|
= |
|
|
|
|
|
(7.6) |
||||
|
|
|
x1 |
x2 |
xk |
y. |
|
|
|
|||||||||||
Вычтем из (7.1) уравнение (7.6), получим |
|
|
|
|
||||||||||||||||
Y − |
|
= b1(x1 − |
|
) + b2(x2 − |
|
) + . . . + bk(xk − |
|
). |
(7.7) |
|||||||||||
y |
x1 |
x2 |
xk |
|||||||||||||||||
Помножим каждый из членов уравнения (7.7) и одновременно разделим на одну и ту же соответствующую данному члену средне-
квадратическую величину: |
|
|
|
|
|
|
|
|
|||||||||||||
|
Sy(Y − |
|
|
) |
|
Sx1 b1(x1 − |
|
|
) |
|
Sxk bk(xk − |
|
) |
. |
|
||||||
|
y |
= |
x1 |
+ . . . + |
xk |
(7.8) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
Sy |
Sx1 |
Sxk |
|
|
|||||||||||||||
Учитывая, что по формулам нормирования (7.2) |
|
|
|||||||||||||||||||
|
|
Y − |
|
= y0; |
x1 − |
|
= x0; . . . ; |
xk − |
|
= x0 |
|
|
|||||||||
|
y |
x1 |
xk |
, |
|
||||||||||||||||
|
|
|
|
|
|||||||||||||||||
|
|
Sy |
1 |
|
|
|
|
k |
|
|
|||||||||||
|
|
Sx1 |
Sxk |
|
|
||||||||||||||||
уравнение (7.8) принимает вид |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
Y 0 = a1x10 + a2x20 + . . . + akxk0 , |
|
(7.9) |
|||||||||||||
где aj = (Sxj /Sy)bj. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Коэффициенты уравнения (7.9) находятся из условия |
|
||||||||||||||||||||
|
|
|
|
|
|
|
Ф = |
N |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
(Yj0 − Yj0 )2 = min . |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i=1
Условия min функции определяются так же, как и в случае
зависимости от одной переменной: |
|
|
|
|||||
|
∂Ф |
= 0; |
∂Ф |
= 0; |
∂Ф |
= 0; · · · ; |
∂Ф |
= 0, |
|
∂a1 |
∂a2 |
∂a3 |
∂ak |
||||
