|
tony1v |
Forum-Newbie
|
|
Beiträge: 1
|
|
|
|
Anmeldedatum: 22.05.09
|
|
|
|
Wohnort: wuppertal
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 20.06.2009, 12:36
Titel: Aitken-Neville Schema
|
|
|
|
|
Aitken-Neville Schema
Das Schema von Aitken-Neville zur Auswertung eines Interpolationspolynoms am Punkt t
bei gegebenen Stützstellen (xi, yi), i = 0, 1, . . . , n lautet
Pi,0(t) := yi, i = 0, . . . , n
Pi,k(t) := Pi,k−1(t) +( t − xi)/(xi − xi−k)*(Pi,k−1(t) − Pi−1,k−1(t)) i >= k, k = 1, . . . , n.
Leiten Sie durch Differenzieren dieser Rekursionsvorschrift ein analoges Schema zur Auswertung
der Ableitung des Interpolationspolynoms am Punkt t her.
Schreiben Sie eine MATLAB-Funktion, die das Schema von Aitken-Neville für die Auswertung
von Interpolationspolynom und dessen Ableitung durchführt. Beachten Sie, dass man das
Schema nicht als Rekursion programmiert, sondern von den gegebenen Werten Pi,0 aus mit
dem Index k fortschreitet. Die Funktion soll folgenden Aufruf besitzen:
function [p,dp] = aitnev(x,y,t)
% Aitken-Neville Schema
% Eingabeparameter:
% x x-Werte der Stuetzstellen als Spaltenvektor
% y y-Werte der Stuetzstellen als Spaltenvektor
% t Stelle an der ausgewertet wird
% Ausgabeparameter:
% p Auswertung des Interpolationspolynoms bei t
% dp Auswertung der Ableitung des Int.pol. bei t
ich habe dieses program programiert aber ich weiss echt nicht ob das schema nicht als Rekursion programiert oder von den gegebenen werte Pi,o aus mit dem Index k fortschreitet was ist mit dp wie kann ich das machen
function p = aitken_neville(x,y,t)
% function p = aitken_neville(x,y,t)
% polynomial interpolant p(t) to data (x,f)
p = y;
for k=1:length(f)-1;
w = (x(1+k:end)-t) ./ (x(1+k:end)-x(1:end-k));
p = w.*p(1:end-1) + (1-w).*p(2:end);
end;
ich bin für jede Antwort sehr dankbar.
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
| RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|