You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
20 SEMINAR 1. CALCUL MATRICEAL ELEMENTAR32: vi=0;33: for i=1:n34: sigma=0;35: for j=i:n36: sigma =sigma + abs(U(i,j));37: end38: if vi < sigma39: vi = sigma;40: end41: end42: vF=0;43: for i=1:n44: for j=i:n45: vF=vF + U(i,j)*U(i,j);46: end47: end48: vF=sqrt(vF);1: function [ Q ]=gramschmidt(B)2:3: %-------------------------------------------------------------------------4: % Algoritmul 1.55: % Functia calculeaza p vectori ortogonali care genereaza acelasi subspatiu6: % ca p vectori liniar independenti dati.7: % Vectorii dati sunt memorati in coloanele matricei B, iar vectorii8: % calculati sunt memorati in coloanele matricei Q.9: % Algoritmul este cunoscut sub numele de Algoritmul Gram-Schmidt10: % Apelul: [ Q ]=gramschmidt(B)11: %12: % Buta Valentin, aprilie, 200613: %-------------------------------------------------------------------------14:15: [ n,p ]=size(B);16: for k=1:n17: Q(k,1)=B(k,1);18: end19: for j=1:p-120: beta(j)=0;21: for k=1:n22: beta(j)=beta(j)+Q(k,j)*Q(k,j);23: Q(k,j+1)=B(k,j+1);24: end25: for i=1:j26: alfa(i,j)=0;27: for k=1:n28: alfa(i,j)=alfa(i,j)+Q(k,i)*B(k,j+1);29: end30: alfa(i,j)=-alfa(i,j)/beta(i);31: for k=1:n32: Q(k,j+1)=Q(k,j+1)+alfa(i,j)*Q(k,i);