有限差分纵波波场模拟(8)

j=i-1;

U2(:,cellnum(2)-j)=U1(:,cellnum(2)-j)+U1(:,cellnum(2)-j-1)-U0(:,cellnum(2)-j-1)-... (V_model(:,cellnum(2)-j)*dt/dx).*((U1(:,cellnum(2)-j)-U1(:,cellnum(2)-j-1))-... (U0(:,cellnum(2)-j-1)-U0(:,cellnum(2)-j-2)));%右边界条件

U2(:,1+j)=U1(:,1+j)+U1(:,2+j)-U0(:,2+j)-(V_model(:,1+j)*dt/dx).*((U1(:,1+j)-... U1(:,2+j))-(U0(:,2+j)-U0(:,3+j)));%左边界条件

U2(cellnum(1)-j,:)=U1(cellnum(1)-j,:)+U1(cellnum(1)-j-1,:)-U0(cellnum(1)-j-1,:)-... (V_model(cellnum(1)-j,:)*dt/dz).*((U1(cellnum(1)-j,:)-U1(cellnum(1)-j-1,:))-... (U0(cellnum(1)-j-1,:)-U0(cellnum(1)-j-2,:)));%下边界条件

U2(1+j,:)=U1(1+j,:)+U1(2+j,:)-U0(2+j,:)-(V_model(1+j,:)*dt/dz).*((U1(1+j,:)-... U1(2+j,:))-(U0(2+j,:)-U0(3+j,:)));%上边界条件 end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% U0=U1; U1=U2; end

colormap(gray) imagesc(U2)

32