Molemmat Newtonia Jaettu differnce muoto ja Lagrangen muoto luoda toimintoja tuottavat saman polynomi, joka vastaa annetuilla panoksilla.
Nämä skriptit näyttää vain lopullinen vastaus. Sillä skriptejä, jotka näyttävät polynomi luotu tuottamaan vastauksen:
>>> Lagrangen ja Newton Divided Ero
Voit luoda visuaalinen kuva perustettuun:
tulo arvot ovat: (1,16), (2,18), (3,21),
>> X = 1: 6;
>> Y = [16 18 21 17 15 12];
>> x = 1: 1/100: 6;
>> y = Lagrangen (X, Y, x); tai >> y = new_div_diff (X, Y, X);
>> tontti (x, y);
Jos haluat laskea yhden pisteen:
>> X = 2.
5
funktio y = Lagrangen (X, Y, X)% LAGRANGE arvioida interpoloimalla polynomi käyttäen Lagrangen lomakkeella.% y = LAGRANGE (X, Y, X) palauttaa y = P (x), jossa P on interpoloiva% polynomi pisteiden kautta määritelty X: llä ja Y: x voi olla vektori,%: jolloin y on myös vektori, y (i) on yhtä suuri P (x (i)) .
n = pituus (X);
jos n ~ = pituus (Y) virhe ("X ja Y tulee olla sama pituus."); päähän
y = 0; % Alusta sumfor i = 1: n% silmukka yli summa indeksi L = 1; % Alustaa tuote J = [1: i-1 i + 1: n]% loop yli tuotteen indeksi L = L.
* (Xx (J)) /(X (i) -X (j)); % Multiply seuraava tekijä lopussa y = y + L * Y (i); % Lisätä seuraava termend
funktio y = new_div_diff (X, Y, X)% Newtonia Divided Ero muoto: Arvioi interpoloimalla polynomi käyttämällä% Newtonia Divided Ero muoto. y = new_div_diff (X, Y, x) palauttaa y = P (x),% jossa P on interpoloiva polynomi kautta pisteiden määritelty X% ja Y x voi olla vektori, jolloin y on myös vektori, Y (i)% vastaa P (x (i)).
n = pituus (X);
jos n ~ = pituus (Y) virhe (X ja Y rypäleen olla samanpituisia.
"); päähän
y = Y
(1); p = 1, i = 1: (n-1), kun j = 1: (ni) Y (j) = ( Y (j + 1) - Y (j)) /(X (j + i) - X (j)); end k = i p = p. * (xx (i)); end y = y + s. * Y
(1); lopussa
>>> Romberg numeerinen integrointi - Matlab Script
> >> Simpsonin sääntö ja trapetsoidisäännöllä Numerical Integration - Matlab Käsikirjoitukset
>>> Newtonin menetelmä löytää Roots - Matlab Script
>>> bisection menetelmä löytää Roots - Matlab Script
>>> sekanttimenetelmää löytää Roots - Matlab Script
com