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

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

存在唯一性定理

0. 对于没有解析解的微分方程，我们可能有什么考察的角度？
1. 叙述皮卡定理中“首先达到左右边界”和“首先达到上下边界”分别对应什么样的解存在闭区间。李氏条件用在了哪里？
2. 叙述与初值问题$y^{\prime }=f(x,y),y(0)=y_0$等价的积分方程。怎么把该方程改写成皮卡序列迭代所用到的公式？
3. 迭代有意义前提是不会超出上下边界。所以需要证明$y_n(x)-y_0$关于$|x-x_0|$有什么结论？
4. 此处证明函数项级数一致收敛用的是M判别法，其中考察的收敛的数项级数是什么？反设有两个解$u(x),v(x)$，证明唯一性时，考察的数项级数和前面的有何异同？
5. 证明皮卡定理时，一致收敛用在了什么地方？
6. 说明李氏条件是Osgood条件的特例（提示：${\int }_0^{\delta }dr/Lr=\infty$）。用李氏条件或Osgood条件证明唯一性的核心思想是考察什么的积分？
7. 不满足李氏条件（从而Osgood条件）下解一定不唯一嘛？皮卡序列一定收敛吗？
8. 利用考察6.中积分类似的方法，还能证明$f$关于$y$单调时，解具有什么样的唯一性？
9. 回忆解的存在唯一性定理证明过程，如何估计$y^{\prime }=x^2-y^2,y(-1)=0,|x+1|\le 1,|y|\le 1$的解的存在区间？
10. 用类似于皮卡定理证明过程的归纳法，可以得到皮卡序列和真解的一个差上界$|y_n-y^{\ast }|\le \frac{M}{L}\frac{L^{n+1}|x-x_0{|}^{n+1}}{(n+1)!}$。由此简要说出如何提高解的精度。

答案

0. 解的结构（线性空间？流形？），数值解，存在唯一性，几何直观等。
1. $[x_0-a,x_0+a],[x_0-b/M,x_0+b/M](M>ma{x}_{x,y}|f(x,y)|)$. 李氏常数用途：对迭代不同次数的结果之间（如$y_1$和$y_2$之间），估计对应的$x$处，两个$f(x,y)$的差。（注：若再积分，即可估计$y_1$与$y_2$的差）
   我们直接使用了一些课本上的记号，下同。
2. $\varphi (x)=y_0+{\int }_{x_0}^{x}f(x,\varphi (x))dx(x在一定范围内)$，不断“近似认为”右侧$\varphi =y_k$，就是迭代公式${\varphi }_{k+1}(x)=y_0+{\int }_{x_0}^{x}f(x,{\varphi }_k(x))dx$.
3. $|y_n(x)-y_0|\le M|x-x_0|$.
4. $\frac{M}{L}\frac{L^n|x-x_0{|}^n}{n!}$. 证明唯一性：初始是一个常数界$K$，表示两解之差上界。证明迭代收敛：初始是$M|x-x_0|$. （当然，证明的表达形式很多，言之有理即可）
5. 交换积分和极限顺序。
6. 略。都要考察$f(x,y_1(x))-f(x,y_2(x))$积分，从而得到两解之差（的变化量）。
7. 不一定，不一定。（如“米勒之例”）
8. 单侧。
9. 最简单的：直接考虑 m i n { a , b / M } = 1 / 4 min\{a,b/M\}=1/4 min{a,b/M}=1/4. 即得解得存在区间为$[-1.25,-0.75]$或更大。
   稍精确：右侧 m i n { a , b / M } = 1 min\{a,b/M\}=1 min{a,b/M}=1显然。考察左侧，为了使得$a=b/M$，列方程$h=1/(1+h{)}^2$，即$h^3+2h^2+h-1:=f(h)=0$.
   $f$在$h>0$递增。$f(0.5)>0,f(0.4)<0$，所以解的存在区间为$[-1.4,0]$或更大。（此处并不能断言比$[-1.5,0]$更小。为什么？）
