108 年 國立中央大學天文研究所碩士班一般生《應用數學》

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

第 (1) 題20 分

(a) (10 points) Does the equation y=x3−3x+2=(x+2)(x−1)2y = x^3 - 3x + 2 = (x+2)(x-1)^2 has a maximum or a minimum? If yes, find out their locations and values. With these information and the location of the root(s), sketch the plot yy vs xx of the equation.
(b) (6 points) Similarly, sketch the curve y2=x3−3x+2y^2 = x^3 - 3x + 2. (Note that the left hand side is y2y^2.)
(c) (4 points) Without going through lengthy calculation, roughly sketch the curve y2=x3−3x+by^2 = x^3 - 3x + b for (i) bb is a little bit smaller than 2, say 1.9; and (ii) bb is a little bit larger than 2, say 2.1.

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

這一題的完整詳解

本題主要考驗微積分基本概念,包括求極值、畫函數圖形,以及對曲線變形的理解。

(a) 求極值與繪製函數圖形

首先,我們需要找出函數 y=x3−3x+2y = x^3 - 3x + 2 的極值。這需要計算一階導數並令其為零。
y′=dydx=3x2−3y' = \frac{dy}{dx} = 3x^2 - 3
令 y′=0y' = 0,得 3x2−3=0  ⟹  x2=1  ⟹  x=±13x^2 - 3 = 0 \implies x^2 = 1 \implies x = \pm 1。

接著,我們使用二階導數來判斷這些點是極大值還是極小值。
y′′=d2ydx2=6xy'' = \frac{d^2y}{dx^2} = 6x

當 x=1x = 1 時,y′′(1)=6(1)=6>0y''(1) = 6(1) = 6 > 0,所以 x=1x=1 是局部極小值點。
對應的 yy 值為 y(1)=(1)3−3(1)+2=1−3+2=0y(1) = (1)^3 - 3(1) + 2 = 1 - 3 + 2 = 0。
極小值點為 (1,0)(1, 0)。

當 x=−1x = -1 時,y′′(−1)=6(−1)=−6<0y''(-1) = 6(-1) = -6 < 0,所以 x=−1x=-1 是局部極大值點。
對應的 yy 值為 y(−1)=(−1)3−3(−1)+2=−1+3+2=4y(-1) = (-1)^3 - 3(-1) + 2 = -1 + 3 + 2 = 4。
極大值點為 (−1,4)(-1, 4)。

接下來,我們找出函數的根(即 y=0y=0 的點)。
題目已給出 y=(x+2)(x−1)2y = (x+2)(x-1)^2。
令 y=0y=0,則 (x+2)(x−1)2=0(x+2)(x-1)^2 = 0,根為 x=−2x = -2 和 x=1x = 1(重根)。
重根 x=1x=1 代表函數在 x=1x=1 處與 x 軸相切。

結合極值點 (−1,4)(-1, 4)、(1,0)(1, 0),以及根 x=−2x=-2 和 x=1x=1(切點),我們可以在 x-y 平面上繪製函數圖形。
函數在 x→∞x \to \infty 時 y→∞y \to \infty,在 x→−∞x \to -\infty 時 y→−∞y \to -\infty。
圖形大致會從左下角上升,經過 (−2,0)(-2, 0),在 (−1,4)(-1, 4) 達到局部極大值,然後下降經過 (1,0)(1, 0) 處與 x 軸相切,之後再上升。

【答案】
函數有局部極大值點 (−1,4)(-1, 4) 和局部極小值點 (1,0)(1, 0)。
函數的根為 x=−2x = -2 和 x=1x = 1(重根)。
繪製的圖形應包含上述特徵。

(b) 繪製曲線 y2=x3−3x+2y^2 = x^3 - 3x + 2

這個曲線是 y=x3−3x+2y = x^3 - 3x + 2 的對稱變形。
首先,觀察函數 f(x)=x3−3x+2f(x) = x^3 - 3x + 2。我們已經知道它的根是 x=−2x=-2 和 x=1x=1(重根),極大值點是 (−1,4)(-1, 4),極小值點是 (1,0)(1, 0)。

由於曲線是 y2=f(x)y^2 = f(x),我們有以下性質:

  1. 對稱性: 由於 yy 只出現在 y2y^2 的形式中,所以如果 (x,y)(x, y) 在曲線上,那麼 (x,−y)(x, -y) 也一定在曲線上。這表示曲線關於 x 軸對稱。
  2. 定義域: y2≥0y^2 \ge 0,所以我們必須有 f(x)=x3−3x+2≥0f(x) = x^3 - 3x + 2 \ge 0。從 (a) 的分析,我們知道 f(x)≥0f(x) \ge 0 的區域是 x∈[−2,1]∪[1,∞)x \in [-2, 1] \cup [1, \infty),也就是 x≥−2x \ge -2。
  3. y=0 的點: 當 y=0y=0 時,f(x)=0f(x)=0,即 x3−3x+2=0x^3 - 3x + 2 = 0。這發生在 x=−2x=-2 和 x=1x=1。所以曲線與 x 軸的交點是 (−2,0)(-2, 0) 和 (1,0)(1, 0)。
  4. 極值:
    • 當 f(x)f(x) 有局部極大值時,例如在 x=−1x=-1,f(−1)=4f(-1)=4。此時 y2=4y^2 = 4,所以 y=±2y = \pm 2。這產生了兩個點 (−1,2)(-1, 2) 和 (−1,−2)(-1, -2)。
    • 當 f(x)f(x) 有局部極小值時,例如在 x=1x=1,f(1)=0f(1)=0。此時 y2=0y^2 = 0,所以 y=0y=0。這產生了點 (1,0)(1, 0),與 x 軸相切。

