章绍辉数学建模第一章
- 格式:pdf
- 大小:539.62 KB
- 文档页数:12
捕食者与被捕食者模型——Logistic-Volterra模型摘要Logistic模型是最常用的模型之一,在其基础上又可以发展出许多其他数学模型,其重要性不言而喻,而Volterra模型则是经典的被捕食者与捕食者模型之一。
本文尝试结合两者,建立一个Logistic-Volterra模型,并做出数值解和分析。
关键词:Logistic模型 Volterra模型数值解一、问题的提出Volterra模型显示的被捕食者与捕食者系统存在着显著的周期振荡,而实际上,多数的捕食者与捕食者系统都是观察不到的。
尝试建立模型,描述这种现象。
二、符号说明r:被捕食者固有增长率d:捕食者固有死亡率a:捕食者掠取被捕食者的能力b:被捕食者供养捕食者的能力N1:被捕食者的最大环境容纳量N2:捕食者的最大环境容纳量三、模型假设1.在没有天敌的情况下,被捕食者数量增加的固有速度与被捕食者数量x和阻滞作用因子(1-x/N1)成正比,即dxdt =rx(1−xN1)2.在没有食物的情况下,捕食者数量减少的固有速度与捕食者数量y和阻滞作用因子(1+y/N2)成正比,即dydt =−dy(1+yN2)3.捕食者与被捕食者在同一环境下生存,它们的种群变化速度互相影响,影响因子应与它们相遇的频率成正比,即捕食导致被捕食者数量减少的速度为-axy,捕食导致捕食者数量增加的速度为bxy四、模型建立与求解1.Volterra模型的分析意大利数学家Volterra在上世纪20年代提出的Volterra模型:dxdt=rx−axydydt=−dy+bxy取r=1 d=0.5 a=0.1 b=0.02,运用matlab的ode45功能函数,做出数值解,并绘图分析。
图1被捕食者与捕食者随时间变化图图2捕食者与被捕食者相图从图形可以看出,捕食者与被捕食者共同生存,数量随时间作周期变化。
2.建立Logistic-Volterra模型在Volterra模型中的物种自身增长率中,考虑自身阻滞作用,即加入Logistic项,得到以下模型:dx dt =rx(1−xN1)−axydy dt =−dy(1+yN2)+bxy取r=1 d=0.5 a=0.1 b=0.02 N1=100 N2=25,运用matlab的ode45功能函数,做出数值解,并绘图分析。
习题3参考解答4. 某成功人士向学院捐献20万元设立优秀本科生奖学金,学院领导打算将这笔捐款以整存整取一年定期的形式存入银行,第二年一到期就支取,取出一部分作为当年的奖学金,剩下的继续以整存整取一年定期的形式存入银行……请你研究这个问题,并向学院领导写一份报告.解答 假设整存整取一年定期的年利率保持不变,记为r ,假设一到期就支取,取出b 元作为当年的奖学金,剩下的继续以整存整取一年定期的形式存入银行……记捐款存入银行之后第k 年一年定期到期日奖学金捐款账户余额为k x 万元,0x =20万元,则列式得1(1), 0,1,2,k k x r x b k +=+-=⋅⋅⋅.其解为()0(1), 0,1,2,k k x r x b r b r k =+-+=⋅⋅⋅ 平衡点为x b r =.因为r >0,所以如果0x b r >,即00b rx <<,k x 就会单调增加趋于无穷大,并且增加得越来越快;如果0x b r <,即0b rx >,k x 就会单调衰减(到零为止),并且减少得越来越快;如果0x b r =,即0b rx =,k x 就会保持不变,即0k x x ≡.如果取r =0.025,则b 的临界值为00.025200.5rx =⨯=(万元). 进一步,可编程分别计算当b =0.4、0.5、0.6、1以及2万元时账户总额k x 的具体变化过程,并绘图.程序:r=0.025; x=[20,20,20,20,20];b=[.4,.5,.6,1,2]; n=20;for k=1:nx(k+1,:)=x(k,:).*(1+r)-b;endplot(0:n,x(:,1),'k.',0:n,x(:,2),'kx',...0:n,x(:,3),'k^',0:n,x(:,4),'ks',0:n,x(:,5),'kp')axis([-1,n+1,0,25])legend('每年用0.4万元','每年用0.5万元',...'每年用0.6万元','每年用1万元','每年用2万元',3)title('奖学金捐款账户余额的演变,年利率2.5%')xlabel('第 k 年'), ylabel('账户余额(万元)')绘得的图形:第 k 年账户余额(万元)奖学金捐款账户余额的演变,年利率2.5%(略去给学院领导的报告,该报告要用非数学语言陈述上述分析)5. 有一位老人60岁时将养老金10万元以整存零取方式(指本金一次存入,分次支取本金的一种储蓄)存入,从第一个月开始每月支取1000元,银行每月初按月利率0.3%把上月结余额孳生的利息自动存入养老金. 请你计算老人多少岁时将把养老金用完?如果想用到80岁,问60岁时应存入多少钱?解答 假设月利率保持不变,记为r ,记养老金存入银行之后第k 月末账户总额为k x 元,从第一个月开始每月支取b 元. 则列式得1(1), 0,1,2,k k x r x b k +=+-=⋅⋅⋅.解得()0(1), 0,1,2,k k x r x b r b r k =+-+=⋅⋅⋅依题意有r =0.003,b =1000,0x =100000. 因为r >0,且0x b r <,所以k x 就会单调衰减(到零为止),并且减少得越来越快;若要0k x ≤,可以解得只需要()()()0log log log 1b r b r x k r --≥+ 所以令()(){}()0log log log 1K b r b r x r ⎡⎤=--+⎢⎥(上取整),则养老金在第K 个月恰好用完. 把具体数据代入,执行以下程序,算得K =120,即10万养老金恰好10年用:x=100000; r=0.003; b=1000;K=ceil((log(b/r)-log(b/r-x))/log(1+r))也可以用条件循环语句按差分方程迭代计算,直到0k x ≤停止. 程序如下:x=100000; r=0.003; b=1000; n=0;while x(n+1)>0n=n+1;x(n+1)=(1+r).*x(n)-b;endn如果养老金想用到80岁,即240x =0,那么()()()2400240111709081b r x r r +-==+.。
高等数学建模系列教材目录第一章数列与数学归纳法1.1 等差数列与等差数列的性质1.2 等比数列与等比数列的性质1.3 数列的求和公式1.4 数学归纳法的基本原理与应用第二章极限与连续2.1 极限的概念与性质2.2 极限运算的基本法则2.3 无穷小与无穷大2.4 一致连续与间断点的分类第三章导数与微分3.1 导数的概念与性质3.2 常见函数的导数及其应用3.3 高阶导数与导数的几何意义3.4 微分的概念与性质第四章微分中值定理与导数的应用4.1 罗尔中值定理与拉格朗日中值定理4.2 洛必达法则与洛必达不定型4.3 泰勒公式与泰勒展开式4.4 导数在几何与物理问题中的应用第五章不定积分与定积分5.1 不定积分的概念与性质5.2 基本积分公式及其运算法则5.3 定积分的概念与性质5.4 定积分的几何与物理应用第六章微分方程6.1 微分方程的概念与基本术语6.2 可分离变量方程与一阶线性方程6.3 高阶线性齐次与非齐次微分方程6.4 常微分方程的初值问题与边值问题第七章无穷级数与傅里叶级数7.1 级数的概念与性质7.2 数项级数、函数项级数与幂级数7.3 收敛判别法与级数的运算7.4 傅里叶级数的概念与性质第八章多元函数微分学8.1 多元函数的极限与连续性8.2 偏导数与全微分8.3 隐函数与参数方程求导8.4 多元函数的极值与条件极值第九章重积分与曲线曲面积分9.1 二重积分的概念与性质9.2 二重积分的计算与应用9.3 三重积分的概念与性质9.4 三重积分的计算与应用第十章曲线与曲面积分10.1 第一类曲线积分的概念与计算10.2 第二类曲线积分的概念与计算10.3 曲面积分的概念与计算10.4 斯托克斯公式与高斯公式第十一章偏微分方程与特殊函数11.1 偏微分方程的基本概念与分类11.2 常见偏微分方程及其解法11.3 常见特殊函数及其性质11.4 特殊函数在物理问题中的应用第十二章线性代数初步12.1 行列式的概念与性质12.2 矩阵的概念与运算12.3 矩阵的逆与初等变换12.4 线性方程组与矩阵的秩第十三章线性空间与线性映射13.1 线性空间的基本概念与性质13.2 线性子空间与线性独立性13.3 线性映射与线性变换13.4 基变换与线性映射的矩阵表示以上为《高等数学建模系列教材》的目录,涵盖了数学建模中所需的基础知识和方法。
第1章数学建模概论随着电子计算机的出现和科学技术的迅猛发展,数学的应用已不再局限于传统的物理领域,而正以空前的广度和深度逐步渗透到人类活动的各个领域。
生物、医学、军事、社会、经济、管理等各学科、各行业都涌现出大量的实际课题,亟待人们去研究、去解决。
利用数学知识研究和解决实际问题,遇到的第一项工作就是要建立恰当的数学模型,简称数学建模,数学建模正在越来越广泛地受到人们的重视。
从这一意义上讲,数学建模被看成是科学研究和技术开发的基础。
没有一个较好的数学模型就不可能得到较好的研究结果,所以,从这一意义上讲,建立一个较好的数学模型乃是解决实际问题的关键步骤之一。
1.1 数学模型与数学建模1.1.1 模型的概念在日常生活和工作中,人们经常会遇到或用到各种模型,如飞机模型、水坝模型、火箭模型、人造卫星模型、大型水电站模型等实物模型;也有文字、符号、图表、公式、框图等描述客观事物的某些特征和内在联系的模型,如模拟模型、数学模型等抽象模型。
模型是客观事物的一种简化的表示和体现,它应具有如下的特点:1.它是客观事物的一种模仿或抽象;它的一个重要作用就是加深人们对客观事物如何运行的理解,为了使模型成为帮助人们合理进行思考的一种工具,因此要用一种简化的方式来表现一个复杂的系统或现象。
2.为了能协助人们解决问题,模型必须具备所研究系统的基本特征和要素。
此外,还应包括决定其原因和效果的各个要素之间的相互关系。
有了这样的一个模型,人们就可以在模型内实际处理一个系统的所有要素,并观察它们的效果。
模型可以分为实物(形象)模型和抽象模型,抽象模型又可以分为模拟模型和数学模型。
对我们来说,最感兴趣的是数学模型。
与上述的各种各样的模型相对应的是它们在现实世界中的原型(原始参照物)。
所谓原型,是指人们研究或从事生产、管理的实际对象,也就是系统科学中所说的实际系统,如电力系统、生态系统、社会经济系统等。
而模型则是指为了某个特定目的,将原型进行适当地简化、提炼而构造的一种原型替代物。
1.(1) n=101;x1=linspace(-1,1,n); x2=linspace(-2,2,n);y1=[sqrt(1-x1.^2);-sqrt(1-x1.^2)];y2=[sqrt(4-x2.^2);-sqrt(4-x2.^2);sqrt(1-(x2.^2)/4);-sqrt(1-(x2.^2)/4)];plot(x1,y1) hold on; plot(x2,y2)title('椭圆x^2/4+y^2=1的内切圆和外切圆') axis equal-2.5-2-1.5-1-0.500.51 1.52 2.5-2-1.5-1-0.500.511.52椭圆x 2/4+y 2=1的内切圆和外切圆(2)x1=linspace(-2,2,101); x2=linspace(-2,8); axis equalplot(exp(x1),x1,x1,exp(x1),x2,x2)title('指数函数y=exp(x)和对数函数y=ln(x)关于y=x 对称')-2-112345678-2-1012345678指数函数y=exp(x)和对数函数y=ln(x)关于y=x 对称(3) hold onq=input('请输入一个正整数q;') for i=1:q for j=1:i if rem(j,i)plot(j/i,1/i) end end end0.10.20.30.40.50.60.70.80.9100.050.10.150.20.250.30.350.40.450.53.代码如下:n=input('请输入实验次数n=') k=0;for i=1:nx=ceil(rand*6)+ceil(rand*6); if x ==3|x==11 k=k+1; elseif x~=2&x~=7&x~=12y= ceil(rand*6)+ceil(rand*6); while y~=x&y~=7y=ceil(rand*6)+ceil(rand*6); end if y==7k=k+1; end end end从上表可看出打赌者赢的概率大约为。
数学建模章绍辉答案【篇一:第三次数学建模作业】数科院105 刘镜韶 20102201092 数科院105 蔡秋荣 20102201166 数科院104 梁浩坤 201022011004、不妨令第k年取出奖学金后,继续存在银行的捐款余额为xk,且银行的整存整取的利率为r,奖学金的金额为d万元,则由已知可得:xk+1 =(1+r)xk-d 故:其解为数列:xk =(x0-d/r)+d/r,且x0=20万元;①奖学金金额d0.6万元,让存在银行的捐款余额每年逐步增加;②奖学金金额d=0.6万元,让存在银行的捐款余额每年保持不变;③奖学金金额d0.6万元,让存在银行的捐款余额每年逐步减少;故对于不同的情况,不妨通过编程对比xk的变化趋势;程序:n=20;r=[0.03,0.03,0.03];x=[20,20,20];d=[0.45,0.6,0.75]; fork=1:nx(k+1,:)=x(k,:).*(1+r)-d; enddisp(本金为20万时不同的奖学金下余额的变化)disp(年 0.45万元0.6万元0.75万元) disp([(0:n),x]);plot(0:n,x(:,1),k^,0:n,x(:,2),ko,0:n,x(:,3),kv) axis([-1,n+1,14,25]) legend(d=0.45,d=0.6,d=0.75,2)title(本金为20万时不同的奖学金下余额的变化) xlabel(第k年),ylabel(余额) 其命令窗口显示结果为:年 0.45万元0.6万元0.75万元 020.000020.000020.00001.000020.150020.000019.85002.000020.304520.000019.69553.000020.463620.000019.53644.000020.627520.000019.37255.000020.796420.000019.20366.000020.970320.000019.02977.000021.149420.000018.8506 8.000021.333920.000018.66619.000021.523920.000018.4761本金为20万时不同的奖学金下余额的变化10.000021.719620.000018.2804 11.000021.921220.000018.078812.000022.128820.000017.871213.000022.342720.000017.657314.000022.562920.000017.4371 15.000022.789820.000017.210216.000023.023520.000016.976517.000023.264220.000016.735818.000023.512220.000016.4878 19.000023.767520.000016.2325 第k年20.000024.030620.000015.9694当利率r=3%时,且以整存整取一年定期的形式来存入银行时;由上述图像可知:①奖学金金额d≤0.6万元时,可以永久持续下去,实现可持续发展,即用20万元本金所得的利息作为奖学金。