1、为介电常数，为真空中的磁导率简记为,则正演模型的离散形式为b5E2RGbCAP=初值条件离散为 。 %时间步长dx = 0.25*10（-2 %空间步长dy = 0.25*10（-2nx = 200。 %网格剖分ny = 300。% wavefield initial value for two timeu1 = zeros（nx,ny %u1初始化u2 = u1。u3 = u1。c = u1。c（：,: = 2*108。%波速 sigma =0.1*10（-9% subwave, ricker waveletf = 5000*106。 %频率s,tw = ricker（f,dt%波源方程sn
2、 = length（s%s的长度T = 350。%迭代次数for k=1:Tfor i=2:nx-1for j=2:ny-1if（kss = （s（k+1else ss = 0。endu3（i,j=（c（i,j*dt/dx2*（u2（i+1,j-2*u2（i,j+u2（i-1,j+.p1EanqFDPw +（c（i,j*dt/dy2*（u2（i,j+1+u2（i,j-1+.DXDiTa9E3d ss*（i=50&（j=ny/2+. +2*u2（i,j-u1（i,j c（i,j2/2*dt*sigma*u1（i,j/（1+c（i,j2/2*dt*sigma；RTCrpUDGiT u3（1,: =
3、0。%边界条件为0 u3（nx,: = 0。 u3（:,1,ny up（k,: u1 = u2。%迭代赋值 u2 = u3。 imagesc（u3 %画出u3pause（0.01 %暂停0.01s矩形子波接收的信号编码： %time grid %space grid% c = 1000。 % wave velocity %grid numberc（1:50,:c（50:70,1:10070,100:200 = 3*108。70,200:300c（70:200,:up=zeros（T,nyf = 7000*106。 %frequencysigma =0.1*10（-9 k ss = （s（k+1 u3（i,j+.5PCzVD7HxA+.jLBHrnAILg ss*（i=2（j=j xHAQX74J0X =u3（2,: = u3（nx-1,: =u3（:,2,ny-1=u3（2,: imagesc（up pause（0.01圆形子波代码c（:f = 4000*106。sigma =0.1*10（-8T = 300。 ss = s（k+1+.LDAYtRyKfE+.Zzz6ZB2LtkdvzfvkwMI1申明：所有资料为本人收集整理，仅限个人学习使用，勿做商业用途。
升级会员 