1. 常微分方程是什么

主要考虑微分方程的解法，而不考虑理论问题。在数学系需要有一个学期的课程进行专题学习。——柳斌

1.1. 数学史引入

1.1.1. Newton第二定律：

md2xdt2=F(t,x→,x˙→)m\frac{\mathrm d^2x}{\mathrm dt^2}=F(t,\overrightarrow{x},\overrightarrow{\dot{x}}) mdt2d2x​=F(t,x,x˙)
就是一个关于xxx的微分方程

1.1.2. 万有引力

利用Newton第二定律
地球所受的万有引力可以近似表示成：
Gmmme∣xm−xe∣2xm−xe∣xm−xe∣+Gmsme∣xs−xe∣2xs−xe∣xs−xe∣=med2xdt2G\frac{m_mm_e}{|x_m-x_e|^2}\frac{x_m-x_e}{|x_m-x_e|}+G\frac{m_sm_e}{|x_s-x_e|^2}\frac{x_s-x_e}{|x_s-x_e|}=m_e\frac{\mathrm d^2 x}{\mathrm dt^2} G∣xm​−xe​∣2mm​me​​∣xm​−xe​∣xm​−xe​​+G∣xs​−xe​∣2ms​me​​∣xs​−xe​∣xs​−xe​​=me​dt2d2x​

1.1.3. 常微分方程的贡献者

• Bernoulli家族
• Euler
• Lagarange
• Gauss

1.1.4. 一个突破

Liouville发现，y′=x2+y2y'=x^2+y^2y′=x2+y2这样一个右端很好的函数无初等函数解。
所以我们的微分方程求解转化到通过方程得出解的一些好性质的阶段。
我们的讨论仅限于之前的可计算的问题。

1.2. 定义

1.2.1. 微分方程

