Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- SavedA=[170,70;156,56;180,80;165,65;172,72;171,65;190,80;192,110;160,96];
- A=[170,70;156,56;180,80;165,65;172,72;171,65;190,80;192,110;160,96];
- m1=mean(A(:,1));
- m2=mean(A(:,2));
- A(:,1)=A(:,1)-m1;
- A(:,2)=A(:,2)-m2;
- C=A'*A;
- v=rand(2,1)
- for k = 1:20
- w=C*v;
- lambda=v'*w;
- v=w/norm(w);
- end
- v
- v=v/norm(v);
- u=A*v;
- u=u/norm(u);
- lambda = v'*(C*v)/(v'*v);
- s=sqrt(lambda);
- P=s*u*v';
- P(:,1)=P(:,1)+m1;
- P(:,2)=P(:,2)+m2;
- plot(P(1:9),P(10:18))
- scatter(SavedA(1:9),SavedA(10:18))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement