Fituici MNCPD
Implicatii in rez. sist. de ec. Liniare Compatibil det (sist patrate): Ax=b, A R nesing. isel_a.m n=size(a); a1=a;b1=b; for k=1:n-1 sig=0; for i=k:n sig=sig+a(i,k)^2; end; sig=sign(a(k,k))*sqrt(sig); for i=1:n if i<=k-1 v(i)=0; elseif i==k v(i)=sig+a(k,k); else v(i)=a(i,k); end; end; beta=sig*(sig+a(k,k)); q=eye(n)-(v'*v)/beta; a=q*a; b=q*b; end; b(n)=b(n)/a(n,n); for i=n-1:-1:1 s=0; for j=i+1:n s=s+a(i,j)*b(j); end; b(i)=(b(i)-s)/a(i,i); end; Ax=b, A R ,m>n...
Domeniu: Calculatoare