112 年 國立清華大學工程與系統科學系碩士班丙組《工程數學》

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

第 1 題

  1. Solve the differential equations.
    (a) dydx+3y=9x\frac{dy}{dx} + 3y = 9x, obtain y(x)y(x) that subjects to y(0)=3y(0) = 3.
    (b) 2d2ydx2−4dydx+4y=exsec⁡x2\frac{d^2y}{dx^2} - 4\frac{dy}{dx} + 4y = e^x \sec x.
    (c) dydx=2y2+3x2xy\frac{dy}{dx} = \frac{2y^2+3x}{2xy}

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

這一題的完整詳解

核心觀念

本題分別考查:

  • 一階線性微分方程:使用積分因子法。
  • 二階常係數非齊次微分方程:先求齊次解,再以變數變換與變參數法求特解。
  • 可化為一階線性方程的非線性微分方程:令 v=y2v=y^2 降階。

(a)dydx+3y=9x,y(0)=3\displaystyle \frac{dy}{dx}+3y=9x,\quad y(0)=3

解題方法

這是一階線性微分方程

y′+P(x)y=Q(x),y'+P(x)y=Q(x),

其中

P(x)=3,Q(x)=9x.P(x)=3,\qquad Q(x)=9x.

積分因子為

μ(x)=e∫3 dx=e3x.\mu(x)=e^{\int 3\,dx}=e^{3x}.

方程兩邊乘上 e3xe^{3x}:

e3xy′+3e3xy=9xe3x.e^{3x}y'+3e^{3x}y=9xe^{3x}.

左側可寫成乘積微分:

ddx(e3xy)=9xe3x.\frac{d}{dx}\left(e^{3x}y\right)=9xe^{3x}.

積分得

e3xy=∫9xe3x dx+C.e^{3x}y=\int 9xe^{3x}\,dx+C.

利用分部積分:

∫9xe3x dx=3xe3x−e3x.\int 9xe^{3x}\,dx =3xe^{3x}-e^{3x}.

因此

e3xy=3xe3x−e3x+C,e^{3x}y=3xe^{3x}-e^{3x}+C,

所以通解為

y=3x−1+Ce−3x.y=3x-1+Ce^{-3x}.

代入初始條件 y(0)=3y(0)=3:

3=−1+C,3=-1+C,

得到

C=4.C=4.

因此

y(x)=3x−1+4e−3x.\boxed{y(x)=3x-1+4e^{-3x}}.

解題技巧與驗算

對特解可直接設為一次式 yp=ax+by_p=ax+b。代入:

a+3(ax+b)=9x.a+3(ax+b)=9x.

比較係數可得

a=3,b=−1.a=3,\qquad b=-1.

再加上齊次解 Ce−3xCe^{-3x},即可快速得到一般解。

驗算:

y′=3−12e−3x,y'=3-12e^{-3x},

因此

y′+3y=3−12e−3x+9x−3+12e−3x=9x.y'+3y =3-12e^{-3x}+9x-3+12e^{-3x} =9x.

且

y(0)=−1+4=3.y(0)=-1+4=3.

(b)2y′′−4y′+4y=exsec⁡x\displaystyle 2y''-4y'+4y=e^x\sec x

核心觀念

先除以 22,得

y′′−2y′+2y=12exsec⁡x.y''-2y'+2y=\frac12 e^x\sec x.

解法分為兩部分:

  1. 求齊次方程的通解。
  2. 求非齊次特解。

由於齊次方程的特徵根為 1±i1\pm i,含有 excos⁡xe^x\cos x、exsin⁡xe^x\sin x;右側也含有 exe^x,因此令

y=exuy=e^xu

可大幅簡化計算。

第一步:求齊次解

考慮齊次方程

y′′−2y′+2y=0.y''-2y'+2y=0.

其特徵方程為

r2−2r+2=0.r^2-2r+2=0.

解得

r=1±i.r=1\pm i.

因此齊次解為

yh=ex(C1cos⁡x+C2sin⁡x).y_h=e^x(C_1\cos x+C_2\sin x).

第二步:變數變換

令

y=exu.y=e^xu.

則