繪製步驟:

  • 標出 x 軸交點 (−2,0)(-2, 0) 和 (1,0)(1, 0)。
  • 標出關於 x 軸對稱的點 (−1,2)(-1, 2) 和 (−1,−2)(-1, -2)。
  • 由於 x≥−2x \ge -2,我們從 x=−2x=-2 開始繪製。
  • 在 x=−2x=-2 到 x=−1x=-1 之間,f(x)f(x) 從 0 上升到 4。這意味著 yy 從 0 上升到 ±2\pm 2。
  • 在 x=−1x=-1 到 x=1x=1 之間,f(x)f(x) 從 4 下降到 0。這意味著 yy 從 ±2\pm 2 下降到 0。
  • 在 x>1x > 1 之間,f(x)f(x) 從 0 上升到 ∞\infty。這意味著 y2y^2 從 0 上升到 ∞\infty,所以 yy 從 0 上升到 ±∞\pm \infty。
  • 由於曲線關於 x 軸對稱,我們繪製上半部 (y≥0y \ge 0) 然後鏡射到下半部 (y≤0y \le 0)。

【答案】
曲線關於 x 軸對稱,定義域為 x≥−2x \ge -2。
主要特徵點包括 (−2,0)(-2, 0), (1,0)(1, 0), (−1,2)(-1, 2), (−1,−2)(-1, -2)。
在 x=1x=1 處,曲線與 x 軸相切。
繪製的圖形應反映上述特徵。

🔒

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

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

免費註冊

第 (2) 題20 分

The relation between Cartesian coordinates (x,y,z)(x, y, z) and spherical coordinates (r,θ,ϕ)(r, \theta, \phi) are x=rsin⁡θcos⁡ϕx = r \sin \theta \cos \phi, y=rsin⁡θsin⁡ϕy = r \sin \theta \sin \phi, z=rcos⁡θz = r \cos \theta. Let {e^x,e^y,e^z}\{\hat{e}_x, \hat{e}_y, \hat{e}_z\} be the orthonormal basis (i.e., orthogonal unit vectors) of the Cartesian coordinate system, and {e^r,e^θ,e^ϕ}\{\hat{e}_r, \hat{e}_\theta, \hat{e}_\phi\} be the orthonormal basis in spherical coordinate system in the direction of increasing r,θ,ϕr, \theta, \phi, respectively.

(a) (4 points) Show the relationship of the two coordinate systems in a figure.
(b) (6 points) Express {e^r,e^θ,e^ϕ}\{\hat{e}_r, \hat{e}_\theta, \hat{e}_\phi\} in terms of {e^x,e^y,e^z}\{\hat{e}_x, \hat{e}_y, \hat{e}_z\} and the spherical coordinates.
(c) (10 points) Express the position vector of a point particle in spherical coordinate system. Find the velocity and acceleration components of the particle in spherical coordinates (r,θ,ϕ)(r, \theta, \phi).

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

這一題的完整詳解

核心觀念

本題考查球座標與直角座標的轉換,以及球座標的局部正交單位向量。採用慣例:

  • rr:原點至粒子的距離。
  • θ\theta:與正 zz 軸的夾角,稱為極角。
  • ϕ\phi:粒子在 xyxy 平面投影與正 xx 軸的夾角,稱為方位角。

球座標的單位向量會隨位置改變,因此求速度與加速度時,除了微分 r,θ,ϕr,\theta,\phi,也必須微分 e^r,e^θ,e^ϕ\hat e_r,\hat e_\theta,\hat e_\phi。


(a) 兩種座標系的關係

粒子位置為 P(x,y,z)P(x,y,z),其在 xyxy 平面的投影為 QQ。定義

OP=r,OQ=rsin⁡θ,z=rcos⁡θ.OP=r,\qquad OQ=r\sin\theta,\qquad z=r\cos\theta.

示意圖如下:

                         z
                         ↑
                         │
                         │        P(r,θ,φ)
                         │       ●
                         │      /│
                         │     / │ r
                         │    /θ │
                         │   /   │
                         O──●─────┘────────→ x
                           Q
                          ╱
                         ╱ φ
                        ↙
                       y

其中:

  • e^r\hat e_r 指向半徑增加的方向,即由原點指向粒子。
  • e^θ\hat e_\theta 指向 θ\theta 增加的方向,即沿經線朝遠離正 zz 軸的方向。
  • e^ϕ\hat e_\phi 指向 ϕ\phi 增加的方向,即繞 zz 軸的切線方向。

三個球座標單位向量彼此正交,但其方向隨著 θ,ϕ\theta,\phi 改變。


(b) 球座標單位向量

徑向單位向量 e^r\hat e_r

位置向量可寫成

r⃗=xe^x+ye^y+ze^z.\vec r=x\hat e_x+y\hat e_y+z\hat e_z.

