Автоматизация геофизических исследований в высоких широтах.

4 . Соотношение между длиной ряда во временной области и разрешением по частоте л f. Если во временной области зада­ но N . узлов, отстоящих друг от друга через Т , то расстояние между узлами в частотной области д £ = 1 /N T . ( 6 ) 5. Выполняется только для вещественных функций. Спектр ве­ щественных функций имеет комплексно-сопряженную симметрию: :(0 * / \ х ( - f ) , (7 ) т.е. мнимая часть спектра - нечетная функция, а вещественная часть четная функция. Поэтому достаточно для вещественных функ­ ций вычислять лишь половину спектра на положительных частотах. Пакет программ, реализующий БПФ при вычислениях с плава­ ющей точкой,написан на языке ФОРТРАН. Для этого использовались тексты, приведенные в формуле ( 3 ) . Назначение модулей дано в табл . 2 . Таблица 2 Пакет программ БПФ с плавающей точкой Модуль Функции модуля S R F T 3 Вычисление спектра на положительных частотах для действительной функции S R F T 1 БПФ для двух действительных функций S R F T S Умножение функции на фактор. i/U S R F T E БПФ для комплексной функции S R F T 4 Вычисление ОБПФ по половине спектра на положи­ тельных частотах S R F T 2 Получение двух действительных функций из БПФ S R F 'T 'l Двоичная инверсия отсчетов дискретной функции Для выполнения БПФ от комплексной функции необходимо вызвать два модуля: CALL S R F T S ( A , R ) ; CALL S R F T E ( A , R ) , Обратное быстрое преобразование Фурье (ОБПФ) для комплек­ сной функции: CALL S R F T 1 ( A ,R ) ; CALL S R F T E ( A , R ) . 59