单摆模型MATLAB程序.docx

《单摆模型MATLAB程序.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、%Simplependulummodel%Allrightsreserved%说明，本程序演示了钟摆摆动过程，并同时给出简写振动方程%代码由董灿生编写，所有权利保留%时间2012年11月7日%%%%%%%%%%%%%%%%%%%%%%%%%Start%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear%重力加速度g=9.81;%Kgm/S^2%摆长l=1.8;%m%初相位phi_0=pi/2;%摆的质量m=1;%Kgtheta_m=pi/3;%最大摆角%圆频率w=sqrt(g/l);%周期T=2*pi/w;%时间向量t=0:0.05:

2、m*T;%简谐方程theta=theta_m*cos(w*t+phi_0);%小球x坐标x=l*sin(theta);%小球y坐标y=l-l*cos(theta);%以下为演示程序代码figureholdon[m,n]=size(x);plot(t,theta);plot(x,y);xlabel('lengthandtime');ylabel('thetaandheight');title('simplependulummodel');plot(0,l,'b+');%gridonfori=1:nh(1)=plot([0,x(i)],[l,y(i)]);h(2)=plot(x(i)

3、,y(i),'ro');h(3)=plot(t(i),theta(i),'gd');pause(0.1)ifi~=ndelete(h);endpause(0.1)endholdoff%%%%%%%%%%%%Theend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%