danbai.m

来自「天津大学matlab讲义MATLAB是一种用于算法开发,数据可视化,数据分析以及」· M 代码 · 共 30 行

%横梁
plot([-.2;0.2],[0;0],'color','y','linestyle','-','linewidth',10);
%初始位置的单摆
g=9.8;
l=1;
theta0=pi/6;x0=l*sin(theta0);
y0=-l*cos(theta0);
axis([-0.75,0.75,-1.25,0]);
axis('off');
%摆锤
head=line(x0,y0,'color','r','marker','.',... 
 'erasemode','xor','markersize',80);
%摆杆
body=line([0;x0],[0;y0],'color','b','linestyle','-', ...
   'erasemode','xor','linewidth',2.5);
%摆动
t=0;
dt=0.005;
while 1
   t=t+dt;
   theta=theta0*cos(sqrt(g/l)*t);
   x=l*sin(theta);
   y=-l*cos(theta);
   set(head,'xdata',x,'ydata',y);
   set(body,'xdata',[0;x],'ydata',[0;y]);
   drawnow;
end