A practical and simple stable method for calculating Fourier integrals is proposed, effective both at low and at high frequencies. An approach based on the fruitful idea of Levin, to use of the collocation method to approximate the slowly oscillating part of the antiderivative of the desired integral, allows reducing the calculation of the integral of a highly oscillating function (with a linear phase) to solving a system of linear algebraic equations with a three-diagonal triangular or five-diagonal band Hermitian matrix. The choice of Gauss-Lobatto grid nodes as collocation points makes it possible to use the properties of discrete "orthogonality" of Chebyshev differentiation matrices in physical and spectral spaces. This is realized by increasing the efficiency of the numerical algorithm for solving the problem. The system pre-conditioning procedure leads to significantly less cumbersome and more economical calculation formulas. To avoid possible numerical instability of the algorithm, we proceed to the solution of a normal system of linear algebraic equations.
翻译:在低频和高频下,建议采用一个实用、简单、稳定的方法来计算Fourier综合体。根据Levin的富有成果的想法,采用合用法,以近似理想综合体抗分体缓慢振动部分,从而降低高度振动功能(线性阶段)的集成计算,以用三对角三角或五对角带Hermitian矩阵解决线性代数方程系统。选择Gaus-Lobatto网格节点作为合用点,使得在物理和光谱空间使用Chebyshev差异矩阵的离散性“正方位性”特性成为可能。这可以通过提高解决问题的数字算法效率来实现。系统预调程序使计算公式明显不那么麻烦、更经济。为避免算法可能的数字不稳定,我们着手解决普通线性平方程式系统。