对一个m*n的矩阵
1. 当 m > n 时
A = rand(m,n)
[U,D,V] = svd(A)
D = D(1:n,1:n)
D = D\eye(n)
U = U(:,1:n)
pinv1 = V*D*U'
2. 当 m < n 时
A = rand(m,n)
[U,D,V] = svd(A)
D = D(1:m,1:m)
D = D\eye(m)
V = V(:,1:m)
pinv1 = V*D*U'
3. 当 m = n 时
A = rand(m)
[U,D,V] = svd(A)
D = D\eye(m)
pinv1 = V*D*U'