常系数线性齐次方程之所以能用代数求解,是因为指数函数求导后只差一个常数倍。代入 erx 后,微分算子就变成关于 r 的多项式。
特征方程
对
a0y(n)+a1y(n−1)+⋯+any=0,a0=0,
试 y=erx,得到
P(r)=a0rn+a1rn−1+⋯+an=0.
每个特征根按照类型生成基解。
不同实根
根 r1,…,rn 不同时,
y=C1er1x+⋯+Cnernx.
实重根
若 r 是 m 重根,对应
erx,xerx,…,xm−1erx.
乘 xk 是为了补出线性无关的解,不能只把同一个指数项重复写 m 次。
复根
若 α±iβ 是一对复根,对应两个实解
eαxcosβx,eαxsinβx.
若这对复根重数为 m,还要依次乘 1,x,…,xm−1。
用算子因式分解理解
令 D=d/dx,方程写作
P(D)y=0.
若 P(r)=(r−r1)m1⋯(r−rk)mk,则
P(D)=(D−r1)m1⋯(D−rk)mk.
这个表示也解释了为什么重根产生 xkerx,并为非齐次方程中的共振规则做准备。
Euler–Cauchy 方程
典型形式
x2y′′+axy′+by=0,x=0.
可以试幂函数 y=xm,得到指标方程
m(m−1)+am+b=0.
在 x>0 上也可令
t=lnx,Y(t)=y(et).
利用
xy′=Y′,x2y′′=Y′′−Y′,
把 Euler 方程变成常系数方程。对 x<0 可用 t=ln∣x∣ 分区间处理,不能让解跨过奇点 x=0。
若指标根为重根 m,基解为
xm,xmln∣x∣.
复根 α±iβ 对应
∣x∣αcos(βln∣x∣),∣x∣αsin(βln∣x∣).
边值问题与特征值
例如
y′′+λy=0,y(0)=y(L)=0
并非任意 λ 都有非零解。根据 λ 的正、零、负三种情形分别写通解,边界条件会筛出离散的特征值
λn=(Lnπ)2,n=1,2,…
及对应特征函数 sin(nπx/L)。这说明边值条件不是简单“解两个常数”,有时会决定参数本身。
快速检查
- 特征多项式次数应等于方程阶数;
- 根的代数重数总和应为 n;
- 每个实根贡献其重数个解,每对复根贡献两倍重数个实解;
- Euler 方程的解区间不能跨 x=0;
- 写完通解后数一遍独立常数,并代回低阶示例检查符号。