单摆模型MATLAB程序

%Simple pendulum model %All rights reserved %说明，本程序演示了钟摆摆动过程，并同时给出简写振动方程 %代码由董灿生编写，所有权利保留 %时间2012年11月7日 %%%%%%%%%%%%%%%%%%%%%%%%%Start%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%% clear %重力加速度 g=9.81;%Kgm/S^2 %摆长 l=1.8;%m %初相位 phi_0=pi/2; %摆的质量 m=1;%Kg theta_m=pi/3;%最大摆角 %圆频率 w=sqrt(g/l); %周期 T=2*pi/w; %时间向量 t=0:0.05:m*T; %简谐方程 theta=theta_m*cos(w*t+phi_0); %小球x坐标 x=l*sin(theta); %小球y坐标 y=l-l*cos(theta); %以下为演示程序代码 figure hold on [m,n]=size(x); plot(t,theta); plot(x,y); xlabel('length and time'); ylabel('\theta and height'); title('simple pendulum model'); plot(0,l,'b+'); % grid on for i=1:n
h(1)=plot([0,x(i)],[l,y(i)]); h(2)=plot(x(i),y(i),'ro'); h(3)=plot(t(i),theta(i),'gd');
pause(0.1) if i~=n
d源自文库lete(h); end
pause(0.1) end hold off %%%%%%%%%%%%The end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%