常微分方程机敏问答[2] #20210619

---

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

一阶线性ODE

0. 一阶齐次线性ODE的解如果有一点等于0，则（）（注：可以通过解的唯一性和通解的表达式说明）。这与非齐次方程中解的唯一性有何联系？请验证$Dy:=dy/dx+p(x)y$中的$D$是线性算子，即$D(ay+bz)=aDy+bDz$.
1. 若$y^{\prime }+p(x)y=q(x)$中$p,q$都是周期为$T$的函数，已知非齐次方程解$f(x)$，那么$f(x+T)$也是解。这个结论和自治系统有何联系和区别？对$f(0)=f(T)$的情况，有什么结论？（提示：现在知道了两个非齐次方程的解，你想做什么？）
2. 若1.中$q(x)=0$，请直接由通解的形式得到方程非零解以$T$为周期的充要条件。
3. 利用2.，当$q(x)$不恒为0时，如果有两个不同的$T$周期解，则有什么结论？
4. 由于通解公式中出现$e^{-\int pdx}$，所以把$tanx,1/(1-x^2)$等函数作为$p$就成了出题的合理选择。请背诵上述两者的不定积分。
5. 将$y^{\prime }=1/cosy+xtany$换元化为一阶线性ODE. $3x{y}^2{y}^{\prime }+y^3+x^3$除了直接用$u=y^3,u^{\prime }=3y^2{y}^{\prime }$换元外，还可以化成什么样的一阶线性ODE？（写出一种）
6. 简述一阶线性ODE的积分因子解法怎么推广到不等式的证明。
7. $f$在$\mathbb{R}$上有界时，$e^{-x}{\int }_{-\infty }^{x}f(s)e^{s}ds={\int }_0^{e^x}f(lny)dy/e^x$有界。为了考察这样的$f$对应的$y^{\prime }+y=f$的所有可能的有界解，我们可以把一阶线性ODE的通解公式改写为$e^{-{\int }_{x_1}^{x}p(s)ds}({\int }_{x_0}^{x}q(s)e^{{\int }_{x_1}^{s}p(t)dt}ds)$. 试解释之前的常数$C$体现在哪里了？

答案

0. 该解恒为0. 注意两个非齐次方程的解相减等于一个齐次方程的解。
1. 自治系统$T$可以为任何实数。提示：相减，得到$f(0)=f(T)$时$f$为周期为$T$的周期函数。
2. $p$在一个周期平均值为0.
3. 一系列结论，比如：非齐次方程任意解都是周期解。（因为非齐次方程通解=齐次方程通解+非齐次方程特解）
4. $-ln|cosx|+C,\frac{1}{2}ln|\frac{1+x}{1-x}|+C$
5. $(siny{)}^{\prime }=cosy{y}^{\prime }$提示我们两边乘$cosy$凑出$(siny{)}^{\prime }=xsiny+1$. 这是一阶线性ODE.
   提示：对于$u^{\prime }+u/x+x^2=0$直接解出$u=\frac{C}{x}-\int x^{3}dx/x=\frac{C}{x}-\frac{x^3}{4}$，从而$v=x{y}^3+\frac{x^4}{4}$时我们有$v^{\prime }=0$是一阶线性ODE.（当然，这个角度有点刁钻）
6. 提示：凑全微分之后，以前是某表达式求导直接等于某确定的关于$x$的函数。如果改为某表达式求导大于等于0，则该表达式递增，可用于证明不等式。
7. $x_0$是待定常数（当然，$x_0$取遍一切值时不一定得到所有的解）。注意$x_1$没有本质影响。（思考：为什么？）

换元法

0. 提取出$y^{\prime }=f(x+y)$换元的核心思想（什么样的系统？）。并由此求解$cosy{y}^{\prime }+cosx=(siny+sinx+x{)}^2$
1. $y^{\prime }=f(y/x)$和0.有什么异同？
2. 什么叫关于$x,y$的$m$次齐次函数？为什么两个关于$x,y$的$m$次齐次函数之商能化为$f(y/x)$？
3. 换元$u=ax+b,v=cy+d,a\ne 0,c\ne 0$时$\frac{dv/du}{dy/dx}=$（），借此严格地求解$y^{\prime }=\frac{2y-x+5}{2x-y-4}$.（不需显式表达$y$）
4. 刚才3.中你解出的解曲线没有显式表达$y$，请解释你如何保证线素场竖直的点被挖去。
5. 换元$x=A{u}^B,y=C{v}^D$，可以给$y^{\prime }+p(x)y=q(x)$造成什么效果？给$y^{\prime }+y^2=x^m$呢？给$y^{\prime }=\frac{x^3+x{y}^2+x}{x^{2}y+y^3+y}$呢？借此严格求解$(x^2+y^2+3)y^{\prime }=2x(2y-x^2/y)$.
6. 回忆抛物线的光学性质，并解微分方程证明满足该光学性质的曲线是抛物线。

