1 2 3 4 5 6 7 8 9 10 11 12 13
//332 clear; close; clc; A=[1 -2;3 -1]; disp(A,'A='); [U diag1 V]=svd(A); Q=U*V'; S=[2 1;1 3]; disp(Q,'Q='); disp(S,'S=') disp(S'*Q,'A=S''Q=') //end