常微分方程机敏问答[4] #20210622

---

本专栏主要作个人复习自测，有相关知识预备的同学也可作复习用。不保证无相应基础的人士能看明白。
万一考试考到了，或者对你的学习有较大帮助，一键三连不过分吧（斜眼笑）

多维系统

0. 四个角上四只乌龟分别盯紧”下一只“，以匀速运动：直觉上是（）螺线。为了说明解存在唯一（从而对称），我们可以考虑（）维复数值系统或者（）维实数值系统，再考察解的存在唯一性。
1. 三维空间中，物体加速度始终和速度、位矢同平面（线性相关），则物体是否在同一平面内运动？如果$\ddot{\vec{r}}(t)=\alpha (t)\dot{\vec{r}}(t)+\beta (t)\dot{\vec{r}}(t)$，且$\vec{r}\times \dot{\vec{r}}$在某时刻附近不为0呢？（提示：在同一平面内运动等价于位移垂直于什么向量？）
2. 高阶和多维有何关系？
3. 自治方程$F(y,y^{\prime },y^{\prime \prime })=0$怎么降为一阶、一维的方程？
4. $x^{\prime \prime }=f(x)$的物理意义是什么？由此怎么作相图？相图上”逆时针转“物理意义是什么？
5. $x^{\prime \prime }=a(x-x^3/6)$（单摆展开到三阶项）时相图如何？（注：这能反映单摆一次沿一个方向摆动角度超过$2\pi$的情况，而只看一阶项不行）
6. 二维情况，皮卡存在定理的解存在区域形如什么？李氏条件怎么表达？

答案

0. 等角，4，8
1. 不一定（如：先停再走）。一定。设$\vec{n}=\vec{r}\times \dot{\vec{r}}$（“法向量”，其在某时刻附近不为0），$\vec{n}$初值${\vec{n}}_0$，则在同一平面内运动等价于$\vec{r}\perp {\vec{n}}_0$. 现直接将$\vec{r}\cdot {\vec{n}}_0$作为时间的函数$f$，则$f^{\prime \prime }(t)=\alpha (t)f(t)+\beta (t)f^{\prime }(t)$，再根据存在唯一性定理，$f(t)=0$是唯一可能的解。
2. 把$y^{(n)}$看成变量，则可将高阶方程升维，并降至一阶。
3. $y^{\prime }=z,F(y,z,dz/dx)=F(y,z,z_y^{\prime }z)=0$
4. 提示：画出势函数，并特别注意正负。略。
5. 设$x^{\prime }=y$，再在二维系统中画使得$y^{\prime }=0$的曲线、画使得$x^{\prime }=0$的曲线，标记出各个区域$x^{\prime },y^{\prime }$符号，定性考察。
6. 根据范数选择不同，可能是各种柱形。表达李氏条件时也要说明范数选择。

连续依赖性

0. 考察解对参数连续依赖性。“$y$一致满足李氏条件”中一致是什么意思？结论中所谓“解连续”的闭区域是什么？
1. 接上，皮卡序列的某一项形如${\varphi }_n(x,\lambda )$，即以$\lambda$为参数的$x$的一元函数。关于$n$维向量参数的连续性实际上相当于（）元函数的连续性。
2. 接上，皮卡序列一致收敛的闭区域是什么？
3. 接上，为什么初值也是参数？把初值看成参数对结论的闭区域有何影响？
4. 解的存在唯一性说明映射$(x,y)\mapsto (x,\varphi (x;\eta =y))$（此记号意为初值$\eta$取$y$时，解在$x$处取值是$\varphi (x;\eta =y)$）是连续的（）射，从而这个映射的逆映射就相当于把解曲线族“拉直”。
5. 大范围对初值连续依赖性：已知初值问题有解$\varphi$，满足局部（一致）李氏条件，在（）（填区间的性质）区间内考察，在任何一点$(x_0,\varphi (x_0))$作小扰动范围内（即初值变成$(x_0,\varphi (x_0)+\delta )$），解在闭区域上（关于参数）连续。为什么不是（一维）闭区间？
6. 接上，利用（）定理把局部李氏条件拓展到整体。
7. 接上，此时作皮卡序列时，第0项和证明局部情况连续依赖性中的第0项有什么异同？
8. 接上，为了防止皮卡序列迭代过程中超出满足（）条件的范围，需要取比（）细得多的”管子“。
9. $y^{\prime }=f(x,y)$，$f$在区域上连续，且经过区域$R$内任一点解曲线存在唯一，则微分方程的解对初值（局部）连续依赖。证明用反证法，在闭区间上找“坏点”的聚点，并利用（）引理找一致收敛子列（并利用唯一性）得矛盾。

答案

0. 关于$x,\lambda$一致。$|x|\le h,|\lambda -{\lambda }_0|\le c$，其中$\lambda$可取的范围就是前述“一致”的范围。
1. $n+1$
2. 同0.中闭区域。（回忆事实：一致收敛保持连续性）
3. 把初值不为0的初值问题平移一下（做一个减法），从而变成函数$f(x,y)$不同了，但初值为0的初值问题，则其中减的量是“参数”。
   结论的闭区域缩小。举例：可以留一半余地给初值的变化，即$|x-x_0|\le \frac{h}{2}=\frac{min(a,b/M)}{2}$而$|\eta -y_0|\le b/2$.
