博士深造矩阵理论考试题目与答案分享

博士深造矩阵理论考试题目与答案分享考试时间：______分钟总分：______分姓名：______一、填空题1.若矩阵A=[(1,2),(3,4)]，则A的转置矩阵A^T=________。2.矩阵B=[(0,1),(2,0)]的特征值为________，相应的特征向量分别为________和________。3.若方阵A可逆，则其逆矩阵A^(-1)满足关系式AA^(-1)=________。4.实对称矩阵的特征值一定是________，且其一定可以________对角化。5.矩阵C=[(2,0),(0,3)]的特征值为________，它是一个________矩阵（请填写“正定”、“负定”、“半正定”、“半负定”、“不定”）。6.矩阵范数||A||_2表示矩阵A的________范数，它等于矩阵A的奇异值中的________。7.若矩阵X满足AX=λX，其中A是n阶方阵，λ是常数，X是非零列向量，则称X是矩阵A的________，λ是A的________。8.矩阵(A+A^T)(A-A^T)是________矩阵（请填写“对称”、“反对称”、“单位”）。9.若矩阵A的特征值分别为λ1,λ2,...,λn，则矩阵A^T的特征值为________。10.奇异值分解（SVD）将任意m×n矩阵A分解为A=UΣV^T，其中Σ是一个m×n的________矩阵，其对角线元素称为A的________。二、计算题1.计算矩阵A=[(1,2,3),(0,4,5),(0,0,6)]的行列式det(A)。2.求矩阵B=[(2,1),(1,2)]的特征值和特征向量。3.计算矩阵C=[(3,1),(1,3)]的逆矩阵C^(-1)（若存在）。4.对矩阵D=[(1,2),(3,4)]进行QR分解（要求Q为正交矩阵，R为上三角矩阵）。5.给定矩阵E=[(2,1),(1,2)]，求矩阵函数f(E)=e^E（提示：可利用对角化方法）。6.求矩阵F=[(1,0),(0,-1),(2,0)]的奇异值分解（SVD）。三、证明题1.证明：若矩阵A是正定矩阵，则A的所有特征值都是正数。2.证明：任何实对称矩阵都可以正交相似于一个对角矩阵。3.证明：若矩阵A是可逆的，则其奇异值分解中的奇异值均为正数。4.证明：矩阵范数满足三角不等式，即||A+B||≤||A||+||B||。四、综合应用题1.考虑线性方程组Ax=b，其中A=[(1,2),(3,4)],b=[(5),(6)]。(1)判断该线性方程组是否有解。(2)若有解，求其最小二乘解（即求解阿达玛方程AX=B，其中B=[b,0]）。(3)利用奇异值分解（SVD）方法求解该线性方程组的最小二乘解，并分析其解的性质。2.设V是n维欧几里得空间中的一组向量，其格拉姆矩阵（GramMatrix）G=V^TV。证明：G必为半正定矩阵。3.在机器学习中，主成分分析（PCA）涉及协方差矩阵的特征值分解。简述PCA的基本思想，并说明协方差矩阵的特征值和特征向量在PCA中的作用。试卷答案一、填空题1.[(2,3),(0,4)]2.-1,2;[(1,1)],[(1,-1)]3.I(单位矩阵)4.实数，正交相似5.2,3;正定6.2范数（或谱范数），最大奇异值7.特征向量，特征值8.反对称9.λ1,λ2,...,λn10.对角，奇异值二、计算题1.解：det(A)=1*(4*6-5*0)-2*(0*6-5*0)+3*(0*4-0*0)=1*24-2*0+3*0=24。2.解：特征方程|B-λI|=|(2-λ,1),(1,2-λ)|=(2-λ)^2-1=λ^2-4λ+3=(λ-3)(λ-1)=0。特征值为λ1=3,λ2=1。对λ1=3，(B-3I)x=0→[(-1,1),(1,-1)]x=0→x1=x2。取特征向量x1=[(1),(1)]。对λ2=1，(B-1I)x=0→[(1,1),(1,1)]x=0→x1=-x2。取特征向量x2=[(1),(-1)]。特征值为3，特征向量可为[(1),(1)]；特征值为1，特征向量可为[(1),(-1)]。3.解：det(C)=3*3-1*1=9-1=8≠0，故C可逆。C^(-1)=(1/det(C))*adj(C)=(1/8)*[(3,-1),(-1,3)]=[(3/8,-1/8),(-1/8,3/8)]。4.解：计算D的特征值|D-λI|=|(1-λ,2),(3,4-λ)|=(1-λ)(4-λ)-6=λ^2-5λ-2=0。解得λ1≈5.414,λ2≈-0.414。对λ1，(D-λ1I)x=0→[(-4.414,2),(3,4-5.414)]x=0→[(-4.414,2),(3,-1.414)]x=0。解得单位特征向量u1≈[0.447,0.894]。对λ2，(D-λ2I)x=0→[(1.414,2),(3,4.414)]x=0→[(1.414,2),(3,4.414)]x=0。解得单位特征向量u2≈[-0.894,0.447]。令Q=[u1,u2]=[(0.447,-0.894),(0.894,0.447)]，则Q是正交矩阵。R=Q^TD=[(1,3),(-2,-4)]。QR分解为[(1,3),(-2,-4)]=[(0.447,-0.894),(0.894,0.447)][(5.414),(0)]。5.解：E是实对称矩阵，可对角化。求E的特征值|E-λI|=|(2-λ,1),(1,2-λ)|=λ^2-4λ+3=(λ-3)(λ-1)=0。特征值为λ1=3,λ2=1。对λ1=3，(E-3I)x=0→[(-1,1),(1,-1)]x=0→x1=x2。取特征向量v1=[(1),(1)]。对λ2=1，(E-1I)x=0→[(1,1),(1,1)]x=0→x1=-x2。取特征向量v2=[(1),(-1)]。单位化特征向量：u1=v1/||v1||=[(1/sqrt(2)),(1/sqrt(2))]，u2=v2/||v2||=[(1/sqrt(2)),(-1/sqrt(2))].令P=[u1,u2]=[(1/sqrt(2),1/sqrt(2)),(1/sqrt(2),-1/sqrt(2))]，则P是正交矩阵，且P^TEP=Λ=[(3),(1)]。f(E)=e^E=Pe^(Λ)P^T=P[(e^3),(e^1)]P^T=[(1/sqrt(2),1/sqrt(2)),(1/sqrt(2),-1/sqrt(2))][(e^3),(e)][(1/sqrt(2)),(1/sqrt(2))][(1/sqrt(2)),(-1/sqrt(2))]。计算得f(E)=[(e^3+e)/2,(e^3-e)/2],[(e^3-e)/2,(e^3+e)/2]。6.解：计算F^TF=[(1,0,2),(0,1,0)][(1,0),(0,-1),(2,0)]=[(1,0,2),(0,-1,0),(2,0,4)]。对矩阵F^TF进行特征值分解。计算特征值|(F^TF-λI)|=|(1-λ,0,2),(0,-1-λ,0),(2,0,4-λ)|。按第二行展开=-(-1-λ)|(1-λ,2),(2,4-λ)|=-(λ+1)[(1-λ)(4-λ)-4]=-(λ+1)(λ^2-5λ)=-(λ+1)λ(λ-5)。特征值为λ1=0,λ2=5,λ3=0。对λ1=0，(F^TF-0I)x=0→[(1,0,2),(0,-1,0),(2,0,4)]x=0。解得x=[0,1,0]^T。对λ2=5，(F^TF-5I)x=0→[(-4,0,2),(0,-6,0),(2,0,-1)]x=0。解得x=[(1/2),(0),(-1)]^T。单位化得v2=[(1/sqrt(2)),(0),(-1/sqrt(2))]^T。对λ3=0，同λ1，解得x=[0,1,0]^T。需与v2正交，取x=[(1/2),(0),(1/2)]^T。单位化得v3=[(1/sqrt(2)),(0),(1/sqrt(2))]^T。令V=[v1,v2,v3]=[(0,1/sqrt(2),1/sqrt(2)),(1,0,0),(0,-1/sqrt(2),1/sqrt(2))]。Σ=diag(σ1,σ2,σ3)=diag(sqrt(5),0,0)=[(sqrt(5),0,0),(0,0,0),(0,0,0)]。U是F的左奇异向量矩阵，由F的奇异向量构成。F^TFv2=σ2^2v2=5v2，故σ2=5，对应的右奇异向量为v2。F^TFv3=σ3^2v3=0，故σ3=0，对应的右奇异向量为v3。F的奇异向量为V的列向量，奇异值为sqrt(5),0,0。SVD分解为F=UΣV^T=[(0,1/sqrt(2),1/sqrt(2)),(1,0,0),(0,-1/sqrt(2),1/sqrt(2))][(sqrt(5),0,0),(0,0,0),(0,0,0)][(0,1/sqrt(2),1/sqrt(2)),(1,0,0),(0,-1/sqrt(2),1/sqrt(2))]。三、证明题1.证明：设λ是正定矩阵A的任意一个特征值，对应的特征向量为x(x≠0)。由A的正定性，对任意非零向量x，有x^TAx>0。令x=v，则v^TAv>0。取x=x/||x||，则(x/||x||)^TA(x/||x||)=(1/||x||^2)x^TAx>0。即x^TAx=||x||^2*(x/||x||)^TA(x/||x||)=||x||^2*(x^TAx)/||x||^2=x^TAx>0。由特征值定义，x^TAx=λ||x||^2。因此，λ||x||^2>0。由于||x||^2>0，故λ>0。由于λ的任意性，正定矩阵A的所有特征值均为正数。2.证明：设A是n阶实对称矩阵。对任意的实向量x，有x^TAx是实数。令x=Au+v，其中u,v是实向量，Au是A的实向量像，v是任意实向量。x^TAx=(Au+v)^TA(Au+v)=u^TA^TAu+2u^TA^Tv+v^TAv=u^TA^2u+2u^TAv+v^TAv。由于A是对称的，A^T=A，A^2也是对称的，Av是实向量。所以u^TA^2u,v^TAv是实数，2u^TAv也是实数。因此，x^TAx是实数。这表明A是厄米矩阵。根据谱定理，厄米矩阵（即实对称矩阵）可正交相似于一个对角矩阵。即存在正交矩阵Q(Q^TQ=I)，使得Q^TAQ=Λ，其中Λ是对角矩阵，其对角线元素是A的特征值。3.证明：设A是m×n矩阵，其奇异值分解为A=UΣV^T，其中Σ=diag(σ1,...,σn)且σ1≥σ2≥...≥σn≥0。A的奇异值σi是A^TA的特征值的平方根。因为A^TA=VΣU^T(UΣV^T)^T=VΣ^2V^T。Σ^2是n阶对角矩阵，其对角线元素为σi^2。因此，A^TA的特征值为σ1^2,...,σn^2。它们都是非负实数。由于A可逆当且仅当其秩为n且σ1,...,σn均不为0。当A可逆时，其奇异值σ1,...,σn均为正数。它们的平方σ1^2,...,σn^2也是正数。故A的奇异值分解中的奇异值均为非负数。若A可逆，则均为正数。4.证明：设A,B是m×n矩阵。根据矩阵范数的定义，对任意向量x，有||Ax||≤||A||||x||和||Bx||≤||B||||x||。考虑向量x=A^Ty+B^Tz，其中y,z是任意实数。||A+B||=||(A+B)(x)||=||A(A^Ty+B^Tz)+B(A^Ty+B^Tz)||。=||AA^Ty+AB^Tz+BA^Ty+BB^Tz||。≤||AA^Ty||+||AB^Tz||+||BA^Ty||+||BB^Tz||（根据三角不等式）≤||A||||A^Ty||+||A||||B^Tz||+||B||||A^Ty||+||B||||B^Tz||（根据向量范数定义）≤||A||||y||+||A||||z||+||B||||y||+||B||||z||（根据向量范数定义，||A^Ty||≤||y||,||B^Tz||≤||z||）=(||A||+||B||)(||y||+||z||)。由于||y||+||z||是任意实数的绝对值之和，可以取y=1,z=1，此时上式达到最大值||A||+||B||。因此，||A+B||≤||A||+||B||。四、综合应用题1.解：(1)系数矩阵A=[(1,2),(3,4)],增广矩阵(A|b)=[(1,2,5),(3,4,6)]。高斯消元：(A|b)→[(1,2,5),(0,-2,-9)]。第二行对应-2x2=-9，即x2=9/2。代入第一行：x1+2(9/2)=5，即x1+9=5，得x1=-4。方程组有唯一解(x1,x2)=(-4,9/2)。(2)最小二乘问题为求解AX=B，其中A=[(1,2),(3,4)],X=[(x1),(x2)],B=[(5),(6)]。阿达玛方程为(A^TA)X=A^TB。A^TA=[(1,3),(2,4)][(1,2),(3,4)]=[(10,14),(14,20)]。A^TB=[(1,3),(2,4)][(5),(6)]=[(23),(34)]。解方程组：(A^TA)X=A^TB→[(10,14),(14,20)][(x1),(x2)]=[(23),(34)]。对增广矩阵进行行变换：[(10,14,23),(14,20,34)]→[(10,14,23),(0,-6/5,-4)]。→[(10,14,23),(0,1,20/3)]。→[(10,0,-47/3),(0,1,20/3)]。→[(1,0,-47/30),(0,1,20/3)]。得解X=[(-47/30),(20/3)]。即最小二乘解为(x1,x2)=(-47/30,20/3)。(3)SVD方法：A的SVD分解为A=UΣV^T。已计算F的SVD为F=UΣV^T，其中U=[(0,1/sqrt(2),1/sqrt(2)),(1,0,0),(0,-1/sqrt(2),1/sqrt(2))]，Σ=[(sqrt(5),0,0),(0,0,0),(0,0,0)],V^T=[(0,1/sqrt(2),1/sqrt(2)),(1,0,0),(0,-1/sqrt(2),1/sqrt(2))]。对于Ax=b，最小二乘解x是方程X=Σ^(-1)U^Tb的解（当Σ的最大奇异值非零时）。Σ^(-1)=diag(1/sqrt(5),0,0)。U^Tb=[(0,1/sqrt(2),1/sqrt(2))][(5),(6),(0)]=[(6/sqrt(2)),(5/sqrt(2)),(5/sqrt(2))]。Σ^(-1)U^Tb=[(6/sqrt(2))/sqrt(5),0,0]=[(6/sqrt(10)),0,0]。X=Σ^(-1)U^Tb=[(6/sqrt(10)),0,0]^T。V^TX=[(0,1/sqrt(2),1/sqrt(2))][(6/sqrt(10)),0,0]^T=[(0),(6/sqrt(20)),(6/sqrt(20))]=[(0),(3/sqrt(5)),(3/sqrt(5))]。A最小二乘解x=VX=[(0),(3/sqrt(5)),(3/sqrt(5))][(0,1/sqrt(2),1/sqrt(2)),(1,0,0),(0,-1/sqrt(2),1/sqrt(2))]。x=[(3/sqrt(5)*(-1/sqrt(2))),(3/sqrt(5)*(1/sqrt(2))),(3/sqrt(5)*(1/sqrt(2)))]=[(-3/sqrt(10)),(3/sqrt(10)),(3/sqrt(10))]。由于UΣV^T中的Σ包含零奇异值，最小二乘解x=VΣ^(-1)U^Tb=[(-47/30),(20/3)]与阿达玛方程解相同。解的性质：此解是最小范数解，即满足所有约束条件且x的欧氏范数最小的解。2.证明：设V是n维欧几里得空间中的一组向量，其格拉姆矩阵G=V^TV。证明G是半正定矩阵。证明方法一：对任意向量x∈R^n，有x^TGx=x^T(V^TV)x=(x^TV)(V^Tx)=||V^Tx||^2。向量范数的平方总是非负的，即||V^Tx||^2≥0。因此，对任意x∈R^n，都