10. 增加迭代次数。或减小区间，以获得更小的$M,L$。

压缩映射

0. 考察$y=x$和$y=f(x)$两曲线在$xOy$平面的图像，说出怎么迭代求$f$的根，并由此给出压缩映射原理的几何直观。
1. 用压缩映射证明皮卡定理，则距离度量和完备性分别是什么？
2. 接上，压缩映射证法得到的解存在区间相比上一节的证法更小，其原因是哪个常数可能较大？
3. $|f(y)-f(x)|<|y-x|$在有界闭区间上成立，则可以用确界存在定理考察$F(x):=|f(x)-x|$作为$x$的函数，最小值点能在闭区间上取得，且$F(x)$最小值为（）。如果改为考察无界区间，那么不动点还存在吗？
4. $g(x)=x^2-4$，则迭代公式$x_{n+1}:=x_n-g(x_n)/g^{\prime }(x_n)$（在一定区域内）为何是压缩映射？其收敛速度如何？

答案

0. 求根方法：交替画横线和竖线。几何直观：考察斜率。
1. “一致”的度量（“差的绝对值最大值”），一致收敛的柯西准则。
2. 李氏常数。
3. 0，不一定，如把$x-e^x+C$（$C$为常数）的一段和射线接起来，得到分段函数。
4. $G(x):=x-g/g^{\prime }$作为关于$x$的函数，在足够小的邻域内有$|G^{\prime }(x)|\le \alpha <1$.
   所以，这里其实用到了$g$一阶导非零，二阶导存在。
   对于本题，$g^{\prime }=2x,x_{n+1}=x_n/2+2/x_n,|x_n-2|=|(x_n-2{)}^2/2x_n|$.
   一定条件下$|(x_n-2{)}^2/2x_n|<|(x_n-2{)}^2|$，从而误差依次为平方，四次方，八次方……量级，即幂次关于$n$以指数增大。收敛速度很快。

解的存在性

0. 解释欧拉折线落在“角形区域”（图像：一个“叉”把矩形切割成4部分），并画出两种可能的简图。
1. 为什么一系列“欧拉折线”满足Ascoli引理的题设？
2. 怎么统一地表达出欧拉折线？（即：形式上表达出欧拉折线对应的函数${\varphi }_n(x)$）
3. 欧拉折线代入上一节2.中积分方程，怎么估计误差？
4. Ascoli引理推出的（）性质用在了哪里？
5. 一致有界和等度连续推出一致收敛成立的区间是且仅是（）区间。
6. 为了构造使得奇数段情况和偶数段情况分别递增和递减的欧拉折线，函数$xcos(\pi /x)$起到什么作用？
7. 如果$f(x,y)$不连续，解可能不存在，举出例子。
8. 皮卡序列是否满足Ascoli引理题设？它能不能用来证明皮亚诺存在定理？
9. 由物理中伯努利方程，桶开小洞漏水时$\dot{h}=-C\sqrt{h}$. 解不唯一表示什么情况下无法反推之前情况？
10. 对于一维自治系统$x^{\prime }=f(t,x)=f(x)$（只与$x$有关）和非自治系统，分别估计欧拉折线法走一步（导致折线和真实解曲线）的误差：其大致是什么常数乘以$\Delta t^2$（步长平方）？
11. 如果改进欧拉折线法：先根据$f(0,0)$走一步$\Delta t$得到点$(\Delta t,f(0,0)\Delta t):=(t_1,x_1)$，再进行一定修正，得到点$(\Delta t,(f(0,0)+f(t_1,x_1))\Delta t/2):=(t_2,x_2)$. 则对自治系统，$(t_2,x_2)$处折线对应的取值和真实解曲线对应的取值误差在什么量级？

答案

