Algoritmo di minimizzazione
x0 ← ph ← (b - x0)n
Stima iniziale
Per k = 1, 2... diamo una stima iniziale y0 del vettore y = (y1, y2, ..., yn), ponendo yk = y0,k. Quindi calcoliamo gli incrementi lJ e mk.
Incrementi di Hooke e Jeeves
Gli incrementi lk sono quelli di Hooke e Jeeves, mentre mk sono quelli di coordinate discendenti, che sono definiti nel modo seguente:
- Determiniamo un incremento l = lk nel vettore yk, tale da minimizzare il vettore yk + l
- mk ← yk - min θ > 0
Fine algoritmo
Programma di riempimento e disegno
begin
- input (rett, rect)
- fill (rett, rect)
- draw (rett, g, g)
- for k=0 to n-1 do
- continue
- end
Inversione della matrice triangolare superiore
Si desidera invertire la seguente matrice triangolare superiore:
U = 0 1 2 0 1 2 3 0 0 2 2 3 0 1 0 0
A tal fine, si devono effettuare le seguenti iterazioni:
- Z(l-i,l-i) = U(i,i)
- Z(I, i-1:l) = -1/l(i-j)
- Z(1:i-1,1:i-1) = Z(1:i-1,1:l)
- Z(l:l)+U(i:i-1,1:l)
Fine programma!
Calcolo del vettore kn
kn = b - Ank ∈ Rn, b - A1 k = ...[b - A0]
I vettori (o scalari) di detti vettori u, v ∈ Rn; dk (che vale 0 o 1), e ovale, 1.