pendulumstateeq.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

function xdot=pendulumstateeq(t, x, flag, par)
% 考虑空气阻力的单摆系统的状态方程函数
% 输入: t 当前计算时刻,  flag 此处不用
% x为2#1矩阵, x(1)为当前摆锤的切向速度; x(2)为当前角位移
% par为4#1系统参数矩阵
% par(1) 为重力加速度参数 g
% par(2) 为空气阻力比例系数 k
% par(3) 为摆锤质量 m
% par(4) 为摆线长度 L
xdot=zeros(2,1);                    % 状态变量矩阵初始化
g=par(1);k=par(2);m=par(3);L=par(4);% 系统参数设置
xdot(1)=g*sin(x(2))-k/m*x(1);       % 状态方程组
xdot(2)=-1/L*x(1);