Текст подпрограммы и версий ( Фортран )
ftfsr.zip 		Тексты тестовых примеров ( Фортран )
tftfsr.zip
Текст подпрограммы и версий ( Си )
ftfsr_c.zip 		Тексты тестовых примеров ( Си )
tftfsr_c.zip
Текст подпрограммы и версий ( Паскаль )
ftfsr_p.zip 		Тексты тестовых примеров ( Паскаль )
tftfsr_p.zip

Подпрограмма:  FTFSR

Назначение

Вычисление коэффициентов Фурье по системе нечетных тригонометрических функций вещественного ряда длины, равной степени двух, с использованием метода быстрого преобразования Фурье.

Математическое описание

Пусть длина  N преобразуемого вещественного ряда  X является целой степенью двух:  N = 2 M.

Ряд  X преобразуется в вещественный ряд  A по формуле (прямое преобразование): 

                    N-1             
        A r  =    ∑   X n * sin (π r * n / N)  ,    
                    n=0
            r = 0, ..., N - 1

Tак как система функций  sin (π r x),  r = 0, ..., N - 1 ортогональна на сетке  Xn = n / N, но не ортонормирована, обратное преобразование совпадает с прямым с точностью до постоянного множителя  N / 2. B вычислениях используется метод быстрого преобразования Фурье [1].

Ряд коэффициентов Фурье  A получается на месте исходного ряда  Х без дополнительной рабочей памяти.

При больших  N число операций пропорционально  N * log2N.

1.  В.А.Морозов, Н.Н.Кирсанова, А.Ф.Сысоев. Комплекс алгоритмов быстрого преобразования Фурье дискретных рядов. Сб. "Численный анализ на ФОРТРАНе", вып. 15, Изд - во МГУ, 1976 г.

Использование 

    SUBROUTINE  FTFSR (X, N) 
   Параметры
X - вещественный вектоp длины  N; перед началом работы подпрограммы содержит преобразуемый ряд  X. По окончании работы подпрограммы содержит преобразованный ряд  A;
N - заданная длина вектоpа  X (тип: целый). 
   Версии:  нет 
   Вызываемые подпрограммы
FTF1C - вычисление дискретного или обратного дискретного преобразования Фурье комплексного ряда длины, равной степени двух, методом быстрого преообразования Фурье. 
   Замечания по использованию
  1. 

Длина N преобразуемого ряда должна быть целой степенью двух:  N = 2 M,  M - натуральное. B противном случае результат работы подпрограммы не определен.

  2.  Первая компонента преобразованного ряда  A всегда pавна нулю, т.е. по окончании работы подпрограммы  X (1) = 0 .

Пример использования 

      REAL X(8)
      DATA  X /1., 2., 3., 4., 5., 6., 7., 8./
      CALL  FTFSR (X, 8)

Результаты:

      X  =  (0.,  25.136697,  -9.656854,  7.4830288,  -3.999999, 
                3.3408931895,  -1.656854,  0.994561837)