只有本站会员才能查看附件，请 登录
这一步是怎么转换的，我眼拙看不出来。
clear
clc
clf
t_0 = 0; %起始时间
t_f = 2*pi; %终止时间
x0 = [0; 0]; %初始值
[t, x] = ode23('f', [t_0,t_f], x0);%结果
[b]y = x(:, 1);
dy = x(:, 2);这是什么意思
plot(t, y, 'r-.', t, dy, 'b--')
legend('y', "y'")
%另建m函数文件f
function dx = f(t, x)
dx = [0, 1; -1, 1]*x + [0; 3*cos(t)];
end