Calculul valorilor si vectorilor proprii

5.3. ALGORITMUL DVS 409
5. Dacă opt 2 = ′ da ′ şi β ≠ 0 atunci
V = Ā(1,1)
f
6. Return
2. % Reducerea la forma bidiagonală
[f,g,U,V ] = JQc(A,opt 1 ,opt 2 )
3. % Faza iterativă
1. p = 0, q = 0
2. C^at timp q < n−1
1. Pentru i = 1 : n−q −1
1. Dacă |g i | ≤ tol(|f i |+|f i+1 |) atunci
g i ← 0
2. % Determinarea parametrului q
1. C^at timp g(n−q −1) = 0
q ← q +1
2. Dacă q = n−1 atunci break
3. % Terminarea fazei iterative
1. Dacă q = n−1 atunci break
4. % Determinarea parametrului p
1. p = n−q −1
2. C^at timp g(p) = 0
1. p ← p−1
2. Dacă p = 0 atunci break
5. k = p+1, l = n−q
6. % Tratarea unui zero diagonal (dacă există)
1. j = 0
2. Pentru i = l−1 : −1 : k
1. Dacă f i = 0 atunci
1. j = i−k +1
2. break
3. Dacă j > 0
1. [f(k : l),g(k : l−1),Z] =
= TZDc(f(k : l),g(k : l−1),j,I l−k+1 ,opt 1 )
2. Dacă opt 1 = ′ da ′ atunci
U(:,k : l) ← U(:,k : l)Z
altfel
1. % Iteraţia curentă
Dacă k < l−1 atunci
[f(k : l),g(k : l−1),U(:,k : l),V(:,k : l)] =
= IT DVSc(f(k : l),g(k : l−1),U(:,k : l),
V(:,k : l),opt 1 ,opt 2 )
altfel