Part 9 常微分方程基础与常见方程求解法

文章目录

1. 常微分方程是什么

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

1.1. 数学史引入

1.1.1. Newton第二定律:

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

1.1.2. 万有引力

利用Newton第二定律
地球所受的万有引力可以近似表示成:
G m m m e ∣ x m − x e ∣ 2 x m − x e ∣ x m − x e ∣ + G m s m e ∣ x s − x e ∣ 2 x s − x e ∣ x s − x e ∣ = m e d 2 x d t 2 G\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} Gxmxe2mmmexmxexmxe+Gxsxe2msmexsxexsxe=medt2d2x

1.1.3. 常微分方程的贡献者

  • Bernoulli家族
  • Euler
  • Lagarange
  • Gauss

1.1.4. 一个突破

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

1.2. 定义

1.2.1. 微分方程

含有未知函数及其至 n n n阶导数的方程
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))xR(1)
称为一个 n n n阶常微分方程。其中各个函数都在对应相同的 x x x处取值。

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

例:

  1. d y d x = y ( x − 1 ) \frac{\mathrm dy}{\mathrm dx}=y(x-1) dxdy=y(x1)不是微分方程,因为其中出现了非 x x x的变量(式中 y ( x − 1 ) y(x-1) y(x1) y y y x − 1 x-1 x1处的取值)
  2. d y d x = 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 φ有直到 n n n阶导数,且 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 = c 1 + c 2 t − 1 2 g t 2 , c 1 c 2 y=c_1+c_2t-\frac{1}{2}gt^2,c_1c_2 y=c1+c2t21gt2,c1c2为常数。

1.2.3. 通解

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

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

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

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

Jacobi矩阵
( ∂ φ c 1 ∂ φ c 2 ⋯ ∂ φ ∂ c n ∂ φ ′ ∂ c 1 ∂ φ ′ c 2 ⋯ ∂ φ ′ ∂ c n ⋮ ⋮ ⋱ ⋮ ∂ φ ( n − 1 ) c 1 ∂ φ ( n − 1 ) φ c 2 ⋯ ∂ φ ( n − 1 ) ∂ c n ) \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φ(n1)c2φc2φφc2φ(n1)cnφcnφcnφ(n1)
不为零可以判别微分满秩。从而得知这些常数不相关。

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

例如:对方程 y ′ ′ + y = 0 y''+y= 0 y+y=0我们可以写出特解 y = sin ⁡ x y=\sin x y=sinx y = cos ⁡ x y=\cos x y=cosx,通解 y = c 1 sin ⁡ x + c 2 cos ⁡ x y=c_1\sin x+c_2\cos x y=c1sinx+c2cosx
但反智的一点是,存在不能通过特解中常数取特殊值将特值表出的情况。这些解叫做“奇解”。

1.2.5. 初值问题

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

d 2 x d t 2 = − g , d x d y = − g t + c 2 , x = − 1 2 g t 2 + c 2 t + c 1 \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_1 dt2d2x=g,dydx=gt+c2,x=21gt2+c2t+c1
如果有一个初态 x = h , d x d t = 0 x=h,\frac{\mathrm dx}{\mathrm dt}=0 x=h,dtdx=0,那么 h = c 1 , d x d t = − g t + c 2 h=c_1,\frac{\mathrm dx}{\mathrm dt}=-gt+c_2 h=c1,dtdx=gt+c2解得: c 2 = 0 c_2=0 c2=0

y ( x 0 ) = y 0 , y ′ ( x 0 ) = y 0 ′ ⋯ y ( n − 1 ) = y 0 ( 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)=y0y(n1)=y0(n1)

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

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 )   d x = G ( x )   d x F(x)\,\mathrm dx=G(x)\,\mathrm dx F(x)dx=G(x)dx的形式,其中方程右侧是 x , y x,y xy分别的函数的乘积。

d y d x = f ( x ) g ( y ) d y g ( y ) = f ( x )   d x ∫ 1 g ( y )   d y = ∫ f ( x )   d x \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)dxg(y)1dy=f(x)dx
这一步先不考虑特殊情况即能不能除,找准主要思路。记
G ( y ) = ∫ 1 g ( y )   d y F ( x ) = ∫ f ( x )   d x G ( 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)1dyF(x)=f(x)dxG(y)=F(x)
这是一个代数方程

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

