卡在這題?用 AI 解開它
Calculus · real student question

解微分方程式 3(1 + t²) dy/dt = 2ty(y³ − 1)。

Question

解

3(1+t2)dydt=2ty(y3−1)3\left(1+t^2\right)\frac{dy}{dt}=2ty\left(y^3-1\right)

Step-by-step solution

  1. 先把平衡解找出來。 右邊在 y=0y=0 或 y3=1y^3=1 時為零,而常數的 yy 滿足 dydt=0\tfrac{dy}{dt}=0,所以

    y≡0與y≡1y\equiv 0\qquad\text{與}\qquad y\equiv 1

    都是原方程式的解。一開始就把它們找出來很重要,因為下面的代換會除以 yy。

  2. 認出這是白努利方程式的形式。 兩邊除以 3(1+t2)3(1+t^2) 並把右邊展開:

    dydt+2t3(1+t2) y=2t3(1+t2) y4\frac{dy}{dt}+\frac{2t}{3(1+t^2)}\,y=\frac{2t}{3(1+t^2)}\,y^{4}

    這是指數 n=4n=4 的白努利方程式,標準解法就是令 v=y1−n=y−3v=y^{1-n}=y^{-3}。

  3. 代入 v=y−3v=y^{-3}。 於是 dvdt=−3y−4dydt\dfrac{dv}{dt}=-3y^{-4}\dfrac{dy}{dt},再用 dydt=2ty(y3−1)3(1+t2)\dfrac{dy}{dt}=\dfrac{2ty(y^3-1)}{3(1+t^2)}:

    dvdt=−3y−4⋅2ty(y3−1)3(1+t2)=−2t(1−y−3)1+t2\frac{dv}{dt}=-3y^{-4}\cdot\frac{2ty(y^3-1)}{3(1+t^2)}=-\frac{2t\left(1-y^{-3}\right)}{1+t^2}

    所以

    dvdt−2t1+t2 v=−2t1+t2\frac{dv}{dt}-\frac{2t}{1+t^2}\,v=-\frac{2t}{1+t^2}

    注意正負號:vv 的係數是負的。在這裡把它寫成正的,是這道題最常見的失誤,而且會讓最後的答案整個翻掉。

  4. 求積分因子。 取 P(t)=−2t1+t2P(t)=-\dfrac{2t}{1+t^2},

    μ(t)=exp⁡ ⁣(∫P dt)=exp⁡(−ln⁡(1+t2))=11+t2\mu(t)=\exp\!\left(\int P\,dt\right)=\exp\left(-\ln\left(1+t^2\right)\right)=\frac{1}{1+t^2}

    兩邊同乘之後,左邊就併成單一個導數:

    ddt(v1+t2)=−2t(1+t2)2\frac{d}{dt}\left(\frac{v}{1+t^2}\right)=-\frac{2t}{\left(1+t^2\right)^2}

  5. 兩邊積分。 右邊本身又是一次代換(令 u=1+t2u=1+t^2):

    ∫−2t(1+t2)2 dt=11+t2+C\int-\frac{2t}{\left(1+t^2\right)^2}\,dt=\frac{1}{1+t^2}+C

    所以 v1+t2=11+t2+C\dfrac{v}{1+t^2}=\dfrac{1}{1+t^2}+C,也就是

    v=1+C(1+t2)v=1+C\left(1+t^2\right)

  6. 換回 yy 並驗算。 因為 v=y−3v=y^{-3},

    1y3=1+C(1+t2)⟹y=[1+C(1+t2)]−1/3\frac{1}{y^{3}}=1+C\left(1+t^{2}\right)\quad\Longrightarrow\quad y=\left[1+C\left(1+t^{2}\right)\right]^{-1/3}

    取 C=0.7C=0.7、t=1.7t=1.7 做數值檢查:左邊 3(1+t2)y′=−1.604485673(1+t^2)y'=-1.60448567,右邊 2ty(y3−1)=−1.604485672ty(y^3-1)=-1.60448567。✓\checkmark(取 C=0C=0 就還原成平衡解 y≡1y\equiv1。)

Answer

y=[1+C(1+t2)]−1/3,並且還有 y≡0y=\left[1+C\left(1+t^{2}\right)\right]^{-1/3},\qquad\text{並且還有 }y\equiv 0

Need to solve a different problem like this? Open the solver →