Octave - Sistemas lineares por Jacobi
Publicado por Daniel Moreira dos Santos (última atualização em 14/07/2010)
[ Hits: 13.455 ]
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 da Secante
Matriz de Hilbert e resolução de sistemas lineares
Crivo de Eratóstenes Simples em Fortran 95
Octave - Método de Runge-Kutta
Octave - Método de Gauss-Legendre
Nenhum comentário foi encontrado.
Melhorando o tempo de boot do Fedora e outras distribuições
Como instalar as extensões Dash To Dock e Hide Top Bar no Gnome 45/46
E a guerra contra bots continua
Tradução do artigo do filósofo Gottfried Wilhelm Leibniz sobre o sistema binário
Conheça o firewall OpenGFW, uma implementação do (Great Firewall of China).
Instalando o FreeOffice no LMDE 6
Anki: Remover Tags de Estilo HTML de Todas as Cartas
Colocando uma opção de redimensionamento de imagem no menu de contexto do KDE
Não consigo acessar os modos de desempenho no mint 2.3 (9)
Problema com alias usando locate (0)
[Shell Script] Script para desinstalar pacotes desnecessários no OpenSuse
[Shell Script] Script para criar certificados de forma automatizada no OpenVpn
[Shell Script] Conversor de vídeo com opção de legenda
[C/C++] BRT - Bulk Renaming Tool
[Shell Script] Criação de Usuarios , Grupo e instalação do servidor de arquivos samba