y ′ = α y ( α y'=\alpha y(\alpha y=αy(α是常数)
1 y   d y = α   d x ln ⁡ ∣ y ∣ = α x + C ∣ y ∣ = e c ⋅ e α x = C 1 ⋅ e α x y = 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) y1dy=αdxlny=αx+Cy=eceαx=C1eαxy=Ceαx(c=0)
加上特解 y ≡ 0 y\equiv0 y0,原方程的解可以表示成 y = C ⋅ e α x y=C\cdot e^{\alpha x} y=Ceαx

2.1.2. 可化为变量分离型

1 。 1^。 1形如
d y d x = f ( a x + b y + c ) \frac{\mathrm dy}{\mathrm dx}=f(ax+by+c) dxdy=f(ax+by+c)
的方程:
引入 u = a x + b y + c , u=ax+by+c, u=ax+by+c,
d u d x = a + b d y d x = a + b f ( a x + b y + c ) = a + b f ( 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+bdxdya+bf(ax+by+c)a+bf(u)
显然 u u u是可解的。 u = u ( x , c ~ ) u=u(x,\tilde{c}) u=u(x,c~)
这又是一个代数方程。

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

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

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

u = y x , y = x u , d y d x = d d x ( x u ) = u + x d u d x = 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+xdxduy=f(u) f(u)
由于 x ≠ 0 x\not=0 x=0 d u d x = f ( u ) − u x \frac{\mathrm du}{\mathrm dx}=\frac{f(u)-u}{x} dxdu=xf(u)u

3 ∘ 3^\circ 3另一类齐次
y = f ( a 1 x + b 1 y + c 1 a x + b y + c ) y=f\left(\frac{a_1x+b_1y+c_1}{ax+by+c}\right) y=f(ax+by+ca1x+b1y+c1)

  1. c = c 1 = 0 , y ′ = f ( a 1 x + b 1 y a x + b y ) = f ( a 1 + b 1 y x a + b y x ) = g ( y x ) 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+bya1x+b1y)=f(a+bxya1+b1xy)=g(xy)
  2. 若不全为0.那么做平移变换。

