Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function autovv()
- n = input('N° ordem da matriz M: '); //numero de elementos do vetor x0
- m = ones(n); //matriz que deseja encontrar o auto valor e o auto vetor
- for(i = 1 : n)
- for(j = 1: n)
- printf("Informe o valor do elemento A(%d,%d)", i, j);
- m(i,j) = input("");
- end
- end
- clc
- x0 = []; //auto vetor incial Obs.: N° de linhas = n° de colunas da matriz M
- for(i = 1 : n)
- if(i == 2) then
- x0(i,1) = 1;
- else
- x0(i,1) = 0;
- end
- end
- k = 1; //auto valor
- ka = 0; //auto valor anterior
- y = x0; // auto vetor
- ka = 0; // auto valor anterior
- l = 1; // contador de iteração
- epsilon = 10^-6;//tolerância
- while(abs(ka-k) > epsilon) do
- x0 = m * y;
- ka = k;
- k = max(abs(x0));
- printf("Iteração %d\n", l);
- printf("Auto valor: %lf ; Auto vetor: (", k);
- for(i = 1 : n)
- y(i) = (1/k)*x0(i);
- printf("%lf ", y(i));
- end
- printf(")\n");
- l = l + 1;
- end
- endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement