算法大全第13章_微分方程建模（以火箭为数学模型案例）

普大帝

你好！我是陪你一起进阶人生的普大帝！愿你成才！祝你成长！ 7 m$ h3 o! s5 R为大家更新一些算法类的辅助资料，大家在想学习时，或者比赛急需时就可以按照对应的名字找到对应的算法，加以应用了。大家按照下图所示箭头处点击主题，就可以查看到其他算法类内容了，或者搜索框内输入算法大全，也可以搜索到对应的贴子，本篇为第13篇。

微分方程建模是数学建模的重要方法，因为许多实际问题的数学描述将导致求解微分方程的定解问题。把形形色色的实际问题化成微分方程的定解问题，大体上可以按以下几步：
7 u# y# A# m- c6 E6 Q1. 根据实际要求确定要研究的量(自变量、未知函数、必要的参数等)并确定坐标系。
2. 找出这些量所满足的基本规律(物理的、几何的、化学的或生物学的等等)。
. U7 d/ S- F# s3. 运用这些规律列出方程和定解条件。
/ x3 M: V: _6 `0 q列方程常见的方法有：
（i）按规律直接列方程
1 e$ z: G2 x5 q. Q8 M在数学、力学、物理、化学等学科中许多自然现象所满足的规律已为人们所熟悉，并直接由微分方程所描述。如牛顿第二定律、放射性物质的放射性规律等。我们常利用这些规律对某些实际问题列出微分方程。
3 c2 P; d$ G' X7 m（ii）微元分析法与任意区域上取积分的方法自然界中也有许多现象所满足的规律是通过变量的微元之间的关系式来表达的。对于这类问题，我们不能直接列出自变量和未知函数及其变化率之间的关系式，而是通过微元分析法，利用已知的规律建立一些变量（自变量与未知函数）的微元之间的关系式，然后再通过取极限的方法得到微分方程，或等价地通过任意区域上取积分的方法来建立微分方程。
（iii）模拟近似法
- y) ?, \% p6 k. y: Y在生物、经济等学科中，许多现象所满足的规律并不很清楚而且相当复杂，因而需要根据实际资料或大量的实验数据，提出各种假设。在一定的假设下，给出实际现象所满足的规律，然后利用适当的数学方法列出微分方程。在实际的微分方程建模过程中，也往往是上述方法的综合应用。不论应用哪种方法，通常要根据实际情况，作出一定的假设与简化，并要把模型的理论或计算结果与实际情况进行对照验证，以修改模型使之更准确地描述实际问题并进而达到预测预报的目的。
0 N4 e6 _  y" [/ X1 E) q" q, w: @本章将利用上述方法讨论具体的微分方程的建模问题。
5 J+ N$ Z8 y/ L7 N' T§1 发射卫星为什么用三级火箭采用运载火箭把人造卫星发射到高空轨道上运行，为什么不能用一级火箭而必须用
, t; ]5 t0 z. G7 D3 j多级火箭系统？
下面通过建立运载火箭有关的数学模型来回答上述问题。
火箭是一个复杂的系统，为了使问题简单明了，我们只从动力系统和整体结构上分析，并且假设引擎是足够强大的。
1.1 为什么不能用一级火箭发射人造卫星
. Y/ F3 J6 n% y7 ^6 o  R下面用三个数学模型回答这个问题
1.1.1 卫星进入 600km 高空轨道时，火箭必须的最低速度首先将问题理想化，假设：
（i）卫星轨道是以地球中心为圆心的某个平面上的圆周，卫星在此轨道上以地球
) I  n2 {0 L4 H5 F) ^引力作为向心力绕地球作平面匀速圆周运动；
0 x: q. E% _0 ~  N8 j: K# h（ii）地球是固定于空间中的一个均匀球体，其质量集中于球心；
$ w1 `4 k7 e5 C（iii）其它星球对卫星的引力忽略不计。建模与求解：设地球半径为 R ，质量为 M ；卫星轨道半径为 r ，卫星质量为 m 。
' V; S8 O/ n/ _* r. G根据假设（ii）和（iii），卫星只受到地球的引力，由牛顿万有引力定律可知其引力大小为

9 H' _* W, |. B1 P( n, Z1 }  o9 m
1.3 多级火箭卫星系统
理想火箭是设想把无用结构质量连续抛弃以达到最佳的升空速度，虽然这在目前的技术条件下办不到，但它确为发展火箭技术指明了奋斗目标。目前已商业化的多级火箭卫星系统便是朝着这种目标迈进的第一步。多级火箭是从末级开始，逐级燃烧，当第 i级燃料烧尽时，第 1 + i 级火箭立即自动点火，并抛弃已经无用的第 i 级。我们用im 表示第 i 级火箭质量，pm 表示有效负载。为了简单起见，先作如下假设：
& Y" }2 y3 ]& s
% T$ m4 v6 B" ^& k8 ?! ?下载资源查看完整版内容
% r& G1 @$ A- j# M. h
, [# M7 b7 }: A- v
3 @; L, o# {3 Q9 y) B6 I3 N
* R# q2 f) b: M. }

数模老哥

发表回复666666

数模老哥

666666666666
: ~3 b, D, ~8 j7 ^$ v$ m; V- R