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 %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(t,x) %绘制x曲线 hold on plot(t,y,'r') %绘制y曲线 xlabel('time') %设置x轴的标签 ylabel('value') %设置y轴的标签 legend({'x','y'}) %设置图例 title('time evolution plot') %设置图形的标题 figure plot(x,y) %绘制phase图 title('phase plane plot') %设置标题 xlabel('x') %设置x轴的标签 ylabel('y') %设置y轴的标签 msgbox('MATLAB编程答疑，请加QQ: 1530497909','MATLAB答疑','help') web http://url.cn/TKcdXk -browser