113 年 國立中正大學機械工程學系碩士班乙組《工程數學》

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

第 1 題5 分

  1. (5%) Please verify v⃗1\vec{v}_1、v⃗2\vec{v}_2、v⃗3\vec{v}_3 are linear independent, where v⃗1=[123]\vec{v}_1 = \begin{bmatrix} 1 \\ 2 \\ 3 \end{bmatrix}, v⃗2=[2−13]\vec{v}_2 = \begin{bmatrix} 2 \\ -1 \\ 3 \end{bmatrix}, v⃗3=[01−1]\vec{v}_3 = \begin{bmatrix} 0 \\ 1 \\ -1 \end{bmatrix}.

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

這一題的完整詳解

核心觀念

向量 v⃗1,v⃗2,v⃗3\vec v_1,\vec v_2,\vec v_3 線性獨立的定義是:若

c1v⃗1+c2v⃗2+c3v⃗3=0⃗,c_1\vec v_1+c_2\vec v_2+c_3\vec v_3=\vec 0,

則必須只有唯一解

c1=c2=c3=0.c_1=c_2=c_3=0.

將三個向量作為矩陣的欄向量:

A=[1202−1133−1].A= \begin{bmatrix} 1&2&0\\ 2&-1&1\\ 3&3&-1 \end{bmatrix}.

對於方陣而言,

det⁡(A)≠0\det(A)\neq 0

等價於矩陣可逆,也等價於其欄向量線性獨立。

解題方法

計算矩陣 AA 的行列式:

🔒

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

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

免費註冊

第 2 題25 分

  1. (25%) There is a simultaneous homogeneous differential equations X′=AXX' = AX, where X=[x1x2x3]X = \begin{bmatrix} x_1 \\ x_2 \\ x_3 \end{bmatrix} and A=[−1101−1000−2]A = \begin{bmatrix} -1 & 1 & 0 \\ 1 & -1 & 0 \\ 0 & 0 & -2 \end{bmatrix}.
    (1) (20%) Please find the eigenvalues and eigenvectors of matrix A.
    (2) (5%) Please find the general solution of simultaneous homogeneous differential equations X′=AXX' = AX.

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

這一題的完整詳解

核心觀念

對線性齊次常微分方程組

X′=AXX'=AX

若矩陣 AA 有特徵值 λ\lambda 及對應特徵向量 v≠0v\neq 0,則

X(t)=eλtvX(t)=e^{\lambda t}v

是方程組的一組解。若 AA 有足夠多個線性獨立的特徵向量,通解就是各組特徵解的線性組合。

本題重點為:

  1. 求特徵方程 det⁡(A−λI)=0\det(A-\lambda I)=0。
  2. 對每個特徵值求特徵向量。
  3. 利用 eλtve^{\lambda t}v 組合出通解。

解題方法

矩陣為

A=[−1101−1000−2].A= \begin{bmatrix} -1&1&0\\ 1&-1&0\\ 0&0&-2 \end{bmatrix}.

矩陣的前兩列、前兩行形成一個 2×22\times 2 區塊,而第三個變數獨立對應到特徵值 −2-2,因此可直接求特徵值與特徵向量。


(1) 特徵值與特徵向量

求特徵值

計算特徵方程:

det⁡(A−λI)=∣−1−λ101−1−λ000−2−λ∣.\det(A-\lambda I) = \begin{vmatrix} -1-\lambda&1&0\\ 1&-1-\lambda&0\\ 0&0&-2-\lambda \end{vmatrix}.

沿第三欄或第三列展開,得到

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

因此

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

整理為

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

令特徵方程等於零:

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

所以特徵值為

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

其中 λ=−2\lambda=-2 是重根,代數重數為 22。


特徵值 λ=0\lambda=0 的特徵向量

解

(A−0I)v=Av=0.(A-0I)v=Av=0.

令

v=[v1v2v3],v= \begin{bmatrix} v_1\\v_2\\v_3 \end{bmatrix},

則

[−1101−1000−2][v1v2v3]=[000].\begin{bmatrix} -1&1&0\\ 1&-1&0\\ 0&0&-2 \end{bmatrix} \begin{bmatrix} v_1\\v_2\\v_3 \end{bmatrix} = \begin{bmatrix} 0\\0\\0 \end{bmatrix}.

由方程組可得

−v1+v2=0,v1−v2=0,−2v3=0.-v_1+v_2=0,\qquad v_1-v_2=0,\qquad -2v_3=0.

因此

v2=v1,v3=0.v_2=v_1,\qquad v_3=0.