答案

0. 一维自治系统$u^{\prime }=f(u)$.
1. 都是化成$f(g(x,y))$再换元$u=g(x,y)$，但$u=x+y$时，$u^{\prime }=1+y^{\prime }=1+f(u)$直接只和$u$有关了。$y/x$求导后虽然不只显含$u$，但是分子和分母都只显含一个变量。
   另外注意对于$u=y/x$，可能需要单独考虑$x=0$情况。
2. $P(tx,ty)=t^{m}P(x,y)$.
3. $c/a$. $u=x-1,v=y+2,v_u^{\prime }=\frac{2v-u}{2u-v}$
   $r=v/u,r_u^{\prime }=\frac{v^{\prime }u-v}{u^2}=\frac{f(r)-r}{u}=\frac{\frac{2r-1}{2-r}-r}{u}$
   注：原始的方程等价于：在$u\ne 0$时$dr/du=\frac{r^2-1}{(2-r)u}$，且在$u=0$时满足原始方程。
   因上述注，我们首先考虑曲线是原方程解的必要条件。
   若解曲线上存在任何$r=1,-1$的点，则根据唯一性定理，这样的解只能是$y=x-3(x\ne 1)$或$y=-x-1(x\ne 1)$.
   否则$du/u=\frac{(2-r)dr}{(r^2-1)}=\frac{(1.5-1.5r+0.5+0.5r)dr}{(r+1)(r-1)}=\frac{-1.5dr}{r+1}+\frac{0.5dr}{r-1},r\ne 2,1,-1$
   $ln|u|=-1.5ln|r+1|+0.5ln|r-1|+C$
   $|x-1|=A|\frac{y+2}{x-1}+1{|}^{-1.5}|\frac{y+2}{x-1}-1{|}^{0.5}$
   $1=A|y+x+1{|}^{-1.5}|y-x+3{|}^{0.5},2x-y-4\ne 0,y+x+1\ne 0,y-x+3\ne 0,A>0$
   当然，以上可以总结成$(y+x+1{)}^3=C(y-x+3),x\ne 1,2x-y-4\ne 0,C\in \mathbb{R}$或$y-x+3=0,x\ne 1,2x-y-4\ne 0$.
   注：$C=0$的情况就是$y+x+1=0$.
   根据$x=1$处的连续性，解曲线是原方程解的一个必要条件是$(y+x+1{)}^3=C(y-x+3),C\in \mathbb{R},2x-y-4\ne 0$或$y-x+3=0,2x-y-4\ne 0$. 而这也是充分条件，这就严格进行了求解。
4. 始终保证右侧分母$2x-y-4$不为0即可。
5. $\frac{CD{v}^{D-1}{v}^{\prime }}{AB{u}^{B-1}}+p^{\ast }(u)C{v}^D=q^{\ast }(u),u^{1-B}{v}^{\prime }+p^{\ast \ast }(u)v=q^{\ast \ast }(u)v^{1-D}$.
   $\frac{CD{v}^{D-1}{v}^{\prime }}{AB{u}^{B-1}}+C^2{v}^{2D}=A^m{u}^{Bm}$，考虑$1-B=Bm,D-1=2D$则$K_1{v}^{\prime }+K_2{u}^{B-1}=K_3{v}^{-2D},K^1{v}^{\prime }+K_2{u}^{-m/(m+1)}=K_3{v}^2$.
   略。
   $y\ne 0$，换元$u=x^2+2,v=y^2+1$得到$(u+v)dv=2(2v-u)du$.
   令$t=u/v$，则$v\ne 0$时$t_v^{\prime }=\frac{\frac{t+1}{4-2t}-t}{v}$，下略。
6. $2arctan{y}^{\prime }=arctan(-x/y),-x/y=2y^{\prime }/(1-y^{\prime 2}),y^{\prime 2}x/y-2y^{\prime }-x/y=0$
   $$y^{\prime }=\frac{1\pm \sqrt{1+x^2/y^2}}{x/y}=y/x\pm \sqrt{y^2/x^2+1}:=f_{\pm }(y/x)=f_{\pm }(u)$$
   $$u_x^{\prime }=\frac{y^{\prime }-u}{x},dx/x=du/(\pm \sqrt{u^2+1}),ln|x|=\pm ln(u+\sqrt{u^2+1})+C$$
   $$x=A(\sqrt{u^2+1}\pm u),(x\pm Ay/x{)}^2=A^2(y^2/x^2+1),x^2\pm 2Ay=A^2$$
   注：根据实际意义，可能需要决定正负号，此处略。

积分因子法