除以長度 rr,得到徑向單位向量:

e^r=sin⁡θcos⁡ϕ e^x+sin⁡θsin⁡ϕ e^y+cos⁡θ e^z.\hat e_r = \sin\theta\cos\phi\,\hat e_x + \sin\theta\sin\phi\,\hat e_y + \cos\theta\,\hat e_z.

極角方向單位向量 e^θ\hat e_\theta

對 e^r\hat e_r 對 θ\theta 微分:

∂e^r∂θ=cos⁡θcos⁡ϕ e^x+cos⁡θsin⁡ϕ e^y−sin⁡θ e^z.\frac{\partial \hat e_r}{\partial\theta} = \cos\theta\cos\phi\,\hat e_x + \cos\theta\sin\phi\,\hat e_y - \sin\theta\,\hat e_z.

此向量的長度為 11,且方向正是 θ\theta 增加的方向,因此

e^θ=cos⁡θcos⁡ϕ e^x+cos⁡θsin⁡ϕ e^y−sin⁡θ e^z.\hat e_\theta = \cos\theta\cos\phi\,\hat e_x + \cos\theta\sin\phi\,\hat e_y - \sin\theta\,\hat e_z.

方位角方向單位向量 e^ϕ\hat e_\phi

對 e^r\hat e_r 對 ϕ\phi 微分:

∂e^r∂ϕ=−sin⁡θsin⁡ϕ e^x+sin⁡θcos⁡ϕ e^y.\frac{\partial \hat e_r}{\partial\phi} = -\sin\theta\sin\phi\,\hat e_x + \sin\theta\cos\phi\,\hat e_y.

其長度為 sin⁡θ\sin\theta,所以除以 sin⁡θ\sin\theta 後得到單位向量:

e^ϕ=−sin⁡ϕ e^x+cos⁡ϕ e^y.\hat e_\phi = -\sin\phi\,\hat e_x + \cos\phi\,\hat e_y.

因此三個單位向量為

e^r=sin⁡θcos⁡ϕ e^x+sin⁡θsin⁡ϕ e^y+cos⁡θ e^z\boxed{ \hat e_r = \sin\theta\cos\phi\,\hat e_x + \sin\theta\sin\phi\,\hat e_y + \cos\theta\,\hat e_z } e^θ=cos⁡θcos⁡ϕ e^x+cos⁡θsin⁡ϕ e^y−sin⁡θ e^z\boxed{ \hat e_\theta = \cos\theta\cos\phi\,\hat e_x + \cos\theta\sin\phi\,\hat e_y - \sin\theta\,\hat e_z } e^ϕ=−sin⁡ϕ e^x+cos⁡ϕ e^y\boxed{ \hat e_\phi = -\sin\phi\,\hat e_x + \cos\phi\,\hat e_y }

快速驗算:

e^r⋅e^θ=e^θ⋅e^ϕ=e^ϕ⋅e^r=0,\hat e_r\cdot\hat e_\theta = \hat e_\theta\cdot\hat e_\phi = \hat e_\phi\cdot\hat e_r =0,

且

∣e^r∣=∣e^θ∣=∣e^ϕ∣=1.|\hat e_r|=|\hat e_\theta|=|\hat e_\phi|=1.

(c) 位置、速度與加速度

位置向量

由直角座標轉換式,

r⃗=xe^x+ye^y+ze^z.\vec r = x\hat e_x+y\hat e_y+z\hat e_z.

代入

x=rsin⁡θcos⁡ϕ,y=rsin⁡θsin⁡ϕ,z=rcos⁡θ,x=r\sin\theta\cos\phi,\qquad y=r\sin\theta\sin\phi,\qquad z=r\cos\theta,

可得

r⃗=r(sin⁡θcos⁡ϕ e^x+sin⁡θsin⁡ϕ e^y+cos⁡θ e^z).\vec r = r\left( \sin\theta\cos\phi\,\hat e_x + \sin\theta\sin\phi\,\hat e_y + \cos\theta\,\hat e_z \right).

因此

r⃗=re^r.\boxed{\vec r=r\hat e_r}.

球座標基底向量的時間微分

由 e^r\hat e_r 對 θ,ϕ\theta,\phi 的偏微分,

∂e^r∂θ=e^θ,∂e^r∂ϕ=sin⁡θ e^ϕ.\frac{\partial\hat e_r}{\partial\theta} =\hat e_\theta, \qquad \frac{\partial\hat e_r}{\partial\phi} =\sin\theta\,\hat e_\phi.

使用鏈鎖律:

de^rdt=∂e^r∂θθ˙+∂e^r∂ϕϕ˙,\frac{d\hat e_r}{dt} = \frac{\partial\hat e_r}{\partial\theta}\dot\theta + \frac{\partial\hat e_r}{\partial\phi}\dot\phi,

所以

de^rdt=θ˙ e^θ+ϕ˙sin⁡θ e^ϕ.\boxed{ \frac{d\hat e_r}{dt} = \dot\theta\,\hat e_\theta + \dot\phi\sin\theta\,\hat e_\phi }.

同理可得

🔒

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

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

免費註冊

第 (3) 題20 分

Consider the matrix

A=(104040101)A = \begin{pmatrix} 1 & 0 & 4 \\ 0 & 4 & 0 \\ 1 & 0 & 1 \end{pmatrix}

