19 lines
415 B
Mathematica
19 lines
415 B
Mathematica
|
% projection de Rayleigh-Ritz
|
||
|
% Données
|
||
|
% A : matrice dont on cherche des couples propres
|
||
|
% V : ensemble de m vecteurs orthonormés
|
||
|
|
||
|
% Résultats
|
||
|
% W : vecteur contenant les approximations des valeurs propres
|
||
|
% V : matrice des vecteurs propres correspondant
|
||
|
function [ W, V ] = rayleigh_ritz_projection( A, V )
|
||
|
|
||
|
H = V'*(A*V);
|
||
|
|
||
|
[VH, DH] = eig(H);
|
||
|
[W, indice] = sort(diag(DH), 'descend');
|
||
|
|
||
|
V = V*VH(:, indice);
|
||
|
|
||
|
end
|