研究生矩阵理课后答案4,5章习题.pptVIP

习题4-1(1) 4-1:求 A= 的满秩分解. 解1: A ?? ?? ?? ?? = C ∴ A=BC, B=(A5,A3,A1)= 习题4-1(1) 4-1:求 A= 的满秩分解. 解2: A ?? ?? ?? ?? = C ∴ A=BC, B=(A1,A2,A3)= 习题4-1(2) 4-1(2):求 A= 的满秩分解. 解: A ?? ?? = C ∴ A=BC, B=(A1,A3)= 习题4-2 求 A= 的奇异值分解. 解: A的奇异值是:?2,1; ?=diag(?2,1) AA*的对应于特征值2,1的单位特征向量是 (1/?2,1/?2,0)T, (1,0,0)T 习题4*1A与B酉等价?A与B奇异值相同 ?必要性: A=UBV ? AA*=UBVV*B*U*=UBB*U* ? BB* ∴ AA*与BB*有相同的特征值集,得证A与B有相同的奇异值集. ?充分性:作A,B的奇异值分解 A=UDV*,B=U1DV1*,D=diag(?,0), 其中,?是由它们的全部正奇异值组成的正对角矩阵.于是 U*AV=D=U1*BV1 ? A=(UU1*)B(V1V*) 因酉矩阵的乘积 UU1*,V1V* 仍为酉矩阵,故上式表明A酉等价于B. 习题4*2 4*2: 设A?Crm?n,U?Um?m,V?Un?n使 B=U*AV=diag(?,0),?=diag(b1,…,br), (*) 则|b1|,…,|br|为A的全部正奇异值. 证: U*AA*U=BB*=diag(??*,0) ? 写成?2不对！ =diag(|b1|2,…,|br|2,0,…,0) ∽ AA* ∴ |b1|,…,|br|为A的全部正奇异值. 奇异值分解定理另一(更强)表述 定理: 令?1,…,?r为A?Crm?n的全部正奇异值; ?=diag(?1,…,?r),则有U?Um?m,V?Un?n使 U*AV= =D?Crm?n (*) 反之,若有U?Um?m,V?Un?n使(*)成立,其中?=diag(d1,…,dr),?i,di0,则d1,…,dr为A的全部正奇异值.(奇异值分解的某种唯一性) 证: AA*=U V*V U*=U U* ∽ diag(d12,…,dr2,0,…,0) ∴ d1,…,dr为A的全部正奇异值. 注:后半部等价于补充题4*2. 4*3已知A奇异值求AT,A*,A-1的奇异值 补充题4*3: 令?1,…,?r为A?Crm?n的全部正奇异值; ?=diag(?1,…,?r),则有U?Um?m,V?Un?n使 A=U V*=Udiag(?,0)V* (*) 易见 A*=Vdiag(?,0)U* AT=(Udiag(?,0)V*)T=(V*)Tdiag(?,0)UT ∴ ?1,…,?r为A*,AT, 的全部正奇异值(利用奇异值分解定理的更强表述). A-1=(U?V*)-1=V?-1U*=Vdiag(?1-1,…,?n-1)U* ∴ ?1-1,…,?n-1为A-1的全部正奇异值. 习题#5-1(2) 试证:?x,y?V,‖x?y‖?|‖x‖-‖y‖|. 证：首先‖x‖=‖(x-y)+y‖?‖x-y‖+‖y‖ ? ‖x-y‖?‖x‖-‖y‖. 其次‖x-y‖=‖-(y-x)‖=‖y-x‖ ?‖y‖-‖x‖= -(‖x‖-‖y‖) ∴ ‖x-y‖?|‖x‖-‖y‖|. 此外 ‖x+y‖=‖x-(-y)‖ ?|‖x‖-‖-y‖|=|‖x‖-‖y‖| ∴ ‖x?y‖?|‖x‖-‖y‖|. 习题#5-2试证‖A‖= n maxi,j|aij|是矩阵范数 ?A=(aij)?Cn?n 证: ?非负性,齐次性显然 ?三角不等式: ‖A+B‖= n maxi,j|aij+bij| ? n maxi,j|aij|+n maxi,j|bij|=‖A‖+‖B‖ ?相容性: ‖AB‖= n maxi,j|ai1b1j+…+ainbnj| ? n2 maxi,t|ait| maxtj|btj| =n maxi,j|aij|(n maxi,j|bij|)=‖A‖‖