{ x = u + α y = v + β s . t .   a 1 α + b 1 β + c 1 = 0 a α + 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

d v d u = f ( a 1 u + b 1 v a u + b v ) \frac{\mathrm dv}{\mathrm du}=f\left(\frac{a_1u+b_1v}{au+bv}\right) dudv=f(au+bva1u+b1v)
利用Cramer法则,
∣ a 1 b 1 a b ∣ ≠ 0 \begin{vmatrix} a_1&b_1\\ a&b \end{vmatrix}\not=0 a1ab1b=0
时, α , β \alpha, \beta α,β存在且唯一,可解。
∣ a 1 b 1 a b ∣ = 0 \begin{vmatrix} a_1&b_1\\ a&b \end{vmatrix}=0 a1ab1b=0
时,设 a 1 a = b 1 b = k \frac{a_1}{a}=\frac{b_1}{b}=k aa1=bb1=k
y ′ = f ( a 1 x + b 1 y + c 1 a x + b y + c ) = f ( k ( a x + b y ) + c 1 a x + b y + 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+ca1x+b1y+c1)=f(ax+by+ck(ax+by)+c1)
现令
d u d x = a + b d y d x = a + b f ( k u + c 1 u + 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)
变量分离,又是可解的。
d y d x = x + y + 2 x − y + 1 \frac{\mathrm d y}{\mathrm dx}=\frac{x+y+2}{x-y+1} dxdy=xy+1x+y+2

2.2. 一阶线性方程

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

2.2.1. 一阶线性方程的判别

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

2.2.2. 一阶线性方程的计算

2.2.2.1. 1 ∘ 1^\circ 1 最简单的情况: q ( x ) = 0 q(x)=0 q(x)=0(齐次)

这个“齐次"和我们前述的齐次函数的齐次并不相同。
y ′ + p ( x ) y = 0 y'+p(x)y=0 y+p(x)y=0
(2)
d y d x = − p ( x ) y \frac{\mathrm dy}{\mathrm dx}=-p(x)y dxdy=p(x)y
首先我们发现 y = 0 y=0 y=0是一个奇解。然后移项积分:
l n ∣ y ∣ = − ∫ x 0 x P ( t )   d t + C ln|y|=-\int_{x_0}^xP(t)\,\mathrm dt+C lny=x0xP(t)dt+C
y = ± e c e ( − ∫ x 0 x P ( t )   d t ) y=\pm e^c e^{\left(-\int_{x_0}^xP(t)\,\mathrm dt\right)} y=±ece(x0xP(t)dt)
结合奇解
y = C ⋅ e ( − ∫ x 0 x P ( t )   d t ) , ∀ C y=C\cdot e^{\left(-\int_{x_0}^xP(t)\,\mathrm dt\right)},\forall C y=Ce(x0xP(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)的两个解,那么 ∀ k 1 , k 2 ∈ R , k 1 φ 1 + k 2 φ 2 \forall k_1,k_2\in\R, k_1\varphi_1+k_2\varphi_2 k1,k2R,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 + y 2 , 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^\circ 2非齐次
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=\Phi y=Φ是(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 ) + C e − ∫ x 0 x P ( t )   d t y=\varphi(x)+Ce^{-\int_{x_0}^xP(t)\,\mathrm dt} y=φ(x)+Cex0xP(t)dt
  • 由上述,这个方程的求解转化成求解(1)的一个特解。
2.2.2.2.3. 方法:常数变易法

将齐次方程的解中的常数变成一个函数表达式,待定为 u ( x ) u(x) u(x)
y ′ = u ′ ( x ) ⋅ e − ∫ x 0 x P ( t )   d t − u ( x ) P ( x ) e − ∫ x 0 x P ( t )   d t = Q − P ( x ) y y'=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)ex0xP(t)dtu(x)P(x)ex0xP(t)dt=QP(x)y
化简得:
u ′ ( x ) ⋅ e − ∫ x 0 x P ( t )   d t = Q u'(x)\cdot e^{-\int_{x_0}^xP(t)\,\mathrm dt}=Q u(x)ex0xP(t)dt=Q
u ′ ( x ) = Q ( x ) ⋅ e ∫ x 0 x P ( t )   d t u'(x)=Q(x)\cdot e^{\int_{x_0}^xP(t)\,\mathrm dt} u(x)=Q(x)ex0xP(t)dt
这个式子可不可积其实还有待讨论……但总之把它积起来就可以得到我们需要的 u ( x ) u(x) u(x)
y = [ ∫ x 0 x Q ( t ) e ∫ x 0 t P ( s )   d s   d t + C ] e − ∫ x 0 x P ( t )   d t y=\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=[x0xQ(t)ex0tP(s)dsdt+C]ex0xP(t)dt

2.2.2.2.4. 步骤整理:
  • step1:求齐次方程通解 C ⋅ y 1 ( x ) C\cdot y_1(x) Cy1(x)
  • step2: u ′ = Q ( x ) y 1 ( x ) u'=\frac{Q(x)}{y_1(x)} u=y1(x)Q(x)
  • step3:移项积分,整理得解

f ∈ C ( a , + ∞ ) f\in C(a, +\infty) fC(a,+) f ( x ) + f ′ ( x ) → 0 f(x)+f'(x)\to0 f(x)+f(x)0 f ( x ) → 0 f(x)\to0 f(x)0

另例:
f ( x ) ∈ C 1 [ 0 , + ∞ ) , a ( x ) ≥ C 0 > 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))=0 x+lim(f(x)+a(x)f(x))=0 lim ⁡ x → + ∞ f ( x ) = 0 \lim\limits_{x\to+\infty}f(x)=0 x+limf(x)=0

b ( x ) = f ′ ( x ) + a ( x ) f ( x ) , b ( x ) ∈ C 0 [ 0 , + ∞ ) , lim ⁡ x → + ∞ b ( x ) = 0 b(x)=f'(x)+a(x)f(x), b(x)\in C^0[0,+\infty),\lim\limits_{x\to+\infty} b(x)=0 b(x)=f(x)+a(x)f(x),b(x)C0[0,+),x+limb(x)=0,从 b b b中解出 f ( x ) f(x) f(x)
f ( x ) = C e − ∫ x 0 x a ( t )   d t + e − ∫ x 0 x a ( t )   d t ⋅ ∫ 0 x b ( s ) e ∫ 0 s a ( t )   d t   d s f(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)=Cex0xa(t)dt+ex0xa(t)dt0xb(s)e0sa(t)dtds
其中取初值条件,可以得 C = f ( 0 ) C=f(0) C=f(0)
f ( x ) = f ( 0 ) + ∫ 0 x b ( s ) e ∫ 0 s a ( t )   d t   d s e ∫ 0 x a ( t )   d t f(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)=e0xa(t)dtf(0)+0xb(s)e0sa(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)}=0 x+limf(x)=x+lima(x)b(x)=0

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

