Octave - Decomposição LU
Publicado por Daniel Moreira dos Santos (última atualização em 10/07/2010)
[ Hits: 21.110 ]
Homepage: http://www.danielmoreira.wordpress.com
Em álgebra linear, a decomposição LU é uma forma de fatoração de uma matriz A não singular como o produto de uma matriz triangular inferior L e uma matriz triangular superior U, isto é A=LU.
Às vezes se deve pré-multiplicar a matriz a ser decomposta por uma matriz de permutação. Esta decomposição se usa em análise numérica para resolver sistemas de equações (mais eficientemente) ou encontrar as matrizes inversas.
function [A, Det, Pivot] = decomLU(n, A) % n, A = ordem e matriz a ser decomposta %A, Det, Pivot = matriz decomposta A = U+L-I, determinante, pivôs for i=1:n Pivot(i) =i; end Det = 1; for j=1:n-1 p = j; Amax = abs(A(j,j)); for k=j+1:n if abs(A(k,j))>Amax Amax = abs(A(k,j)); p =k; end end if p~=j for k=1:n t = A(j,k); A(j,k) = A(p,k); A(p,k) = t; end m = Pivot(j); Pivot(j)=Pivot(p); Pivot(p)=m; Det = -Det; end Det = Det*A(j,j); if abs(A(j,j))~=0 r = 1/A(j,j); for i=j+1:n Mult = A(i,j)*r; A(i,j) = Mult; for k=j+1:n A(i,k) = A(i,k)-Mult*A(j,k); end end end end Det = Det*A(n,n)
Octave - Método de Euler-Modificado para solução de EDO
Função Fatorial Recursiva em Fortran 95
Matriz de Hilbert e resolução de sistemas lineares
Octave - Calcular raiz pelo método de Newton
Nenhum comentário foi encontrado.
Desktop Linux ganha fôlego enquanto distros refinam experiência e segurança
Wayland em alta, IA em debate e novos ventos para distros e devices
Qt 6.10 e Python 3.14 puxam o ritmo do software livre nesta quarta
Atualizações de Apps, Desktop e Kernel agitam o ecossistema Linux nesta terça-feira
Como mudar o nome de dispositivos Bluetooth via linha de comando
Adicionando o repositório backports no Debian 13 Trixie
Meu Amigo não consegue entrar em meu mundo, via tailscale, hamachi e z... (1)
Erro ao fazer remaster usando Penguin eggs (1)
Warcraft II Remastered no Linux? (0)
queria saber como instalar o android power tools igual ao winux no lin... (7)