*   >> läser Utbildning artiklar >> education >> college and university

Lagrange Metod och Newton Divided Skillnad Metod - Matlab Scripts

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

Lagrange Omdömen

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

Newton Divided Skillnad 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

Andra Matlab skript jag har skapat inkluderar:

>>> 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

Copyright © 2008 - 2016 läser Utbildning artiklar,https://utbildning.nmjjxx.com All rights reserved.