横波的matlab程序

x=0:.2:10; %离散波线
 n=length(x)
 y=zeros(1,n)
 u=1 %波速
 plot(x,y,'or',x,y) %波线的各个质元都在平衡位置
 axis([0,10,-1,1])
 grid on 
 pause %按任意键继续，波线上各个质元从波源开始由近及远依次振动
 
 
 for t=0:0.1:20    %时间离散   
  y=0.5*cos(t-x-pi/2).*(t>x/u);%波动方程，波源的振动状态传到的质元才能振动   
  plot(x,y,'o-.r') %依次画出t时刻的波形图
  xlabel('x轴')
  ylabel('y轴')
  title('横波的动画')
  axis([0,10,-1,1])
  grid on
  pause(0.05)   %暂停0.05秒，视觉暂留
 end