Advertenties
Als p = n-1, X = V
.
n
Als p < n-1, verwijder dan kolommen p+2, ..., n-1, n uit V
gevormd.
Als p > n-1, voeg dan kolommen n+1, ..., p-1, p+1, toe aan V
matrix X wordt gevormd.
In stap 3 uit deze lijst moeten we erop letten dat de kolom i (i= n+1, n+2, ...,
i
i
p+1) de vector [x
x
... x
1
2
in plaats van een vector gebruiken, dus x = { x
volgorde eenvoudig berekenen { x
omzetten in een vector en het menu COL gebruiken om deze kolommen aan
de matrix V
toe te voegen totdat X voltooid is.
n
Als X klaar is en als de vector y beschikbaar is, dan is de berekening van de
coëfficiëntvector b hetzelfde als de meervoudige lineaire aanpassing (de
vorige matrixtoepassing). We kunnen dus een programma schrijven voor het
berekenen van de polynomiale aanpassing die voordeel heeft van het
programma dat al voor de meervoudige lineaire aanpassing was ontwikkeld.
We moeten de hierboven vermelde stappen 1 tot en met 3 aan het
programma toevoegen.
Het algoritme voor het programma kan dus als volgt worden geschreven:
Voer de vectoren x en y, met dezelfde afmeting, in als lijsten. (Opmerking:
omdat de functie VANDERMONDE een lijst als invoer gebruikt, is het
handiger de gegevens van (x,y) als een lijst in te voeren.) Voer ook de
waarde van p in.
•
Bepaal n = grootte vector x.
•
Gebruik de functie VANDERMONDE om de Vandermonde-matrix V
te genereren voor de ingevoerde lijst x.
•
Als p = n-1, dan
X = V
,
n
Anders, als p < n-1
Verwijder de kolommen p+2, ..., n uit V
(Gebruik een FOR-lus en COL-)
i
] is. Als we een lijst met gegevenswaarden voor x
n
x
1
2
i
i
i
x
... x
}. We kunnen deze lijst dan
1
2
n
zodat X wordt
n
zodat de
n
... x
}, dan kunnen we de
n
om X te vormen
n
Blz. 18-64
n
Advertenties
