MATLAB使用欧拉Euler法求解微分方程组

MATLAB使用欧拉Euler法求解微分方程组 部分源码 clear;clc c=2/3; %设置c的值 x(1)=0.1; %设置x初值为0.1 y(1)=0.3; %设置y初值为0.3 h=0.05; %设置步长为0.05

clear;clc
c=2/3;  %设置c的值
x（1）=0.1;  %设置x初值为0.1
y（1）=0.3;  %设置y初值为0.3
h=0.05;    %设置步长为0.05

%for循环：根据Euler法，求解微分方程组
for i=1:1000
    x（i+1）=x（i）+h*（x（i）*（c-x（i）/y（i）））;
    y（i+1）=y（i）+h*（y（i）*（1-y（i））-x（i）*y（i））;
end

t=0:h:1000*h;   %计算时间
plot（tx）       %绘制x曲线
hold on       
plot（ty‘r‘）   %绘制y曲线
xlabel（‘time‘）  %设置x轴的标签
ylabel（‘value‘） %设置y轴的标签
legend（{‘x‘‘y‘}）  %设置图例
title（‘time evolution plot‘）  %设置图形的标题
figure   
plot（xy）   %绘制phase图
title（‘phase plane plot‘）  %设置标题
xlabel（‘x‘）   %设置x轴的标签
ylabel（‘y‘）   %设置y轴的标签

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         637  2018-04-21 18:38  MATLAB使用欧拉Euler法求解微分方程组 源程序代码\chengxu.m
     目录           0  2018-04-21 18:38  MATLAB使用欧拉Euler法求解微分方程组 源程序代码\