Exemplo

-

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