取 v1=1v_1=1,可得一組特徵向量

v1=[110].\boxed{ v_1= \begin{bmatrix} 1\\1\\0 \end{bmatrix} }.

特徵值 λ=−2\lambda=-2 的特徵向量

解

(A+2I)v=0.(A+2I)v=0.

其中

A+2I=[110110000].A+2I= \begin{bmatrix} 1&1&0\\ 1&1&0\\ 0&0&0 \end{bmatrix}.

因此

[110110000][v1v2v3]=[000],\begin{bmatrix} 1&1&0\\ 1&1&0\\ 0&0&0 \end{bmatrix} \begin{bmatrix} v_1\\v_2\\v_3 \end{bmatrix} = \begin{bmatrix} 0\\0\\0 \end{bmatrix},

得到

v1+v2=0.v_1+v_2=0.

所以

v2=−v1,v_2=-v_1,

而 v3v_3 可任意取值。特徵向量的一般形式為

🔒

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

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

免費註冊

第 3 題20 分

  1. (20%) There is a periodic function f(x)=x2f(x) = x^2, (−π<x<π)(-\pi < x < \pi) and f(x)=f(x+2π)f(x) = f(x+2\pi).
    (1) (10%) Please find the Fourier series of f(x)f(x).
    (2) (10%) Please use the Parseval equality to prove the following equality.
    112+122+132+142+⋯=π26\frac{1}{1^2} + \frac{1}{2^2} + \frac{1}{3^2} + \frac{1}{4^2} + \dots = \frac{\pi^2}{6}

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

這一題的完整詳解

(1) 求解 f(x)f(x) 之傅立葉級數

由於 f(x)=x2f(x) = x^2 於 (−π,π)(-\pi, \pi) 為偶函數(Even function),其正弦項係數 bn=0b_n = 0 (n=1,2,3,…n = 1, 2, 3, \dots)。

1. 計算常數項 a0a_0

a0=1π∫−ππx2 dx=2π∫0πx2 dx=2π[x33]0π=2π23a_0 = \frac{1}{\pi} \int_{-\pi}^{\pi} x^2 \, dx = \frac{2}{\pi} \int_{0}^{\pi} x^2 \, dx = \frac{2}{\pi} \left[ \frac{x^3}{3} \right]_0^{\pi} = \frac{2\pi^2}{3}
故直流分量(Constant term)為 a02=π23\frac{a_0}{2} = \frac{\pi^2}{3}。

2. 計算餘弦項係數 ana_n (n≥1n \ge 1)

an=1π∫−ππx2cos⁡(nx) dx=2π∫0πx2cos⁡(nx) dxa_n = \frac{1}{\pi} \int_{-\pi}^{\pi} x^2 \cos(nx) \, dx = \frac{2}{\pi} \int_{0}^{\pi} x^2 \cos(nx) \, dx

利用分部積分法(Integration by parts):
∫0πx2cos⁡(nx) dx=[x2sin⁡(nx)n]0π−∫0π2xsin⁡(nx)n dx\int_{0}^{\pi} x^2 \cos(nx) \, dx = \left[ x^2 \frac{\sin(nx)}{n} \right]_0^{\pi} - \int_{0}^{\pi} 2x \frac{\sin(nx)}{n} \, dx
=0−2n([−xcos⁡(nx)n]0π−∫0π−cos⁡(nx)n dx)= 0 - \frac{2}{n} \left( \left[ -x \frac{\cos(nx)}{n} \right]_0^{\pi} - \int_{0}^{\pi} -\frac{\cos(nx)}{n} \, dx \right)
=−2n(−π(−1)nn+0)=2π(−1)nn2= -\frac{2}{n} \left( -\frac{\pi (-1)^n}{n} + 0 \right) = \frac{2\pi (-1)^n}{n^2}

代回 ana_n 可得:
an=2π⋅2π(−1)nn2=4(−1)nn2a_n = \frac{2}{\pi} \cdot \frac{2\pi (-1)^n}{n^2} = \frac{4(-1)^n}{n^2}

🔒

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

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

免費註冊

第 4 題25 分

  1. (25%) Consider the following differential equation
    dxdt+3x=e−2t,x(0)=1\frac{dx}{dt} + 3x = e^{-2t}, \quad x(0) = 1
    (1) (10%) Please find the homogeneous solution.
    (2) (15%) Please find the exact solution for x(t)x(t) without using Laplace transform.

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

這一題的完整詳解

核心觀念

本題考查一階線性微分方程

dxdt+p(t)x=q(t)\frac{dx}{dt}+p(t)x=q(t)

