Suponer es una matriz definida positiva simétrica de dimensiones. Matriz es una matriz de valor real de rango completo con estrictamente menor que , es decir,
A través de los experimentos de Montecarlo, me di cuenta de que
Me pregunto de dónde viene esta desigualdad.
Para más ilustración, se puede ejecutar el siguiente código corto en Matlab.
n=30;
m=29;
c=0;
for i=1:1000;
Q = randn(n,n);
eigen_mean = 0.1; %can be made anything,
A = Q' * diag(abs(eigen_mean+randn(n,1))) * Q; %A random symmetric positive definite
B=randn*randn(m,n);
if max(eig(A))< max(eig(B * A * pinv(B)))
c = c +1 ;
end
end
c
Cada vez se devuelve zeor, ya que aparentemente nunca es más pequeño que
Dejar . Tenga en cuenta que y son simétrico matrices y, por tanto, su producto es diagonalizable y sólo tiene valores propios.
Más precisamente, es similar a la siguiente matriz simétrica
.
Para cada vector , dónde .
Entonces dónde y hemos terminado.
fimpellizzeri
Amir Amini
fimpellizzeri
Amir Amini