не могут меняться существенно, и функция ф(£) должна удовлетворять условиям стационарности [8]. Проверка на стационарность может быть проведена, например, по дисперсии с использованием критерия масштаба Кокре- на [23]. Для всех N 2 интервалов (выборок объемом N n ) *2 процесса ф вычисляют оценки дисперсии о г - и статисти ку G = Omax/2X- 2 которая сравнивается с критическим (=1 значением Z 0 (y, N n , P d }, vjifi v =N n — 1. Если G0,9). Затем программа GAMMA выполняет вещественный гармонический анализ временного ряда ф«. Дискретное прямое преобразование Фурье производится с помощью процедуры FFT (метод БПФ), имеющейся в математиче ском обеспечении ЕС ЭВМ на языке ПЛ/1 [22]. Числовой ряд должен при этом содержать 2 м элементов (М — целое число), а для хранения коэффициентов Фурье .требуется 2 м + 2 полей памяти. Вычисляется параметр М и из ряда фі формируется массив А1 с ближайшим к Ао + 2 числом элементов М +2 = 2 ЛІ + 2^А'' 0 + 2, при чем добавленные сверх числа N o элементы заполняют нулями. Обращение к процедуре FFT содержит вектор А1, параметр М и признак ОРТ, определяющий харак тер операции. При ОПТ=0 выполняется вещественный анализ, 1 — комплексный анализ, 2 — вещественный синтез, 3 — комплексный синтез. В результате работы FFT (ОРТ=0) на место входного массива А1 помеща ются коэффициенты Фурье Х 3 - (j=l, 2, 3, ..., /Ѵ ( ) ,/2а о- 6 о = О.«Р * ) ..... a N _ u b N _ x , l/2a N ,b N =0 (^=1/2^). Фиксируется временной интервал Т — (N) — 1)Н с уче «9