Matlab-Moore-Penrose广义逆矩阵

>> clear
A=[16 2 3 13;
5 11 10 8;
9 7 6 12;
4 14 15 1]
A =
    16     2     3    13
     5    11    10     8
     9     7     6    12
     4    14    15     1
>> [U,S,V]=svd(A)
U =
   -0.5000    0.6708    0.5000   -0.2236
   -0.5000   -0.2236   -0.5000   -0.6708
   -0.5000    0.2236   -0.5000    0.6708
   -0.5000   -0.6708    0.5000    0.2236
S =
   34.0000         0         0         0
         0   17.8885         0         0
         0         0    4.4721         0
         0         0         0    0.0000
V =
   -0.5000    0.5000    0.6708   -0.2236
   -0.5000   -0.5000   -0.2236   -0.6708
   -0.5000   -0.5000    0.2236    0.6708
   -0.5000    0.5000   -0.6708    0.2236
>> for i=1:3
S(i,i)=1/S(i,i);
end
>> S
S =
    0.0294         0         0         0
         0    0.0559         0         0
         0         0    0.2236         0
         0         0         0    0.0000
>> pinv(A)
a