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 |
Factorizare QR/Householder Matlab
Last Updated: Jun 05 2010 18:54, Started by
no-limit
, Jun 05 2010 16:52
·
0
#1
Posted 05 June 2010 - 16:52
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); endex2) 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
Posted 05 June 2010 - 18:54
fa si tu o proba, ce mare lucru ? sau astepti invitatie speciala ?
|
Anunturi
▶ 0 user(s) are reading this topic
0 members, 0 guests, 0 anonymous users