Método da Interpolação de Newton em Scilab

-
Versão do Scilab: 5.5.2.

clc
clear

 function [polinomio] = polinomioNewton(vetX, vetY)
    n  = size(vetX, 1) - 1
    x = poly(0,"x")
    cpy = vetY
    coeficientes = []
    polinomio = vetY(1)
    for i=1:n
        diferences = []
        for j=1:(n - i+1)
            diference = (cpy(j+1) - cpy(j) ) / (vetX(j+1) - vetX(j))
            diferences = [diferences ; diference]
        end
        coeficientes(i) = diferences(1);
        cpy = diferences
    end
    cpy = x -vetX
    for i=1: n 
        polinomio = polinomio + prod(cpy(1:i)) * coeficientes(i)
    end
endfunction

 //ENTRADA DE DADOS
vetX = [0; 1822; 2994; 3640; 8848]
vetY = [100; 94; 91; 87; 71]

 polinomio = polinomioNewton(vetX, vetY)
disp("Polinomio de Newton: ")
disp(polinomio)

 //calcula o valor de um ponto no polinomio
result = horner(polinomio,2892)
disp(result)

Comentários

Postar um comentário

Obrigado pelo contato, retornarei o mais breve possível!

Postagens mais visitadas deste blog

Método da Interpolação por Eliminação Gaussiana em Scilab

-
Versão do Scilab: 5.5.2 Observação: Só calcula equação linear e quadrática; clear clc //Eliminação Gaussiana //Variavel matriz é a matriz incompleta //Variavel matrizResposta é a matriz que contem as respostas dos sistemas lineares function  matrizCoeficientes= eliminacaoGauss(matriz, matrizResposta)     matrizCompleta=[matriz matrizResposta];     n=size(matriz,1);     disp("Resolução do Sistema")     disp("Matriz Ampliada = ")     disp(matrizCompleta)     for j=1:(n-1)         for i=(j+1):n            mij=matrizCompleta(i,j)/matrizCompleta(j,j)            matrizCompleta(i,:)=matrizCompleta(i,:)-mij*matrizCompleta(j,:)         end     end     disp("Matriz Escalonada")     disp(matrizCompleta)     //Fazer Retrosubstituição     matrizCoeficientes(n)=matr...
Leia mais

Cálculo do Erro para Interpolação em Scilab

-
Versão do Scilab: 5.5.2 clear clc function [orderValues] = getOrder(lastOrder, vectorX, orderNumber)     for i = 1 : size(lastOrder, 2) - 1 // size 1 = quantas linhas, size 2 = quantas colunas        orderValues(1,i) = (lastOrder(1,i+1) - lastOrder(1,i)) / (vectorX(1,orderNumber+i-1) - vectorX(1,i))     end endfunction function [f, polyF] = newtonInterpol(vectorY, vectorX)     lastOrder = vectorY     polyF = poly([vectorY(1,1)],"x", "coeff")     for i = 2 : size(vectorX, 2) // size 1 = quantas linhas, size 2 = quantas colunas         lastOrder = getOrder(lastOrder, vectorX, i)         polyA = poly([-vectorX(1,1) 1],"x","coeff")         for j = 2 : i-1             polyA = polyA * poly([-vectorX(1,j) 1],"x","coeff")         end         polyA = polyA * lastOrder(1,1) ...
Leia mais

Método de Interpolação de Lagrange em Scilab

-
Versão do Scilab : 5.5.2 clc  clear //funcao que retorna o polinomio de lagrange function y=polinomioLagrange(tab); // calcula o polinomio interpolador pelo metodo de lagrange;     [l,n]=size(tab);     x=tab(1,:);     for i=1:n           if i>size(x) then            error(" i incompativel com dimensao da tabela");         end;         [l,n]=size(x);         for k=1:i-1             r(k)=x(k)         end;         for k=i:n-1              r(k)=x(k+1)         end;         P=1;         for k=1:(n-1)              P=P*(x(i)-r(k))         end;         if P==0 then    ...
Leia mais