像这样
x0 = [0;
0);
tspan = [0, 10];
(t, x) =数值(@func tspan, x0);
次要情节(2,1,1)
情节(t) x (: 1));
网格(“上”);
%轴([0 10 1 1]);
包含(“时间[s]”);
ylabel (“x [m]”);
次要情节(2,1,2)
情节(t) x (:, 2));
网格;
%轴([0 10 1 1]);
包含(“时间[s]”);
ylabel (“dx / dt (m / s)”);
函数xdot = func (x ~)
M = 1.0;
K = 10;
C = 1;
l = 0.5;
g = 9.8;
xdot = [x (2);
- k * (x (1) - l) * (1) / mc * x (2) / M + g);%(1)不仅仅是x。不需要dt
结束