(a) (10 points) Find the eigenvalues and the corresponding eigenvectors of the matrix.
(b) (2 points) Are there any pairs of eigenvectors perpendicular to each other? If yes, what are they?
(c) (8 points) Find the inverse of the matrix.

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

這一題的完整詳解

核心觀念

本題評量線性代數中矩陣特徵理論與矩陣求逆運算的三大基本主題:

  1. 特徵值與特徵向量(Eigenvalues and Eigenvectors):
    • 滿足特徵方程式 det⁡(A−λI)=0\det(A - \lambda I) = 0 的標量 λ\lambda 為特徵值。
    • 對應特定特徵值 λ\lambda 的非零向量 v\mathbf{v},滿足 (A−λI)v=0(A - \lambda I)\mathbf{v} = \mathbf{0},構成該特徵值的特徵空間(Eigenspace)。
  2. 正交性(Orthogonality):
    • 在標準內積空間中,兩實向量 u,v\mathbf{u}, \mathbf{v} 若滿足內積 u⋅v=uTv=0\mathbf{u} \cdot \mathbf{v} = \mathbf{u}^T \mathbf{v} = 0,則兩向量互相垂直(正交)。
  3. 反矩陣(Inverse Matrix):
    • 若矩陣 AA 之行列式值 det⁡(A)≠0\det(A) \neq 0,則 AA 可逆。
    • 可利用古典伴隨矩陣法(Adjugate Formula):A−1=1det⁡(A)adj⁡(A)A^{-1} = \frac{1}{\det(A)} \operatorname{adj}(A),或列運算高斯-喬登消去法(Gauss-Jordan Elimination)求得。

解題方法與推導過程

(a) 求特徵值與對應之特徵向量

Step 1:求解特徵多項式

計算特徵方程式 det⁡(A−λI)=0\det(A - \lambda I) = 0:

det⁡(A−λI)=∣1−λ0404−λ0101−λ∣=0\det(A - \lambda I) = \begin{vmatrix} 1 - \lambda & 0 & 4 \\ 0 & 4 - \lambda & 0 \\ 1 & 0 & 1 - \lambda \end{vmatrix} = 0

沿著第二列(第二列僅有一個非零元素 4−λ4 - \lambda)進行拉普拉斯展開:

det⁡(A−λI)=(4−λ)∣1−λ411−λ∣=0\det(A - \lambda I) = (4 - \lambda) \begin{vmatrix} 1 - \lambda & 4 \\ 1 & 1 - \lambda \end{vmatrix} = 0

展開二階行列式:

(4−λ)[(1−λ)2−4]=0(4 - \lambda) \left[ (1 - \lambda)^2 - 4 \right] = 0 (4−λ)(λ2−2λ−3)=0(4 - \lambda) \left( \lambda^2 - 2\lambda - 3 \right) = 0 (4−λ)(λ−3)(λ+1)=0(4 - \lambda)(\lambda - 3)(\lambda + 1) = 0

解得特徵值為:

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

Step 2:求各特徵值對應之特徵向量

  1. 當 λ1=4\lambda_1 = 4 時:
    解方程式 (A−4I)v1=0(A - 4I)\mathbf{v}_1 = \mathbf{0}:

    (−30400010−3)(x1x2x3)=(000)\begin{pmatrix} -3 & 0 & 4 \\ 0 & 0 & 0 \\ 1 & 0 & -3 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}

    由第一式 −3x1+4x3=0-3x_1 + 4x_3 = 0 與第三式 x1−3x3=0x_1 - 3x_3 = 0 可得聯立解:

    x1=0,x3=0x_1 = 0, \quad x_3 = 0

    變數 x2x_2 為自由變數(Free variable)。令 x2=1x_2 = 1(或任意非零實數 c1c_1),對應之特徵向量為:

    v1=(010)\mathbf{v}_1 = \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix}
  2. 當 λ2=3\lambda_2 = 3 時:
    解方程式 (A−3I)v2=0(A - 3I)\mathbf{v}_2 = \mathbf{0}:

    (−20401010−2)(x1x2x3)=(000)\begin{pmatrix} -2 & 0 & 4 \\ 0 & 1 & 0 \\ 1 & 0 & -2 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}

    由第二列直接得 x2=0x_2 = 0;
    由第一列及第三列得 −2x1+4x3=0  ⟹  x1=2x3-2x_1 + 4x_3 = 0 \implies x_1 = 2x_3。
    令 x3=1x_3 = 1(或任意非零實數 c2c_2),對應之特徵向量為:

    v2=(201)\mathbf{v}_2 = \begin{pmatrix} 2 \\ 0 \\ 1 \end{pmatrix}
  3. 當 λ3=−1\lambda_3 = -1 時:
    解方程式 (A−(−1)I)v3=(A+I)v3=0(A - (-1)I)\mathbf{v}_3 = (A + I)\mathbf{v}_3 = \mathbf{0}:

    (204050102)(x1x2x3)=(000)\begin{pmatrix} 2 & 0 & 4 \\ 0 & 5 & 0 \\ 1 & 0 & 2 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}

    由第二列得 5x2=0  ⟹  x2=05x_2 = 0 \implies x_2 = 0;
    由第一列與第三列得 x1+2x3=0  ⟹  x1=−2x3x_1 + 2x_3 = 0 \implies x_1 = -2x_3。
    令 x3=1x_3 = 1(或任意非零實數 c3c_3),對應之特徵向量為:

    v3=(−201)\mathbf{v}_3 = \begin{pmatrix} -2 \\ 0 \\ 1 \end{pmatrix}

