微分方程数值解Euler法

最新推荐文章于 2023-04-23 15:10:27 发布

weixin_30861797

最新推荐文章于 2023-04-23 15:10:27 发布

阅读量178

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/wdfrog/p/5872742.html

本文通过具体的MATLAB代码演示了如何使用欧拉方法求解一阶微分方程dy/dt=1+y的数值解，并与解析解进行对比。此外，还介绍了Logistic增长模型的微分方程形式及其数值解法。

微分方程:dy/dt=1+y; 解是y=2exp(x)-1;

clc
clear
figure(2)
dx=0.1;
x=0:dx:8;
y=zeros(size(x));

x(1)=0;
y(1)=1;
for i=1:length(x)-1
%   z= y(i)+ (1+y(i)) * dx;
%   y(i+1)= y(i)+ ( (1+y(i) + 1+ z) /2) * dx;
  y(i+1)= y(i)+ (1+y(i)) * dx;
end

plot(x,y)
grid on
hold on
ezplot('2*exp(x)-1',[0,8])

View Code

logistic增长模型

%微分方程入下
%dp/dt=0.001*(100-p)*p;
clc
clear
figure(2)
dx=0.1;
x=0:dx:100;
y=zeros(size(x));

x(1)=0;
y(1)=10;
for i=1:length(x)-1
%   z= y(i)+ (1+y(i)) * dx;
%   y(i+1)= y(i)+ ( (1+y(i) + 1+ z) /2) * dx;
%   y(i+1)= y(i)+ (1+y(i)) * dx;
z=y(i) +0.001 *(100-y(i))*y(i)*dx;
y(i+1)=y(i) + ((0.001 *(100-y(i))*y(i) + 0.001*(100-z)*z )/2  ) *dx;
end

plot(x,y)
grid on
hold on