109 年 國立成功大學電腦與通信工程研究所丙組《工程數學》

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

第 1 題20 分

  1. Solve the first – order differential equations (20%)
    (a) (2x−4y+5)y′+x−2y+3=0(2x - 4y + 5)y' + x - 2y + 3 = 0
    (b) xy′+y+4=0xy' + y + 4 = 0

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

這一題的完整詳解

核心觀念

本題考查一階微分方程的變數代換與可分離變數法:

  • 若方程中反覆出現線性組合 x−2yx-2y,可令
    u=x−2y.u=x-2y.
  • 若左側含有
    xy′+y,xy'+y,
    可利用乘積微分公式
    (xy)′=xy′+y.(xy)'=xy'+y.

(a)(2x−4y+5)y′+x−2y+3=0(2x-4y+5)y'+x-2y+3=0

解題方法

觀察
2x−4y+5=2(x−2y)+5,2x-4y+5=2(x-2y)+5,
因此令
u=x−2y.u=x-2y.

對 xx 微分得
dudx=1−2y′,\frac{du}{dx}=1-2y',
所以
y′=1−u′2.y'=\frac{1-u'}{2}.

原方程可改寫為
(2u+5)y′+u+3=0.(2u+5)y'+u+3=0.

代入 y′=1−u′2y'=\dfrac{1-u'}{2}:

(2u+5)1−u′2+u+3=0.(2u+5)\frac{1-u'}{2}+u+3=0.

兩邊乘以 22:

(2u+5)(1−u′)+2u+6=0.(2u+5)(1-u')+2u+6=0.

整理得

2u+5−(2u+5)u′+2u+6=0,2u+5-(2u+5)u'+2u+6=0, (2u+5)u′=4u+11.(2u+5)u'=4u+11.

因此

dxdu=2u+54u+11.\frac{dx}{du}=\frac{2u+5}{4u+11}.

將分式拆解:

2u+54u+11=12−12(4u+11).\frac{2u+5}{4u+11} =\frac12-\frac{1}{2(4u+11)}.

積分得到

x=∫(12−12(4u+11)) du,x=\int\left(\frac12-\frac{1}{2(4u+11)}\right)\,du, x=u2−18ln⁡∣4u+11∣+C.x=\frac{u}{2}-\frac18\ln|4u+11|+C.

代回 u=x−2yu=x-2y,得一般解

x−2y2−18ln⁡∣4x−8y+11∣=x+C\boxed{ \frac{x-2y}{2} -\frac18\ln|4x-8y+11| =x+C }

也可整理為

x−2y2−x−18ln⁡∣4x−8y+11∣=C.\boxed{ \frac{x-2y}{2}-x -\frac18\ln|4x-8y+11| =C }.

特殊解

在推導過程中曾除以 4u+114u+11,因此需檢查

4u+11=0.4u+11=0.

若

u=−114,u=-\frac{11}{4},

則

x−2y=−114,x-2y=-\frac{11}{4},

即

y=x2+118.\boxed{ y=\frac{x}{2}+\frac{11}{8} }.

代入原方程:

2x−4y+5=−12,x−2y+3=14,2x-4y+5=-\frac12,\qquad x-2y+3=\frac14,

且 y′=12y'=\dfrac12,所以

🔒

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

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

免費註冊

第 2 題15 分

  1. Use the calculations of residues, calculate the following integral P=∫02πdθ(1−bcos⁡θ)2P=\int_{0}^{2\pi} \frac{d\theta}{(1-b \cos\theta)^2} where b<1b <1 (15%)

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

這一題的完整詳解

核心觀念

本題考查複變函數中的留數定理,並利用單位圓上的參數代換

z=eiθ,dθ=dziz,cos⁡θ=12(z+1z).z=e^{i\theta},\qquad d\theta=\frac{dz}{iz},\qquad \cos\theta=\frac12\left(z+\frac1z\right).

當 ∣b∣<1|b|<1 時,原積分可轉換為沿單位圓 ∣z∣=1|z|=1 的複數線積分,再計算圓內極點的留數:

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

題目寫作 b<1b<1,通常隱含 bb 為非負實數,因此以下以 0≤b<10\le b<1 為條件;一般而言,積分有限的充分條件為 ∣b∣<1|b|<1。

解題方法

令

P=∫02πdθ(1−bcos⁡θ)2.P=\int_0^{2\pi}\frac{d\theta}{(1-b\cos\theta)^2}.

先處理 b=0b=0:

P=∫02π1 dθ=2π.P=\int_0^{2\pi}1\,d\theta=2\pi.

以下考慮 0<b<10<b<1。取

z=eiθ,z=e^{i\theta},

則當 θ\theta 從 00 變化到 2π2\pi 時,zz 逆時針繞單位圓 C:∣z∣=1C:|z|=1 一周。由

cos⁡θ=12(z+1z),dθ=dziz,\cos\theta=\frac12\left(z+\frac1z\right),\qquad d\theta=\frac{dz}{iz},

可得

1−bcos⁡θ=1−b2(z+1z)=2z−b(z2+1)2z.1-b\cos\theta =1-\frac b2\left(z+\frac1z\right) =\frac{2z-b(z^2+1)}{2z}.

因此

dθ(1−bcos⁡θ)2=4z2(2z−b(z2+1))2⋅dziz.\frac{d\theta}{(1-b\cos\theta)^2} = \frac{4z^2}{\left(2z-b(z^2+1)\right)^2}\cdot\frac{dz}{iz}.

整理為

P=∮C4zi(2z−b(z2+1))2 dz.P=\oint_C \frac{4z}{i\left(2z-b(z^2+1)\right)^2}\,dz.

將二次式因式分解。令

s=1−b2,s=\sqrt{1-b^2},

則

2z−b(z2+1)=−b(z−1−sb)(z−1+sb).2z-b(z^2+1) =-b\left(z-\frac{1-s}{b}\right) \left(z-\frac{1+s}{b}\right).

定義

α=1−sb,β=1+sb.\alpha=\frac{1-s}{b},\qquad \beta=\frac{1+s}{b}.

由於

αβ=1,\alpha\beta=1,

且 0<b<10<b<1,所以

0<α<1,β>1.0<\alpha<1,\qquad \beta>1.

因此只有 z=αz=\alpha 位於單位圓內,而 z=βz=\beta 位於單位圓外。原積分化為

P=∮C4zib2(z−α)2(z−β)2 dz.P=\oint_C \frac{4z}{ib^2(z-\alpha)^2(z-\beta)^2}\,dz.

點 z=αz=\alpha 是二階極點。令

🔒

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

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

免費註冊

第 3 題15 分

  1. Let f(x)={0−π≤x<0x20≤x<πf(x) = \begin{cases} 0 & -\pi \leq x < 0 \\ x^2 & 0 \leq x < \pi \end{cases}
    Find the Fourier series of f(x)f(x) on the given interval (15%)

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

這一題的完整詳解

核心觀念

題目要求在區間 [−π,π)[-\pi,\pi) 上求函數的 Fourier series,採用週期 2π2\pi 的傅立葉級數:

f(x)∼a02+∑n=1∞(ancos⁡nx+bnsin⁡nx),f(x)\sim \frac{a_0}{2}+\sum_{n=1}^{\infty}\left(a_n\cos nx+b_n\sin nx\right),

其中

a0=1π∫−ππf(x) dx,a_0=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\,dx, an=1π∫−ππf(x)cos⁡nx dx,bn=1π∫−ππf(x)sin⁡nx dx.a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\cos nx\,dx,\qquad b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin nx\,dx.

由於

f(x)={0,−π≤x<0,x2,0≤x<π,f(x)= \begin{cases} 0,&-\pi\le x<0,\\ x^2,&0\le x<\pi, \end{cases}

所有積分皆只需計算 [0,π][0,\pi] 上的部分。


解題方法

1. 計算常數項 a0a_0

a0=1π∫0πx2 dx=1π[x33]0π=π23.a_0=\frac{1}{\pi}\int_0^\pi x^2\,dx =\frac{1}{\pi}\left[\frac{x^3}{3}\right]_0^\pi =\frac{\pi^2}{3}.

因此

a02=π26.\frac{a_0}{2}=\frac{\pi^2}{6}.

2. 計算餘弦係數 ana_n

an=1π∫0πx2cos⁡nx dx.a_n=\frac{1}{\pi}\int_0^\pi x^2\cos nx\,dx.

先作分部積分:

∫x2cos⁡nx dx=x2sin⁡nxn−2n∫xsin⁡nx dx.\int x^2\cos nx\,dx =\frac{x^2\sin nx}{n}-\frac{2}{n}\int x\sin nx\,dx.

又

∫xsin⁡nx dx=−xcos⁡nxn+sin⁡nxn2.\int x\sin nx\,dx =-\frac{x\cos nx}{n}+\frac{\sin nx}{n^2}.

所以

∫x2cos⁡nx dx=x2sin⁡nxn+2xcos⁡nxn2−2sin⁡nxn3.\int x^2\cos nx\,dx = \frac{x^2\sin nx}{n} +\frac{2x\cos nx}{n^2} -\frac{2\sin nx}{n^3}.

代入 00 與 π\pi:

∫0πx2cos⁡nx dx=2πcos⁡(nπ)n2=2π(−1)nn2.\int_0^\pi x^2\cos nx\,dx = \frac{2\pi\cos(n\pi)}{n^2} = \frac{2\pi(-1)^n}{n^2}.

因此

an=1π⋅2π(−1)nn2=2(−1)nn2.a_n=\frac{1}{\pi}\cdot\frac{2\pi(-1)^n}{n^2} =\frac{2(-1)^n}{n^2}.

3. 計算正弦係數 bnb_n

bn=1π∫0πx2sin⁡nx dx.b_n=\frac{1}{\pi}\int_0^\pi x^2\sin nx\,dx.

分部積分得

∫x2sin⁡nx dx=−x2cos⁡nxn+2n∫xcos⁡nx dx.\int x^2\sin nx\,dx =-\frac{x^2\cos nx}{n} +\frac{2}{n}\int x\cos nx\,dx.

而

∫xcos⁡nx dx=xsin⁡nxn+cos⁡nxn2.\int x\cos nx\,dx =\frac{x\sin nx}{n}+\frac{\cos nx}{n^2}.

因此

∫x2sin⁡nx dx=−x2cos⁡nxn+2xsin⁡nxn2+2cos⁡nxn3.\int x^2\sin nx\,dx = -\frac{x^2\cos nx}{n} +\frac{2x\sin nx}{n^2} +\frac{2\cos nx}{n^3}.

代入 00 與 π\pi:

🔒

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

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

免費註冊

第 4 題15 分

  1. Find the inverse Laplace transform of F(s)=1s(s+2)e−4sF(s) = \frac{1}{s(s+2)} e^{-4s}. (15%)

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

這一題的完整詳解

核心觀念

本題考查兩個拉普拉斯轉換觀念:

  1. 部分分式分解

1s(s+2)\frac{1}{s(s+2)}

  1. 第二位移定理
    若

L{f(t)}=F(s),\mathcal{L}\{f(t)\}=F(s),

則

L−1{e−asF(s)}=u(t−a)f(t−a),\mathcal{L}^{-1}\{e^{-as}F(s)\} =u(t-a)f(t-a),

其中 u(t−a)u(t-a) 為單位階梯函數。


解題方法

先處理未含延遲因子的部分:

1s(s+2)=As+Bs+2.\frac{1}{s(s+2)} =\frac{A}{s}+\frac{B}{s+2}.

因此

1=A(s+2)+Bs=(A+B)s+2A.1=A(s+2)+Bs =(A+B)s+2A.

比較係數得

A+B=0,2A=1,A+B=0,\qquad 2A=1,

所以

A=12,B=−12.A=\frac12,\qquad B=-\frac12.

故

1s(s+2)=12(1s−1s+2).\frac{1}{s(s+2)} =\frac12\left(\frac1s-\frac1{s+2}\right).

利用基本反拉普拉斯轉換公式

L−1{1s}=1,L−1{1s+2}=e−2t,\mathcal{L}^{-1}\left\{\frac1s\right\}=1, \qquad \mathcal{L}^{-1}\left\{\frac1{s+2}\right\}=e^{-2t},

可得

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

原式為

🔒

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

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

免費註冊

第 5 題20 分

  1. Solve the initial value problem below, which involves a Dirac delta function. (20%)
    y′′−4y′+13y=δ(t−3);y(0)=y′(0)=0y'' - 4y' + 13y = \delta(t-3); y(0) = y'(0) = 0

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

這一題的完整詳解

核心觀念

本題考查含 Dirac delta 函數的二階線性常微分方程,核心公式為:

  • 單位階躍函數 u(t−a)u(t-a) 的位移性質:
    L{u(t−a)f(t−a)}=e−asF(s)\mathcal{L}\{u(t-a)f(t-a)\}=e^{-as}F(s)
  • Dirac delta 函數的拉普拉斯變換:
    L{δ(t−a)}=e−as\mathcal{L}\{\delta(t-a)\}=e^{-as}
  • 基本反變換:
L−1{1(s−α)2+β2}=1βeαtsin⁡(βt)\mathcal{L}^{-1}\left\{\frac{1}{(s-\alpha)^2+\beta^2}\right\} =\frac{1}{\beta}e^{\alpha t}\sin(\beta t)

由於外力 δ(t−3)\delta(t-3) 只在 t=3t=3 瞬間作用,因此解在 t<3t<3 為零,t=3t=3 之後才產生響應。


解題方法:拉普拉斯變換

原方程為

y′′−4y′+13y=δ(t−3),y(0)=0,y′(0)=0.y''-4y'+13y=\delta(t-3), \qquad y(0)=0,\quad y'(0)=0.

設

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

利用初始條件,

L{y′′}=s2Y(s)−sy(0)−y′(0)=s2Y(s),\mathcal{L}\{y''\} =s^2Y(s)-sy(0)-y'(0) =s^2Y(s), L{y′}=sY(s)−y(0)=sY(s).\mathcal{L}\{y'\}=sY(s)-y(0)=sY(s).

因此左側拉普拉斯變換為

s2Y(s)−4sY(s)+13Y(s)=(s2−4s+13)Y(s).s^2Y(s)-4sY(s)+13Y(s) =(s^2-4s+13)Y(s).

右側為

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

所以

(s2−4s+13)Y(s)=e−3s,(s^2-4s+13)Y(s)=e^{-3s},

得到

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

將分母配方:

s2−4s+13=(s−2)2+9.s^2-4s+13=(s-2)^2+9.

故

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

先取未含位移因子的反拉普拉斯變換:

L−1{1(s−2)2+32}=13e2tsin⁡(3t).\mathcal{L}^{-1} \left\{\frac{1}{(s-2)^2+3^2}\right\} =\frac{1}{3}e^{2t}\sin(3t).

再利用第二位移定理:

y(t)=u(t−3)⋅13e2(t−3)sin⁡(3(t−3)).y(t) =u(t-3)\cdot \frac{1}{3}e^{2(t-3)} \sin\bigl(3(t-3)\bigr).

因此分段表示為

🔒

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

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

免費註冊

第 6 題15 分

  1. Evaluate ∫−∞∞1x6+64dx\int_{-\infty}^{\infty} \frac{1}{x^6 + 64} dx. (15%)

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

這一題的完整詳解

解題步驟

  1. 因子與極點

x6+64=0  ⟹  x6=−64=64eiπx^{6}+64=0\;\Longrightarrow\;x^{6}= -64 =64e^{i\pi}

第六根為

ak=2 eiπ+2kπ6,k=0,…,5a_k =2\,e^{i\frac{\pi+2k\pi}{6}},\qquad k=0,\dots ,5

上半平面的三個簡單極點

a1=2eiπ/6,  a2=2eiπ/2=2i,  a3=2ei5π/6a_1=2e^{i\pi/6},\; a_2=2e^{i\pi/2}=2i,\; a_3=2e^{i5\pi/6}

  1. 留數
    對 f(z)=1z6+64f(z)=\dfrac{1}{z^{6}+64} ,
Res⁡(f,z=a)=1(z6+64)′∣z=a=16a5.\operatorname{Res}\bigl(f,z=a\bigr)=\frac{1}{(z^{6}+64)'}\Big|_{z=a} =\frac{1}{6a^{5}} .

因 a=2eiθa=2e^{i\theta},a5=32ei5θa^{5}=32e^{i5\theta},故

Res⁡(f,a)=1192 e−i5θ.\operatorname{Res}(f,a)=\frac{1}{192}\,e^{-i5\theta}.

計算三個上半平面極點的留數和:

🔒

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

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

免費註冊

其他考古題