Skip to content

习题4.2


2. 求解下列常系数线性微分方程

(1) x(4)5x+4x=0x^{(4)} - 5x'' + 4x = 0

  • 分析: 这是一个4阶齐次方程。虽然阶数高,但因为它没有一阶和三阶导数项,实际上解起来像解二元一次方程一样简单。
  • 解题步骤:
    1. 写出特征方程:x(n)x^{(n)} 换成 λn\lambda^n

      λ45λ2+4=0\lambda^4 - 5\lambda^2 + 4 = 0

    2. 求解特征根: 这看起来是个4次方程,但我们可以把它看作关于 λ2\lambda^2 的二次方程。 利用十字相乘法(或者求根公式):u25u+4=(u1)(u4)u^2 - 5u + 4 = (u-1)(u-4)。 所以:

      (λ21)(λ24)=0(\lambda^2 - 1)(\lambda^2 - 4) = 0

      这意味着 λ2=1\lambda^2 = 1λ2=4\lambda^2 = 4。 解得四个互不相等的实根:

      λ1=1,λ2=1,λ3=2,λ4=2\lambda_1 = 1, \quad \lambda_2 = -1, \quad \lambda_3 = 2, \quad \lambda_4 = -2

    3. 写出通解: 因为是四个不同的实根,根据规则,直接线性叠加。

      x=c1et+c2et+c3e2t+c4e2tx = c_1 \mathrm{e}^t + c_2 \mathrm{e}^{-t} + c_3 \mathrm{e}^{2t} + c_4 \mathrm{e}^{-2t}


(2) x3ax+3a2xa3x=0x''' - 3a x'' + 3a^2 x' - a^3 x = 0

  • 分析: 这里系数变成了字母 aa,把它当成普通数字处理即可。
  • 解题步骤:
    1. 写出特征方程:

      λ33aλ2+3a2λa3=0\lambda^3 - 3a\lambda^2 + 3a^2\lambda - a^3 = 0

    2. 求解特征根: 观察这个式子,系数是 1, -3, 3, -1。这正是二项式展开公式 (AB)3=A33A2B+3AB2B3(A-B)^3 = A^3 - 3A^2B + 3AB^2 - B^3 的形式。 所以方程可以因式分解为:

      (λa)3=0(\lambda - a)^3 = 0

      这里得出特征值 λ=a\lambda = a注意: 这里的 aa 出现了3次,所以它是三重根(重数 k=3k=3)。
    3. 写出通解: 对于三重根,我们需要补充乘上 ttt2t^2 才能凑够3个独立的解。 基础解:eat\mathrm{e}^{at} 补乘一次:teatt\mathrm{e}^{at} 补乘两次:t2eatt^2\mathrm{e}^{at}通解:

      x=(c1+c2t+c3t2)eatx = (c_1 + c_2 t + c_3 t^2) \mathrm{e}^{at}


(6) x4x+5x2x=2t+3x''' - 4x'' + 5x' - 2x = 2t + 3

  • 分析: 这是一个非齐次方程。

    • 第一步:先算左边等于0的通解(齐次解 xhx_h)。
    • 第二步:利用待定系数法算右边的特解(xpx_p)。
    • 第三步:加起来。
  • 步骤一:求齐次通解 xhx_h

    1. 特征方程:λ34λ2+5λ2=0\lambda^3 - 4\lambda^2 + 5\lambda - 2 = 0
    2. 试根法: 这种高次方程通常试一下 ±1,±2\pm 1, \pm 2。 试 λ=1\lambda = 114+52=01 - 4 + 5 - 2 = 0。成立!说明 (λ1)(\lambda - 1) 是一个因子。
    3. 多项式除法(或综合除法):(λ34λ2+5λ2)(\lambda^3 - 4\lambda^2 + 5\lambda - 2) 除以 (λ1)(\lambda - 1),得到 λ23λ+2\lambda^2 - 3\lambda + 2
    4. 继续分解 λ23λ+2=(λ1)(λ2)\lambda^2 - 3\lambda + 2 = (\lambda - 1)(\lambda - 2)
    5. 综合起来,特征方程为:

      (λ1)2(λ2)=0(\lambda - 1)^2 (\lambda - 2) = 0

      特征根为:λ1=λ2=1\lambda_1 = \lambda_2 = 1 (二重根), λ3=2\lambda_3 = 2 (单根)。
    6. 齐次解:

      xh=(c1+c2t)et+c3e2tx_h = (c_1 + c_2 t)\mathrm{e}^t + c_3 \mathrm{e}^{2t}

  • 步骤二:求特解 xpx_p

    1. 观察右边: f(t)=2t+3f(t) = 2t + 3
      • 这是一个一次多项式 P1(t)P_1(t)
      • 对应指数形式是 e0t\mathrm{e}^{0t}(因为没有显式的 ee,指数就是0)。
    2. 检查共振(确定 kk): 右边的指数是 0。我们检查特征根(1和2),里面没有0。 所以,不是特征根,k=0k=0(不需要补乘 tt)。
    3. 设特解: 因为右边是一次多项式,我们设特解也是一次多项式:

      x~=At+B\widetilde{x} = At + B

    4. 代入求导:

      x~=A,x~=0,x~=0\widetilde{x}' = A, \quad \widetilde{x}'' = 0, \quad \widetilde{x}''' = 0

    5. 代入原方程:

      04(0)+5(A)2(At+B)=2t+30 - 4(0) + 5(A) - 2(At + B) = 2t + 3

      5A2At2B=2t+35A - 2At - 2B = 2t + 3

      整理得到:

      2At+(5A2B)=2t+3-2A t + (5A - 2B) = 2t + 3

    6. 比较系数:
      • tt 的系数:2A=2    A=1-2A = 2 \implies A = -1
      • 常数项:5A2B=35A - 2B = 3。将 A=1A=-1 代入:52B=3    2B=8    B=4-5 - 2B = 3 \implies -2B = 8 \implies B = -4。 由此得到特解:

      x~=t4\widetilde{x} = -t - 4

  • 步骤三:合成通解

    x=(c1+c2t)et+c3e2tt4x = (c_1 + c_2 t)\mathrm{e}^t + c_3 \mathrm{e}^{2t} - t - 4


