（5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？
（6） V(:,1)'*null(A) 得到什么结果 ?
（7） abs(A*null(A)) < tol 得到什么结果 ?
（8） U(:,1:2) = = orth(A) 的运行结果是什么 ?
1) M 文件代码[No3code.m]
clear;
clc;
disp('设 A=[1,2,3;4,5,6;7,8,9]，得')
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)
2) 实验结果[No3result.txt]
设 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
