一种高效的宽带简正波本征值计算方法  

An efficient method for broadband eigenvalue computation

作　　者：杨雪峰[1,2] 王好忠 骆文于[4] 胡长青[1] YANG Xue-feng;WANG Hao-zhong;LUO Wen-yu;HU Chang-qing(Shanghai Acoustic Laboratory,Institute of Acoustics,Chinese Academy of Sciences,Shanghai 201815,China;University of Chinese Academy of Sciences,Beijing 100049,China;Ocean University of China,O＇ngdao 266100,Shandong,China;Institute of Acoustics,Chinese Academy of Sciences,Beijing 100190,China)

机构地区：[1]中国科学院声学研究所东海研究站,上海201815 [2]中国科学院大学,北京100049 [3]中国海洋大学,山东青岛266100 [4]中国科学院声学研究所,北京100190

出　　处：《声学技术》2018年第3期201-204,共4页Technical Acoustics

摘　　要：为提高宽带简正波本征值的计算效率,汉密尔顿方法通过公式变换抵消频率项,将宽带简正波本征值的计算由频率、本征值实部和虚部的三维寻根降低至本征值实部和虚部的二维寻根,可以一次性求解单号简正波所有频率对应的简正波本征值。在已有工作的基础上优化和完善了汉密尔顿宽带简正波本征值算法,并加入并行计算方法进一步提高计算效率。以简正波模型KRAKENC作为对比,通过若干数值算例验证了该算法对于宽带简正波本征值的计算精度和计算效率。数值仿真结果显示,在保证宽带简正波本征值计算精度的前提下,该方法的计算效率相对KRAKENC有着明显的优势;加入并行算法后,该方法的计算效率得到大幅提高。Based on existing Hamiltonian method for eigenvalue calculation of normal mode, an efficient method for broadband eigenvalue calculation is developed. Hamiltonian method gets rid of the frequency term by formula transformation to change the root finding problem from a three-dimensional problem with respect to frequency as well as real and imaginary parts of eigenvalue into a two-dimensional one. As a result, the eigenvalues of singular normal mode corresponding to all frequencies could be solved at one time. A MATLAB parallel computing method and other optimizations are included to improve the computational efficiency of Hamiltonian method. The significant advantage of Hamiltonian method in efficiency is indicated by several numerical simulation comparisons with KRAKENC, while retaining the same accuracy.

关 键 词：宽带 简正波方法 本征值 汉密尔顿方法 高效 

分 类 号：P733.21[天文地球—物理海洋学]