用matlab编写驻波图,Matlab驻波图

clear;clc

A=1;x=50;a=5;b=3;i=1;

t=[0:pi/50:2*pi]/(2*pi)*b;

y1=A*sin(2*pi*(t./b-x/a));

y2=A*sin(2*pi*(t./b+x/a));

% subplot(311);plot(t,y1);

% subplot(312);plot(t,y2);

% subplot(313);plot(t,y1+y2);

b1=[];b2=[];

x10=0;

x20=100;

b=zeros(1,1000);

while 1

y10=y1(i);

y20=y2(i);

b1=[y1(i) b1];

b2=[b2 y2(i)];

n1=length(b1);

b11=[b1 zeros(1,1000-n1)];

n2=length(b2);

b22=[zeros(1,1000-n2) b2];

b=b11+b22;

plot([0:0.1:100-0.1],b,[x10 x20],[b(1) b(1000)],'o')

axis([-1 101 -2.5 2.5])

drawnow

i=i+1;

if i==101

i=1;

end

if length(b1)==1000

b1=b1(1:999);

b2&