跳到主要内容

東北大学 工学研究科 電気・情報系 2017年3月実施 基礎科目 問題6 数学基礎

Author​

祭音Myyura (co-authored with GPT 5.6 SOL)

Description​

日本語版​

(1) n×nn\times n エルミート行列 AA、nn 次元複素列ベクトル x\boldsymbol x、および二次形式 f(x)=x∗Axf(\boldsymbol x)=\boldsymbol x^*A\boldsymbol x を考える。ここで、x∗\boldsymbol x^* は x\boldsymbol x の複素共役転置(随伴行列)を表す。零ベクトルでない全ての x\boldsymbol x に対して f(x)≥0f(\boldsymbol x)\ge0 を満たすとき、AA は半正定値行列であると言う。次の問に答えよ。

(a) ある正方行列 UU が存在して A=U∗UA=U^*U と表せるとき、AA が半正定値行列であることを示せ。

(b) AA の固有値がすべて実数であることを示せ。

(c) AA の異なる固有値に対応する固有ベクトルは直交することを示せ。

(d) AA の固有値がすべて非負であるとき、AA は半正定値行列であることを示せ。必要に応じて、エルミート行列はユニタリー行列によって対角化可能であるという事実を用いてよい。

(e) 次の関係式を満たす 3×33\times3 実対称行列 BB を求めよ。

x12+x22+x32−x1x2−x2x3=(x1x2x3)B(x1x2x3)x_1^2+x_2^2+x_3^2-x_1x_2-x_2x_3=(x_1\quad x_2\quad x_3)B\begin{pmatrix}x_1\\x_2\\x_3\end{pmatrix}

また、BB の固有値をすべて求め、BB が半正定値行列かどうか判定せよ。

(2) 複素変数 zz の関数 g(z)=z(z2−6z+1)2g(z)=\dfrac{z}{(z^2-6z+1)^2} を考える。

CC は z=e2iθz=e^{2i\theta}(0≤θ≤π0\le\theta\le\pi)により表される円周上を正の向きに回る積分路である。ii は虚数単位である。次の問に答えよ。

(a) 関数 g(z)g(z) のすべての孤立特異点とその留数を求めよ。

(b) 複素積分 ∫Cg(z) dz\displaystyle\int_Cg(z)\,dz を求めよ。

(c) 実定積分 ∫0πdθ(1+sin⁡2θ)2\displaystyle\int_0^\pi\frac{d\theta}{(1+\sin^2\theta)^2} を求めよ。

题目描述​

  1. 设 AA 为 n×nn\times n Hermitian 矩阵,f(x)=x∗Axf(\boldsymbol x)=\boldsymbol x^*A\boldsymbol x。若对任意非零复向量 x\boldsymbol x 都有 f(x)≥0f(\boldsymbol x)\ge0,称 AA 半正定。
    • (a) 若存在方阵 UU 使 A=U∗UA=U^*U,证明 AA 半正定。
    • (b) 证明 AA 的全部特征值为实数。
    • (c) 证明属于不同特征值的特征向量正交。
    • (d) 若全部特征值非负,证明 AA 半正定;可用 Hermitian 矩阵可酉对角化。
    • (e) 求实对称矩阵 BB,使 x12+x22+x32−x1x2−x2x3=xTBxx_1^2+x_2^2+x_3^2-x_1x_2-x_2x_3=\boldsymbol x^TB\boldsymbol x;求其全部特征值并判断半正定性。
  2. 设 g(z)=z/(z2−6z+1)2g(z)=z/(z^2-6z+1)^2,C:z=e2iθC:z=e^{2i\theta},0≤θ≤π0\le\theta\le\pi,取正向。
    • (a) 求全部孤立奇点及留数。
    • (b) 求 ∮Cg(z) dz\oint_Cg(z)\,dz。
    • (c) 求 ∫0π(1+sin⁡2θ)−2 dθ\displaystyle\int_0^\pi(1+\sin^2\theta)^{-2}\,d\theta。

Kai​

(1)​

(a) x∗Ax=(Ux)∗(Ux)=∥Ux∥2≥0\boldsymbol x^*A\boldsymbol x=(U\boldsymbol x)^*(U\boldsymbol x)=\|U\boldsymbol x\|^2\ge0。

(b) 若 Av=λvAv=\lambda v,则 λv∗v=v∗Av=v∗Av‾=λˉv∗v\lambda v^*v=v^*Av=\overline{v^*Av}=\bar\lambda v^*v,故 λ∈R\lambda\in\mathbb R。

(c) 若 Au=λu,Av=μvAu=\lambda u,Av=\mu v,则 λu∗v=(Au)∗v=u∗Av=μu∗v\lambda u^*v=(Au)^*v=u^*Av=\mu u^*v。当 λ≠μ\lambda\ne\mu 时,u∗v=0u^*v=0。

(d) 写成 A=Qdiag⁡(λ1,…,λn)Q∗A=Q\operatorname{diag}(\lambda_1,\ldots,\lambda_n)Q^*,令 y=Q∗xy=Q^*x,则 x∗Ax=∑jλj∣yj∣2≥0x^*Ax=\sum_j\lambda_j|y_j|^2\ge0。

(e)

B=(1−1/20−1/21−1/20−1/21),det⁡(λI−B)=(λ−1)((λ−1)2−12).B=\begin{pmatrix}1&-1/2&0\\-1/2&1&-1/2\\0&-1/2&1\end{pmatrix},\qquad \det(\lambda I-B)=(\lambda-1)\left((\lambda-1)^2-\frac12\right).

特征值为 1,1±1/2\boxed{1,1\pm1/\sqrt2},均正,所以 BB 正定,因而半正定。

(2)​

令 a=3−22,b=3+22a=3-2\sqrt2,b=3+2\sqrt2,两者均为二阶极点。

Res⁡(g,a)=[ddzz(z−b)2]z=a=−a+b(a−b)3=32128,\operatorname{Res}(g,a)=\left[\frac{d}{dz}\frac{z}{(z-b)^2}\right]_{z=a} =-\frac{a+b}{(a-b)^3}=\frac{3\sqrt2}{128},
Res⁡(g,b)=−32128.\operatorname{Res}(g,b)=-\frac{3\sqrt2}{128}.

单位圆仅包含 aa,所以

∮Cg(z) dz=32πi64.\boxed{\oint_Cg(z)\,dz=\frac{3\sqrt2\pi i}{64}}.

由 z=e2iθz=e^{2i\theta} 得 dz=2iz dθdz=2iz\,d\theta,且

z2−6z+1=−4z(1+sin⁡2θ).z^2-6z+1=-4z(1+\sin^2\theta).

因而

∮Cg(z) dz=i8∫0πdθ(1+sin⁡2θ)2,\oint_Cg(z)\,dz=\frac{i}{8}\int_0^\pi\frac{d\theta}{(1+\sin^2\theta)^2},

故 ∫0πdθ(1+sin⁡2θ)2=3π42\boxed{\displaystyle\int_0^\pi\frac{d\theta}{(1+\sin^2\theta)^2}=\frac{3\pi}{4\sqrt2}}。