0. 注意导函数$f$有界。略。
1. 一致有界显然。由导函数$f$有界得到等度连续性。
2. ${\varphi }_n(x)=y_0+{\sum }_{k=0}^{s-1}f(x_k,y_k)h/n+f(x_s,y_s)(x-x_s)$.
3. 提示：利用$f$关于第二个变量（$y$）的一致连续性（这比李氏条件显然弱）。
4. 一致收敛，交换积分和极限。
5. 有限。
6. 保证原点处连续，且注意奇偶性影响。
7. 提示：一系列可数个同心菱形，半径依次缩小，一圈$f=1$，一圈$f=-1$.
8. 满足。不能。（回忆皮卡定理中需要两边取极限后，才能说明$\varphi$确实是解）
9. 桶漏光水之后。
10. 自治：$\frac{1}{2}f(0)f^{\prime }(0)$. 非自治：全导数$\frac{1}{2}df/dt=\frac{1}{2}({\partial }_{t}f+{\partial }_{x}f\cdot f){|}_{(0,0)}$.（原因：泰勒展开的二阶项）
11. $f(t_1,x_1)=f(x_1)\approx f(0)+f^{\prime }(0)x_1+f^{\prime \prime }(0)x_1^2/2,\frac{f(x_1)+f(0)}{2}\approx f(0)+f^{\prime }(0)f(0)\Delta t/2+f^{\prime \prime }(0)f^2(0)\Delta t^2/4$
    误差$\approx [(f^{\prime \prime }{f}^2+f^{\prime 2}f)/6-f^{\prime \prime }{f}^2/4]\Delta t^3$. （泰勒展开要到三阶）

分岔

0. 含参数微分方程的参数变化时解结构会改变。如果有多个参数，那么参数和解结构的关系图从原理上类似于（）（说出一种大学物理中学到的事物）。注：含多个参数的时候，如果系统参数只在某条线上变化，那又相当于退化到一维，仍然可能出现类似于单参数情况的分岔。
1. 分岔反映在解的几何意义和物理意义上，可能有突变、相变等。对$y^{\prime }=-ay-y^2(a\ge 0)$说明。
2. 用二维系统的鞍点、结点限制在某一直线上的情形解释“鞍-结点”分岔。
3. 分岔图的横轴是参数$r$，纵轴是（）。导函数图像横纵轴是（）。相图中空心点对应分岔图中（）。分岔图中曲线如何反映到相图上？
4. “鞍-结点”分岔与其它三种（transcritical, subcritical, supercritical分岔）在考察导函数图像时的几何意义显然不同（提示：平移与旋转），试直观解释。
5. 解释三种“临界”分岔中不动点的失稳，比较亚临界和超临界分岔的物理意义，画图解释磁滞（hysteresis）并猜想“磁滞”的词源。
6. 临界分岔$x^{\prime }=x^3-rx$如果右侧加“不完美项”常数$h$使得体系不再具有对称性，则临界分岔变为（）分岔。为了从几何直观考察，应该绘制什么样的三维图？
7. 接上，对于$x^{\prime }=rx-x^2+h$中的$h$呢？（提示：考虑$h$的符号）
8. 回忆0，考察7.（请特别考察两种越过界线的方向）
9. 为了考察$x^{\prime }=rx-sinx$带来的分岔现象，画出（）的图像再考察过原点不同斜率直线与其交点即可。该系统有一个（）分岔和许多（）分岔。

答案

0. （物态变化中的）相图。注意”相图“一词多义。
1. $y^{\prime }=-y^2$则$y=1/(x+C)$. 这种衰减比指数慢得多。
2. 平衡点相碰消失（或凭空产生）。（可以尝试$x^{\prime }=r-x^2,y^{\prime }=-y$）
3. $x$，纵轴$x^{\prime }$横轴$x$，虚线，略。
4. ”鞍-结点“是平移，后三种是旋转。
5. 同一个不动点随参数变化，稳定性改变。亚临界有“灾变”（灾难性后果），图略。回忆物理中铁磁性介质的磁滞回线。
6. 鞍-结点。例如$z(x,r)=x^3-rx$的三维图像，然后用不同高度截取。
7. $h<0$则：在$r=0$处，$x^{\prime }=0$无解，反映在分岔图上即$x=0$处没有曲线。此时是（两个）鞍-结点分岔。
   $h=0$则跨临界分岔。
   $h>0$则可以解出整个$r\in \mathbb{R}$上两条函数曲线$x_1(r),x_2(r)$，并不存在分岔。
