Både Newton Divided differnce Form och Lagrange Form skapa funktioner producera samma polynom som matchar givna ingångar Omdömen
Dessa skript visa endast det slutliga svaret.. För skript som visar polynomet skapats för att producera svaret:
>>> Lagrange och Newton Divided Skillnad
För att skapa en visuell bild av den skapade funktion:
ingångsvärden är: (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); eller >> y = new_div_diff (X, Y, x);
>> plot (x, y), Omdömen
För att beräkna en enda punkt:
>> X = 2,5
funktionen y = lagrange (X, Y, x)% LAGRANGE Utvärdera interpole polynom med hjälp av Lagrange form.
% y = LAGRANGE (X, Y, x) returnerar y = P (x), där P är interpole% polynom genom punkterna som definieras av X och Y. x kan vara en vektor,% i vilket fall y är också en vektor, med y (i) lika med P (x (i)) .
n = längd (X), Omdömen
om n ~ = längd (Y) fel ("X och Y måste vara lika långa. '); slutet Omdömen
y = 0; % Initiera sumfor i = 1: n% slinga över summa index L = 1; % Initiera produkt för j = [1: i-1 i + 1: n].
% Slinga över produktindex L = L * (XX (j)) /(X (i) -X (j)); % Multiplicera Nästa faktor ände y = y + L * y (i); % Lägga nästa termend Omdömen
funktionen y = new_div_diff (X, Y, x)% Newton Divided Skillnad Form: Utvärdera interpole polynoma använder% Newton Divided Skillnad Form. y = new_div_diff (X, Y, x) returnerar y = P (x),% där P är den interpole polynom genom punkterna definieras av X% och Y. x kan vara en vektor, i vilket fall y är också en vektor, med y (i)% är lika med P (x (i)) katalog
n = längd (X),.
Omdömen
om n ~ = längd (Y) fel ("X och Y måste ha samma längd '); slutet
y = Y
(1); p = 1; för i = 1:. (n-1) för j = 1: (ni) Y (j) = ( Y (j + 1) - Y (j)) /(X (j + i) - X (j)); avsluta för k = i p = p * (XX (i)).; avsluta y = y + p * Y
(1), end Omdömen
>>> Romberg Numerisk integration - Matlab Script Omdömen
>.
>> Simpsons regel och trapetsregeln för numerisk Integration - Matlab skript
>>> Newtons metod för att hitta Roots - Matlab Script Omdömen
>>> bisektionsmetoden att finna Roots - Matlab Script Omdömen
>>> sekantmetoden att finna Roots - Matlab Script Omdömen