(b) 檢驗特徵向量之間的正交性

分別計算三組特徵向量之間的兩兩點積(內積):

  1. v1\mathbf{v}_1 與 v2\mathbf{v}_2:

    v1⋅v2=(0)(2)+(1)(0)+(0)(1)=0\mathbf{v}_1 \cdot \mathbf{v}_2 = (0)(2) + (1)(0) + (0)(1) = 0

    因此 v1⊥v2\mathbf{v}_1 \perp \mathbf{v}_2。

  2. v1\mathbf{v}_1 與 v3\mathbf{v}_3:

    v1⋅v3=(0)(−2)+(1)(0)+(0)(1)=0\mathbf{v}_1 \cdot \mathbf{v}_3 = (0)(-2) + (1)(0) + (0)(1) = 0

    因此 v1⊥v3\mathbf{v}_1 \perp \mathbf{v}_3。

  3. v2\mathbf{v}_2 與 v3\mathbf{v}_3:

    v2⋅v3=(2)(−2)+(0)(0)+(1)(1)=−4+0+1=−3≠0\mathbf{v}_2 \cdot \mathbf{v}_3 = (2)(-2) + (0)(0) + (1)(1) = -4 + 0 + 1 = -3 \neq 0
🔒

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

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

免費註冊

第 (4) 題20 分

(a) (5 points) A force can be expressed as the gradient of a potential, i.e., F=−∇ϕ\mathbf{F} = -\nabla \phi is called a conservative force. Show that ∇×F=0\nabla \times \mathbf{F} = 0.
(b) (10 points) The centre of the Earth is located at the origin of a Cartesian coordinate system. The centre of the Moon is located on the z-axis at a fixed distance R from the origin. The tidal force exerted by the Moon on a point mass at the surface of the Earth (x,y,z)(x, y, z) can be approximated by (when x,y,zx, y, z are much smaller than R)

Fx=−GMmR3x,Fy=−GMmR3y,Fz=2GMmR3zF_x = -\frac{GMm}{R^3}x, \quad F_y = -\frac{GMm}{R^3}y, \quad F_z = \frac{2GMm}{R^3}z

where mm is the mass of the Moon and MM is the mass of the Earth. Show that the tidal force is a conservative force, and work out the corresponding potential.
(c) (5 points) Sketch the contours of the potential surfaces on the x−zx-z plane. Don't forget to label the relative levels of the contours (e.g., which one corresponds to is higher potential and which one to lower).

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

這一題的完整詳解

本題主要考驗向量微積分中的保守力、旋度、梯度,以及重力位能的概念。

(a) 證明保守力旋度為零

一個力場 F\mathbf{F} 被定義為一個勢能函數 ϕ\phi 的梯度負值,即 F=−∇ϕ\mathbf{F} = -\nabla \phi,則稱 F\mathbf{F} 為保守力。
我們需要證明其旋度 ∇×F=0\nabla \times \mathbf{F} = \mathbf{0}。

設 F=Fxe^x+Fye^y+Fze^z\mathbf{F} = F_x \hat{e}_x + F_y \hat{e}_y + F_z \hat{e}_z。
根據定義, Fx=−∂ϕ∂xF_x = -\frac{\partial \phi}{\partial x}, Fy=−∂ϕ∂yF_y = -\frac{\partial \phi}{\partial y}, Fz=−∂ϕ∂zF_z = -\frac{\partial \phi}{\partial z}。

旋度的定義為:

∇×F=∣e^xe^ye^z∂∂x∂∂y∂∂zFxFyFz∣\nabla \times \mathbf{F} = \begin{vmatrix} \hat{e}_x & \hat{e}_y & \hat{e}_z \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ F_x & F_y & F_z \end{vmatrix} =(∂Fz∂y−∂Fy∂z)e^x+(∂Fx∂z−∂Fz∂x)e^y+(∂Fy∂x−∂Fx∂y)e^z= \left( \frac{\partial F_z}{\partial y} - \frac{\partial F_y}{\partial z} \right) \hat{e}_x + \left( \frac{\partial F_x}{\partial z} - \frac{\partial F_z}{\partial x} \right) \hat{e}_y + \left( \frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y} \right) \hat{e}_z

我們來計算每一項的係數:

e^x\hat{e}_x 的係數:∂Fz∂y−∂Fy∂z\frac{\partial F_z}{\partial y} - \frac{\partial F_y}{\partial z}
=∂∂y(−∂ϕ∂z)−∂∂z(−∂ϕ∂y)= \frac{\partial}{\partial y} \left( -\frac{\partial \phi}{\partial z} \right) - \frac{\partial}{\partial z} \left( -\frac{\partial \phi}{\partial y} \right)
=−∂2ϕ∂y∂z+∂2ϕ∂z∂y= -\frac{\partial^2 \phi}{\partial y \partial z} + \frac{\partial^2 \phi}{\partial z \partial y}
由於二階偏導數在連續區域是連續的,所以 ∂2ϕ∂y∂z=∂2ϕ∂z∂y\frac{\partial^2 \phi}{\partial y \partial z} = \frac{\partial^2 \phi}{\partial z \partial y}。
因此,e^x\hat{e}_x 的係數為 0。

