is there any correction..?

Question

Lalit Patil 2013-2-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/63905-is-there-any-correction

I am doing Singular value decomposition using this.

 Xi = [1 2;3 4]; Yi = [1 2;3 4]; Zi = [1 2;3 4];
 A=[Xi,Yi,Zi,ones(length(Xi),1)];
 [U S V] = svd(A)

Now, i want to do the same without function. so, i written little as below.

 Xi = [1 2;3 4]; Yi = [1 2;3 4]; Zi = [1 2;3 4];
 A=[Xi,Yi,Zi,ones(length(Xi),1)];
 f = eig(A'*A);
 for l = 1:length(f);
 p(l) = sqrt(f(l));
 end
 S = diag(p);
 [U,D] = eig(A*A');
 [V,D] = eig(A'*A);
 V = V';
 U
 S

But in answer i am getting only some of its elements same in matrix. And remaining are different. So, please suggest any modification if there or solution..?