瑞雷波在n层介质中传播的频散曲线的程序附程序:三层介质瑞雷波程序程序思想:1、首先根据凡友华的标量传递算法求出瑞雷波的频散函数，即令x5=0，它是关于波速和频率的非线性方程;2、然后固定频率的数值，使波速从1500到3500取值，给定某一搜索步长，搜索有根区间，只要找到一个有根区间，就调用二分法的求根子函数，求得零点;3、最后改变频率的数值重复2步，得到不同频率时候的零点值。主函数:[reslut.m]functionresult求根的主函数i=1;%
记录
混凝土 养护记录下载土方回填监理旁站记录免费下载集备记录下载集备记录下载集备记录下载
根的个数n=9.9;%n为搜索步长forf=0:10:3000;%频率分别在0到3000变化forct=1700:n:3500%波速在1700到3500之间变化时搜素有根区间a=feval('rayleigh',ct,f);b=feval('rayleigh',ct+9.9,f);ifa==0nn(f/10+1,i)=ct;i=i+1;elseifa*b<0jie(f/10+1,i)=erfen('rayleigh',ct,ct+9.9,f);%i=i+1;p=i-1;elsecontinueendendi=1;调用求根函数endfid=fopen('test.txt','wt');forj=1:(3000/10)+1fprintf(fid,'%3.3e%3.3e',10*(j-1),jie(j,1:p));fprintf(fid,'\n');endfclose(fid);二分法求根子函数:[erfen.m]functionc=erfen(fun,a,b,t)%二分法求解频散方程的子函数%传递参数时fun是一个函数文件%a,b区间两端点%delta计算精度%c所求方程的根delta=0.0001;ya=feval(fun,a,t);yb=feval(fun,b,t);max1=1+round((log(b-a)-log(delta))/log(2));%max1表示最大迭代次数fork=1:max1c=(b+a)/2;yc=feval(fun,c,t);ifyc==0a=c;b=c;elseifyb*yc>0b=c;yb=yc;elsea=c;ya=yc;endifb-a
本文档为【瑞雷波在n层介质中传播的频散曲线的程序】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483