e^y\hat{e}_y 的係數:∂Fx∂z−∂Fz∂x\frac{\partial F_x}{\partial z} - \frac{\partial F_z}{\partial x}
=∂∂z(−∂ϕ∂x)−∂∂x(−∂ϕ∂z)= \frac{\partial}{\partial z} \left( -\frac{\partial \phi}{\partial x} \right) - \frac{\partial}{\partial x} \left( -\frac{\partial \phi}{\partial z} \right)
=−∂2ϕ∂z∂x+∂2ϕ∂x∂z= -\frac{\partial^2 \phi}{\partial z \partial x} + \frac{\partial^2 \phi}{\partial x \partial z}
同樣,由於二階偏導數的連續性,此項也為 0。

e^z\hat{e}_z 的係數:∂Fy∂x−∂Fx∂y\frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}
=∂∂x(−∂ϕ∂y)−∂∂y(−∂ϕ∂x)= \frac{\partial}{\partial x} \left( -\frac{\partial \phi}{\partial y} \right) - \frac{\partial}{\partial y} \left( -\frac{\partial \phi}{\partial x} \right)
=−∂2ϕ∂x∂y+∂2ϕ∂y∂x= -\frac{\partial^2 \phi}{\partial x \partial y} + \frac{\partial^2 \phi}{\partial y \partial x}
此項也為 0。

所以,∇×F=0e^x+0e^y+0e^z=0\nabla \times \mathbf{F} = 0 \hat{e}_x + 0 \hat{e}_y + 0 \hat{e}_z = \mathbf{0}。
這證明了任何保守力場的旋度都為零。

【答案】
通過計算,所有旋度的分量均為零,因此 ∇×F=0\nabla \times \mathbf{F} = \mathbf{0}。

(b) 證明潮汐力是保守力並求勢能

給定的潮汐力分量為:
Fx=−GMmR3xF_x = -\frac{GMm}{R^3}x
Fy=−GMmR3yF_y = -\frac{GMm}{R^3}y
Fz=2GMmR3zF_z = \frac{2GMm}{R^3}z

首先,我們需要證明這個力場是保守的,即證明 ∇×F=0\nabla \times \mathbf{F} = \mathbf{0}。
設 k=GMmR3k = \frac{GMm}{R^3}。則 F=(−kx,−ky,2kz)\mathbf{F} = (-kx, -ky, 2kz)。

計算旋度:
∂Fz∂y=∂(2kz)∂y=0\frac{\partial F_z}{\partial y} = \frac{\partial (2kz)}{\partial y} = 0
∂Fy∂z=∂(−ky)∂z=0\frac{\partial F_y}{\partial z} = \frac{\partial (-ky)}{\partial z} = 0
所以 e^x\hat{e}_x 的係數為 0−0=00 - 0 = 0。

∂Fx∂z=∂(−kx)∂z=0\frac{\partial F_x}{\partial z} = \frac{\partial (-kx)}{\partial z} = 0
∂Fz∂x=∂(2kz)∂x=0\frac{\partial F_z}{\partial x} = \frac{\partial (2kz)}{\partial x} = 0
所以 e^y\hat{e}_y 的係數為 0−0=00 - 0 = 0。

∂Fy∂x=∂(−ky)∂x=0\frac{\partial F_y}{\partial x} = \frac{\partial (-ky)}{\partial x} = 0
∂Fx∂y=∂(−kx)∂y=0\frac{\partial F_x}{\partial y} = \frac{\partial (-kx)}{\partial y} = 0
所以 e^z\hat{e}_z 的係數為 0−0=00 - 0 = 0。

因此,∇×F=0\nabla \times \mathbf{F} = \mathbf{0}。這證明了潮汐力是一個保守力場。

接下來,我們需要找到對應的勢能 ϕ\phi。根據保守力的定義,F=−∇ϕ\mathbf{F} = -\nabla \phi。
Fx=−∂ϕ∂x=−kx  ⟹  ∂ϕ∂x=kxF_x = -\frac{\partial \phi}{\partial x} = -kx \implies \frac{\partial \phi}{\partial x} = kx
Fy=−∂ϕ∂y=−ky  ⟹  ∂ϕ∂y=kyF_y = -\frac{\partial \phi}{\partial y} = -ky \implies \frac{\partial \phi}{\partial y} = ky
Fz=−∂ϕ∂z=2kz  ⟹  ∂ϕ∂z=−2kzF_z = -\frac{\partial \phi}{\partial z} = 2kz \implies \frac{\partial \phi}{\partial z} = -2kz

對第一個方程積分:
ϕ(x,y,z)=∫kx dx=12kx2+f(y,z)\phi(x, y, z) = \int kx \, dx = \frac{1}{2} kx^2 + f(y, z)
其中 f(y,z)f(y, z) 是積分常數,它可能是 yy 和 zz 的函數。

