利用matlab实现驻波仿真

%stationary wave
A=6;k=2;w=0.05;
x=0:0.01:15;
u_0=pi/8;

for t=0:250
    E_1=A*cos(k*x+w*t)
    E_2=A*cos(k*x-w*t+u_0)
    E=E_1+E_2
    plot(x,E_1,'r')
    hold on
    plot(x,E_2,'c')
    hold on
    plot(x,E,'k')
    hold off
    ylim([-15,15])
    title('驻波合成')
    drawnow
end

ylim处要限制振幅，否则y轴会不停抖动，得不到想要的驻波波形。（这步我也不知道为什么，有理解的可以告诉我）