Смысл уравнения состоит в том, что без сезонных колебаний все уровни были бы равны среднемесячному, т. е. у¯; колебания же в равной мере разнесены на sin t и cos t. В первом квадранте (т. е. от января до апреля) косинус является положительной величиной и снижается от 1 до 0, синус тоже положителен и возрастает от 0 до 1. Во втором квадранте (апрель — июль) косинус отрицателен и снижается от 0 до -1, синус положителен и снижается от 1 до 0. В третьем квадранте (июль — октябрь) косинус отрицателен, но возрастает от -1 до 0, а синус снижается от 0 до -1. В четвертом квадранте косинус возрастает от 0 до 1 (к декабрю до 0,866), а синус возрастает от -1 до 0 (к декабрю до -0,5). Цикл завершается новым январем. За счет комбинации изменений косинуса и синуса при разных значениях параметров b1 и Ь2 удается отобразить, как показывает табл. 6.4 (графа y^i), любое синусоидальное колебание уровней временного ряда. Имеем: b1 = -484/6 = -80,7; Ь2 = 284/6 = 47,3. Уравнение сезонных колебаний продуктивности коров имеет вид:
y^i = 310 — 80,7∙cos ti + 47,3 sin ti,
где ti = 0° в январе, а месяц = 30° дуги.
Отклонения фактических уровней (но усредненных за ряд лет) от расчетных по ряду Фурье очень малы: максимальное отклонение 7, среднее (по модулю) 3,33, что составляет лишь 1,07 %. Такая точность вполне достаточна для прогнозов и других расчетов. Если же отклонения оказались значительными, следует на основании ряда отклонений повторить расчет, т. е. рассчитать вторую гармонику, и тогда окончательные уровни модели (ряда Фурье) будут представлять собой сумму всех гармоник:
где т — число гармоник;
к — номер гармоники.
Однако если колебания явно не имеют синусоидальной формы, то требуется много гармоник, расчет становится трудоемким и гораздо проще применить метод, описанный в разд. 6.3.2.
6.4. Измерение тренда колеблемости
Неоднократно указывалось на большое значение мониторинга колебаний. Как правило, производство, экономика заинтересованы в уменьшении колеблемости. Чтобы измерить изменение абсолютного показателя силы колебаний S(t), проще всего рассчитать эту величину за последовательные отрезки времени, а затем по полученным значениям S(t)1, S(t)2 и т. д. до S(t)n провести аналитическое выравнивание, т. е. вычислить тренд того или другого типа. Однако для более надежного вычисления меры колеблемости необходимо как минимум 7–9 уровней первичного временного ряда, а для вычисления тренда по этим мерам колеблемости — опять 7–9 таких же частных мер S(t). А для этого первичный ряд должен содержать примерно 8 x 8 = 64 уровня. Такие ряды анализируются нечасто, а значит, пет и условий для расчета тренда мер колеблемости.
Положение отчасти спасает то, что для вычисления тренда колеблемости вовсе необязательно, чтобы за весь изучаемый период существовал единый тренд уровней показателя. Вполне допустимо для расчета тренда колеблемости объединить отрезки времени с разными по типу трендами или с кусочно-линейным трендом. От изменения скорости роста или даже типа роста, или направления тенденции динамики колеблемость зависит мало или совсем не зависит. Но и с учетом этой ее особенности измерить тренд колеблемости по ряду отдельных отрезков времени сложно. При длине первичного ряда в 15–20 уровней получается всего два значения S(t), чего явно не хватает для расчета тренда.
Не вполне корректными с математической точки зрения являются расчет скользящих показателей колеблемости со сдвигом в один период времени и последующее их аналитическое выравнивание. Конечно, скользящие показатели уже зависят друг от друга, но выявить общую тенденцию изменения силы колебаний и приближенно измерить тренд S(t) все же возможно. Покажем применение этого метода на примере временного ряда урожайности зерновых культур во Франции (см. разд. 5.1). В приложении 1 вычислены отклонения уровней от тренда, с которых и начинается измерение тренда среднего квадратического отклонения (табл. 6.5).
Скользящие показатели колеблемости S(t)i будем рассчитывать по 11-летним подпериодам, т. е. первый за 1970–1980 гг., второй ~ за 1971–1981 гг. и т. д. Первая величина S(t) будет относиться к середине подпериода, т. е. 1975 г. и т. д., последняя скользящая средняя за 1985–1995 гг. относится к 1990 г. Итого получаем 16 скользящих значений показателей колеблемости, которые и выравниваем по уравнению прямой.