114 年 國立中央大學統計研究所碩士班在職生《基礎數學》

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

第 1 題16 分

  1. (16%)
    (a) Let g(x)g(x) be a nondecreasing nonnegative function for 0<x<∞0 < x < \infty and 0≤f(x)≤10 \le f(x) \le 1 for 0<x<∞0 < x < \infty. Given g(ϵ)>0g(\epsilon) > 0 with some ϵ>0\epsilon > 0, show that
    ∫g(ϵ)∞f(y)dy=∫ϵ∞g(y)f(y)dy\int_{g(\epsilon)}^\infty f(y) dy = \int_\epsilon^\infty g(y) f(y) dy
    (b) Let α\alpha be a positive integer and λ>0\lambda > 0. Show that
    ∫ϵ∞yα−1e−λydy<(α+1)!λα+2\int_\epsilon^\infty y^{\alpha-1} e^{-\lambda y} dy < \frac{(\alpha+1)!}{\lambda^{\alpha+2}}

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

這一題的完整詳解

核心觀念

本題依照題面所給條件,兩小題的結論皆無法成立。原因如下:

  • (a) 僅知道 gg 非遞減、非負,以及 0≤f≤10\le f\le 1,不足以推出兩個積分相等。
  • (b) 右側的估計未加入 λ\lambda 的限制,因此不可能對所有 λ>0\lambda>0 成立。

(a) 題面結論不成立

取

g(x)=x,f(y)=e−y,ϵ=1.g(x)=x,\qquad f(y)=e^{-y},\qquad \epsilon=1.

則 gg 為非遞減非負函數,且對所有 y>0y>0 有 0<f(y)≤10<f(y)\le 1,同時 g(ϵ)=1>0g(\epsilon)=1>0。

左側為

∫g(ϵ)∞f(y) dy=∫1∞e−y dy=e−1.\int_{g(\epsilon)}^\infty f(y)\,dy =\int_1^\infty e^{-y}\,dy =e^{-1}.

右側為

∫ϵ∞g(y)f(y) dy=∫1∞ye−y dy.\int_\epsilon^\infty g(y)f(y)\,dy =\int_1^\infty ye^{-y}\,dy.

分部積分得

∫1∞ye−y dy=[−ye−y]1∞+∫1∞e−y dy=e−1+e−1=2e.\int_1^\infty ye^{-y}\,dy =\left[-ye^{-y}\right]_1^\infty+\int_1^\infty e^{-y}\,dy =e^{-1}+e^{-1} =\frac{2}{e}.

因此

∫g(ϵ)∞f(y) dy=1e≠2e=∫ϵ∞g(y)f(y) dy.\int_{g(\epsilon)}^\infty f(y)\,dy =\frac1e \ne \frac2e =\int_\epsilon^\infty g(y)f(y)\,dy.

故 (a) 的等式不由題目條件推出。

解題技巧

遇到聲稱「對所有函數皆成立」的積分恆等式,可先取最簡單的函數測試,例如

g(x)=x,f(x)=e−x.g(x)=x,\qquad f(x)=e^{-x}.

若兩側分別出現不同的權重,通常需要額外的變數代換條件或函數關係,不能僅靠單調性推出。


(b) 題面不等式不成立

取

α=1,λ=10,ϵ=0.001.\alpha=1,\qquad \lambda=10,\qquad \epsilon=0.001.

則左側為

∫0.001∞e−10y dy=e−0.0110≈0.099005.\int_{0.001}^{\infty}e^{-10y}\,dy =\frac{e^{-0.01}}{10} \approx 0.099005.

右側為

🔒

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

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

免費註冊

第 2 題24 分

  1. (24%)
    (a) (8%) Show that
    ∫0∞e−x2dx=π2\int_0^\infty e^{-x^2} dx = \frac{\sqrt{\pi}}{2}
    (b) (8%) Show that for all y≥0y \ge 0
    ∫0∞e−y2/xx2dx≤1y\int_0^\infty \frac{e^{-y^2/x}}{x^2} dx \le \frac{1}{y}
    (c) (8%) Show that
    ∫0∞∫0∞e−x2+y2dxdy≤π\int_0^\infty \int_0^\infty e^{-x^2+y^2} dx dy \le \sqrt{\pi}

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

