规则区域上Helmholtz方程的一种快速算法
A Fast Algorithm for Solving Helmholtz Equation on Regular Domain
-
摘要: 采用有限差分法对Helmholtz方程进行五点差分离散,在规则区域上引入快速傅里叶变换(FFT),将差分方程变换成一组三对角方程,使求解规则子区域上Helmholtz方程的计算量降为O(PlgP),最后的数值结果证明了文中的算法是一种快速算法Abstract: The complicated structure can be divided into several almost regular domains using decomposition domain method (DDM). In this paper, finite difference method(FD)is used to discretize the Helmholtz equation applied with five-points difference discretization. On the regular domain, the fast Fourier transform(FFT)is presented to transform the difference equations to a group of trigonal equations. The theory analysis and numerical results show that the computation complex of this algorithm is about O(PlgP).