Sign In
Create Account
Back to top |
Chirurgia spinală minim invazivă
Chirurgia spinală minim invazivă oferă pacienților oportunitatea unui tratament eficient, permițându-le o recuperare ultra rapidă și nu în ultimul rând minimizând leziunile induse chirurgical. Echipa noastră utilizează un spectru larg de tehnici minim invazive, din care enumerăm câteva: endoscopia cu variantele ei (transnazală, transtoracică, transmusculară, etc), microscopul operator, abordurile trans tubulare și nu în ultimul rând infiltrațiile la toate nivelurile coloanei vertebrale. www.neurohope.ro |
Factorizare QR/Householder Matlab
Last Updated: Jun 05 2010 18:54, Started by
no-limit
, Jun 05 2010 16:52
·
0
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
-
- Change Theme
- English
- Mark Community Read
- Termene & conditii
- Confidentialitate
- Consimtamant
- Cookies
- Help
© 2001-2024 Softpedia. All rights reserved. Softpedia® and the Softpedia logo are registered trademarks of SoftNews Net SRL.
Ora implicita a Forumului Softpedia este ora standard a Romaniei (GMT+2). Mai multe informatii →
⚠ Postarile afisate pe Forumul Softpedia reprezinta o opinie subiectiva a membrilor care le-au publicat si nu a SoftNews Net SRL.