Математическое моделирование комплексных процессов / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Вычисл. центр ; [редкол.: В. С. Мингалев (отв. ред.) и др.]. - Апатиты : Кольский филиал АН СССР, 1982.

плазмы был разработан комплекс программ M ERCA N I , по зволяющий проводить обработку каждым из указанных методов. Этот комплекс программ включает в себя расчет элементов траектории ИСЗ, параметров дифракционной картины в топоцентрической и сле дящей системе координат, а также расчет высот, масштабов и флук туаций концентрации электронной плотности Д N в неоднородности. Расчет орбитального движения ИСЗ в гринвичской системе координат осуществляется соответствующей частью программы КАДР, любезно предоставленной нам сотрудниками ИКИ АН СССР / 6 /. В топоцентрической системе, т.е. в прямоугольной системе координат с центром, расположенным в пункте приема радиосигналов (H’qj Л 0 - географические широта и долгота пункта приема), расчет координат и скорости спутника производится с учетом сжатия Земли (формула 1 1 . 0 6 0 в книге / 7 / ) по известным соотношениям, например гла ва П книги / 8 /. По соответствующим формулам рассчитываются Ас - азимут, 8 С - угол места, г с - дальность и другие т^аекторные параметры ИСЗ, необходимые для определения геометрии распростра нения радиоволны. В се эти величины рассчитываются на момент времени, соответствующий середине восьмисекундной зоны регистра ции мерцаний радиосигнала. Входными элементами для расчета траектории ИСЗ являются шесть оскулирующих элементов: а - большая полуось, е - эксцентри ситет, і - наклонение, £1 - долгота восходящего узла, со - аргумент перигея, и - параметр, определяющий положение спут ника на орбите. Исходным материалом для статистического анализа мерцаний являются три временных ряда, представляющие собой три сечения наземной дифракционной картины и полученные во время ее движе ния над приемными антеннами. Первым этапом анализа является нахождение автокорреляционной и трех кросскорреляционных функций между парами антенн; они рассчитываются по стандартным выраже ниям / 9 / - . Л » (1) N - k с . . ( к ) =%7" Г ( х . . - х . ) ( х . . . - X . ) , О < к < L —1 , ( 2 ) Ч N it i /v jt+k j ' 1 N ~ k c . . ( - k ) = — У ( x . . . - x . ) ( x . . - x . ) , 0 < к < L—1 , (3 ) i j v ' N £-> 4 i t +к j jt j ' где N - число точек в массиве выборки временного ряда, а Ь = 0 .1 N . Метод полного корреляционного анализа используется в соот ветствии с общим приближением, данным в работах / 3 , 4 / . При про ведении анализа предполагается, что исследуемый статистический процесс стационарен во времени и однороден в пространстве; прост ранственный контур постоянной корреляции амплитуды может быть аппроксимирован квадратичной формой; вид пространственной корре ляционной функции р(Е,,т?,о') подобен виду временной корреляцион- 135