問題出現(xiàn)在給的 tspan 參數(shù)錯(cuò)誤, 0 和 5 之間添個(gè)空格就行.

% 轉(zhuǎn)換為方程組：y1 = y, y2 = y'
ode_fun = @(x, y) [y(2); -3*y(2) - 2*y(1)];
% 數(shù)值求解
[x_ode, y_ode] = ode45(ode_fun, [0 5], [1; 0]);
% 解析解
y_exact = @(x) 2*exp(-x) - exp(-2*x);
% 繪圖
figure;
plot(x_ode, y_ode(:,1), 'b--', x_ode, y_exact(x_ode), 'r-');
xlabel('x');
ylabel('y');
legend('數(shù)值解', '解析解');
title('二階齊次方程解對(duì)比');