ФУРЬЕ ИНТЕГРАЛОВ СПОСОБЫ ВЫЧИСЛЕНИЯ.
Интеграл Фурье (и. Ф.) функции
принадлежащей пространствам
или
Пространство абстрактное в функциональном анализе), имеет вид
И.Ф. широко используется в автоматического управления теории, в теории дифф. ур-ний, в спектральном анализе, распознавании образов, при анализе речевых сигналов и т. д. Рассмотрим некоторые способы вычисления
в зависимости от гладкости подынтегральной ф-ции.
1.
- непрерывная ф-ция, приближенно равная нулю вне отрезка
и заданная табл. своих значений в точках
. Тогда
Если
и погрешности, с которыми задается ф-ция в точках
взаимно независимы и распределены равновероятно на отрезке
, то справедлива следующая оценка полной погрешности А (см. Погрешностей вычислений теория) приведенного алгоритма:
где модуль непрерывности
— число двоичных цифр в мантиссах чисел для ЦВМ с плавающей запятой,
- к-во арифм. операций для вычисления
— округления погрешность при вычислении
. Для вычисления
применяют алгоритм быстрого преобразования Фурье. Отличительной особенностью этого алгоритма является высокая эффективность по сравнению с известными до сих пор алгоритмами. Он основывается на возможности вычисления коэфф.
итерационным методом, что приводит к значительной экономии вычислительного (машинного) времени. Напр., при
требуемое время сокращается прибл. в 100 раз. С ростом
преимущество быстрого преобразования. Фурье становится еще более ощутимым. Быстрое преобразование Фурье можно с успехом применять и для вычисления интегралов типа свертки, автокорреляционной функции, двумерного преобразования Фурье и т. д.
2. Увеличение степени гладкости
не ведет к существенному увеличению точности
и поэтому может оказаться более целесообразным следующий способ, который изложен для прибл. вычисления обратного Фурье преобразования. Пусть
достаточно гладкая ф-ция, заданная на оси и обращающаяся в нуль на бесконечности. Сделаем замену
и положим
Аппроксимируем
многочленом интерполяции
где
тогда
Для вычисления коэфф.
целесообразно
-менять алгоритм быстрого преобразования Фурье. При этом нужно учитывать, что если
вещественна, то
вещественны. Положим
Тогда имеет место следующее соотношение:
где
— многочлены Лягерра, ортогональные на полуоси
с весом
Справедлива оценка
, где
погрешность метода в метрике
граница единичной окружности,
величина наилучшего приближения
многочленами вида
метрике Чебышева
Аппроксимация функций равномерная
).
Если
имеет ограниченную
то
Для того, чтобы получить аналогичную оценку в пространстве
введем ф-цию
и построим и
Положим
Тогда
Справедлива оценка абсолютной погрешности метода в метрике
Если
имеет
то
Коэффициенты
действительны, что дает возможность существенно упростить программы вычисления их. В более общем случае, когда в окрестности бесконечности
следует
применить к
. В результате получают:
где
дельта-функция,
— коэфф. интерполяционного многочлена для
В случае, когда подынтегральная ф-ция принадлежит
гильбертовому пространству
комплекснозначных ф-ций
отличающихся более чем на константу,
производные которых (в обобщенном смысле) квадратично интегрируемы, со скалярным произведением
построен функционал
где
, который является оптимальной (в смысле минимизации погрешности метода) аппроксимацией функционала
в пространстве линейных функционалов вида
При вычислении функционала
целесообразно использовать алгоритм быстрого преобразования Фурье.
Лит.: Иванов В. В. Теория приближенных методов и ее применение к численному решению сингулярных интегральных уравнений. К., 1968 [библиогр. с. 281—285]; Бабушка И., Витасек Э., Прагер М. Численные процессы решения дифференциальных уравнений. Пер. с англ. М., 1969 [библиогр. с. 354—358]. В. К. Задирака.