这是一个将非线性方程化成线性方程的好例子。
d y d x + 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 − α d y d x + P ⋅ y 1 − α = Q ( x ) y^{-\alpha}\frac{\mathrm dy}{\mathrm dx}+P\cdot y^{1-\alpha}=Q(x) yαdxdy+Py1α=Q(x)
z = 1 − α z=1-\alpha z=1α,
1 1 − α d z d x + P ( x ) z = Q ( x ) \frac{1}{1-\alpha}\frac{\mathrm dz}{\mathrm dx}+P(x)z=Q(x) 1α1dxdz+P(x)z=Q(x)
这是一个线性方程。

2.3. 全微分方程

2.3.1. 全微分方程的形式

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

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

∃ F ( x , y ) ∈ C 1 ( D ) ,   d F ( x , y ) = P   d x + Q   d y ⇒ d F ( 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+QdydF(x,y)=0

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

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

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

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

偏积分得
F = ∫ x 0 x P ( s , y )   d s + φ ( y ) F=\int_{x_0}^xP(s,y)\,\mathrm ds+\varphi(y) F=x0xP(s,y)ds+φ(y)
(1)
Q = ∂ ∂ y ∫ x 0 x P ( s , y )   d s + φ ′ ( y ) Q=\frac{\partial}{\partial y}\int_{x_0}^xP(s,y)\,\mathrm ds+\varphi'(y) Q=yx0xP(s,y)ds+φ(y)
由于光滑,求导积分可交换。
Q = ∫ x 0 x ∂ ∂ y P ( s , y )   d s + φ ′ ( y ) Q=\int_{x_0}^x\frac{\partial}{\partial y}P(s,y)\,\mathrm ds+\varphi'(y) Q=x0xyP(s,y)ds+φ(y)
再由充分条件得:
∫ x 0 x ∂ ∂ x Q ( s , y )   d s + φ ′ ( y ) = Q ( x , y ) \int_{x_0}^x\frac{\partial}{\partial x}Q(s,y)\,\mathrm ds+\varphi'(y)=Q(x,y) x0xxQ(s,y)ds+φ(y)=Q(x,y)
Q ( x , y ) − Q ( x 0 , 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 ( x 0 , y ) \varphi'(y)=Q(x_0,y) φ(y)=Q(x0,y)
φ = ∫ y 0 y Q ( x 0 , t )   d t \varphi=\int_{y_0}^yQ(x_0,t)\,\mathrm dt φ=y0yQ(x0,t)dt
结合(1)得:
F ( x , y ) = ∫ x 0 x P ( s , y )   d s + ∫ y 0 y Q ( x 0 , t )   d t F(x,y) =\int_{x_0}^xP(s,y)\,\mathrm ds+\int_{y_0}^yQ(x_0,t)\,\mathrm dt F(x,y)=x0xP(s,y)ds+y0yQ(x0,t)dt
说人话,就是说明了二型线与路径无关的条件是全微分方程直接用原函数是解的条件。

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

接下来要解决的问题是,如果 ∂ P ∂ y \frac{\partial P}{\partial y} yP ∂ Q ∂ x \frac{\partial Q}{\partial x} xQ不相等呢?

2.3.3. 积分因子

对以上的问题,基本想法是:
寻找一个(就够的) μ = μ ( x , y ) ≠ 0 \mu=\mu(x,y)\not=0 μ=μ(x,y)=0s.t.
μ P   d x + μ Q   d y = 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=xF,μQ=yF

结合光滑:
∂ ∂ 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) PyμQxμ=μ(xQyP)

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

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

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

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

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

如果方程
{ y ′ = f ( x , y ) y ( x 0 ) = y 0 (1) \begin{cases} y'=f(x,y)\tag1\\ y(x_0)=y_0 \end{cases} {y=f(x,y)y(x0)=y0(1)
f f f连续,且满足Lipschitz条件(斜率不大于某个值):
∃ L > 0 , ∀ ( x , y 1 ) , ( x , y 0 ) ∈ D \exist L>0,\forall(x,y_1),(x,y_0)\in D L>0,(x,y1),(x,y0)D
∣ f ( x , y 0 ) − f ( x , y 1 ) ∣ ≤ L ∣ y 1 − y 0 ∣ |f(x,y_0)-f(x,y_1)|\leq L|y_1-y_0| f(x,y0)f(x,y1)Ly1y0
那么其在 ∣ x − x 0 ∣ ≤ h |x-x_0|\leq h xx0h解存在且唯一
h = min ⁡ { a , b M } 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}Dmaxf(x,y)

  • 注1: h h h是越界保护。其中 M M M通过(1)中代换得到的: y ′ ≤ M y'\leq M yM

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

Step3:
序列一致收敛,故存在(证明略

Step4:
解的唯一性判定
∣ φ ( x ) − ψ ( x ) ∣ ≤ ∣ ∫ x 0 x f ( s , φ ( s ) ) − f ( s , ψ ( s ) )   d s ∣ ≤ L ∫ x 0 x ∣ φ ( s ) − ψ ( s ) ∣   d s |\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)x0xf(s,φ(s))f(s,ψ(s))dsLx0xφ(s)ψ(s)ds
∣ φ ( x ) − ψ ( x ) ∣ ≤ L M 0 ∣ x − x 0 ∣ |\varphi(x)-\psi(x)|\leq LM_0|x-x_0| φ(x)ψ(x)LM0xx0
故(2)代入(1)得
∣ φ ( x ) − ψ ( x ) ∣ ≤ L 2 M 0 2 ∣ x − x 0 ∣ 2 |\varphi(x)-\psi(x)|\leq \frac{L^2M_0}{2}|x-x_0|^2 φ(x)ψ(x)2L2M0xx02
累次代入,得:
∣ φ ( x ) − ψ ( x ) ∣ ≤ L n M 0 n ! h n |\varphi(x)-\psi(x)|\leq\frac{L^nM_0}{n!}h^n φ(x)ψ(x)n!LnM0hn
显然,
φ ( x ) ≡ ψ ( x ) \varphi(x)\equiv\psi(x) φ(x)ψ(x)

注1(Cauchy条件)
f ∈ C 1 ( D ) f\in C^1(D) fC1(D)
f ( x , y 1 ) − f ( x , y 2 ) = ∂ f ∂ y ( x , ξ ) ( y 1 − y 0 ) 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)=yf(x,ξ)(y1y0)
显然这是满足Lipschitz条件的。

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