這一題的完整詳解

核心觀念

本題主要考查:

  • 高斯積分 ∫0∞e−x2 dx\displaystyle \int_0^\infty e^{-x^2}\,dx 的極座標計算。
  • 變數代換在廣義積分中的應用。
  • 非負函數二重積分的收斂性判斷,以及 Tonelli 定理。
  • 由積分值進行不等式與發散性的檢查。

需注意:題目第 (b)、(c) 小題依照目前文字敘述皆不成立,無法證明。以下先給出正確計算,並指出矛盾所在。


(a) 證明 ∫0∞e−x2 dx=π2\displaystyle \int_0^\infty e^{-x^2}\,dx=\frac{\sqrt{\pi}}{2}

令

I=∫0∞e−x2 dx.I=\int_0^\infty e^{-x^2}\,dx.

因為 I>0I>0,考慮其平方:

I2=(∫0∞e−x2 dx)(∫0∞e−y2 dy).I^2 = \left(\int_0^\infty e^{-x^2}\,dx\right) \left(\int_0^\infty e^{-y^2}\,dy\right).

因此

I2=∫0∞∫0∞e−(x2+y2) dx dy.I^2 = \int_0^\infty\int_0^\infty e^{-(x^2+y^2)}\,dx\,dy.

在第一象限使用極座標代換

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

其中

0≤r<∞,0≤θ≤π2,0\le r<\infty,\qquad 0\le\theta\le\frac{\pi}{2},

且 Jacobian 為 rr,所以

I2=∫0π/2∫0∞e−r2r dr dθ.I^2 = \int_0^{\pi/2}\int_0^\infty e^{-r^2}r\,dr\,d\theta.

令 u=r2u=r^2,則 du=2r drdu=2r\,dr,故

∫0∞e−r2r dr=12∫0∞e−u du=12.\int_0^\infty e^{-r^2}r\,dr = \frac12\int_0^\infty e^{-u}\,du = \frac12.

因此

I2=∫0π/212 dθ=12⋅π2=π4.I^2 = \int_0^{\pi/2}\frac12\,d\theta = \frac12\cdot\frac{\pi}{2} = \frac{\pi}{4}.

由於 I>0I>0,取正平方根得

∫0∞e−x2 dx=π2.\boxed{ \int_0^\infty e^{-x^2}\,dx=\frac{\sqrt{\pi}}{2} }.

(b) 題目敘述不成立

題目聲稱對所有 y≥0y\ge0,

∫0∞e−y2/xx2 dx≤1y.\int_0^\infty \frac{e^{-y^2/x}}{x^2}\,dx\le\frac1y.

首先,當 y=0y=0 時,右側 1y\frac1y 沒有定義;而左側為

∫0∞1x2 dx,\int_0^\infty \frac{1}{x^2}\,dx,

此積分在 x=0x=0 附近發散。因此 y=0y=0 時不成立。

對 y>0y>0,令

u=y2x,x=y2u,dx=−y2u2 du.u=\frac{y^2}{x}, \qquad x=\frac{y^2}{u}, \qquad dx=-\frac{y^2}{u^2}\,du.

此外,

1x2=u2y4.\frac1{x^2}=\frac{u^2}{y^4}.

當 x→0+x\to0^+ 時,u→∞u\to\infty;當 x→∞x\to\infty 時,u→0+u\to0^+。因此

🔒

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

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

免費註冊

第 3 題16 分

  1. (16%) Give constants 0<α1<10 < \alpha_1 < 1, 0<α2<10 < \alpha_2 < 1, 0<β1<10 < \beta_1 < 1, 0<β2<10 < \beta_2 < 1, and 0<μ<10 < \mu < 1. Minimize β1w1+β2w2\beta_1 w_1 + \beta_2 w_2 with respect to w1w_1 and w2w_2 under two constraints w1+w2=1w_1 + w_2 = 1 and α1w1+α2w2=μ\alpha_1 w_1 + \alpha_2 w_2 = \mu with 0≤w1≤10 \le w_1 \le 1 and 0≤w2≤10 \le w_2 \le 1.
    (a) (8%) Find the minimizer w1w_1 and w2w_2.
    (b) (8%) Find the minimal value of β1w1+β2w2\beta_1 w_1 + \beta_2 w_2.

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

