Método das potências

Páginas: 2 (268 palavras) Publicado: 28 de junho de 2012
function [autovalor, autovetor,eixox,eixoy1,eixoy2]=metodo_das_potencias(A,b,n)

t= input('Digite o valor máximo de iterações:');
e= input ('Digite a precisão desejada(na forma decimal): ');
y{1}=b;
z{2}=A*y{1};
p{2}=abs(max(abs(z{2}(:,1))));y{2}=z{2}*(1/p{2});
z{3}=A*y{2};
l{2}=zeros(n,1);
for c=1:n;
l{2}(c,1)=z{3}(c,1)*(1/y{2}(c,1));
end

for r=2:t;z{r+1}=A*y{r};
p{r+1}=abs(max(abs(z{r+1}(1:n,1))));
y{r+1}=z{r+1}*(1/p{r+1});
z{r+2}=A*y{r+1};l{r+1}=zeros(n,1);
for h=1:n;
l{r+1}(h,1)=z{r+2}(h,1)*(1/y{r+1}(h,1));
enderro=zeros(n,1);
for g=1:n;
erro(g,1)=abs((l{r+1}(g,1)-l{r}(g,1))*(1/l{r+1}(g,1)));
end
ifmin(erro(:,1))<=e;
[Q,w]=min(erro(:,1));
autovalor= l{r+1}(w,1);
autovetor=y{r+1};eixoy1(r)=l{r+1}(w,1);
eixoy2(r)=max(eig(A));
eixox(r)=r;
break
end
erro=zeros(n,1);
forg=1:n;
erro(g,1)=abs(l{r+1}(g,1)-l{r}(g,1))*(1/l{r}(g,1));
end
[T,j]=min(erro(:,1));autovalor=l{r+1}(j,1);
autovetor=y{r+1};
eixoy1(r)=l{r+1}(j,1);
eixoy2(r)=max(eig(A));
eixox(r)=r;
end
end
Ler documento completo

Por favor, assinar para o acesso.

Estes textos também podem ser interessantes

  • Métodos da Potência
  • Métodos de compensação da potência reactiva
  • determinação do fator de amortecimento pelo metodo da meia potencia
  • Método de Pontos Interiores Aplicados ao Problema de Fluxo de Potência Ótimo com Restrições de Reserva de Potência...
  • Fluxo De Potencia Pelo Metodo De Gauss Seidel Aula Dia 23 04 2015
  • Potencias
  • potencia
  • Potencia

Seja um membro do Trabalhos Feitos

CADASTRE-SE AGORA!