This preview shows pages 1–2. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.View Full Document
Unformatted text preview: val = ( x(:,1,i) - xm ) .^2 + ( x (:,2,i) - ym ) .^2 ; val = sum (val); srval = sqrt(val); G = cat (1, G , srval); xj = ( x(:,1,i) - xm ) ./ srval; yj = ( x(:,2,i) - ym ) ./ srval; x(:,1,i) = xj; x(:,2,i) = yj; end for k =1:40 z = corr( x(:,:,1), x(:,:,k)); [U S V] = svd (z); R = V*U'; if (k == 2) R end M = sum( diag(S) ); ScF = M * G(1) / G(k); if (k ==2) ScF end T = mean(x(:,:,1)) - ScF * mean(x(:,:,k)) * R; if (k == 2) T end PM = ScF * x(:,:,k) * R ; for i = 1:64 PM(i,:) = PM(i,:) +T; end PMS = cat(3, PMS , PM); end z =PMS; end...
View Full Document
This note was uploaded on 11/06/2010 for the course CSE 527 taught by Professor Ab during the Fall '09 term at Cornell University (Engineering School).
- Fall '09