第9讲-差分方程.pptVIP

数学建模与数学实验 后勤工程学院数学教研室 差分方程 实验目的 实验内容 2、掌握用数学软件求解差分方程问题。 1、直观了解差分方程基本内容。 1、差分方程的基本理论。 3、实验作业。 2、用数学软件求解差分方程问题。 差分方程的基本理论 1 差分方程模型 解：在第一个月时，只有一对小兔子，过了一个月，那对兔子成熟了，在第三个月时便生下一对小兔子，这时有两对兔子。再过一个月，成熟的兔子再生一对小兔子，而另一对小兔子长大，有三对小兔子。如此推算下去，我们便发现一个规律：  55 34 21 10 34 21 13 9 21 13 8 8 13 8 5 7 8 5 3 6 5 3 2 5 3 2 1 4 2 1 1 3 1 1 0 2 1 0 1 1 兔子总数(对) 成熟兔子(对) 初生兔子(对) 时间(月) 2 差分方程的解法 1 ) 常系数线性齐次差分方程的解法 解 差分方程的特征方程: 2) 常系数线性非齐次差分方程的解法 应用举例 例8 如图正方形ABCD的四个顶点各有一人。在某一时刻，四人同时出发以匀速v=1米/秒按顺时针方向追逐下一人，如果他们始终保持对准目标，则最终按螺旋状曲线于中心点O.试求出这种情况下每个人的行进轨迹。 O C D A B 图1 求解过程: 1. 建立平面直角坐标系:A(x1,y1), B(x2,y2), C(x3,y3), D(x4,y4)。 MATLAB程序如下： v=1; dt=0.05; x=[0 0 10 10]; y=[0 10 10 0]; for i=1:4 plot(x(i),y(i),.);hold on end d=20; while(d0.1) x(5)=x(1);y(5)=y(1); for i=1:4 d=sqrt((x(i+1)-x(i))^2+(y(i+1)-y(i))^2); x(i)=x(i)+v*dt*(x(i+1)-x(i))/d; y(i)=y(i)+v*dt*(y(i+1)-y(i))/d; plot(x(i),y(i),.);hold on end end 所得到得轨迹模拟图如图1 To Matlab （ cx1.m ） 例8 最优捕鱼策略 生态学表明，对可再生资源的开发策略应在事先可持续收获的前提下追求最大经济效益.考虑具有4个年龄组：1龄鱼，...4龄鱼的某种鱼.该鱼类在每年后4个月季节性集中产卵繁殖.而按规定，捕捞作业只允许在前8个月进行，每年投入的捕捞能力固定不变，单位时间捕捞量与各年龄组鱼群条数的比例称为捕捞强度系数.使用只能捕捞3、4龄鱼的13mm网眼的拉网，其两个捕捞强度系数比为0.42：1.渔业上称这种方法式为固定努力量捕捞. 有如下问题需要解决: 1) 分析如何实现可持续捕获(即每年开始捕捞时渔场中各年龄组鱼群不变)，并在此前提下得到最高收获量； 2) 合同要求某渔业公司在5年合同期满后鱼群的生产能力不能受到太大的破坏,承包时各年龄组鱼群数量为122，29.7，10.1，3.29(×109条),在固定努力量的捕捞方式下,问该公司应采取作怎样的捕捞策略,才能使总收获量最高. 这里只讨论问题2),即可持续捕获策略模型. 以一年为一个离散化的单位时间. 因为3龄鱼与4龄鱼捕捞强度系数比为0.42：1