%求解變換矩陣T和T'
function [T,Q]=transform(P,n)
%P為n×2的矩陣 每列表示一點坐標值

%求點集的形心 并將形心移到坐標原點
ave=mean(P');
P(1,:)=P(1,:)-ave(1);
P(2,:)=P(2,:)-ave(2);

%對點進行縮放使它們到原點的平均距離等于2^(1/2)
d_sum=0;

for i=1:n
    d_sum=d_sum+(P(1,i).^2+P(2,i).^2);
end
    d_ave=d_sum/n;

%times為縮放的倍數
s=d_ave/2;
s=sqrt(s);
%T為變換矩陣
T=[1/s 0 -ave(1)/s;0 1/s -ave(2)/s;0 0 1];
Q=T*[P;ones(1,n)];