*   >> Leitura Educação Artigos >> education >> colégio e universidade

Método Lagrange e de Newton Dividido Diferença Método - Matlab Scripts

Ambos Newtons Dividido differnce Forma e Lagrange Formulário de criar funções produzem o mesmo polinômio que corresponde dadas entradas

Esses scripts só exibem a resposta final.. Para os scripts que exibem o polinômio criado para produzir a resposta:

>>> Lagrange e de Newton Diferença dividida

Para criar uma imagem visual da função criada:

A Os valores de entrada são: (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); ou >> y = new_div_diff (X, Y, x);

>> plot (x, y);

Para calcular um único ponto:

>> X = 2,5

Lagrange

função y = Lagrange (X, Y, x)% LAGRANGE Avaliar interpolação polinomial usando Lagrange formulário.

% y = LAGRANGE (X, Y, x) retorna y = P (x), onde P é a interpolação polinomial% através dos pontos definidos por X e Y. X pode ser um vector,%, caso em que Y é também um vector, com y (i) igual a P (X (I)) .

n = comprimento (X);

Se o erro n ~ = comprimento (Y) ('X e Y devem ter o mesmo comprimento.'); final

y = 0; % Inicializar sumfor i = 1: n% loop sobre índice soma L = 1; % Inicializar produto para j = [1: 1 i-i + 1: n].

% Loop sobre o índice de produtos L = L * (XX (j)) /(X (i) -X (j)); % Multiplicam próximo fim fator y = y + L * Y (i); % Adicionar próxima termend

Newton Diferença dividida

função y = new_div_diff (X, Y, X) Forma Diferença% Newtons Dividida: Avaliar interpolação polinomiais usando% Newtons Dividido Formulário diferença. y = new_div_diff (X, Y, X) retorna Y = P (x), onde% P é a interpolação polinomial definido pelos pontos em X% e Y. X pode ser um vector, em que Y é caso também um vector, com y (i)% igual a P (x (i))

n = comprimento (X);.


Se n ~ = comprimento (Y) erro ('X e Y must ter o mesmo comprimento ');

final

= Y Y

(1); p = 1; para i = 1:. (n-1) para j = 1: (ni) Y (J) = ( Y (J + 1) - Y (J)) /(X (J + i) - X (J)); terminar para k = i P = P * (XX (i)).; acabar y = y + p * Y

(1);

finais

Outros Scripts Matlab eu criei incluem:

>>> Romberg Integração Numérica - Matlab script

>.

>> Regra de Simpson e trapezoidal Estado de Integração Numérica - Scripts Matlab

>>> Método Newtons de encontrar Roots - Matlab script

>>> Método Bisection de encontrar Roots -

Matlab Script

>>> Método da Secante de encontrar Roots -

Matlab Script

Page   <<       [1] [2] >>

Copyright © 2008 - 2016 Leitura Educação Artigos,https://artigos.nmjjxx.com All rights reserved.