Математические модели управления производственно-экономическими системами. Учебное пособие
.pdf
Экономическая значимость параметров α и τ в модели ЛНТ следующая. На первом периоде значения α – 0,68 и 0,70 – соответствуют ранее рассмотренным моделям ЛКД и ФКД. Величина τ имеет очень малое значение порядка 10-3, оно отрицательно, и в соответствии с t-статистикой недостоверно. Поэтому судить о вкладе НТП на основе модели ЛНТ невозможно. Аналогичное заключение относится к модели на всём временном интервале и к варианту несглаженных данных на втором интервале.
Модель для второго интервала в варианте сглаженных данных достоверна по всем параметрам и характеризует отрицательный вклад от использования основного капитала в анализируемом комплексе в 1990–2001 гг. – α < 0. Это негативное обстоятельство частично компенсируется положительной ролью НТП в эти годы – τ > 0.
Так же как и в предыдущих ситуациях, на основе идентифици- рованных в модели ЛТП коэффициентов проанализируем аппрокси- мативные свойства однородной ПФ Кобба – Дугласа с научно-техни- ческим прогрессом в форме модели ОНТ, отвечающей зависимости
(5.14).
Графики моделируемых и реальных траекторий развития отрасли на основе моделей ОНТ представлены на рис. 5.10, 5.11.
6.625 ×10 37000
6000
|
Y q |
|
|
5000 |
||||||||
|
KD q |
|
|
|
||||||||
|
KD1 |
|
|
4000 |
||||||||
|
|
q1 |
|
|||||||||
|
|
|
|
|||||||||
|
KD2 q2 |
3000 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2000 |
1.179 ×10 31000
1965 |
|
1970 |
1975 |
1980 |
1985 |
1990 |
1995 |
2000 |
2005 |
|
1.965 |
×10 |
3 |
|
|
t q , t q , t q1 |
, t q2 |
|
2.001 |
×10 |
3 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
Рис. 5.10. Траектории выпуска продукции (модель ОНТ, фактические данные) (– – реальный выпуск продукции Y (тыс. руб. в сопоставимых ценах);
- - - расчёт по модели ОНТ на всём интервале исследования; расчёты по ОНТ для первого и второго интервалов)
111
6.632 ×10 37000 |
|
|
|
|
|
|
|
|
|
|
||
|
|
6000 |
|
|
|
|
|
|
|
|
|
|
Y q |
|
5000 |
|
|
|
|
|
|
|
|
|
|
KD q |
|
|
|
|
|
|
|
|
|
|
|
|
KD1 |
|
4000 |
|
|
|
|
|
|
|
|
|
|
q1 |
|
|
|
|
|
|
|
|
|
|
|
|
KD2 q2 |
3000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
|
|
1.233 ×10 31000 |
|
|
|
|
|
|
|
|
|
|
||
|
|
1965 |
|
1970 |
1975 |
1980 |
1985 |
1990 |
1995 |
2000 |
2005 |
|
|
|
1.965 |
×10 |
3 |
|
|
t q , t q , t q1 |
, t q2 |
|
2.001 |
×10 |
3 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
||||
Рис. 5.11. Траектории выпуска продукции (модель ОНТ, сглаженные данные) (– – реальный выпуск продукции Y (тыс. руб. в сопоставимых ценах);
- - - расчёт по модели ОНТ на всём интервале исследования; расчёты по ОНТ для первого и второго интервалов)
Из рис. 5.10, 5.11 видно хорошее согласование расчётных и фактических значений. Имеет место весьма высокое соответствие результатов для сглаженных данных первого периода. Для второго периода, как и в предыдущих случаях, построенная модель требует существенной доработки.
В табл. 5.7 представлены параметры данных ПФ и выборочные, качественные характеристики модели ОНТ.
Таблица 5.7
Параметры и характеристики моделей однородной ПФ Кобба – Дугласа с учетом НТП (модели ОНТ)
Период: 1965–2001 гг.
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
A |
1,3953 |
1,3829 |
|
|
|
α |
0,7569 |
0,7557 |
|
|
|
τ |
-0,0142 |
-0,0136 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
R2 |
0,9191 |
0,9348 |
F |
193,0815 |
243,8519 |
|
|
|
DW |
0,6123 |
0,2799 |
|
|
|
112
|
|
|
|
Окончание табл. 5.7 |
|
|
|
|
|
Период: 1965–1990 гг. |
|
|
|
|
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
|
Сглаженные данные: |
|
|
|
|
|
|
A |
|
1,1830 |
1,2189 |
|
|
|
|
|
|
α |
|
0,6791 |
0,7071 |
|
|
|
|
|
|
τ |
|
-1,6042×10-3 |
-3,3729×10-3 |
|
|
|
|
|
|
Качественные |
характе- |
|
|
|
ристики модели: |
|
|
|
|
|
|
|
|
|
R2 |
|
0,9965 |
|
0,9985 |
|
|
|
|
|
F |
|
3253,9833 |
|
7834,5450 |
|
|
|
|
|
DW |
|
1,3727 |
|
0,7808 |
|
|
|
|
|
Период: 1990–2001 гг. |
|
|
|
|
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
|
Сглаженные данные: |
|
|
|
|
|
|
A |
|
0,4274 |
|
0,2311 |
|
|
|
|
|
α |
|
-0,4387 |
|
-0,9148 |
|
|
|
|
|
τ |
|
0,0638 |
|
0,0999 |
|
|
|
|
|
Качественные |
характе- |
|
|
|
ристики модели: |
|
|
|
|
|
|
|
|
|
R2 |
|
0,4689 |
|
0,7017 |
|
|
|
|
|
R2sc |
|
0,4158 |
|
0,6719 |
|
|
|
|
|
F |
|
3,9727 |
|
10,5850 |
|
|
|
|
|
DW |
|
1,1246 |
|
0,7312 |
|
|
|
|
|
Из данных табл. 5.7 следует, что наиболее адекватной модель ОНТ является на периоде с 1965 по 1990 гг. в варианте сглаженных данных, где коэффициент детерминации возрос до значения R2=0,9985. Для второго периода значения скорректированного коэффициента детерминации R2sc в модели ОНТ формально уменьшились по сравнению с моделью ЛНТ, но при этом одновременно снизился уровень достоверности значений этих величин – F-критерий.
Сохранилась общая тенденция улучшения аппроксимативных качеств моделей с использованием сглаженных данных, при одновременном уменьшении прогнозных качеств моделей и снижении в 1,5–2 раза значений статистики DW, что свидетельствует
113
о росте степени автокорреляции остатков в моделях и неустранении этого явления введением фактора НТП. В целом, модели с учётом научно-технического прогресса – ЛНТ и ОНТ – несмотря на некоторое улучшение показателей качества, являются недостаточно пригодными для адекватного моделирования протекания процессов в исследуемом комплексе.
Попытаемся усовершенствовать модельные конструкции и рассмотрим неоднородную ПФ Кобба – Дугласа с учётом влияния НТП. Запишем модель в виде:
Y (t) = A × K (t)α × L(t) β × eτ×t , α + β ¹ 1. |
(5.16) |
Приведём её к линейной зависимости в логарифмической форме:
ln( )= ln( |
) + α × ln( |
|
)+ β × ln( |
|
) + τ × |
|
. |
(5.17) |
Y |
A |
K |
|
L |
|
t |
|
Четыре неизвестных параметра уравнения (2.17) – модели ФНТ – отыщем методом МНК и проведём статистический анализ качества полученных моделей. Результаты расчётов представлены в табл. 5.8.
Таблица 5.8
Параметры и характеристики линеаризованных моделей неоднородной ПФ Кобба – Дугласа с учетом НТП (модели ФНТ)
Период: 1965–2001 гг.
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
ln (А) |
-2,8765 |
-2,7920 |
|
|
|
α |
0,4740 |
0,4804 |
|
|
|
β |
0,8762 |
0,8601 |
|
|
|
λ (α+β) |
1,3502 |
1,3405 |
|
|
|
τ |
0,0085 |
0,0082 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
σ |
0,0785 |
0,0633 |
|
|
|
DA |
0,5152 |
0,3348 |
|
|
|
Da |
0,0124 |
9,7075×10-3 |
Db |
0,0274 |
0,0210 |
Dt |
7,3020×10-5 |
5,5980×10-5 |
SA |
0,7178 |
0,5786 |
|
|
|
114
|
|
Продолжение табл. 5.8 |
|
|
|
Sα |
0,1113 |
0,0985 |
Sβ |
0,1655 |
0,1450 |
Sτ |
8,5452×10-3 |
7,4820×10-3 |
tA |
4,0074 (tтабл,= 2,04) |
4,8252 (tтабл,= 2,04) |
tα |
4,2577 (tтабл,= 2,04) |
4,8756 (tтабл,= 2,04) |
tβ |
5,2940 (tтабл,= 2,04) |
6,4469 (tтабл,= 2,04) |
tτ |
0,9987 (tтабл,= 2,04) |
1,1012 (tтабл,= 2,04) |
R2 |
0,9782 |
0,9855 |
|
|
|
F |
494,0575 |
747,2597 |
|
|
|
DW |
0,9825 |
0,3701 |
|
|
|
Период: 1965–1990 гг. |
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
ln (А) |
-1,3351 |
-1,7686 |
|
|
|
α |
0,5940 |
0,5964 |
|
|
|
β |
0,5807 |
0,6333 |
|
|
|
λ (α+β) |
1,1747 |
1,2297 |
|
|
|
τ |
0,0029 |
0,0016 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
σ |
0,0351 |
0,0227 |
|
|
|
DA |
1,0235 |
0,6702 |
|
|
|
Dα |
4,8525×10-3 |
3,1699×10-3 |
Dβ |
0,0311 |
0,0211 |
Dτ |
1,9698×10-5 |
1,1195×10-5 |
SA |
1,0117 |
0,8186 |
|
|
|
Sα |
0,0697 |
0,0563 |
Sβ |
0,1765 |
0,1451 |
Sτ |
4,4382×10-3 |
3,3459×10-3 |
tA |
1,3197 (tтабл,= 2,07) |
2,7014 (tтабл,= 2,07) |
tα |
8,5266 (tтабл,= 2,07) |
10,5935 (tтабл,= 2,07) |
tβ |
3,2907 (tтабл,= 2,07) |
4,3644 (tтабл,= 2,07) |
tτ |
0,6571(tтабл,= 2,07) |
0,4695 (tтабл,= 2,07) |
R2 |
0,9966 |
0,9985 |
F |
2146,7231 |
4990,6490 |
|
|
|
DW |
1,1745 |
0,7057 |
|
|
|
115
|
|
Окончание табл. 5.8 |
|
|
|
Период: 1990–2001 гг. |
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
ln (А) |
84,8244 |
50,1088 |
|
|
|
α |
-11,3012 |
-7,2611 |
|
|
|
β |
1,9892 |
2,0611 |
|
|
|
λ (α+β) |
-9,3120 |
-5,200 |
|
|
|
τ |
0,3177 |
0,2442 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
σ |
0,0865 |
0,0542 |
|
|
|
DA |
841,2737 |
261,8357 |
|
|
|
Dα |
11,5730 |
4,0485 |
Dβ |
0,2073 |
0,0932 |
Dτ |
6,7254×10-3 |
2,4729×10-3 |
SA |
29,0047 |
16,1813 |
|
|
|
Sα |
3,4019 |
2,0121 |
Sβ |
0,4553 |
0,3052 |
Sτ |
0,0820 |
0,0497 |
tA |
2,9245 (tтабл,= 2,31) |
3,0967 (tтабл,= 2,31) |
tα |
3,3220 (tтабл,= 2,31) |
3,6087 (tтабл,= 2,31) |
tβ |
4,3689 (tтабл,= 2,31) |
6,7528 (tтабл,= 2,31) |
tτ |
3,8744(tтабл,= 2,31) |
4,9104 (tтабл,= 2,31) |
R2 |
0,7552 |
0,8662 |
|
|
|
F |
0,6634 |
0,8160 |
|
|
|
DW |
8,2251 |
17,2634 |
|
|
|
DA |
2,6167 |
2,3209 |
|
|
|
Данные табл. 5.8 показывают, что для общего и первого периодов расчёты на основе неоднородной ПФ с учётом НТП (модель ФНТ) близки к результатам расчётов по неоднородной ПФ без учёта НТП (модель ФКД) и несколько лучше расчётов по однородной ПФ с учётом НТП (модель ЛНТ).
116
Значения среднеквадратической ошибки СКВО, коэффициентов детерминации R2 для всех вариантов фактических и сглаженных данных зачастую отличаются на доли процентов, а то и полностью совпадают. Так R2 для общего периода равно: 0,9776 и 0,9850 для ФКД; 0,9605 и 0,9701 для ЛКТ и 0,9782 и 0,9865 для ФНТ. Для первого периода соответствующие значения равны: 0,9965 и 0,9985
для ФКД; 0,9924 и 0,9967 для ЛНТ и 0,9966 и 0,9985 для ФНТ.
Практическое совпадение расчётов по модели ФНТ и ФКД можно объяснить тем, что параметр НТП – τ – в модели ФНТ также как и в модели ЛНТ на общем и первом интервалах оказывается незначимым по t-критерию (табл. 5.8, 5.6), и его учёт не приводит к совершенствованию моделей ПФ. Низкое значение критерия DW дополнительно указывает на существование важных факторов, неучтённых в модели ФНТ на общем и первом интервалах.
Совершенно иная ситуация имеет место при использовании модели ФНТ для описания процессов на втором интервале. Здесь показатели качества резко улучшаются. Значения СКВО уменьшаются с 13,83 % и 10,23 % для ФКД и с 13,18 % и 7,58 % для ЛНТ до 8,65 % и 5,42 % для ФНТ. Значения коэффициента детерминации R2 возрастают с 0,2959 и 0,4630 для ФКД и с 0,5228 и 0,8195 для ЛНТ до 0,7552 и 0,8662 для ФНТ. Соответственно, значения скорректированного коэффициента детерминации R2sc возрастают с 0,1395 и 0,3437 для ФКД и с 0,4751 и 0,8015 для ЛНТ до
0,6634 и 0,8160 для ФНТ.
Такое резкое улучшение показателей качества модели ФНТ на интервале 1990–2001 гг. соответствует тому, что параметр НТП – τ – впервые оказывается значимым: расчётные значения t-статистики для него существенно превышают табличные (табл. 5.8). На этом втором интервале оказались значимыми и все остальные параметры модели ФКТ. Эти оценки значимы на уровне α=0,01 для фактических данных и при α=0,001 в случае сглаженных данных.
117
По F-критерию достоверными являются значения R2. Величина критерия Дарбина – Уотсона DW свидетельствует о практическом отсутствии автокорреляции остатков, и о том, что все учтённые в модели факторы – A, α, β, τ – достаточно полно описывают существо протекающих процессов. Таким образом, впервые из всех рассмот- ренных ранее видов ПФ, модель ФНТ оказалась конструктивной при описании процессов на этапе структурных перестроек 1990–2001 гг.
При этом лучшее совпадение расчётов с исходными данными имеется для усреднённых значений. Соответствующую этому случаю величину СКВО σ=0,0542 можно признать удовлетворительной при описании процессов на втором интервале.
Значения параметров α, β и λ в модели ФНТ аналогичны соответствующим величинам модели ФКД (табл. 5.4), и они имеют практически ту же содержательную экономическую сущность, кото- рая была выявлена ранее при анализе модели ФКД. Значения пара- метров τ на общем и первом интервалах являются недостоверными, и содержательной интерпретации не имеют. Параметры α, β и τ на втором интервале, являются достоверными и также, как и в модели ЛНТ, они характеризуют, в целом, отрицательную эффективность использования капитала и трудовых ресурсов, и положительную, компенсирующую роль НТП.
Результаты расчётов с идентифицированными, приведёнными в табл. 5.8 параметрами иллюстрируются графиками на рис. 5.12, 5.13, описывающими выпуск продукции неоднородной ПФ Кобба – Дугласа с учётом НТП в форме зависимости (5.16) – модель ННТ.
На графиках наглядно видно достаточно высокая степень совпадения между модельными и реальными данными. Наиболее чётко согласованность прослеживается для варианта сглаженных исходных данных – рис. 5.13.
118
6.823×1038000 |
|
|
|
|
|
|
|
|
|
6000 |
|
|
|
|
|
|
|
|
|
Yq |
|
|
|
|
|
|
|
|
|
KDq |
|
|
|
|
|
|
|
|
|
4000 |
|
|
|
|
|
|
|
|
|
KD1q1 |
|
|
|
|
|
|
|
|
|
KD2q2 |
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
|
1.179×103 |
0 |
|
|
|
|
|
|
|
|
|
1965 |
1970 |
1975 |
1980 |
1985 |
1990 |
1995 |
2000 |
2005 |
|
1.965×103 |
|
|
|
tq , tq , tq1 ,tq2 |
|
|
|
2.001×103 |
Рис. 5.12. Траектории выпуска продукции (модель ННТ, фактические данные) (– – реальный выпуск продукции Y (тыс. руб. в сопоставимых ценах);
- - - расчёт по модели ННТ на всём интервале исследования; расчёты по ННТ для первого и второго интервалов)
6.7×1037000 |
|
|
|
|
|
|
|
|
|
|
6000 |
|
|
|
|
|
|
|
|
Yq |
5000 |
|
|
|
|
|
|
|
|
KDq |
|
|
|
|
|
|
|
|
|
|
4000 |
|
|
|
|
|
|
|
|
KD1q1 |
|
|
|
|
|
|
|
|
|
KD2q2 |
|
|
|
|
|
|
|
|
|
|
3000 |
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
1.233×10 |
3 |
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
|
|
1965 |
1970 |
1975 |
1980 |
1985 |
1990 |
1995 |
2000 |
2005 |
|
1.965×103 |
|
|
|
tq , tq , tq1 ,tq2 |
|
|
|
2.001×103 |
Рис. 5.13. Траектории выпуска продукции (модель ННТ, сглаженные данные) (– – реальный выпуск продукции Y (тыс. руб. в сопоставимых ценах);
- - - расчёт по модели ННТ на всём интервале исследования; расчёты по ННТ для первого и второго интервалов)
Параметры и основные качественные характеристики модели ННТ представлены в табл. 5.9.
119
Таблица 5.9
Параметры и характеристики моделей неоднородной ПФ Кобба – Дугласа с учетом НТП (модели ННТ)
Период: 1965–2001 гг.
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
A |
0,0563 |
0,0613 |
|
|
|
α |
0,4740 |
0,4804 |
|
|
|
β |
0,8762 |
0,8601 |
|
|
|
λ (α+β) |
1,3502 |
1,3405 |
|
|
|
τ |
0,0085 |
0,0082 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
R2 |
0,9483 |
0,9621 |
|
|
|
F |
201,9241 |
279,1595 |
|
|
|
DW |
0,9119 |
0,3605 |
|
|
|
Период: 1965–1990 гг. |
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
A |
0,2631 |
0,1706 |
|
|
|
α |
0,5940 |
0,5964 |
|
|
|
β |
0,5807 |
0,6333 |
|
|
|
λ (α+β) |
1,1747 |
1,2297 |
|
|
|
τ |
0,0029 |
0,0016 |
|
|
|
Качественные характе- |
|
|
ристики модели: |
|
|
|
|
|
R2 |
0,9953 |
0,9980 |
|
|
|
F |
1560,4193 |
3743,9141 |
|
|
|
DW |
1,2237 |
0,6059 |
|
|
|
Период: 1990–2001 гг. |
|
|
|
|
|
Параметры модели: |
Несглаженные данные: |
Сглаженные данные: |
|
|
|
A |
6,8986×10+36 |
5,7805×10+21 |
|
|
|
α |
-11,3012 |
-7,2611 |
|
|
|
β |
1,9892 |
2,0611 |
|
|
|
λ (α+β) |
-9,3120 |
-5,200 |
|
|
|
τ |
0,3177 |
0,2442 |
|
|
|
120