含有未知函数及其至nnn阶导数的方程
F(x,y,y′,⋯,y(n))（x∈R）(1)F(x,y,y',\cdots, y^{(n)})（x\in\R）(1) F(x,y,y′,⋯,y(n))（x∈R）(1)
称为一个nnn阶常微分方程。其中各个函数都在对应相同的xxx处取值。

阶数等于方程中最高阶导数的阶数。

例：

1. dydx=y(x−1)\frac{\mathrm dy}{\mathrm dx}=y(x-1)dxdy​=y(x−1)不是微分方程，因为其中出现了非xxx的变量（式中y(x−1)y(x-1)y(x−1)是yyy在x−1x-1x−1处的取值）
2. dydx=y(y(x))\frac{\mathrm dy}{\mathrm dx}=y(y(x))dxdy​=y(y(x))不是微分方程。因为函数关系不明确。

1.2.2. 方程的解

称y=φ(x),x∈(a,b)y=\varphi(x), x\in(a,b)y=φ(x),x∈(a,b)是(1)的解，如果φ\varphiφ有直到nnn阶导数，且F(x,φ(x),φ′(x),⋅φn(x))≡0,x∈(a,b)F(x,\varphi(x), \varphi'(x), \cdot\varphi^{n}(x))\equiv0, x\in(a, b)F(x,φ(x),φ′(x),⋅φn(x))≡0,x∈(a,b)。

说人话：将函数代入微分方程使得方程成为恒等式，这个函数就是这个方程的解。换言之，微分方程的求解，是利用函数和函数导数关系，求得一个函数作解。

解y=c1+c2t−12gt2,c1c2y=c_1+c_2t-\frac{1}{2}gt^2,c_1c_2y=c1​+c2​t−21​gt2,c1​c2​为常数。

1.2.3. 通解

称y=φ(x,c1,c2,⋯c,n)y=\varphi(x,c_1,c_2,\cdots c,_n)y=φ(x,c1​,c2​,⋯c,n​)为方程的通解，如果

1. y=φ(x,c1,c2⋯,cn)y=\varphi(x,c_1,c_2\cdots, c_n)y=φ(x,c1​,c2​⋯,cn​)是解。
2. c1,⋯,cnc_1,\cdots, c_nc1​,⋯,cn​是独立的常数（大致理解，这些常数通过积分产生）。

有时候通解不能通过显函数的方式得到，那么利用隐函数形式表达的通解叫做通积分。

1.2.4. 方程特解中常数独立性的判别

Jacobi矩阵
(∂φc1∂φc2⋯∂φ∂cn∂φ′∂c1∂φ′c2⋯∂φ′∂cn⋮⋮⋱⋮∂φ(n−1)c1∂φ(n−1)φc2⋯∂φ(n−1)∂cn)\left( \begin{matrix} \frac{\partial \varphi}{c_1} &\frac{\partial \varphi}{c_2} &\cdots & \frac{\partial \varphi}{\partial c_n}\\\\ \frac{\partial \varphi'}{\partial c_1} & \frac{\partial \varphi'}{c_2} & \cdots & \frac{\partial \varphi'}{\partial c_n}\\\\ \vdots & \vdots & \ddots & \vdots\\\\ \frac{\partial \varphi^{(n-1)}}{c_1}&\frac{\partial \varphi^{(n-1)}}{\varphi c_2}&\cdots&\frac{\partial \varphi^{(n-1)}}{\partial c_n}\\ \end{matrix} \right)⎝⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎛​c1​∂φ​∂c1​∂φ′​⋮c1​∂φ(n−1)​​c2​∂φ​c2​∂φ′​⋮φc2​∂φ(n−1)​​⋯⋯⋱⋯​∂cn​∂φ​∂cn​∂φ′​⋮∂cn​∂φ(n−1)​​⎠⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎞​
不为零可以判别微分满秩。从而得知这些常数不相关。

反向理解：如果其中任意两个常数存在关联，那么它们所在的两列线性相关。从而为0
留待思考：如果求出的Jacobi式是一个和主元有关的式子，那么是否会影响通解各常数之间的独立性。

例如：对方程y′′+y=0y''+y= 0y′′+y=0我们可以写出特解y=sin⁡xy=\sin xy=sinx或y=cos⁡xy=\cos xy=cosx，通解y=c1sin⁡x+c2cos⁡xy=c_1\sin x+c_2\cos xy=c1​sinx+c2​cosx
但反智的一点是，存在不能通过特解中常数取特殊值将特值表出的情况。这些解叫做“奇解”。

1.2.5. 初值问题

初始条件与特解的关系
初始条件确定特解：确定了通解中的常数之后，解成为了特解。
研究添加初始条件后被赋予特定意义、解中的常数得到确定的问题，就是初值问题。

例 d2xdt2=−g,dxdy=−gt+c2,x=−12gt2+c2t+c1\frac{\mathrm d^2x}{\mathrm dt^2}=-g,\frac{\mathrm dx}{\mathrm dy}=-gt+c_2,x=-\frac{1}{2}gt^2+c_2t+c_1dt2d2x​=−g,dydx​=−gt+c2​,x=−21​gt2+c2​t+c1​
如果有一个初态x=h,dxdt=0x=h,\frac{\mathrm dx}{\mathrm dt}=0x=h,dtdx​=0，那么h=c1,dxdt=−gt+c2h=c_1,\frac{\mathrm dx}{\mathrm dt}=-gt+c_2h=c1​,dtdx​=−gt+c2​解得：c2=0c_2=0c2​=0

y(x0)=y0,y′(x0)=y0′⋯y(n−1)=y0(n−1)y(x_0)=y_0,y'(x_0)=y'_0\cdots y^{(n-1)}=y_0^{(n-1)}y(x0​)=y0​,y′(x0​)=y0′​⋯y(n−1)=y0(n−1)​

补充：常微分方程的解的图形叫做积分曲线。

1.3. 常微分方程的求解综述

常微分方程的求解是我们这里的核心问题。
我们以一阶线性常系数的齐次、非齐次方程说明基本求解方法，另外以二阶为例，说明高阶方程的一般求解法。

对于非齐次方程求解，要充分利用分治定理，将其分解成齐次方程通解和非齐次特解的和。
求解齐次方程，使用特征根法。
非齐次的特解，使用待定系数法。

2. 一阶微分方程的求解

任何研究都应该从最基本，最直观的情形开始。一开始也不要用特例情况自我束缚、对思维节外生枝。不妨在解决主要问题之后，再来“修补”这些。

重点也只有三大类：即可变量分离的方程、齐次微分方程、一阶线性微分方程。

2.1. 变量分离的方程

2.1.1. 典型变量分离型

形如y′=f(x)g(y)y'=f(x)g(y)y′=f(x)g(y)的方程可以化成F(x)dx=G(x)dxF(x)\,\mathrm dx=G(x)\,\mathrm dxF(x)dx=G(x)dx的形式，其中方程右侧是x，yx，yx，y分别的函数的乘积。
解
dydx=f(x)g(y)dyg(y)=f(x)dx∫1g(y)dy=∫f(x)dx\frac{\mathrm dy}{\mathrm dx}=f(x)g(y)\\ \frac{\mathrm dy}{g(y)}=f(x)\,\mathrm dx\\ \int\frac{1}{g(y)}\,\mathrm dy=\int f(x)\,\mathrm dx dxdy​=f(x)g(y)g(y)dy​=f(x)dx∫g(y)1​dy=∫f(x)dx
这一步先不考虑特殊情况即能不能除，找准主要思路。记
G(y)=∫1g(y)dyF(x)=∫f(x)dxG(y)=F(x)G(y)=\int\frac{1}{g(y)}\,\mathrm dy\\ F(x)=\int f(x)\,\mathrm dx\\ G(y)=F(x) G(y)=∫g(y)1​dyF(x)=∫f(x)dxG(y)=F(x)
这是一个代数方程

现在重新考虑不能相除的情况，即∃bs.t.g(b)=0\exist b\,s.t.g(b)=0∃bs.t.g(b)=0，那么y≡by\equiv by≡b是特解（千万不要忘记啊！）

例 y′=αy(αy'=\alpha y(\alphay′=αy(α是常数)
1ydy=αdxln⁡∣y∣=αx+C∣y∣=ec⋅eαx=C1⋅eαxy=C⋅eαx(c≠0)\frac{1}{y}\,\mathrm dy=\alpha\,\mathrm dx\\ \ln|y|=\alpha x+C\\ |y|=e^c\cdot e^{\alpha x}=C_1\cdot e^{\alpha x}\\ y=C\cdot e^{\alpha x}(c\not=0) y1​dy=αdxln∣y∣=αx+C∣y∣=ec⋅eαx=C1​⋅eαxy=C⋅eαx(c​=0)
加上特解y≡0y\equiv0y≡0，原方程的解可以表示成y=C⋅eαxy=C\cdot e^{\alpha x}y=C⋅eαx

2.1.2. 可化为变量分离型

1。1^。1。形如
dydx=f(ax+by+c)\frac{\mathrm dy}{\mathrm dx}=f(ax+by+c) dxdy​=f(ax+by+c)
的方程：
引入u=ax+by+c,u=ax+by+c,u=ax+by+c,
dudx=a+bdydx=a+bf(ax+by+c)=a+bf(u)\begin{aligned} &\frac{\mathrm du}{\mathrm dx} =a+b\frac{\mathrm dy}{\mathrm dx}\\ =&a+bf(ax+by+c)\\ =&a+bf(u) \end{aligned}==​dxdu​=a+bdxdy​a+bf(ax+by+c)a+bf(u)​
显然uuu是可解的。u=u(x,c~)u=u(x,\tilde{c})u=u(x,c~)
这又是一个代数方程。

特别注意在求解过程当中，要将f(u)=0f(u)=0f(u)=0的情况考虑成一个奇解。

2∘2^\circ2∘
“一阶导数导数=零次齐次函数”的齐次方程

nnn次齐次函数定义为：φ(tytx)=tn(yx)\varphi\left(\frac{ty}{tx}\right)=t^n\left(\frac{y}{x}\right)φ(txty​)=tn(xy​)
y′=f(yx)(x≠0)y'=f\left(\frac{y}{x}\right)(x\not=0) y′=f(xy​)(x​=0)

令u=yx,y=xu,dydx=ddx(xu)=u+xdudx=y′=f(u)f(u)u=\frac{y}{x},y=xu,\frac{\mathrm dy}{\mathrm dx}=\frac{\mathrm d}{\mathrm dx}(xu)=u+x\frac{\mathrm du}{\mathrm dx}\xlongequal{y'=f(u)}f(u)u=xy​,y=xu,dxdy​=dxd​(xu)=u+xdxdu​y′=f(u)f(u)
由于x≠0x\not=0x​=0，dudx=f(u)−ux\frac{\mathrm du}{\mathrm dx}=\frac{f(u)-u}{x}dxdu​=xf(u)−u​

3∘3^\circ3∘另一类齐次
y=f(a1x+b1y+c1ax+by+c)y=f\left(\frac{a_1x+b_1y+c_1}{ax+by+c}\right) y=f(ax+by+ca1​x+b1​y+c1​​)

1. 若c=c1=0,y′=f(a1x+b1yax+by)=f(a1+b1yxa+byx)=g(yx)c=c_1=0, y'=f\left(\frac{a_1x+b_1y}{ax+by}\right)=f\left(\frac{a_1+b_1\frac{y}{x}}{a+b\frac{y}{x}}\right)=g(\frac{y}{x})c=c1​=0,y′=f(ax+bya1​x+b1​y​)=f(a+bxy​a1​+b1​xy​​)=g(xy​)
2. 若不全为0.那么做平移变换。

{x=u+αy=v+βs.t.a1α+b1β+c1=0aα+bβ+c=0\begin{cases} x=u+\alpha\\ y=v+\beta \end{cases}\\ s.t.\, a_1\alpha+b_1\beta+c_1=0\\ a\alpha+b\beta+c=0 {x=u+αy=v+β​s.t.a1​α+b1​β+c1​=0aα+bβ+c=0
则
dvdu=f(a1u+b1vau+bv)\frac{\mathrm dv}{\mathrm du}=f\left(\frac{a_1u+b_1v}{au+bv}\right) dudv​=f(au+bva1​u+b1​v​)
利用Cramer法则，
∣a1b1ab∣≠0\begin{vmatrix} a_1&b_1\\ a&b \end{vmatrix}\not=0∣∣∣∣​a1​a​b1​b​∣∣∣∣​​=0
时，α,β\alpha, \betaα,β存在且唯一，可解。
∣a1b1ab∣=0\begin{vmatrix} a_1&b_1\\ a&b \end{vmatrix}=0∣∣∣∣​a1​a​b1​b​∣∣∣∣​=0
时，设a1a=b1b=k\frac{a_1}{a}=\frac{b_1}{b}=kaa1​​=bb1​​=k
y′=f(a1x+b1y+c1ax+by+c)=f(k(ax+by)+c1ax+by+c)y'=f\left(\frac{a_1x+b_1y+c_1}{ax+by+c}\right)=f\left(\frac{k(ax+by)+c_1}{ax+by+c}\right) y′=f(ax+by+ca1​x+b1​y+c1​​)=f(ax+by+ck(ax+by)+c1​​)
现令
dudx=a+bdydx=a+bf(ku+c1u+c)\frac{\mathrm du}{\mathrm dx}=a+b\frac{\mathrm dy}{\mathrm dx}=a+bf(\frac{ku+c_1}{u+c}) dxdu​=a+bdxdy​=a+bf(u+cku+c1​​)
变量分离，又是可解的。
dydx=x+y+2x−y+1\frac{\mathrm d y}{\mathrm dx}=\frac{x+y+2}{x-y+1} dxdy​=x−y+1x+y+2​

2.2. 一阶线性方程

y′+p(x)y=q(x),p,q∈C0(a,b)y'+p(x)y=q(x), p,q\in C^0(a,b) y′+p(x)y=q(x),p,q∈C0(a,b)
(1)
其中的C0C^0C0表示连续，C1C^1C1表示连续可微。

2.2.1. 一阶线性方程的判别

关于未知函数yyy与其导数dydx\frac{\mathrm dy}{\mathrm dx}dxdy​都是一次的。如果q(x)=0q(x)=0q(x)=0成为一阶齐次线性微分方程，否则成为一阶非齐次。

2.2.2. 一阶线性方程的计算

2.2.2.1. 1∘1^\circ1∘ 最简单的情况：q(x)=0q(x)=0q(x)=0（齐次）

这个“齐次"和我们前述的齐次函数的齐次并不相同。
y′+p(x)y=0y'+p(x)y=0 y′+p(x)y=0
(2)
dydx=−p(x)y\frac{\mathrm dy}{\mathrm dx}=-p(x)y dxdy​=−p(x)y
首先我们发现y=0y=0y=0是一个奇解。然后移项积分：
ln∣y∣=−∫x0xP(t)dt+Cln|y|=-\int_{x_0}^xP(t)\,\mathrm dt+C ln∣y∣=−∫x0​x​P(t)dt+C
y=±ece(−∫x0xP(t)dt)y=\pm e^c e^{\left(-\int_{x_0}^xP(t)\,\mathrm dt\right)} y=±ece(−∫x0​x​P(t)dt)
结合奇解
y=C⋅e(−∫x0xP(t)dt),∀Cy=C\cdot e^{\left(-\int_{x_0}^xP(t)\,\mathrm dt\right)},\forall C y=C⋅e(−∫x0​x​P(t)dt),∀C
这个结果告诉我们，线性齐次方程的解仅有两种情况：≡0\equiv0≡0或≢0\not\equiv0​≡0，
讨论：

• φ1(x),φ2(x)\varphi_1(x), \varphi_2(x)φ1​(x),φ2​(x)是（2）的两个解，那么∀k1,k2∈R,k1φ1+k2φ2\forall k_1,k_2\in\R, k_1\varphi_1+k_2\varphi_2∀k1​,k2​∈R,k1​φ1​+k2​φ2​是(2)(2)(2)的解。
• 若φ(x)\varphi(x)φ(x)是（2）的解，则φ(x)∈C(a,b)\varphi(x)\in C(a,b)φ(x)∈C(a,b)
• 一般非线性方程可能存在解的长度比定义域长度短的情况。例：y′=1+y2,y=tan⁡x,x∈(−π2,π2)y'=1+y^2, y=\tan x, x\in(-\frac{\pi}{2}, \frac{\pi}{2})y′=1+y2,y=tanx,x∈(−2π​,2π​)，对于线性，一定是一样长的。

2.2.2.2. 2∘2^\circ2∘非齐次

2.2.2.2.1. 解的性质：

• 若y=φ(x),y=ψ(x)y=\varphi(x),y=\psi(x)y=φ(x),y=ψ(x)是（1）的解，则y=φ(x)−ψ(x)y=\varphi(x)-\psi(x)y=φ(x)−ψ(x)是（2）的解
• 若y=φ(x)y=\varphi(x)y=φ(x)是（1）的解，y=Φy=\Phiy=Φ是（2）的解，那么y=φ(x)+Φ(x)y=\varphi(x)+\Phi(x)y=φ(x)+Φ(x)是（1）的解。

2.2.2.2.2. 求解思路：

• 从而推知若y=φ(x)y=\varphi(x)y=φ(x)是（1）的一个解。那么（1）的所有解可以写成y=φ(x)+Ce−∫x0xP(t)dty=\varphi(x)+Ce^{-\int_{x_0}^xP(t)\,\mathrm dt}y=φ(x)+Ce−∫x0​x​P(t)dt
• 由上述，这个方程的求解转化成求解(1)的一个特解。

2.2.2.2.3. 方法：常数变易法

将齐次方程的解中的常数变成一个函数表达式，待定为u(x)u(x)u(x)
y′=u′(x)⋅e−∫x0xP(t)dt−u(x)P(x)e−∫x0xP(t)dt=Q−P(x)yy'=u'(x)\cdot e^{-\int_{x_0}^xP(t)\,\mathrm dt}-u(x)P(x)e^{-\int_{x_0}^xP(t)\,\mathrm dt}=Q-P(x)y y′=u′(x)⋅e−∫x0​x​P(t)dt−u(x)P(x)e−∫x0​x​P(t)dt=Q−P(x)y
化简得：
u′(x)⋅e−∫x0xP(t)dt=Qu'(x)\cdot e^{-\int_{x_0}^xP(t)\,\mathrm dt}=Q u′(x)⋅e−∫x0​x​P(t)dt=Q
u′(x)=Q(x)⋅e∫x0xP(t)dtu'(x)=Q(x)\cdot e^{\int_{x_0}^xP(t)\,\mathrm dt} u′(x)=Q(x)⋅e∫x0​x​P(t)dt
这个式子可不可积其实还有待讨论……但总之把它积起来就可以得到我们需要的u(x)u(x)u(x)
y=[∫x0xQ(t)e∫x0tP(s)dsdt+C]e−∫x0xP(t)dty=\left[\int_{x_0}^xQ(t)e^{\int_{x_0}^tP(s)\,\mathrm ds}\,\mathrm dt+C\right]e^{-\int_{x_0}^xP(t)\,\mathrm dt} y=[∫x0​x​Q(t)e∫x0​t​P(s)dsdt+C]e−∫x0​x​P(t)dt

2.2.2.2.4. 步骤整理：

• step1：求齐次方程通解C⋅y1(x)C\cdot y_1(x)C⋅y1​(x)
• step2：u′=Q(x)y1(x)u'=\frac{Q(x)}{y_1(x)}u′=y1​(x)Q(x)​
• step3：移项积分，整理得解

若f∈C(a,+∞)f\in C(a, +\infty)f∈C(a,+∞)且f(x)+f′(x)→0f(x)+f'(x)\to0f(x)+f′(x)→0则f(x)→0f(x)\to0f(x)→0

另例：
f(x)∈C1[0,+∞),a(x)≥C0>0,a(x)∈C[0,+∞)f(x)\in C^1[0,+\infty), a(x)\geq C_0>0, a(x)\in C[0,+\infty)f(x)∈C1[0,+∞),a(x)≥C0​>0,a(x)∈C[0,+∞)且lim⁡x→+∞(f′(x)+a(x)f(x))=0\lim\limits_{x\to+\infty}(f'(x)+a(x)f(x))=0x→+∞lim​(f′(x)+a(x)f(x))=0则lim⁡x→+∞f(x)=0\lim\limits_{x\to+\infty}f(x)=0x→+∞lim​f(x)=0

令b(x)=f′(x)+a(x)f(x),b(x)∈C0[0,+∞),lim⁡x→+∞b(x)=0b(x)=f'(x)+a(x)f(x), b(x)\in C^0[0,+\infty),\lim\limits_{x\to+\infty} b(x)=0b(x)=f′(x)+a(x)f(x),b(x)∈C0[0,+∞),x→+∞lim​b(x)=0，从bbb中解出f(x)f(x)f(x)
f(x)=Ce−∫x0xa(t)dt+e−∫x0xa(t)dt⋅∫0xb(s)e∫0sa(t)dtdsf(x)=Ce^{-\int_{x_0}^xa(t)\,\mathrm dt}+e^{-\int_{x_0}^xa(t)\,\mathrm dt}\cdot \int_0^xb(s)e^{\int_0^sa(t)\,\mathrm dt}\,\mathrm ds f(x)=Ce−∫x0​x​a(t)dt+e−∫x0​x​a(t)dt⋅∫0x​b(s)e∫0s​a(t)dtds
其中取初值条件，可以得C=f(0)C=f(0)C=f(0)
f(x)=f(0)+∫0xb(s)e∫0sa(t)dtdse∫0xa(t)dtf(x)=\frac{f(0)+\int_0^xb(s)e^{\int_0^sa(t)\,\mathrm dt}\,\mathrm ds}{e^{\int_0^xa(t)\,\mathrm dt}} f(x)=e∫0x​a(t)dtf(0)+∫0x​b(s)e∫0s​a(t)dtds​
分母趋向于无穷（且求导后无穷，即满足洛必达法则），
即有lim⁡x→+∞f(x)=lim⁡x→+∞b(x)a(x)=0\lim\limits_{x\to+\infty}f(x)=\lim\limits_{x\to+\infty}\frac{b(x)}{a(x)}=0x→+∞lim​f(x)=x→+∞lim​a(x)b(x)​=0

2.2.2.3. 3∘3^\circ3∘Bernoulli方程

这是一个将非线性方程化成线性方程的好例子。
dydx+P(x)y=Q(x)yα\frac{\mathrm dy}{\mathrm dx}+P(x)y=Q(x)y^\alpha dxdy​+P(x)y=Q(x)yα
y−αdydx+P⋅y1−α=Q(x)y^{-\alpha}\frac{\mathrm dy}{\mathrm dx}+P\cdot y^{1-\alpha}=Q(x) y−αdxdy​+P⋅y1−α=Q(x)
令z=1−αz=1-\alphaz=1−α,
11−αdzdx+P(x)z=Q(x)\frac{1}{1-\alpha}\frac{\mathrm dz}{\mathrm dx}+P(x)z=Q(x) 1−α1​dxdz​+P(x)z=Q(x)
这是一个线性方程。

2.3. 全微分方程

2.3.1. 全微分方程的形式

dydx=−PQ\frac{\mathrm dy}{\mathrm dx}=-\frac{P}{Q} dxdy​=−QP​
对称形式
P(x,y)dx+Q(x,y)dy=0P(x,y)\,\mathrm dx+Q(x,y)\,\mathrm dy=0 P(x,y)dx+Q(x,y)dy=0

P,Q∈C1(D)P,Q\in C^1(D)P,Q∈C1(D)(这里要可微，因为如果不可导的话，我们的方法实在有限)

若∃F(x,y)∈C1(D),dF(x,y)=Pdx+Qdy⇒dF(x,y)=0\exist F(x,y)\in C^1(D), \,\mathrm dF(x,y)=P\,\mathrm dx+Q\,\mathrm dy\Rightarrow \mathrm dF(x,y)=0∃F(x,y)∈C1(D),dF(x,y)=Pdx+Qdy⇒dF(x,y)=0

从而我们求这样的方程的解的问题，转化成：

1. 判断FFF的存在性
2. 求出FFF

在二型线和路径无关的条件时，我们写出了几个这个FFF存在的等价条件。
在这里我们常用
∂P∂y=∂Q∂x\frac{\partial P}{\partial y}=\frac{\partial Q}{\partial x} ∂y∂P​=∂x∂Q​
必要性利用连续函数的混合导数相等说明。

2.3.2. 判别条件的充分性证明：

偏积分得
F=∫x0xP(s,y)ds+φ(y)F=\int_{x_0}^xP(s,y)\,\mathrm ds+\varphi(y) F=∫x0​x​P(s,y)ds+φ(y)
(1)
Q=∂∂y∫x0xP(s,y)ds+φ′(y)Q=\frac{\partial}{\partial y}\int_{x_0}^xP(s,y)\,\mathrm ds+\varphi'(y) Q=∂y∂​∫x0​x​P(s,y)ds+φ′(y)
由于光滑，求导积分可交换。
Q=∫x0x∂∂yP(s,y)ds+φ′(y)Q=\int_{x_0}^x\frac{\partial}{\partial y}P(s,y)\,\mathrm ds+\varphi'(y) Q=∫x0​x​∂y∂​P(s,y)ds+φ′(y)
再由充分条件得：
∫x0x∂∂xQ(s,y)ds+φ′(y)=Q(x,y)\int_{x_0}^x\frac{\partial}{\partial x}Q(s,y)\,\mathrm ds+\varphi'(y)=Q(x,y) ∫x0​x​∂x∂​Q(s,y)ds+φ′(y)=Q(x,y)
Q(x,y)−Q(x0,y)+φ′(y)=Q(x,y)Q(x,y)-Q(x_0,y)+\varphi'(y)=Q(x,y) Q(x,y)−Q(x0​,y)+φ′(y)=Q(x,y)
得
φ′(y)=Q(x0,y)\varphi'(y)=Q(x_0,y) φ′(y)=Q(x0​,y)
φ=∫y0yQ(x0,t)dt\varphi=\int_{y_0}^yQ(x_0,t)\,\mathrm dt φ=∫y0​y​Q(x0​,t)dt
结合（1）得:
F(x,y)=∫x0xP(s,y)ds+∫y0yQ(x0,t)dtF(x,y) =\int_{x_0}^xP(s,y)\,\mathrm ds+\int_{y_0}^yQ(x_0,t)\,\mathrm dt F(x,y)=∫x0​x​P(s,y)ds+∫y0​y​Q(x0​,t)dt
说人话，就是说明了二型线与路径无关的条件是全微分方程直接用原函数是解的条件。

例 (xcos⁡y+2xy2)dx+(−12x2sin⁡y+2x2y+y2)=0(x\cos y+2xy^2)\,\mathrm dx+(-\frac{1}{2}x^2\sin y+2x^2y+y^2)=0(xcosy+2xy2)dx+(−21​x2siny+2x2y+y2)=0
∂P∂y=∂Q∂x\frac{\partial P}{\partial y}=\frac{\partial Q}{\partial x}∂y∂P​=∂x∂Q​然后使用偏积分求出通积分。

接下来要解决的问题是，如果∂P∂y\frac{\partial P}{\partial y}∂y∂P​和∂Q∂x\frac{\partial Q}{\partial x}∂x∂Q​不相等呢？

2.3.3. 积分因子

对以上的问题，基本想法是：
寻找一个（就够的）μ=μ(x,y)≠0\mu=\mu(x,y)\not=0μ=μ(x,y)​=0s.t.
μPdx+μQdy=0\mu P\,\mathrm dx+\mu Q\,\mathrm dy=0 μPdx+μQdy=0
是全微分方程。

μ\muμ满足μP=∂F∂x,μQ=∂F∂y\mu P=\frac{\partial F}{\partial x}, \mu Q=\frac{\partial F}{\partial y}μP=∂x∂F​,μQ=∂y∂F​

结合光滑：
∂∂y(μP)=∂∂x(μQ)\frac{\partial }{\partial y}(\mu P)=\frac{\partial }{\partial x}(\mu Q) ∂y∂​(μP)=∂x∂​(μQ)
转化成：

P∂μ∂y−Q∂μ∂x=μ(∂Q∂x−∂P∂y)P\frac{\partial \mu}{\partial y}-Q\frac{\partial \mu}{\partial x}=\mu\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right) P∂y∂μ​−Q∂x∂μ​=μ(∂x∂Q​−∂y∂P​)

这是一个一阶偏微分方程。
但我们仍只需要找到一个就够。

若存在μ=μ(x)\mu=\mu(x)μ=μ(x)
μ′μ=−∂Q∂x−∂P∂yQ(x,y)\frac{\mu'}{\mu}=-\frac{\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}}{Q(x,y)} μμ′​=−Q(x,y)∂x∂Q​−∂y∂P​​
若右端只依赖xxx那么就是变量分离的方程，可求这个积分因子。
μ=e∫x0xf(s)ds\mu=e^{\int_{x_0}^xf(s)\,\mathrm ds} μ=e∫x0​x​f(s)ds
实际求解过程中，我们只需要看P,QP,QP,Q二者谁和右端相除可以化成单变量式。

附：另一个常用的表达式；
μ′μ=∂Q∂x−∂P∂xP(x,y)\frac{\mu'}{\mu}=\frac{\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial x}}{P(x,y)} μμ′​=P(x,y)∂x∂Q​−∂x∂P​​

3. 微分方程解的存在、唯一性定理

这是节外生枝的一部分，但也具有相当高的理论意义。

如果方程
{y′=f(x,y)y(x0)=y0(1)\begin{cases} y'=f(x,y)\tag1\\ y(x_0)=y_0 \end{cases}{y′=f(x,y)y(x0​)=y0​​(1)
若fff连续，且满足Lipschitz条件（斜率不大于某个值）：
∃L>0,∀(x,y1),(x,y0)∈D\exist L>0,\forall(x,y_1),(x,y_0)\in D∃L>0,∀(x,y1​),(x,y0​)∈D
∣f(x,y0)−f(x,y1)∣≤L∣y1−y0∣|f(x,y_0)-f(x,y_1)|\leq L|y_1-y_0| ∣f(x,y0​)−f(x,y1​)∣≤L∣y1​−y0​∣
那么其在∣x−x0∣≤h|x-x_0|\leq h∣x−x0​∣≤h解存在且唯一
h=min⁡{a,bM}M=max⁡D∣f(x,y)∣\begin{aligned} h=&\min\{a,\frac{b}{M}\}\\ M=&\max\limits_D|f(x,y)| \end{aligned}h=M=​min{a,Mb​}Dmax​∣f(x,y)∣​

• 注1：hhh是越界保护。其中MMM通过(1)中代换得到的：y′≤My'\leq My′≤M。

证明：
Step1：
原式等价于
y(x)−y0=∫x0xf(s,y(s))dsy(x)-y_0=\int_{x_0}^xf(s,y(s))\,\mathrm ds y(x)−y0​=∫x0​x​f(s,y(s))ds
Step2：
构造Picard序列：
φn+1(x)=y0+∫x0xf(s,φn(s))ds\varphi_{n+1}(x)=y_0+\int_{x_0}^xf(s,\varphi_n(s))\,\mathrm ds φn+1​(x)=y0​+∫x0​x​f(s,φn​(s))ds
这也是近似求解的方法。逐次逼近为我们提供了一种求近似解的计算方法。

Step3：
序列一致收敛，故存在（证明略

Step4：
解的唯一性判定
∣φ(x)−ψ(x)∣≤∣∫x0xf(s,φ(s))−f(s,ψ(s))ds∣≤L∫x0x∣φ(s)−ψ(s)∣ds|\varphi(x)-\psi(x)|\leq \left|\int_{x_0}^xf(s,\varphi(s))-f(s,\psi(s))\,\mathrm ds\right| \leq L\int_{x_0}^x|\varphi(s)-\psi(s)|\,\mathrm ds ∣φ(x)−ψ(x)∣≤∣∣∣∣​∫x0​x​f(s,φ(s))−f(s,ψ(s))ds∣∣∣∣​≤L∫x0​x​∣φ(s)−ψ(s)∣ds
∣φ(x)−ψ(x)∣≤LM0∣x−x0∣|\varphi(x)-\psi(x)|\leq LM_0|x-x_0| ∣φ(x)−ψ(x)∣≤LM0​∣x−x0​∣
故（2）代入（1）得
∣φ(x)−ψ(x)∣≤L2M02∣x−x0∣2|\varphi(x)-\psi(x)|\leq \frac{L^2M_0}{2}|x-x_0|^2 ∣φ(x)−ψ(x)∣≤2L2M0​​∣x−x0​∣2
累次代入，得：
∣φ(x)−ψ(x)∣≤LnM0n!hn|\varphi(x)-\psi(x)|\leq\frac{L^nM_0}{n!}h^n ∣φ(x)−ψ(x)∣≤n!LnM0​​hn
显然，
φ(x)≡ψ(x)\varphi(x)\equiv\psi(x) φ(x)≡ψ(x)

注1（Cauchy条件）
f∈C1(D)f\in C^1(D)f∈C1(D)
f(x,y1)−f(x,y2)=∂f∂y(x,ξ)(y1−y0)f(x,y_1)-f(x,y_2)=\frac{\partial f}{\partial y}(x,\xi)(y_1-y_0) f(x,y1​)−f(x,y2​)=∂y∂f​(x,ξ)(y1​−y0​)
显然这是满足Lipschitz条件的。

注2 满足Lipschitz条件的未必处处可微。例如f(x,y)=∣y∣f(x,y)=|y|f(x,y)=∣y∣

推论：对y′=f(x,y),(x,y)∈Dy'=f(x,y),(x,y)\in Dy′=f(x,y),(x,y)∈D，若函数f(x,y)f(x,y)f(x,y)和fyf_yfy​在区域上连续，那么过D内任一点(x0,y0)(x_0,y_0)(x0​,y0​)，有且仅有一条方程的积分曲线通过。（利用Cantor定理可说一致连续）

注3 虽然Picard定理所求的范围有限。但是我们可以通过不断的延拓，仍能获得很好的结果

4. 高阶举例——二阶方程的求解

y′′(x)+p(x)y′(x)+q(x)y(x)=f(x)(3)y''(x)+p(x)y'(x)+q(x)y(x)=f(x)\tag3 y′′(x)+p(x)y′(x)+q(x)y(x)=f(x)(3)
y′′(x)+p(x)y′+q(x)y=0(4)y''(x)+p(x)y'+q(x)y=0\tag4 y′′(x)+p(x)y′+q(x)y=0(4)
性质1：显然（4）的解具有线性性质。

对于二阶方程，我们首先有这样一类容易解决的问题。

4.1. 可降阶的方程

这是解决高阶方程的一个很简便、清晰的思路。但条件有限制。

常用 不显含xxx的二阶方程（通常设题比较刁钻，只能用此法解）
F(y,y′,y′′)=0F(y,y',y'')=0 F(y,y′,y′′)=0
令
P=y′=dydxP=y'=\frac{\mathrm dy}{\mathrm dx} P=y′=dxdy​
y′′=ddxy′=ddyy′dydx=dPdyPy''=\frac{\mathrm d}{\mathrm dx}y'=\frac{\mathrm d}{\mathrm dy}y'\frac{\mathrm dy}{\mathrm dx}=\frac{\mathrm dP}{\mathrm dy}P y′′=dxd​y′=dyd​y′dxdy​=dydP​P
原式可以转化为
F1(y,P,dPdy)=0F_1(y,P,\frac{\mathrm dP}{\mathrm dy})=0 F1​(y,P,dydP​)=0
这是一阶线性方程。
所得解P=dydxP=\frac{\mathrm dy}{\mathrm dx}P=dxdy​又是变量分离的。

另外还有F(x,y′,y′′)=0F(x,y',y'')=0F(x,y′,y′′)=0型（将y′y'y′记为PPP）。

这是一个快捷化简。

4.2. 函数的线性相关

4.2.1. 判据——W行列式

其他的方程在求解上具有一定的复杂性，不能通过初等方法求解，但我们可以通过一定的方式对其解的结构进行研究。

设g1(x),g2(x)∈D(a,b),∃c1,c2,s.t.┐(c1=0∧c2=0)∧c1g1+c2g2≡0g_1(x),g_2(x)\in D(a,b),\exist c_1,c_2,s.t.\,\urcorner(c_1=0\wedge c_2=0)\wedge c_1g_1+c_2g_2\equiv0g1​(x),g2​(x)∈D(a,b),∃c1​,c2​,s.t.┐(c1​=0∧c2​=0)∧c1​g1​+c2​g2​≡0则称g1,g2g_1,g_2g1​,g2​线性相关。
结合性质1，我们提出以下引理：
设φ1(x),φ2(x)\varphi_1(x), \varphi_2(x)φ1​(x),φ2​(x)是（4）的两个解。那么它们在(a,b)(a,b)(a,b)上线性相关的充要条件是：它们确定的Wronski行列式
W=∣φ1φ2φ1′φ2′∣≡0W=\begin{vmatrix} \varphi_1&\varphi_2\\ \varphi_1'&\varphi_2' \end{vmatrix}\equiv0 W=∣∣∣∣​φ1​φ1′​​φ2​φ2′​​∣∣∣∣​≡0
证明：必要性对（4）和（4）的一阶导数式利用Cramer法则。充分性利用解的存在唯一性。详见高数下P180

4.2.2. Wronski行列式的计算

w(x)=φ1(x)φ2′(x)−φ2(x)φ1′(x)w′(x)=φ1′φ2′+φ1φ2′′−φ2′φ1′−φ2φ1′′=φ1φ2′′−φ2φ1′′w(x)=\varphi_1(x)\varphi_2'(x)-\varphi_2(x)\varphi'_1(x)\\ w'(x)=\varphi'_1\varphi'_2+\varphi_1\varphi_2''-\varphi'_2\varphi_1'-\varphi_2\varphi''_1=\varphi_1\varphi''_2-\varphi_2\varphi_1'' w(x)=φ1​(x)φ2′​(x)−φ2​(x)φ1′​(x)w′(x)=φ1′​φ2′​+φ1​φ2′′​−φ2′​φ1′​−φ2​φ1′′​=φ1​φ2′′​−φ2​φ1′′​
方程组的秩为3，三个变量，可解。
φ1′′+pφ1′+qφ1=0φ2′′+pφ2′+qφ2=0\varphi''_1+p\varphi'_1+q\varphi_1=0\\ \varphi_2''+p\varphi'_2+q\varphi_2=0 φ1′′​+pφ1′​+qφ1​=0φ2′′​+pφ2′​+qφ2​=0
代入消元可得：
w′=p⋅(−w)w'=p\cdot(-w) w′=p⋅(−w)
解这个变量分离的方程有：
w(x)=w(x0)e−∫x0xp(s)ds(Liouville公式)w(x)=w(x_0)e^{-\int_{x_0}^xp(s)\,\mathrm ds} (Liouville公式) w(x)=w(x0​)e−∫x0​x​p(s)ds(Liouville公式)
这个公式告诉我们，在待求区间(a,b)(a,b)(a,b)上如果有一点的Wronski行列式为0，那么整个区间上都为零。

4.2.3. 二阶线性齐次微分方程

前面已经讨论了函数线性无关的条件。

对于方程（4），我们只需要找到两个线性无关的解，这两个解的初值分别为(0,1)(0,1)(0,1)和(1,0)(1,0)(1,0)，
那么所有解都可以经由这两个解表出。

为什么是两个？二阶方程看似是三个不定元，但实际上方程会隐式地确定xxx，所以只要将带有积分常数CCC的y,y′y,y'y,y′的初值确定下来，那y′′y''y′′也就是唯一确定的了。所以这个方程如果展成方程组，对于y(n)y^{(n)}y(n)结构来看，显示的秩只有n−1n-1n−1。因而我们只需求解两个初值即可。

4.3. 求解二阶线性非齐次方程

y′′+p(x)y′+q(x)y=f(x)(3)y''+p(x)y'+q(x)y=f(x)\tag3 y′′+p(x)y′+q(x)y=f(x)(3)
设y=φ1(x),y=φ2(x)y=\varphi_1(x), y=\varphi_2(x)y=φ1​(x),y=φ2​(x)是（4）的线性无关解。
y=φ0(x)y=\varphi_0(x)y=φ0​(x)是（3）的一个特解。

则（3）的所有解都可以表示成：
y=c1ψ1(x)+c2ψ(x)+φ0(x)y=c_1\psi_1(x)+c_2\psi(x)+\varphi_0(x) y=c1​ψ1​(x)+c2​ψ(x)+φ0​(x)
这仍然是由方程的线性性质所保证的。
求解过程就退化成找特解。

为了扩大可解的范围，我们给出如下的

4.3.1. 分治定理

设y=φ1(x)y=\varphi_1(x)y=φ1​(x)是
y′′+py+qy=f1y''+py+qy=f_1 y′′+py+qy=f1​
的解，
设y=f2(x)y=f_2(x)y=f2​(x)是
y′′+py′+qy=f2y''+py'+qy=f_2 y′′+py′+qy=f2​
的解。

那么φ1+φ2\varphi_1+\varphi_2φ1​+φ2​是原方程的解。
这为分治的思路提供了指导。

在正式研究求解之前，我们给出关于可解性的一般规律：如果p,qp,qp,q与xxx有关，一般不可解。
例如：
y′′+a(x)y=0a(x)=1+εcos⁡xy''+a(x)y=0\\ a(x)=1+\varepsilon\cos x y′′+a(x)y=0a(x)=1+εcosx
除了y≡0y\equiv0y≡0以外我们找不到其他的有用的解。所以接下来，我们着重讨论二阶方程一种比较简单的形式：

4.3.2. 求通解——二阶常系数微分方程

y′′+py′+qy=f(x)(5)y''+py'+qy=f(x)\tag5 y′′+py′+qy=f(x)(5)
核心思想是尝试，方法是待定系数法。找到线性无关的解基，我们就赢了。

首先解决齐次情形：
这里考虑到eee指数函数的求导性质优良，我们优先考虑形如y=eλxy=e^{\lambda x}y=eλx的解。代入得：
(λ2+pλ+q)eλx=0(\lambda^2+p\lambda+q)e^{\lambda x}=0 (λ2+pλ+q)eλx=0
eλxe^{\lambda x}eλx是解⇔λ2+pλ+q=0\Leftrightarrow \lambda^2+p\lambda+q=0⇔λ2+pλ+q=0

case 1：特征根相异时y1=eλ1x,y2=eλ2xy_1=e^{\lambda_1 x},y_2=e^{\lambda_2x}y1​=eλ1​x,y2​=eλ2​x是两个线性无关的解
W(x)=∣eλ1xeλ2xλ1eλ1xλ2eλ2x∣=(λ2−λ1)e(λ1+λ2)xW(x)=\begin{vmatrix} e^{\lambda_1x}&e^{\lambda_2x}\\ \lambda_1e^{\lambda_1x}&\lambda_2e^{\lambda_2x} \end{vmatrix}=(\lambda_2-\lambda_1)e^{(\lambda_1+\lambda_2)x} W(x)=∣∣∣∣​eλ1​xλ1​eλ1​x​eλ2​xλ2​eλ2​x​∣∣∣∣​=(λ2​−λ1​)e(λ1​+λ2​)x

case2：特征根为二重根时。第一个根是容易找到的。常用方法是对另一根进行乘xnx^nxn.
这里乘xxx就够用了y2=xeλxy_2=xe^{\lambda x}y2​=xeλx经验证，也是一个解。
且y1=eλ1x,y2=xeλ2xy_1=e^{\lambda_1x},y_2=xe^{\lambda_2x}y1​=eλ1​x,y2​=xeλ2​x可得Wronski行列式不为零。

有了以上的讨论，我们就可以解特征方程Δ>0\Delta >0Δ>0的问题了。

case3：特征方程有共轭复根
λ1,2=α±iβ,β>0\lambda_{1,2}=\alpha\pm i\beta, \beta>0λ1,2​=α±iβ,β>0

y1∗(x)=e(α+iβ)x=e(αx(cos⁡βx+isin⁡βx)y2∗(x)=e(α−iβ)x=eαx(cos⁡βx−isin⁡βx)y^*_1(x)=e^{(\alpha+i\beta )x}=e^{(\alpha x}(\cos\beta x+i\sin\beta x)\\ y^*_2(x)=e^{(\alpha-i\beta) x}=e^{\alpha x}(\cos\beta x-i\sin\beta x) y1∗​(x)=e(α+iβ)x=e(αx(cosβx+isinβx)y2∗​(x)=e(α−iβ)x=eαx(cosβx−isinβx)
消虚数，取得两个实函数解基：
y1(x)=12(y1∗+y2∗)=eαxcos⁡βxy2(x)=12i(y1∗−y2∗)=eαxsin⁡βxy_1(x)=\frac{1}{2}(y_1^*+y_2^*)=e^{\alpha x}\cos\beta x\\ y_2(x)=\frac{1}{2i}(y^*_1-y_2^*)=e^{\alpha x}\sin \beta x y1​(x)=21​(y1∗​+y2∗​)=eαxcosβxy2​(x)=2i1​(y1∗​−y2∗​)=eαxsinβx
由于β>0\beta >0β>0故而Wronski行列式为βe2αx≠0\beta e^{2\alpha x}\not=0βe2αx​=0

于是我们可以将通解形式总结为：

特征方程判别式	特征根	通解形式
Δ>0\Delta >0Δ>0	两个相异实根λ1,λ2\lambda_1, \lambda_2λ1​,λ2​	C1eλ1x+C2eλ2xC_1e^{\lambda_1 x}+C_2e^{\lambda_2x}C1​eλ1​x+C2​eλ2​x
Δ=0\Delta = 0Δ=0	二重根λ1\lambda_1λ1​	(C1+C2x)⋅eλ1x(C_1+C_2x)\cdot e^{\lambda_1}x(C1​+C2​x)⋅eλ1​x
Δ<0\Delta <0Δ<0	λ1,2=α±iβ\lambda_{1,2}=\alpha\pm i\betaλ1,2​=α±iβ	eαx(C1cos⁡βx+C2sin⁡βx)e^{\alpha x}(C_1\cos\beta x+C_2\sin\beta x)eαx(C1​cosβx+C2​sinβx)

• 从某个角度上说，非重根的所有情况都是第三种情形的退化情况。
• 注意这个βx\beta xβx是放在一起写的。这是eixe^{i x}eix的De Moivre展开。

其实在常系数的情况下，高阶的也相对容易求解。

特征根	通解
单重实根λ\lambdaλ	eλxe^{\lambda x}eλx
k重实根(k>1k>1k>1)λ\lambdaλ	eλx,xeλx,⋯,xk−1eλxe^{\lambda x},xe^{\lambda x}, \cdots, x^{k-1}e^{\lambda x}eλx,xeλx,⋯,xk−1eλx
单共轭复根λ1,2=α+iβ\lambda_{1,2}=\alpha+i\betaλ1,2​=α+iβ	eαxcos⁡βx,eαxsin⁡βxe^{\alpha x}\cos\beta x,e^{\alpha x}\sin\beta xeαxcosβx,eαxsinβx
nnn重共轭复根λ1,2=α+iβ(n>1)\lambda_{1,2}=\alpha+i\beta(n>1)λ1,2​=α+iβ(n>1)	eαxcos⁡βx,eαxsin⁡βx,⋯,xn−1eαxcos⁡βx,xm−1sin⁡βxe^{\alpha x}\cos\beta x,e^{\alpha x}\sin\beta x, \cdots,x^{n-1}e^{\alpha x}\cos\beta x,x^{m-1}\sin\beta xeαxcosβx,eαxsinβx,⋯,xn−1eαxcosβx,xm−1sinβx

4.3.3. 求特解——二阶常系数线性非齐次方程

虽然已经研究了关于线性齐次方程的问题。但是，解决大部分非齐次仍然是困难的。
在这里，我们着重讨论以下几种：

1. f(x)=Pn(x)f(x)=P_n(x)f(x)=Pn​(x)
2. f(x)=a⋅eαxf(x)=a\cdot e^{\alpha x}f(x)=a⋅eαx
3. f(x)=acos⁡βx+bsin⁡βxf(x)=a\cos\beta x+b\sin\beta xf(x)=acosβx+bsinβx

它们的方法都相对固定：核心方法是待定系数，目标是找通解，退化时乘xxx

我们从最简单的情形开始：

4.3.3.1. 多项式型

利用待定系数法，设解为：
Qn(x)=b0xn+b1xn−1+⋯+bn,b0≠0Q_n(x)=b_0x^n+b_1x^{n-1}+\cdots +b_n, b_0\not=0 Qn​(x)=b0​xn+b1​xn−1+⋯+bn​,b0​​=0
（0不是特征根）case 1:
q≠0q\not=0q​=0时，左右都是nnn次多项式。比较各次项系数，解出bib_ibi​

（0是单特征根）case 2:
考虑Q(x)=xQn(x)Q(x)=xQ_n(x)Q(x)=xQn​(x)，代入得：
Q′(x)=Qn(x)+xQn′(x)Q′′(x)=2Qn′(x)+xQn′′(x)Q'(x)=Q_n(x)+xQ_n'(x)\\ Q''(x)=2Q'_n(x)+xQ_n''(x) Q′(x)=Qn​(x)+xQn′​(x)Q′′(x)=2Qn′​(x)+xQn′′​(x)
整理出nnn次项作为首项的式子：
p(xQn′+Qn)+(xQn′′+2Qn′)=Pnp(xQ_n'+Q_n)+(xQ_n''+2Q_n')=P_n p(xQn′​+Qn​)+(xQn′′​+2Qn′​)=Pn​
如果p≠0p\not=0p​=0退化成case 1中的讨论情形。

（0是重特征根）case 3:
对case 3中的p=0p=0p=0情形，考虑Qn(x)=x2QnQ_n(x)=x^2Q_nQn​(x)=x2Qn​可以直接代入对比系数（北大P188）。

史逸老师提供了一种方便的理解，case3包含条件p=0,q=0p=0,\ q=0p=0, q=0，原式即：
y′′=Pn(x)y''=P_n(x) y′′=Pn​(x)
两次积分得：
y=∫Pn+1dx+C1x+C2y=\int P_{n+1}\,\mathrm dx+C_1x+C_2 y=∫Pn+1​dx+C1​x+C2​

简要(emmm)总结：q!=0? y : (q!=0 ? x.y : xx.y)

常微分方程部分应该多加练习。寻找一种求解的直觉。

4.3.3.2. eee指数型

核心思路仍然是待定系数：
f(x)=a⋅eαx,a,α∈R,a≠0f(x)=a\cdot e^{\alpha x}, a, \alpha\in\R,a\not=0f(x)=a⋅eαx,a,α∈R,a​=0
设y(x)=A⋅eαxy(x)=A\cdot e^{\alpha x}y(x)=A⋅eαx
原方程转化为：
(Aα2+pAα+qA)eαx=a⋅eαx(A\alpha^2+pA\alpha+qA)e^{\alpha x}=a\cdot e^{\alpha x} (Aα2+pAα+qA)eαx=a⋅eαx
（α\alphaα是特征根）case 1:
若α2+pα+q≠0\alpha^2+p\alpha+q\not=0α2+pα+q​=0，则
A∗=aα2+pα+q(1)A^*=\frac{a}{\alpha^2+p\alpha +q}\tag1 A∗=α2+pα+qa​(1)
A∗⋅eαxA^*\cdot e^{\alpha x}A∗⋅eαx为一个特解。

（α\alphaα是单特征根）case 2:
若α2+pα+q=0\alpha^2+p\alpha+q=0α2+pα+q=0，则。。。没那么简单

同样使用乘xxx的方法，y=Axeαxy=Axe^{\alpha x}y=Axeαx，原式即：
Aeαx[(α2+pα+q)x+(2α+p)]Aeαx(2α+p)=aeαxAe^{\alpha x}\left[(\alpha^2+p\alpha+q)x+(2\alpha+p)\right]\\ Ae^{\alpha x}\left(2\alpha+p\right)=ae^{\alpha x} Aeαx[(α2+pα+q)x+(2α+p)]Aeαx(2α+p)=aeαx
（α\alphaα不是重根）再分两种情况：2α+p≠02\alpha+p\not=02α+p​=0直接找到一个特解
a2α+p(2)\frac{a}{2\alpha+p}\tag22α+pa​(2)

（α\alphaα是重根）若2α+p=02\alpha+p=02α+p=0，则α\alphaα为λ2+pλ+q=0\lambda^2+p\lambda+q=0λ2+pλ+q=0的重根。此时AxeαxAxe^{\alpha x}Axeαx仍然不是满足条件的特解。我们还需要构建一个新的形式的解。

设Ax2eαxAx^2e^{\alpha x}Ax2eαx，代入得
Aeαx[(α2+pα+q)x2+(4α+2p)x+2]=2Aeαx=aeαxAe^{\alpha x}\left[(\alpha^2+p\alpha +q)x^2+(4\alpha+2p)x+2\right]\\ =2Ae^{\alpha x}=ae^{\alpha x} Aeαx[(α2+pα+q)x2+(4α+2p)x+2]=2Aeαx=aeαx

即
y=a2x2eαx(3)y=\frac{a}{2}x^2e^{\alpha x}\tag3y=2a​x2eαx(3)
为一个特解.

4.3.3.3. 三角型

设y=Acos⁡βx+Bsin⁡βxy=A\cos\beta x+B\sin\beta xy=Acosβx+Bsinβx

代入得：
(−Aβ2+pBβ+qA)cos⁡βx+(−Bβ2−pAβ+Bq)sin⁡βx=acos⁡βx+bsin⁡βx(-A\beta^2+pB\beta+qA)\cos\beta x+(-B\beta^2-pA\beta+Bq)\sin\beta x=a\cos\beta x+b\sin \beta x (−Aβ2+pBβ+qA)cosβx+(−Bβ2−pAβ+Bq)sinβx=acosβx+bsinβx

转化为方程
(q−β2)A+pβB=a−pβA+(q−β2)B=b(q-\beta^2)A+p\beta B=a\\ -p\beta A+(q-\beta^2)B=b (q−β2)A+pβB=a−pβA+(q−β2)B=b
由Cramer法则
Δ=(q−β)2+p2β2≠0\Delta =(q-\beta)^2+p^2\beta^2\not=0 Δ=(q−β)2+p2β2​=0
case 1：
如果Δ≠0\Delta \not=0Δ​=0，那么有唯一解。记作
A∗cos⁡βx+B∗sin⁡βxA^*\cos\beta x+B^*\sin\beta x A∗cosβx+B∗sinβx
case2：
（βi\beta iβi是特征根）如果q−β2q-\beta^2q−β2与pβp\betapβ同时为0
设y=x(Acos⁡βx+Bsin⁡βx)y=x(A\cos\beta x+B\sin\beta x)y=x(Acosβx+Bsinβx)
代入得：
{pA+2βB=a−2βA+pB=b\begin{cases} pA+2\beta B=a\\ -2\beta A+pB=b \end{cases} {pA+2βB=a−2βA+pB=b​
由Cramer法则，方程有唯一解。
y=x(A∗cos⁡βx+B∗sin⁡βx)y=x(A^*\cos\beta x+B^*\sin\beta x) y=x(A∗cosβx+B∗sinβx)
就是特解。

4.3.3.4. 其他组合

除了以上三种，我们有以下可解的组合情况
Pn(x)eαx(αcos⁡βx+bsin⁡βx)P_n(x)e^{\alpha x}(\alpha \cos\beta x+b\sin\beta x) Pn​(x)eαx(αcosβx+bsinβx)
如果α±iβ\alpha\pm i\betaα±iβ不是特征根，eαx[Qn(x)cos⁡βx+Rn(x)sin⁡βx]e^{\alpha x}[Q_n(x)\cos\beta x+R_n(x)\sin\beta x]eαx[Qn​(x)cosβx+Rn​(x)sinβx]，如果是特征根（退化情况就乘xxx即可），即：eαx[Qn(x)cos⁡βx+Rn(x)sin⁡βx]e^{\alpha x}[Q_n(x)\cos\beta x+R_n(x)\sin\beta x]eαx[Qn​(x)cosβx+Rn​(x)sinβx]

这个结合分治定理便可以解决大量的非齐次问题了。

比如这个例子
y′′+y=xcos⁡2x+sin⁡xy''+y=x\cos 2x+\sin x y′′+y=xcos2x+sinx
利用分治定理我们就可以化成我们熟悉的形式了。

补充：
为什么最后的总结我们将判别条件利用是否为特征根表达? 做题真的很方便啊！

4.4. 其他重要方法

4.4.1. 常数变易法

当右式不是以上讨论的几种类型的时候，我们需要使用先哲发现的伟大通解法：常数变易法。

常数变易法是解线性微分方程行之有效的一种方法。它是拉格朗日十一年的研究成果，我们所用仅是他的结论，并无过程。

常数变易，顾名思义就是将本来是常数的部分变成函数。即：
设
y(x)=C1(x)φ1(x)+C2(x)⋅φ2(x)y(x)=C_1(x)\varphi_1(x)+C_2(x)\cdot\varphi_2(x)y(x)=C1​(x)φ1​(x)+C2​(x)⋅φ2​(x)
是
y′′(x)+p(x)y′(x)+q(x)y(x)=f(x)(*)y''(x)+p(x)y'(x)+q(x)y(x)=f(x)\tag{*} y′′(x)+p(x)y′(x)+q(x)y(x)=f(x)(*)
的一个解。
y′(x)=C1(x)φ1′(x)+C2(x)φ2′(x)+C1′(x)φ1(x)+C2′(x)φ2(x)\begin{aligned} y'(x)&=C_1(x)\varphi_1'(x)+C_2(x)\varphi'_2(x)\\ &+C_1'(x)\varphi_1(x)+C_2'(x)\varphi_2(x) \end{aligned}y′(x)​=C1​(x)φ1′​(x)+C2​(x)φ2′​(x)+C1′​(x)φ1​(x)+C2′​(x)φ2​(x)​

令：
C1′(x)φ(x)+C2′(x)φ2(x)=0(1)C_1'(x)\varphi(x)+C_2'(x)\varphi_2(x)=0\tag1 C1′​(x)φ(x)+C2′​(x)φ2​(x)=0(1)
便有
y′(x)=C1(x)φ1′(x)+C2(x)φ2′(x)y'(x)=C_1(x)\varphi_1'(x)+C_2(x)\varphi_2'(x) y′(x)=C1​(x)φ1′​(x)+C2​(x)φ2′​(x)
再求导
y′′=C1(x)φ1′′(x)+C2(x)φ2′′(x)+C1′(x)φ1′(x)+C2′(x)φ2′(x)y''=C_1(x)\varphi_1''(x)+C_2(x)\varphi_2''(x)+C_1'(x)\varphi_1'(x)+C_2'(x)\varphi_2'(x) y′′=C1​(x)φ1′′​(x)+C2​(x)φ2′′​(x)+C1′​(x)φ1′​(x)+C2′​(x)φ2′​(x)
注意到φ1(x)\varphi_1(x)φ1​(x)与φ2(x)\varphi_2(x)φ2​(x)是齐次方程的解：
C1′(x)φ1′(x)+C2′(x)φ2′(x)=f(x)(2)C_1'(x)\varphi'_1(x)+C_2'(x)\varphi_2'(x)=f(x)\tag2 C1′​(x)φ1′​(x)+C2′​(x)φ2′​(x)=f(x)(2)
联立(1)(2)，是一个Wronski行列式不为零的方程组，可解C1′,C2′C_1',C_2'C1′​,C2′​
关于常数变易法，我们特别注意的公式是：
{C1′(x)φ1(x)+C2′(x)φ2(x)=0C1′(x)φ1′(x)+C2′(x)φ2′(x)=f(x)\begin{cases} C_1'(x)\varphi_1(x)+C_2'(x)\varphi_2(x)=0\\ C_1'(x)\varphi_1'(x)+C_2'(x)\varphi_2'(x)=f(x) \end{cases} {C1′​(x)φ1​(x)+C2′​(x)φ2​(x)=0C1′​(x)φ1′​(x)+C2′​(x)φ2′​(x)=f(x)​

4.4.2. Euler方程

一类可以化成常系数方程的线性变系数方程。（但方法固定，考察并不多）
a0xny(n)+a1xn−1y(n−1)+⋯+an−1xy′+any=0a_0x^ny^{(n)}+a_1x^{n-1}y^{(n-1)}+\cdots+a_{n-1}xy'+a_ny=0 a0​xny(n)+a1​xn−1y(n−1)+⋯+an−1​xy′+an​y=0
其中ai(i=0,1,2⋯)a_i(i=0,1,2\cdots)ai​(i=0,1,2⋯)为常数

方法：
代入x=etx=e^tx=et，然后就可以化成常系数线性微分方程，通解公式为：
xky(k)=D(D−1)⋯(D−k+1)yx^ky^{(k)}=D(D-1)\cdots(D-k+1)y xky(k)=D(D−1)⋯(D−k+1)y
其中DDD为微分算符ddt\frac{\mathrm d}{\mathrm dt}dtd​，DDD的次数直接对应为微分方程中对应项的次数。

理解过程可以使用算子法：
∂∂x=1x∂∂t\frac{\partial}{\partial x}=\frac{1}{x}\frac{\partial}{\partial t} ∂x∂​=x1​∂t∂​

4.4.3. 常系数线性微分方程组

4.4.3.1. 基本概念

• 常数的独立性利用Jacobi行列式检验。
• 特解（代入后使得每一个方程成为恒等式的函数组），通解（带有n个独立的任意常数的函数组），向量函数
• 唯一性定理：只要n个初值给定，就存在唯一

4.4.3.2. 几种低阶的解法

4.4.3.2.1. 二阶

代入法，保留一个，消去一个从而成为二阶线性方程

4.4.3.2.2. 三阶

讨论是否有两个元的系数成比例：
利用Wronski行列式

• 有两个变量的系数成比例：先解另一个变量的二阶线性方程
• 均不成比例：不妨设x,y,zx,y,zx,y,z，那么先用xxx及其导数表示y,zy,zy,z然后得到xxx的三阶方程。

Part 9 常微分方程基础与常见方程求解法相关推荐

1. 已知基础解系反求有效方程(矩阵)

   已知基础解系反求有效方程(矩阵) @(数学) 这个是很有趣的推导过程,原理需要弄清楚. 即:已知Ax = 0的基础解系,由Ax = 0的系数行向量与解向量的关系可以反过来求解A. 具体推导如下: 齐次 ...

2. 常见的求最大公约数的方法

   最大公约数 最大公约数和最小公倍数求解,常用的方法是短除法进行因式分解,然后最大公约数是所有公共因子的乘积,最小公倍数是所有因子的乘积. 本质上求最小公倍数就是求最大公倍数:x=m*a, y=m*b: ...

3. 三个点在同一个半圆的概率_【国际数学竞赛】列方程求概率

   在国际数学竞赛中概率的考察往往也是古典概型和几何概型,但有些题目却找不到样本空间,只有事件的交错与转化,看着很是复杂,这里就介绍一种方法--通过事件间的关系列出方程组求出概率值.下面通过2017年AM ...

4. UA PHYS515 电磁理论I 麦克斯韦方程组基础3 麦克斯韦方程的势能形式

   UA PHYS515 电磁理论I 麦克斯韦方程组基础3 麦克斯韦方程的势能形式 定义电磁场的potential 改写Maxwell方程 上一讲我们基于实验定律导出了真空中电磁场的Maxwell方程: ...

5. 关于方程求根的解决方案

   对于方程求根主要的思想主要采取迭代的思想,通过条件判断,循环执行直到满足条件以后直接跳出循环输出 下面以x-cos(x)=0:为例  采用do-while 循环,输出Root: #include &q ...

6. 1.5编程基础之循环控制 32 求分数序列和 python

   http://noi.openjudge.cn/ch0105/32/ """ 1.5编程基础之循环控制 32 求分数序列和 http://noi.openjudge.cn ...

7. 1.5编程基础之循环控制 01 求平均年龄

   题目链接 http://noi.openjudge.cn/ch0105/01/ #include<iostream> using namespace std; int main() {in ...

8. 拼题a答案c语言计算整数各位数字之和,这是小白刷的基础题，跪求大佬帮忙！！！题目:读入一个自然数n...

   该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 这是小白刷的基础题,跪求大佬帮忙!!! 题目:读入一个自然数n,计算其各位数字之和,用汉语拼音写出和的每一位数字. 输入格式:每个测试输入包含1个测试用例 ...

9. 【机器学习】数值分析02——任意方程求根

   任意方程求根 全文目录 (博客园)机器学习 (Github)MachineLearning Math 1.简介 方程和函数是代数数学中最为重要的内容之一,从初中直到大学,我们都在研究着方程与函数,甚至 ...

最新文章

1. Python中fnmatch模块的使用
2. MS/OR国际期刊排名
3. 左耳朵耗子：程序员如何把控自己的职业？
4. linux命令的导入,[导入]Linux基本命令
5. python的8种标准数据类型有哪些_Python的八种数据类型
6. [转载] Python——函数练习（包括简单递归）
7. jbpm的流程变量与任务变量的存取与更新
8. 在linux中安装Jupyter notebook
9. 服装设计与工程_百度百科
10. R语言威尔科克森(Wilcoxon)分布
11. Layer Tree 绘制
12. 用python实现相声、评书自动下载（六）
13. JS中数组(Array)、Json对象长度(length)获取方法
14. windows进注册表快捷键
15. cocos2d-x分享系统
16. linux shell 在后台运行命令
17. 小米笔记本适合计算机专业吗,小米电脑笔记本好吗应该选哪个,小米入手使用感受...
18. cocoapods——更新
19. 基于STC89C52单片机的智能车控制系统设计
20. 简单的爱，简单的幸福

热门文章

1. WinXP系统怎么重装？
2. 胖哈勃 web--NewSql Mysql 8 注入
3. CRM销售管理系统：用数据共享提升客户满意度
4. CUDA RuntimeError: CUDA error: an illegal memory access was encountered
5. scada系统web服务器,基于IEC61970的Web-SCADA系统服务器后台的设计与实现
6. Clickhouse单机部署以及从mysql增量同步数据
7. c罗图片带字经典语言,c罗足球名言短句 c罗的那些霸气名言
8. OpenCV 图像梯度 ：cv2.Sobel()，cv2.Schar()，cv2.Laplacian() + 数据类型设置：cv2.CV_8U，cv2.CV_16S，cv2.CV_64F
9. Apache Tomcat 下载网站惊现校内（xiaonei）镜像
10. 医疗机构如何确保符合HIPAA标准 保障网络安全