Algorithmensammlung: Numerik: Dividierte Differenzen

Aus Wikibooks

Algorithmensammlung: Numerik

• Dividierte Differenzen
• Hermiteinterpolation
• Horner-Schema
• Quadratur

[Bearbeiten] 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.

[Bearbeiten] Matlab

Siehe Octave.

[Bearbeiten] 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;