
книги из ГПНТБ / Пивоваров, С. Э. Моделирование процессов прогнозирования в приборостроении
.pdfКоэффициенты а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) = а0• 1 + 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