m_chap2_ex3_sub.m

来自「simulink建模与仿真 基础入门源代码」· M 代码 · 共 17 行

function y=f(t,y)

m=70;
k=20;
a1=1;   a2=1;
g=10;
y0=[-30;0];

x=y(1);       %  x代表蹦极者的位移（在以桥梁面为基准面、以指向地面方向为正方向的坐标系中）
x_dot=y(2);   %  x_dot 代表 x 的速度
if x>0
    y=[0,1;-k/m,-a1/m-(a2/m)*abs(x_dot)]*y+[0;g];
else
    y=[0,1;0,-a1/m-(a2/m)*abs(x_dot)]*y+[0;g];
end