Ngày tải lên :
10/04/2013, 00:31
... centroid =mean(m,2); T1(1:d-1,d) =-centroid(1:d-1); else T1(1:d-1,d) =-orig(1:d-1); end m = T 1*m; % Scaling scale =rms / sqrt(mean(sum(m(1:d-1,:).J\2,1»); T2 = [scale*eye(d-1)zeros(d-1,1);zeros(1,d-1) ... []; A(1:n,:) = [zeros(n,4) A(n+l:2*n,:) = [(mi(3,:)'*ones(1,4».*Mi' zeros(n,4) -(mi(1,:)'*ones(l,4».*Mi']; A(2*n+ 1:3*n,:) =[-(mi(2,:)'*ones(1,4» *Mi' (mi(1,:)'*ones(l,4» *Mi' zeros(n,4 )]; PI -(mi(3,:)'*ones(1,4» ... tr(in d6ng anh bie'n diim anh VaG{e1,e2,e3,e4} % A: ma tr(in 3x3 hI = cross(m(:,2), m(:,3)); h2 = cross(m(:,3), m(:,I)); h3 = cross(m(:,I), m(:,2)); % L(ip h~ thudn nha't Aa A =0 =[ dot(h2,m(:,4))...