瑞雷波在n层介质中传播的频散曲线的程序

附程序：(1)三层介质瑞雷波程序程序思想：1、首先根据凡友华的标量传递算法求出瑞雷波的频散函数，即令x5=0，它是关于波速和频率的非线性方程；2、然后固定频率的数值，使波速从1500到3500取值，给

附程序： (1) 三层介质瑞雷波程序 程序思想： 1x5=0 、首先根据凡友华的标量传递算法求出瑞雷波的频散函数，即令，它是关 于波速和频率的非线性方程； 215003500 、然后固定频率的数值，使波速从到取值，给定某一搜索步长，搜 索有根区间，只要找到一个有根区间，就调用二分法的求根子函数，求得零点； 32 、最后改变频率的数值重复步，得到不同频率时候的零点值。 主函数： [reslut.m] functionresult % 求根的主函数 i=1;% 记录根的个数 n=9.9;%n 为搜索步长 forf=0:10:3000;%03000 频率分别在到变化 forct=1700:n:3500%17003500 波速在到之间变化时搜素有根区间 a=feval('rayleigh',ct,f); b=feval('rayleigh',ct+9.9,f); ifa==0 nn(f/10+1,i)=ct; i=i+1; elseifa*b<0 jie(f/10+1,i)=erfen('rayleigh',ct,ct+9.9,f);% 调用求根函数 i=i+1; p=i-1; else continue end end i=1;