Файл: Пивоваров, С. Э. Моделирование процессов прогнозирования в приборостроении.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)

скалярно

на

вектор

 

 

 

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 С, Э, Пивоваров