Octave - Calcular raiz pelo método de Newton
Publicado por Daniel Moreira dos Santos (última atualização em 15/07/2010)
[ Hits: 22.930 ]
Homepage: http://www.danielmoreira.wordpress.com
Em análise numérica, o método de Newton (ou método de Newton-Raphson) tem o objetivo de estimar as raízes de uma função. Para isso, toma-se um ponto qualquer da função, calcula-se a equação da tangente (derivada) da função nesse ponto, calcula-se o intercepto da tangente ao eixo das abcissas, calcula-se o valor da função nesse ponto, e repete-se o processo, que deve tender a uma das raízes da função rapidamente, ou não tender a nada, deixando isso claro logo.
function [Raiz,Iter,Erro] = Newton_Raiz(funcao,x0,Toler,IterMax); %parametros de entrada: funcao,x0,Toler,IterMax -> funcao, valor inicial, %tolerancia e num. max. de iteracoes %parametros de saida: Raiz,Iter,Erro-> raiz, num. gasto de iteracoes e %cond. de erro x = x0; Fx = eval(funcao); DFx = eval(diff(funcao)); Iter = 0; disp(' Iter x DFx Fx DeltaX'); disp([Iter x DFx Fx ]); while (1) DeltaX = -Fx/DFx; x = x+DeltaX; Fx = eval(funcao); DFx = eval(diff(funcao)); Iter = Iter+1; disp([Iter x DFx Fx DeltaX]); if(abs(DeltaX)<Toler && abs(Fx)<Toler)||abs(DFx)==0 || Iter>= IterMax break; end end Raiz = x; if abs(Fx)<Toler Erro = 0; else Erro =1; end y = []; for x=0.001:0.001:10 y = [y,eval(funcao)]; end x = 0.001:0.001:10; plot(x,y); end
Octave - Sistemas lineares por Gauss-Seidel
Octave - Método de Runge-Kutta
Octave - Sistemas lineares por Jacobi
Octave - Método de Euler-Modificado para solução de EDO
Nenhum coment�rio foi encontrado.
Conciliando o uso da ZRAM e SWAP em disco na sua máquina
Servidor de Backup com Ubuntu Server 24.04 LTS, RAID e Duplicati (Dell PowerEdge T420)
Visualizar câmeras IP ONVIF no Linux sem necessidade de instalar aplicativos
Realizar overclock no Miyoo Mini (plus ou normal)
Otimização de memória para máquinas modestas
Unbuntu não atualiza o firmware (1)
linux mint reconhece microfone de lapela como fone de ouvido sem micro... (0)
Dúvidas sobre a originalidade de conteúdos online (10)
Erro de interface de Rede no Virt Manager dentro Debian 13 KDE (12)