常微分方程机敏问答[1] #20210611

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

解的含义、通解

0. 针对$y=\varphi (x,C_1,C_2)=C_{1}cosx+C_{2}sinx$，计算雅可比行列式$\frac{D[\varphi ,{\varphi }^{\prime }]}{D[C_1,C_2]}$. 注：$x$是一族函数的自变量，$C_i$是参数。
1. 对于$F(x,y,y^{\prime },y^{\prime \prime })=0$和它的通解$y=\varphi (x,A,B)$，写出全导数$dF/dA$的表达式，据此推出矩阵$(\begin{array}{ccc}\partial \varphi /\partial A& \partial {\varphi }^{\prime }/\partial A& \partial {\varphi }^{\prime \prime }/\partial A\\ \partial \varphi /\partial B& \partial {\varphi }^{\prime }/\partial B& \partial {\varphi }^{\prime \prime }/\partial B\end{array})$
   的第三列在任意满足一定非退化条件的$x,A,B$处，都可以被前两列线性表出。此处线性表出的系数一般来说是关于$x,A,B$恒定的吗？
2. $n$阶微分方程的解$y=\varphi (x,C_1,\cdots ,C_n)$中$n$个常数不独立的含义是什么为0？这个为0是恒为0吗？请举例说明你的判断。将你的判断和1.作联系，提取出一定的相同之处。
3. 平面$xOy$上过一定点的所有圆构成的曲线族有几个自由度？该曲线族至少是几阶微分方程的通解？
4. 对于$(x-a{)}^2+(y-b{)}^2=a^2+b^2$，把$y$看成$x$的函数，求一阶导后结果是什么？由原式和求一阶导后式子共两式，能得到$a,b$关于谁的表达式？把这个表达式代入什么式子可以得到不出现$a,b$的几阶微分方程？
5. 在$P:x=\xi ,C_1=c_1,\cdots ,C_n=c_n$邻域考察$F(x,y,\cdots ,y^{(n-1)})$的通解$y=\varphi (x,C_1,\cdots ,C_n)$。设${\eta }_k={\varphi }^{(k)}(\xi ,c_1,\cdots ,c_n)$，$Q$点表示$({\eta }_0,\cdots ,{\eta }_{n-1})$。
   现在为了得到以在$Q$点邻域中任取的某点$({\eta }_0+\Delta {\eta }_0,\cdots ,{\eta }_{n-1}+\Delta {\eta }_{n-1})$作为初值的满足$F(x,y,\cdots ,y^{(n-1)})=0$的特解，使用反函数存在定理，我们对哪个行列式有什么要求？
