MATLAB程序分享实现单摆在外力矩作用下的动画源程序-MATLAB实现单摆在外力矩作用下的动画 源程序代码.rar

MATLAB程序分享实现单摆在外力矩作用下的动画源程序-MATLAB实现单摆在外力矩作用下的动画 源程序代码.rar 程序代码见附件

clear;clc;close all
l=4;
m=1; 
g=9.81;
zeta=0.2; 
I=m*l^2;
w0=sqrt（g/l）;
c=zeta*2*w0*I; 
dt=0.2/w0;

theta（1）=0.6;
dtheta（1）=0; 

Ek=0.5*m*（l*dtheta（1））^2; 
Ep=l*（1-cos（theta（1）））*m*g;

figure          
h1=subplot（211）;
set（h1‘XLim‘[-8 8]‘YLim‘[-5 0]）
set（gca‘drawmode‘‘fast‘）
line_handle1=plot（[0 l*sin（theta（1））][0 -l*cos（theta（1））]‘-o‘）; 
axis equal
title（‘Pendulum Simulation‘）
blank=[‘                ‘]; 
parameters={[‘Length‘ blank（1:14） num2str（l）];[‘Mass‘ blank（1:16） num2str（m）]; ...
            [‘Damping ratio‘ blank（1:4） num2str（zeta）];[‘Gravity‘ blank（1:14） num2str（g）]};
text（3-1.5parameters）
text_show={[‘Time‘ blank（1:10） ‘0‘];[‘Theta‘ blank（1:9） num2str（theta（1））]; ...
           [‘dTheta/dt‘ blank（1:4） num2str（dtheta（1））];[‘Energy‘ blank（1:7） num2str（Ek+Ep）]};
text_handle=text（-7.5-3.5text_show）;

h2=subplot（212）;
set（h2‘XLim‘[-0.8 0.8]‘YLim‘[-0.8 0.8]） 
set（gca‘drawmode‘‘fast‘）
line_handle2=plot（theta（1）dtheta（1））;
axis equal
title（‘theta-dtheta curve‘） 
xlabel（‘theta （rad）‘） 
ylabel（‘dtheta （rad/s）‘）

N=80;
for k=1:N
    Tw=5*sin（2*（k-1）*dt）;
    ddtheta（k）=（Tw/I）-（g/l）*sin（theta（k））-（c/（m*l^2））*dtheta（k）;
    dtheta（k+1）=dtheta（k）+dt*ddtheta（k）; 
    theta（k+1）=theta（k）+dt*dtheta（k）; 
    
    Ek=0.5*m*（l*dtheta（k+1））^2;
    Ep=l*（1-cos（theta（k+1）））*m*g;
    
    set（line_handle1‘xdata‘[0 l*sin（theta（k+1））]‘ydata‘[0 -l*cos（theta（k+1））]）
    set（h1‘XLim‘[-8 8]‘YLim‘[-5 0]） 
    text_show={[‘Time‘ blank（1:10） num2str（k*dt）];[‘Theta‘ blank（1:9） num2str（theta（k+1））]; ...
               [‘dTheta/dt‘ blank（1:4） num2str（dtheta（k+1））];[‘Energy‘ blank（1:7） num2str（Ek+Ep）]};
    set（text_handle‘string‘text_show）
    
    set（line_handle2‘xdata‘theta‘ydata‘dtheta）
    set（h2‘XLim‘[-0.8 0.8]‘YLim‘[-0.8 0.8]）
    drawnow
    pause（0.001）
end

msgbox（‘MATLAB编程答疑，请加QQ: 1530497909‘‘MATLAB答疑‘‘help‘）
web http://url.cn/NSFcAs -browser

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1980  2014-05-17 21:02  MATLAB实现单摆在外力矩作用下的动画 源程序代码\chengxu.m

     文件        333  2014-05-17 20:48  MATLAB实现单摆在外力矩作用下的动画 源程序代码\Read Me.url

     文件        294  2014-02-06 11:28  MATLAB实现单摆在外力矩作用下的动画 源程序代码\联系我们.txt

     目录          0  2014-05-17 21:01  MATLAB实现单摆在外力矩作用下的动画 源程序代码

----------- ---------  ---------- -----  ----

                 2607                    4