這一題的完整詳解

核心觀念

本題為「受約束的最佳化問題(Constrained Optimization Problem)」,本質上屬於線性聯立方程組的可行性分析與線性規劃(Linear Programming)。

  1. 約束條件的維度分析:
    變數僅有兩個(w1,w2w_1, w_2),但等式約束有兩個: {w1+w2=1α1w1+α2w2=μ\begin{cases} w_1 + w_2 = 1 \\ \alpha_1 w_1 + \alpha_2 w_2 = \mu \end{cases} 當 α1≠α2\alpha_1 \neq \alpha_2 時,這是兩個二元一次線性方程,其解為唯一點(Unique Solution)。
  2. 可行域(Feasible Region)的退化:
    在 0≤w1≤10 \le w_1 \le 1 與 0≤w2≤10 \le w_2 \le 1 的限制下,若求得的唯一解落在 [0,1][0, 1] 區間內,則可行域僅包含「單一點」。此時目標函數沒有多餘的自由度可供極小化,該唯一定點即為唯一的極值點(Minimizer)。
  3. 退化情形(α1=α2\alpha_1 = \alpha_2):
    若 α1=α2=α\alpha_1 = \alpha_2 = \alpha:
    • 當 μ=α\mu = \alpha 時,第二式化為 α(w1+w2)=α\alpha(w_1 + w_2) = \alpha,與第一式重合,可行域為整條線段 w1+w2=1w_1 + w_2 = 1(w1,w2∈[0,1]w_1, w_2 \in [0, 1])。此時極小化問題轉為一般的一維線性規劃。
    • 當 μ≠α\mu \neq \alpha 時,聯立方程組無解,可行域為空集合。

解題方法

步驟一:解約束條件的線性聯立方程組

由題意,變數 w1,w2w_1, w_2 需同時滿足:

w1+w2=1⋯⋯(1)α1w1+α2w2=μ⋯⋯(2)\begin{aligned} w_1 + w_2 &= 1 \quad \cdots\cdots (1) \\ \alpha_1 w_1 + \alpha_2 w_2 &= \mu \quad \cdots\cdots (2) \end{aligned}

將 (1)(1) 式代入 (2)(2) 式,利用 w2=1−w1w_2 = 1 - w_1:

α1w1+α2(1−w1)=μ\alpha_1 w_1 + \alpha_2 (1 - w_1) = \mu

展開並移項整理:

(α1−α2)w1=μ−α2(\alpha_1 - \alpha_2) w_1 = \mu - \alpha_2

步驟二:分情況討論極值點與可行性

情形一:一般情形(α1≠α2\alpha_1 \neq \alpha_2)

由於 α1−α2≠0\alpha_1 - \alpha_2 \neq 0,可直接除以 (α1−α2)(\alpha_1 - \alpha_2) 解出 w1w_1:

w1=μ−α2α1−α2w_1 = \frac{\mu - \alpha_2}{\alpha_1 - \alpha_2}

代回 w2=1−w1w_2 = 1 - w_1:

w2=1−μ−α2α1−α2=α1−μα1−α2w_2 = 1 - \frac{\mu - \alpha_2}{\alpha_1 - \alpha_2} = \frac{\alpha_1 - \mu}{\alpha_1 - \alpha_2}

可行性檢驗(Feasibility Check):
題目要求 0≤w1≤10 \le w_1 \le 1 且 0≤w2≤10 \le w_2 \le 1。

  • 若 α1>α2\alpha_1 > \alpha_2: 0≤μ−α2α1−α2≤1  ⟺  α2≤μ≤α10 \le \frac{\mu - \alpha_2}{\alpha_1 - \alpha_2} \le 1 \iff \alpha_2 \le \mu \le \alpha_1
  • 若 α1<α2\alpha_1 < \alpha_2: 0≤μ−α2α1−α2≤1  ⟺  α1≤μ≤α20 \le \frac{\mu - \alpha_2}{\alpha_1 - \alpha_2} \le 1 \iff \alpha_1 \le \mu \le \alpha_2

因此,只要 μ\mu 介於 α1\alpha_1 與 α2\alpha_2 之間(即 min⁡(α1,α2)≤μ≤max⁡(α1,α2)\min(\alpha_1, \alpha_2) \le \mu \le \max(\alpha_1, \alpha_2)),滿足全部約束條件的點 (w1,w2)(w_1, w_2) 是唯一的。

因為可行解集合為單元素集合,目標函數在該集合上的極小值點即為此唯一解。

將求得的 (w1,w2)(w_1, w_2) 代入目標函數 f(w1,w2)=β1w1+β2w2f(w_1, w_2) = \beta_1 w_1 + \beta_2 w_2:

🔒

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

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

免費註冊

第 4 題16 分

  1. (16%) Given the two by two matrix
    G=(λ00μ)G = \begin{pmatrix} \lambda & 0 \\ 0 & \mu \end{pmatrix}
    with λ>0\lambda > 0 and μ>0\mu > 0.
    (a) (8%) Find the eigendecomposition G=BAB−1G = BAB^{-1}.
    (b) (8%) Evaluate eG=∑n=0∞Gnn!e^G = \sum_{n=0}^\infty \frac{G^n}{n!} using G=BAB−1G = BAB^{-1}. Note that G0=IG^0 = I where II is the identity matrix.

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

這一題的完整詳解

核心觀念

本題考查兩個觀念:

  1. 矩陣的特徵值分解

    若矩陣 GG 可對角化,則可寫成

    G=BAB−1G=BAB^{-1}

    其中 AA 是由特徵值組成的對角矩陣,BB 的欄向量是對應的特徵向量。

  2. 矩陣指數函數

    矩陣指數定義為

    eG=∑n=0∞Gnn!e^G=\sum_{n=0}^{\infty}\frac{G^n}{n!}

    若 G=BAB−1G=BAB^{-1},則

    Gn=BAnB−1G^n=BA^nB^{-1}

    因此

    eG=BeAB−1e^G=Be^AB^{-1}

    這可將矩陣指數問題化為對角矩陣的指數計算。

解題方法

題目給定

G=(λ00μ),λ>0, μ>0G= \begin{pmatrix} \lambda&0\\ 0&\mu \end{pmatrix}, \qquad \lambda>0,\ \mu>0

由於 GG 本身已經是對角矩陣,因此可直接由對角元素讀出特徵值與特徵向量。

(a)求特徵值分解 G=BAB−1G=BAB^{-1}

特徵值由特徵方程決定:

det⁡(G−tI)=det⁡(λ−t00μ−t)\det(G-tI) = \det \begin{pmatrix} \lambda-t&0\\ 0&\mu-t \end{pmatrix}

因此

det⁡(G−tI)=(λ−t)(μ−t)\det(G-tI) = (\lambda-t)(\mu-t)

故特徵值為

t1=λ,t2=μt_1=\lambda,\qquad t_2=\mu

對應的特徵向量分別可取

v1=(10),v2=(01)v_1= \begin{pmatrix} 1\\ 0 \end{pmatrix}, \qquad v_2= \begin{pmatrix} 0\\ 1 \end{pmatrix}

將特徵向量排列成矩陣:

B=(1001)=IB= \begin{pmatrix} 1&0\\ 0&1 \end{pmatrix} =I

特徵值組成的對角矩陣為

A=(λ00μ)A= \begin{pmatrix} \lambda&0\\ 0&\mu \end{pmatrix}

由於 B=IB=I 且 B−1=IB^{-1}=I,所以

BAB−1=IAI=A=(λ00μ)=GBAB^{-1} = IAI = A = \begin{pmatrix} \lambda&0\\ 0&\mu \end{pmatrix} =G

因此

B=(1001),A=(λ00μ)\boxed{ B= \begin{pmatrix} 1&0\\ 0&1 \end{pmatrix}, \qquad A= \begin{pmatrix} \lambda&0\\ 0&\mu \end{pmatrix} }

當 λ=μ\lambda=\mu 時,G=λIG=\lambda I,任意可逆矩陣都可以作為特徵向量矩陣;取 B=IB=I 仍是最簡潔且正確的選擇。

(b)利用 G=BAB−1G=BAB^{-1} 求 eGe^G

由

G=BAB−1G=BAB^{-1}

可得

Gn=(BAB−1)n=BAnB−1G^n=(BAB^{-1})^n=BA^nB^{-1}
🔒

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

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

免費註冊

第 5 題28 分

  1. (28%) Give a full rank matrix
    A=(a11a12a13a14a21a22a23a24a31a32a33a34a41a42a43a44a51a52a53a54)A = \begin{pmatrix} a_{11} & a_{12} & a_{13} & a_{14} \\ a_{21} & a_{22} & a_{23} & a_{24} \\ a_{31} & a_{32} & a_{33} & a_{34} \\ a_{41} & a_{42} & a_{43} & a_{44} \\ a_{51} & a_{52} & a_{53} & a_{54} \end{pmatrix}
    and H=A(ATA)−1ATH = A(A^T A)^{-1} A^T where ATA^T is the transpose of AA and A−1A^{-1} is the inverse of AA.
    (a) (7%) Show that H5H^5 is a positive semi-definite matrix.
    (b) (7%) If rank[H5]=r<5[H^5] = r < 5, show that it has rr eigenvalues equal to unity and 5−r5 - r eigenvalues equal to zero.
    (c) (7%) Show trace[H]=rank[H5][H] = \text{rank}[H^5].
    (d) (7%) Find the eigenvalues for I5−H5I_5 - H^5 where I5I_5 is the 5×55 \times 5 identity matrix.

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

這一題的完整詳解

此題考察矩陣的投影矩陣性質,以及特徵值、秩、跡等概念。矩陣 HH 是一個投影矩陣。

矩陣 AA 是 5×45 \times 4 的矩陣,且是 full rank。這意味著 AA 的秩 rank(A)=4(A)=4。
矩陣 ATAA^T A 是 4×44 \times 4 的矩陣。由於 AA 是 full rank,rank(ATA)=rank(A)=4(A^T A) = \text{rank}(A) = 4。
因此,(ATA)−1(A^T A)^{-1} 存在。

矩陣 H=A(ATA)−1ATH = A(A^T A)^{-1} A^T 是一個 5×55 \times 5 的矩陣。
HH 是 AA 的列空間上的投影矩陣。
性質:

  1. HT=(A(ATA)−1AT)T=(AT)T((ATA)−1)TAT=A((ATA)T)−1AT=A(ATA)−1AT=HH^T = (A(A^T A)^{-1} A^T)^T = (A^T)^T ((A^T A)^{-1})^T A^T = A ((A^T A)^T)^{-1} A^T = A (A^T A)^{-1} A^T = H。
    所以 HH 是對稱矩陣。
  2. H2=HH=(A(ATA)−1AT)(A(ATA)−1AT)=A(ATA)−1(ATA)(ATA)−1ATH^2 = H H = (A(A^T A)^{-1} A^T)(A(A^T A)^{-1} A^T) = A(A^T A)^{-1} (A^T A) (A^T A)^{-1} A^T
    H2=A(ATA)−1I(ATA)−1AT=A(ATA)−1AT=HH^2 = A(A^T A)^{-1} I (A^T A)^{-1} A^T = A(A^T A)^{-1} A^T = H。
    所以 HH 是冪等矩陣(idempotent matrix)。

由 H2=HH^2 = H,我們可以推斷 Hn=HH^n = H 對於所有 n≥1n \ge 1。
因此,H5=HH^5 = H。

(a) 證明 H5H^5 是正半定矩陣(positive semi-definite matrix)。
由於 H5=HH^5 = H,我們需要證明 HH 是正半定矩陣。
對於任意向量 xx,我們需要證明 xTHx≥0x^T H x \ge 0。
xTHx=xT(A(ATA)−1AT)x=(xTA)(ATA)−1(ATx)x^T H x = x^T (A(A^T A)^{-1} A^T) x = (x^T A) (A^T A)^{-1} (A^T x)。
令 y=ATxy = A^T x。則 yy 是一個 4×14 \times 1 的向量。
xTHx=yT(ATA)−1yx^T H x = y^T (A^T A)^{-1} y。
由於 ATAA^T A 是一個對稱矩陣 (因為 (ATA)T=AT(AT)T=ATA(A^T A)^T = A^T (A^T)^T = A^T A),且 rank(ATA)=4(A^T A) = 4,所以 ATAA^T A 是正定的。
因此,(ATA)−1(A^T A)^{-1} 也是正定的。
對於任意向量 y≠0y \ne 0, yT(ATA)−1y>0y^T (A^T A)^{-1} y > 0。
如果 y=0y = 0 (即 ATx=0A^T x = 0), 則 yT(ATA)−1y=0y^T (A^T A)^{-1} y = 0。
所以 xTHx≥0x^T H x \ge 0 對於所有向量 xx。
這證明了 HH 是正半定矩陣。
由於 H5=HH^5 = H, H5H^5 也是正半定矩陣。

【答案】 H5=HH^5 = H。由於 HH 是投影矩陣,它是對稱且冪等的。xTHx=xTA(ATA)−1ATx=(ATx)T(ATA)−1(ATx)x^T H x = x^T A (A^T A)^{-1} A^T x = (A^T x)^T (A^T A)^{-1} (A^T x)。令 y=ATxy = A^T x。則 xTHx=yT(ATA)−1yx^T H x = y^T (A^T A)^{-1} y。矩陣 ATAA^T A 是 4×44 \times 4 的正定對稱矩陣,因此 (ATA)−1(A^T A)^{-1} 也是正定的。所以 yT(ATA)−1y≥0y^T (A^T A)^{-1} y \ge 0。故 HH 是正半定矩陣。

(b) 如果 rank[H5]=r<5[H^5] = r < 5,證明它有 rr 個特徵值等於 1,以及 5−r5-r 個特徵值等於 0。
由於 H5=HH^5 = H,rank[H5]=rank[H][H^5] = \text{rank}[H]。
矩陣 HH 是 AA 的列空間上的投影矩陣。
HH 的像空間(image space)是 AA 的列空間,其維度是 rank(H)(H)。
AA 是 5×45 \times 4 的矩陣,rank(A)=4(A)=4。
HH 的像空間是 AA 的列空間,其維度是 rank(A)=4(A) = 4。
所以,rank[H]=4[H] = 4。
這意味著 r=4r = 4。
題目說 rank[H5]=r<5[H^5] = r < 5,這與 rank[H]=4[H]=4 矛盾,除非 r=4r=4。
如果 r=4r=4,則 5−r=15-r = 1。
這表示 HH 有 4 個特徵值等於 1,1 個特徵值等於 0。

為什麼 rank[H][H] 是 AA 的列空間維度?
Hx=A(ATA)−1ATxH x = A (A^T A)^{-1} A^T x。
令 z=(ATA)−1ATxz = (A^T A)^{-1} A^T x。則 Hx=AzHx = Az。
這表明 HH 的像空間是 AA 的列空間的子空間。
由於 AA 的列空間維度是 4,而 HH 是 5×55 \times 5 矩陣,rank(H)≤5(H) \le 5。
因為 HH 是投影矩陣,其特徵值只能是 0 或 1。
HH 是 5×55 \times 5 矩陣,rank(H)(H) 是非零特徵值的個數。
由於 HH 是 AA 的列空間上的投影,其像空間的維度就是 AA 的列空間維度,即 rank(A)=4(A)=4。
所以 rank[H]=4[H] = 4。
因此,r=4r=4。
HH 的特徵值是 1 (重數為 rank(H)(H)) 和 0 (重數為 5−rank(H)5 - \text{rank}(H))。

🔒

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

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

免費註冊

其他考古題