推论:对 y ′ = f ( x , y ) , ( x , y ) ∈ D y'=f(x,y),(x,y)\in D y=f(x,y),(x,y)D,若函数 f ( x , y ) f(x,y) f(x,y) f y f_y fy在区域上连续,那么过D内任一点 ( x 0 , y 0 ) (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. 可降阶的方程

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

常用 不显含 x x x的二阶方程(通常设题比较刁钻,只能用此法解)
F ( y , y ′ , y ′ ′ ) = 0 F(y,y',y'')=0 F(y,y,y)=0

P = y ′ = d y d x P=y'=\frac{\mathrm dy}{\mathrm dx} P=y=dxdy
y ′ ′ = d d x y ′ = d d y y ′ d y d x = d P d y P y''=\frac{\mathrm d}{\mathrm dx}y'=\frac{\mathrm d}{\mathrm dy}y'\frac{\mathrm dy}{\mathrm dx}=\frac{\mathrm dP}{\mathrm dy}P y=dxdy=dydydxdy=dydPP
原式可以转化为
F 1 ( y , P , d P d y ) = 0 F_1(y,P,\frac{\mathrm dP}{\mathrm dy})=0 F1(y,P,dydP)=0
这是一阶线性方程。
所得解 P = d y d x P=\frac{\mathrm dy}{\mathrm dx} P=dxdy又是变量分离的。

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

这是一个快捷化简。

4.2. 函数的线性相关

4.2.1. 判据——W行列式

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

g 1 ( x ) , g 2 ( x ) ∈ D ( a , b ) , ∃ c 1 , c 2 , s . t .   ┐ ( c 1 = 0 ∧ c 2 = 0 ) ∧ c 1 g 1 + c 2 g 2 ≡ 0 g_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\equiv0 g1(x),g2(x)D(a,b),c1,c2,s.t.(c1=0c2=0)c1g1+c2g20则称 g 1 , g 2 g_1,g_2 g1,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 ′ ∣ ≡ 0 W=\begin{vmatrix} \varphi_1&\varphi_2\\ \varphi_1'&\varphi_2' \end{vmatrix}\equiv0 W=φ1φ1φ2φ20
证明:必要性对(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 ( x 0 ) e − ∫ x 0 x p ( s )   d s ( L i o u v i l l e 公 式 ) w(x)=w(x_0)e^{-\int_{x_0}^xp(s)\,\mathrm ds} (Liouville公式) w(x)=w(x0)ex0xp(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)
那么所有解都可以经由这两个解表出。

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

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 = c 1 ψ 1 ( x ) + c 2 ψ ( 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 ′ ′ + p y + q y = f 1 y''+py+qy=f_1 y+py+qy=f1
的解,
y = f 2 ( x ) y=f_2(x) y=f2(x)
y ′ ′ + p y ′ + q y = f 2 y''+py'+qy=f_2 y+py+qy=f2
的解。

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

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

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

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

首先解决齐次情形:
这里考虑到 e e e指数函数的求导性质优良,我们优先考虑形如 y = e λ x y=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 λ x e^{\lambda x} eλx是解 ⇔ λ 2 + p λ + q = 0 \Leftrightarrow \lambda^2+p\lambda+q=0 λ2+pλ+q=0

case 1:特征根相异时 y 1 = e λ 1 x , y 2 = e λ 2 x y_1=e^{\lambda_1 x},y_2=e^{\lambda_2x} y1=eλ1x,y2=eλ2x是两个线性无关的解
W ( x ) = ∣ e λ 1 x e λ 2 x λ 1 e λ 1 x λ 2 e λ 2 x ∣ = ( λ 2 − λ 1 ) e ( λ 1 + λ 2 ) x W(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λ1xλ1eλ1xeλ2xλ2eλ2x=(λ2λ1)e(λ1+λ2)x

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

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

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

y 1 ∗ ( x ) = e ( α + i β ) x = e ( α x ( cos ⁡ β x + i sin ⁡ β x ) y 2 ∗ ( x ) = e ( α − i β ) x = e α x ( cos ⁡ β x − i sin ⁡ β 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βxisinβx)
消虚数,取得两个实函数解基:
y 1 ( x ) = 1 2 ( y 1 ∗ + y 2 ∗ ) = e α x cos ⁡ β x y 2 ( x ) = 1 2 i ( y 1 ∗ − y 2 ∗ ) = e α x sin ⁡ β x y_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(y1y2)=eαxsinβx
由于 β > 0 \beta >0 β>0故而Wronski行列式为 β e 2 α x ≠ 0 \beta e^{2\alpha x}\not=0 βe2αx=0

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

特征方程判别式特征根通解形式
Δ > 0 \Delta >0 Δ>0两个相异实根 λ 1 , λ 2 \lambda_1, \lambda_2 λ1,λ2 C 1 e λ 1 x + C 2 e λ 2 x C_1e^{\lambda_1 x}+C_2e^{\lambda_2x} C1eλ1x+C2eλ2x
Δ = 0 \Delta = 0 Δ=0二重根 λ 1 \lambda_1 λ1 ( C 1 + C 2 x ) ⋅ e λ 1 x (C_1+C_2x)\cdot e^{\lambda_1}x (C1+C2x)eλ1x
Δ < 0 \Delta <0 Δ<0 λ 1 , 2 = α ± i β \lambda_{1,2}=\alpha\pm i\beta λ1,2=α±iβ e α x ( C 1 cos ⁡ β x + C 2 sin ⁡ β x ) e^{\alpha x}(C_1\cos\beta x+C_2\sin\beta x) eαx(C1cosβx+C2sinβx)
  • 从某个角度上说,非重根的所有情况都是第三种情形的退化情况。
  • 注意这个 β x \beta x βx是放在一起写的。这是 e i x e^{i x} eix的De Moivre展开。

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

特征根通解
单重实根 λ \lambda λ e λ x e^{\lambda x} eλx
k重实根( k > 1 k>1 k>1) λ \lambda λ e λ x , x e λ x , ⋯   , x k − 1 e λ x e^{\lambda x},xe^{\lambda x}, \cdots, x^{k-1}e^{\lambda x} eλx,xeλx,,xk1eλx
单共轭复根 λ 1 , 2 = α + i β \lambda_{1,2}=\alpha+i\beta λ1,2=α+iβ e α x cos ⁡ β x , e α x sin ⁡ β x e^{\alpha x}\cos\beta x,e^{\alpha x}\sin\beta x eαxcosβx,eαxsinβx
n n n重共轭复根 λ 1 , 2 = α + i β ( n > 1 ) \lambda_{1,2}=\alpha+i\beta(n>1) λ1,2=α+iβ(n>1) e α x cos ⁡ β x , e α x sin ⁡ β x , ⋯   , x n − 1 e α x cos ⁡ β x , x m − 1 sin ⁡ β x e^{\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 x eαxcosβx,eαxsinβx,,xn1eαxcosβx,xm1sinβx

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

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

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

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

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

4.3.3.1. 多项式型

利用待定系数法,设解为:
Q n ( x ) = b 0 x n + b 1 x n − 1 + ⋯ + b n , b 0 ≠ 0 Q_n(x)=b_0x^n+b_1x^{n-1}+\cdots +b_n, b_0\not=0 Qn(x)=b0xn+b1xn1++bn,b0=0
(0不是特征根)case 1:
q ≠ 0 q\not=0 q=0时,左右都是 n n n次多项式。比较各次项系数,解出 b i b_i bi

(0是单特征根)case 2:
考虑 Q ( x ) = x Q n ( x ) Q(x)=xQ_n(x) Q(x)=xQn(x),代入得:
Q ′ ( x ) = Q n ( x ) + x Q n ′ ( x ) Q ′ ′ ( x ) = 2 Q n ′ ( x ) + x Q n ′ ′ ( 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)
整理出 n n n次项作为首项的式子:
p ( x Q n ′ + Q n ) + ( x Q n ′ ′ + 2 Q n ′ ) = P n p(xQ_n'+Q_n)+(xQ_n''+2Q_n')=P_n p(xQn+Qn)+(xQn+2Qn)=Pn
如果 p ≠ 0 p\not=0 p=0退化成case 1中的讨论情形。

(0是重特征根)case 3:
对case 3中的 p = 0 p=0 p=0情形,考虑 Q n ( x ) = x 2 Q n Q_n(x)=x^2Q_n Qn(x)=x2Qn可以直接代入对比系数(北大P188)。

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

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

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

4.3.3.2. e e e指数型

核心思路仍然是待定系数:
f ( x ) = a ⋅ e α x , a , α ∈ R , a ≠ 0 f(x)=a\cdot e^{\alpha x}, a, \alpha\in\R,a\not=0 f(x)=aeαx,a,αR,a=0
y ( x ) = A ⋅ e α x y(x)=A\cdot e^{\alpha x} y(x)=Aeαx
原方程转化为:
( A α 2 + p A α + q A ) 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=aeα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 α x A^*\cdot e^{\alpha x} Aeαx为一个特解。

α \alpha α是单特征根)case 2:
α 2 + p α + q = 0 \alpha^2+p\alpha+q=0 α2+pα+q=0,则。。。没那么简单

同样使用乘 x x x的方法, y = A x e α x y=Axe^{\alpha x} y=Axeαx,原式即:
A e α x [ ( α 2 + p α + q ) x + ( 2 α + p ) ] A e α x ( 2 α + p ) = a e α x Ae^{\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 ≠ 0 2\alpha+p\not=0 2α+p=0直接找到一个特解
a 2 α + p (2) \frac{a}{2\alpha+p}\tag2 2α+pa(2)

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

A x 2 e α x Ax^2e^{\alpha x} Ax2eαx,代入得
A e α x [ ( α 2 + p α + q ) x 2 + ( 4 α + 2 p ) x + 2 ] = 2 A e α x = a e α x Ae^{\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 = a 2 x 2 e α x (3) y=\frac{a}{2}x^2e^{\alpha x}\tag3 y=2ax2eαx(3)
为一个特解.

4.3.3.3. 三角型

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

代入得:
( − A β 2 + p B β + q A ) cos ⁡ β x + ( − B β 2 − p A β + B q ) sin ⁡ β x = a cos ⁡ β x + b sin ⁡ β 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β2pAβ+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=apβA+(qβ2)B=b
由Cramer法则
Δ = ( q − β ) 2 + p 2 β 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 ⁡ β x A^*\cos\beta x+B^*\sin\beta x Acosβx+Bsinβx
case2:
β i \beta i βi是特征根)如果 q − β 2 q-\beta^2 qβ2 p β p\beta pβ同时为0
y = x ( A cos ⁡ β x + B sin ⁡ β x ) y=x(A\cos\beta x+B\sin\beta x) y=x(Acosβx+Bsinβx)
代入得:
{ p A + 2 β B = a − 2 β A + p B = b \begin{cases} pA+2\beta B=a\\ -2\beta A+pB=b \end{cases} {pA+2βB=a2βA+pB=b
由Cramer法则,方程有唯一解。
y = x ( A ∗ cos ⁡ β x + B ∗ sin ⁡ β x ) y=x(A^*\cos\beta x+B^*\sin\beta x) y=x(Acosβx+Bsinβx)
就是特解。

4.3.3.4. 其他组合

除了以上三种,我们有以下可解的组合情况
P n ( x ) e α x ( α cos ⁡ β x + b sin ⁡ β 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 [ Q n ( x ) cos ⁡ β x + R n ( 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],如果是特征根(退化情况就乘 x x x即可),即: e α x [ Q n ( x ) cos ⁡ β x + R n ( 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 = x cos ⁡ 2 x + sin ⁡ x y''+y=x\cos 2x+\sin x y+y=xcos2x+sinx
利用分治定理我们就可以化成我们熟悉的形式了。

补充:
为什么最后的总结我们将判别条件利用是否为特征根表达? 做题真的很方便啊!
在这里插入图片描述

4.4. 其他重要方法

4.4.1. 常数变易法

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

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

常数变易,顾名思义就是将本来是常数的部分变成函数。即:

y ( x ) = C 1 ( x ) φ 1 ( x ) + C 2 ( 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 ) = C 1 ( x ) φ 1 ′ ( x ) + C 2 ( x ) φ 2 ′ ( x ) + C 1 ′ ( x ) φ 1 ( x ) + C 2 ′ ( 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)

令:
C 1 ′ ( x ) φ ( x ) + C 2 ′ ( 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 ) = C 1 ( x ) φ 1 ′ ( x ) + C 2 ( 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 ′ ′ = C 1 ( x ) φ 1 ′ ′ ( x ) + C 2 ( x ) φ 2 ′ ′ ( x ) + C 1 ′ ( x ) φ 1 ′ ( x ) + C 2 ′ ( 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)是齐次方程的解:
C 1 ′ ( x ) φ 1 ′ ( x ) + C 2 ′ ( 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行列式不为零的方程组,可解 C 1 ′ , C 2 ′ C_1',C_2' C1,C2
关于常数变易法,我们特别注意的公式是:
{ C 1 ′ ( x ) φ 1 ( x ) + C 2 ′ ( x ) φ 2 ( x ) = 0 C 1 ′ ( x ) φ 1 ′ ( x ) + C 2 ′ ( 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方程

一类可以化成常系数方程的线性变系数方程。(但方法固定,考察并不多)
a 0 x n y ( n ) + a 1 x n − 1 y ( n − 1 ) + ⋯ + a n − 1 x y ′ + a n y = 0 a_0x^ny^{(n)}+a_1x^{n-1}y^{(n-1)}+\cdots+a_{n-1}xy'+a_ny=0 a0xny(n)+a1xn1y(n1)++an1xy+any=0
其中 a i ( i = 0 , 1 , 2 ⋯   ) a_i(i=0,1,2\cdots) ai(i=0,1,2)为常数

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

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

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 , z x,y,z x,y,z,那么先用 x x x及其导数表示 y , z y,z y,z然后得到 x x x的三阶方程。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值