orthonormal vectors in generalized eigenvalue problem

Question

Consider the generalized eigenvalue problem $$A\textbf{v}=\lambda B \textbf{v},$$ where $A$ is assumed to be symmetric, nonsingular with distinct eigenvalues and $B$ symmetric and positive definite and let $(\lambda_i,\textbf{v}_i)$ be the $i$-th eigenpair of the problem. I would like to show that the $\textbf{v}_i$ can be chosen such that $$\textbf{v}_i^TB\textbf{v}_j=\delta_{ij}.$$ I have managed to show that for $i \neq j$ the vectors $\textbf{v}_i$ and $\textbf{v}_j$ are orthogonal with respect to $B$ (i.e. $\textbf{v}_i^TB\textbf{v}_j= 0$ for $i \neq j$), but I am not sure how to show that in the case $i=j$ I can have $$\textbf{v}_i^TB\textbf{v}_i=1.$$

Answer 1

I'm assuming your matrices are real.

The generalized eigenvectors $v$ correspond to eigenvectors $w$ of $B^{-1/2} A B^{-1/2}$ by $v = B^{1/2} w$. Since $B^{-1/2} A B^{-1/2}$ is symmetric, it has an orthonormal basis of eigenvectors $w_j$. Then $v_i^T B v_j = w_i^T w_j = \delta_{ij}$.