y′=ex(u′+u),y'=e^x(u'+u), y′′=ex(u′′+2u′+u).y''=e^x(u''+2u'+u).

代入左側:

\begin{align*}
y''-2y'+2y
&=e^x(u''+2u'+u)-2e^x(u'+u)+2e^xu\
&=e^x(u''+u).
\end{align*}

所以原方程化為

ex(u′′+u)=12exsec⁡x,e^x(u''+u)=\frac12 e^x\sec x,

即

u′′+u=12sec⁡x.u''+u=\frac12\sec x.

第三步:用變參數法求 uu 的特解

對齊次方程

u′′+u=0,u''+u=0,

取基本解

u1=cos⁡x,u2=sin⁡x.u_1=\cos x,\qquad u_2=\sin x.

其 Wronskian 為

W=∣cos⁡xsin⁡x−sin⁡xcos⁡x∣=cos⁡2x+sin⁡2x=1.W= \begin{vmatrix} \cos x&\sin x\\ -\sin x&\cos x \end{vmatrix} =\cos^2x+\sin^2x=1.

令

f(x)=12sec⁡x.f(x)=\frac12\sec x.

變參數法的特解公式為

up=−u1∫u2fW dx+u2∫u1fW dx.u_p=-u_1\int \frac{u_2f}{W}\,dx +u_2\int \frac{u_1f}{W}\,dx.

代入各項:

up=−cos⁡x∫12sin⁡xsec⁡x dx+sin⁡x∫12cos⁡xsec⁡x dx.u_p =-\cos x\int \frac12\sin x\sec x\,dx +\sin x\int \frac12\cos x\sec x\,dx.

因為

sin⁡xsec⁡x=tan⁡x,cos⁡xsec⁡x=1,\sin x\sec x=\tan x,\qquad \cos x\sec x=1,

故

up=−12cos⁡x∫tan⁡x dx+12sin⁡x∫1 dx.u_p =-\frac12\cos x\int\tan x\,dx +\frac12\sin x\int 1\,dx.

又

🔒

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

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

免費註冊

第 2 題

  1. Find the series solution of the following differential equation about x=0x = 0.
    9x2d2ydx2−9xdydx−4y=09x^2 \frac{d^2y}{dx^2} - 9x\frac{dy}{dx} - 4y = 0.
    You have to express the solution in the form of y(x)=C1y1(x)+C2y2(x)y(x) = C_1y_1(x) + C_2y_2(x). To save time, you can only show the first four terms of y1(x)y_1(x) and y2(x)y_2(x).

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

這一題的完整詳解

核心觀念

  1. 正則奇異點(Regular Singular Point)之判別:
    對於二階線性齊次常微分方程式 y′′+P(x)y′+Q(x)y=0y'' + P(x)y' + Q(x)y = 0,若在 x=x0x = x_0 處 P(x)P(x) 或 Q(x)Q(x) 不具解析性(Analytic),但極限 lim⁡x→x0(x−x0)P(x)=p0\lim_{x \to x_0} (x - x_0)P(x) = p_0 與 lim⁡x→x0(x−x0)2Q(x)=q0\lim_{x \to x_0} (x - x_0)^2 Q(x) = q_0 均存在且為有限值,則 x=x0x = x_0 稱為正則奇異點。
  2. 佛羅貝尼烏斯方法(Frobenius Method):
    在正則奇異點 x=0x = 0 展開的級數解形式為: y(x)=∑n=0∞anxn+r,a0≠0y(x) = \sum_{n=0}^{\infty} a_n x^{n+r}, \quad a_0 \neq 0
  3. 指標方程式(Indicial Equation):
    由級數最低次項係數為零所得之代數方程式 r(r−1)+p0r+q0=0r(r-1) + p_0 r + q_0 = 0。兩指標根 r1,r2r_1, r_2 若滿足 r1−r2∉Zr_1 - r_2 \notin \mathbb{Z}(差值非整數),則保證存在兩個線性獨立的 Frobenius 級數解 y1(x)y_1(x) 與 y2(x)y_2(x)。
  4. 柯西-尤拉方程式(Cauchy-Euler Equation)特質:
    本題為標準柯西-尤拉方程式,代入級數展開時無平移項,其遞迴關係式將顯示高階項係數全為零。

解題方法

步驟一:判別展開點 x=0x = 0 之奇異性

將原微分方程式化為標準形式:

y′′−1xy′−49x2y=0y'' - \frac{1}{x} y' - \frac{4}{9x^2} y = 0

其中係數函數為 P(x)=−1xP(x) = -\frac{1}{x}, Q(x)=−49x2Q(x) = -\frac{4}{9x^2}。
計算極限:

p0=lim⁡x→0xP(x)=−1p_0 = \lim_{x \to 0} x P(x) = -1 q0=lim⁡x→0x2Q(x)=−49q_0 = \lim_{x \to 0} x^2 Q(x) = -\frac{4}{9}

由於 p0p_0 與 q0q_0 均為有限常數,故 x=0x = 0 為方程式之正則奇異點,適用 Frobenius 級數解法。


步驟二:設定 Frobenius 級數解並求導

設解為:

y(x)=∑n=0∞anxn+r(a0≠0)y(x) = \sum_{n=0}^{\infty} a_n x^{n+r} \quad (a_0 \neq 0)

對 xx 逐項微分可得:

y′(x)=∑n=0∞(n+r)anxn+r−1y'(x) = \sum_{n=0}^{\infty} (n+r) a_n x^{n+r-1} y′′(x)=∑n=0∞(n+r)(n+r−1)anxn+r−2y''(x) = \sum_{n=0}^{\infty} (n+r)(n+r-1) a_n x^{n+r-2}

步驟三:代回原微分方程式整理同次項

將級數代入原方程式 9x2y′′−9xy′−4y=09x^2 y'' - 9x y' - 4y = 0:

9x2∑n=0∞(n+r)(n+r−1)anxn+r−2−9x∑n=0∞(n+r)anxn+r−1−4∑n=0∞anxn+r=09x^2 \sum_{n=0}^{\infty} (n+r)(n+r-1) a_n x^{n+r-2} - 9x \sum_{n=0}^{\infty} (n+r) a_n x^{n+r-1} - 4 \sum_{n=0}^{\infty} a_n x^{n+r} = 0

將 xx 的冪次項整併:

∑n=0∞[9(n+r)(n+r−1)−9(n+r)−4]anxn+r=0\sum_{n=0}^{\infty} \left[ 9(n+r)(n+r-1) - 9(n+r) - 4 \right] a_n x^{n+r} = 0

定義多項式函數 F(s)=9s(s−1)−9s−4=9s2−18s−4F(s) = 9s(s-1) - 9s - 4 = 9s^2 - 18s - 4,方程式可簡記為:

∑n=0∞F(n+r)anxn+r=0\sum_{n=0}^{\infty} F(n+r) a_n x^{n+r} = 0

步驟四:求解指標方程式與指標根

令最低次項(n=0n = 0)的係數為零:

F(r)a0=0F(r) a_0 = 0

因為 a0≠0a_0 \neq 0,故得指標方程式:

9r2−18r−4=09r^2 - 18r - 4 = 0

利用二次方程式公式求解:

🔒

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

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

免費註冊

第 3 題

  1. Use the Laplace transform to solve the problem and obtain y(t)y(t).
    d2ydt2+4y=4tu(t−2)\frac{d^2y}{dt^2} + 4y = 4t u(t - 2),
    where y(0)=y′(0)=0y(0) = y'(0) = 0 and u(t)u(t) is the unit step function.

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

這一題的完整詳解

核心觀念

本題考查:

  • 單位階梯函數 u(t−a)u(t-a) 的拉普拉斯轉換。
  • 第二位移定理:
    L−1{e−asF(s)}=u(t−a)f(t−a).\mathcal{L}^{-1}\{e^{-as}F(s)\}=u(t-a)f(t-a).
  • 微分方程的拉普拉斯轉換:
    L{y′′(t)}=s2Y(s)−sy(0)−y′(0).\mathcal{L}\{y''(t)\}=s^2Y(s)-sy(0)-y'(0).
  • 部分分式分解。

由於右側在 t=2t=2 才開始作用,使用拉普拉斯轉換能直接處理階梯函數。

解題方法

原方程為

y′′+4y=4t u(t−2),y''+4y=4t\,u(t-2),

且

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

1. 右側的拉普拉斯轉換

將 tt 改寫成 (t−2)+2(t-2)+2:

4t u(t−2)=4[(t−2)+2]u(t−2).4t\,u(t-2) =4\big[(t-2)+2\big]u(t-2).

令

g(t)=4(t+2),g(t)=4(t+2),

則

4t u(t−2)=u(t−2)g(t−2).4t\,u(t-2)=u(t-2)g(t-2).

因此由第二位移定理,

L{4t u(t−2)}=e−2sL{4(t+2)}.\mathcal{L}\{4t\,u(t-2)\} =e^{-2s}\mathcal{L}\{4(t+2)\}.

計算得

L{4(t+2)}=4(1s2+2s).\mathcal{L}\{4(t+2)\} =4\left(\frac{1}{s^2}+\frac{2}{s}\right).

所以

L{4t u(t−2)}=4e−2s(1s2+2s).\mathcal{L}\{4t\,u(t-2)\} =4e^{-2s}\left(\frac{1}{s^2}+\frac{2}{s}\right).

2. 對微分方程取拉普拉斯轉換

設

Y(s)=L{y(t)}.Y(s)=\mathcal{L}\{y(t)\}.

由初始條件 y(0)=y′(0)=0y(0)=y'(0)=0,

L{y′′}=s2Y(s).\mathcal{L}\{y''\}=s^2Y(s).

因此

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

整理:

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

即

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

3. 部分分式分解

分解

4+8ss2(s2+4)=As+Bs2+Cs+Ds2+4.\frac{4+8s}{s^2(s^2+4)} = \frac{A}{s}+\frac{B}{s^2} +\frac{Cs+D}{s^2+4}.

兩邊同乘 s2(s2+4)s^2(s^2+4):

4+8s=As(s2+4)+B(s2+4)+(Cs+D)s2.4+8s =As(s^2+4)+B(s^2+4)+(Cs+D)s^2.

比較各次方係數:

A=2,B=1,C=−2,D=−1.A=2,\qquad B=1,\qquad C=-2,\qquad D=-1.

故

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

因此

🔒

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

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

免費註冊

第 4 題

  1. Consider the matrix
    M=[233−2−1−2201]M = \begin{bmatrix} 2 & 3 & 3 \\ -2 & -1 & -2 \\ 2 & 0 & 1 \end{bmatrix}
    (a) Find the determinant of MM and obtain the inverse matrix M−1M^{-1}.
    (b) Estimate the eigenvalues and eigenvectors of MM.

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

這一題的完整詳解

核心觀念

本題考查:

  • 3×33\times 3 矩陣的行列式計算。

  • 伴隨矩陣法求反矩陣:

    M−1=1det⁡Madj⁡(M)M^{-1}=\frac{1}{\det M}\operatorname{adj}(M)

  • 特徵值與特徵向量的定義:

    det⁡(M−λI)=0\det(M-\lambda I)=0

    求得特徵值 λ\lambda 後,再由

    (M−λI)v=0(M-\lambda I)\mathbf v=\mathbf 0

    求對應的特徵向量 v\mathbf v。


(a)行列式與反矩陣

矩陣為

M=[233−2−1−2201]M= \begin{bmatrix} 2&3&3\\ -2&-1&-2\\ 2&0&1 \end{bmatrix}

沿第一列展開行列式:

det⁡M=2∣−1−201∣−3∣−2−221∣+3∣−2−120∣=2(−1)−3((−2)(1)−(−2)(2))+3((−2)(0)−(−1)(2))=−2−3(2)+3(2)=−2.\begin{aligned} \det M &=2 \begin{vmatrix} -1&-2\\ 0&1 \end{vmatrix} -3 \begin{vmatrix} -2&-2\\ 2&1 \end{vmatrix} +3 \begin{vmatrix} -2&-1\\ 2&0 \end{vmatrix}\\ &=2(-1)-3((-2)(1)-(-2)(2)) +3((-2)(0)-(-1)(2))\\ &=-2-3(2)+3(2)\\ &=-2. \end{aligned}

因此

det⁡M=−2≠0\boxed{\det M=-2\neq 0}

所以 MM 可逆。

餘因子矩陣

各元素的代數餘因子為

C=[−1−22−3−46−3−24]C= \begin{bmatrix} -1&-2&2\\ -3&-4&6\\ -3&-2&4 \end{bmatrix}

因此伴隨矩陣為餘因子矩陣的轉置:

adj⁡(M)=CT=[−1−3−3−2−4−2264]\operatorname{adj}(M)=C^T = \begin{bmatrix} -1&-3&-3\\ -2&-4&-2\\ 2&6&4 \end{bmatrix}

故

M−1=1−2[−1−3−3−2−4−2264]=[123232121−1−3−2].\begin{aligned} M^{-1} &=\frac{1}{-2} \begin{bmatrix} -1&-3&-3\\ -2&-4&-2\\ 2&6&4 \end{bmatrix}\\ &= \begin{bmatrix} \frac12&\frac32&\frac32\\ 1&2&1\\ -1&-3&-2 \end{bmatrix}. \end{aligned}

所以

M−1=[123232121−1−3−2]\boxed{ M^{-1}= \begin{bmatrix} \frac12&\frac32&\frac32\\ 1&2&1\\ -1&-3&-2 \end{bmatrix}}

(b)特徵值與特徵向量

特徵方程為

det⁡(λI−M)=0\det(\lambda I-M)=0

其中

λI−M=[λ−2−3−32λ+12−20λ−1].\lambda I-M= \begin{bmatrix} \lambda-2&-3&-3\\ 2&\lambda+1&2\\ -2&0&\lambda-1 \end{bmatrix}.

計算可得

det⁡(λI−M)=λ3−2λ2−λ+2.\det(\lambda I-M) =\lambda^3-2\lambda^2-\lambda+2.

因式分解:

λ3−2λ2−λ+2=(λ−1)(λ2−λ−2)=(λ−1)(λ−2)(λ+1).\lambda^3-2\lambda^2-\lambda+2 =(\lambda-1)(\lambda^2-\lambda-2) =(\lambda-1)(\lambda-2)(\lambda+1).

因此三個特徵值為

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

對應 λ=1\lambda=1

解

(M−I)v=0(M-I)\mathbf v=\mathbf 0

即

🔒

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

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

免費註冊

第 5 題

  1. (a) Use the divergence theorem to determine the surface integral of the vector function F=xyzi+(yz+3)j+(z−2)kF = xyz\mathbf{i} + (yz + 3)\mathbf{j} + (z-2)\mathbf{k} on the surface SS of the region 0≤z≤50 \le z \le 5, x2+y2≤9x^2 + y^2 \le 9.
    (b) Use Stokes' theorem to determine the line integral of the vector function F=x2ln⁡xi+ysin⁡(y)j+tan⁡(z2)kF = x^2 \ln x \mathbf{i} + y\sin(y)\mathbf{j} + \tan(z^2)\mathbf{k} along a close curve CC: r(θ)=cos⁡θi+sin⁡θj+(cos⁡(2θ)+1)k\mathbf{r}(\theta) = \cos\theta \mathbf{i} + \sin\theta \mathbf{j} + (\cos(2\theta) + 1)\mathbf{k}, 0≤θ≤2π0 \le \theta \le 2\pi.

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

這一題的完整詳解

第 5 題

核心觀念

本題分別考查:

  • 散度定理(Divergence Theorem)

    對封閉曲面 SS 及其所包圍的體積區域 VV,

∬SF⋅n dS=∭V∇⋅F dV.\iint_S \mathbf F\cdot \mathbf n\,dS =\iiint_V \nabla\cdot\mathbf F\,dV.
  • Stokes 定理(Stokes’ Theorem)

    對閉曲線 CC 及其所張成的曲面 Σ\Sigma,

∮CF⋅dr=∬Σ(∇×F)⋅n dS.\oint_C \mathbf F\cdot d\mathbf r =\iint_\Sigma (\nabla\times\mathbf F)\cdot\mathbf n\,dS.

(a) 使用散度定理

給定

F=xyz i+(yz+3) j+(z−2) k.\mathbf F=xyz\,\mathbf i+(yz+3)\,\mathbf j+(z-2)\,\mathbf k.

區域為

V={(x,y,z):0≤z≤5, x2+y2≤9},V=\{(x,y,z):0\le z\le 5,\ x^2+y^2\le 9\},

即半徑 33、高度 55 的圓柱體,且 SS 為其封閉表面。

計算散度

∇⋅F=∂∂x(xyz)+∂∂y(yz+3)+∂∂z(z−2).\nabla\cdot\mathbf F =\frac{\partial}{\partial x}(xyz) +\frac{\partial}{\partial y}(yz+3) +\frac{\partial}{\partial z}(z-2).

因此

∇⋅F=yz+z+1.\nabla\cdot\mathbf F=yz+z+1.

由散度定理,

∬SF⋅n dS=∭V(yz+z+1) dV.\iint_S\mathbf F\cdot\mathbf n\,dS =\iiint_V(yz+z+1)\,dV.

改用柱座標

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

其中

0≤r≤3,0≤θ≤2π,0≤z≤5,0\le r\le 3,\qquad 0\le\theta\le 2\pi,\qquad 0\le z\le 5,

且 dV=r dz dr dθdV=r\,dz\,dr\,d\theta。

故

∭Vyz dV=0,\iiint_V yz\,dV=0,

因為圓盤 x2+y2≤9x^2+y^2\le 9 關於 y=0y=0 對稱,而被積函數 yzyz 對 yy 為奇函數。

接著,

∭Vz dV=(∬x2+y2≤9dA)(∫05z dz).\iiint_V z\,dV =\left(\iint_{x^2+y^2\le 9}dA\right) \left(\int_0^5z\,dz\right).

圓盤面積為 9π9\pi,所以

∭Vz dV=9π[z22]05=9π⋅252=225π2.\iiint_V z\,dV =9\pi\left[\frac{z^2}{2}\right]_0^5 =9\pi\cdot\frac{25}{2} =\frac{225\pi}{2}.

另外,

∭V1 dV=圓柱體體積=π(3)2(5)=45π.\iiint_V1\,dV =\text{圓柱體體積} =\pi(3)^2(5)=45\pi.

因此

∬SF⋅n dS=0+225π2+45π=315π2.\iint_S\mathbf F\cdot\mathbf n\,dS =0+\frac{225\pi}{2}+45\pi =\frac{315\pi}{2}.

(b) 使用 Stokes 定理

給定

🔒

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

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

免費註冊

第 6 題

  1. Solve the Laplace equation ∂2u∂x2+∂2u∂y2+∂2u∂z2=0\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} = 0 in a rectangular parallelepiped shown in the figure with the right side (y=by = b) kept at temperature g(x,z)g(x,z) and the left side (y=0y = 0) kept at f(x,z)f(x,z). The remaining sides are kept at temperature zero.
    🖼️【此處有附圖,請對照原卷】
🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 3 頁

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

這一題的完整詳解

核心觀念

本題求解長方體區域內的三維 Laplace 方程:

uxx+uyy+uzz=0.u_{xx}+u_{yy}+u_{zz}=0.

由圖中座標可視為

0<x<a,0<y<b,0<z<c.0<x<a,\qquad 0<y<b,\qquad 0<z<c.

邊界條件為

u(x,0,z)=f(x,z),u(x,b,z)=g(x,z),u(x,0,z)=f(x,z),\qquad u(x,b,z)=g(x,z),

而其餘四個側面溫度為零:

u(0,y,z)=u(a,y,z)=0,u(0,y,z)=u(a,y,z)=0, u(x,y,0)=u(x,y,c)=0.u(x,y,0)=u(x,y,c)=0.

由於 x,zx,z 方向的邊界值為零,適合使用雙重 Fourier 正弦級數。


解題方法:分離變數

令

u(x,y,z)=X(x)Y(y)Z(z).u(x,y,z)=X(x)Y(y)Z(z).

代入 Laplace 方程得

X′′X+Y′′Y+Z′′Z=0.\frac{X''}{X}+\frac{Y''}{Y}+\frac{Z''}{Z}=0.

針對 xx 與 zz 方向的零邊界條件,分別有

X(0)=X(a)=0,X(0)=X(a)=0, Z(0)=Z(c)=0.Z(0)=Z(c)=0.

