% Check the Jacobian of the projection x' = N(PX) px = [vecP;X] y = N(x') Wolfgang Förstner wfoerstn@uni-bonn.de last changes: Susanne Wenzel 04/18 wenzel@igg.uni-bonn.de
0001 %% Check the Jacobian of the projection 0002 % x' = N(PX) 0003 % 0004 % px = [vecP;X] 0005 % y = N(x') 0006 % 0007 % Wolfgang Förstner 0008 % wfoerstn@uni-bonn.de 0009 % 0010 % last changes: Susanne Wenzel 04/18 0011 % wenzel@igg.uni-bonn.de 0012 % 0013 function y = projection_n(px) 0014 0015 y = reshape( px(1:12),3,4 ) * px(13:16); 0016 y = y/norm(y); 0017