112 年 國立聯合大學電機工程學系研究所《工程數學》

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

第 1 題

The matrix A is as below.

A=[40−2010101]A = \begin{bmatrix} 4 & 0 & -2 \\ 0 & 1 & 0 \\ 1 & 0 & 1 \end{bmatrix}

(1). Find the eigenvalues λ1,λ2,λ3\lambda_1, \lambda_2, \lambda_3. (4%)
(2). Find the corresponding eigenvectors X1,X2,X3X_1, X_2, X_3 with matrix X=[X1X2X3]X = [X_1 X_2 X_3]. (6%)

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 1 頁

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

這一題的完整詳解

核心觀念

本題考查矩陣的特徵值與特徵向量:

  • 特徵值 λ\lambda 滿足
    det⁡(A−λI)=0\det(A-\lambda I)=0
  • 對應於特徵值 λ\lambda 的特徵向量 XX 滿足
    AX=λXAX=\lambda X
    或等價地
    (A−λI)X=0(A-\lambda I)X=0

解題方法

給定

A=[40−2010101]A= \begin{bmatrix} 4&0&-2\\ 0&1&0\\ 1&0&1 \end{bmatrix}

(1) 求特徵值

計算特徵方程:

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

沿著第二列或第二欄展開:

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

因此

det⁡(A−λI)=(1−λ)[(4−λ)(1−λ)+2]\det(A-\lambda I) = (1-\lambda)\left[(4-\lambda)(1-\lambda)+2\right]

整理括號內:

(4−λ)(1−λ)+2=λ2−5λ+6=(λ−2)(λ−3)(4-\lambda)(1-\lambda)+2 = \lambda^2-5\lambda+6 = (\lambda-2)(\lambda-3)

所以

det⁡(A−λI)=(1−λ)(λ−2)(λ−3)\det(A-\lambda I) = (1-\lambda)(\lambda-2)(\lambda-3)

令其等於零,得到三個特徵值:

λ1=1,λ2=2,λ3=3\boxed{\lambda_1=1,\qquad \lambda_2=2,\qquad \lambda_3=3}

(2) 求對應的特徵向量

令

X=[xyz]X= \begin{bmatrix} x\\y\\z \end{bmatrix}

對應 λ1=1\lambda_1=1

(A−I)X=0(A-I)X=0

即

[30−2000100][xyz]=[000]\begin{bmatrix} 3&0&-2\\ 0&0&0\\ 1&0&0 \end{bmatrix} \begin{bmatrix} x\\y\\z \end{bmatrix} = \begin{bmatrix} 0\\0\\0 \end{bmatrix}

由第三列得 x=0x=0,再由第一列得 z=0z=0,而 yy 可自由取值。因此可取

X1=[010]X_1= \begin{bmatrix} 0\\1\\0 \end{bmatrix}
🔒

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

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

免費註冊

第 2 題

