用幂级数求解微分方程2y''+y'=0时的疑问
嘿,我来帮你捋清楚这个问题!咱们一步步拆解,看看你用幂级数求解时的误区在哪里,以及怎么和官方解对应上。
首先,先回忆一下幂级数解法的标准步骤:
假设微分方程的解可以表示为幂级数形式:y = Σ(n=0到∞) aₙxⁿ
然后求导得到:
y' = Σ(n=1到∞) n aₙxⁿ⁻¹y'' = Σ(n=2到∞) n(n-1)aₙxⁿ⁻²
把这两个代入原方程2y'' + y' = 0,为了合并级数,我们统一下标:
- 对
y'',令k = n-2,则n = k+2,转化为2Σ(k=0到∞) (k+2)(k+1)aₖ₊₂xᵏ - 对
y',令k = n-1,则n = k+1,转化为Σ(k=0到∞) (k+1)aₖ₊₁xᵏ
合并后每一项的系数必须为0,所以:2(k+2)(k+1)aₖ₊₂ + (k+1)aₖ₊₁ = 0(k≥0)
因为k+1≠0,两边除以k+1得到递推关系:aₖ₊₂ = -aₖ₊₁/[2(k+2)]
接下来我们用这个递推关系展开系数:
a₀和a₁是任意常数(二阶微分方程有两个线性无关解)a₂ = -a₁/(2×2),a₃ = -a₂/(2×3) = a₁/(2²×2×3),a₄ = -a₃/(2×4) = -a₁/(2³×2×3×4),以此类推
把这些系数代回y的展开式:
y = a₀ + a₁x + a₂x² + a₃x³ + ... = a₀ + a₁[ x - x²/(2×2!) + x³/(2²×3!) - x⁴/(2³×4!) + ... ]
现在重点来了:你觉得“多了x项”,是因为没把这个括号里的级数和已知的指数函数展开式对应上。我们知道e^t = Σ(n=0到∞) tⁿ/n!,那e^(-x/2)的展开式是:e^(-x/2) = 1 - x/(2×1!) + x²/(2²×2!) - x³/(2³×3!) + ...
对比括号里的级数,我们可以做个变形:
括号里的级数 = 2[ 1 - e^(-x/2) ]
(因为1 - e^(-x/2) = x/(2×1!) - x²/(2²×2!) + x³/(2³×3!) - ...,乘以2正好就是括号里的内容)
把这个代回y的表达式:y = a₀ + a₁×2[1 - e^(-x/2)] = (a₀ + 2a₁) + (-2a₁)e^(-x/2)
最后做常数替换:令c₁ = a₀ + 2a₁,c₂ = -2a₁,就得到了官方解的形式:y = c₁ + c₂e^(-x/2)
所以你所谓的“额外的x”其实是幂级数转化为指数函数时的中间形式,通过合并常数项就能和官方解完全等价——本质上只是常数的表达方式不同,你的解和官方解是同一个结果,只是换了个常数代号而已!
可能你出错的地方是:直接把自己得到的带x的级数当成了e^(-x/2)的展开,但实际上它是1 - e^(-x/2)的倍数,只要调整常数项就能统一成标准形式啦。
备注:内容来源于stack exchange,提问作者timetravel_0

