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
MATLAB使用欧拉Euler法求解微分方程组 源程序代码.zip
版权申诉
169 浏览量
2022-11-16
21:16:29
上传
评论
收藏 1009B ZIP 举报
芝麻粒儿
- 粉丝: 6w+
- 资源: 2万+
最新资源
- Linux系统mysql安装配置教程.md
- Linux系统mysql安装配置教程.md
- 宜搭低代码开发时(高级)题库
- 51单片机设计-交通灯与串口模拟.zip,可作参考
- OPENCV学习总结,包含图像处理常用的数据类型和库函数
- [MacOS (AppleSilicon)] Proxyman-5.2.0
- 20240414-开源证券-开源证券量化评论-形态识别,均线的收敛与发散.txt
- 2023-04-06-项目笔记 - 第一百六十一阶段 - 4.4.2.159全局变量的作用域-159 -2024.06.11
- 数字图像处理-基于Qt 5.8.0.zip
- 2023-04-06-项目笔记 - 第一百六十一阶段 - 4.4.2.159全局变量的作用域-159 -2024.06.11
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