8.2. Автономная подпрограмма spid1

Данная подпрограмма написала на алгоритмическом языке ФОРТРАН-4 и предназначается для проведения операции вычисления нормированной спектральной плотности одномерного /скалярного/ СП по его дискретизированной корреляционной функции.

Обращение к подпрограмме имеет вид:

CALL SPID1 (R0, NKORF, DTAU, S , FI, NSPEK).

Здесь:

R0 (NКОRF) - массив значений статистической нормированной кор­реляционной функции , соответствующих равномерному разбиению области определения с шагом ;

NKORF- число узлов дискретизации нормированной корреляционной функции СП;

DTAU - шаг разбиения области определения функции ,c;

S(NSРЕК) - массив оценок значений функции нормированной спе­ктральной плотности , построенных на заданной шкале аргумента этой функции;

FI(NSPEK) - шкала аргумента функции спектральной плотности, Гц;

NSPEK - число узлов дискретизации функции нормированной спектральной плотности.

Операция вычисления оценки значений функции нормированной спе­ктральной плотности в подпрограмме производится по формуле:

, .

8.3. Автономная подпрограмма spid2

Данная подпрограмма написана на алгоритмическом языке ФОРТРАН-4 и предназначается для проведения операции вычисления оценки нор­мированных авто-и взаимных функций спектральных плотностей двух одномерных /скалярных/ СП по их дискретизированным корреляционным функциям.

Обращение к подпрограмме имеет вид:

CALL SPID2 (R0X, R0У, R0ХУ, R0XY, NКОRF, DTAU,

SX, SY, MODSXY, ARGSXY, FI, NSPEK) .

Здесь:

R0X, R0У(NКОRF)- массивы значений статистических нормирован­ных автокорреляционных функций , СП Х(Т) и У(Т), со­ответствующих равномерному разбиению областей их определения с ша­гом ;

R0ХУ, R0YX(NKORF)- массивы значений статистических нормирован­ных взаимных корреляционных функции , СП Х (Т) и УСТ), соответствующих равномерному разбиению областей их определения с шагом ;

NKORF - число узлов дискретизации нормированных авто-и взаимных корреляционных функций;

DTAU - шаг разбиения областей определения авто-и взаимных корреляционных функций, с;

SX, SУ (NSPEK)- массивы значений статистических нормированных функций спектральной плотности , СП X (T) и Y (Т), построенных на заданной шкале частоты ;

MODSXY, ARGSXY(NSPEK) массивы значений модуля и углового аргумента /в радианах/ статистической функции нормированной взаимной спектральной плотности СП Х (Т) и У (Т ) , построенных на заданной шкале частоты ;

FI(NSPEK) - шкала частоты , Гц;

NSPEK - число узлов дискретизации функций спектральной плотности.

Операция вычисления оценки значений функций нормированной спектральной плотности в подпрограмме производится по формулам:

8.4. Автономная подпгограмма spidm

Данная подпрограмма написана на алгоритмическом языке ФОРТРАН-4 и предназначается для проведения операции вычисления оценки нор­мированных авто-и взаимных функций спектральных плотностей составляющих многомерного /векторного/ СП по их дискретизированным корреляционным функциям.

Обращение к подпрограмме имеет вид:

CALL SPIDM (R0, NKORF, M, DTAU, MODS, ARGS, FI, NSPEK) .

Здесь:

R0 (NKORF,M,M) - матрица значений авто-и взаимных корреляционных функции составляющих СП Х (Т) с номерами "J" и "к" соответствующих равномерному разбиению области их определения с шагом ;

NKORF - число узлов дискретизации корреляционной функции СП;

М - составляющих векторного СП Х(Т);

DTAU - шаг разбиения области определения корреляционной функции, с;

MODS (NSPEK,M,M) - матрица значений модуля нормированных статистических авто-и взаимных функций спектральной плотности составляющих СП Х (Т) с номерами "J" и "к" ;

соответствующих заданной шкале частоты ;

ARGS (NSPEK,M,M) - матрица значений углового аргумента статистических функций спектральной плотности составляющих СП Х(Т) с номерами "J" и "к" , соответствующих заданной шкале частоты , рад;

FI (NSPEK) - шкала частоты , Гц;

NSРЕК - число , узлов дискретизации функции спектральной плотности.

Операция вычисления оценки значений функции нормированной спектральной плотности в подпрограмме производится по формулам:

,

,