二阶线性微分方程幂级数解的递推关系验证及通解推导问询
关于(a)部分递推公式的正确性
你的推导思路完全正确!不过有个小细节需要修正:当 ( m=1 ) 时,代入递推公式应该是 ( (1+2)(1+1)a_{3} + \left(p + \frac{1}{2}\right)a_1 = 0 ),也就是 ( 6a_3 + \left(p + \frac{1}{2}\right)a_1 = 0 ),你之前写的 ( 3a_3 ) 是系数计算失误,但整体的推导逻辑没问题。
再帮你梳理一遍核心步骤,确认递推关系的合理性:
- 假设解的幂级数形式 ( y = \sum_{n=0}^{\infty} a_n x^n ),求导得到 ( y' = \sum_{n=1}^{\infty} n a_n x^{n-1} ),( y'' = \sum_{n=2}^{\infty} n(n-1)a_n x^{n-2} )
- 将这些代入原方程 ( y'' + \left(p + \frac{1}{2} - \frac{1}{4}x^2\right)y = 0 ),拆分后得到三个求和项:
[
\sum_{n=2}^{\infty} n(n-1)a_n x^{n-2} + \left(p + \frac{1}{2}\right)\sum_{n=0}^{\infty} a_n x^n - \frac{1}{4}\sum_{n=0}^{\infty} a_n x^{n+2} = 0
] - 统一所有求和项的幂次为 ( x^m ):
- 对第一个求和项,令 ( n = m+2 ),替换后变为 ( \sum_{m=0}^{\infty} (m+2)(m+1)a_{m+2}x^m )
- 对第三个求和项,令 ( n = m-2 ),替换后变为 ( \sum_{m=2}^{\infty} \frac{1}{4}a_{m-2}x^m )
- 合并同幂次的系数:
- 当 ( m=0 ):只有前两项有贡献,得到 ( 2a_2 + \left(p + \frac{1}{2}\right)a_0 = 0 )
- 当 ( m=1 ):只有前两项有贡献,得到 ( 6a_3 + \left(p + \frac{1}{2}\right)a_1 = 0 )
- 当 ( m \geq 2 ):三项都有贡献,得到 ( (m+1)(m+2)a_{m+2} + \left(p + \frac{1}{2}\right)a_m - \frac{1}{4}a_{m-2} = 0 )
(b)部分:幂级数形式的通解推导
原方程是二阶线性齐次微分方程,通解是两个线性无关解的线性组合。我们分两种情况,分别推导偶次幂主导的解 ( y_1(x) ) 和奇次幂主导的解 ( y_2(x) ):
情况1:偶次幂解 ( y_1(x) )(取 ( a_0 \neq 0, a_1 = 0 ))
此时所有奇数项系数 ( a_{2k+1}=0 ),只需要计算偶数项:
m=0:
[
2a_2 + \left(p + \frac{1}{2}\right)a_0 = 0 \implies a_2 = -\frac{1}{2}\left(p + \frac{1}{2}\right)a_0
]m=2:
代入递推公式 ( 12a_4 + \left(p + \frac{1}{2}\right)a_2 - \frac{1}{4}a_0 = 0 ),将 ( a_2 ) 代入后整理得:
[
a_4 = \frac{1}{24}\left(p + \frac{1}{2}\right)^2 a_0 + \frac{1}{48}a_0
]m=4:
代入递推公式 ( 30a_6 + \left(p + \frac{1}{2}\right)a_4 - \frac{1}{4}a_2 = 0 ),代入 ( a_4 ) 和 ( a_2 ) 并整理得:
[
a_6 = -\frac{1}{720}\left(p + \frac{1}{2}\right)^3 a_0 + \frac{1}{288}\left(p + \frac{1}{2}\right)a_0
]所以 ( y_1(x) ) 的前7项(到 ( x^6 ))为:
[
y_1(x) = a_0\left[1 - \frac{1}{2}\left(p + \frac{1}{2}\right)x^2 + \left(\frac{1}{24}\left(p + \frac{1}{2}\right)^2 + \frac{1}{48}\right)x^4 + \left(-\frac{1}{720}\left(p + \frac{1}{2}\right)^3 + \frac{1}{288}\left(p + \frac{1}{2}\right)\right)x^6 + \dots\right]
]
情况2:奇次幂解 ( y_2(x) )(取 ( a_1 \neq 0, a_0 = 0 ))
此时所有偶数项系数 ( a_{2k}=0 ),只需要计算奇数项:
m=1:
代入递推公式 ( 6a_3 + \left(p + \frac{1}{2}\right)a_1 = 0 ),得:
[
a_3 = -\frac{1}{6}\left(p + \frac{1}{2}\right)a_1
]m=3:
代入递推公式 ( 20a_5 + \left(p + \frac{1}{2}\right)a_3 - \frac{1}{4}a_1 = 0 ),代入 ( a_3 ) 后整理得:
[
a_5 = \frac{1}{120}\left(p + \frac{1}{2}\right)^2 a_1 + \frac{1}{80}a_1
]m=5:
代入递推公式 ( 42a_7 + \left(p + \frac{1}{2}\right)a_5 - \frac{1}{4}a_3 = 0 ),代入 ( a_5 ) 和 ( a_3 ) 并整理得:
[
a_7 = -\frac{1}{5040}\left(p + \frac{1}{2}\right)^3 a_1 + \frac{1}{1440}\left(p + \frac{1}{2}\right)a_1
]所以 ( y_2(x) ) 的前7项(到 ( x^7 ))为:
[
y_2(x) = a_1\left[x - \frac{1}{6}\left(p + \frac{1}{2}\right)x^3 + \left(\frac{1}{120}\left(p + \frac{1}{2}\right)^2 + \frac{1}{80}\right)x^5 + \left(-\frac{1}{5040}\left(p + \frac{1}{2}\right)^3 + \frac{1}{1440}\left(p + \frac{1}{2}\right)\right)x^7 + \dots\right]
]
最终通解
将两个线性无关解组合,得到原方程的通解:
[
y(x) = C_1 y_1(x) + C_2 y_2(x)
]
其中 ( C_1, C_2 ) 为任意常数。
备注:内容来源于stack exchange,提问作者esmeraldabrown

