Файл: Пивоваров, С. Э. Моделирование процессов прогнозирования в приборостроении.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.10.2024
Просмотров: 78
Скачиваний: 0
Умножим обе части уравнения (3.4.9) на f, ^/ = 0, 1, . . . .
скалярно. Так как (h, | / ) = 0 при 1 Ф /, то
(У, Ь ) = |
С ;(|/, 1 ) |
+ (А, |
Ь ) , |
(/ = 0, |
(3.4.10) |
Отсюда получаем: |
|
|
|
|
|
|
С’,= |
|
|
& . |
(3.4.11) |
|
(th h) |
(h, h) |
|
||
Определим математические |
ожидания |
случайных |
величин |
||
с |
(Д , %i) |
____ . |
. Зга- f - n 2 |
|
|
б' “ |
С В |
для |
|
|
|
Обозначим математическое ожидание случайной величины симво лом М, получим
(3-412>
так как МД = 0 . Таким образом,
а .
(£/> £//
Отсюда получим формулу для оценки коэффициентов
с; {1 = о, 1......... |
Щ * \ . |
|
|
(у, Г,) |
(3.4.13) |
||
(£/» |
ъ) |
||
|
Погрешности оценки коэффициентов С\ можно определить, если вы числить дисперсии случайных величин:
* & Гг) .
|
|
{р = 1 |
[ ( I |
р = 1 |
< З Л 1 4 ) |
|
|
6, = C f-C ,.
&
Строим для оценки значений дисперсии а2 новый А^-мерный
вектор |
р по формуле |
|
|
3п+ п2 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
( Y , |
Q |
h. |
|
|
|
|
|
|
|
|
|
|
|
’ - |
у - |
1=2о (lh ll) |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Преобразуя |
вектор |
р по формулам (3.4.9) и (3.4.11), получим: |
||||||||||||||
|
3п +п2 |
|
|
|
Зп + п1 |
|
|
|
|
3п+ я 2 |
|
|
|
|||
|
|
V |
(уЛ,) |
|
|
|
|
|
|
|
|
|
(у, id |
|
||
Р = |
1^ Q |
|
ll) |
s- 2 № +*- 2 М б - |
||||||||||||
|
|
|
|
|
1= 0 |
(|/. ъ) |
|
|
1=0 |
(lh ll) |
|
|||||
|
|
Зп+ п2 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
= А - |
|
|
|
3я + я 2 |
|
|
3л + я 2 |
|
|
|||||||
1=02 |
(lh |
ll) |
|
|
2 |
|
|
|
2 ст-саб. |
|||||||
|
|
V |
(Д. |
ь) |
|
|
|
|
|
|
|
|
|
|
|
|
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3n + |
я 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
p = A - 2 |
6,f,. |
|
|
|
|
|
(3.4.15) |
|||
|
|
|
|
|
|
|
|
1=0 |
|
|
|
|
|
|
|
|
Домножив обе |
части |
выражения |
(3.4.15) |
скалярно |
на |
вектор |
||||||||||
1т |
|
|
|
2 |
|
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3/1 + |
я2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
(р, D |
= |
(A, |
L ) - \ |
2 |
6&, |
L / |
= (A,fm) - ( A , f m) = 0. |
(3.4.16) |
||||||||
|
|
|
|
|
|
/= |
0 |
|
|
|
|
|
|
|
|
|
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
Зп + |
я 2 |
|
|
|
|
|
|
|
|
(р. Р) = |
(р, |
А) = |
(А, |
А) - |
21 |
(&• |
А); |
|
|
||||
|
|
|
|
|
|
|
|
|
|
г = 1 |
|
|
|
|
||
Поэтому |
|
|
|
(С/. А) = б, (Ь, Q . |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
Зя + |
я 2 |
|
|
|
|
3/1 + я 2 |
|
|
|||
(?, |
Д - ( д , |
4 ) - |
2 «КБ, § ) - |
N |
|
|
2 |
|
«?©,£). |
|||||||
2 Д !р - |
|
|||||||||||||||
|
|
|
|
|
1= 1 |
|
Р= 1 |
|
/ = 1 |
|
|
|||||
Вводим случайные величины |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
Ч, |
( / = 1, |
2, |
|
2 |
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V(£ь £/)• |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
Ч/ — |
|
|
|
|
|
|
62
Тогда параметры нормальной величины тр суть
Мчр = |
0; |
|
|
Dr\i = а2. |
|
||
Отсюда |
|
3п + п 2 |
|
|
|
|
|
(р, р) = 2 J Д2Р “ |
2 |
Ль |
|
Р = 1 |
|
;= i |
|
где |
|
|
|
Л/ = 8/ К ( & . 9 = |
2 |
Ар |
(3‘4Л7) |
(£/. I/ |
р |
= 1 |
1/ 2 |
|
|||
|
|
|
' р=| |
Продолжим ортонормированное |
|
преобразование для (/ = |
— Зп + ? ! 1,
yv (любым способом). Тогда все г]р (р = 1, 2,
N) суть нормальные случайные величины с параметрами (0, а2). Поскольку при ортонормированием преобразовании сохраняется инвариантность расстояния (поворот координатных осей в прост ранстве), имеет место соотношение
|
|
Зп+ |
п2 |
^ |
N |
2 |
N |
|
(р> р) = |
2 |
Лр - 2 л? = 2 |
Лр, |
|
||
|
р = 1 |
/ = 1 |
|
_ _ Зя + л 2 , , |
|
||
|
|
|
|
Р |
2 |
i |
|
|
|
3n + /i2 |
\ |
|
|
|
|
так • как сумма -квадратов {^N ■---- 2— |
) одинаково распределенных |
||||||
нормальных |
случайных |
величин |
распределена |
по |
г_!±±£! и |
||
М [(?, » ] = |
( « - |
|
|
|
|
|
|
М(р. р )
УЗя+ Ра
2 |
|
|
Следовательно, |
|
|
(р. Р) |
~ О2. |
(3.4.18) |
З я + п 2 |
||
д/------- |
|
|
В результате можно обосновать метод отбора значимых факторов, существенно влияющих на выходную функцию. Для этого строим
3п-\-п2
— дисперсии:
(рь Pi) |
о?; |
3 n + |
rfi \ |
(3.4.19) |
N — I |
2 |
) ’ |
63
где
|
(3.4.20) |
Если при переходе от I к I + |
1 факторам, т. е. при добавлении |
(/ + 1)-го фактора, величины а) |
и 07 + 1 несущественно отличаются |
друг от друга, то это означает, что (/ -+- 1)-й фактор не оказывает на выходную функцию непосредственного влияния и может быть отброшен.
Для проверки статистических гипотез о значимости расхожде ния двух дисперсий о] и crf |_ 1, взятых из нормальных генеральных совокупностей, используем распределение отношений несмещенных оценок дисперсий.
(3.4.21)
Распределение случайной величины F (так называемое распре деление Фишера) зависит только от числа степеней свободы /q = = к2 = N — 1, если выборка, т. е. число наблюдений, имеет объем, равный N. Чтобы проверить гипотезу о равенстве сх| и о'{+ ь построим критическую область для критерия F. После сравнения эмпириче ских значений F с величинами этой области делаем окончательный вывод об их равенстве.
Будем задавать коэффициент доверия р и фиксировать два зна чения Fy и F2, таких, что при соответствующих степенях свободы кг и к2
(3.4.22}
Если выборочное значение F оказывается в критической области, т. е. вне области значений (Flt F2), то гипотеза ст| = crf+i отвергается, так как вероятность того, что F не содержится в интервале (Flf F2)
равна |
2 ( 1 — р) и весьма незначительна. |
|
|
Если |
|
В качестве коэффициента доверия р принимаем р — 0,95. |
|||||
значение F не выходит за пределы интервала (Е,, F2), т. е. FYsg |
F eg: |
||||
C F2, считаем, что (/ + 1)-й фактор оказывает несущественное влия |
|||||
ние и может не учитываться. При F =sg Fx или F ^ |
F2 расхождение |
||||
между о'} и (Tf+ i считаем существенным, т. е. фактор (/ + |
1), |
как |
|||
оказывающий наибольшее влияние, необходимо учитывать. |
= |
1 |
до |
||
В |
результате, оценивая изменения дисперсий |
aj от / |
|||
I |
+ п%' можно отобрать все значимые факторы уг. |
|
|
|
Отбор влияющих факторов методом сравнения остаточных дис персий позволяет объективно оценить значимость факторов с по мощью имеющихся статистических данных, причем трудоемкость вычислений значительно уменьшается по сравнению с вычислением обычным методом наименьших квадратов.
64
По формуле (3.4.11) рассчитываем значения С\ только для зна чимых факторов. Для получения итоговых значений Ct применяется рекуррентное соотношение вида
Ст = С'т\ |
|
|
|
Cm—к— Ст—к ' |
* |
(gm_ K, V/ ) |
(3.4.23) |
|
|||
. } =т2—к 1 |
к> \ т —к) |
|
где т — число наиболее влияющих факторов, отобранных с уче том попарных произведений нелинейных составляющих.
Если учесть, что случайная компонента А есть мера случайного влияния множества случайных факторов, характеризующих вы ходную функцию У, то вероятность того, что истинное значение У отклонится от вычисленного по формуле (3.4.5) не более чем на ko, определяется по формуле
( |
ш, |
m, |
\ |
k |
_ /> |
Со+ |
2 |
Cfti —ko <;М <гС0+ ^ |
Cfli -f- ko 1 = у ^ - |
^ е 2 dt. |
|
|
|
|
|
° |
(3.4.24) |
Эта формула определяется из центральной предельной теоремы Ля пунова.
Оценки коэффициентов С), определяемые по формуле (3.4.13), являются приближенными, а следовательно, и оценки, получаемые линейной комбинацией значений С}, также приближенные. Отсюда целесообразно оценить меру точности определения коэффициентов С;. Среднее квадратичное отклонение коэффициентов С/ приближенно равно
О |
£i |
(3.4.25) |
а (Ci) |
h ) C’l |
|
V |
|
После оценки коэффициентов зависимостей выходных функций У (основных технико-экономических показателей) от наборов наибо
лее влияющих факторов xt (i — 1, т) и выбора окончательного вида этих зависимостей определяем вариантные значения У в прогнози руемом периоде. При этом задаем различные значения для влияю щих факторов в диапазоне граничных значений, отражающих воз можности функционирования исследуемой системы. Для определе ния граничных значений по факторам используются директивные указания и различные вариантные проработки, учитывающие пессимистические и оптимистические значения темпов и пропор ций экономического развития.
3 С, Э, Пивоваров