的解法,以及齊次解與初始條件的使用。

  • 齊次方程:令右側 q(t)=0q(t)=0。
  • 積分因子:
μ(t)=e∫p(t) dt\mu(t)=e^{\int p(t)\,dt}
  • 初始條件 x(0)=1x(0)=1 用來決定積分常數。

解題方法

(1)齊次解

將非齊次項 e−2te^{-2t} 設為 00,得到齊次方程:

dxhdt+3xh=0\frac{dx_h}{dt}+3x_h=0

整理為:

dxhdt=−3xh\frac{dx_h}{dt}=-3x_h

分離變數:

1xh dxh=−3 dt\frac{1}{x_h}\,dx_h=-3\,dt

兩邊積分:

ln⁡∣xh∣=−3t+C\ln|x_h|=-3t+C

因此:

xh(t)=Ce−3tx_h(t)=Ce^{-3t}

其中 CC 為任意常數。


(2)不用 Laplace transform 求 x(t)x(t)

原方程為:

dxdt+3x=e−2t\frac{dx}{dt}+3x=e^{-2t}

這是一階線性微分方程,其中 p(t)=3p(t)=3,因此積分因子為:

μ(t)=e∫3 dt=e3t\mu(t)=e^{\int 3\,dt}=e^{3t}

方程兩側同乘 e3te^{3t}:

e3tdxdt+3e3tx=e3te−2te^{3t}\frac{dx}{dt}+3e^{3t}x=e^{3t}e^{-2t}

左側可寫成乘積微分:

ddt(e3tx)=et\frac{d}{dt}\left(e^{3t}x\right)=e^t

兩邊對 tt 積分:

🔒

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

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

免費註冊

第 5 題25 分

  1. (25%) Consider the following differential equation
    2dx1dt+dx2dt+3x1=1;dx1dt+2dx2dt+4x2=e−t,x1(0)=1;x2(0)=02 \frac{dx_1}{dt} + \frac{dx_2}{dt} + 3x_1 = 1; \quad \frac{dx_1}{dt} + 2 \frac{dx_2}{dt} + 4x_2 = e^{-t}, \quad x_1(0) = 1; x_2(0) = 0
    (1) (5%) Let the Laplace transform of x1(t)x_1(t) be X1(s)X_1(s). Please express the Laplace transform of dx1dt\frac{dx_1}{dt} in terms of X1(s)X_1(s).
    (2) (20%) Please find the Laplace transform of x1(t)x_1(t) and x2(t)x_2(t) directly from the differential equation.

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

這一題的完整詳解

(1) 求 L{dx1dt}\mathcal{L}\left\{\frac{dx_1}{dt}\right\} 之表示式

由拉普拉斯轉換之微分性質:
L{dx1dt}=sX1(s)−x1(0)\mathcal{L}\left\{\frac{dx_1}{dt}\right\} = s X_1(s) - x_1(0)

代入初始條件 x1(0)=1x_1(0) = 1:
L{dx1dt}=sX1(s)−1\mathcal{L}\left\{\frac{dx_1}{dt}\right\} = s X_1(s) - 1


(2) 直接由微分方程求 x1(t)x_1(t) 與 x2(t)x_2(t) 之拉普拉斯轉換

設 L{x1(t)}=X1(s)\mathcal{L}\{x_1(t)\} = X_1(s),L{x2(t)}=X2(s)\mathcal{L}\{x_2(t)\} = X_2(s)。對原微分方程組兩端取拉普拉斯轉換並代入初始條件 x1(0)=1x_1(0) = 1、x2(0)=0x_2(0) = 0:

  1. 第一式:
    2[sX1(s)−1]+[sX2(s)−0]+3X1(s)=1s2[sX_1(s) - 1] + [sX_2(s) - 0] + 3X_1(s) = \frac{1}{s}
    (2s+3)X1(s)+sX2(s)=1s+2=2s+1s(2s + 3)X_1(s) + sX_2(s) = \frac{1}{s} + 2 = \frac{2s + 1}{s}

  2. 第二式:
    [sX1(s)−1]+2[sX2(s)−0]+4X2(s)=1s+1[sX_1(s) - 1] + 2[sX_2(s) - 0] + 4X_2(s) = \frac{1}{s + 1}
    sX1(s)+(2s+4)X2(s)=1s+1+1=s+2s+1sX_1(s) + (2s + 4)X_2(s) = \frac{1}{s + 1} + 1 = \frac{s + 2}{s + 1}

寫成線性方程組之矩陣型態:

🔒

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

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

免費註冊

其他考古題