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