Matlab 欧拉法求解常微分方程组 - du/dt = -4u + 6v, dv/dt = 3u - 7v

% 定义函数 function [u, v, t] = euler_method() % 初始化参数 h = 0.01; % 步长 t = 0:h:1; % 时间步 u = zeros(size(t)); % u的解 v = zeros(size(t)); % v的解 u(1) = 1; % 初始条件 v(1) = 2; % 初始条件

% 欧拉法求解 for i = 1:length(t)-1 u(i+1) = u(i) + h*(-4u(i) + 6v(i)); % du/dt = -4u + 6v v(i+1) = v(i) + h*(3u(i) - 7v(i)); % dv/dt = 3u - 7v end

% 绘制图像 subplot(2,1,1); plot(t,u,'r'); xlabel('t'); ylabel('u'); title('du/dt = -4u + 6v');

subplot(2,1,2); plot(t,v,'b'); xlabel('t'); ylabel('v'); title('dv/dt = 3u - 7v'); end

% 调用函数 [u, v, t] = euler_method();