8. 提示：全平面$rOh$分为两个区域。
   注意：固定$h$考察$r$的分岔，与固定$r$考察$h$的分岔，本质上是两个东西。
9. $sinx$，亚临界，鞍-结点。

分岔的应用

0. 阻尼力$-b\dot{\varphi }$（$\varphi$是角度）中$b$单位是什么？负号是什么意思？
1. “环”上运动的珠子中，怎么用无量纲化定义“阻尼很大”？
2. $\dot{x}=sinx(\gamma cosx-1)$的分岔怎么考察？
3. “环”上运动的珠子中，导函数图像体现在数轴上是怎样的？（提示：注意阐述周期性）
4. 请在二维相图上表示奇异极限。此时线素场如何？
5. 对于$ϵx^{\prime \prime }+x^{\prime }+x=0$的奇异极限，怎么在$x(t)$图像上体现？（参考4.）
6. 分岔的稳定性与物理中的势如何结合？（提示：构造势函数）

答案

0. $N\cdot s$. 阻碍运动（力与运动方向相反）
1. 略。
2. 泰勒展开到三阶（不要忘了$sinx$三阶项）
3. 提示：“环”结构中，数轴上无穷个点表示实际的一个点。导函数有周期性，故实际上只用在一个周期内考察和画图。注意分岔导致有几种可能的图。
4. 提示：先极快到一条曲线附近，然后跟着曲线走。除了一条曲线附近，其余处线素场近似竖直。解曲线上水平。（即：偏离就很快拉回来）
5. 导数很快从初值变成某个和初值无关的定值。（注：且这过程其实是”指数的“趋近）
6. $x^{\prime }$积分取相反数就是势函数。回忆物理中$力=-\nabla 势$.
   尝试画出势函数图，发现各类分岔的物理意义非常明显。

幂级数解法

0. 对于$y^{\prime }=x^2+y^2,y(0)=0$分别用皮卡序列和待定系数法求解的幂级数展开到$x^7$项（不用严格说明）
1. 什么叫函数在$G\in {\mathbb{R}}^2$内解析？
2. 解释求幂级数形式解中，$C_2=\frac{1}{2!}{y}^{\prime \prime }(x_0)=\frac{f_x^{\prime }(x_0,y_0)+f_y^{\prime }(x_0,y_0)y^{\prime }(x_0)}{2}=\frac{a_{10}+a_{01}{C}_1}{2}$是什么意思。
3. 进行形式运算：对于$y^{\prime }=ry-y^2+a{y}^3+O(y^4)$，若$y=x-b{x}^3$，则代换得（）。上式中对于$r\ne 0$可适当取（）使得$x^{\prime }$关于$x$表达式三阶项消失。

答案

0. 皮卡：$y_0=0,y_1=x^3/3,y_2=\int (x^6/9+x^2)dx=x^3/3+x^7/63$
   待定系数：$a_0=0,a_1=a_0^2=0,2a_2=2a_{0}2a_1=0,\cdots$这样计算下去。
1. 提示：展开式每一项是二元单项式。展开收敛，且等于原函数。
2. 提示：$C_i$是解的（一元）幂级数展开系数，$a_{ij}$是$f(x,y)$的（二元）幂级数展开系数。注意多元微分链式法则。
3. $x^{\prime }-3b{x}^2{x}^{\prime }=rx-rb{x}^3-x^2+a{x}^3+O(x^4),x^{\prime }=rx+rx3b{x}^2-rb{x}^3-x^2+a{x}^3+O(x^4)=rx-x^2+(a+2br)x^3+O(x^4)$.
   $b$.

