Calculul valorilor si vectorilor proprii

270 CAPITOLUL 4. VALORI ŞI VECTORI PROPRII
3. % Terminarea normală a fazei iterative
1. Dacă q ≥ n−2 atunci break
4. % Terminarea prin eşec a algoritmului
1. Dacă cont it > 30 atunci
1. Tipăreşte ’S-au consumat 30 iteraţii QR pentru
evidenţierea unui bloc diagonal fără a se atinge
acest obiectiv. Este posibil ca, pentru aceste date
de intrare, algoritmul QR să nu fie convergent.’
2. Return
5. % Determinarea lui p
1. p = n−q −1
2. C^at timp a p+1,p ≠ 0
1. p = p−1
2. Dacă p = 0 atunci break
6. % Iteraţia curentă
1. k = p+1, l = n−q
2. w = VD2(A(k:l,k:l))
3. % Calculul deplasării implicite modificate
1. Dacă cont it = 10 sau cont it = 20 atunci
1. s = 1.5(|a l,l−1 |+|a l−1,l−2 |)
2. p = (|a l,l−1 |+|a l−1,l−2 |) 2
⎡
3. w = ⎣ a2 kk +a k,k+1a k+1,k −sa kk +p
a k+1,k (a kk +a k+1,k+1 −s) ⎦
a k+1,k a k+2,k+1
4. [A(k : l,k : l),V,b] = IT QR2(A(k:l,k:l),w)
5. t = 3
6. Pentru i = 1 : l−k
1. Dacă i = l−k atunci t = 2
2. r = min(p+i+2,l)
3. Dacă k > 1 atunci
1. A(1 : p,p+i : r) =
= Hrd(A(1 : p,p+i : r),V(1:t,i),b i )
4. Dacă l < n atunci
1. A(p+i : r,l+1 : n) =
= Hrs(V(1:t,i),b i ,A(p+i : r,l+1 : n))
7. cont it = cont it+1
7. Dacă opt =’da’ atunci
1. t = 3
2. Pentru i = 1 : l−k
1. Dacă i = l−k atunci t = 2
2. r = min(p+i+2,l)
3. Q(:, p+i : r) = Hrd(Q(:, p+i : r),V(1:t,i),b i ))
⎤