對第二個方程積分:
ϕ(x,y,z)=∫ky dy=12ky2+g(x,z)\phi(x, y, z) = \int ky \, dy = \frac{1}{2} ky^2 + g(x, z)
其中 g(x,z)g(x, z) 是積分常數,它可能是 xx 和 zz 的函數。

將第一個積分結果代入第二個方程的條件:
∂∂y(12kx2+f(y,z))=ky\frac{\partial}{\partial y} (\frac{1}{2} kx^2 + f(y, z)) = ky
∂f(y,z)∂y=ky\frac{\partial f(y, z)}{\partial y} = ky
對 f(y,z)f(y, z) 進行積分:

🔒

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

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

免費註冊

第 (5) 題20 分

Given an inhomogeneous ODE (ordinary differential equation)

x¨+x˙t−xt2=δ(t−t0)\ddot{x} + \frac{\dot{x}}{t} - \frac{x}{t^2} = \delta(t - t_0)

where x˙=dx/dt\dot{x} = dx/dt, x¨=d2x/dt2\ddot{x} = d^2x/dt^2, δ(t)\delta(t) is the Dirac delta function and t0t_0 is a constant.

(a) (5 points) If you are seeking continuous solution for x(t)x(t), what kind of matching condition(s) will you suggest for x˙\dot{x} at time t=t0t = t_0. Please state your reasons.
(b) (15 points) Consider t0>0t_0 > 0. If the initial conditions for xx is x(0)=0x(0) = 0 and x˙(0)=1\dot{x}(0) = 1, solve the ODE for xx in the two domains (i) 0<t<t00 < t < t_0 and (ii) t0<tt_0 < t.

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

這一題的完整詳解

本題考驗常微分方程的求解,特別是包含狄拉克 δ 函數的非齊次方程,以及如何處理匹配條件。

(a) 連續解的匹配條件

考慮方程 x¨+x˙t−xt2=δ(t−t0)\ddot{x} + \frac{\dot{x}}{t} - \frac{x}{t^2} = \delta(t - t_0)。
狄拉克 δ 函數 δ(t−t0)\delta(t-t_0) 在 t=t0t=t_0 處有一個衝激。
我們尋求一個連續的解 x(t)x(t)。

首先,觀察方程的齊次部分:x¨+x˙t−xt2=0\ddot{x} + \frac{\dot{x}}{t} - \frac{x}{t^2} = 0。
這個方程的結構類似於 Euler-Cauchy 方程,但多了 x˙t\frac{\dot{x}}{t} 這一項。
如果我們假設解的形式是 x(t)=trx(t) = t^r,則 x˙=rtr−1\dot{x} = rt^{r-1},x¨=r(r−1)tr−2\ddot{x} = r(r-1)t^{r-2}。
代入齊次方程:
r(r−1)tr−2+rtr−1t−trt2=0r(r-1)t^{r-2} + \frac{rt^{r-1}}{t} - \frac{t^r}{t^2} = 0
r(r−1)tr−2+rtr−2−tr−2=0r(r-1)t^{r-2} + rt^{r-2} - t^{r-2} = 0
tr−2[r(r−1)+r−1]=0t^{r-2} [r(r-1) + r - 1] = 0
r2−r+r−1=0r^2 - r + r - 1 = 0
r2−1=0  ⟹  r=±1r^2 - 1 = 0 \implies r = \pm 1。
所以,齊次方程的兩個獨立解是 xh(t)=c1t+c2t−1x_h(t) = c_1 t + c_2 t^{-1}。
這意味著,在 t=0t=0 附近,齊次解 c2t−1c_2 t^{-1} 是不連續的(趨於無窮大)。
然而,題目中給定的初始條件是 x(0)=0x(0)=0 和 x˙(0)=1\dot{x}(0)=1。這意味著解在 t=0t=0 處必須是有限的。
如果 x(t)=c1t+c2t−1x(t) = c_1 t + c_2 t^{-1},那麼 x˙(t)=c1−c2t−2\dot{x}(t) = c_1 - c_2 t^{-2}。
為了使 x(0)x(0) 有限, c2c_2 必須為 0。
為了使 x˙(0)\dot{x}(0) 有限, c2c_2 也必須為 0。
所以,在 t=0t=0 附近,齊次解必須是 xh(t)=c1tx_h(t) = c_1 t。
這與初始條件 x(0)=0x(0)=0 和 x˙(0)=1\dot{x}(0)=1 相符,因為這會得到 x(t)=tx(t)=t (當沒有 δ 函數時)。

現在考慮 δ 函數的影響。狄拉克 δ 函數表示在 t=t0t=t_0 處有一個無限大的瞬時「力」或「衝激」。
對於一個連續的解 x(t)x(t),我們要求 x(t)x(t) 在 t=t0t=t_0 處是連續的。
然而,δ 函數會導致 x¨\ddot{x} 在 t=t0t=t_0 處有一個「跳躍」。
我們可以通過對原方程積分來分析這種跳躍。

對方程兩邊從 t0−ϵt_0 - \epsilon 到 t0+ϵt_0 + \epsilon 積分,其中 ϵ→0+\epsilon \to 0^+。

