微分方程建模是数学建模的重要方法,因为许多实际问题的数学描述将导致求解微分方程的定解问题。在高教杯数学建模竞赛中每年都会有一道微分方程建模问题,大体上可以按以 下几步:
1. 根据实际要求确定要研究的量(自变量、未知函数、必要的参数等)并确定坐标系。
2. 找出这些量所满足的基本规律(物理的、几何的、化学的或生物学的等等)。
3. 运用这些规律列出方程和定解条件。
常见的列方程方法有:
(i)按规律直接列方程
在数学、力学、物理、化学等学科中许多自然现象所满足的规律已为人们所熟悉, 并直接由微分方程所描述。如牛顿第二定律、放射性物质的放射性规律等。我们常利用 这些规律对某些实际问题列出微分方程。
(ii)微元分析法与任意区域上取积分的方法
自然界中也有许多现象所满足的规律是通过变量的微元之间的关系式来表达的。对于这类问题,我们不能直接列出自变量和未知函数及其变化率之间的关系式,而是通过微元分析法,利用已知的规律建立一些变量(自变量与未知函数)的微元之间的关系式,然后再通过取极限的方法得到微分方程,或等价地通过任意区域上取积分的方法来建立微分方程。
(iii)模拟近似法
在生物、经济等学科中,许多现象所满足的规律并不很清楚而且相当复杂,因而需要根据实际资料或大量的实验数据,提出各种假设。在一定的假设下,给出实际现象所满足的规律,然后利用适当的数学方法列出微分方程。
示例:发射卫星为什么用三级火箭问题
采用运载火箭把人造卫星发射到高空轨道上运行,为什么不能用一级火箭而必须用多级火箭系统? 下面通过建立运载火箭有关的数学模型来回答上述问题。
火箭是一个复杂的系统,为了使问题简单明了,我们只从动力系统和整体结构上分析,并且假设引擎是足够强大的。
卫星进入 600km 高空轨道时,火箭必须的最低速度,首先将问题理想化,假设:
(i)卫星轨道是以地球中心为圆心的某个平面上的圆周,卫星在此轨道上以地球引力作为向心力绕地球作平面匀速圆周运动;
(ii)地球是固定于空间中的一个均匀球体,其质量集中于球心;
(iii)其它星球对卫星的引力忽略不计。
建模与求解:
设地球半径为 R ,质量为 M ;卫星轨道半径为 r ,卫星质量为m 。 根据假设,卫星只受到地球的引力,由牛顿万有引力定律可知其引力大小为:
其中G 为引力常数。
为消去常数G ,把卫星放在地球表面:
代入得:
其中 g 为重力加速度,一般为10或者9.81m/s。
根据假设(i),若卫星围绕地球作匀速圆周运动的速度为v ,则其向心力为 , 因为卫星所受的地球引力就是它作匀速运动的向心力,故有:
由此便推得卫星距地面为(r − R)km ,必须的最低速度的数学模型为 :
即要把卫星送入离地面 600km 高的轨道,火箭的末速度最低应为 7.6km/s。
火箭的简单模型是由一台发动机和一个燃料仓组成。燃料燃烧产生大量气体从火箭 末端喷出,给火箭一个向前的推力。火箭飞行要受地球引力、空气阻力、地球自转与公转等的影响,使火箭升空后作曲线运动。为使问题简化,假设:
(i)火箭在喷气推动下作直线运动,火箭所受的重力和空气阻力忽略不计。
(ii)在t 时刻火箭质量为m(t) ,速度为v(t) ,且均为时间t 的连续可微函数;
(iii)从火箭末端喷出气体的速度(相对火箭本身)为常数u 。
火箭—卫星系统的质量可分为三部分:mp (有效负载,如卫星),mF(燃料质量), ms(结构质量,如外壳、燃料容器及推进器)。一级火箭末速度上限主要是受目前技术条件的限制,假设:
从前面对问题的假设和分析可以看出:火箭推进力自始至终在加速着整个火箭,然而随着燃料的不断消耗,所出现的无用结构质量也在随之不断加速,作了无用功,故效益低,浪费大。 所谓理想火箭,就是能够随着燃料的燃烧不断抛弃火箭的无用结构。下面建立它的数学模型。
假设:在(t,t + Δt) 时段丢弃的结构质量与烧掉的燃料质量以α 与1−α 的比例同时进行。
建模与分析:由动量守恒定律,有
多级火箭是从末级开始,逐级燃烧,当第i 级燃料烧尽时,第i +1级火箭立即自动点火,并抛弃已经无用的第i 级。我们用mi 表 示第i 级火箭质量,mp 表示有效负载。为了简单起见,先作如下假设:
实际上,由于受技术条件的限制,采用四级或四级以上的火箭,经济效益是不合算 的,因此采用三级火箭是最好的方案。
最佳结构设计:下面将考虑当用n 级火箭发射卫星时的最佳结构,即使m0 / mp最小的结构 。