南京邮电大学matlab软件设计(超详细).doc(4)

njupt

A=[1,2,3;4,5,6;7,8,9] %一个3*3矩阵 disp('设 [U,S,V]=svd(A)，得')

[U,S,V]=svd(A) %奇异值分解(SVD)USV:U m*m酉矩阵 V n*n酉矩阵 S 对角阵

disp('设 tol=S(1,1)*3*eps，得') tol=S(1,1)*3*eps

disp('设 a=sum(diag(S)>tol) - rank(A)，求 a') a=sum(diag(S)>tol) - rank(A)

disp('设 b=S(1,1) - norm(A)，求 b') b=S(1,1) - norm(A)

disp('设 c=sqrt(sum(diag(S*S))) - norm(A, fro )，求 c') c=sqrt(sum(diag(S*S))) - norm(A,'fro') disp('设 d=S(1,1)/S(3,3) - cond(A)，求 d') d=S(1,1)/S(3,3) - cond(A)

disp('设 e=S(1,1)*S(2,2)*S(3,3) - det(A) < tol ，求 e') e=S(1,1)*S(2,2)*S(3,3) - det(A) < tol disp('设 f=V(:,1) *null(A)，求 f') f=V(:,1)'*null(A)

disp('设 g=abs(A*null(A)) < tol，求 g') g=abs(A*null(A)) < tol

disp('设 h=U(:,1:2) == orth(A)，求 h') h=U(:,1:2) == orth(A) 结果：

设 A=[1,2,3;4,5,6;7,8,9]，得 A =

1 2 3 4 5 6 7 8 9

设 [U,S,V]=svd(A)，得 U =

-0.2148 0.8872 0.4082 -0.5206 0.2496 -0.8165 -0.8263 -0.3879 0.4082 S =

16.8481 0 0 0 1.0684 0

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库南京邮电大学matlab软件设计(超详细).doc(4)在线全文阅读。