matlab怎么编写驻波的动画显示代码

%function f=prob5_11(K,ki)
%prob5_10函数产生驻波动画，K控制驻波动画动态变化的次数，ki控制曲线动态变化的快慢
nargin=0
clf
t=0;k=0;
x=(0:100)/100*3*pi;
n=length(x);
x0=[0 pi 2*pi 3*pi];
y0=zeros(1,4);
y=sin(t).*sin(x);
h=line(x,y,'LineWidth',4,'Color',[0 0 1]);
h1=line(x0,y0,'LineWidth',4,'Marker', '.','MarkerSize',50,'Color',[1 0 0]);
axis off;
% set(h,'erasemode','background') %h是需要执行动画图像的句柄，一般都是由line或者plot创建
if nargin==0
K=1;ki=1;
elseif nargin==1
ki=1;
end
while 1
t=t+0.5;
y=sin(t).*sin(x);
axis([ 0 10 -1.1 1.1])
set(h,'xdata',x,'ydata',y) %更新图像的坐标数据
drawnow %刷新屏幕
set(h1,'xdata',x0,'ydata',y0,'erasemode','background')%更新图像的坐标数据
drawnow %刷新屏幕
if round(t/48)==K
break
end
pause(0.1*ki)