|
|
| Line 1: |
Line 1: |
| In mathematics, the '''Smith normal form''' is a [[Canonical form|normal form]] that can be defined for any matrix (not necessarily square) with entries in a [[principal ideal domain]] (PID). The Smith normal form of a matrix is [[Diagonal matrix|diagonal]], and can be obtained from the original matrix by multiplying on the left and right by [[invertible]] square matrices. In particular, the integers are a PID, so one can always calculate the Smith normal form of an integer matrix. The Smith normal form is very useful for working with finitely generated modules over a PID, and in particular for deducing the structure of a quotient of a [[free module]].
| | Bryan is really a celebrity during the generating along with the profession progress initially [http://lukebryantickets.sgs-suparco.org rihanna concert] next to his third theatre record, And , is [http://lukebryantickets.citizenswebcasting.com miley cyrus concert tickets] the proof. He broken on the scene in 2011 along with his amazing blend of down-home availability, motion picture celebrity wonderful appears and words, is scheduled t within a key way. The brand new recor in the nation graph and #2 about the pop maps, generating it the 2nd maximum debut during that time of 2005 for a nation musician. <br><br> |
|
| |
|
| ==Definition==
| | The boy of the , knows perseverance and perseverance are important elements in terms of a successful occupation- . His initially record, Stay Me, produced the very best hits “All My Friends “Country and Say” Gentleman,” although his energy, luke bryan tickets indiana [[http://lukebryantickets.asiapak.net previous]] Doin’ Thing, discovered the singer-about three directly No. 1 single people: [http://www.ffpjp24.org tickets luke bryan] Different Contacting Is really a Wonderful Point.”<br><br>In the fall of 2002, Concert tours: Luke And that have a remarkable set of , such as Downtown. “It’s almost like you are obtaining a authorization to look to another level, affirms those designers which were a part of the Concert tourabove in to a bigger degree of musicians.” It wrapped among the most successful tours in their 10-year history.<br><br>My blog post ... [http://okkyunglee.com prince tickets] |
| | |
| Let ''A'' be a nonzero ''m''×''n'' matrix over a [[principal ideal domain]] ''R''. There exist invertible <math>m \times m</math> and <math>n \times n</math>-matrices ''S, T'' so that the product ''S A T'' is
| |
| | |
| :<math>
| |
| \begin{pmatrix}
| |
| \alpha_1 & 0 & 0 & & \cdots & & 0 \\
| |
| 0 & \alpha_2 & 0 & & \cdots & & 0 \\
| |
| 0 & 0 & \ddots & & & & 0\\
| |
| \vdots & & & \alpha_r & & & \vdots \\
| |
| & & & & 0 & & \\
| |
| & & & & & \ddots & \\
| |
| 0 & & & \cdots & & & 0
| |
| \end{pmatrix}.
| |
| </math>
| |
| | |
| and the diagonal elements <math>\alpha_i</math> satisfy <math>\alpha_i \mid \alpha_{i+1}\;\forall\;1 \le i < r</math>. This is the Smith normal form of the matrix ''A''. The elements <math>\alpha_i</math> are unique [[up to]] multiplication by a [[Unit (ring theory)|unit]] and are called the ''elementary divisors'', ''invariants'', or ''invariant factors''. They can be computed (up to multiplication by a unit) as
| |
| : <math>\alpha_i = \frac{d_i(A)}{d_{i-1}(A)},</math>
| |
| where <math>d_i(A)</math> (called ''i''-th ''determinant divisor'') equals the [[greatest common divisor]] of all <math>i\times i</math> [[Minor (linear algebra)|minors]] of the matrix ''A''.
| |
| | |
| == Algorithm ==
| |
| Our first goal will be to find invertible square matrices ''S'' and ''T'' such that the product ''S A T'' is diagonal. This is the hardest part of the algorithm and once we have achieved diagonality it becomes relatively easy to put the matrix in Smith normal form. Phrased more abstractly, the goal is to show that, thinking of ''A'' as a map from <math>R^n</math> (the free ''R''-[[Module (mathematics)|module]] of rank ''n'') to <math>R^m</math> (the free ''R''-[[Module (mathematics)|module]] of rank ''m''), there are isomorphisms <math>S:R^m \to R^m</math> and <math>T:R^n \to R^n</math> such that <math>S \cdot A \cdot T</math> has the simple form of a [[diagonal matrix]]. The matrices ''S'' and ''T'' can be found by starting out with identity matrices of the appropriate size, and modifying ''S'' each time a row operation is performed on ''A'' in the algorithm by the same row operation, and similarly modifying ''T'' for each column operation performed. Since row operations are left-multiplications and column operations are right-multiplications, this preserves the invariant <math>A'=S'\cdot A\cdot T'</math> where <math>A',S',T'</math> denote current values and ''A'' denotes the original matrix; eventually the matrices in this invariant become diagonal. Only invertible row and column operations are performed, which ensures that ''S'' and ''T'' remain invertible matrices.
| |
| | |
| For ''a'' in ''R'' \ {0}, write δ(''a'') for the number of prime factors of ''a'' (these exist and are unique since any PID is also a [[unique factorization domain]]). In particular, ''R'' is also a [[Bézout domain]], so it is a [[gcd domain]] and the gcd of any two elements satisfies a [[Bézout's identity]].
| |
| | |
| To put a matrix into Smith normal form, one can repeatedly apply the following, where ''t'' loops from 1 to ''m''.
| |
| | |
| ===Step I: Choosing a pivot===
| |
| Choose ''j''<sub>''t''</sub> to be the smallest column index of ''A'' with a non-zero entry, starting the search at column index ''j''<sub>''t''-1</sub>+1 if ''t'' > 1.
| |
| | |
| We wish to have <math>a_{t,j_t}\neq0</math>; if this is the case this step is complete, otherwise there is by assumption some ''k'' with <math>a_{k,j_t} \neq 0</math>, and we can exchange rows <math>t</math> and ''k'', thereby obtaining <math>a_{t,j_t}\neq0</math>.
| |
| | |
| Our chosen pivot is now at position (''t'', ''j''<sub>''t''</sub>).
| |
| | |
| ===Step II: Improving the pivot===
| |
| If there is an entry at position (''k'',''j''<sub>''t''</sub>) such that <math>a_{t,j_t} \nmid a_{k,j_t}</math>, then, letting <math>\beta =\gcd\left(a_{t,j_t}, a_{k,j_t}\right)</math>, we know by the Bézout property that there exist σ, τ in ''R'' such that
| |
| | |
| :<math>
| |
| a_{t,j_t} \cdot \sigma + a_{k,j_t} \cdot \tau=\beta.
| |
| </math>
| |
| | |
| By left-multiplication with an appropriate invertible matrix ''L'', it can be achieved that row ''t'' of the matrix product is the sum of σ times the original row ''t'' and τ times the original row ''k'', that row ''k'' of the product is another linear combination of those original rows, and that all other rows are unchanged. Explicitly, if σ and τ satisfy the above equation, then for <math>\alpha=a_{t,j_t}/\beta</math> and <math>\gamma=a_{k,j_t}/\beta</math> (which divisions are possible by the definition of β) one has
| |
| | |
| :<math>
| |
| \sigma\cdot \alpha + \tau \cdot \gamma=1,
| |
| </math>
| |
| | |
| so that the matrix
| |
| | |
| :<math> L_0= | |
| \begin{pmatrix}
| |
| \sigma & \tau \\
| |
| -\gamma & \alpha \\
| |
| \end{pmatrix}
| |
| </math>
| |
| | |
| is invertible, with inverse
| |
| | |
| :<math> | |
| \begin{pmatrix}
| |
| \alpha & -\tau \\
| |
| \gamma & \sigma \\
| |
| \end{pmatrix}
| |
| .</math>
| |
| | |
| Now ''L'' can be obtained by fitting <math>L_0</math> into rows and columns ''t'' and ''k'' of the identity matrix. By construction the matrix obtained after left-multiplying by ''L'' has entry β at position (''t'',''j''<sub>''t''</sub>) (and due to our choice of α and γ it also has an entry 0 at position (''k'',''j''<sub>''t''</sub>), which is useful though not essential for the algorithm). This new entry β divides the entry <math>a_{t,j_t}</math> that was there before, and so in particular <math>\delta(\beta) < \delta(a_{t,j_t})</math>; therefore repeating these steps must eventually terminate. One ends up with a matrix having an entry at position (''t'',''j''<sub>''t''</sub>) that divides all entries in column ''j''<sub>''t''</sub>.
| |
| | |
| ===Step III: Eliminating entries===
| |
| Finally, adding appropriate multiples of row ''t'', it can be achieved that all entries in column ''j''<sub>''t''</sub> except for that at position (''t'',''j''<sub>''t''</sub>) are zero. This can be achieved by left-multiplication with an appropriate matrix. However, to make the matrix fully diagonal we need to eliminate nonzero entries on the row of position (''t'',''j''<sub>''t''</sub>) as well. This can be achieved by repeating the steps in Step II for columns instead of rows, and using multiplication on the right. In general this will result in the zero entries from the prior application of Step III becoming nonzero again.
| |
| | |
| However, notice that the [[ideal (ring theory)|ideals]] generated by the elements at position (''t'',''j''<sub>''t''</sub>) form an [[ascending chain condition|ascending chain]], because entries from a later step always divide entries from a previous step. Therefore, since ''R'' is a [[Noetherian ring]] (it is a [[principal ideal domain|PID]]), the ideals eventually become stationary and do not change. This means that at some stage after Step II has been applied, the entry at (''t'',''j''<sub>''t''</sub>) will divide all nonzero row or column entries before applying any more steps in Step II. Then we can eliminate entries in the row or column with nonzero entries while preserving the zeros in the already-zero row or column. At this point, only the block of ''A'' to the lower right of (''t'',''j''<sub>''t''</sub>) needs to be diagonalized, and conceptually the algorithm can be applied recursively, treating this block as a separate matrix. In other words, we can increment ''t'' by one and go back to Step I.
| |
| | |
| ===Final step===
| |
| Applying the steps described above to the remaining non-zero columns of the resulting matrix (if any), we get an <math>m \times n</math>-matrix with column indices <math>j_1 < \ldots < j_r</math> where <math>r \le \min(m,n)</math>. The matrix entries <math>(l,j_l)</math> are non-zero, and every other entry is zero.
| |
| | |
| Now we can move the null columns of this matrix to the right, so that the nonzero entries are on positions <math>(i,i)</math> for <math>1 \le i\le r</math>. For short, set <math>\alpha_i</math> for the element at position <math>(i,i)</math>.
| |
| | |
| The condition of divisibility of diagonal entries might not be satisfied. For any index <math>i<r</math> for which <math>\alpha_i\nmid\alpha_{i+1}</math>, one can repair this shortcoming by operations on rows and columns <math>i</math> and <math>i+1</math> only: first add column <math>i+1</math> to column <math>i</math> to get an entry <math>\alpha_{i+1}</math> in column ''i'' without disturbing the entry <math>\alpha_i</math> at position <math>(i,i)</math>, and then apply a row operation to make the entry at position <math>(i,i)</math> equal to <math>\beta=\gcd(\alpha_i,\alpha_{i+1})</math> as in Step II; finally proceed as in Step III to make the matrix diagonal again. Since the new entry at position <math>(i+1,i+1)</math> is a linear combination of the original <math>\alpha_i,\alpha_{i+1}</math>, it is divisible by β.
| |
| | |
| The value <math>\delta(\alpha_1)+\cdots+\delta(\alpha_r)</math> does not change by the above operation (it is δ of the determinant of the upper <math>r\times r</math> submatrix), whence that operation does diminish (by moving prime factors to the right) the value of
| |
| :<math>\sum_{j=1}^r(r-j)\delta(\alpha_j).</math>
| |
| So after finitely many applications of this operation no further application is possible, which means that we have obtained <math>\alpha_1\mid\alpha_2\mid\cdots\mid\alpha_r</math> as desired.
| |
| | |
| Since all row and column manipulations involved in the process are invertible, this shows that there exist invertible <math>m \times m</math> and <math>n \times n</math>-matrices ''S, T'' so that the product ''S A T'' satisfies the definition of a Smith normal form. In particular, this shows that the Smith normal form exists, which was assumed without proof in the definition.
| |
| | |
| == Applications ==
| |
| | |
| The Smith normal form is useful for computing the [[homology (mathematics)|homology]] of a [[chain complex]] when the chain modules of the chain complex are [[Finitely generated module|finitely generated]]. For instance, in [[topology]], it can be used to compute the homology of a [[simplicial complex]] or [[CW complex]] over the integers, because the boundary maps in such a complex are just integer matrices. It can also be used to prove the well known [[structure theorem for finitely generated modules over a principal ideal domain]].
| |
| | |
| == Example ==
| |
| As an example, we will find the Smith normal form of the following matrix over the integers.
| |
| | |
| :<math>
| |
| \begin{pmatrix}
| |
| 2 & 4 & 4 \\
| |
| -6 & 6 & 12 \\
| |
| 10 & -4 & -16 | |
| \end{pmatrix}
| |
| </math>
| |
| | |
| The following matrices are the intermediate steps as the algorithm is applied to the above matrix.
| |
| | |
| :<math>
| |
| \to
| |
| \begin{pmatrix}
| |
| 2 & 0 & 0 \\
| |
| -6 & 18 & 24 \\
| |
| 10 & -24& -36
| |
| \end{pmatrix}
| |
| \to
| |
| \begin{pmatrix}
| |
| 2 & 0 & 0 \\
| |
| 0 & 18 & 24 \\
| |
| 0 & -24& -36
| |
| \end{pmatrix}
| |
| </math> | |
| | |
| :<math>
| |
| \to
| |
| \begin{pmatrix}
| |
| 2 & 0 & 0 \\
| |
| 0 & 18 & 24 \\
| |
| 0 & -6 & -12
| |
| \end{pmatrix}
| |
| \to
| |
| \begin{pmatrix}
| |
| 2 & 0 & 0 \\
| |
| 0 & 6 & 12 \\
| |
| 0 & 18 & 24
| |
| \end{pmatrix}
| |
| </math>
| |
| | |
| :<math>
| |
| \to
| |
| \begin{pmatrix}
| |
| 2 & 0 & 0 \\
| |
| 0 & 6 & 12 \\
| |
| 0 & 0 & -12
| |
| \end{pmatrix}
| |
| \to
| |
| \begin{pmatrix}
| |
| 2 & 0 & 0 \\
| |
| 0 & 6 & 0 \\
| |
| 0 & 0 & 12
| |
| \end{pmatrix}
| |
| </math>
| |
| | |
| So the Smith normal form is
| |
| | |
| :<math>
| |
| \begin{pmatrix}
| |
| 2 & 0 & 0 \\
| |
| 0 & 6 & 0 \\
| |
| 0 & 0 & 12
| |
| \end{pmatrix}
| |
| </math>
| |
| | |
| and the elementary divisors are 2, 6 and 12.
| |
| | |
| == Similarity ==
| |
| The Smith normal form can be used to determine whether or not matrices with entries over a common field are [[similar (linear algebra)|similar]]. Specifically two matrices ''A'' and ''B'' are similar if and only if the [[characteristic matrix|characteristic matrices]] <math>xI-A</math> and <math>xI-B</math> have the same Smith normal form.
| |
| | |
| For example, with
| |
| :<math>
| |
| \begin{align}
| |
| A & {} =\begin{bmatrix}
| |
| 1 & 2 \\
| |
| 0 & 1
| |
| \end{bmatrix}, & & \mbox{SNF}(xI-A) =\begin{bmatrix}
| |
| 1 & 0 \\
| |
| 0 & (x-1)^2
| |
| \end{bmatrix} \\
| |
| B & {} =\begin{bmatrix}
| |
| 3 & -4 \\
| |
| 1 & -1
| |
| \end{bmatrix}, & & \mbox{SNF}(xI-B) =\begin{bmatrix}
| |
| 1 & 0 \\
| |
| 0 & (x-1)^2
| |
| \end{bmatrix} \\
| |
| C & {} =\begin{bmatrix}
| |
| 1 & 0 \\
| |
| 1 & 2
| |
| \end{bmatrix}, & & \mbox{SNF}(xI-C) =\begin{bmatrix}
| |
| 1 & 0 \\
| |
| 0 & (x-1)(x-2)
| |
| \end{bmatrix}.
| |
| \end{align}
| |
| </math>
| |
| | |
| ''A'' and ''B'' are similar because the Smith normal form of their characteristic matrices match, but are not similar to ''C'' because the Smith normal form of the characteristic matrices do not match.
| |
| | |
| == See also ==
| |
| * [[Canonical form]]
| |
| * [[Elementary divisors]]
| |
| * [[Frobenius normal form]]
| |
| * [[Hermite normal form]]
| |
| * [[Invariant factor]]
| |
| * [[Henry John Stephen Smith]] (1826–1883), [[eponym]] of the Smith normal form
| |
| * [[Structure theorem for finitely generated modules over a principal ideal domain]]
| |
| | |
| ==References==
| |
| * {{cite journal |last=Smith |first=Henry J. Stephen |authorlink=Henry John Stephen Smith |year=1861 |title=On systems of linear indeterminate equations and congruences |journal=[[Philosophical Transactions of the Royal Society of London|Phil. Trans. R. Soc. Lond.]] |volume=151 |issue=1 |pages=293–326 |url=http://www.jstor.org/stable/108738 |doi=10.1098/rstl.1861.0016 }} Reprinted (pp. [http://archive.org/stream/collectedmathema01smituoft#page/366/mode/2up 367–409]) in [http://archive.org/details/collectedmathema01smituoft ''The Collected Mathematical Papers of Henry John Stephen Smith'', Vol. I], edited by [[James Whitbread Lee Glaisher|J. W. L. Glaisher]]. Oxford: Clarendon Press (1894), ''xcv''+603 pp.
| |
| * {{PlanetMath |urlname=GausssAlgorithmForPrincipalIdealDomains |title=Smith normal form}}
| |
| * {{PlanetMath |urlname=ExampleOfSmithNormalForm |title=Example of Smith normal form}}
| |
| * K. R. Matthews, [http://www.numbertheory.org/courses/MP274/smith.pdf Smith normal form]. MP274: Linear Algebra, Lecture Notes, University of Queensland, 1991.
| |
| | |
| [[Category:Matrix theory]]
| |
| [[Category:Matrix normal forms]]
| |
Bryan is really a celebrity during the generating along with the profession progress initially rihanna concert next to his third theatre record, And , is miley cyrus concert tickets the proof. He broken on the scene in 2011 along with his amazing blend of down-home availability, motion picture celebrity wonderful appears and words, is scheduled t within a key way. The brand new recor in the nation graph and #2 about the pop maps, generating it the 2nd maximum debut during that time of 2005 for a nation musician.
The boy of the , knows perseverance and perseverance are important elements in terms of a successful occupation- . His initially record, Stay Me, produced the very best hits “All My Friends “Country and Say” Gentleman,” although his energy, luke bryan tickets indiana [previous] Doin’ Thing, discovered the singer-about three directly No. 1 single people: tickets luke bryan Different Contacting Is really a Wonderful Point.”
In the fall of 2002, Concert tours: Luke And that have a remarkable set of , such as Downtown. “It’s almost like you are obtaining a authorization to look to another level, affirms those designers which were a part of the Concert tourabove in to a bigger degree of musicians.” It wrapped among the most successful tours in their 10-year history.
My blog post ... prince tickets