0. 对于$x{y}^{\prime }+ay$结构，两项的$x$次数差一提示我们考虑形如（）的积分因子。求解出待定系数即可。
1. 直接写出$Pdx+Qdy=0$恰当的充要条件和$\mu Pdx+\mu Qdy=0$恰当的充要条件。若代入$\mu =\mu (x)$（即$\mu$与$y$无关）则如何？
2. 把上面的$\mu (x)$改为$\mu (\varphi (x,y))$（复合函数）则如何？
3. $\mu (Pdx+Qdy)=dF,\lambda (Rdx+Sdy)=dG$，那么怎么处理$(P+R)dx+(Q+S)dy$？
4. $P,Q$是齐$n$次函数则$Pdx+Qdy=0$有积分因子（）。解释证明过程用到的$n{t}^{n-1}P(x,y)=(t^n{)}^{\prime }P(x,y)=dP(tx,ty)/dt=x{\partial }_{x}P{|}_{(tx,ty)}+y{\partial }_{y}P{|}_{(tx,ty)}=t^{n-1}(x{\partial }_{x}P{|}_{(x,y)}+y{\partial }_{y}P{|}_{(x,y)})$。（下标表示求偏导）
5. 假设ODE的一族解是$F(x,y)=C$，若$C$表示”高度“，此时解曲线的直观意义是什么？

答案

0. $x^b$
1. $P_y=Q_x,(\mu P{)}_y=(\mu Q{)}_x$（下标表示偏导），$Q{\mu }_x=\mu (P_y-Q_x),(ln|\mu |{)}_x={\mu }_x/\mu =(P_y-Q_x)/Q$，如果右侧只和$x$有关就找到了积分因子$\mu$.
2. 略。
3. $\mu$是积分因子则$\Phi (F)\mu$根据链式法则也是。所以考虑$\Phi ,\Psi$使得$\mu \Phi (F)=\lambda \Psi (G)$.
4. $(xP+yQ{)}^{-1}$. 第二个等号：直接考虑$P((t+\Delta t)x,(t+\Delta t)y)=(t+\Delta t{)}^{n}P(x,y)$. 第三个等号：多元微分链式法则。第四个等号：直接根据定义验证。
5. 等高线。（等高，即$dF=0$. 因此直观上我们相信有积分因子）

应用举例

0. 回忆由单参数曲线族$\Phi (x,y;C)=0$如何决定对称形式的微分方程$Pdx+Qdy=0$，以及原始形式的微分方程$y^{\prime }=H(x,y)$，并请由此直接写出“正交轨线”满足什么微分方程。
1. $2ydx-xdy=0$的正交轨线族满足什么对称形式的微分方程？
2. 求与$xy=C$相交始终成$\pi /4$的曲线族。
3. 已知$x^{\prime }=x(-1+y),y^{\prime }=y(1-x)$的解曲线是某个第一象限点附近的一族闭轨，利用$(lnx{)}^{\prime }$和$(lny{)}^{\prime }$表达式说明一个周期内$y$的平均值$\bar{y}$对一切闭轨为定值。
4. 请形式上把$y^{\prime }=y/(x+y^2)$化为$x^{\prime }=f(x,y)$.
   平面上，$P$从原点出发，以匀速1向右运动，$Q$从$(0,1)$出发以匀速2追赶$P$，速度方向永远指向$P$，则考察$Q$相对于$P$坐标$(x,y)$，有$(x_t^{\prime },y_t^{\prime })=$（），现请利用前面的提示，考察$x_y^{\prime }$，并利用积分因子法一节的0.进行求解。

答案

0. 求导后联立消$C$. $y^{\prime }=-1/H(x,y)$.
   注：如果使用对称形式，则无需考虑分母为0的特例。
1. $xdx+2ydy=0$（同心椭圆）
2. 对$xdy+ydx=0$，形式上$y^{\prime }=-y/x$，则旋转$\pi /4$后$y^{\prime }=tan(arctan{y}^{\prime }\pm \pi /4)=\frac{-y\pm x}{x\pm y}$，$(x\pm y)dy+(y\mp x)dx=0$
   对于第一种情况，$d(xy+\frac{y^2}{2}-\frac{x^2}{2})=0$，第二种情况略。
3. 一个周期中，$0=\Delta (lnx)=\int x^{\prime }/x=\int (-1+y)=-T+T\bar{y}$
4. $x^{\prime }=x/y+y$（为一阶线性ODE）
   $(-1-\frac{2x}{\sqrt{x^2+y^2}},-\frac{2y}{\sqrt{x^2+y^2}})$
   $x_y^{\prime }=\frac{-r-2x}{-2y},-r-2x+2y{x}_y^{\prime }=0,(x/y{)}_y^{\prime }{y}^2=y{x}_y^{\prime }-x=r/2=\frac{\sqrt{y^2(1+(x/y{)}^2)}}{2}$，此时$x/y$可以换元为$u$.