研究生数学建模竞赛优秀论文B

《研究生数学建模竞赛优秀论文B》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、一.问题的重述考虑航天器在仅受到地球万有引力、航天器自身发动机作用力的作用下作平面运动，将地球和航天器视为质点，建立航天器运动的数学模型。显然这样的数学模型在精度上是远远不能满足实际需要的，在其他耍求糟确制导等有关高科技的实际问题屮，我们都而临着类似的问题：我们必须建立高精度的数学模型，必须高精度地估计模型中的大批参数，因为只有这样的数学模型才能解决实际问题，而不会出现差之毫厘，结果却失之千里的情况。由于航天器的问题太复杂，本题仅考虑较简单的确定高精度参数问题。假设有一个生态系统，其中含有两种生物，即：A生物和B生物，其中A生物是捕食者，B生物是

2、被捕食者。假设/时刻捕食者A的数口为x(r),被捕食者B数目为y(r),它们Z间满足以下变化规律:⑴国+也)心)]初始条件为:y，(t)=y(t)[_a3^a4x(ty”仏)=巾其屮^,(1<^<6)为模型的待定参数。通过对此生态系统的观测，可以得到相关的观测数据。要利用冇关数据，解决以下问题：1)在观测数据无误差的情况下，若已知勺，求其它5个参数乙伙=1,3,4,5,6)?2)若勺也未知，至少需要多少组观测数据，才能确定参数^(1<^<6)?3)在观测资料有误差(时间变量不含有误差)的情况下，确定参数^.(1<^<6)在某种意义下的最优解，并

3、与仿真结果比较，进而改进数学模型。4)假设连观测资料的吋间变量也含有误差，确定参数色在某种意义下的最优解。二、航天器运动模型的建立考虑航天器在仅受到地球力有引力、航天器自身发动机作用力的作用下作平面运动，将地球和航天器视为质点，由理论力学可知，一个刚体在空间的运动可以看作质心的移动，因此可以应用质心运动定理來研究刚体质心的移动规律。以地球中心为原点，建立直角坐标系，航犬器绕地球飞行，可以出现在该直角坐标系屮四个彖限的任意一个Z内。平而直角坐标系如图lo符号说明如下：航天器在X方向的速度航天器在y方向的速度F2——航天器发动机作用力，为控制变量6Z

4、——万有引力与X轴正方向的夹角0——航天器发动机作用力与X轴正方向的夹角G——初始时刻兀0，儿—航天器初始位置Vv0航天器X方向初速度vv0航天器y方向初速度航天器受的万有引力存方向指向地球中心（原点），航犬器受推力坊的方向与X轴正方向成0角。将斥和尸2投影到该直角处标系上，见图1图1航天器受力分解图初始条件为Mo)=儿叫(G)=叫0Vv(r0)=Vv0x(r),y(f)都是关于时间的位置函数，航天器在x方向的分速度即兀“)对时间t求导，航天器在y方向的分速度即y(/)对吋间t求导，航天器在x方向的加速度即兀(/)对时间t求二阶导，航天器在y方向

5、的加速度即y(/)对时间t求二阶导，根据牛顿第二定律有方程(3)和(4)。由此建立的航天器模型如下:兀d/(Z)=Vyvf(t}=^=Z卩_g—乙一-『E)AmmD+m)jm)+m)()上=空空_M也—2mmm)+y2(“jm)+m)显然这样的数学模型在精度上是远远不能满足实际需要的，在其他要求精确制导等有关高科技的实际问题屮，我们都面临着类似的问题：我们必须建立高精度的数学模型，必须高精度地估计模型中的人批参数，因为只冇这样的数学模型才能解决实际问题，而不会出现茅Z毫厘，结果却失Z千里的情况。这时所建立数学模型的精度就成了数学模型的生命线。例如

6、上述问题屮的航天器述耍受到地球质量分布不均匀所引起的摄动力，大气阻力，日、刀及其它星球的摄动引力的影响，以及航天器发动机为调整航天器自身姿态运作吋作用力的影响。这样不但数学模型十分复杂，而且在这些数学模型中还要涉及到许多重要的参数，如地球的引力场模型就有许多待定参数。不仅如此，在对航大器进行测量时，还涉及到观测站的地理位置以及设备的系统误差等参数。为此人们耍设法利用长期积累的丰富的观测资料，高精度确定这些重要的参数。由于航天器的问题太复杂，下面木题仅考虑较简单的确定高精度参数问题。三、捕食者与被捕食者生态系统问题的分析题中假设有一个生态系统，含有

7、两种生物，A生物和B生物，A生物是捕食者，B生物是被捕食者。假设r时刻捕食者A的数目为x(r),被捕食者B数目为y(f),它们之间满足以下变化规律:x/(r)=x(z)[aI+3(“]⑴[禺+*")]初始条件为:Go)=^5该模型屮，e捕食者独自存在时死广率，«,<0；也被捕食者对捕食者的供养能力«2>0；勺是被捕食者的独立生存增长率，巾>0；夠是捕食者掠取被捕食者的能力,勺<0。〔2]这个方程就是生态系统中被捕食者与捕食者的volterra模型，务（】彳比6）为模型的待定参数。对于该模型理论上不存在解析解，因此我们不能通过参数拟合确定模型的参数

8、。Volterra模型在给定参数和初始值的情形下可以采用数值积分获得任意时间点的数值解。根据volterra模型进行一些公式推导如下：两