Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Пивоваров, С. Э. Моделирование процессов прогнозирования в приборостроении

.pdf
Скачиваний:
13
Добавлен:
19.10.2023
Размер:
6.89 Mб
Скачать

Коэффициенты а0, аи .... аг вычисляются из условия: i

Q = 2 [(а0+ аМ + . . . + aru'r) - Р’\ = min,

/=1

где щ — значение

г-го фактора, Р'

— значение

потребности в /-м

наблюдении. Величины а0,

аг,

..., аг, при которых величина Q до­

стигает

минимума,

определяется

из системы т уравнений

т. е.

 

1

=

0

0 - 0

.

1......... *

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

^

( ао + а1х{ + --- + аг х 1 - Р ') - 1= 0 ;

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

Y ^ =

2

+

QlAr‘ +

' ' • +

arX>r ~

p/)*i =

° :

 

dQ

^

(ао +

а

+

• ■• +

Я'*/ -

РО

=

0.

 

2 даг; =

 

 

/=1

 

 

 

 

 

 

 

 

 

 

Введя

обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ьц 1;

Ьь г+2 2

Р^\

 

 

 

 

 

 

 

 

 

 

 

/= I

 

 

 

 

 

 

 

 

/= 1

 

 

(А =

2,

r +

l);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

bvi=

2

'

 

 

 

 

 

 

 

 

 

 

 

/=i

 

 

 

(и = 2, г + 1);

 

 

 

 

 

 

 

 

 

 

 

*4 /-+2 =

 

 

 

 

 

 

 

 

 

 

I

 

 

/= 1

_______

 

_______

 

 

 

_

 

 

 

6 **=

 

 

 

(у = 2,

г + 1 ;

6 = 2,

/г +

1).

 

/ = 1

 

 

 

 

 

 

 

 

 

 

 

получим систему

 

 

 

 

 

 

 

 

 

 

 

 

 

^11^0 +

^12а1 +

• • +

^1, r+l^r — ^1, л+2>

 

 

^21а0+ ^22а1 + • • + ^2, /-+1а/- = ^2, /Ч-2>

 

 

b r+i, i«o +

++*, га г +

• +

+ + J, /-++(- — br+l, r+p

50

для решения которой можно применить один из методов решения систем линейных алгебраических уравнений.

Прогнозирование потребности в каком-либо виде продукции на длительный срок должно сопровождаться, с нашей точки зре­ ния, построением нескольких вариантов моделей и отбором опти­ мального варианта, наиболее точно отражающего развитие потре­ бления за базовый период и учитывающего возможные структур­ ные изменения в прогнозируемом периоде.

Ниже изложена методика выбора оптимального варианта мо­ дели. Рассмотрим несколько зависимостей из перечисленных выше (можно выбрать и другие виды зависимости, например вместо

факторов xlt

х2,

хг использовать их степени х* x*t

xKrt где

к ± 1, ± 2,

...).

 

корреляции

Следует

отметить, что в моделях множественной

по сравнению с моделями парной корреляции выбор вида зависи­

мости много сложнее, так как действия факторов

переплетаются

и отсутствует возможность графического контроля.

В данном слу­

чае большое значение имеет исследование зависимости каждогр фактора и выходной функции.

Далее сравним численные значения потребности в годы базо­ вого и прогнозируемого периодов. Последние получим из формулы

 

р(Л = ф (дл/^ х(р t

х</>)

 

где / — го£

прогнозируемого

периода,

считая значения х(р, х</>,

..., х® известными.

 

 

 

С помощью ^-статистики

Стьюдента проверим,

принадлежат

ли эти две

независимые

частичные

совокупности

объемов

(длина базового периода) и /2 (длина прогнозируемого периода) одной и той же нормально^ распределенной общей совокупности,

имеющей среднее значение р = 0 и дисперсию а2.

Пусть средние значения рассматриваемых частичных совокуп­ ностей равны рг и р2, а для оценок дисперсий справедливо

s;- 2

/.=1

Г л ( h ) - P il2

и - 1

: i i .

[Рг (/2) — Рг\2

/2- 1

/.=i

и пусть проверяемая гипотеза верна. Основой проверки гипотезы является разность рх р2, дисперсия которой равна [51]

<7а

/ + А

0^- •Ра т + £ -

hk

51

Так как оценки SJ и 5 | дисперсии ст2 имеют вес 4 — 1 и 4 — 1,

то полная оценка дисперсии

о2 будет равна:

It

 

U

02 _ (4—1) ■St + (4—О *S|

2 (Pi ( i i — Pil2 +

(P'i (/2) —Pal2

1=j____________b=i__________

(h~ !) + (/,-!)

U+ 4 2

В данном случае для /-статистики справедливо [51]

i Pi Р2

S

Таким образом, можем воспользоваться таблицей распределения /-статистики, взяв число степеней свободы

v = 4 + 4-

Далее, из зависимостей, удовлетворяющих этому критерию, выбираем наилучшую в смысле наименьшей суммы квадратов отклонений расчетных значений от фактических, т. е. такую зави­ симость

Р = ф (*i> *2, .... хг),

для которой

2 [р{ - фМ> .... *0 Г = min-

/= 1

Модель множественной корреляции может иметь также вид

р = «О + a j (лу) + . . . + arf (хг) + А,

где f — какая-нибудь из функций в перечисленных выше видах зависимостей;

А — случайная нормально распределенная величина с нуле­ вым средним и дисперсией а2, характеризующая влия­ ние не поддающихся учету случайных факторов.

Реализация этой модели представляется более сложной, чем описанных выше.

Исследование

метода получения коэффициентов

щ (г = 0, г)

данной модели с

одновременной оценкой дисперсии

о2 описано

в следующем параграфе.

 

После выбора модели и получения её коэффициентов целесооб­ разным представляется рассмотрение частных коэффициентов эла­ стичности. Если зависимость потребности от влияющих на неё факторов имеет вид

Р = ф(*1, *2. •••. Хг),

то для коэффициента эластичности в зависимости от изменения фактора Xi (i = 0, г) справедливо

чдР

Кх , = р дх[

52

Коэффициенты эластичности характеризуют частные скорости роста потребности в зависимости от изменения рассматриваемого фак­ тора. Таким образом, получаем возможность сопоставить разно­ родные (измеряемые в разных единицах) факторы и расположить их по степени их влияния на величину потребности.

3.4. Прогнозирование экономического развития отрасли

Основные цели экономического прогнозирования заключаются в определении тенденций развития экономических процессов, в ра­ скрытии основных закономерностей этих тенденций, в выявлении возможных альтернатив развития отрасли, в определении темпов и пропорций расширенного воспроизводства.

При этом перед прогнозированием экономического развития от­ расли ставится задача исследования важнейших технико-экономи­ ческих показателей, характеризующих экономическое развитие отрасли, таких как объем валовой (товарной) продукции, себе­ стоимость валовой (товарной продукции), производительность труда, балансовая прибыль и др.

В процессе прогнозирования необходимо выяснить характер воздействия на названные выше показатели экономического раз­ вития различных факторов, в частности таких, как: результаты научно-технического прогресса, уровень организации и специали­ зации, капитальные затраты, трудовые и материальные ресурсы и др.

При этом исследуемый показатель следует рассматривать как выходную функцию совместного влияния факторов. Весь процесс экономического прогнозирования отрасли состоит из следующих этапов:

отбор наиболее значимых факторов; составление математической модели зависимости;

прогнозирование основных технико-экономических показате­ лей.

При решении задачи отбора наиболее значимых факторов к ним предъявляется ряд существенных требований, а именно: факторы должны быть выражены в количественных показателях (даже если они носят качественный характер); ни один из факторов не должен находиться в точной (функциональной) зависимости от другого или группы других.

После выбора выходного значения и установления набора фак­ торов составляется таблица (матрица) наблюдений, в которой строки представляют собой фиксированные значения выбран­ ных факторов в отчетные годы исследуемого периода.

Рассматривая выделенные основные технико-экономические по­ казатели как выходные функции от совместного влияния ряда факторов, можно при наличии матрицы наблюдений построить статистическую модель зависимости выходной функции от этих факторов, которая в известной степени будет решать задачу полу­

53

чения значений этих показателей в любой период времени. Таким образом, на основании статистической модели, используя методы математической статистики, отбираем факторы, наиболее влияю­ щие на выходную функцию, и строим математическую модель зави­ симости выходной функции от уже отобранных факторов.

Определение зависимости выходной функции от набора влияю­ щих факторов возможно при использовании метода наименьших квадратов, сущность которого заключается в следующем.

Пусть имеются результаты N независимых наблюдений,

в которых каждому набору значений М факторов xt (i = 1, М) соответствует значение функции У. Из соображений, отражаю­ щих специфику экономического развития исследуемой отрасли, выбираем зависимость

У= Ф(*Ъ х2, ..., хм),

содержащую ряд неизвестных параметров аъ а2, .... ап. Требуется

выбрать эти параметры

так, чтобы

кривая

 

У = Ф {хъ

х2, .... хм,

аг, а2,

ап)

наилучшим образом отображала полученную путем статистических наблюдений зависимость.

Метод наименьших квадратов основан на требовании свести к минимуму квадраты отклонений полученных статистических

значений от сглаживающей кривой,

т. е.

величина

N

 

 

 

5 = 2 [У7 — ф М ,

х 2,' ...,

х!м, аъ а2, .... ал)]2

/= 1

 

 

 

должна быть минимальной.

фактора;

У1— значение выходной

Здесь х[ — значение t'-ro

функции в /-м наблюдении.

Величины аи а2, ..., ап, при которых величина S достигает мини­

мума, определяются из системы

п уравнений:

£ = - 0 ,

(к = Т Г я ).

Основным недостатком метода наименьших квадратов является трудность решения задачи отбора наиболее значимых факторов. Так как известно, что на величину У оказывают совокупное влия­

ние М факторов xlt х2,

..., хм, то представляет большой интерес

не только построение

математической модели зависимости Y =

= ф (хг, х2, ..., хм), но и оценка степени значимости каждого из факторов в отдельности (отбор наиболее влияющих факторов).

Другим возможным способом определения зависимости выход­ ных функций от отобранных факторов, также часто встречающимся в практике статистического исследования экономических показа­ телей, является метод множественной корреляции, сущность ко­

54

торого заключается в том, что ставится задача определения зави­ симости вида

Y = F(x1, х2........

хм).

При этом функция F (модель Брандона) является произведе­ нием некоторых функций отдельных факторов:

Y = kh (Jfi) /s (хг) ...?м(хм)>

причем каждая из (х,) принимается для упрощения расчетов линейной

fi(xt) = ai + bixl.

Таким образом, задача сводится к определению. коэффициен­

тов ah bi и k (i = 1, М). Исходными данными для расчета служит таблица наблюдений, как и в методе наименьших квадратов.

В начале расчета строим график зависимости Y от хх. На гра­ фике получим корреляционное поле точек, из которых методом наименьших квадратов или методом средних определяется наиболее вероятная линия регрессии

fi(x1) = a 1+ b1xl

и находятся коэффициенты ах и bt. После этого рассчитывается зна­ чение нового приведенного фактора

Ух ~ ^

~

kf% (х%) ■• -fм (хм)

и строится график зависимости

от х2, по которому определяется

уравнение линии регрессии / 2 (х2). Далее рассчитывается значение следующего приведенного фактора

v -

Ух

У г

к ( ч )

Такой расчет продолжается до тех пор, пока не будут опре­ делены все функции fi (xi) и коэффициент к.

Недостатками модели Брандона являются трудность определе­ ния остаточной дисперсии погрешности оценки степени значимости факторов с помощью функции

Y = F(xlt х2 ... , хм)

и, как следствие, невозможность проведения отбора значимых фак­ торов.

Отбор факторов может быть проведен с помощью коэффициен­ тов корреляции. Для этого рассмотрим влияние факторов друг на

друга, рассчитывая коэффициенты парной корреляции (между любой парой из них) по формуле

 

 

Гщх! =

X t X j X i X j

 

 

а ,

а_

 

 

 

xi

xj

 

N

 

 

 

где Xi =

У xW

_ среднее

значение факторов хг;

к= 1

N

1 V (*) по — среднее значение произведения пары фак-

X#J = -N

2 i Xl х‘

торов W

 

 

 

Л'= 1

 

 

 

 

[xf0]2— (х,)2— дисперсия случайной величины хг.

Коэффициент парной корреляции гх.х.

характеризует тесноту

связи между факторами х,- и Ху, причем 0

| гх,х.| sg 1.

Чем |тх.х/|

ближе к 1,

тем теснее связь между х( и Ху,

при этом,

если гх/х/ =

± 1, между xt и Xj предполагается прямая функциональная зави­ симость. Наличие такой тесной связи между двумя факторами ука­ зывает на их полную взаимную зависимость, т. е. один из них нужно исключить из исследуемого набора факторов, но какой именно — определяется при последующем логико-априорном анализе специа­ листами.

Далее исследуем влияние каждого фактора на выходную функ­ цию, для чего определяем коэффициенты парной корреляции по формуле

 

 

 

_

Yxj-Yxi

 

 

 

Т Yx.

 

 

N

 

 

 

где У

к=1

 

— среднее значение выходной функции;

 

 

 

 

77,-12к = 1

Ymx^

— среднее значение произведения выход­

ной

функции и влияющего фактора;

 

 

 

 

 

N

 

— дисперсия случайной величины

ау =

1

V [К(к>]2 — (F)2

 

VN

КId= 1

 

выходной функции.

 

 

 

Если Туч близко или равно 0, то корреляционная связь между

У и xt слаба, т. е. такие факторы надо исключать из дальнейшего анализа.

Однако коэффициент корреляции rYx, не может являться надеж­ ной характеристикой тесноты связи между величинами У и Х(, если

нет совместного нормального распределения этих величин. Если совместное распределение не является нормальным, то из равен­ ства Гуч = 0 ещё не следует стохастической независимости перемен­

ных Y и xt. Поскольку коэффициент гУх{ является линейным коэф­ фициентом корреляции, это означает, что при rYxi = 0 отсутствует

линейная корреляционная связь, но при этом может существовать нелинейная корреляционная связь. При наличии же нелинейной корреляционной связи гУщ даёт лишь ориентировочное представле­

ние о её тесноте. В то же время гУч позволяет характеризовать сте­

пень приближения корреляционной зависимости к функциональ­ ной, прямолинейной зависимости. Кроме того, за счет влияния не­ учтенных факторов может возникнуть ложная корреляция и даже при | гп ^. | = 1 нельзя окончательно считать, что xt и Xj — зависимые

величины. Поэтому при последующем логико-априорном анализе особое внимание следует обратить на те факторы, для которых rxl%j

близки к 1, а Гуч — к 0. Отбрасывать эти факторы можно лишь после

тщательного анализа их связей.

Далее необходимо исследовать метод отбора значимых факторов (существенно влияющих на выходную функцию) и построения ста­ тистической модели зависимости от них в специфических условиях прогнозирования экономического развития отрасли, учитывая, что помимо переменных величин заданных факторов на значение вы­ ходной функции оказывает влияние большое количество не поддаю­ щихся учету случайных факторов. Это случайное влияние необхо­ димо фиксировать в математической модели в качестве некоторой нормально распределенной случайной величины со средним b и дисперсией а2. "Нормальность распределения обосновывается зако­ ном больших чисел. Учет случайной компоненты в математической модели определяется самой природой исследуемых экономических

процессов.

 

функция

описывается выражением:

Исследуемая выходная

 

Y = + 2

a‘xi ~Ь 2

.S

aijxixj

(3.4.1),

 

i = 1

1=1 / —1

 

 

где *i, х2,

..., хп — переменные факторы, определяющие

значение

а0,

выходной функции

Y;

 

ач — численные значения

коэффициентов;

 

 

Д —нормально распределенная по закону

 

 

Р{-Х)

aVto

С

2а‘

(3.4.2)

случайная величина со средним b и дисперсией о2. Значение b можно без ограничения общности выражения (3.4.1) прибавить к величине коэффициента а0 и тем самым свести пара­ метры А к (0, о2). Среднее значение Д будем считать равным нулю.

57

Пусть имеется таблица наблюдений значений исследуемой вы­ ходной функции:

(3.4.3)

y ( N ) x ( N ) ' Д.(М )> _ _ _ ^ X ( N )

Пусть число значений исследуемой выходной функции будет больше числа факторов, влияющих на неё, т. е. N > п. Индекс снизу означает номер параметра, а индекс сверху — номер года исследуемого периода. Необходимо получить приближение мате­

матического описания зависимости У от факторов xt (i = 1, п). Коэффициенты а0, а;, а,-у и значение а2неизвестны и подлежат оценке.

Запишем значения Y в следующем виде:

пп п

У(1) = а0 - 1+ 2 о м 1 +2 2 aijx‘i"x'i' + Лъ

пп п

у<ло = flo • 1 + 2 а,л:\N) + 2 2 al]Xf ^ + Л„.

Обозначим:

Тогда предыдущие уравнения примут вид:

П2

У(1) = а01 + 2 aiX? + 2 Ькгк'" + дь

 

 

П

2

+Ajy.

 

У(Л/) = й0•1+2 aix\N) + 2

Введем обозначения:

 

_

_

единичный N-мерный

вектор — 1;

вектор

N наблюдений

значений

У — К;

вектор

N наблюдений значений

х<(/=ТГ~п) —*г;

58

Значения случайной величины А могут быть обозначены Димер­ ным вектором А.

Окончательно получим записанное в векторной форме выраже­ ние:

п*+ п

 

 

 

п

 

 

2

 

 

 

 

 

 

У = «0 1 + 2

a-iX-i + 2

 

+ А.

 

 

(3-4-4)

 

 

 

г=1

 

 

к- 1

 

 

 

 

Обозначим через

 

 

 

 

 

 

 

 

 

 

 

йо

 

 

при

/ =

();

 

 

 

 

с ,=

а/ (t = 1, п)

при

1 = \ , П;

 

 

 

 

 

!

 

п2-\-п\

 

 

 

 

 

, /л 2+ л\

 

Ьк [к ==1,

2 )

при

/ =

П-\- 1 "+(

I

)•

 

 

 

1 = 0,

»+№)'

 

 

 

 

 

 

п

при

/ =

0 ;

 

 

 

 

 

 

 

н

Xi

при

/ =

1,

 

 

^/г2+ л \

 

 

 

ZK

при

/ = п + 1,

« +

 

 

 

2

) ’

 

 

 

 

/ = 0 , л + (- 2+ я )■

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У= с , . т +

 

2

 

< А + £

 

 

(3.4.5)

Обозначим:

 

 

 

 

 

 

 

 

 

 

 

Получим

3п + п‘

3п + п»

 

У = Со+ 2

Cfli ) 1 2 Ciyl + A,

(3.4.6)

/*= 1

г- i

 

где уг представляет собой A-мерный вектор, полученный вычитанием из всех значений координат вектора уг величины уг, т. е.

7г = Уг —УгЛ

(Т, у/) = 0 для всех

59

Соседние файлы в папке книги из ГПНТБ