Matlab学习笔记（三）

单摆（二）

一、单摆的运动方程

for n=1:iterations
    % Iteratively solve equations of motion using Euler's Method
    theta(n+1,:)=theta(n,:)+omega(n,:)*tStep; % new angular position
    omega(n+1,:)=omega(n,:)+alpha(n,:)*tStep; % new angular velocity
    alpha(n+1,:)=(-g*sin(theta(n+1,:)))/l-c*omega(n+1,:); % new angular acceleration
    
    T(n+1,:)=m*g*cos(theta(n+1,:))+m*l*(omega(n+1,:)).^2; % tension in cable
end

二、单摆的运动方程

% g_pendu.m
function dx=dynamics_pendulum(t,x)
mass=1; % define the mass of the bob
l=1; % define the length of the cord
g=9.81; % gravitational acceleration
c = 0.5;
% define the state-space model
% dx=[x(2,1); -(g/l)*sin(x(1,1))];% 不考虑阻尼
dx=[x(2,1); -(g/l)*sin(x(1,1))-c*x(2,1)];
end

% G_pendulum.m
[time_vector2,solution]=ode45(@g_pendu,time_vector,x0);