Matriz de Hilbert e resolução de sistemas lineares

Publicado por Daniel Moreira dos Santos (última atualização em 06/07/2010)

[ Hits: 10.485 ]

Homepage: http://www.danielmoreira.wordpress.com

Download sistema_hilbert.m




Uma matriz de Hilbert é uma matriz quadrada com a seguinte forma:

H_{ij} = \frac{1}{i+j-1}

Aqui vamos mostrar como resolver um sistema linear HX=b utilizando o octave, onde H é uma matriz de Hilbert.

  



Esconder código-fonte

%Resolve o sistema linear AX=b onde A é a matriz de Hilbert de ordem n
%Isso não é portugol e sim octave, mas não existe essa categoria de contribuição de scripts.

function X=sistema_hilbert(n)

H=hilb(n);
 %cria H como uma matriz de Hilbert de ordem n
b(n)=zeros;
b=b';
for i=1 : n
  for j=1 : n
     b(i)=b(i)+H(i,j); %calcula b
  end
end

X=H\b;  %resolve o sistema (equivalente a X=inv(H)*b) e retorna o vetor solução X

Scripts recomendados

Numero Primo

Octave - Sistemas lineares por Gauss-Seidel

Octave - Decomposição LU

Octave - Método de Runge-Kutta

Teste de hardware e outros


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor RedeHost.
Linux banner
Linux banner
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts