%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轴会不停抖动,得不到想要的驻波波形。(这步我也不知道为什么,有理解的可以告诉我)