Octave - Sistemas lineares por Jacobi
Publicado por Daniel Moreira dos Santos (última atualização em 14/07/2010)
[ Hits: 13.534 ]
Homepage: http://www.danielmoreira.wordpress.com
O método de Jacobi trata-se dum algoritmo para determinar a solução de um sistema de equações lineares com os maiores valores absolutos em cada linha e coluna dominados pelo elemento da sua diagonal. Trata-se duma versão simplificada do algoritmo de valores próprios de Jacobi.
function [x,Iter,CondErro] = Jacobi(n,A,b,Toler,IterMax) %n,A,b,Toler,IterMax = ordem,matriz,vetor independente,tolerancia,numero %maximo de iteracoes %x,Iter,CondErro = vetor solucao,numero de iteracoes e condicao de erro for i=1:n r=1/A(i,i); for j=1:n if i~=j A(i,j)=A(i,j)*r; end end b(i)=b(i)*r; x(i)=b(i); end Iter =0; fprintf('\nIter = %d ' ,Iter); x while 1 Iter = Iter+1; for i=1:n Soma=0; for j=1:n if i~=j Soma=Soma+A(i,j)*x(j); end end v(i)=b(i)-Soma; end NormaNum =0; NormaDen=0; for i=1:n t=abs(v(i)-x(i)); if t>NormaNum NormaNum =t; end if abs(v(i))>NormaDen NormaDen = abs(v(i)); end x(i)=v(i); end NormaRel = NormaNum/NormaDen; fprintf('\nIter = %d ' ,Iter); x fprintf('NormaRel = %f \n',NormaRel); if NormaRel<=Toler | Iter>=IterMax break; end end if NormaRel<=Toler CondErro =0; else CondErro =1; end
Octave - Calcular raiz pelo método de Newton
Octave - Sistemas lineares por Gauss-Seidel
Octave - Método de Euler-Melhorado para solução de EDO
Octave - Método de Gauss-Legendre
Nenhum comentário foi encontrado.
Como escolher o melhor escalonador de CPU para melhorar o desempenho da máquina
Curiosidade sobre DOOM Guy e Isabelle de Animal Crossing
Inicializando servidor Ubuntu na AWS e rodando apache em Container
Conheça o Octopi, outro frontend para o Pacman com acesso ao ARU (Arch Linux e derivados)
Terminal transparente no Debian 12 com interface i3wm usando Xfce4-Terminal e Compton
Como compilar o DOOM RETRO no Linux
Stacer - um otimizador gráfico básico de Linux
Corrigindo o erro de "WARNING: Possibly missing firmware for module" no Arch Linux
O que siginifica esse texto escrito em binário! (5)
Ubuntu simplesmente morreu (1)
horário adiantado no sistema (9)