因此特徵函數為

Xm(x)=sin⁡mπxa,m=1,2,…,X_m(x)=\sin\frac{m\pi x}{a},\qquad m=1,2,\ldots, Zn(z)=sin⁡nπzc,n=1,2,….Z_n(z)=\sin\frac{n\pi z}{c},\qquad n=1,2,\ldots.

相對應的特徵值為

λmn2=(mπa)2+(nπc)2.\lambda_{mn}^2 = \left(\frac{m\pi}{a}\right)^2 + \left(\frac{n\pi}{c}\right)^2.

故

λmn=(mπa)2+(nπc)2.\lambda_{mn} = \sqrt{ \left(\frac{m\pi}{a}\right)^2 + \left(\frac{n\pi}{c}\right)^2 }.

此時 yy 方向滿足

Y′′−λmn2Y=0,Y''-\lambda_{mn}^2Y=0,

其通解為

Y(y)=Asinh⁡(λmny)+Bcosh⁡(λmny).Y(y)=A\sinh(\lambda_{mn}y)+B\cosh(\lambda_{mn}y).

因此整體解可寫成

u(x,y,z)=∑m=1∞∑n=1∞Ymn(y)sin⁡mπxasin⁡nπzc.u(x,y,z) = \sum_{m=1}^{\infty} \sum_{n=1}^{\infty} Y_{mn}(y) \sin\frac{m\pi x}{a} \sin\frac{n\pi z}{c}.

