Daarna krijgt u de vector van de coëfficiënten uit b = (X
vector y = [y
y
1
In hoofdstuk 10 hebben de Vandermonde-matrix die overeenkomt met een
vector x = [x
x
1
matrix X die interessant is voor de polynomiale aanpassing, maar heeft alleen
n in plaats van (p+1) kolommen.
We kunnen de functie VANDERMONDE gebruiken om de matrix X te maken
als we ons aan de volgende regels houden:
Als p = n-1, X = V
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, ...,
p+1) de vector [x
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.
T
... y
]
is.
2
n
... x
] gedefinieerd. De Vandermonde-matrix lijkt op de
2
m
.
n
i
i
i
x
... x
] is. Als we een lijst met gegevenswaarden voor x
1
2
n
x
1
i
i
i
x
... x
}. We kunnen deze lijst dan
1
2
n
T
-1
T
⋅X)
⋅X
⋅y, waarbij y
zodat X wordt
n
zodat de
n
... x
}, dan kunnen we de
2
n
n
Blz. 18-59