4. 证明题

题目: 给定 x+5x+6x=f(t)x''' + 5x'' + 6x' = f(t)φ1,φ2\varphi_1, \varphi_2 是解,证明 limt[φ1(t)φ2(t)]\lim\limits_{t \to \infty} [\varphi_1(t) - \varphi_2(t)] 存在.

  • 思路引导: 这道题看起来很抽象,其实是在考查非齐次方程通解的结构。 核心定理: 非齐次方程的任意两个解之差,一定是对应齐次方程的解。

  • 证明步骤:

    1. y(t)=φ1(t)φ2(t)y(t) = \varphi_1(t) - \varphi_2(t)
    2. 因为 φ1,φ2\varphi_1, \varphi_2 都是原方程的解,代入方程相减:

      (φ1+5φ1+6φ1)(φ2+5φ2+6φ2)=f(t)f(t)=0(\varphi_1''' + 5\varphi_1'' + 6\varphi_1') - (\varphi_2''' + 5\varphi_2'' + 6\varphi_2') = f(t) - f(t) = 0

      整理得:

      y+5y+6y=0y''' + 5y'' + 6y' = 0

      这说明 差函数 y(t)y(t) 满足对应的齐次方程
    3. 现在问题转化为:求齐次方程 y+5y+6y=0y''' + 5y'' + 6y' = 0 的通解,并看它在 tt \to \infty 时的极限。
    4. 特征方程:

      λ3+5λ2+6λ=0\lambda^3 + 5\lambda^2 + 6\lambda = 0

      提取公因式 λ\lambda

      λ(λ2+5λ+6)=0\lambda(\lambda^2 + 5\lambda + 6) = 0

      λ(λ+2)(λ+3)=0\lambda(\lambda + 2)(\lambda + 3) = 0

      特征根为:λ1=0,λ2=2,λ3=3\lambda_1 = 0, \quad \lambda_2 = -2, \quad \lambda_3 = -3
    5. 齐次方程的解(即 φ1φ2\varphi_1 - \varphi_2 的表达式):

      y(t)=c1e0t+c2e2t+c3e3ty(t) = c_1 \mathrm{e}^{0t} + c_2 \mathrm{e}^{-2t} + c_3 \mathrm{e}^{-3t}

      y(t)=c1+c2e2t+c3e3ty(t) = c_1 + c_2 \mathrm{e}^{-2t} + c_3 \mathrm{e}^{-3t}

    6. 求极限:

      limty(t)=limt(c1+c2e2t+c3e3t)\lim_{t \to \infty} y(t) = \lim_{t \to \infty} (c_1 + \frac{c_2}{\mathrm{e}^{2t}} + \frac{c_3}{\mathrm{e}^{3t}})

      tt \to \infty 时,e2t\mathrm{e}^{2t}e3t\mathrm{e}^{3t} 趋向于无穷大,倒数趋向于0。

      limt[φ1(t)φ2(t)]=c1+0+0=c1\lim_{t \to \infty} [\varphi_1(t) - \varphi_2(t)] = c_1 + 0 + 0 = c_1

      这是一个常数,说明极限存在。证毕。

