数学建模报告46.docVIP

数学建模竞赛 一、问题简述： 问题一中要求检测A种鱼的数量随时间的变化的变化情况；问题二中给出了A,B两种鱼的微分方程组，要求我们根据已知条件求出A,B两种鱼随时间变化而产生的数量上的变化；问题三中给出了两种鱼在一段时间内数量的观测数据，要求我们根据问题二求出微分方程中的各个常数以及在特定时间点的两种鱼的数量。根据题目的要求，我们对问题一采用三次曲线拟合的方法，对问题二采用ODE45函数的方法，对问题三也采用ODE45函数的方法。 二、问题分析： 对问题一的分析，问题属于一元三次方程的数学问题，我们采用三次曲线拟合的方法，x=[25.81 35.47 51.08 60.45 65.87 92.82 120.64 149.20 184.54 222.62 257.68 292.53 325.26 352.38 385.63 406.92 427.30 433.87 439.42 451.38]，这些x的数据都随着时间的增大而递增，所以我们可以试着用方程来求出它的三次曲线。 对问题二的分析：第一步先在matlab中建立M文件C，编写一段函数，第二步在MatLab主窗口中调用ode45函数解微分方程组，第三步绘制曲线图，计算出t=10,20,30......100的A、B两种鱼的数量值。 问题三的分析：第一步，在已给的数据中任意抽取4组数据求取a0，第二步，MatLab中建立M文件L，编写一下函数，第三步，用已求的a0初始值调用lsqnonlin函数进行最小二乘求解，第四步，根据以求得的a0*，在MatLab中建立M文件C3，编写一下函数，第五步，根据A、B鱼类的两曲线循环周期，设定步长，调用ode45函数画出数值解曲线图，第六步，采用搜索算法在曲线上搜索出（x,y）的对应时间，第七步，抽取数据计算时间间隔 三：模型的建立与求解 （1）在某水域生长的一种鱼类（A种鱼），根据表格记录了一段时间内所观测到的该鱼类在若干观测观测的数量。建立数学模型，描述A种鱼的数量随时间变化的情况。 ①样本数据输入： x=[25.81 35.47 51.08 60.45 65.87 92.82 120.64 149.20 184.54 222.62 257.68 292.53 325.26 352.38 385.63 406.92 427.30 433.87 439.42 451.38]; t=0:1:19; ②综合比较，用3次曲线拟合样本数据效果较好，有MatLab编程得到3次方程系数： p=polyfit(t,x,3); px=poly2str(p,x) px =-0.14419 x^3 + 4.1232 x^2 - 4.6016 x + 33.8034 ③MatLab绘制拟合曲线和原始数据曲线： x=[25.81 35.47 51.08 60.45 65.87 92.82 120.64 149.20 184.54 222.62 257.68 292.53 325.26 352.38 385.63 406.92 427.30 433.87 439.42 451.38]; t=0:1:19; y=-0.14419*t.^3 + 4.1232*t.^2 - 4.6016*t + 33.8034;plot(t,x,c+,t,y,r) 图（1） ④最后得到的拟合方程： y=-0.14419*t.^3 + 4.1232*t.^2 - 4.6016*t + 33.8034 （2）根据给出的微分方程，计算t=10，20,30....100 时，两种鱼的数量。 ①MatLab中建立M文件C，编写一下函数： function x=C(t,x) x=[2*x(1)-0.01*x(1)*x(2);-x(2)+0.01*x(1)*x(2)]; ②在MatLab主窗口中调用ode45函数解微分方程组： t_final=100;x0=[10,100];[t,x]=ode45(C,[0:10:t_final],x0); ③绘制曲线图，计算出t=10,20,30......100的A、B两种鱼的数量值： plot(t,x(:,1),r,t,x(:,2),b,t,x(:,1),ro,t,x(:,2),bo) 根据在一段时间内两种鱼的数量的观测数据，两种鱼的数量变化符合问题二中的微分方程组，求其中的常数以及初始时刻两种鱼的数量。 定理一：对于参数完全相同的生态系统，其系统常数与观测时间单位长度存在反比例关系，即。 证明：在volterra模型中： 把，代