Take the 2-minute tour ×
MathOverflow is a question and answer site for professional mathematicians. It's 100% free, no registration required.

I have the equation $\Sigma_k(M_k{p_k})V=EV$, where the $M_k$ are n*n real Hermitian matrices, $V$ is a n*n eigenvector matrix, $E$ a dim-n energy eigenvector and the $p_k$ scalar parameters. The $M_k$ are known exactly, and I get m energies $e$ from experiment, where (usually) $k\lt{m}\lt{n}$.
I now want to find the set of $p_k$ such that $\Sigma_m(e_m-E_m)^2$ is minimized (where $E_m$ are the "theoretical" eigenvalues following from the above equation, and $e_m$ my experimental ones - knowing beforehand which $E_m$ and $e_m$ "belong" together is another can of worms).
Since 30 years I solve this problem with INVEIX, a FORTRAN kludge which starts from an initial set $p_k$ and converges to a miminum by a steepest gradient algorithm. Kludge or not, it works splendid since 30 years :-)
Still, I'm curious: Can this minimax problem be solved more analytically, maybe even as closed form? I'm very reluctant to differentiate this matrix equation after the $p_k$...

share|improve this question
    
What does "*" mean here? If just multiplication, it would be clearer to omit the symbol. –  David Ketcheson Dec 6 '12 at 8:35
    
It is just multiplication. Edited. –  Hauke Reddmann Dec 7 '12 at 12:30

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.