5. 欧拉方程 (Euler Equation)

这部分是全新的知识点,我们先搞懂它的原理。 核心难点: 变量代换中的求导法则(链式法则)。

定义与证明 题目给出的定义是:xndnydxn+=0x^n \frac{\mathrm{d}^n y}{\mathrm{d} x^n} + \cdots = 0。这里的自变量是 xx,因变量是 yy。 变换公式:x=etx = \mathrm{e}^t (这意味着 t=lnxt = \ln x,其中 tt 是新的自变量)。

证明变换过程(从零推导): 我们要把关于 xx 的导数 dydx\frac{\mathrm{d}y}{\mathrm{d}x} 变成关于 tt 的导数 dydt\frac{\mathrm{d}y}{\mathrm{d}t}

  1. 一阶导数转换: 根据链式法则:

    dydx=dydtdtdx\frac{\mathrm{d} y}{\mathrm{d} x} = \frac{\mathrm{d} y}{\mathrm{d} t} \cdot \frac{\mathrm{d} t}{\mathrm{d} x}

    因为 t=lnxt = \ln x,所以 dtdx=1x\frac{\mathrm{d} t}{\mathrm{d} x} = \frac{1}{x}

    dydx=1xdydt\frac{\mathrm{d} y}{\mathrm{d} x} = \frac{1}{x} \frac{\mathrm{d} y}{\mathrm{d} t}

    两边同乘 xx,得出一个重要结论:

    xdydx=dydt(记作 Dy)x \frac{\mathrm{d} y}{\mathrm{d} x} = \frac{\mathrm{d} y}{\mathrm{d} t} \quad \text{(记作 D}y)

  2. 二阶导数转换: 我们需要算 d2ydx2=ddx(dydx)\frac{\mathrm{d}^2 y}{\mathrm{d} x^2} = \frac{\mathrm{d}}{\mathrm{d} x} \left( \frac{\mathrm{d} y}{\mathrm{d} x} \right)。 把刚才的一阶结果代进去:

    ddx(1xdydt)\frac{\mathrm{d}}{\mathrm{d} x} \left( \frac{1}{x} \frac{\mathrm{d} y}{\mathrm{d} t} \right)

    利用乘法求导法则 (uv)=uv+uv(uv)' = u'v + uv'

    =(1x)dydt+1xddx(dydt)= \left( \frac{1}{x} \right)' \frac{\mathrm{d} y}{\mathrm{d} t} + \frac{1}{x} \frac{\mathrm{d}}{\mathrm{d} x} \left( \frac{\mathrm{d} y}{\mathrm{d} t} \right)

    =1x2dydt+1x[ddt(dydt)dtdx](再次用链式法则)= -\frac{1}{x^2} \frac{\mathrm{d} y}{\mathrm{d} t} + \frac{1}{x} \left[ \frac{\mathrm{d}}{\mathrm{d} t} \left( \frac{\mathrm{d} y}{\mathrm{d} t} \right) \cdot \frac{\mathrm{d} t}{\mathrm{d} x} \right] \quad \text{(再次用链式法则)}

    =1x2dydt+1xd2ydt21x= -\frac{1}{x^2} \frac{\mathrm{d} y}{\mathrm{d} t} + \frac{1}{x} \frac{\mathrm{d}^2 y}{\mathrm{d} t^2} \cdot \frac{1}{x}

    =1x2(d2ydt2dydt)= \frac{1}{x^2} \left( \frac{\mathrm{d}^2 y}{\mathrm{d} t^2} - \frac{\mathrm{d} y}{\mathrm{d} t} \right)

    两边同乘 x2x^2,得出第二个重要结论:

    x2d2ydx2=d2ydt2dydtx^2 \frac{\mathrm{d}^2 y}{\mathrm{d} x^2} = \frac{\mathrm{d}^2 y}{\mathrm{d} t^2} - \frac{\mathrm{d} y}{\mathrm{d} t}