展開兩個非齊次邊界條件

將 f(x,z)f(x,z) 與 g(x,z)g(x,z) 展開成雙重正弦級數:

f(x,z)=∑m=1∞∑n=1∞Fmnsin⁡mπxasin⁡nπzc,f(x,z) = \sum_{m=1}^{\infty} \sum_{n=1}^{\infty} F_{mn} \sin\frac{m\pi x}{a} \sin\frac{n\pi z}{c}, g(x,z)=∑m=1∞∑n=1∞Gmnsin⁡mπxasin⁡nπzc.g(x,z) = \sum_{m=1}^{\infty} \sum_{n=1}^{\infty} G_{mn} \sin\frac{m\pi x}{a} \sin\frac{n\pi z}{c}.

由正交性,係數為

🔒

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

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

免費註冊

第 7 題

  1. Apply Laplace transform to solve the displacement u(x,t)u(x,t) of a string which is driven by an external force:
    ∂2u∂t2=c2∂2u∂x2+sin⁡(πx)sin⁡(ωt),0<x<1,t>0\frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2} + \sin(\pi x) \sin(\omega t), \quad 0 < x < 1, t > 0
    u(0,t)=0,u(1,t)=0,t>0u(0, t) = 0, \quad u(1, t) = 0, \quad t > 0
    u(x,0)=0,∂u∂t(x,0)=0,0<x<1u(x, 0) = 0, \quad \frac{\partial u}{\partial t}(x, 0) = 0, \quad 0 < x < 1

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

