基于MATLAB的简正模正常谱峰分裂计算

基于MINEOS软件、MATLAB的ode45工具箱和样条插值函数interp1，设计计算简正模正常谱峰分裂的接口。采用静力水准重力势能J2、J4系数计算地球椭率一阶、二阶项，考虑PREM和1066A两种起始模型，计算0S2模态的正常谱峰分裂参数和单谱频率，并与Dahlen的结果比较，证明了MATLAB接口的正确性。研究发现，椭率二阶项对简正模谱峰分裂参数影响较小，对单谱频率影响可忽略。最后以₀S₂为例，计算其耦合矩阵，评估耦合强度，为谱峰分裂函数反演地球内部密度结构奠定基础。

Computation of Normal Splitting of Normal Modes Based on MATLAB

In this paper, an interface is built to compute the normal modes splitting parameters, based on the MINEOS software, ode45 toolbox of MATLAB, and interpolation function. The hydrochloric coefficient of J2 and J4 are used to compute the first and second order term of Earth’s ellipticity. The splitting parameters and singlets frequencies of 0S2 are generated based on the PREM and 1066A. It is proven to be correct when compared with the results of Dahlen. Our research shows that the second order ellipticity has a little effect on the splitting parameters, but almost no influence on the singlet frequencies. Finally, the coupling matrix and coupling strength between ₀S₂ and selected modes are given. This lays a good base to use the splitting function to inverse the interior Earth structure.