结论: 这个变换把前面的变系数 xxx2x^2 完美消掉了,剩下的全是常数系数!这就变成了我们学过的常系数方程。


特别注意: 在下面的两道小题中,题目把字母换了! 题目写的是:t2x+txx=0t^2 x'' + t x' - x = 0。 这里:自变量是 tt,因变量是 xx。 为了不混淆,我们令变换为 t=est = \mathrm{e}^s (也就是引进新变量 s=lnts = \ln t)。 根据刚才的推导,我们有直接替换公式:

  • txt x' 替换为 dxds\frac{\mathrm{d}x}{\mathrm{d}s}
  • t2xt^2 x'' 替换为 d2xds2dxds\frac{\mathrm{d}^2x}{\mathrm{d}s^2} - \frac{\mathrm{d}x}{\mathrm{d}s}

(1) t2x+txx=0t^2 x'' + t x' - x = 0

  1. 代入变换: 将上面两个替换公式代入原方程:

    (d2xds2dxds)+dxdsx=0\left( \frac{\mathrm{d}^2x}{\mathrm{d}s^2} - \frac{\mathrm{d}x}{\mathrm{d}s} \right) + \frac{\mathrm{d}x}{\mathrm{d}s} - x = 0

  2. 化简:

    d2xds2x=0\frac{\mathrm{d}^2x}{\mathrm{d}s^2} - x = 0

    化简后变成了一个简单的常系数方程。
  3. 求解(关于 ss): 特征方程:λ21=0    λ=±1\lambda^2 - 1 = 0 \implies \lambda = \pm 1。 解为:x=c1es+c2esx = c_1 \mathrm{e}^s + c_2 \mathrm{e}^{-s}
  4. 回代变量(从 ss 变回 tt): 因为 t=est = \mathrm{e}^s,所以 es=t\mathrm{e}^s = tes=1t\mathrm{e}^{-s} = \frac{1}{t}最终答案:

    x=c1t+c2tx = c_1 t + \frac{c_2}{t}


(2) t2x4tx+6x=tt^2 x'' - 4t x' + 6x = t

  1. 代入变换(令 t=est = \mathrm{e}^s): 左边:

    (d2xds2dxds)4dxds+6x\left( \frac{\mathrm{d}^2x}{\mathrm{d}s^2} - \frac{\mathrm{d}x}{\mathrm{d}s} \right) - 4\frac{\mathrm{d}x}{\mathrm{d}s} + 6x

    右边:tt 变成了 es\mathrm{e}^s
  2. 整理方程:

    d2xds25dxds+6x=es\frac{\mathrm{d}^2x}{\mathrm{d}s^2} - 5\frac{\mathrm{d}x}{\mathrm{d}s} + 6x = \mathrm{e}^s

    这是一个非齐次常系数方程。
  3. 解齐次部分: 特征方程:λ25λ+6=0    (λ2)(λ3)=0\lambda^2 - 5\lambda + 6 = 0 \implies (\lambda-2)(\lambda-3)=0。 齐次解:xh=c1e2s+c2e3sx_h = c_1 \mathrm{e}^{2s} + c_2 \mathrm{e}^{3s}
  4. 解特解部分: 右边是 es\mathrm{e}^s(指数是1)。 检查:1 不是特征根(特征根是2和3)。 设特解 x~=Aes\widetilde{x} = A \mathrm{e}^s。 代入方程:

    (A5A+6A)es=es(A - 5A + 6A)\mathrm{e}^s = \mathrm{e}^s

    2Aes=es    A=122A \mathrm{e}^s = \mathrm{e}^s \implies A = \frac{1}{2}

    所以特解是 12es\frac{1}{2}\mathrm{e}^s
  5. 写出关于 ss 的通解:

    x=c1e2s+c2e3s+12esx = c_1 \mathrm{e}^{2s} + c_2 \mathrm{e}^{3s} + \frac{1}{2}\mathrm{e}^s

  6. 回代变量(变回 tt): 利用 es=t\mathrm{e}^s = t,则 e2s=t2\mathrm{e}^{2s} = t^2e3s=t3\mathrm{e}^{3s} = t^3最终答案:

    x=c1t2+c2t3+12tx = c_1 t^2 + c_2 t^3 + \frac{1}{2}t


書體

本站所載,間有由 AI 所生成者。其辭義真偽,請君自審之。