這一題的完整詳解

核心觀念

本題考查:

  1. 對時間變數 tt 作 Laplace transform。
  2. 利用初始條件處理時間微分:
L{∂2u∂t2}=s2U(x,s)−su(x,0)−ut(x,0).\mathcal{L}\left\{\frac{\partial^2u}{\partial t^2}\right\} =s^2U(x,s)-su(x,0)-u_t(x,0).
  1. 辨認空間外力 sin⁡(πx)\sin(\pi x) 正好是固定端弦的特徵函數,因此可令
    u(x,t)=a(t)sin⁡(πx).u(x,t)=a(t)\sin(\pi x).
  2. 注意共振情形 ω=cπ\omega=c\pi,此時一般公式的分母為零,必須單獨處理。

令

α=cπ.\alpha=c\pi.

解題方法一:先作空間模態分離

由於外力項為

sin⁡(πx)sin⁡(ωt),\sin(\pi x)\sin(\omega t),

且邊界條件要求 u(0,t)=u(1,t)=0u(0,t)=u(1,t)=0,可設

u(x,t)=a(t)sin⁡(πx).u(x,t)=a(t)\sin(\pi x).

此時

utt=a′′(t)sin⁡(πx),u_{tt}=a''(t)\sin(\pi x),

以及

uxx=−π2a(t)sin⁡(πx).u_{xx}=-\pi^2a(t)\sin(\pi x).

