跳到主要内容

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

Author

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

Description

日本語版

(1) n×nn\times n エルミート行列 AAnn 次元複素列ベクトル x\boldsymbol x、および二次形式 f(x)=xAxf(\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=UUA=U^*U と表せるとき、AA が半正定値行列であることを示せ。

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

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

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

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

x12+x22+x32x1x2x2x3=(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(z26z+1)2g(z)=\dfrac{z}{(z^2-6z+1)^2} を考える。

CCz=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+sin2θ)2\displaystyle\int_0^\pi\frac{d\theta}{(1+\sin^2\theta)^2} を求めよ。

题目描述

  1. AAn×nn\times n Hermitian 矩阵,f(x)=xAxf(\boldsymbol x)=\boldsymbol x^*A\boldsymbol x。若对任意非零复向量 x\boldsymbol x 都有 f(x)0f(\boldsymbol x)\ge0,称 AA 半正定。
    • (a) 若存在方阵 UU 使 A=UUA=U^*U,证明 AA 半正定。
    • (b) 证明 AA 的全部特征值为实数。
    • (c) 证明属于不同特征值的特征向量正交。
    • (d) 若全部特征值非负,证明 AA 半正定;可用 Hermitian 矩阵可酉对角化。
    • (e) 求实对称矩阵 BB,使 x12+x22+x32x1x2x2x3=xTBxx_1^2+x_2^2+x_3^2-x_1x_2-x_2x_3=\boldsymbol x^TB\boldsymbol x;求其全部特征值并判断半正定性。
  2. g(z)=z/(z26z+1)2g(z)=z/(z^2-6z+1)^2C: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+sin2θ)2dθ\displaystyle\int_0^\pi(1+\sin^2\theta)^{-2}\,d\theta

Kai

(1)

(a) xAx=(Ux)(Ux)=Ux20\boldsymbol x^*A\boldsymbol x=(U\boldsymbol x)^*(U\boldsymbol x)=\|U\boldsymbol x\|^2\ge0

(b)Av=λvAv=\lambda v,则 λvv=vAv=vAv=λˉvv\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,则 λuv=(Au)v=uAv=μuv\lambda u^*v=(Au)^*v=u^*Av=\mu u^*v。当 λμ\lambda\ne\mu 时,uv=0u^*v=0

(d) 写成 A=Qdiag(λ1,,λn)QA=Q\operatorname{diag}(\lambda_1,\ldots,\lambda_n)Q^*,令 y=Qxy=Q^*x,则 xAx=jλjyj20x^*Ax=\sum_j\lambda_j|y_j|^2\ge0

(e)

B=(11/201/211/201/21),det(λIB)=(λ1)((λ1)212).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=322,b=3+22a=3-2\sqrt2,b=3+2\sqrt2,两者均为二阶极点。

Res(g,a)=[ddzz(zb)2]z=a=a+b(ab)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=2izdθdz=2iz\,d\theta,且

z26z+1=4z(1+sin2θ).z^2-6z+1=-4z(1+\sin^2\theta).

因而

Cg(z)dz=i80πdθ(1+sin2θ)2,\oint_Cg(z)\,dz=\frac{i}{8}\int_0^\pi\frac{d\theta}{(1+\sin^2\theta)^2},

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