wave2.m

四面体有限元程序,单元刚度的计算及总刚度矩征集成

%wave2.m
t=0:pi/20:4*pi;
x=0:0.1*pi:4*pi;
h=plot(x,y);
y=sin(x);
axis([0,12,-1,1]);
set(h,'EraseMode','xor')
for i=2:60
   y=sin(x-t(i));
   set(h,'XData',x,'YData',y)
   drawnow;
   pause(0.001)
end