7. 级数法求解微分方程 Power Series Solutions

  • 二阶变系数齐次 ODE 的通式 Px+Qx+Rx=0Px''+Qx'+Rx = 0

    • 首先通过令P0P \not = 0 同时除以 PP 得到 x+px+qx=0x'' + px' +qx = 0

    • 级数收敛定理

      • 对于函数 p,q 通过级数展开,拥有收敛半径 ρ1,ρ2\rho _1, \rho_2,则解集 xx 也拥有大于等于 min(ρ1,ρ2)\min(\rho_1,\rho_2) 的级数展开

      • 也就是说,取两个收敛的min以内的部分的 t 一定满足收敛,即可以展开

      • 或者说,通过级数收敛得到的答案大概在 (ρmin,+ρmin)(-\rho_{min}, +\rho_{min}) 区间上是可以展开的,外部就会 \to \infty

        • 这里说的是至少,也就是说并不是在 ρmin\rho_{min} 处就会立即发散
  • 欧拉公式 Euler’s Equation

    • 定义:t2x+αtx+βx=0t^2x''+\alpha tx'+\beta x = 0

    • 奇点位置 0,我们主要看 t>0t>0 的部分,令 x=trx =t^r

    • 得到等式 r=α12±12(α1)24βr = -\frac{\alpha -1 }{2}\pm \frac 12\sqrt{(\alpha - 1)^2-4\beta} , 一个很眼熟的二次方程的解

      • Δ>0\Delta > 0 我们有 x=c1tr1+c2tr2x = c_1 t^{r_1} +c_2 t^{r_2}
      • Δ=0\Delta = 0 我们通过 reduction of ratio 得到一个 x2=tr1lntx_2 = t^{r_1} \ln t
  • 奇点 Singular Points

    • x+p(t)x+q(t)=0x'' + p(t) x' +q(t) = 0 , 其中 p,qp,q 是复数平面的函数

    • 定义:如果 p,q 在 t0t_0 附近可分析 (有限),则称这个点为 ordinary point, 反之为 singular point

    • regular:对于奇点,令 p 在 t0t_0 有一个pole 且对 t0t_0重数不超过 1, q 在t0t_0 有一个pole 且重数不超过 2, 则称这个点为 regular;反之为 irregular

    • Frobenius ansatz: 对于一般的 $x’’ + px’ +qx = 0 $ 我们变换为欧拉形式 t2x+t(tp(t))x+t2q(t)x=0t^2 x'' + t(tp(t))\cdot x' +t^2 q(t)\cdot x = 0

      • 这样我们得到 试根为 x(t)=trΣaktkx(t) = t^r\Sigma a_kt ^k

      • 首先将试根带回原方程得到 F(x):=x(x1)+p0x+q0F(x):= x(x-1)+p_0x + q_0, 其中 p0p_0 表示的是 p(t)p(t) 的泰勒展开的第一项,也就是 p(t0)p(t_0) , q 同理

        • 这其实是一个等式 F(x)=0F(x) = 0
        • 这是一个二次方程,但是在这里我们只考虑一个根
      • 再通过爆算得到等式 $ ((r+ m )(r+m-1)+q_0+(r+m)p_0)a_m+\Sigma_{k=0}^{m-1}(q_{m-k}+(r+k)p_{m-k})a_k$

      • 那么我们对 $F(x) $ 在 rr 附近展开,得到 indicial equation F(r)=0,amF(r+m)=Σk=0m1(qmk+(r+k)pmk)akF(r) = 0, a_mF(r+m) = -\Sigma_{k=0}^{m-1}(q_{m-k}+(r+k)p_{m-k})a_k

    • 第一类贝塞尔方程 JvJ_v Bessel Function of First Kind

      • t2x+tx+(t2v2)x=0t^2x'' + tx'+(t^2-v^2)x = 0
      • t=0t = 0 处得到 x(t)=0x(t) = 0
    • 贝塞尔方程的根

      • tt\to \infty, 则 有 w=txw = \sqrt t x, 那么我们有等式 w+w+1x2(14v2)w=0w'' + w' + \frac 1{x^2}(\frac 1 4 - v^2)w = 0 成立, 其中 w0w\to 0tt\to \infty 时候成立,那么这个结果就变成了一个类似一阶振动的微分方程,也就是说解 xccos(ωt)x \approx c \cdot \cos (\omega t)
      • 那么也就是说第一类贝塞尔方程具有无数的解,且 在 tt\to \infty 的时候根的间距大概为 π\pi