matlab二维平面波场绘制

%二维平面简谐波绘制
%参数定义
nx=201;     %x方向网格点数
nz=201;     %z方向网格点数
nt=201;      %z方向网格点数
h=10;       %空间网格大小
dt=0.002;   %时间网格大小
A=100;      %振幅
f=35;       %频率
w=2*pi*f;   %圆频率
v=3000;     %地震波波速
alpha=pi/4; %地震波传播方向与x轴夹角
k=w/v;
% time=0:0.002:1; %t以0.002为步长从0到1
%%计算视波速
kx=k*cos(alpha);%x方向视波速
kz=k*sin(alpha);%z方向视波速
%%8计算波函数
wave=zeros(nz,nx); %初始化二维平面波场
for m=1:nt       %时间网格坐标
    for i=1:nz   %z轴网格坐标
        for j=1:nx %x轴网格坐标
            x=(j-1)*h; %x坐标
            z=(i-1)*h;  %z坐标
            t=(m-1)*dt; %t坐标
             wave(i,j)=A*exp(1i*(kx*x+kz*z-w*t)); %二维平面波场
        end
    end
hold on
imagesc(real(wave));
axis([0 200 0 200]);
xlabel('x/grid');
ylabel('z/grid');
title('二维平面波场');
drawnow;
end

% figure;
% plot(real(wave));
% figure;
% plot(wave);
% figure;
% plot(real(wave(101,:)));