MATLAB使用欧拉Euler法求解微分方程组 源程序代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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轴的标签 复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：MATLAB atlab matla Euler 微分方程 legend figure 源程序 title hold

谢谢分享。

好贴啊

非常感谢分享

学了matlab的确可以秒做很多数学题，特别是线性代数，哈哈