奇解和包络

0. 一阶隐式微分方程相比显式情况，“线素场”怎么样了？
1. 对$y=x{y}^{\prime }-\frac{1}{4}{y}^{\prime 2}$两边求微分得到（）为0或（）为0，这是否是解的充要条件？
2. 证明上式中“抛物线”解每一点处切线都是解。
3. 参数法求隐式方程的解：请解释$dy=y^{\prime }dx\Rightarrow g_{u}du+g_{v}dv=h(u,v)(f_{u}du+f_{v}dv)$的各字母含义。
4. 对于一阶隐式微分方程$F(x,y,y^{\prime })=0$一个解$y=\varphi (x)$和通解$y=f(x;C)$，当（）不属于（），且（）解曲线上任一点出发都得到（）中的一个解，则（）是奇解。奇解唯一嘛？
5. $F(x,y,y^{\prime })$在区域内连续且对第二、第三个变量可偏导时，奇解上每一点$F$处对第几个变量偏导为0？证明思路：反设一点附近（），则由（）元、（）维情况隐函数存在定理，得到局部有（）元函数$p(x,y)$使得（）为0，这样形式上可以把隐式微分方程化为显式，然后根据（）（人名）定理推出矛盾。由上可以得到奇解的（）条件。
6. 为什么$y^{\prime }=\sqrt{y}$不能用5.考察？
7. 什么是包络？（一阶连续可微的）通解对应曲线族的包络线是奇解对应的解曲线如何直观解释？
8. 给出包络满足的“C-判别式”。为了证明它，我们可以形式上写出包络线上$x,y,C$用$C$表示的参数方程（），由$V(x,y,C)$恒为0，求全微分得（），再根据向量（）和向量$(V_x,V_y)$关系即得结果。
9. 8.的几何直观解释：设$C$表示高度，每个高度上有一条曲线$V(x,y,C_0)$，如果$V_C^{\prime }\ne 0$则局部有二元一维反函数（），设$V$一阶连续可微（且$V_x^{\prime },V_y^{\prime }$不全为0），则局部类似于”山坡的等高线“，不可能有包络。

答案

0. 不一定唯一（不一定良定义）
1. $y^{\prime \prime }=0$或$x-y^{\prime }/2=0$. 必要不充分。
2. 提示：由上，抛物线满足$y^{\prime }=2x$，直接写出切线方程，验证是解即可。
3. 把$x,y,y^{\prime }$都用两个参数$u,v$表示，分别为$x=f(u,v),y=g(u,v),y^{\prime }=h(u,v)$.
4. $\varphi (x)$， { f ( x ; C ) ∣ C 取 一 切 可 能 值 } \{f(x;C)|C取一切可能值\} {f(x;C)∣C取一切可能值}，$\varphi (x)$， { f ( x ; C ) ∣ C 取 一 切 可 能 值 } \{f(x;C)|C取一切可能值\} {f(x;C)∣C取一切可能值}，$\varphi (x)$. 不一定。
5. 第三个。${\partial }_{y^{\prime }}F\ne 0$，2，1，2，$F(x,y,p)$，皮卡，必要不充分。
6. 对第二个变量不可偏导，于是皮卡定理条件不满足。
7. 对曲线族$V(x,y,C)=0$和不在其中的曲线$U(x,y)$，$U$上任一点都和某条曲线族$V$中的曲线相切，则$U(x,y)$称为“包络线”。提示：相切表示导数相等，而点点都有切线也呼应了奇解的定义。
8. $V(x,y,C)=V_C^{\prime }(x,y,C)=0$. $\{\begin{array}{c}x=f(C)\\ y=g(C)\\ C=C\end{array}$，$V_x{f}_C+V_y{g}_C+V_C=0$，$(f_C,g_C)$.（下标表示偏导）
9. $C(x,y)$.