4. 单
5. 有界闭。注意关于参数也要连续，所以考察的当然不是一元的连续性。
6. 有限覆盖
7. 不同：不是常数，而是“管”的“中心线”。一个相同点：两者中都出现参数，因为考察的是对参数的连续依赖性！
8. （整体）李氏，满足整体李氏条件的“管子”
9. Arzela-Ascoli.
   此处我们再次强调：得到一致收敛的子列不需要李氏条件。

多维线性微分方程组

0. 写出多维线性微分方程组的（列）向量形式。它的满足任何可能初值条件的解在整个区间上存在且唯一，该事实用到对于任意确定的$x$，都有$\dot{\vec{y}}$关于$\vec{y}$（）。
1. （）微分方程组的解的集合构成线性空间。如何考察它的维数？线性空间总有零元素此处说明什么？什么是一组函数的线性无关？在齐次线性微分方程组中一组解线性无关和初值的线性无关有什么关系？
2. 基本解组和通解有何关系？如果向量都是列向量，那么如何利用矩阵书写？
3. 接上，朗斯基行列式$W(x)$是关于$x$的函数，它和2.有何关系？$W$满足什么一阶线性微分方程？为什么？
4. 接上，讨论$W(x)$为0与线性无关的关系。
5. 齐次线性微分方程组解矩阵的导数是什么？
6. 考察解矩阵左乘或右乘常数非奇异矩阵。
7. 直接背诵非齐次线性微分方程组的特解公式，分析怎么退化到一维。该公式相比一维情况，实用意义有何不同？
8. $y^{\prime }=y+z,z^{\prime }=z/x$的求解：这是（）微分方程组，所以可以求出两个（）的特解得到通解。比如先用第二式求解$z$的所有可能解，再代入另一式。
9. 为什么$(1,0,0),(x,0,0),(x^2,0,0)$不可能同时满足任何一个三阶齐次线性微分方程组？
10. 为了证明基解矩阵完全决定齐次线性微分方程组（当然可能多个基解矩阵对应同一个方程组），回忆线性无关情况下（）（人名）行列式恒不为零，就可以考察任一点$x_0$处情况，并用（）（方程组种类）方程组有$n$个无关解证明（）恒为零矩阵。
11. 用常数变易法考察$\dot{\vec{y}}=A\vec{y}+f(x,\vec{y})$.

答案

0. $\dot{\vec{y}}=A(x)\vec{y}+\vec{e}(x)$，注意带$(x)$的记号表示$x$的函数构成的矩阵或向量。是仿射变换从而李氏连续。
1. 齐次线性。构造初值所处的线性空间到它的线性同构。齐次线性微分方程组总有零解（当且仅当任意一点初值为0）。能否非平凡表出恒为0的函数。等价。
2. $通解=\sum C_i{\varphi }_i(x)$（此处${\varphi }_i$均是$\vec{n}$维（列）向量值函数）。可以定义基解矩阵$\Phi =({\varphi }_1,\cdots ,{\varphi }_n)$，则$通解=\Phi (x)\vec{C}$，$\vec{C}$是列（常值）向量。
3. $W(x)$就是基解矩阵在某$x$处行列式。$W^{\prime }=tr(A(x))W$. 提示：行列式的多重线性。（注：可以用一维情形记忆）
4. $W$只可能恒为0或恒不为0，且恒为0当且仅当线性相关。（注意根据1.的线性同构，可以只考察$W(x_0)$就能得到相关性）
5. $A\Phi$（是“矩阵值函数”）
6. $\Phi B$仍是基解矩阵，而$B\Phi$导数为$BA\Phi$，所以是齐次线性微分方程组$\dot{\vec{y}}=BA{B}^{-1}\vec{y}$的基解矩阵。
7. $\Phi (x){\int }_{x_0}^x{\Phi }^{-1}(s)f(s)ds$，$\Phi$就是一维的$e^{-\int pdx}$（注意反号），$f$就是$q$. 注意此处$f$是列向量，$\Phi$恒非奇异。
   高维情况不好算（不能像一维一样显式表达）。
8. 齐次线性，线性无关。（注：有趣的是，$x=0$处，朗斯基行列式可能等于0. 但注意$x=0$不在定义域）
9. 线性无关（即不能非平凡组合出恒为0的函数）但朗斯基行列式可能为0.（实际上是恒为0）
10. 朗斯基，齐次线性（对于数值的，不是对于函数的），$A_1-A_2$（这里设我们考察的两组齐次线性微分方程组分别是$\dot{\vec{y}}=A_i\vec{y}(i=1,2)$）
11. 方法还是代入$\vec{y}=\Phi (x)\vec{C}$（千万注意求导的时候不要忘记$\vec{C}$也要求）。