原文筆記

本課程筆記數學 · 4 篇

4-2 作业


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

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

  • 分析: 这是一个4阶齐次方程。虽然阶数高,但因为它没有一阶和三阶导数项,实际上解起来像解二元一次方程一样简单。
  • 解题步骤:
    1. 写出特征方程: 把 x(n)x^{(n)} 换成 λn\lambda^n: λ4−5λ2+4=0\lambda^4 - 5\lambda^2 + 4 = 0
    2. 求解特征根: 这看起来是个4次方程,但我们可以把它看作关于 λ2\lambda^2 的二次方程。 利用十字相乘法(或者求根公式):u2−5u+4=(u−1)(u−4)u^2 - 5u + 4 = (u-1)(u-4)。 所以: (λ2−1)(λ2−4)=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+c2e−t+c3e2t+c4e−2tx = c_1 \mathrm{e}^t + c_2 \mathrm{e}^{-t} + c_3 \mathrm{e}^{2t} + c_4 \mathrm{e}^{-2t}

(2) x′′′−3ax′′+3a2x′−a3x=0x''' - 3a x'' + 3a^2 x' - a^3 x = 0

  • 分析: 这里系数变成了字母 aa,把它当成普通数字处理即可。
  • 解题步骤:
    1. 写出特征方程:λ3−3aλ2+3a2λ−a3=0\lambda^3 - 3a\lambda^2 + 3a^2\lambda - a^3 = 0
    2. 求解特征根: 观察这个式子,系数是 1, -3, 3, -1。这正是二项式展开公式 (A−B)3=A3−3A2B+3AB2−B3(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. 写出通解: 对于三重根,我们需要补充乘上 tt 和 t2t^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) x′′′−4x′′+5x′−2x=2t+3x''' - 4x'' + 5x' - 2x = 2t + 3

  • 分析: 这是一个非齐次方程。
    • 第一步:先算左边等于0的通解(齐次解 xhx_h)。
    • 第二步:利用待定系数法算右边的特解(xpx_p)。
    • 第三步:加起来。
  • 步骤一:求齐次通解 xhx_h
    1. 特征方程:λ3−4λ2+5λ−2=0\lambda^3 - 4\lambda^2 + 5\lambda - 2 = 0。
    2. 试根法: 这种高次方程通常试一下 ±1,±2\pm 1, \pm 2。 试 λ=1\lambda = 1:1−4+5−2=01 - 4 + 5 - 2 = 0。成立!说明 (λ−1)(\lambda - 1) 是一个因子。
    3. 多项式除法(或综合除法): 用 (λ3−4λ2+5λ−2)(\lambda^3 - 4\lambda^2 + 5\lambda - 2) 除以 (λ−1)(\lambda - 1),得到 λ2−3λ+2\lambda^2 - 3\lambda + 2。
    4. 继续分解 λ2−3λ+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. 代入原方程:0−4(0)+5(A)−2(At+B)=2t+30 - 4(0) + 5(A) - 2(At + B) = 2t + 35A−2At−2B=2t+35A - 2At - 2B = 2t + 3 整理得到: −2At+(5A−2B)=2t+3-2A t + (5A - 2B) = 2t + 3
    6. 比较系数:
      • tt 的系数:−2A=2  ⟹  A=−1-2A = 2 \implies A = -1
      • 常数项:5A−2B=35A - 2B = 3。将 A=−1A=-1 代入:−5−2B=3  ⟹  −2B=8  ⟹  B=−4-5 - 2B = 3 \implies -2B = 8 \implies B = -4。 由此得到特解: x~=−t−4\widetilde{x} = -t - 4
  • 步骤三:合成通解x=(c1+c2t)et+c3e2t−t−4x = (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 是解,证明 lim⁡t→∞[φ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 的通解,并看它在 t→∞t \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+c2e−2t+c3e−3ty(t) = c_1 \mathrm{e}^{0t} + c_2 \mathrm{e}^{-2t} + c_3 \mathrm{e}^{-3t}y(t)=c1+c2e−2t+c3e−3ty(t) = c_1 + c_2 \mathrm{e}^{-2t} + c_3 \mathrm{e}^{-3t}
    6. 求极限: lim⁡t→∞y(t)=lim⁡t→∞(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}}) 当 t→∞t \to \infty 时,e2t\mathrm{e}^{2t} 和 e3t\mathrm{e}^{3t} 趋向于无穷大,倒数趋向于0。 lim⁡t→∞[φ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=ln⁡xt = \ln x,其中 tt 是新的自变量)。

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

  1. 一阶导数转换: 根据链式法则: dydx=dydt⋅dtdx\frac{\mathrm{d} y}{\mathrm{d} x} = \frac{\mathrm{d} y}{\mathrm{d} t} \cdot \frac{\mathrm{d} t}{\mathrm{d} x} 因为 t=ln⁡xt = \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)′=u′v+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+1xd2ydt2⋅1x= -\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(d2ydt2−dydt)= \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=d2ydt2−dydtx^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}

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


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

  • tx′t x' 替换为 dxds\frac{\mathrm{d}x}{\mathrm{d}s}
  • t2x′′t^2 x'' 替换为 d2xds2−dxds\frac{\mathrm{d}^2x}{\mathrm{d}s^2} - \frac{\mathrm{d}x}{\mathrm{d}s}

(1) t2x′′+tx′−x=0t^2 x'' + t x' - x = 0

  1. 代入变换: 将上面两个替换公式代入原方程: (d2xds2−dxds)+dxds−x=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. 化简:d2xds2−x=0\frac{\mathrm{d}^2x}{\mathrm{d}s^2} - x = 0 化简后变成了一个简单的常系数方程。
  3. 求解(关于 ss): 特征方程:λ2−1=0  ⟹  λ=±1\lambda^2 - 1 = 0 \implies \lambda = \pm 1。 解为:x=c1es+c2e−sx = c_1 \mathrm{e}^s + c_2 \mathrm{e}^{-s}。
  4. 回代变量(从 ss 变回 tt): 因为 t=est = \mathrm{e}^s,所以 es=t\mathrm{e}^s = t,e−s=1t\mathrm{e}^{-s} = \frac{1}{t}。 最终答案:x=c1t+c2tx = c_1 t + \frac{c_2}{t}

(2) t2x′′−4tx′+6x=tt^2 x'' - 4t x' + 6x = t

  1. 代入变换(令 t=est = \mathrm{e}^s): 左边: (d2xds2−dxds)−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. 整理方程:d2xds2−5dxds+6x=es\frac{\mathrm{d}^2x}{\mathrm{d}s^2} - 5\frac{\mathrm{d}x}{\mathrm{d}s} + 6x = \mathrm{e}^s 这是一个非齐次常系数方程。
  3. 解齐次部分: 特征方程:λ2−5λ+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。 代入方程: (A−5A+6A)es=es(A - 5A + 6A)\mathrm{e}^s = \mathrm{e}^s2Aes=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^2,e3s=t3\mathrm{e}^{3s} = t^3。 最终答案:x=c1t2+c2t3+12tx = c_1 t^2 + c_2 t^3 + \frac{1}{2}t