2025年大学《数学与应用数学》专业题库- 数学在航空航天领域的应用

2025年大学《数学与应用数学》专业题库——数学在航空航天领域的应用考试时间：______分钟总分：______分姓名：______一、已知飞行器在忽略空气阻力和地球自转影响时，其运动轨迹满足微分方程$\vec{r}''=-\frac{GM}{r^3}\vec{r}$，其中$\vec{r}$为飞行器相对地心的位矢，$G$为万有引力常数，$M$为地球质量，$r=|\vec{r}|$。设飞行器在$t=0$时刻位于点$(r_0,0,0)$，初始速度为$\vec{v}(0)=(0,v_{0y},0)$。1.将该矢量微分方程转化为关于$r$和$\theta$(极坐标中的角)的常微分方程组（$r=|\vec{r}|$,$\theta$为极角，假设$z=0$平面为初始运动平面）。2.推导出轨道方程$r=\frac{h}{1+e\cos\theta}$，并说明式中$h$和$e$的物理意义。其中$h=r_0v_{0y}$为特定角动量，$e$为轨道离心率。二、飞行器绕地球做椭圆轨道运动，其近地点距离为$r_p$，远地点距离为$r_a$。利用题一推导出的轨道方程相关结论，计算该飞行器在近地点和远地点时的速度大小。三、考虑一个由四个质点组成的简单模型，代表航天器四个关键部件。质点质量分别为$m_1,m_2,m_3,m_4$，位置向量分别为$\vec{r}_1,\vec{r}_2,\vec{r}_3,\vec{r}_4$。定义航天器的质心位矢$\vec{R}=\frac{1}{M}\sum_{i=1}^4m_i\vec{r}_i$，其中$M=m_1+m_2+m_3+m_4$为航天器总质量。1.推导出描述该航天器整体运动的方程（牛顿第二定律应用于质心）。2.若已知$\vec{r}_1(t),\vec{r}_2(t),\vec{r}_3(t),\vec{r}_4(t)$，如何计算该航天器在时刻$t$的总动量$\vec{P}(t)$和总动能$T(t)$？四、在飞行器姿态控制中，常用欧拉角描述飞行器的姿态。设飞行器绕其质心的旋转运动由三个刚体转动惯量$I_1,I_2,I_3$决定，角速度矢量为$\vec{\omega}=(\omega_x,\omega_y,\omega_z)$。陀螺仪测得的总角速度（相对于惯性空间）为$\vec{\omega}_e=(\omega_{ex},\omega_{ey},\omega_{ez})$。假设存在一个由陀螺仪和飞轮组成的控制回路，产生的控制力矩$\vec{M}_c$与角速度偏差$\vec{\omega}_e-\vec{\omega}$成正比，即$\vec{M}_c=K(\vec{\omega}_e-\vec{\omega})$，其中$K$为比例常数。1.建立描述该简化姿态动力学和控制系统的方程（运动方程）。2.分析该控制系统的稳定性（定性分析，如判断特征值实部的符号），并说明比例常数$K$对系统稳定性的影响。五、空气动力学中，飞行器表面某点处的压力$p$可以用流函数$\psi$表示，满足拉普拉斯方程$\nabla^2\psi=0$。对于层流，速度势$\phi$也满足拉普拉斯方程$\nabla^2\phi=0$。假设在飞行器附近某区域，流场可近似为二维，速度分量$u=\frac{\partial\phi}{\partialx}$,$v=-\frac{\partial\phi}{\partialy}$。1.推导出压力$p$与速度势$\phi$的关系（伯努利方程）。2.若已知该区域的速度分布为$u=U_0(1+\frac{y}{a})e^{-x/a}$，其中$U_0$和$a$为常数，且$v=0$，求该区域的速度势$\phi(x,y)$和压力分布$p(x,y)$。假设远场压力为$p_\infty$，且$\phi=0$时$x\to\infty$。六、航天任务设计常涉及轨道机动。假设一航天器从圆形轨道（半径$r_0$）进行霍曼转移，前往另一个半径为$r_1$的圆形轨道。忽略非引力因素。1.求在第一个Burn（加速）中，航天器需要增加的速度$\Deltav_1$。2.求在第二个Burn（减速）中，航天器需要减少的速度$\Deltav_2$。3.求整个霍曼转移过程所需的总$\Deltav$以及航天器在转移轨道远地点和近地点的速度大小。七、卫星接收地面站发送的信号会受到大气层的影响。假设信号在穿过大气层时，其衰减服从指数规律，即信号强度$I$随高度$h$的变化关系为$I(h)=I_0e^{-\kappah}$，其中$I_0$为地面（$h=0$）的信号强度，$\kappa$为大气衰减系数。若已知地面信号强度为$I_0=1$瓦米${}^{-2}$，大气衰减系数$\kappa=0.1$公里${}^{-1}$，求卫星在高度$h=400$公里处接收到的信号强度。试卷答案一、1.将矢量方程$\vec{r}''=-\frac{GM}{r^3}\vec{r}$投影到极坐标系（$r,\theta$），有$\vec{r}=r\hat{e}_r$，$\vec{r}''=\ddot{r}\hat{e}_r+r\ddot{\theta}\hat{e}_\theta+2\dot{r}\dot{\theta}\hat{e}_\theta$。其中$\hat{e}_r,\hat{e}_\theta$是径向和角向单位矢量，满足$\hat{e}_r\cdot\hat{e}_r=1,\hat{e}_r\cdot\hat{e}_\theta=0,\hat{e}_\theta\cdot\hat{e}_\theta=1,\dot{\hat{e}}_r=\dot{\theta}\hat{e}_\theta,\ddot{\hat{e}}_r=\ddot{\theta}\hat{e}_\theta-\dot{\theta}^2\hat{e}_r,\dot{\hat{e}}_\theta=-\dot{\theta}\hat{e}_r,\ddot{\hat{e}}_\theta=-\ddot{\theta}\hat{e}_r-2\dot{\theta}^2\hat{e}_\theta$。代入原方程：$\ddot{r}\hat{e}_r+r\ddot{\theta}\hat{e}_\theta+2\dot{r}\dot{\theta}\hat{e}_\theta=-\frac{GM}{r^3}(r\hat{e}_r)$。整理得：$\ddot{r}=-\frac{GM}{r^3}r=-\frac{GM}{r^2}$；$r\ddot{\theta}+2\dot{r}\dot{\theta}=0$。消去$\dot{\theta}$，得到：$r^2\ddot{\theta}+r\dot{r}\dot{\theta}'+2r\dot{r}^2=0$。或者写成$\frac{1}{r}\frac{d}{dt}(r^2\dot{\theta})=0$，即$r^2\dot{\theta}=h$，$h$为常量。因此，转化后的常微分方程组为：$\dot{r}=v_r,\quad\ddot{r}=-\frac{GM}{r^2}$$\dot{\theta}=\frac{h}{r^2},\quad\dot{h}=0$（$h$为常量）其中$v_r=\frac{dr}{dt}$。2.由$\dot{r}=\frac{dr}{d\theta}\dot{\theta}=\frac{dr}{d\theta}\frac{h}{r^2}$，得到$\frac{dr}{d\theta}=\frac{r^3}{h}$。分离变量积分：$\int_{r_0}^{r}r^{-3}dr=\int_{0}^{\theta}\frac{1}{h}d\theta$。得到：$-\frac{1}{2r^2}\Big|_{r_0}^{r}=\frac{\theta}{h}$，即$\frac{1}{2r^2}-\frac{1}{2r_0^2}=\frac{\theta}{h}$。整理得：$\frac{1}{r^2}=\frac{1}{r_0^2}+\frac{2\theta}{h}$。进一步整理，令$e=\frac{h}{r_0v_{0y}}=\frac{h}{r_0\cdot\frac{h}{r_0}}=1$（对于$t=0$时刻位于$(r_0,0,0)$，$\vec{v}(0)=(0,v_{0y},0)$的初始条件，$e=1$对应抛物线轨道，更一般的椭圆轨道需考虑其他初始条件，但此形式常表示椭圆方程，需理解为$h=r_0v_{0y}$是特定角动量）。最终得到轨道方程：$r=\frac{h}{1+e\cos\theta}$。其中，$h=r_0v_{0y}$是特定角动量，表示飞行器相对于地心的角动量大小。$e$是轨道离心率，决定轨道形状：$e=0$为圆轨道，$0<e<1$为椭圆轨道，$e=1$为抛物线轨道，$e>1$为双曲线轨道。二、由题一，近地点$r=r_p$，远地点$r=r_a$。根据机械能守恒定律，飞行器在任意位置的机械能$E$相同，即$\frac{1}{2}mv^2-\frac{GMm}{r}=E$。在近地点，$v=v_p$，$r=r_p$，有$\frac{1}{2}mv_p^2-\frac{GMm}{r_p}=E$。在远地点，$v=v_a$，$r=r_a$，有$\frac{1}{2}mv_a^2-\frac{GMm}{r_a}=E$。两式相减，消去$E$：$\frac{1}{2}m(v_p^2-v_a^2)=\frac{GMm}{r_p}-\frac{GMm}{r_a}=GMm\left(\frac{1}{r_p}-\frac{1}{r_a}\right)$。$v_p^2-v_a^2=\frac{2GMm}{m}\left(\frac{1}{r_p}-\frac{1}{r_a}\right)=2GM\left(\frac{1}{r_p}-\frac{1}{r_a}\right)$。由角动量守恒定律，近地点和远地点的角动量相等，$mr_pv_p=mr_av_a$，即$r_pv_p=r_av_a$。将$v_p=\sqrt{v_p^2}$和$v_a=\sqrt{v_a^2}$代入角动量守恒式，平方得$r_p^2v_p^2=r_a^2v_a^2$。代入$v_p^2-v_a^2$的表达式：$r_p^2\left(\frac{2GM}{r_p}-\frac{2GM}{r_a}\right)=r_a^2v_a^2$。整理得：$2GMr_a-2GMr_p=r_a^2v_a^2$。$v_a^2=\frac{2GM(r_a-r_p)}{r_a^2}$。$v_a=\sqrt{\frac{2GM(r_a-r_p)}{r_a^2}}$。同理，将$v_p^2-v_a^2$的表达式代入$r_p^2v_p^2=r_a^2v_a^2$，得到$r_a^2v_a^2=r_p^2\left(\frac{2GM}{r_p}-\frac{2GM}{r_a}\right)$。整理得：$2GMr_p-2GMr_a=r_p^2v_p^2$。$v_p^2=\frac{2GM(r_p-r_a)}{r_p^2}$。$v_p=\sqrt{\frac{2GM(r_p-r_a)}{r_p^2}}$。三、1.航天器整体运动方程即牛顿第二定律应用于质心。质心$\vec{R}$的加速度为$\vec{R}''=\frac{d^2\vec{R}}{dt^2}$。总外力为$\vec{F}_{ext}$。根据牛顿第二定律：$\vec{F}_{ext}=M\vec{R}''$。其中$M=m_1+m_2+m_3+m_4$，$\vec{R}=\frac{1}{M}(m_1\vec{r}_1+m_2\vec{r}_2+m_3\vec{r}_3+m_4\vec{r}_4)$。因此，描述航天器整体运动的方程为：$M\vec{R}''=\vec{F}_{ext}$。2.总动量$\vec{P}(t)$定义为航天器各部件动量的矢量和：$\vec{P}(t)=m_1\vec{v}_1+m_2\vec{v}_2+m_3\vec{v}_3+m_4\vec{v}_4$。其中$\vec{v}_i=\frac{d\vec{r}_i}{dt}$。由于质心速度$\vec{V}=\frac{d\vec{R}}{dt}$，有$\vec{v}_i=\vec{V}+\vec{v}_{iR}$，其中$\vec{v}_{iR}$是第$i$个部件相对于质心的速度。总动量可写为：$\vec{P}(t)=M\vec{V}+\sum_{i=1}^4m_i\vec{v}_{iR}$。由质心定义$\vec{R}=\frac{1}{M}\sum_{i=1}^4m_i\vec{r}_i$，对时间求导得$\vec{V}=\frac{1}{M}\sum_{i=1}^4m_i\vec{v}_i$。将$\vec{v}_i=\vec{V}+\vec{v}_{iR}$代入上式：$\vec{P}(t)=M\vec{V}+\sum_{i=1}^4m_i(\vec{V}+\vec{v}_{iR})=M\vec{V}+M\vec{V}+\sum_{i=1}^4m_i\vec{v}_{iR}=2M\vec{V}+\sum_{i=1}^4m_i\vec{v}_{iR}$。由于$\vec{R}=\frac{1}{M}\sum_{i=1}^4m_i\vec{r}_i$，质心加速度$\vec{R}''=\vec{V}$，意味着$\sum_{i=1}^4m_i\vec{v}_{iR}=0$（质心平动参考系下，内力之和为零，或对质心求和）。因此，总动量$\vec{P}(t)=2M\vec{V}=2M\frac{d\vec{R}}{dt}$。总动能$T(t)$定义为航天器各部件动能之和：$T(t)=\frac{1}{2}m_1v_1^2+\frac{1}{2}m_2v_2^2+\frac{1}{2}m_3v_3^2+\frac{1}{2}m_4v_4^2$。其中$v_i=|\vec{v}_i|=|\vec{V}+\vec{v}_{iR}|$。$T(t)=\frac{1}{2}\sum_{i=1}^4m_i(\vec{V}+\vec{v}_{iR})\cdot(\vec{V}+\vec{v}_{iR})$$=\frac{1}{2}\sum_{i=1}^4m_i(\vec{V}\cdot\vec{V}+2\vec{V}\cdot\vec{v}_{iR}+\vec{v}_{iR}\cdot\vec{v}_{iR})$$=\frac{1}{2}\sum_{i=1}^4m_iV^2+\sum_{i=1}^4\frac{1}{2}m_iv_{iR}^2+\sum_{i=1}^4m_i\vec{V}\cdot\vec{v}_{iR}$。由于$\sum_{i=1}^4m_i\vec{v}_{iR}=0$，且$\vec{V}\cdot\vec{v}_{iR}$是相对速度与质心速度的点积，其总和不一定为零，但通常在质心系下讨论时，此项代表整体旋转动能，若只考虑整体平动，则此项为零或被包含在旋转动能中。若假设部件相对质心速度很小或只考虑整体平动动能，则$\sum_{i=1}^4m_i\vec{V}\cdot\vec{v}_{iR}=0$。因此，总动能近似为$T(t)=\frac{1}{2}MV^2+\sum_{i=1}^4\frac{1}{2}m_iv_{iR}^2=\frac{1}{2}MV^2+T_{rot}$，其中$T_{rot}$为相对质心的旋转动能。四、1.航天器绕质心旋转的运动方程为$I\vec{\omega}''=\vec{M}_{ext}$，其中$I$是惯量张量，$\vec{M}_{ext}$是外部总力矩。总力矩$\vec{M}_{ext}=\vec{M}_c+\vec{M}_{ext,other}$，其中$\vec{M}_c=K(\vec{\omega}_e-\vec{\omega})$是控制力矩，$\vec{M}_{ext,other}$是其他外力矩（如科里奥利力矩、重力梯度力矩等，题目未明确给出，可视为零或包含在$\vec{M}_{ext}$中）。运动方程可写为：$I\vec{\omega}''=K(\vec{\omega}_e-\vec{\omega})+\vec{M}_{ext,other}$。将$\vec{\omega}''=\dot{\vec{\omega}}$代入，得到：$I\dot{\vec{\omega}}=K(\vec{\omega}_e-\vec{\omega})+\vec{M}_{ext,other}$。这是一个关于$\vec{\omega}$的一阶线性非齐次微分方程组。写成分量形式，假设$I=\text{diag}(I_1,I_2,I_3)$，$\vec{\omega}=(\omega_x,\omega_y,\omega_z)$，$\vec{\omega}_e=(\omega_{ex},\omega_{ey},\omega_{ez})$：$I_1\dot{\omega}_x=K(\omega_{ex}-\omega_x)+M_{ext,x}$$I_2\dot{\omega}_y=K(\omega_{ey}-\omega_y)+M_{ext,y}$$I_3\dot{\omega}_z=K(\omega_{ez}-\omega_z)+M_{ext,z}$2.为分析稳定性，考虑系统在平衡点的线性化方程。平衡点是指$\vec{\omega}=\vec{\omega}_e$时系统保持不变的状态。设$\vec{\omega}=\vec{\omega}_e+\vec{\delta\omega}$，其中$\vec{\delta\omega}$是小扰动角速度。当$\vec{\omega}=\vec{\omega}_e$时，方程为$I\dot{\vec{\omega}}=K(\vec{\omega}_e-\vec{\omega}_e)+\vec{M}_{ext,other}=\vec{M}_{ext,other}$。如果$\vec{M}_{ext,other}=0$，则平衡点是$\vec{\omega}=\vec{\omega}_e$。将$\vec{\omega}=\vec{\omega}_e+\vec{\delta\omega}$代入运动方程$I\dot{\vec{\omega}}=K(\vec{\omega}_e-\vec{\omega})+\vec{M}_{ext,other}$，得到：$I\dot{\vec{\delta\omega}}=K(\vec{\omega}_e-(\vec{\omega}_e+\vec{\delta\omega}))+\vec{M}_{ext,other}$$I\dot{\vec{\delta\omega}}=-K\vec{\delta\omega}+\vec{M}_{ext,other}$。若$\vec{M}_{ext,other}=0$，则扰动方程为$I\dot{\vec{\delta\omega}}=-K\vec{\delta\omega}$。该方程的特征方程为$\det(I\lambda+K)=0$。特征值$\lambda$为$-K/I$。惯量张量$I$的对角元素$I_1,I_2,I_3$均为正，比例常数$K$也为正（代表控制力矩与速度偏差成正比，且起抑制作用的预期）。因此，特征值$\lambda_1=-K/I_1,\lambda_2=-K/I_2,\lambda_3=-K/I_3$均为负实数。根据线性稳定性理论，所有特征值都具有负实部意味着系统在平衡点$\vec{\omega}=\vec{\omega}_e$是李雅普诺夫意义下稳定的。比例常数$K$的大小影响特征值的负值大小，即影响系统恢复平衡的快慢。$K$越大，特征值越负，系统越快恢复稳定。但$K$过大可能导致系统振荡或响应过快引起饱和，工程上需选择合适的$K$值。五、1.速度势$\phi$满足拉普拉斯方程$\nabla^2\phi=0$。速度分量$u=\frac{\partial\phi}{\partialx}$,$v=-\frac{\partial\phi}{\partialy}$。压力$p$满足伯努利方程：$p+\frac{1}{2}\rhou^2+\rhogh=\text{常数}$，其中$\rho$为流体密度，$g$为重力加速度，$h$为高度。在二维水平流动问题中，通常不考虑重力项或假设$gh$为常数可吸收进常数项，简化为$p+\frac{1}{2}\rhou^2=\text{常数}$。在流函数$\psi$满足的拉普拉斯方程$\nabla^2\psi=0$的流场中，速度与流函数的关系为$u=\frac{\partial\psi}{\partialy}$,$v=-\frac{\partial\psi}{\partialx}$。结合伯努利方程，有$p+\frac{1}{2}\rho\left(\left(\frac{\partial\psi}{\partialy}\right)^2+\left(-\frac{\partial\psi}{\partialx}\right)^2\right)=p+\frac{1}{2}\rho|\nabla\psi|^2=\text{常数}$。这表明在流函数满足拉普拉斯方程的流场中，压力$p$只与速度大小有关，与速度方向无关。沿流线（$\psi=\text{常数}$），压力变化与速度平方的变化有关。2.已知$u=U_0(1+\frac{y}{a})e^{-x/a}$,$v=0$。首先，验证该速度场是否满足流函数方程$\frac{\partial\psi}{\partialy}=u=U_0(1+\frac{y}{a})e^{-x/a}$和$\frac{\partial\psi}{\partialx}=-v=0$。积分$\frac{\partial\psi}{\partialy}=U_0(1+\frac{y}{a})e^{-x/a}$，得到$\psi=\intU_0(1+\frac{y}{a})e^{-x/a}dy=U_0\inte^{-x/a}dy+\frac{U_0}{a}\intye^{-x/a}dy$。$\psi=U_0ye^{-x/a}+\frac{U_0}{a}\left(y\frac{a}{a}e^{-x/a}-\int\frac{a}{a}e^{-x/a}dy\right)=U_0ye^{-x/a}+\frac{U_0}{a}(ye^{-x/a}+ae^{-x/a})=(U_0y+U_0a)e^{-x/a}+C(x)$。由于$\frac{\partial\psi}{\partialx}=-v=0$，要求$\frac{\partial}{\partialx}\left[(U_0y+U_0a)e^{-x/a}+C(x)\right]=0$。$-U_0(y+a)\frac{1}{a}e^{-x/a}+C'(x)=0$。要使此式恒等于零，需$C'(x)=U_0(y+a)\frac{1}{a}e^{-x/a}$。由于$C'(x)$只能是$x$的函数，而右侧含有$y$，要使等式恒成立，必须有$y+a=0$，这不可能。因此，积分常数$C(x)$必须包含所有与$x$相关的部分，使得$\frac{\partialC(x)}{\partialx}=U_0(y+a)\frac{1}{a}e^{-x/a}$。这暗示积分形式可能需要调整。更常见的形式是直接写出速度势。另一种处理方式是直接假设速度势$\phi$，使得其梯度给出给定速度场。设$\phi=f(x,y)$。$\frac{\partial\phi}{\partialx}=-v=0\implies\phi$只与$y$有关，或$\phi=\phi(y)$。$\frac{\partial\phi}{\partialy}=u=U_0(1+\frac{y}{a})e^{-x/a}\implies\phi(y)=\intU_0(1+\frac{y}{a})e^{-x/a}dy=U_0ye^{-x/a}+U_0ae^{-x/a}+C$。令$C=0$（在$x\to\infty$时$\phi=0$的简化条件下确定），得到$\phi(y)=U_0(y+a)e^{-x/a}$。现在求压力分布。使用简化伯努利方程$p+\frac{1}{2}\rhou^2=p_\infty$（假设远场压力为$p_\infty$）。$p=p_\infty-\frac{1}{2}\rhou^2$。$u=U_0(1+\frac{y}{a})e^{-x/a}$。$u^2=U_0^2(1+\frac{y}{a})^2e^{-2x/a}$。$p=p_\infty-\frac{1}{2}\rhoU_0^2(1+\frac{y}{a})^2e^{-2x/a}$。六、1.霍曼转移轨道由两个共面、连续的圆轨道（或椭圆轨道）连接而成，其中一个是初始轨道，另一个是目标轨道。第一个Burn（加速）发生在初始轨道上，使航天器进入一条椭圆转移轨道的近地点，该近地点与初始圆轨道相同。第二个Burn（减速）发生在椭圆转移轨道的远地点，使航天器进入目标圆轨道。初始条件：航天器在半径为$r_0$的圆形轨道上，速率为$v_0=\sqrt{\frac{GM}{r_0}}$。第一个Burn（加速）：在初始圆轨道上，沿速度方向施加$\Deltav_1$，使航天器进入椭圆转移轨道。转移轨道的近地点为$r_0$，远地点为$r_1$。根据机械能守恒和角动量守恒，在椭圆轨道的近地点和远地点，速率满足：$\frac{1}{2}mv_{p,trans}^2-\frac{GMm}{r_0}=\frac{1}{2}mv_{a,trans}^2-\frac{GMm}{r_1}$。$r_0v_{p,trans}=r_1v_{a,trans}$。解得：$v_{p,trans}=v_0=\sqrt{\frac{GM}{r_0}}$。$v_{a,trans}=v_0\