MATLAB 常微分方程数值解法:显式欧拉法
% 显式欧拉法求解常微分方程 clear all; clc;
% 设置时间步长和时间区间 dt = 0.01; t = 0:dt:1;
% 初始化u和v的值 u = zeros(size(t)); v = zeros(size(t)); u(1) = 1; v(1) = 2;
% 显式欧拉法求解常微分方程 for i = 2:length(t) u(i) = u(i-1) + dt*(-4u(i-1) + 6v(i-1)); v(i) = v(i-1) + dt*(3u(i-1) - 7v(i-1)); end
% 绘制图像 plot(t,u,'-r',t,v,'-b'); xlabel('t'); ylabel('u,v'); legend('u','v');
原文地址: https://www.cveoy.top/t/topic/oRjY 著作权归作者所有。请勿转载和采集!