%% 学习目标:欧拉法求解微分方程组
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轴的标签
matlab数理统计和数据分析及优化求解:40 扩展参考: matlab欧拉Euler法求解微分方程组.zip
版权申诉
33 浏览量
2023-08-06
22:01:46
上传
评论
收藏 867B ZIP 举报
skyJ
- 粉丝: 2477
- 资源: 2038
最新资源
- 使用 C 语言实现的计算非负整数的阶乘
- 2011-2021最新版本北京大学数字普惠金融指数(PKU-DFIIC).xlsx
- 县域数字乡村指数2018-2020(1).xlsx
- Docker容器配置进阶
- tensorflow-gpu-2.7.4-cp37-cp37m-manylinux2010-x86-64.whl
- 多段线、 圆、弧转多段线(仅我可见)
- tensorflow-2.7.2-cp38-cp38-manylinux2010-x86-64.whl
- 李慧琴C语言基础部分.zip
- yeyue-p8Yi4-ve4a83792.apk
- tensorflow-gpu-2.7.3-cp38-cp38-manylinux2010-x86-64.whl
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