114 年 國立臺灣大學數學系碩士班《常微分方程式》

📄 試題原卷 免費註冊後即可對照原始考卷 PDF免費註冊

第 1 題25 分

Solve the following equations.
(a) dydx=1+(3−x+y)3\frac{dy}{dx} = 1 + (3-x+y)^3.
(b) ey2dydx+(1+y2)=0e^{y^2} \frac{dy}{dx} + (1+y^2) = 0.
(c) d2ydx2+xdydx+y=0\frac{d^2y}{dx^2} + x\frac{dy}{dx} + y = 0, y(0)=3y(0) = 3, y′(0)=0y'(0) = 0.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題分別考查:

  • (a) 變數代換與可分離變數法。
  • (b) 可分離變數法,以及非初等積分的表示方式。
  • (c) 尋找微分方程中的全導數,將二階方程降階,再利用初始條件求解。

本題不是選擇題,因此沒有需逐一判斷的選項。


(a) dydx=1+(3−x+y)3\displaystyle \frac{dy}{dx}=1+(3-x+y)^3

解題方法

右側的複合式 3−x+y3-x+y 直接出現在三次方中,令

u=3−x+y.u=3-x+y.

對 xx 微分:

dudx=−1+dydx.\frac{du}{dx}=-1+\frac{dy}{dx}.

