Beide Newtons Verdeeld differnce Vorm en Lagrange Vorm van het creëren van functies leveren hetzelfde polynoom die overeenkomt gegeven ingangen
Deze scripts alleen het definitieve antwoord te geven.. Voor de scripts die de polynoom geschapen om het antwoord te produceren geven:
>>> Lagrange en Newton Verdeeld Verschil
Om een visueel beeld van de gecreëerde functie te maken:
De ingang waarden zijn: (1,16) (2,18) (3,21),
>> X = 1: 6;
>> Y = [16 18 21 17 15 12];
>> x = 1: 1/100: 6;
>> y = Lagrange (X, Y, x); of >> y = new_div_diff (X, Y, x);
>> plot (x, y);
Om een enkel punt te berekenen:
>> X = 2,5
functie y = Lagrange (X, Y, x)% LAGRANGE Evalueer interpoleren veelterm met behulp van Lagrange formulier.
% y = LAGRANGE (X, Y, x) geeft y = P (x), waarbij p de interpolerende% polynoom met de door X en Y. X punten kan een vector,% waarbij y een vector met y (i) gelijk aan P (x (i)) .
n = lengte (X);
als n ~ = lengte (Y) fout ('X en Y moeten even lang zijn.'); einde
y = 0; % Initialiseren sumfor i = 1: n% lus over som index L = 1; % Initialiseren product voor j = [1: i-1 i + 1: n].
% Lus over product index L = L * (XX (j)) /(X (i) -X (j)); % Meervoudig volgende factor end y = y + L * Y (i); % Toe te voegen volgende termend
functie y = new_div_diff (X, Y, x)% Newtons Verdeeld Verschil Vorm: Evalueer interpoleren veelterm met behulp% Newtons Verdeeld Verschil Vorm. y = new_div_diff (X, Y, x) geeft y = P (x), waarbij p% van de interpolerende polynoom door de punten gedefinieerd door X en Y.% x kan een vector zijn, waarbij y een vector, met y (i)% gelijk aan P (x (i))
n = lengte (X).
als n ~ = lengte (Y) fout ('X en Y must worden dezelfde lengte ") uiteinde
y = y
(1); p = 1, i = 1:. (n-1) voor j = 1: (ni) Y (j) = ( Y (j + 1) - Y (j)) /(X (j + i) - X (j)); end voor k = i p = p * (xX (i)).; end y = y + p * Y
(1); einde
>>> Romberg Numerieke Integratie - Matlab Script
>.
>> Simpson's Rule en Trapezium Regel van numerieke integratie - Matlab Scripts
>>> Newton methode om Roots - Matlab Script
>>> Bisectie Wijze van het vinden van Roots - Matlab Script
>>> Secant Wijze van het vinden van Roots - Matlab Script