Solve the initial value problem of the following differential equations.
(1). 2y(x2−1)y′+(y2+1)ln⁡∣y2+1∣=02y(x^2-1)y'+(y^2+1)\ln|y^2+1| = 0, y(3)=1y(3) =1 (10%)
(2). x2y′′−xy′+4y=cos⁡(ln⁡∣x∣)+xsin⁡(ln⁡∣x∣)x^2y''-xy'+4y = \cos(\ln|x|) + x \sin(\ln|x|), y(1)=1,y′(1)=2y(1) = 1, y'(1) = 2 (10%)
(3). y′+y′+(x+1)y=ex2y3y'+ y'+ (x+1)y = e^{x^2}y^3, y(0)=4y(0) = 4 (10%)
(4). y′′′−9y′′+27y′−27y=54sin⁡3xy'''-9y''+27y'-27y = 54\sin 3x, y(0)=1,y′(0)=2,y′′(0)=3y(0) = 1, y'(0) = 2, y''(0) = 3 (10%)

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

這一題的完整詳解

(1)

2y(x2−1)y′+(y2+1)ln⁡(y2+1)=0,y(3)=12y(x^{2}-1)y'+(y^{2}+1)\ln (y^{2}+1)=0,\qquad y(3)=1

2y(y2+1)ln⁡(y2+1) dy=−dxx2−1\frac{2y}{(y^{2}+1)\ln (y^{2}+1)}\,dy=-\frac{dx}{x^{2}-1}

左側以 u=y2+1u=y^{2}+1 代換得 ∫duuln⁡u=ln⁡∣ln⁡u∣\displaystyle\int\frac{du}{u\ln u}= \ln|\ln u|。
右側 ∫−dxx2−1=−12ln⁡∣x−1x+1∣\displaystyle\int-\frac{dx}{x^{2}-1}= -\frac12\ln\Big|\frac{x-1}{x+1}\Big|。

ln⁡∣ln⁡(y2+1)∣=−12ln⁡∣x−1x+1∣+C\ln|\ln (y^{2}+1)|=-\frac12\ln\Big|\frac{x-1}{x+1}\Big|+C

取正號(y>0y>0),指數化後

ln⁡(y2+1)=Kx+1x−1,K=ln⁡22  (x=3,  y=1)\ln (y^{2}+1)=K\sqrt{\frac{x+1}{x-1}},\qquad K=\frac{\ln2}{\sqrt2}\;(x=3,\;y=1)

 y(x)=exp⁡ ⁣(ln⁡22x+1x−1)−1 \boxed{\,y(x)=\sqrt{\exp\!\left(\frac{\ln2}{\sqrt2}\sqrt{\frac{x+1}{x-1}}\right)-1}\,}

【答案】y(x)=eln⁡22x+1x−1−1y(x)=\sqrt{e^{\frac{\ln2}{\sqrt2}\sqrt{\frac{x+1}{x-1}}}-1}


(2)

x2y′′−xy′+4y=cos⁡(ln⁡∣x∣)+xsin⁡(ln⁡∣x∣),y(1)=1,  y′(1)=2x^{2}y''-xy'+4y=\cos(\ln|x|)+x\sin(\ln|x|),\qquad y(1)=1,\;y'(1)=2

令 t=ln⁡∣x∣  (x=et)t=\ln|x|\;(x=e^{t})。則

ytt−2yt+4y=cos⁡t+etsin⁡ty_{tt}-2y_{t}+4y=\cos t+e^{t}\sin t

齊次解:yh=et[C1cos⁡(3 t)+C2sin⁡(3 t)]y_{h}=e^{t}\big[C_{1}\cos(\sqrt3\,t)+C_{2}\sin(\sqrt3\,t)\big]。

特解

yp=Acos⁡t+Bsin⁡t+12etsin⁡t,A=35,  B=−25y_{p}=A\cos t+B\sin t+\tfrac12e^{t}\sin t, \quad A=\frac35,\;B=-\frac25

合成之

y(t)=et[C1cos⁡(3t)+C2sin⁡(3t)]+35cos⁡t−25sin⁡t+12etsin⁡ty(t)=e^{t}\big[C_{1}\cos(\sqrt3t)+C_{2}\sin(\sqrt3t)\big] +\frac35\cos t-\frac25\sin t+\frac12e^{t}\sin t

初值 t=0t=0(x=1x=1)給
C1=2/5,  C2=3/2C_{1}=2/5,\;C_{2}=\sqrt3/2。

回代 t=ln⁡xt=\ln x 並以 et=xe^{t}=x 得

🔒

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

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

免費註冊

第 3 題

Find the solutions of the following differential equations by use of Laplace transforms.
(1). y′′+4y′+3y=sin⁡(t−1)u(t−1)+e−2tδ(t−2)y''+4y'+3y = \sin(t-1)u(t-1)+e^{-2t}\delta(t-2), y(0)=0,y′(0)=0y(0) = 0, y'(0) = 0 (10%)
(2). y(t)+∫0ty(τ)cos⁡(t−τ)dτ=tu(t−1)+e−2ty(t) + \int_0^t y(\tau) \cos(t-\tau)d\tau = t u(t-1)+e^{-2t} (10%)
(3).

{y1′=4y1+8y2+2cos⁡t−16sin⁡ty2′=6y1+2y2+cos⁡t−14sin⁡t,y1(0)=15,y2(0)=13(10%)\begin{cases} y_1' = 4y_1 + 8y_2 + 2 \cos t - 16 \sin t \\ y_2' = 6y_1 + 2y_2 + \cos t - 14 \sin t \end{cases} , \quad y_1(0) = 15, y_2(0) = 13 \quad (10\%)

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

這一題的完整詳解

核心觀念

本題統一考查拉普拉斯轉換的下列性質:

  • 微分:
    L{y′}=sY(s)−y(0)\mathcal{L}\{y'\}=sY(s)-y(0)
    L{y′′}=s2Y(s)−sy(0)−y′(0)\mathcal{L}\{y''\}=s^2Y(s)-sy(0)-y'(0)

  • 位移函數:
    L{f(t−a)u(t−a)}=e−asF(s)\mathcal{L}\{f(t-a)u(t-a)\}=e^{-as}F(s)

  • Dirac 脈衝:
    L{δ(t−a)}=e−as\mathcal{L}\{\delta(t-a)\}=e^{-as}
    g(t)δ(t−a)=g(a)δ(t−a)g(t)\delta(t-a)=g(a)\delta(t-a)

  • 卷積定理:
    L{∫0tf(τ)g(t−τ) dτ}=F(s)G(s)\mathcal{L}\left\{\int_0^t f(\tau)g(t-\tau)\,d\tau\right\}=F(s)G(s)

  • 聯立微分方程式可利用矩陣或逐式拉普拉斯轉換求解。


(1)

題目為

y′′+4y′+3y=sin⁡(t−1)u(t−1)+e−2tδ(t−2),y(0)=0,y′(0)=0.y''+4y'+3y=\sin(t-1)u(t-1)+e^{-2t}\delta(t-2), \qquad y(0)=0,\quad y'(0)=0.

解題方法

先對兩邊取拉普拉斯轉換。

左側因初始條件皆為零:

L{y′′+4y′+3y}=(s2+4s+3)Y(s).\mathcal{L}\{y''+4y'+3y\} =(s^2+4s+3)Y(s).

右側第一項為第二平移:

L{sin⁡(t−1)u(t−1)}=e−s1s2+1.\mathcal{L}\{\sin(t-1)u(t-1)\} =e^{-s}\frac{1}{s^2+1}.

第二項利用脈衝函數性質:

e−2tδ(t−2)=e−4δ(t−2),e^{-2t}\delta(t-2)=e^{-4}\delta(t-2),

因此

L{e−2tδ(t−2)}=e−4e−2s.\mathcal{L}\{e^{-2t}\delta(t-2)\} =e^{-4}e^{-2s}.

所以

(s2+4s+3)Y(s)=e−ss2+1+e−4e−2s.(s^2+4s+3)Y(s) =\frac{e^{-s}}{s^2+1}+e^{-4}e^{-2s}.

因為

s2+4s+3=(s+1)(s+3),s^2+4s+3=(s+1)(s+3),

得到

Y(s)=e−s1(s2+1)(s+1)(s+3)+e−4e−2s1(s+1)(s+3).Y(s) =e^{-s}\frac{1}{(s^2+1)(s+1)(s+3)} +e^{-4}e^{-2s}\frac{1}{(s+1)(s+3)}.

第一部分反轉換

作部分分式分解:

1(s2+1)(s+1)(s+3)=14(s+1)−120(s+3)−s5(s2+1)+110(s2+1).\frac{1}{(s^2+1)(s+1)(s+3)} = \frac{1}{4(s+1)} -\frac{1}{20(s+3)} -\frac{s}{5(s^2+1)} +\frac{1}{10(s^2+1)}.

因此

f(t)=L−1{1(s2+1)(s+1)(s+3)}f(t)=\mathcal{L}^{-1} \left\{ \frac{1}{(s^2+1)(s+1)(s+3)} \right\}

為

f(t)=14e−t−120e−3t−15cos⁡t+110sin⁡t.f(t) =\frac14e^{-t}-\frac1{20}e^{-3t} -\frac15\cos t+\frac1{10}\sin t.

由第一平移性質:

L−1{e−sF(s)}=f(t−1)u(t−1).\mathcal{L}^{-1} \left\{ e^{-s}F(s) \right\} =f(t-1)u(t-1).

故第一部分為

[14e−(t−1)−120e−3(t−1)−15cos⁡(t−1)+110sin⁡(t−1)]u(t−1).\left[ \frac14e^{-(t-1)} -\frac1{20}e^{-3(t-1)} -\frac15\cos(t-1) +\frac1{10}\sin(t-1) \right]u(t-1).

第二部分反轉換

又有

1(s+1)(s+3)=12(1s+1−1s+3),\frac{1}{(s+1)(s+3)} =\frac12\left(\frac1{s+1}-\frac1{s+3}\right),

所以

L−1{1(s+1)(s+3)}=12(e−t−e−3t).\mathcal{L}^{-1} \left\{ \frac{1}{(s+1)(s+3)} \right\} =\frac12(e^{-t}-e^{-3t}).

考慮平移量 22 與係數 e−4e^{-4}:

L−1{e−4e−2s1(s+1)(s+3)}=e−42[e−(t−2)−e−3(t−2)]u(t−2).\mathcal{L}^{-1} \left\{ e^{-4}e^{-2s}\frac{1}{(s+1)(s+3)} \right\} = \frac{e^{-4}}2 \left[ e^{-(t-2)}-e^{-3(t-2)} \right]u(t-2).

因此

y(t)=[14e−(t−1)−120e−3(t−1)−15cos⁡(t−1)+110sin⁡(t−1)]u(t−1)+e−42[e−(t−2)−e−3(t−2)]u(t−2).\boxed{ \begin{aligned} y(t) ={}& \left[ \frac14e^{-(t-1)} -\frac1{20}e^{-3(t-1)} -\frac15\cos(t-1) +\frac1{10}\sin(t-1) \right]u(t-1) \\ &+\frac{e^{-4}}2 \left[ e^{-(t-2)}-e^{-3(t-2)} \right]u(t-2). \end{aligned} }

解題技巧

脈衝項的係數必須先將指數函數在 t=2t=2 取值:

e−2tδ(t−2)=e−4δ(t−2).e^{-2t}\delta(t-2)=e^{-4}\delta(t-2).

不能誤寫成 e−2δ(t−2)e^{-2} \delta(t-2)。另外,e−ase^{-as} 對應的時間位移是 t−at-a,並且必須乘上 u(t−a)u(t-a)。


(2)

題目為

y(t)+∫0ty(τ)cos⁡(t−τ) dτ=tu(t−1)+e−2t.y(t)+\int_0^t y(\tau)\cos(t-\tau)\,d\tau =tu(t-1)+e^{-2t}.

解題方法

積分項是卷積:

∫0ty(τ)cos⁡(t−τ) dτ=(y∗cos⁡t)(t).\int_0^t y(\tau)\cos(t-\tau)\,d\tau =(y*\cos t)(t).

因此拉普拉斯轉換為

Y(s)+Y(s)ss2+1=L{tu(t−1)}+1s+2.Y(s)+Y(s)\frac{s}{s^2+1} =\mathcal{L}\{tu(t-1)\}+\frac1{s+2}.

先處理右側的平移項。因為

tu(t−1)=[(t−1)+1]u(t−1),tu(t-1)=[(t-1)+1]u(t-1),

所以

L{tu(t−1)}=e−s(1s2+1s).\mathcal{L}\{tu(t-1)\} =e^{-s}\left(\frac1{s^2}+\frac1s\right).

因此

Y(s)(1+ss2+1)=e−s(1s2+1s)+1s+2.Y(s)\left(1+\frac{s}{s^2+1}\right) =e^{-s}\left(\frac1{s^2}+\frac1s\right)+\frac1{s+2}.

整理左側:

1+ss2+1=s2+s+1s2+1.1+\frac{s}{s^2+1} =\frac{s^2+s+1}{s^2+1}.

故

Y(s)=e−s(s2+1)(s+1)s2(s2+s+1)+s2+1(s+2)(s2+s+1).Y(s) = e^{-s}\frac{(s^2+1)(s+1)}{s^2(s^2+s+1)} + \frac{s^2+1}{(s+2)(s^2+s+1)}.

第一部分反轉換

部分分式分解得

(s2+1)(s+1)s2(s2+s+1)=1s2+s−1s2+s+1.\frac{(s^2+1)(s+1)}{s^2(s^2+s+1)} = \frac1{s^2}+\frac{s-1}{s^2+s+1}.

注意

s2+s+1=(s+12)2+34.s^2+s+1 =\left(s+\frac12\right)^2+\frac34.

因此

L−1{1s2}=t,\mathcal{L}^{-1}\left\{\frac1{s^2}\right\}=t,

且

s−1s2+s+1=s+12(s+12)2+(32)2−32(s+12)2+(32)2.\frac{s-1}{s^2+s+1} = \frac{s+\frac12}{\left(s+\frac12\right)^2+\left(\frac{\sqrt3}{2}\right)^2} - \frac{\frac32}{\left(s+\frac12\right)^2+\left(\frac{\sqrt3}{2}\right)^2}.

反轉換後為

L−1{s−1s2+s+1}=e−t/2cos⁡3t2−3e−t/2sin⁡3t2.\mathcal{L}^{-1} \left\{ \frac{s-1}{s^2+s+1} \right\} = e^{-t/2}\cos\frac{\sqrt3t}{2} -\sqrt3e^{-t/2}\sin\frac{\sqrt3t}{2}.

令

🔒

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

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

免費註冊

第 4 題

Calculate the following integral of complex functions using Residue theorem
(1). ∫02π2π2+cos⁡θdθ\int_0^{2\pi} \frac{2\pi}{2+\cos\theta} d\theta (10%)
(2). pr.v.∫−∞∞1(1+x)(x−2)2dx\text{pr.v.} \int_{-\infty}^{\infty} \frac{1}{(1+x)(x-2)^2} dx (10%)

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

這一題的完整詳解

核心觀念

本題考查複變數積分中的兩項基本技巧:

  1. 利用 z=eiθz=e^{i\theta} 將三角積分轉換為單位圓上的複積分。
  2. 利用留數定理計算封閉路徑積分;若積分路徑上有極點,必須判斷主值是否存在。

留數定理為

∮Cf(z) dz=2πi∑Res⁡(f;zk),\oint_C f(z)\,dz=2\pi i\sum \operatorname{Res}(f;z_k),

其中 zkz_k 為路徑 CC 內部的孤立奇點。


(1)計算 ∫02π2π2+cos⁡θ dθ\displaystyle \int_0^{2\pi}\frac{2\pi}{2+\cos\theta}\,d\theta

解題方法

令

z=eiθ,dθ=dziz,z=e^{i\theta},\qquad d\theta=\frac{dz}{iz},

且

cos⁡θ=12(z+1z).\cos\theta=\frac{1}{2}\left(z+\frac{1}{z}\right).

當 θ\theta 從 00 變化至 2π2\pi 時,zz 沿單位圓 ∣z∣=1|z|=1 逆時針繞行一次。因此

I=∮∣z∣=12π2+12(z+1z)dziz=∮∣z∣=14πi(z2+4z+1) dz.\begin{aligned} I &=\oint_{|z|=1} \frac{2\pi}{2+\frac12\left(z+\frac1z\right)} \frac{dz}{iz} \\ &=\oint_{|z|=1} \frac{4\pi}{i(z^2+4z+1)}\,dz. \end{aligned}

因為

z2+4z+1=0z^2+4z+1=0

的兩根為

z=−2±3.z=-2\pm\sqrt3.

其中

z1=−2+3z_1=-2+\sqrt3

滿足 ∣z1∣<1|z_1|<1,而

z2=−2−3z_2=-2-\sqrt3

滿足 ∣z2∣>1|z_2|>1。所以單位圓內只有 z1z_1 一個極點。

令

f(z)=4πi(z2+4z+1).f(z)=\frac{4\pi}{i(z^2+4z+1)}.

在 z1=−2+3z_1=-2+\sqrt3 處的留數為

Res⁡(f;z1)=4πi(2z1+4)=4πi(23)=−2πi3.\begin{aligned} \operatorname{Res}(f;z_1) &=\frac{4\pi}{i(2z_1+4)} \\ &=\frac{4\pi}{i(2\sqrt3)} =-\frac{2\pi i}{\sqrt3}. \end{aligned}

由留數定理,

I=2πi(−2πi3)=4π23.\begin{aligned} I &=2\pi i\left(-\frac{2\pi i}{\sqrt3}\right) \\ &=\frac{4\pi^2}{\sqrt3}. \end{aligned}

(2)計算 p.v.⁡∫−∞∞1(1+x)(x−2)2 dx\displaystyle \operatorname{p.v.}\int_{-\infty}^{\infty}\frac{1}{(1+x)(x-2)^2}\,dx

核心判斷

被積函數

f(x)=1(x+1)(x−2)2f(x)=\frac{1}{(x+1)(x-2)^2}

在實軸上有兩個奇點:

  • x=−1x=-1:一階極點;
  • x=2x=2:二階極點。

一階極點可以用柯西主值處理,但二階極點附近的主要行為為

f(x)∼13(x−2)2.f(x)\sim \frac{1}{3(x-2)^2}.

而

🔒

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

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

免費註冊

其他考古題