∫t0−ϵt0+ϵ(x¨+x˙t−xt2)dt=∫t0−ϵt0+ϵδ(t−t0)dt\int_{t_0-\epsilon}^{t_0+\epsilon} \left( \ddot{x} + \frac{\dot{x}}{t} - \frac{x}{t^2} \right) dt = \int_{t_0-\epsilon}^{t_0+\epsilon} \delta(t - t_0) dt ∫t0−ϵt0+ϵx¨dt+∫t0−ϵt0+ϵx˙tdt−∫t0−ϵt0+ϵxt2dt=1\int_{t_0-\epsilon}^{t_0+\epsilon} \ddot{x} dt + \int_{t_0-\epsilon}^{t_0+\epsilon} \frac{\dot{x}}{t} dt - \int_{t_0-\epsilon}^{t_0+\epsilon} \frac{x}{t^2} dt = 1

當 ϵ→0+\epsilon \to 0^+ 時:
∫t0−ϵt0+ϵx¨dt=[x˙]t0−ϵt0+ϵ=x˙(t0+)−x˙(t0−)\int_{t_0-\epsilon}^{t_0+\epsilon} \ddot{x} dt = [\dot{x}]_{t_0-\epsilon}^{t_0+\epsilon} = \dot{x}(t_0^+) - \dot{x}(t_0^-)
∫t0−ϵt0+ϵx˙tdt≈x˙(t0)t0∫t0−ϵt0+ϵdt=x˙(t0)t0(2ϵ)→0\int_{t_0-\epsilon}^{t_0+\epsilon} \frac{\dot{x}}{t} dt \approx \frac{\dot{x}(t_0)}{t_0} \int_{t_0-\epsilon}^{t_0+\epsilon} dt = \frac{\dot{x}(t_0)}{t_0} (2\epsilon) \to 0
∫t0−ϵt0+ϵxt2dt≈x(t0)t02∫t0−ϵt0+ϵdt=x(t0)t02(2ϵ)→0\int_{t_0-\epsilon}^{t_0+\epsilon} \frac{x}{t^2} dt \approx \frac{x(t_0)}{t_0^2} \int_{t_0-\epsilon}^{t_0+\epsilon} dt = \frac{x(t_0)}{t_0^2} (2\epsilon) \to 0

所以,主要的貢獻來自 x¨\ddot{x} 的積分:
x˙(t0+)−x˙(t0−)=1\dot{x}(t_0^+) - \dot{x}(t_0^-) = 1

這意味著,在 t=t0t=t_0 處,一階導數 x˙\dot{x} 有一個大小為 1 的跳躍。
題目要求的是「連續解 for x(t)x(t)」,這意味著 x(t)x(t) 本身在 t=t0t=t_0 處是連續的,即 x(t0+)=x(t0−)x(t_0^+) = x(t_0^-)。
但是,對於 x˙\dot{x},我們發現它有一個跳躍。

然而,題目問的是 "what kind of matching condition(s) will you suggest for x˙\dot{x} at time t=t0t = t_0."
如果我們尋求一個「廣義解」 (weak solution),那麼 x˙\dot{x} 的跳躍是允許的。
但如果我們尋求的是一個「經典解」 (classical solution),那麼 x¨\ddot{x} 必須是良態的 (well-behaved),並且 x˙\dot{x} 必須連續。
但是,δ 函數的引入本身就暗示了我們需要處理某種形式的不連續性。

通常,對於帶有 δ 函數的 ODE,我們允許解的某些導數有跳躍。
對於 x¨+P(t)x˙+Q(t)x=f(t)δ(t−t0)\ddot{x} + P(t)\dot{x} + Q(t)x = f(t)\delta(t-t_0) 類型的方程:
如果 f(t0)≠0f(t_0) \ne 0 且 P(t0)≠0P(t_0) \ne 0,那麼 x˙\dot{x} 會有跳躍。
如果 P(t0)=0P(t_0) = 0 且 f(t0)≠0f(t_0) \ne 0,那麼 x¨\ddot{x} 會有跳躍。

在我們的方程 x¨+1tx˙−1t2x=δ(t−t0)\ddot{x} + \frac{1}{t}\dot{x} - \frac{1}{t^2}x = \delta(t - t_0) 中,
P(t)=1/tP(t) = 1/t
Q(t)=−1/t2Q(t) = -1/t^2
f(t)=1f(t) = 1
在 t=t0t=t_0 處,P(t0)=1/t0P(t_0) = 1/t_0。
如果 t0≠0t_0 \ne 0,那麼 P(t0)≠0P(t_0) \ne 0。
因此,我們預期 x˙\dot{x} 會有跳躍,其大小為 f(t0)=1f(t_0) = 1。
即 x˙(t0+)−x˙(t0−)=1\dot{x}(t_0^+) - \dot{x}(t_0^-) = 1。

所以,匹配條件是 x˙(t)\dot{x}(t) 在 t=t0t=t_0 處有一個大小為 1 的跳躍。

理由:
狄拉克 δ 函數代表一個無限大的瞬時作用。當它出現在方程的右側時,它會導致解的某些導數發生跳躍。對於形如 x¨+P(t)x˙+Q(t)x=f(t)δ(t−t0)\ddot{x} + P(t)\dot{x} + Q(t)x = f(t)\delta(t-t_0) 的方程,如果 P(t0)≠0P(t_0) \ne 0 且 f(t0)≠0f(t_0) \ne 0,則

🔒

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

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

免費註冊

其他考古題

108 年中央大學的其他科目