声波有限差分法正演模拟c语言程序.docx

《声波有限差分法正演模拟c语言程序.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、#include#include#definefm30#definedt0.001#definePI3.#defineNt401#defineNx200#defineNz200//---------------加载震源，雷克子波-----------------------------voidfun(floatsource[]){FILE*fp;intit,i;floatt1,t2,t0;for(i=0;i

2、");for(it=0;it

3、,a4,a5,a6;intflag;floatdx,dz;floatUn[Nx][Nz],Unm[Nx][Nz],Unp[Nx][Nz],v[Nx][Nz],record[Nx][Nt];//Un表示n时刻，Unm表示n-1时刻，Unp表示n+1时刻波场floatsource[Nt];nsx=Nx/2;nsz=0.0;//震源位置N=6;order=2*N;//精度if(order==2){a0=-2.00000;a1=1.00000;a2=0;a3=0;a4=0;a5=0;a6=0;}if(order==4){a0=-2.50000;a1=1.33333;a2=-8.33333e-2;a3

4、=0;a4=0;a5=0;a6=0;}if(order==6){a0=-2.72222;a1=1.50000;a2=-1.50000e-1;a3=1.11111e-2;a4=0;a5=0;a6=0;}if(order==8){a0=-2.84722;a1=1.60000;a2=-2.00000e-1;a3=2.53968e-2;a4=-1.78571e-3;a5=0;a6=0;}if(order==10){a0=-2.92722;a1=1.66667;a2=-2.38095e-1;a3=3.96825e-2;a4=-4.96032e-3;a5=3.17460e-4;a6=0;}if(orde

5、r==12){a0=-2.98278;a1=1.71429;a2=-2.67857e-1;a3=5.29101e-2;a4=-8.92857e-3;a5=1.03896e-3;a6=-6.01251e-5;}/*for(ix=0;ix

6、z

7、ix][iz]=2*Un[ix][iz]-Unm[ix][iz]+0.5*v[ix][iz]*v[ix][iz]*dt*dt/(dx*dx)*(a0*Un[ix][iz]+a1*(Un[ix+1][iz]+Un[ix-1][iz])+a2*(Un[ix+2][iz]+Un[ix-2][iz])+a3*(Un[ix+3][iz]+Un[ix-3][iz])+a4*(Un[ix+4][iz]+Un[ix-4][i