pendulumstateeqlinear.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 15 行

function xdot=pendulumstateeqlinear(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);% 系统参数设置
A=[ -k/m, g; -1/L, 0];  % 状态矩阵赋值
B=0; 
f=0;                    % 无输入信号
xdot=A*x+B*f;           % 通用的状态方程