[MATLAB] Calcul valeurs propres matrices
Résolu
snakehill
Messages postés
39
Date d'inscription
Statut
Membre
Dernière intervention
-
555 -
555 -
Bonjour,
toujours dans le probleme du tambour , qui devient de plus en plus un problème ch**nt.
Il s'agit de calculer les valeur propres d'une matrice 300*300 avec la méthode des puissances itérées. J'ai donc conçu un petit algorithme qui fonctionne ... Mais juste si l'on ne dépasse pas la dimension 10*10 , ce qui est vraimennt étrange ... donc voici le code ... Si quelqu'un peut m'éclairer , ou juste le signaler où se trouve l'érreur ce serait cool ^^
toujours dans le probleme du tambour , qui devient de plus en plus un problème ch**nt.
Il s'agit de calculer les valeur propres d'une matrice 300*300 avec la méthode des puissances itérées. J'ai donc conçu un petit algorithme qui fonctionne ... Mais juste si l'on ne dépasse pas la dimension 10*10 , ce qui est vraimennt étrange ... donc voici le code ... Si quelqu'un peut m'éclairer , ou juste le signaler où se trouve l'érreur ce serait cool ^^
function [BD,mat,VP] = vp(Tambour) taille=50; BD=zeros(taille,1); % On prend un Y quelconque Y=eye(taille,1); % On initie mat mat=Tambour; %Pour tout le tambour for j=1:1:taille %On lance la boucle des puissances itérées for i=1:1:15 el = norm(Y); X = Y/el; Y = Tambour*X; end %Récupération de la valeur propre lambda = X'*Y; %Stockage de la valeur propre dans le vecteur BD en sortie BD(j)=lambda; % On cherche le vecteur propre en utilisant jacobi (A-lambda*Id)*VP=B avec % B un vecteur nul B = zeros(taille,1); mat = mat -(lambda*eye(taille,taille)); VP = seidel(mat,B,0.01); % On enlève la matrice du vecteur propre récupéré (VP) Tambour = Tambour - lambda*((VP*VP')/(VP'*VP)); % On doit stocker les valeurs propres et les vecteurs propres %fprintf('Valeur propre %d /n', lambda) %fprintf('vecteur propre :/n') %VP end BD end
A voir également:
- Méthode de la puissance itérée matlab
- Test puissance pc - Guide
- Puissance en c++ - Forum C++
- Puissance en c ✓ - Forum C
- &Quot;pow " pour calculer la puissance en C - Forum C
- Puissance signal wifi - Guide