奇异值.m

SVD的应用

m=input('please iuput m number:')
n=input('please iuput n number:')
A=rand(m,n)
[U,S,V]=svd(A)
S1=S(:,1:n-1);
S2=zeros(m,1);
S1=[S1 S2]
A1=U*S1*V'
n1=norm(A-A1);
a=10000;
t=0;
while(a~=0)
B=rand(m,n);
n2=norm(A-B);
if(n1<n2)
    t=t+1;
end
a=a-1;
end
n2=min(n2)
n1
t