代回原偏微分方程:

a′′(t)sin⁡(πx)=−c2π2a(t)sin⁡(πx)+sin⁡(πx)sin⁡(ωt).a''(t)\sin(\pi x) = -c^2\pi^2a(t)\sin(\pi x) +\sin(\pi x)\sin(\omega t).

消去共同因子 sin⁡(πx)\sin(\pi x),得到常微分方程

a′′(t)+c2π2a(t)=sin⁡(ωt),a''(t)+c^2\pi^2a(t)=\sin(\omega t),

即

a′′(t)+α2a(t)=sin⁡(ωt).a''(t)+\alpha^2a(t)=\sin(\omega t).

初始條件變為

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

對時間作 Laplace transform

設

A(s)=L{a(t)}.A(s)=\mathcal{L}\{a(t)\}.

由初始條件 a(0)=a′(0)=0a(0)=a'(0)=0,

L{a′′(t)}=s2A(s).\mathcal{L}\{a''(t)\}=s^2A(s).

因此

s2A(s)+α2A(s)=L{sin⁡(ωt)}.s^2A(s)+\alpha^2A(s) = \mathcal{L}\{\sin(\omega t)\}.

又

L{sin⁡(ωt)}=ωs2+ω2,\mathcal{L}\{\sin(\omega t)\} = \frac{\omega}{s^2+\omega^2},

故

A(s)=ω(s2+ω2)(s2+α2).A(s) = \frac{\omega} {(s^2+\omega^2)(s^2+\alpha^2)}.

非共振情形:ω≠α\omega\ne\alpha

進行部分分式分解:

ω(s2+ω2)(s2+α2)=ωα2−ω2(1s2+ω2−1s2+α2).\frac{\omega} {(s^2+\omega^2)(s^2+\alpha^2)} = \frac{\omega}{\alpha^2-\omega^2} \left( \frac{1}{s^2+\omega^2} - \frac{1}{s^2+\alpha^2} \right).

利用

L−1{1s2+k2}=1ksin⁡(kt),\mathcal{L}^{-1} \left\{\frac{1}{s^2+k^2}\right\} = \frac{1}{k}\sin(kt),

得到

a(t)=ωα2−ω2(sin⁡(ωt)ω−sin⁡(αt)α).a(t) = \frac{\omega}{\alpha^2-\omega^2} \left( \frac{\sin(\omega t)}{\omega} - \frac{\sin(\alpha t)}{\alpha} \right).

整理為

🔒

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

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

免費註冊

第 8 題

  1. (a) Find the Laurent series of the complex function f(z)=1z(z+4i)f(z) = \frac{1}{z(z+4i)} for ∣z∣>3|z| > 3.
    (b) Determine the integral ∮C(z2exp⁡(1z)+sin⁡zz5+3z3)dz\oint_C (z^2 \exp(\frac{1}{z}) + \frac{\sin z}{z^5+3z^3}) dz where CC is the circle of ∣z−3i2∣=2|z - \frac{3i}{2}| = 2, oriented counterclockwise.

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

這一題的完整詳解

第 8 題

核心觀念

本題考查:

  1. Laurent 級數的收斂環域與幾何級數展開。
  2. 留數定理:
    ∮CF(z) dz=2πi∑Res⁡(F;曲線內奇點)\oint_C F(z)\,dz=2\pi i\sum \operatorname{Res}(F;\text{曲線內奇點})
  3. 本性奇點的留數計算。
  4. 有理函數乘上解析函數後,在單極點與高階極點的留數求法。
  5. 判斷奇點是否位於圓形路徑內。

