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法求解微分方程组 源程序代码.rar
版权申诉
63 浏览量
2023-11-07
18:12:46
上传
评论
收藏 740B RAR 举报
公众号:智慧方案文库
- 粉丝: 2187
- 资源: 1万+
最新资源
- 以下是一些适用于英语六级作文的万能句型模板,涵盖了引言、正文和结论部分的各类表达方式.docx
- MATLAB中的非线性规划
- 进行C语言面试资格确认是招聘过程中一个重要的步骤,目的是确保候选人具备足够的C语言编程能力和知识.docx
- Java 轻量级的集群负载均衡设计
- 纹身师个人网站模板.jpg
- 在C语言中,连接两个字符串(即将一个字符串附加到另一个字符串的末尾)通常可以使用标准库中的 `strcat` 函数.docx
- 数据库管理工具:dbeaver-ce-23.1.1-stable.x86-64.rpm
- 以下是几个具体竞赛题目的详细解答,包括建模思路、方法和步骤 .docx
- 一份关于全国大学生建模大赛的相关教程!!
- 以下是关于计算机网络和现代通信组网的详细教程、案例和相关项目的推荐.docx
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