Speed up the pseudoinverse calculation

Question

Vincent 2015-6-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224404-speed-up-the-pseudoinverse-calculation

评论： Vincent 2015-6-20

Test.zip

Hi,

I am looking for a way to decrease the computational time for some lines called thousands of times. The idea is to improve the computational speed of the calculation of the pseudoinverse (Jplus in the code). Do you think there is a faster way than using the svd decomposition in Matalb ? If you think its fun to try alternatives, you can give a shot (Enclosed is a save of the workspace with two mex functions). Just unzip the folder and run the code below (the mex-files will only work with 64bit Win platform. Sorry for that!! :))

for k=1:1500
            markk=mark(:,k);
            M=MR(1:3,markk,k,s)'; M=M(:);
            markk3D=[markk;markk;markk];
            a=1;er=1;
            while mean(abs(er))>0.005 && a<6
                T=Tags3M(qk,TS);
                T=T(markk3D);
                er=T-M;
                J=TagsJacobian3M(qk,TS);
                J=J(markk3D,:);
                [U,S,V] = svd(J,0);   r = size(S,1);      SMin = S(r,r);
                Jplus=V/S*U';
                qk = qk - Jplus*er; 
                a=a+1;
            end
end