第 8 讲:常系数齐次方程与 Euler 方程

常系数线性齐次方程之所以能用代数求解,是因为指数函数求导后只差一个常数倍。代入 erxe^{rx} 后,微分算子就变成关于 rr 的多项式。

特征方程

对

a0y(n)+a1y(n−1)+⋯+any=0,a0≠0,a_0y^{(n)}+a_1y^{(n-1)}+\cdots+a_ny=0, \qquad a_0\ne0,

试 y=erxy=e^{rx},得到

P(r)=a0rn+a1rn−1+⋯+an=0.P(r)=a_0r^n+a_1r^{n-1}+\cdots+a_n=0.

每个特征根按照类型生成基解。

不同实根

根 r1,…,rnr_1,\dots,r_n 不同时,

y=C1er1x+⋯+Cnernx.y=C_1e^{r_1x}+\cdots+C_ne^{r_nx}.

实重根

若 rr 是 mm 重根,对应

erx,xerx,…,xm−1erx.e^{rx},xe^{rx},\dots,x^{m-1}e^{rx}.

乘 xkx^k 是为了补出线性无关的解,不能只把同一个指数项重复写 mm 次。

复根

若 α±iβ\alpha\pm i\beta 是一对复根,对应两个实解

eαxcos⁡βx,eαxsin⁡βx.e^{\alpha x}\cos\beta x,\qquad e^{\alpha x}\sin\beta x.

若这对复根重数为 mm,还要依次乘 1,x,…,xm−11,x,\dots,x^{m-1}。

用算子因式分解理解

令 D=d/dxD=d/dx,方程写作

P(D)y=0.P(D)y=0.

若 P(r)=(r−r1)m1⋯(r−rk)mkP(r)=(r-r_1)^{m_1}\cdots(r-r_k)^{m_k},则

P(D)=(D−r1)m1⋯(D−rk)mk.P(D)=(D-r_1)^{m_1}\cdots(D-r_k)^{m_k}.

这个表示也解释了为什么重根产生 xkerxx^ke^{rx},并为非齐次方程中的共振规则做准备。

Euler–Cauchy 方程

典型形式

x2y′′+axy′+by=0,x≠0.x^2y''+axy'+by=0,\qquad x\ne0.

可以试幂函数 y=xmy=x^m,得到指标方程

m(m−1)+am+b=0.m(m-1)+am+b=0.

在 x>0x>0 上也可令

t=ln⁡x,Y(t)=y(et).t=\ln x,\qquad Y(t)=y(e^t).

利用

xy′=Y′,x2y′′=Y′′−Y′,xy'=Y',\qquad x^2y''=Y''-Y',

把 Euler 方程变成常系数方程。对 x<0x<0 可用 t=ln⁡∣x∣t=\ln|x| 分区间处理,不能让解跨过奇点 x=0x=0。

若指标根为重根 mm,基解为

xm,xmln⁡∣x∣.x^m,\qquad x^m\ln|x|.

复根 α±iβ\alpha\pm i\beta 对应

∣x∣αcos⁡(βln⁡∣x∣),∣x∣αsin⁡(βln⁡∣x∣).|x|^\alpha\cos(\beta\ln|x|),\qquad |x|^\alpha\sin(\beta\ln|x|).

边值问题与特征值

例如

y′′+λy=0,y(0)=y(L)=0y''+\lambda y=0,\qquad y(0)=y(L)=0

并非任意 λ\lambda 都有非零解。根据 λ\lambda 的正、零、负三种情形分别写通解,边界条件会筛出离散的特征值

λn=(nπL)2,n=1,2,…\lambda_n=\left(\frac{n\pi}{L}\right)^2,\qquad n=1,2,\dots

及对应特征函数 sin⁡(nπx/L)\sin(n\pi x/L)。这说明边值条件不是简单“解两个常数”,有时会决定参数本身。

快速检查

  • 特征多项式次数应等于方程阶数;
  • 根的代数重数总和应为 nn;
  • 每个实根贡献其重数个解,每对复根贡献两倍重数个实解;
  • Euler 方程的解区间不能跨 x=0x=0;
  • 写完通解后数一遍独立常数,并代回低阶示例检查符号。

评论