(a) Laurent 級數

給定

f(z)=1z(z+4i).f(z)=\frac{1}{z(z+4i)}.

其奇點為

z=0,z=−4i.z=0,\qquad z=-4i.

相對於 Laurent 展開中心 z=0z=0,兩個奇點的距離分別為 00 與 44。因此,題目寫成 ∣z∣>3|z|>3 時,區域中包含了奇點 z=−4iz=-4i 的圓周位置,實際上不能用同一個 Laurent 級數涵蓋整個 ∣z∣>3|z|>3。

所以必須分成兩個收斂環域:

3<∣z∣<43<|z|<4

以及

∣z∣>4.|z|>4.

情形一:3<∣z∣<43<|z|<4

將函數部分分式分解:

1z(z+4i)=14i(1z−1z+4i).\frac{1}{z(z+4i)} = \frac{1}{4i}\left(\frac{1}{z}-\frac{1}{z+4i}\right).

因為 ∣z∣<4|z|<4,可將第二項寫成

1z+4i=14i11+z4i.\frac{1}{z+4i} = \frac{1}{4i}\frac{1}{1+\frac{z}{4i}}.

利用幾何級數公式

11+w=∑n=0∞(−w)n,∣w∣<1,\frac{1}{1+w}=\sum_{n=0}^{\infty}(-w)^n,\qquad |w|<1,

得到

1z+4i=14i∑n=0∞(−z4i)n.\frac{1}{z+4i} = \frac{1}{4i} \sum_{n=0}^{\infty} \left(-\frac{z}{4i}\right)^n.

因此

f(z)=14iz+116∑n=0∞(−z4i)n,3<∣z∣<4.f(z) = \frac{1}{4iz} + \frac{1}{16} \sum_{n=0}^{\infty} \left(-\frac{z}{4i}\right)^n, \qquad 3<|z|<4.

前幾項為

f(z)=14iz+116+iz64−z2256−iz31024+⋯ .f(z) = \frac{1}{4iz} +\frac{1}{16} +\frac{iz}{64} -\frac{z^2}{256} -\frac{iz^3}{1024} +\cdots.

情形二:∣z∣>4|z|>4

將原式改寫為

f(z)=1z211+4iz.f(z) = \frac{1}{z^2}\frac{1}{1+\frac{4i}{z}}.

因為 ∣z∣>4|z|>4,所以 ∣4iz∣<1\left|\frac{4i}{z}\right|<1,可使用幾何級數:

11+4iz=∑n=0∞(−4iz)n.\frac{1}{1+\frac{4i}{z}} = \sum_{n=0}^{\infty} \left(-\frac{4i}{z}\right)^n.

因此

f(z)=1z2∑n=0∞(−4iz)n=∑n=0∞(−4i)nz−n−2,∣z∣>4.f(z) = \frac{1}{z^2} \sum_{n=0}^{\infty} \left(-\frac{4i}{z}\right)^n = \sum_{n=0}^{\infty} (-4i)^n z^{-n-2}, \qquad |z|>4.

前幾項為

f(z)=1z2−4iz3−16z4+64iz5+256z6+⋯ .f(z) = \frac{1}{z^2} -\frac{4i}{z^3} -\frac{16}{z^4} +\frac{64i}{z^5} +\frac{256}{z^6} +\cdots.

(b) 複變線積分

要求

∮C(z2e1/z+sin⁡zz5+3z3) dz,\oint_C \left( z^2e^{1/z} + \frac{\sin z}{z^5+3z^3} \right)\,dz,

其中 CC 為

∣z−3i2∣=2\left|z-\frac{3i}{2}\right|=2

且為逆時針方向。

令

F(z)=z2e1/z+sin⁡zz3(z2+3).F(z)=z^2e^{1/z}+\frac{\sin z}{z^3(z^2+3)}.

第一步:判斷曲線內的奇點

奇點包括:

  • z=0z=0;
  • z=i3z=i\sqrt{3};
  • z=−i3z=-i\sqrt{3}。

計算各點到圓心 3i2\frac{3i}{2} 的距離:

∣0−3i2∣=32<2,\left|0-\frac{3i}{2}\right|=\frac32<2,

所以 z=0z=0 在曲線內。

∣i3−3i2∣=3−32<2,\left|i\sqrt3-\frac{3i}{2}\right| = \sqrt3-\frac32<2,

所以 z=i3z=i\sqrt3 在曲線內。

∣−i3−3i2∣=3+32>2,\left|-i\sqrt3-\frac{3i}{2}\right| = \sqrt3+\frac32>2,

所以 z=−i3z=-i\sqrt3 在曲線外。

🔒

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

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

免費註冊

其他考古題