`

由原方程

dydx=1+u3,\frac{dy}{dx}=1+u^3,

因此

dudx=u3.\frac{du}{dx}=u^3.

此方程為可分離變數方程。

計算

當 u≠0u\neq 0 時,

u−3 du=dx.u^{-3}\,du=dx.

積分得

∫u−3 du=∫dx,\int u^{-3}\,du=\int dx, −12u2=x+C.-\frac{1}{2u^2}=x+C.

重新整理常數後可寫成

1u2=C−2x.\frac{1}{u^2}=C-2x.

因此

u=±1C−2x.u=\pm \frac{1}{\sqrt{C-2x}}.

代回 u=3−x+yu=3-x+y:

3−x+y=±1C−2x,3-x+y=\pm \frac{1}{\sqrt{C-2x}},

所以

y=x−3±1C−2x.\boxed{y=x-3\pm \frac{1}{\sqrt{C-2x}}}.

另外,分離變數時除以了 u3u^3,必須補回 u=0u=0 的解。由

3−x+y=0,3-x+y=0,

得到

y=x−3.\boxed{y=x-3}.

解題技巧與檢查

看到 3−x+y3-x+y 這種「一次式整體出現」時,優先令其為新變數。代換後:

u′=y′−1,u'=y'-1,

恰好消去原方程右側的常數 11,形成最簡單的可分離方程 u′=u3u'=u^3。

此外,u=0u=0 是平衡解,不能因分離變數時除以 u3u^3 而遺漏。


(b) ey2dydx+(1+y2)=0\displaystyle e^{y^2}\frac{dy}{dx}+(1+y^2)=0

核心觀念

先將方程整理為

dydx=−(1+y2)e−y2.\frac{dy}{dx}=-(1+y^2)e^{-y^2}.

這是可分離變數方程。將 yy 的部分與 xx 的部分分開即可。

計算

由

ey2dydx=−(1+y2),e^{y^2}\frac{dy}{dx}=-(1+y^2),

可得

ey21+y2 dy=−dx.\frac{e^{y^2}}{1+y^2}\,dy=-dx.

兩邊積分:

∫ey21+y2 dy=−x+C.\int \frac{e^{y^2}}{1+y^2}\,dy=-x+C.

因此通解可表示為

∫ey21+y2 dy+x=C.\boxed{\int \frac{e^{y^2}}{1+y^2}\,dy+x=C}.

若希望將積分寫成定積分形式,可定義

F(y)=∫0yet21+t2 dt.F(y)=\int_0^y \frac{e^{t^2}}{1+t^2}\,dt.

則解為

x+F(y)=C,\boxed{x+F(y)=C},

即

x+∫0yet21+t2 dt=C.\boxed{x+\int_0^y\frac{e^{t^2}}{1+t^2}\,dt=C}.
🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 2 題30 分

Let A=(11−45)A = \begin{pmatrix} 1 & 1 \\ -4 & 5 \end{pmatrix}.
(a) Find eAte^{At}.
(b) Solve the system x′(t)=Ax(t)x'(t) = Ax(t), x(0)=(10)x(0) = \begin{pmatrix} 1 \\ 0 \end{pmatrix}.
(c) Solve the system x′(t)=Ax(t)+(−21)e−2tx'(t) = Ax(t) + \begin{pmatrix} -2 \\ 1 \end{pmatrix} e^{-2t}, x(0)=(10)x(0) = \begin{pmatrix} 1 \\ 0 \end{pmatrix}.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查線性常係數系統

x′(t)=Ax(t)+f(t)x'(t)=Ax(t)+f(t)

的矩陣指數 eAte^{At}、齊次系統解,以及變數變換法(variation of parameters)。

關鍵觀察是將矩陣 AA 分解為

A=3I+B,B=A−3I=(−21−42).A=3I+B, \qquad B=A-3I= \begin{pmatrix} -2&1\\ -4&2 \end{pmatrix}.

直接計算可得

B2=(−21−42)(−21−42)=(0000).B^2= \begin{pmatrix} -2&1\\ -4&2 \end{pmatrix} \begin{pmatrix} -2&1\\ -4&2 \end{pmatrix} = \begin{pmatrix} 0&0\\ 0&0 \end{pmatrix}.

因此 BB 是冪零矩陣,矩陣指數級數只保留一次項:

eBt=I+tB.e^{Bt}=I+tB.

又因為 3I3I 與 BB 可交換,

eAt=e(3I+B)t=e3teBt=e3t(I+tB).e^{At}=e^{(3I+B)t}=e^{3t}e^{Bt} =e^{3t}(I+tB).

解題方法

(a)求 eAte^{At}

由 B2=0B^2=0,

eBt=I+Bt+(Bt)22!+⋯=I+tB.e^{Bt} =I+Bt+\frac{(Bt)^2}{2!}+\cdots =I+tB.

因此

eAt=e3t[I+t(−21−42)].e^{At} =e^{3t} \left[ I+t \begin{pmatrix} -2&1\\ -4&2 \end{pmatrix} \right].

整理矩陣得

eAt=e3t(1−2tt−4t1+2t).\boxed{ e^{At} = e^{3t} \begin{pmatrix} 1-2t&t\\ -4t&1+2t \end{pmatrix} }.

(b)求解齊次系統

題目給定

x′(t)=Ax(t),x(0)=(10).x'(t)=Ax(t), \qquad x(0)= \begin{pmatrix} 1\\ 0 \end{pmatrix}.

常係數齊次系統的解為

x(t)=eAtx(0).x(t)=e^{At}x(0).

代入(a)的結果:

x(t)=e3t(1−2tt−4t1+2t)(10).x(t) = e^{3t} \begin{pmatrix} 1-2t&t\\ -4t&1+2t \end{pmatrix} \begin{pmatrix} 1\\ 0 \end{pmatrix}.

因此

x(t)=e3t(1−2t−4t).\boxed{ x(t)= e^{3t} \begin{pmatrix} 1-2t\\ -4t \end{pmatrix} }.

(c)求解非齊次系統

現在考慮

x′(t)=Ax(t)+(−21)e−2t,x(0)=(10).x'(t)=Ax(t)+ \begin{pmatrix} -2\\ 1 \end{pmatrix}e^{-2t}, \qquad x(0)= \begin{pmatrix} 1\\ 0 \end{pmatrix}.

令

b=(−21).b= \begin{pmatrix} -2\\ 1 \end{pmatrix}.

由常數變易法,解為

x(t)=eAtx(0)+∫0teA(t−s)b e−2s ds.x(t)=e^{At}x(0)+ \int_0^t e^{A(t-s)}b\,e^{-2s}\,ds.

由

eA(t−s)=e3(t−s)[I+(t−s)B],e^{A(t-s)} = e^{3(t-s)} \left[I+(t-s)B\right],

且

Bb=(−21−42)(−21)=(510),Bb= \begin{pmatrix} -2&1\\ -4&2 \end{pmatrix} \begin{pmatrix} -2\\ 1 \end{pmatrix} = \begin{pmatrix} 5\\ 10 \end{pmatrix},

故非齊次積分項為

∫0teA(t−s)b e−2s ds=e3t∫0te−5s[b+(t−s)Bb]ds.\begin{aligned} \int_0^t e^{A(t-s)}b\,e^{-2s}\,ds &= e^{3t}\int_0^t e^{-5s} \left[b+(t-s)Bb\right]ds. \end{aligned}

所需的兩個積分為

∫0te−5s ds=1−e−5t5,\int_0^t e^{-5s}\,ds = \frac{1-e^{-5t}}{5},

以及

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 3 題25 分

Assume that both y1(x)y_1(x) and y2(x)y_2(x) have continuous second derivatives on [0,π][0, \pi] and satisfies
(E1)

{y′′+y=0for x∈[0,π/2]y′′+xy=0for x∈[π/2,π]\begin{cases} y'' + y = 0 & \text{for } x \in [0, \pi/2] \\ y'' + xy = 0 & \text{for } x \in [\pi/2, \pi] \end{cases}

(a) Find the solution y1(x)y_1(x) with y1(π)=2y_1(\pi) = 2.
(b) Is it possible that y2(0)=1y_2(0) = 1?
(c) Let S={y:[0,π]→R∣y has a continuous second derivative and satisfies (E1)}S = \{y: [0, \pi] \rightarrow \mathbb{R} \mid y \text{ has a continuous second derivative and satisfies (E1)}\}. Show that SS is a vector space over R\mathbb{R} and find its dimension.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  1. 二階線性齊次微分方程的通解。
  2. 分段微分方程在分界點的接續條件。
  3. 連續二階導數所帶來的額外限制。
  4. 向量空間的封閉性與維度。
  5. 方程式
y′′+xy=0y''+xy=0

可用 Airy 函數表示其解。

令

a=π2.a=\frac{\pi}{2}.

一、分界點的必要條件

在左半段,

y′′+y=0,y''+y=0,

因此由左側逼近 aa 可得

y′′(a−)=−y(a).y''(a^-)= -y(a).

在右半段,

y′′+xy=0,y''+xy=0,

因此由右側逼近 aa 可得

y′′(a+)=−a y(a).y''(a^+)= -a\,y(a).

題目要求 yy 具有連續的二階導數,所以

y′′(a−)=y′′(a+).y''(a^-)=y''(a^+).

因此

−y(a)=−a y(a),-y(a)=-a\,y(a),

即

(a−1)y(a)=0.(a-1)y(a)=0.

由於 a=π2≠1a=\dfrac{\pi}{2}\ne 1,故必須有

y(π2)=0.\boxed{y\left(\frac{\pi}{2}\right)=0}.

這是本題最關鍵的條件。若只要求 y,y′y,y' 連續,不能得到此限制;正是因為題目要求二階導數連續,才會產生 y(a)=0y(a)=0。


二、右半段的基本解

右半段方程為

y′′+xy=0.y''+xy=0.

Airy 函數 Ai⁡\operatorname{Ai} 與 Bi⁡\operatorname{Bi} 滿足

Ai⁡′′(z)−zAi⁡(z)=0,Bi⁡′′(z)−zBi⁡(z)=0.\operatorname{Ai}''(z)-z\operatorname{Ai}(z)=0, \qquad \operatorname{Bi}''(z)-z\operatorname{Bi}(z)=0.

因此 Ai⁡(−x)\operatorname{Ai}(-x) 與 Bi⁡(−x)\operatorname{Bi}(-x) 滿足

y′′+xy=0.y''+xy=0.

定義函數 Φ\Phi 為右半段的基本解:

Φ(x)=π[Bi⁡(−a)Ai⁡(−x)−Ai⁡(−a)Bi⁡(−x)].\Phi(x) = \pi\left[ \operatorname{Bi}(-a)\operatorname{Ai}(-x) - \operatorname{Ai}(-a)\operatorname{Bi}(-x) \right].

由定義可得

Φ(a)=0.\Phi(a)=0.

利用 Airy 函數的 Wronskian 關係

Ai⁡(z)Bi⁡′(z)−Ai⁡′(z)Bi⁡(z)=1π,\operatorname{Ai}(z)\operatorname{Bi}'(z) - \operatorname{Ai}'(z)\operatorname{Bi}(z) = \frac{1}{\pi},

可得

Φ′(a)=1.\Phi'(a)=1.

所以 Φ\Phi 是滿足

Φ(a)=0,Φ′(a)=1\Phi(a)=0,\qquad \Phi'(a)=1

的唯一解。


(a) 求滿足 y1(π)=2y_1(\pi)=2 的解

左半段

左半段方程

y′′+y=0y''+y=0

的通解為

y(x)=Acos⁡x+Bsin⁡x.y(x)=A\cos x+B\sin x.

由 y(a)=0y(a)=0 得

Acos⁡π2+Bsin⁡π2=B=0.A\cos\frac{\pi}{2}+B\sin\frac{\pi}{2}=B=0.

因此左半段可寫成

y(x)=Acos⁡x.y(x)=A\cos x.

為了使用在 x=ax=a 的導數條件,令

c=y′(a).c=y'(a).

由於 y(a)=0y(a)=0,左半段解可直接寫成

y(x)=csin⁡(x−a).y(x)=c\sin(x-a).

確實有

y(a)=0,y′(a)=c.y(a)=0,\qquad y'(a)=c.

右半段

右半段解必須滿足

y(a)=0,y′(a)=c.y(a)=0,\qquad y'(a)=c.

由 Φ(a)=0\Phi(a)=0、Φ′(a)=1\Phi'(a)=1,可知右半段為

y(x)=cΦ(x).y(x)=c\Phi(x).

再利用邊界條件 y(π)=2y(\pi)=2:

cΦ(π)=2,c\Phi(\pi)=2,

因此

c=2Φ(π).c=\frac{2}{\Phi(\pi)}.

所以

y1(x)={2Φ(π)sin⁡(x−π2),0≤x≤π2,2Φ(x)Φ(π),π2≤x≤π,\boxed{ y_1(x)= \begin{cases} \dfrac{2}{\Phi(\pi)} \sin\left(x-\dfrac{\pi}{2}\right), & 0\le x\le \dfrac{\pi}{2},\\[1.2em] \dfrac{2\Phi(x)}{\Phi(\pi)}, & \dfrac{\pi}{2}\le x\le \pi, \end{cases} }

其中

Φ(x)=π[Bi⁡(−π2)Ai⁡(−x)−Ai⁡(−π2)Bi⁡(−x)].\boxed{ \Phi(x) = \pi\left[ \operatorname{Bi}\left(-\frac{\pi}{2}\right)\operatorname{Ai}(-x) - \operatorname{Ai}\left(-\frac{\pi}{2}\right)\operatorname{Bi}(-x) \right]. }

為何 Φ(π)≠0\Phi(\pi)\ne 0?

由 Φ(a)=0\Phi(a)=0、Φ′(a)=1\Phi'(a)=1,且右半段方程為

Φ′′+xΦ=0.\Phi''+x\Phi=0.

在 [a,π][a,\pi] 上有 x≤πx\le \pi。與方程

z′′+πz=0z''+\pi z=0

比較,其從 aa 出發的第一個正零點為

a+ππ=π2+π>π.a+\frac{\pi}{\sqrt{\pi}} = \frac{\pi}{2}+\sqrt{\pi} > \pi.

因此 Φ\Phi 在 (a,π](a,\pi] 上不會再次變成零,故

Φ(π)≠0.\Phi(\pi)\ne 0.

(b) 是否可能有 y2(0)=1y_2(0)=1?

由分界點條件,所有解都必須滿足

y2(π2)=0.y_2\left(\frac{\pi}{2}\right)=0.

在左半段,

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 4 題20 分

Let yy satisfy
y′′+y3=0.y'' + y^3 = 0.
(a) Show that (y′2)2+y44=constant\left(\frac{y'}{2}\right)^2 + \frac{y^4}{4} = \text{constant}.
(b) Assume y(0)≠0y(0) \neq 0. Show that yy is a periodic function.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  • 將微分方程乘以 y′y',建立能量守恆式。
  • 利用位能 V(y)=y44V(y)=\dfrac{y^4}{4} 的性質,判斷解的運動範圍。
  • 由相平面上的封閉軌道或有限週期積分,證明解為週期函數。

不過,題目第(a)小題的係數與微分方程不一致。


(a)題目式子存在係數錯誤

題目聲稱

(y′2)2+y44=(y′)24+y44\left(\frac{y'}{2}\right)^2+\frac{y^4}{4} =\frac{(y')^2}{4}+\frac{y^4}{4}

為常數。對其微分:

ddx((y′)24+y44)=y′y′′2+y3y′.\frac{d}{dx}\left(\frac{(y')^2}{4}+\frac{y^4}{4}\right) = \frac{y'y''}{2}+y^3y'.

整理得

ddx((y′)24+y44)=y′(y′′2+y3).\frac{d}{dx}\left(\frac{(y')^2}{4}+\frac{y^4}{4}\right) = y'\left(\frac{y''}{2}+y^3\right).

由原方程

y′′+y3=0y''+y^3=0

可得 y′′=−y3y''=-y^3,因此

ddx((y′)24+y44)=12y3y′,\frac{d}{dx}\left(\frac{(y')^2}{4}+\frac{y^4}{4}\right) = \frac{1}{2}y^3y',

一般並不等於 00。所以題目所列的式子並非守恆量。

例如取某一點使得 y=1y=1、y′=1y'=1,則上式導數為

12≠0.\frac{1}{2}\neq 0.

因此第(a)小題照原題文字無法成立。

正確的能量守恆式應為

(y′)22+y44=C.\boxed{\frac{(y')^2}{2}+\frac{y^4}{4}=C}.

證明如下。將方程式乘以 y′y':

y′y′′+y3y′=0.y'y''+y^3y'=0.

注意到

y′y′′=ddx((y′)22),y3y′=ddx(y44),y'y''=\frac{d}{dx}\left(\frac{(y')^2}{2}\right), \qquad y^3y'=\frac{d}{dx}\left(\frac{y^4}{4}\right),

故

ddx((y′)22+y44)=0.\frac{d}{dx} \left( \frac{(y')^2}{2}+\frac{y^4}{4} \right)=0.

因此

(y′)22+y44=C.\boxed{\frac{(y')^2}{2}+\frac{y^4}{4}=C}.

若原方程式應為 y′′+2y3=0y''+2y^3=0,則題目所列的

(y′)24+y44\frac{(y')^2}{4}+\frac{y^4}{4}

才會是守恆量。


(b)證明 yy 為週期函數

以下使用正確的能量式

(y′)22+y44=C.\frac{(y')^2}{2}+\frac{y^4}{4}=C.

令

y(0)=y0≠0,y′(0)=v0.y(0)=y_0\neq 0,\qquad y'(0)=v_0.

則

C=v022+y044>0.C=\frac{v_0^2}{2}+\frac{y_0^4}{4}>0.

由能量式可得

(y′)2=2C−y42.(y')^2=2C-\frac{y^4}{2}.

定義

A=(4C)1/4>0,A=(4C)^{1/4}>0,

則 4C=A44C=A^4,所以

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

其他考古題