6. 反函数存在定理针对的是$n$元$n$维函数（即定义域和值域都是${\mathbb{R}}^n$），为什么上一问中不是这样？
7. 写出一个通解是$y=\{\begin{array}{c}-(x-A{)}^2,-\infty <x<A\\ 0,A\le x\le B\\ (x-B{)}^2,B<x<+\infty \end{array}$的一阶ODE，解释为什么一阶ODE出现了两个待定常数。
8. 写出一个解集中包含所有$y=arctan(x+C)$的一阶ODE，并找出它的所有解。
9. $y=A{x}^B,x>0,y>0$是什么二阶ODE的通解？
10. 设$y=g(x,C_1,\cdots ,C_n)$是一个充分光滑的函数族，$x$是自变量，$C_1,\cdots ,C_n$是$n$个独立参数。回忆并叙述独立参数的含义。借此独立性，我们知道利用一个小区域内的反函数存在定理，可得到$C_0,\cdots ,C_n$关于初值$y(x_0),\cdots ,y^{(n-1)}(x_0)$具有以$x_0$为参数的函数关系（用含参数的反函数组$c_{1,x_0}(\cdot ),\cdots ,c_{n,x_0}(\cdot )$表示）。
    请解释”含参数的反函数组“，借此说明存在一个$n$阶ODE，通解恰好是$g$.

答案

0. $\varphi =C_{1}cosx+C_{2}sinx,{\varphi }^{\prime }=-C_{1}sinx+C_{2}cosx,\partial \varphi /\partial C_1=cosx,\partial \varphi /\partial C_2=sinx,\cdots$，计算得结果为1.
1. $dF/dA=\frac{\partial F}{\partial y}\frac{\partial \varphi }{\partial A}+\frac{\partial F}{\partial y^{\prime }}\frac{\partial {\varphi }^{\prime }}{\partial A}+\frac{\partial F}{\partial y^{\prime \prime }}\frac{\partial {\varphi }^{\prime \prime }}{\partial A}$，把上式中A换为B又得到一个式子。再令前二式都为0，就得到待求的线性表出关系。”非退化“意指$\frac{\partial F}{\partial y^{\prime \prime }}\ne 0$. 不是恒定。
2. 雅可比行列式$\frac{D[\varphi ,\cdots ,{\varphi }^{(n-1)}]}{D[C_1,\cdots ,C_n]}$. 是恒为0. 例如考察$y=sin(x+C)$是$y^2+y^{\prime 2}=1$的通解，此处一个参数$C$显然应当是“独立”的。雅可比行列式$cos(x+C)$有时为0但不恒为0.
   联系前一问：参数不独立等价于雅可比行列式恒为0，或说雅可比矩阵最后一列恒能被前$n-1$列线性表出。这和1.中的$2\times 3$矩阵最后一列恒能被前2列线性表出意义类似。
   注：雅可比矩阵最后一列恒能被前$n-1$列线性表出，但对于不同的$x,C_1,\cdots$，系数可以不同。这和1.也是类似的。这种不一致性并不影响”恒为0“的判断。（这里的系数与$x$有关要区分于行列式值与$x$有关）
3. 两个。二阶。
4. $2(x-a)+2(y-b)y^{\prime }=0$. $x,y,y^{\prime }$. 原式求二阶导得到的$2-2b{y}^{\prime \prime }+2y^{\prime 2}+2y{y}^{\prime \prime }=0$. 二阶。
5. $P$点处的雅可比行列式$\frac{D[\varphi ,\cdots ,{\varphi }^{(n-1)}]}{D[C_1,\cdots ,C_n]}{|}_{x=\xi ,C_i=c_i}$.
   注：$P$点一点处的雅可比行列式不为0可以推出其邻域的该行列式恒不为0.
6. 提示：$x$是”参数“。这里针对每个固定的$x$值例如$x_0$，用的都是$n$元$n$维反函数存在定理，即已知各参数$C_i$在一点$x_0$处决定各阶导$y^{(i)}$后，反求如何由各阶导$y^{(i)}$决定该点处应当取的各参数$C_i$.
7. $y^{\prime }=2\sqrt{|y|}$（注意此处正负号问题，特别是$x-C_1<0$）. 一阶ODE只能有一个独立的待定常数。但，这里两个待定常数考察的甚至都不是一个区间，更没法求雅可比行列式，更无法定义是否独立。
   实际上这里的待定常数反映的是解的不唯一性和“奇解”现象。
8. 求导，得$y^{\prime }=(1+(x+C{)}^2{)}^{-1}$，代入$x+C=tany$，即得到一个可能的方程$y^{\prime }=co{s}^{2}y$.
   找前述方程所有解：如果$cosy=0$则$y=k\pi +\pi /2(k\in \mathbb{Z})$，否则$dy/co{s}^{2}y=dx,tany=x+C,y=arctan(x+C)+k\pi (k\in \mathbb{Z})$.
9. （求二阶导后消去常数，得）$xy{y}^{\prime \prime }=y^{\prime }(x{y}^{\prime }-y)$.
10. 提示：$y-g(x,c_{1,x}(y,y^{\prime },\cdots ,y^{(n-1)}),\cdots ,c_{n,x}(y,y^{\prime },\cdots ,y^{(n-1)}))$从形式上可以看成一个$x,y,y^{\prime },\cdots ,y^{(n-1)}$的函数。

线素场

0. 对于$dy/dx=P(x,y)/Q(x,y)$中，$P\ne 0,Q=0$的点，如何作线素合理？如何用一个方程统一所有$P^2+Q^2>0$的情况？前述方程并不符合ODE原始定义的形式$F(x,y,\cdots ,y^{(n-1)})=0$，那它是什么意思呢？
1. 线素场在$P=0,Q=0$（奇异点）处是否有意义？奇异点在奇异点集中一定孤立嘛？给出证明或反例。
2. 利用线素场研究$y^{\prime }=1+xy$，说出解曲线（这里指“尽量延伸”后的解曲线，即定义域取得极大的区间，下同）的增减性可能有几种情况。
3. 利用线素场研究$y^{\prime }=x^2-y^2$。标记$y^{\prime }=0,y^{\prime }>0,y^{\prime }<0$分别对应的区域。
4. 在3.中记$y$正半轴，$x$负半轴，$y$负半轴，$x$正半轴所在区域分别为$A,B,C,D$. 判断经过且仅经过以下区域组合的解曲线是否可能存在，并说明理由。$CD,ABCD$.
5. 在3.中仅从$y^{\prime }$符号判断，经过$BCD$和经过$BC$的解曲线都是可能出现的。实际上它们确实可能出现。用比较定理说明$x=0,y=-100$的点对应的解曲线经过$BC$，而$x=0,y=-0.01$的点对应的解曲线经过$BCD$。用确界存在定理以此为基础给出一个优美的结论。

答案

0. 垂直于$x$轴。$Qdy-Pdx=0$（“对称形式”）. 提示：在线素场竖直的点附近，满足另一个方程。
1. 无意义。不一定。略。
2. （画线素场，找到单调性成立的区域）3种。即先减后增、先增后减、增。
3. 略。
4. 不可能。$CD$：注意$C$区域，函数递减，则函数向左延伸一定经过$B$区域。$ABCD$：经过$C$则容易说明只可能是$BCD$或$BC$，不可能经过$A$.（证明要点是考察“越过边界”时的斜率是否引出矛盾）
5. 对于$x=0,y=-100$，每次让$x$变大1，在长度为1的区间上用比较定理，把$x$和$y$都放缩成定值，并归纳。
   对于$x=0,y=-0.01$，在$x\in [0,1]$区间用比较定理比较$y^{\prime }=-y^2$的同初值的解即可说明解曲线与$C$和$D$的边界线有交点。
   注意到：若$x=0,y=a<0$对应的解曲线不经过$D$，则$x=0,y=a_1<a$对应的解曲线也不经过$D$.（这两条解曲线初值不同，所以应当平移后用比较定理）
   类似地考察经过$D$的解曲线，从而根据确界存在定理存在一个初值的“分界点”$a^{\ast }$，当初值$x=0,y>a^{\ast }$时解曲线经过$D$，否则不经过。

恰当方程和分离变量法

0. 二阶偏导数在满足一定条件下可交换和“恰当方程”的判断有何联系？
1. 利用分组凑全微分和公式法两种方法分别解$(y{e}^x+2e^x+y^2)dx+(e^x+2xy)dy$。在分组凑全微分时，没有”配对“的$2e^x$具有什么特点使得其可以不“配对”？
2. 小心地求解$\frac{2s-1}{t}ds+\frac{s-s^2}{t^2}dt=0$.
3. 解恰当方程时，非单连通区域可能带来什么问题？举例说明。
4. 写出一个以$x=0$（严格来说，是$x=0,y\ne 0$，下同）为一个特解，$y=0$为另一个特解，且原点是奇异点（无法定义线素场方向）的ODE. （写上节0.中的“对称形式”，下同）
5. 再写出一个满足4.条件且至少在一个非奇异点处解不唯一的ODE.
6. 一维自治系统$y^{\prime }=f(y)$中$f(y)$有唯一零点0，讨论${\int }_0^1\frac{dy}{f(y)}$的敛散性和系统导函数图像，一维相图，线素场，解唯一性的关系。
7. 甲从原点出发沿$x$轴正方向作速度为1的匀速直线运动，乙从$(0,1)$出发，运动方向始终指向甲且和甲保持距离恒定。写出乙轨迹曲线上$y^{\prime }$的表达式并求乙运动轨迹（结果可以表示成$x(y)$形式）。
8. 求出乙在任意时刻的速度。

答案

0. 提示：判断标准是对于$Pdx+Qdy=0$，看是否有${\partial }_{y}P={\partial }_{x}Q$. 这里$P$可以看成某个二元函数关于$x$的偏导。
1. 略。提示：$2e^{x}dx=d(2e^x)$，不涉及$y$. （或者说$\frac{\partial 2e^x}{\partial y}=0$）
2. $\frac{2s-1}{s-s^2}ds+dt/t=0$（此处特解$s=0,t\ne 0$或$s=1,t\ne 0$），若非特解，则
   $(1/(1-s)-1/s)ds+dt/t=d(-ln|s(1-s)|+ln|t|)=0$
   解得$t/s(1-s)=C,C\ne 0$（注意此时$t\ne 0,s(1-s)\ne 0$自然满足）
3. 多值性。例如$\frac{ydx-xdy}{x^2+y^2}=0,x^2+y^2>0$.
4. 如$ydx-xdy=0$.
5. 如$P(x,y)dx+Q(x,y)dy=0$，其中
   $P(x,y)=\{\begin{array}{c}\sqrt{|y|},|x|>2\\ 略（使得光滑连接）,1\le |x|\le 2\\ y,|x|<1\end{array}$
   $Q(x,y)=\{\begin{array}{c}-1,|x|>2\\ 略（使得光滑连接）,1\le |x|\le 2\\ -x,|x|<1\end{array}$
6. 提示：积分发散等价于系统解唯一，一维相图上运动到平衡点的“时间”需要无穷大。注意：导函数图像在$0$附近和$x$轴不正交（成一定角度，相当于李氏连续性）是解唯一的充分不必要条件。（思考：请举例说明不必要）
7. $y^{\prime }=-y/\sqrt{1-y^2},x=\frac{1}{2}ln\frac{1+\sqrt{1-y^2}}{1-\sqrt{1-y^2}}-\sqrt{1-y^2}$
8. $t$时刻，甲处于$(t,0)$，乙处于$(x(y),y)$且$y^{\prime }=-y/\sqrt{1-y^2}=-y/(t-x)$，得
   $$t=\frac{1}{2}ln\frac{1+\sqrt{1-y^2}}{1-\sqrt{1-y^2}},y=\sqrt{1-(\frac{e^{2t}-1}{e^{2t}+1}{)}^2}$$
   $$v=|\sqrt{1+y^{\prime -2}}\frac{dy}{dt}|=略$$

一阶线性ODE：基本知识

0. 直接说出$x^{\prime }=ax$的通解。直接说出$x^{\prime }=b(x-x_{final}),b<0,x(0)=x_{init}$的解并说出物理意义。
1. 直接说出$x^{\prime }=x+cost$保留变限积分的通解。
2. 用欧拉公式和考察对称式两种方式求不定积分$\int e^{-t}costdt$
3. 说出一阶线性ODE一般形式。
4. 一句话概括积分因子法。
5. 一句话概括常数变易法。背诵一阶线性ODE通解的公式，指出变易体现在哪里。

答案

0. $x=A{e}^{at}$（注意$A$可以为任意实数）。$x=x_{final}+(x_{init}-x_{final})e^{bt}$，“指数趋近”。
1. $x=(C+{\int }_{t_0}^t{e}^{-u}cosudu)e^t$。
2. 欧拉公式：${\int }_a^b{e}^{-t}costdt={\int }_a^{b}Re(e^{(-1+i)t})dt$，$a,b$为实数。设$L_{ab}$表示复平面上连接$a,b$两点的线段，于是
   $$\begin{gathered} 原式=Re({\int }_{L_{ab}}{e}^{(-1+i)z}dz)=Re(\frac{e^{(-1+i)z}}{-1+i}{|}_a^b)=Re(\frac{e^{(-1+i)b}-e^{(-1+i)a}}{-1+i}) \\ =(令a=0)Re(\frac{(-1-i)(e^{-b}cosb-1+i{e}^{-b}sinb)}{2})=\frac{-e^{-b}cosb+1+e^{-b}sinb}{2} \end{gathered}$$
   则不定积分结果为$\frac{-cosb+sinb}{2}{e}^{-b}+C$.
   考察对称式：分别对$e^{-t}cost$和$e^{-t}sint$求导，再解方程。
3. $y^{\prime }+p(x)y=q(x)$.
   （注：这里$py$项如果书写在右边也很合理，且容易推广到高维情况）
4. $y^{\prime }+py=q$的两边同乘积分因子$e^{\int pdx}$（严格来说是$e^{{\int }_{x_0}^{x}p(s)ds}$，使用“哑元”或说局部变量$s$）后左边变成全微分。
5. 齐次（$q(x)=0$）通解$y=C{e}^{-\int pdx}$中设$C=C(x)$，代回原方程解出$C(x)$.
   $$y=e^{-\int pdx}(C+\int q{e}^{\int pdx}dx)(注：此处没有用“哑元”书写)$$
   变易体现在括号内多了一项。