Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function testSVD()
- A = rand(10,5);
- [u,s,v]=svd(A,'econ');
- novDokument=rand(15,1);
- novA=[A;zeros(5,5)];
- novA=[novA,novDokument];
- rekDokum=inv(s)*u'*novDokument(1:10,1);
- mojV=[v;rekDokum'];
- rekBesed=novA(11:15,:)*mojV*inv(s);
- mojU=[u;rekBesed];
- [u1,s1,v1]=svd(novA,'econ');
- u1
- mojU
- v
- mojV
- endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement