wave1.m

来自「四面体有限元程序,单元刚度的计算及总刚度矩征集成」· M 代码 · 共 10 行

t=0:pi/20:4*pi;   
x=0:0.1*pi:4*pi;
m=moviein(40)
for i=1:40
   y=sin(x-t(i));
   plot(x,y);
   axis([0,12,-1,1]);
   m(:,i)=getframe;
end
movie(m,3,30)