一种低网格高精度基于数值与符号混合计算的波方程本征值算法.docVIP

  一种低网格高精度基于数值与符号混合计 算的波方程本征值算法 喻志远* （电子科技大学物理电子学院，成都 610054） 摘要：本文给出了一种解波方程问题的新的算法。文中将数值差分方法与 Fourier 变换相结 合并以数值计算和符号计算相混合来求解波方程的本征值问题。与频域差分法中用迭代法求 解矩阵方程不同，本方法用 Matlab 中的符号计算工具，求解几个网格数不同的矩阵方程。 与时域有限差分方法相比，本方法避免了时间空间步迭代与有限 Fourier 反变换所引入的数 值误差，而保留了时域有限差分方法的优点，即一次计算就可以得到全部频率下的信息。论 文采用了一种级数加速收敛技术，极大地提高了本征值的精度。作为例子，一个一维波方程 的本征值问题和一矩形波导本征值问题用本方法进行了求解。其结果与理论值的对比，充分 说明了本方法的有效性。 关键词：符号与数值计算；波方程；本征值；加速收敛 中图分类号：TN011;TP391.41 An Algorithm of Low Mesh and high Accuracy with Symbolic and Numerical Computation For Eigenvalue Problems Yu Zhiyuan (University of Electronic Science and Technology of China School of Physical Electronics, ChengDu 610054) Abstract: A new method for solving the wave equations is presented in this paper. The numerical difference method is combined with Fourier transformation and a symbol calculation is mixed with numerical calculation to solve the eigenvalue problems of wave equations. Instead of using the iteration method to solve the matrix in finite difference frequency domain, it solve the several matrices using the symbolic computation by a simple Matlab? codes. This eliminates the errors caused by both time space iteration and finite numerical Fourier transformation in common FDTD method. But it retains the strong point of FDTD, that is, one can get all frequency information through only one calculation. A scheme of accelerating the convergence a series is carried out. This greatly enhanced the accuracy of the calculation results Examples of the wave equation in one dimension cavity and the cutoff wave number in a rectangular waveguide are given to demonstrate the feasibility of the new method. Keywords: Symbolic and numerical computational; wave equations; eigenvalue problems; acceleration of Convergence 0 引言 在工程与理论研究中，符号与数值混合运算是一强有力的工具。在 E. L. Kaltofen 的论 文中[1]符号与数值混合运算被认为是“七个小矮人”中的第五个。符号与数值混合运算被用 于电磁兼容，天线问题，以及其他电磁基本问题的求解[2-4]。 符号运算还被用于数学处理中， 如线性多变量系统中的分解，矩阵行列式的求解[5-6]。另外，符号与数值混合运算在控制系 统的设计中还有重要应用[7]。 边值问题在电磁理论和相关工程问题中是一基本问题。人们对此问题进行了大量研究， 也出