Jacobi Method For Complex Hermitian Matrices
   HOME
*





Jacobi Method For Complex Hermitian Matrices
In mathematics, the Jacobi method for complex Hermitian matrices is a generalization of the Jacobi iteration method. The Jacobi iteration method is also explained in "Introduction to Linear Algebra" by . Derivation The complex unitary rotation matrices ''R''''pq'' can be used for Jacobi iteration of complex Hermitian matrices in order to find a numerical estimation of their eigenvectors and eigenvalues simultaneously. Similar to the Givens rotation matrices, ''R''''pq'' are defined as: : \begin (R_)_ & = \delta_ & \qquad m,n \ne p,q, \\0pt(R_)_ & = \frac e^, \\0pt(R_)_ & = \frac e^, \\0pt(R_)_ & = \frac e^, \\0pt(R_)_ & = \frac e^ \end Each rotation matrix, ''R''''pq'', will modify only the ''p''th and ''q''th rows or columns of a matrix ''M'' if it is applied from left or right, respectively: : \begin (R_ M)_ & = \begin M_ & m \ne p,q \\ pt\frac (M_ e^ - M_ e^) & m = p \\ pt\frac (M_ e^ + M_ e^) & m = q \end \\ pt(MR_^\dagger)_ & = \begin M_ & n \ne p,q \\ \frac (M_ e^ ...
[...More Info...]      
[...Related Items...]     OR:     [Wikipedia]   [Google]   [Baidu]  


Hermitian Matrices
In mathematics, a Hermitian matrix (or self-adjoint matrix) is a complex square matrix that is equal to its own conjugate transpose—that is, the element in the -th row and -th column is equal to the complex conjugate of the element in the -th row and -th column, for all indices and : or in matrix form: A \text \quad \iff \quad A = \overline . Hermitian matrices can be understood as the complex extension of real symmetric matrices. If the conjugate transpose of a matrix A is denoted by A^\mathsf, then the Hermitian property can be written concisely as Hermitian matrices are named after Charles Hermite, who demonstrated in 1855 that matrices of this form share a property with real symmetric matrices of always having real eigenvalues. Other, equivalent notations in common use are A^\mathsf = A^\dagger = A^\ast, although note that in quantum mechanics, A^\ast typically means the complex conjugate only, and not the conjugate transpose. Alternative characterizations Hermit ...
[...More Info...]      
[...Related Items...]     OR:     [Wikipedia]   [Google]   [Baidu]  


Jacobi Eigenvalue Algorithm
In numerical linear algebra, the Jacobi eigenvalue algorithm is an iterative method for the calculation of the eigenvalues and eigenvectors of a real symmetric matrix (a process known as diagonalization). It is named after Carl Gustav Jacob Jacobi, who first proposed the method in 1846, but only became widely used in the 1950s with the advent of computers. Description Let S be a symmetric matrix, and G=G(i,j,\theta) be a Givens rotation matrix. Then: :S'=G S G^\top \, is symmetric and similar to S. Furthermore, S^\prime has entries: :\begin S'_ &= c^2\, S_ - 2\, s c \,S_ + s^2\, S_ \\ S'_ &= s^2 \,S_ + 2 s c\, S_ + c^2 \, S_ \\ S'_ &= S'_ = (c^2 - s^2 ) \, S_ + s c \, (S_ - S_ ) \\ S'_ &= S'_ = c \, S_ - s \, S_ & k \ne i,j \\ S'_ &= S'_ = s \, S_ + c \, S_ & k \ne i,j \\ S'_ &= S_ &k,l \ne i,j \end where s=\sin(\theta) and c=\cos(\theta). Since G is orthogonal, S and S^\prime have the same Frobenius norm , , \cdot, , _F (the square-root sum of squares ...
[...More Info...]      
[...Related Items...]     OR:     [Wikipedia]   [Google]   [Baidu]  


Unitary Matrices
In linear algebra, a complex square matrix is unitary if its conjugate transpose is also its inverse, that is, if U^* U = UU^* = UU^ = I, where is the identity matrix. In physics, especially in quantum mechanics, the conjugate transpose is referred to as the Hermitian adjoint of a matrix and is denoted by a dagger (†), so the equation above is written U^\dagger U = UU^\dagger = I. The real analogue of a unitary matrix is an orthogonal matrix. Unitary matrices have significant importance in quantum mechanics because they preserve norms, and thus, probability amplitudes. Properties For any unitary matrix of finite size, the following hold: * Given two complex vectors and , multiplication by preserves their inner product; that is, . * is normal (U^* U = UU^*). * is diagonalizable; that is, is unitarily similar to a diagonal matrix, as a consequence of the spectral theorem. Thus, has a decomposition of the form U = VDV^*, where is unitary, and is diagonal and u ...
[...More Info...]      
[...Related Items...]     OR:     [Wikipedia]   [Google]   [Baidu]  


picture info

Rotation Matrix
In linear algebra, a rotation matrix is a transformation matrix that is used to perform a rotation in Euclidean space. For example, using the convention below, the matrix :R = \begin \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end rotates points in the plane counterclockwise through an angle with respect to the positive axis about the origin of a two-dimensional Cartesian coordinate system. To perform the rotation on a plane point with standard coordinates , it should be written as a column vector, and multiplied by the matrix : : R\mathbf = \begin \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end \begin x \\ y \end = \begin x\cos\theta-y\sin\theta \\ x\sin\theta+y\cos\theta \end. If and are the endpoint coordinates of a vector, where is cosine and is sine, then the above equations become the trigonometric summation angle formulae. Indeed, a rotation matrix can be seen as the trigonometric summation angle formulae in matrix form. One w ...
[...More Info...]      
[...Related Items...]     OR:     [Wikipedia]   [Google]   [Baidu]  




Givens Rotation
In numerical linear algebra, a Givens rotation is a rotation in the plane spanned by two coordinates axes. Givens rotations are named after Wallace Givens, who introduced them to numerical analysts in the 1950s while he was working at Argonne National Laboratory. Matrix representation A Givens rotation is represented by a matrix of the form :G(i, j, \theta) = \begin 1 & \cdots & 0 & \cdots & 0 & \cdots & 0 \\ \vdots & \ddots & \vdots & & \vdots & & \vdots \\ 0 & \cdots & c & \cdots & -s & \cdots & 0 \\ \vdots & & \vdots & \ddots & \vdots & & \vdots \\ 0 & \cdots & s & \cdots & c & \cdots & 0 \\ \vdots & & \vdots & & \vdots & \ddots & \vdots \\ 0 & \cdots & 0 & \cdots & 0 & \cdots & 1 \end, where and appear at the int ...
[...More Info...]      
[...Related Items...]     OR:     [Wikipedia]   [Google]   [Baidu]  


Hermitian Matrix
In mathematics, a Hermitian matrix (or self-adjoint matrix) is a complex square matrix that is equal to its own conjugate transpose—that is, the element in the -th row and -th column is equal to the complex conjugate of the element in the -th row and -th column, for all indices and : or in matrix form: A \text \quad \iff \quad A = \overline . Hermitian matrices can be understood as the complex extension of real symmetric matrices. If the conjugate transpose of a matrix A is denoted by A^\mathsf, then the Hermitian property can be written concisely as Hermitian matrices are named after Charles Hermite, who demonstrated in 1855 that matrices of this form share a property with real symmetric matrices of always having real eigenvalues. Other, equivalent notations in common use are A^\mathsf = A^\dagger = A^\ast, although note that in quantum mechanics, A^\ast typically means the complex conjugate only, and not the conjugate transpose. Alternative characterizations Hermit ...
[...More Info...]      
[...Related Items...]     OR:     [Wikipedia]   [Google]   [Baidu]