Matriz de Hilbert e resolução de sistemas lineares

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

[ Hits: 12.397 ]

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

Octave - Calcular raiz pelo método Regula-Falsi

Teste de hardware e outros

Octave - Decomposição LU

Script em Bath "DOS" para Upload e Download de arquivos usando o SSH

Octave - Calcular raiz pelo método de Newton


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor RedeHost.
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts