4.6. Применение МГУА в задачах прогнозирования в макроэкономике.
В результате работы была разработана программа Designer, которая позволяет применить НМГУА для произвольной выборки исходных данных. Для задания исходных данных необходимо сформировать файл-проект текстового формата, в котором должны быть заданы сложность модели опорной функции (то есть наибольшая степень полинома), количество переменных, количество узлов интерполяции, матрица исходных данных, свобода выбора и некоторая дополнительная информация. Программа старается построить модель, адекватную исходной выборке, и выводит результаты в виде таблиц и графиков. Пользователь в процессе работы программы имеет возможность изменять свободу выбора на ряде селекции, весовые коэффициенты внешних критериев оптимальности и другие параметры, сохранять результаты работы программы в текстовом и графическом форматах данных.
Исходя из постановки задачи, решение будет компромисным в силу того, что необходимо достичь высокой точности построения модели при наименьшей ее сложности. Поэтому в принципе может существовать несколько таких моделей, которые удовлетворяют заданным требованиям, а выбор наилучшей из них возлагается на эксперта. Программа дает возможность не только находить наилучшую модель по внешним критериям оптимальности, а также выводит данные по всем моделям на последнем ряде селекции.
Рис.1. Модель, построенная с помощью НМГУА. Вертикальной линией разделены обучающая и тестовая последовательности. Нулю соответствует июль 1995 года, а двадцати восьми - декабрь 1997 года. Среднеквадратичное отклонение - 0.51.
На рисунке 1 приведены результаты построения модели с квадратичной опорной функцией. Исходя из полученных данных можно утверждать, что модель имеет достаточно высокое качество, она близка к реальным данным как на обучающей, так и на тестовой последовательности. Причем качественные показатели полученной модели на тестовой последовательности не ниже аналогичных показателей для обучающей последовательности.
В таблице 1 приведен вид модели, которая соответствует рисунку 1. При свертывании этой модели использовалась геделевская нумерация переменных, в которой индекс каждой переменной должен быть простым числом. Обычно полином имеет достаточно большое количество членов, и большинство из них близки к нулю, поэтому для упрощения данной модели можно воспользоваться процедурой "фильтрации" членов, влияние которых незначительно.
Таблица 1. Модель прогноза ИПЦ |
(2,543001*1)+(1,251916*x6*x6)+(0,290404*x6*x12)+(0,001978*x1*x1)+(0,018901*x1*x6)+ (0,026448*x6*x6*x6*x6)+(0,069067*x6*x6*x6*x12)+(0,024024*x6*x6*x12*x12)+(0,007946*x1*x6*x6)+ (0,001631*x1*x12)+(0,000218*x1*x1*x6*x12)+(0,001322*x1*x6*x6*x12)+(0,000113*x1*x12*x12)+ (5,988685*x1*x6*x12*x12)+(5,962365*x1*x1*x1*x6)+(0,086859*x15)+(0,000618*x12*x12*x12)+ (1,981465*x12*x12*x12*x12)+(0,166181*x6*x6*x15)+(0,004761*x6*x15*x15)+(0,021663*x6*x12*x15)+ (0,000408*x12*x15*x15)+(0,002495*x6*x6*x15*x15)+(0,000305*x6*x12*x12*x15)+ (2,750355*x12*x12*x15*x15)+(0,000521*x6*x15*x15*x15)+(0,000189*x1*x6*x6*x6*x6)+ (0,000724*x6*x6*x6*x6*x12*x12)+(7,027005*x6*x6*x12*x12*x12)+(2,663105*x6*x12*x12*x12*x12)+ (0,004047*x6*x6*x6*x6*x15)+(0,008955*x6*x6*x6*x12*x15)+(0,000149*x6*x6*x12*x15*x15)+ (0,000581*x6*x6*x12*x12*x15)+(2,089095*x6*x12*x12*x15*x15)+(1,345485*x12*x12*x12*x15)+ (1,162365*x12*x15*x15*x15)+(4,512145*x6*x6*x6*x6*x6*x6)+(0,000334*x6*x6*x6*x6*x6*x12)+ (0,000376*x6*x6*x6*x12*x12*x12)+(0,000227*x6*x6*x6*x6*x15*x15)+ (0,000213*x6*x6*x6*x12*x15*x15)+(3,678255*x6*x6*x12*x12*x12*x15)+ (3,447085*x6*x6*x12*x15*x15*x15)+(0,000397*x1*x6*x15)+(2,599605*x1*x15*x15)+ (3,223425*x1*x6*x6*x15)+(1,200305*x1*x12*x15)+(3,942445*x1*x6*x6*x12*x15)+(2,053845*x1*x1*x15)+ (3,543775*x1*x6*x6*x6*x12*x12)+(3,847715*x6*x6*x6*x6*x6*x12*x15)+ (7,192535*x6*x6*x6*x6*x12*x12*x15)+(1,280725*x1*x6*x6*x6*x6*x15) |
То есть приведенная в таблице 1 модель прогноза имеет вид:
При построении модели для квадратичной опорной функции мы столкнулись с проблемой выбора весовых коэффициентов внешних критериев, которые обеспечивают построение адекватной модели для показателя ИПЦ. Для различных исходных данных (имеется в виду выбор разного размера тестовой последовательности) мы получали различный допустимый диапазон изменения веса одного из критериев, который обеспечивает получение адекватной модели. Таким образом, допускается существование эффективной областиі весовых коэффициентов, которые обеспечивают адекватное построение модели. Пример такой области для исследуемых экспериментальных данных приведен на рисунке 2. Эффективная область в данном случае лежит между двумя приведенными зависимостями.
Рис.2. Эффективная область выбора весовых коэффициентов внешних критериев оптимальности для квадратичной опорной функции.
В таблице 2 приведены значения среднеквадратичного отклонения для различных размеров тестовой последовательности для моделей разной сложности. Как видно из таблицы, среднеквадратичное отклонение для случая линейных опорных функций в среднем вдвое превышает соответствующие значения для квадратичной опорной функции, но преимуществом линейных опорных функций является значительно большая скорость расчетов, поэтому для грубой оценки целесообразно использовать линейное моделирование. Однако при неэффективном выборе весовых коэффициентов внешних критериев оптимальности для квадратичной опорной функции мы получим результаты в 9 раз хуже, чем при использовании линейных моделей, и в 18 раз хуже, чем при использовании квадратичных моделей в эффективной области. Необходимо также заметить, что существование такой области доказывает необходимость использования нескольких внешних критериев, что значительно повышает качество моделирования.
Таблица 2. Среднеквадратичные ошибки для построенных линейных и квадратичных моделей для заданных экспериментальных данных.
Размер тестовой последова-тельности |
Среднеквадратичное отклонение |
||
Линейные модели |
Квадратичные модели в эффективной области | Квадратичные модели в неэффективной области | |
1 |
0.821 |
0.386 |
0.506 |
2 |
0.463 |
0.410 |
19.427 |
3 |
1.880 |
0.483 |
4.922 |
4 |
0.922 |
0.562 |
1.632 |
5 |
0.842 |
0.526 |
6.862 |
6 |
1.073 |
0.901 |
1.550 |
7 |
1.039 |
0.703 |
3.044 |
8 |
0.984 |
0.661 |
2.847 |
9 |
1.053 |
0.634 |
3.738 |
10 |
1.212 |
0.511 |
52.523 |