Dividierte Differenzen
Das Schema der dividierten Differenzen wird in der Polynominterpolation zur Bestimmung der Koeffizienten für den Newtonschen Algorithmus verwendet. Für weitere Informationen siehe Dividierte Differenzen.
Die Auswertung der ermittelten Stützstellen erfolgt mithilfe des Horner-Schemas.
Matlab
Siehe Octave.
Octave
Dies ist eine rekursive Implementierung. Für eine Implementierung mit Schleifen siehe Hermiteinterpolation.
function coeff = divdiff(x, y, x0, x1)
if ! exist('x0')
coeff = arrayfun(@(n) divdiff(x, y, 1, n), 1:length(x));
else if x0 == x1
coeff = y(x0);
else
coeff = (divdiff(x, y, x0 + 1, x1) - \
divdiff(x, y, x0, x1 - 1)) / \
(x(x1) - x(x0));
end;
end;
Mithilfe der expliziten Formel für dividierte Differenzen kann man auch ganz auf Schleifen und rekursive Aufrufe verzichten:
function coeff = divdiff(x, y)
divs = repmat(x, length(x), 1);
divs = divs - divs';
divs(logical(eye(size(divs)))) = 1;
coeff = arrayfun(@(n) sum(y(1:n) ./ prod(divs(1:n, 1:n), 1)), 1:length(x));
end
This article is issued from Wikibooks. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.