Математические модели управления производственно-экономическими системами. Учебное пособие
.pdf
|
|
Продолжение табл. 5.4 |
|
|
|
tA |
4,4959 (tтабл,= 2,04) |
5,2748 (tтабл,= 2,04) |
tα |
35,8499 (tтабл,= 2,04) |
44,1654 (tтабл,= 2,04) |
tβ |
11,1470 (tтабл,= 2,04) |
13,2453 (tтабл,= 2,04) |
R2 |
0,9776 |
0,9850 |
|
|
|
F |
740,6457 |
1113,3183 |
|
|
|
DW |
0,9733 |
0,3901 |
|
|
|
Период: 1965–1990 гг. |
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
ln (А) |
-0,9780 |
-1,5190 |
|
|
|
α |
0,6366 |
0,6214 |
|
|
|
β |
0,5002 |
0,5810 |
|
|
|
λ (α+β) |
1,1368 |
1,2024 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
σ |
0,0347 |
0,0223 |
|
|
|
DA |
0,7102 |
0,3744 |
|
|
|
Dα |
6,2446×10-4 |
3,2720×10-4 |
Dβ |
0,0157 |
8,3349×10-3 |
SA |
0,8428 |
0,6119 |
|
|
|
Sα |
0,025 |
0,0181 |
Sβ |
0,1254 |
0,0913 |
tA |
1,1605 (tтабл,= 2,07) |
2,4826 (tтабл,= 2,07) |
tα |
25,4755 (tтабл,= 2,07) |
34,3534 (tтабл,= 2,07) |
tβ |
3,9879 (tтабл,= 2,07) |
6,3636 (tтабл,= 2,07) |
R2 |
0,9965 |
0,9985 |
F |
3301,4448 |
7748,4921 |
|
|
|
DW |
1,2014 |
0,7426 |
|
|
|
Период: 1990–2001 гг. |
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
ln (А) |
-11,9155 |
-21,1617 |
|
|
|
α |
1,4452 |
2,2018 |
|
|
|
β |
0,9435 |
1,2484 |
|
|
|
λ (α+β) |
2,3887 |
3,4502 |
|
|
|
101
|
|
Окончание табл. 5.4 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
σ |
0,1383 |
0,1023 |
|
|
|
DA |
284,4554 |
182,042 |
Dα |
1,918 |
1,1876 |
Dβ |
0,3437 |
0,2347 |
SA |
16,8658 |
13,4923 |
|
|
|
Sα |
1,3849 |
1,0898 |
Sβ |
0,5863 |
0,4845 |
tA |
0,7065 (tтабл,= 2,26) |
1,5684 (tтабл,= 2,26) |
tα |
1,0435 (tтабл,= 2,26) |
2,0204 (tтабл,= 2,26) |
tβ |
1,6093 (tтабл,= 2,26) |
2,5766 (tтабл,= 2,26) |
R2 |
0,2959 |
0,4630 |
|
|
|
F |
0,1395 |
0,3437 |
|
|
|
DA |
1,8914 |
3,8797 |
DW |
0,9921 |
0,4507 |
|
|
|
Из анализа данных табл. 5.4 следует, что на общем периоде идентификации все оценки коэффициентов построенных моделей значимы. На первом периоде незначимыми являются оценки для масштабного коэффициента ln(А). На втором интервале исследования значимой является только оценка параметра β в варианте сглаженных данных. Наблюдается общая тенденция роста расчётных значений t-статистики для моделей, рассчитанных по сглаженным данным.
Аппроксимативные свойства моделей ФКД улучшились по сравнению с моделями ЛКД. Значения среднеквадратической ошибки стали меньше и составили 7,85 % и 6,35 % для общего и 3,47 % и 2,23 % для первого периодов. На втором периоде оценки СКВО практические не изменились.
Коэффициент детерминации для общего и первого периодов увеличился, наилучшее значение данный показатель принимает для модели с 1965 по 1990 гг. в варианте сглаженных данных – R2=0.9985. Эти коэффициенты являются значимыми, F-статистика значительно превышает табличные значения.
102
Модели на интервале с 1990 по 2001 гг. стали иметь меньшие значения скорректированного коэффициента детерминации R2sc. Однако эти показатели нельзя признать значимыми, так как расчётные значения F-статистики составляют 1,89 и 3,88, при табличном значении F2,9=4.2 (при уровне значимости α=0,05).
Устойчиво проявляется тенденция уменьшения ошибки аппроксимации и повышения адекватности при использовании сглаженного варианта первичной статистики на всех интервалах исследования. Значения критерия Дарбина – Уотсона DW изменились незначительно и остаются достаточно низкими для признания отсутствия автокорреляции.
Для содержательной экономической интерпретации модели ФКД помимо факторных эластичностей α и β в табл. 5.4 приведены также значения полной эластичности λ=α+β, характеризующей интегра- льную, масштабную эффективность производства, в целом.
На первом периоде 1965–1990 гг. значения эластичности по фондам α равны 0,64 и 0,62. Эти величины близки к значениям, идентифицированным в модели ЛКД однородной ПФ Кобба – Дугласа (см. табл. 2.1). Эластичность выпуска по труду равна 0,50 и 0,58. Это достаточно высокие оценки вклада трудовых ресурсов. При этом α>β, что отражает больший вклад капитала в конечные продукт. Суммарная эластичность выпуска, в целом, λ равна 1,14 и 1,20. Эти значения характеризуют достаточно высокую суммарную эффективность использования капитальных и трудовых ресурсов в первом периоде.
Во втором периоде значения α чрезмерно высоки 1,45 и 2,20. Однако в соответствии с t-статистикой их следует признать недостоверными, и содержательного экономического смысла они не имеют. Аналогично обстоят дела и со значениями β и λ на этом интервале.
На всём периоде моделирования значения α – 0,58 и 0,59 практически те же, что и у однородной ПФ – модели ЛКД (табл. 5.1),
103
и они достаточно достоверны. Значения β равны 0,72 и 0,71. Это очень высокие значения вклада трудовых ресурсов, существенно превышающие вклад капитала β>α. Соответственно, весьма высокой оказывается и суммарная эластичность λ – 1,31 и 1,30. Значения β и λ в модели ФКД с экономической точки зрения являются завы- шенными, что с математической точки зрения является следствием влияния периода с 1990 по 2001 гг., где модель ФКД оказывается недостоверной по всем параметрам и показателям качества.
В целом, модель ФКД также как модель ЛКД не учитывает важные факторы, определяющие сущность экономических процессов, что приводит к неудовлетворительным характеристикам этих моделей, кроме, частично, первого периода с 1965 по 1990 гг.
Проанализируем качество описания исследуемых процессов на основе модели неоднородной ПФ Кобба – Дугласа – НКД в форме (5.12). Численные значения параметров модели возьмём по результатам идентификации, представленным в табл. 5.4.
На рис. 5.8, 5.9 представлены графики, иллюстрирующие пове- дение расчётных траекторий моделей НКД.
6.625×1037000 |
|
|
|
|
|
|
|
|
|
|
6000 |
|
|
|
|
|
|
|
|
Yq |
5000 |
|
|
|
|
|
|
|
|
KDq |
|
|
|
|
|
|
|
|
|
KD1 |
4000 |
|
|
|
|
|
|
|
|
q1 |
|
|
|
|
|
|
|
|
|
KD2q2 |
|
|
|
|
|
|
|
|
|
|
3000 |
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
1.179×10 |
1000 |
|
|
|
|
|
|
|
|
|
1965 |
1970 |
1975 |
1980 |
1985 |
1990 |
1995 |
2000 |
2005 |
|
3 |
|
|
|
tq,tq,tq1 |
,tq2 |
|
|
3 |
|
1.965×10 |
|
|
|
|
|
2.001×10 |
||
|
|
|
|
|
|
|
|
||
Рис. 5.8. Траектории выпуска продукции (модель НКД, фактические данные) (– – реальный выпуск продукции Y (тыс. руб. в сопоставимых ценах);
- - - расчёт по модели НКД на всём интервале исследования; расчёты по НКД для первого и второго интервалов)
104
6.587×1037000 |
|
|
|
|
|
|
|
|
|
|
6000 |
|
|
|
|
|
|
|
|
Yq |
5000 |
|
|
|
|
|
|
|
|
KDq |
|
|
|
|
|
|
|
|
|
KD1 |
4000 |
|
|
|
|
|
|
|
|
q1 |
|
|
|
|
|
|
|
|
|
KD2q2 |
|
|
|
|
|
|
|
|
|
|
3000 |
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
1.233×10 |
1000 |
|
|
|
|
|
|
|
|
|
1965 |
1970 |
1975 |
1980 |
1985 |
1990 |
1995 |
2000 |
2005 |
|
1.965×103 |
|
|
|
tq,tq,tq1,tq2 |
|
|
|
2.001×103 |
Рис. 5.9. Траектории выпуска продукции (модель НКД, сглаженные данные) (– – реальный выпуск продукции Y (тыс. руб. в сопоставимых ценах);
- - - расчёт по модели НКД на всём интервале исследования; расчёты по НКД для первого и второго интервалов)
Визуально наблюдается удовлетворительное совпадение расчётных и фактических значений характеристик для первого временного интервала. Особенно хорошее соответствие имеет место для сглаженных данных. На втором интервале сходимость расчётов с фактическими значениями остаётся низкой.
Значения параметров и основных качественных характеристик моделей НКД представлены в табл. 5.5.
Таблица 5.5
Параметры и характеристики моделей неоднородной ПФ Кобба – Дугласа (модели НКД)
Период: 1965–2001 гг.
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
A |
0,0910 |
0,0976 |
|
|
|
α |
0,5839 |
0,5879 |
|
|
|
β |
0,7241 |
0,7117 |
|
|
|
λ (α+β) |
1,3080 |
1,2996 |
|
|
|
105
|
|
|
Окончание табл. 5.5 |
|
|
|
|
Качественные |
характе- |
|
|
ристики модели: |
|
|
|
|
|
|
|
R2 |
|
0,9468 |
0,9603 |
F |
|
302,3056 |
411,6964 |
|
|
|
|
DW |
|
0,8837 |
0,3666 |
|
|
|
|
Период: 1965–1990 гг. |
|
|
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
|
|
A |
|
0,3761 |
0,2189 |
|
|
|
|
α |
|
0,6366 |
0,6214 |
|
|
|
|
β |
|
0,5002 |
0,5810 |
|
|
|
|
λ (α+β) |
|
1,1368 |
1,2024 |
|
|
|
|
Качественные |
характе- |
|
|
ристики модели: |
|
|
|
|
|
|
|
R2 |
|
0,9958 |
0,9982 |
|
|
|
|
F |
|
2759,3471 |
6531,9131 |
|
|
|
|
DW |
|
1,3201 |
0,6632 |
|
|
|
|
Период: 1990–2001 гг. |
|
|
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
|
|
A |
|
0,0000067 |
6,4502×10-10 |
|
|
|
|
α |
|
1,4452 |
2,2018 |
|
|
|
|
β |
|
0,9435 |
1,2484 |
|
|
|
|
λ (α+β) |
|
2,3887 |
3,4502 |
|
|
|
|
Качественные |
характе- |
|
|
ристики модели: |
|
|
|
|
|
|
|
R2 |
|
0,3136 |
0,4600 |
|
|
|
|
R2sc |
|
0,1610 |
0,3400 |
|
|
|
|
F |
|
2,0557 |
3,8339 |
|
|
|
|
DW |
|
0,8851 |
0,4359 |
|
|
|
|
Из сопоставления характеристик неоднородных производствен- ных функций – ФКД (табл. 5.4) и НКД (табл. 5.5) следует, что, в целом, аппроксимативные свойства модели НКД несколько ниже, чем модели ФКД. Уменьшилась величина коэффициентов детерминации R2 и R2sc для всех моделей, кроме варианта несглаженных данных для второго периода. Снизились расчётные значения F-статистик (кроме
106
этой же ситуации). Остались неудовлетворительными величины критериев Дарбина – Уотсона DW для всех моделей.
В целом модели, построенные на основе неоднородных ПФ для общего и первого интервалов, имеют несколько лучшие аппрокси- мативные свойства и адекватность по сравнению с аналогичными моделями однородных ПФ, однако учёт в моделях значимых факторов неудовлетворителен. В случае моделирования поведения процессов на интервале с 1990 по 2001 гг. модели продолжают оставаться непригодными для адекватного описания протекающих процессов.
5.3. СИНТЕЗ МОДЕЛЕЙ С УЧЁТОМ НАУЧНО-ТЕХНИЧЕСКОГО ПРОГРЕССА
Как было отмечено ранее, одним из возможных существенных факторов, неучтённых в предшествующих построениях является научно-технический прогресс. Введём данный фактор в построенные
впредыдущих разделах модели, и проведём сравнительный анализ качественных характеристик таких моделей.
Втеории производственных функций фактор НТП обычно вводится как дополнительный член уравнения, экспоненциально зависящий от времени t в виде (4.51). Будем применять модели ПФ с нейтральным НТП, – нематериализованным ни в одном из рассматри- ваемых факторов производства [89].
Математическую модель в форме ОКД с учётом НТП запишем
ввиде
Y (t ) = A × K (t )α × L(t )1-α × eτ×t , |
(5.14) |
где τ – характеристика темпа роста эффективности технологических
процессов вследствие воздействия НТП [55]. Приведя зависимость (5.14) к виду, аналогичному
линейную форму модели:
Y |
|
K |
|
|||
ln |
|
|
= ln(A) + α × ln |
|
|
+ τ × t . |
|
|
|||||
|
L |
|
L |
|
||
(5.13), получим
(5.15)
107
Методом наименьших квадратов найдём неизвестные параметры: A, α, τ, и проведём качественный анализ полученных характеристик модели. Результаты расчётов для линейной формы (5.15) пред- ставлены в табл. 5.6.
Таблица 5.6
Параметры и характеристики линеаризованных моделей на основе однородной ПФ Кобба – Дугласа с учетом научно-технического прогресса (модели ЛНТ)
Период: 1965–2001 гг.
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
ln (А) |
0,3331 |
0,3242 |
|
|
|
α |
0,7569 |
0,7557 |
|
|
|
τ |
-0,0142 |
-0,0136 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
σ |
0,0972 |
0,0842 |
|
|
|
DA |
0,0217 |
0,0181 |
|
|
|
Dα |
0,0142 |
0,0122 |
Dτ |
8,0520×10-5 |
6,7919×10-5 |
SA |
0,1474 |
0,1347 |
|
|
|
Sα |
0,1192 |
0,1104 |
Sτ |
8,9733×10-3 |
8,2413×10-3 |
tA |
2,2605 (tтабл,= 2,04) |
2,4067 (tтабл,= 2,04) |
tα |
6,3500 (tтабл,= 2,04) |
6,8485 (tтабл,= 2,04) |
tτ |
1,5823 (tтабл,= 2,04) |
1,6558 (tтабл,= 2,04) |
R2 |
0,9605 |
0,9701 |
|
|
|
F |
413,8423 |
550,6695 |
|
|
|
DW |
0,9275 |
0,2930 |
|
|
|
Период: 1965–1990 гг. |
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
ln (А) |
0,1681 |
0,1979 |
|
|
|
α |
0,6791 |
0,7071 |
|
|
|
τ |
-1,6042×10-3 |
-3,3729×10-3 |
|
|
|
108
|
|
Окончание табл. 5.6 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
σ |
0,0378 |
0,0250 |
|
|
|
DA |
4,3428×10-3 |
2,1445×10-3 |
Dα |
2,5594×10-3 |
1,3011×10-3 |
Dτ |
1,7349×10-5 |
8,6010×10-6 |
SA |
0,0659 |
0,0463 |
|
|
|
Sα |
0,0506 |
0,0361 |
Sτ |
4,1652×10-3 |
2,9328×10-3 |
tA |
2,5501 (tтабл,= 2,07) |
4,2741 (tтабл,= 2,07) |
tα |
13,4228 (tтабл,= 2,07) |
19,6021 (tтабл,= 2,07) |
tτ |
0,3852 (tтабл,= 2,07) |
1,1501 (tтабл,= 2,07) |
R2 |
0,9924 |
0,9967 |
|
|
|
F |
1507,9664 |
3442,3455 |
|
|
|
DW |
1,0794 |
0,7219 |
|
|
|
Период: 1990–2001 гг. |
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
ln (А) |
-0,85003 |
-1,4647 |
|
|
|
α |
-0,4387 |
-0,9148 |
|
|
|
τ |
0,0638 |
0,0999 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
σ |
0,1318 |
0,0758 |
|
|
|
DA |
0,7470 |
0,2276 |
Dα |
0,5747 |
0,1741 |
Dτ |
2,7765×10-3 |
8,3744×10-4 |
SA |
0,8643 |
0,4770 |
Sα |
0,7581 |
0,4173 |
Sτ |
0,0527 |
0,0289 |
tA |
0,9835 (tтабл,= 2,26) |
3,0706 (tтабл,= 2,26) |
tα |
0,5787 (tтабл,= 2,26) |
2,1922 (tтабл,= 2,26) |
tτ |
1,2116 (tтабл,= 2,26) |
3,4546 (tтабл,= 2,26) |
R2 |
0,5228 |
0,8195 |
R2sc |
0,4751 |
0,8015 |
F |
4,9308 |
20,4316 |
|
|
|
DW |
1,4964 |
0,7433 |
|
|
|
109
Проанализируем данные табл. 5.6 и сравним их с характеристи- ками линеаризованной однородной ПФ Кобба – Дугласа без учёта НТП – ЛКД (5.3).
Отметим снижение оценок значимости параметров однородной ПФ – A и α – в модели ЛНТ (табл. 5.6) по сравнению с моделью ЛКД (табл. 5.1), кроме одного значения для сглаженных данных второго периода. Введённый в модель ЛНТ фактор НТП – τ – значим лишь в модели на интервале с 1990 года в варианте сглаженных данных. Во всех других моделях фактор НТП является незначимым. Значения среднеквадратических ошибок в моделях ЛНТ и ЛКД близки между собой.
Для моделей второго временного интервала характерно значительное увеличение скорректированного коэффициента детерминации R2sc с 0,3895 до 0,4751 по фактической статистике и с 0,5381 до 0,8015 при сглаженных данных.
Все значения R2 являются значимыми на уровне α=0,001, расчётные значения F-критерия в несколько раз превосходят табличные. Исключение составляет второй интервал, где уравнение регрессии (5.15) значимо только на уровне α=0,05 для исходных данных.
Учёт в моделях НТП фактора времени τ, лишь незначительно повысил оценки DW. В моделях для первого и второго интервалов значения критерия DW превышают единицу в вариантах несгла- женных данных, однако и в этих случаях нельзя признать верной гипотезу об отсутствии автокорреляции остатков.
Как и во всех предыдущих моделях, более удовлетворительные оценки имеют модели ЛНТ на интервале с 1965 по 1990 гг. Эти модели имеют высокие аппроксимативные свойства (коэффициент детерминации R2 возрастает до 0,9967), однако статистическая незначимость параметра τ и наличие автокорреляции не позволяют в полной мере удовлетворительно оценить эти модели.
110
