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

摘要

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