Электроэнергетика была самой значимой и ведущей отраслью советского хозяйства по плану ГОЭЛРО. Поэтому важно понять, что при жёстком ручном управлении в её развитии всё же были колебания. Подробная экспертиза динамики электроэнергетики нашей страны с 1913 года дана в [5–8]. Эта методика [4] применена и при выявлении закономерностей инновационной деятельности в промышленности России. Такой анализ важен и по другим отраслевым направлениям при разработке и реализации национального проекта научно-технологического развития «Сколково», где одним из научных направлений в деятельности станет «анализ трендов технологического развития» [6].
В статье [7] была рассмотрена структура потребления электроэнергии в динамике от 1913 до 1987 г. по отраслям хозяйства царской России и бывшего СССР. В другой статье [5] нами показаны закономерности распределения электроэнергии между отраслями хозяйства. Это позволит в будущем России избежать перекосов в энергетической политике и планах диверсификации всей технологической базы экономики. Масштабная технологическая перестройка в 20-х и 30-х годах ХХ века была обусловлена бурным развитием электроэнергетики. Современная Россия пытается восстановить промышленный потенциал до уровня СССР, но имеет малую по времени предысторию, поэтому в статье [8] был предложен анализ поведения электроэнергетики на инерционное будущее при двух динамических рядах с 1913 года.
В данной статье на примерах из [1] покажем проведение научно-технической экспертизы технологического развития России с 1995 года. По данным [2, 3] на период 1992–2010 гг. будут отдельные статьи. Но сразу же оговоримся, что надежные прогнозные модели можно получать только по динамическим рядам за столетний период развития нашей страны.
Динамика основных показателей инновационной деятельности по статистическим данным из [1] показана в таблице.
Инновационная активность предприятий (добывающие, обрабатывающие производства, производство и распределение электроэнергии, газа и воды)
Год учета |
Время t, лет |
Доля из общего числа организаций, осуществлявшая технологические инновации α, % |
Доля из организаций с техноинновациями, осуществлявшая отдельные виды инноватики αm |
1995 |
3 |
5,5 |
57,9 |
1997 |
5 |
4,7 |
4,2 |
1998 |
6 |
5,0 |
54,3 |
1999 |
7 |
6,2 |
35,1 |
2000 |
8 |
10,6 |
42,7 |
2001 |
9 |
9,6 |
35,5 |
2002 |
10 |
9,8 |
33,6 |
2003 |
11 |
10,3 |
30,2 |
2004 |
12 |
10,5 |
33,0 |
2005 |
13 |
9,3 |
31,6 |
2006 |
14 |
9,4 |
33,4 |
2007 |
15 |
9,4 |
33,5 |
2008 |
16 |
9,6 |
33,2 |
2009 |
17 |
9,4 |
34,6 |
2010 |
18 |
9,3 |
35,7 |
Начало отсчета принимается по данным Росстата [3] t = 0 1992 год.
Без подробных выкладок приведем результаты идентификации устойчивых закономерностей с волновыми составляющими. Из них в дальнейшем можно будет составлять прогнозные модели.
Доля из общего числа организаций, осуществлявшая технологические инновации α, %, определяется (рис. 1) трендом
(1)
Коэффициент корреляции 0,9344 гораздо выше известного предела 0,7. Абсолютная погрешность (или остаток) формулы (1) дана на рис. 1. Сопоставление остатков с данными из таблице показало погрешность моделирования в некоторых точках более 20 %.
Технология моделирования проста: формула, вместе с исходными данными, запускается в программную среду Curve Expert для идентификации связей между количественно измеренными факторами. Совместно с компьютером исследователь выполняет полуавтоматический поиск значений параметров исходно заданной общей модели. Такой процесс поиска значений параметров модели называется – структурно-параметрическая идентификация.
На рис. 2 показан волновой сигнал, полученный по остаткам.
После объединения тренда с волной получена (рис. 3) трехчленная формула вида
α = α1 + α2 + α3, (2)
В формуле (2) и далее приняты следующие условные обозначения: A – амплитуда (половина) колебательного возмущения; p – половина периода колебания системы, лет.
Коэффициент корреляции стал равным единице только из-за малого количества членов статистического ряда. Программная среда способна обрабатывать ряд, содержащий до 3500 членов. Примем концепцию моделирования динамических рядов до максимального выявления волновых составляющих колебательного возмущения системы в процессе функционирования. Следующая волна возмущения показана на рис. 3.
После объединения всех четырех составляющих была получена математическая функция
α = α1 + α2 + α3 + α4; (3)
Тренд по формуле (1) Остатки после тренда (1)
Рис. 1. Динамика доли организаций, выполнявших технологические инновации, из общего количества организаций, по тренду (1): S – дисперсия; r – коэффициент корреляции как мера адекватности модели
Вейвлет-сигнал по остаткам от (1) Тренд и волна по формуле (2)
Рис. 2. Динамика доли организаций, выполнявших технологические инновации
Вейвлет-сигнал по остаткам от (2) График по модели (3) динамики
Волна по формуле (4) Волна по формуле (5)
Рис. 3. Волновая динамика доли организаций с технологическими инновациями
Но возможности программной среды не дают сконструировать общую модель с числом параметров более 19, поэтому вне модели (3) остатки дали еще две волны (рис. 3).
Поэтому нужна специальная программная среда по нашим сценариям, которая позволяла бы одновременно считать по модели со 100 составляющими и не менее 1000 параметрами модели. Причем поиск параметров и идентификации полной модели на простых компьютерах является слишком трудоемким. В итоге нужен суперкомпьютер.
Следующую волну возмущения (рис. 3) нашли вне модели (3) по вейвлет-функции вида
(4)
Даже при малой амплитуде колебания в будущем волна может разрастись до кризисного влияния. Поэтому нужен анализ остатков дальше.
После параметрической идентификации вейвлет-сигнала была получена (рис. 3) шестая составляющая вида
(5)
Как и в предыдущих волнах, АЧХ поведения российской экономической системы сложная – поведение характеризуется переменными значениями амплитуды и периода колебания.
Волна возмущения возрастает на будущее время. Это означает, что в российской действительности с 1999 года (левая граница вейвлета) зреет позитивная (положительный знак пред амплитудой) волна роста доли организаций, занимающихся инновационной деятельностью. Но по периоду она сокращается, начиная от 1,34354∙2 = 2,68708 ≈ 2,7 года, то есть частота возмущения экономической системы нарастет и может наступить тремор. Высокочастотные колебания очень опасны именно из-за того, что не успеет руководство с ручным управлением экономикой.
Суперкомпьютер со специальной программной средой по нашим сценариям смог бы «потрясти» все шесть составляющих и тогда не было бы графика остатков, как показано на рис. 4. Из-за нестыковки без встряски составляющих общей статистической модели получается поправочная формула по биотехническому закону (не показана).
Поправочная составляющая к модели (5) Остатки после поправки
Рис. 4. Динамика остатков от общей модели изучаемого показателя
Остатки после поправки также даны на рис. 4.
На прогноз эти остатки дальше не влияют. Эта волна была в прошлом, поэтому она имеет смысл только для историографического анализа.
Таким образом, мы подробно показали технологию идентификации.
Более кратко продолжим дальше анализ данных из таблицы.
Доля из организаций с техноинновациями, осуществлявшая отдельные виды инноватики αm, %, оценивается (рис. 5) трендом вида
(6)
Колебание было в прошлом, и на прогноз оно не влияет. По первой составляющей формулы (6) наблюдается дальнейший спад по закону экспоненциальной гибели.
Критерии оценки. Они необходимы для сравнительного анализа параметров системы (экономической или иной с участием человека, то есть так называемой эргатической системы) и выбора наиболее адекватного из них для последующего проводимого ежегодно итерационного научно-технологического и социально-экономического мониторинга.
![pic_76.wmf](/fs/i/2013/10-11/pic_76_fmt.jpg)
![pic_77.wmf](/fs/i/2013/10-11/pic_77_fmt.jpg)
Тенденция по модели (6) Остатки после модели (6)
Рис. 5. Динамика из организаций, выполнявших технологические инновации, по отдельным видам технологических инноваций, по тренду (6)
Максимальная относительная погрешность Δmax показывает не только адекватность составной модели, но и характеризует качество процесса управления системой. При этом должно соблюдаться условие, чтобы структурно-параметрическая идентификация исходной закономерности была проведена до максимально возможного количества составляющих (в зависимости от точности измеренных или расчетных исходных данных).
Этот критерий одновременно является важнейшим для оценки замены табличной модели типа таблицы на её математические модели.
Остальные три критерия позволяют выяснить картину поведения в прошлом, то есть в статистически учетный период.
Значимость α1 трендовой (среднестатистической по динамике физической величины, принципиально отличающейся от средних арифметических и иных средних величин) закономерности факторной функции характеризует стабильность и устойчивость будущего функционирования системы. Как говорил В.В. Путин: держать тренд.
Этот оценочный показатель детерминированной части (%) факторной функции (может быть более одной составляющей) вычисляется как отношение расчетных значений, в нашем примере первой или первых двух составляющих, к общему расчетному результату, то есть по математическому выражению
a1 = 100Y1/Y, (13)
где Y1 – первая или первые две составляющие тренда; Y – расчетное значение параметра системы по общей модели.
Коэффициент приспособляемости k изучаемой системы к внешней (в данной статье к экономической) среде показывает устойчивость в прошлом по отношению к разным внешним и внутренним возмущениям.
Этот оценочный показатель исчисляется по формуле
(14)
где Y1 – трендовая закономерность; i – номер составляющей общей закономерности.
Коэффициент динамичности Kд поведения системы показывает характеристику стабильности в прошлом и возможную устойчивость к внешним воздействиям в будущем до горизонта инерционного прогноза.
Указанный критерий самоуправления вычисляется по соотношению
(15)
где yтр – тренд из нескольких не волновых составляющих; ε – абсолютная погрешность (остатки), то есть разница между фактически измеренными значениями параметров системы по данным таблицы и рассчитанными значениями по конкретным моделям.
Вывод
Предлагаемая методология исследования при проведении научно-технической экспертизы, по своей функциональной структуре будущего поведения, вполне позволяет оценивать системы мер в известном документе «Долгосрочный прогноз научно-технологического развития Российской Федерации (до 2025 года)», а затем итерационно проводить мониторинг и уточнять ежегодно программу научно-технологического развития и диверсификации технологической базы экономики России.
Рецензенты:
Салихов М.Г., д.т.н., профессор, заведующий кафедрой автомобильных дорог, ФГБОУ ВПО «Поволжский государственный технологический университет», г. Йошкар-Ола;
Колупаев Б.И., д.б.н., профессор кафедры водных ресурсов, ФГБОУ ВПО «Поволжский государственный технологический университет», г. Йошкар-Ола;
Савин К.Н., д.э.н., профессор кафедры «Экономический анализ и качество», ГОУ ВПО «Тамбовский государственный технический университет», г. Тамбов.
Работа поступила в редакцию 08.11.2013.