Jump to content

SUBIECTE NOI
« 1 / 5 »
RSS
Incalzire in pardoseala etapizata

Suprataxa card energie?!

Cum era nivelul de trai cam din a...

probleme cu ochelarii
 Impozite pe proprietati de anul v...

teava rezistenta panou apa calda

Acces in Curte din Drum National

Sub mobila de bucatarie si sub fr...
 Rezultat RMN

Numar circuite IPAT si prindere t...

Pareri brgimportchina.ro - teapa ...

Lucruri inaintea vremurilor lor
 Discuții despre TVR Sport HD.

Cost abonament clinica privata

Tremura toata, dar nu de la ro...

Renault Android
 

Factorizare QR/Householder Matlab

  • Please log in to reply
1 reply to this topic

#1
no-limit

no-limit

    Member

  • Grup: Members
  • Posts: 800
  • Înscris: 12.12.2007
Buna, am de realizat si eu un program cu interfata grafica in Matlab care sa realizeze factorizarea QR a matricilor folosind metoda Householder.
Pentru descompunerea QR am gasit un exemplu functional:

A=[1 2 3;4 5 6;7 8 9][Q,R]=qr(A)


Insa pentru Householder nu am nicio idee cum as putea continua programul ,am gasit 2 algoritmi insa nu stiu daca sunt buni
ex1)
function [A] = QR_HOUSE(A);
[n,m] = size(A);
for k = 1:min(n-1,m) v(k:n,1) = HOUSEHOLDER(A(k:n,k));
	   A(k:n,k:m) = HOUSEHOLDER_MULT(A(k:n,k:m),v(k:n,1));
	 A(k+1:n,k) = v(k+1:n,1);
end
ex2)
function a = househ(a)
% Turn matrix a into upper- hessenberg with same eigenvalues[m,n]=size(a);
if m~=n, error('Not a square matrix!');
 endfor i=1:n-2		   c=zeros(n,1);
		   c(i+1:n)=a(i+1:m,i);
% set up c		 q=norm?;
				b=zeros(n,1);
		   if a(i+1,i)>=0,				 b(i+1)=-q;
			  else					b(i+1)=q;
			   end			 w=(c-B)/norm(c-B);
 Q=eye(n)-2*w*w';
				a=Q*a*Q;
				a(i+2:n,i)=0;
 %let's correct the roundoff errorend
 

V-as ruga sa imi spuneti daca vreunul dintre cei 2 algoritmi householder este bun, iar in caz contrar unde as putea gasi un exemplu ok si functional cu aceasta metoda de rezolvare .

Edited by MarianG, 14 March 2019 - 02:14.


#2
soarce

soarce

    Active Member

  • Grup: Members
  • Posts: 1,467
  • Înscris: 02.10.2008
fa si tu o proba, ce mare lucru ? sau astepti invitatie speciala ?

Anunturi

Second Opinion Second Opinion

Folosind serviciul second opinion ne puteți trimite RMN-uri, CT -uri, angiografii, fișiere .pdf, documente medicale.

Astfel vă vom putea da o opinie neurochirurgicală, fără ca aceasta să poată înlocui un consult de specialitate. Răspunsurile vor fi date prin e-mail în cel mai scurt timp posibil (de obicei în mai putin de 24 de ore, dar nu mai mult de 48 de ore). Second opinion – Neurohope este un serviciu gratuit.

www.neurohope.ro

0 user(s) are reading this topic

0 members, 0 guests, 0 anonymous users

Forumul Softpedia foloseste "cookies" pentru a imbunatati experienta utilizatorilor Accept
Pentru detalii si optiuni legate de cookies si datele personale, consultati Politica de utilizare cookies si Politica de confidentialitate