这一步是怎么转换的，我眼拙看不出来。


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