Matlab数学建模实战——（Lokta-Volterra掠食者-猎物方程）

最新推荐文章于 2024-10-01 22:14:15 发布

沿途有李

最新推荐文章于 2024-10-01 22:14:15 发布

阅读量3.2k

点赞数 5

分类专栏： Matlab 文章标签：数学建模 matlab 开发语言

本文链接：https://blog.csdn.net/weixin_43475628/article/details/131620498

版权

Matlab 专栏收录该内容

7 篇文章

订阅专栏

文章通过三个步骤展示了如何使用Matlab解决不同的微分方程组问题，涉及变量替换、编程求解和结果分析。通过改变初值，观察到周期变化，例如问题1和问题2中周期分别为5和8。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

1.题目

请添加图片描述

问题1

该数学建模的第一问和第二问主要是用Matlab求解微分方程组，直接编程即可。

求解

Step1改写

y(1)=r
y(2)=f

Step2得y的导数

y(1).=2y(1)-ay(1)*y(2)
y(2).=-y(2)+a*y(1)*y(2)

Step3编程

clear;
a=0.01;
F=@(t,y)[2*y(1)-a*y(1)*y(2);
    -y(2)+a*y(1)*y(2)];
[t,y]=ode45(F,[0,10],[300,150]);
subplot(121);
plot(t,y(:,1));
xlabel('时间');
ylabel('r数量');
subplot(122);
plot(y(:,1),y(:,2));
xlabel('r数量');
ylabel('y数量');

Step4结果
在这里插入图片描述
从图中可以看出对应的周期tp大概是5

问题2

改变初值就行啦，直接进入step3

求解

Step3编程

clear;
a=0.01;
F=@(t,y)[2*y(1)-a*y(1)*y(2);
    -y(2)+a*y(1)*y(2)];
[t,y]=ode45(F,[0,10],[15,22]);
subplot(121);
plot(t,y(:,1));
xlabel('时间');
ylabel('r数量');
subplot(122);
plot(y(:,1),y(:,2));
xlabel('r数量');
ylabel('y数量');

Step4结果
在这里插入图片描述
从图中可以看出对应的周期tp大概是8

问题3

分析

首先列公式，给了u和v的定义，咱就对它求导，然后把它代入捕食者方程，再把r和f用u和v替换，根据题目所说忽略二阶项uv，这样可以算得
在这里插入图片描述

求解

Step1改写

y(1)=u
y(2)=v

Step2得y的导数

y(1).=-y(2)
y(2).=2y(1)

Step3编程

clear;
a=0.01;
F=@(t,y)[-y(2);
    2*y(1)];
[t,y]=ode45(F,[0,10],[0.1,0.1]);
subplot(121);
plot(t,y(:,1));
xlabel('时间');
ylabel('u');
subplot(122);
plot(y(:,1),y(:,2));
xlabel('u');
ylabel('v');