Compute the two sides of the normal equations for linear least squares and divide by :
The first equation gives . Substitute this into the second and solve:
These simple linear regression coefficients are unique provided the observed distances are not all equal, since
The condition alone does not guarantee this. If all , the design matrix has matrix rank one and the valid alternative is the family
In this rank-deficient ordinary least squares case, the fitted predictions are uniquely but the two parameters cannot be separately determined.
A two-index Cartesian tensor has components transforming as under an orthogonal change of orthonormal basis. The tensor order here counts indices, rather than the matrix rank. Its scalar, symmetric-traceless and axial tensor decomposition follows by setting
Then , with antisymmetric and symmetric and traceless. There are three independent antisymmetric components and five independent symmetric traceless components: three off-diagonal entries and two independent diagonal ones in the latter case.
For the symmetric tensor relative to a nonzero vector , put , , and . The symmetric tensor decomposition relative to an axis gives
In particular . Define the orthogonal projection onto the plane perpendicular to , and set
This tensor is symmetric and kills . Its trace is . To verify the full decomposition, set and split into its longitudinal component , two mixed components, and . The mixed terms are exactly , and gives the longitudinal coefficient. Thus all requested constraints and the reconstruction hold.
All singular values are positive, so and are invertible matrices and has matrix rank . The image of a linear map is therefore
with those columns furnishing an orthonormal basis. Also implies and hence . Thus
Finally exactly when . The remaining columns of the full orthogonal matrix span this orthogonal complement:
If , this last null space also equals and has the empty basis.
No. Distinct singular values do not force uniqueness in the spectral norm. For an explicit example take
The singular values of are . Every has matrix rank one, and its error in the spectral norm is
Thus there are infinitely many best rank-one approximations even with a strict gap between the singular values.
The same mechanism works for every under the question's positive singular values. Starting with , define
When the singular values are distinct, , so every has matrix rank , while
The spectral norm measures only the largest remaining singular value, leaving room to change retained components without increasing that maximum. This is the nonuniqueness of spectral-norm low-rank approximation. The exceptional case is unique for the elementary reason given above.
Take the identity matrix
Both and have matrix rank one, and
They are distinct and both attain the minimum spectral norm error from the preceding part. Hence a best low-rank approximation of rank one need not be unique. More generally, every orthogonal projection with gives another such minimizer.
Retain the first terms in the singular value decomposition:
The Eckart-Young-Mirsky theorem says that this truncated singular value decomposition is a best low-rank approximation in the spectral norm. The remaining singular values of the difference are , together with zeros, so
For the displayed matrix has matrix rank exactly , since every retained singular value is positive. For , the only matrix of matrix rank zero is the zero matrix, and the same formula gives error .
A full real singular value decomposition is
where has only diagonal entries, with . These are the singular values. The number of positive singular values is the matrix rank. The columns of are right singular vectors and those of are left singular vectors, giving orthonormal input and output coordinate systems.
The nullspace is the vector subspace . A full real singular value decomposition is
where are orthogonal matrices and for , with for . We may order . The columns of form an orthonormal basis of the input space , and the columns of form one of the output space . The positive singular values number exactly the matrix rank.
Let be the matrix rank, and partition the right singular vectors as , with consisting of the columns associated with zero singular values. From , these columns form an orthonormal basis of the nullspace. Therefore its orthogonal projector is
It acts on the input space . Symmetry and idempotence follow from . At full column rank it is zero, and at rank zero it is the identity.
A real singular value decomposition has the form , where and are orthogonal matrices, and is a rectangular diagonal matrix with entries . These entries are the singular values. Their squares are the eigenvalues of , omitting the additional zero eigenvalues when . The positive singular values number exactly the matrix rank.
By the Eckart-Young-Mirsky theorem, truncate the singular value decomposition:
The same minimizes the spectral norm and Frobenius norm errors over matrices of matrix rank at most . This conventional meaning of rank- low-rank approximation includes the case , when .
Put . Since has columns, the rank bound for a matrix product gives . The Eckart-Young-Mirsky theorem and part (a) then give
One can also see the lower bound directly. The restriction of to the -dimensional linear span of has a nonzero null space, because its matrix rank is at most . For a unit vector in that null space, by the singular value decomposition.
Minimizing gives the normal equations . Since has full column rank, is positive definite and invertible, and the objective has a unique minimum. Set ; then and . The model gives and , so and . For the alternative estimator the corresponding general formulas are and .
Inverting the detector rotation at this time gives
For a representation with one system factor and one environment factor per displayed detector state, each parenthesized system–environment vector would have to be a product state. But the first has coefficient matrix in the indicated orthonormal bases; the second has matrix . Their matrix ranks are two whenever . Equivalently, their two nonzero Schmidt coefficients are , which no independent basis changes can remove. Thus the proposed factorization is impossible when .
The exception matters: at , the detector basis is the original basis up to phases or interchange, and a factorization does exist. In particular, the literal condition is insufficient if it includes . The intended obstruction concerns genuinely mixed detector bases.
The environment has selected robust detector records: its interaction commutes with the original projectors, and orthogonal environmental records suppress their interference in the reduced density matrix. A rotated detector branch instead leaves the system and environment entangled. This is environment-induced pointer-basis selection, the issue analyzed by W. H. Zurek in his apparatus-environment model.[ref]
This addresses the basis ambiguity but does not by itself select one realized outcome from the global superposition. Moreover, a two-state environment has recurrences: returns to unit modulus, so the toy model does not establish irreversible macroscopic collapse. A large environment can make records effectively persistent over experimental times, but that requires further dynamical assumptions. The distinction is part of the measurement problem, not an algebraic error cured by a partial trace.
The elementary row operations preserve the determinant, and the diagonal product of the upper triangular matrix found above is
Consequently
At , all four rows are equal. At , successive rows alternate between one row and its negative. In either case the matrix rank is one, confirming the vanishing determinant.
For an matrix over the real numbers, write its rows as and its columns as . Then
Thus the row space and column space consist of all linear combinations of the corresponding rows and columns; the null space is the kernel of a linear map represented by . The matrix rank counts independent row or column directions, with row rank equal to column rank.
For a real matrix , the row space is the linear span of its rows in . The column space is the linear span of its columns in , equivalently . The null space is
The matrix rank is the dimension of the column space, which equals the dimension of the row space. In particular,
The last equality is the rank-nullity theorem.
The reduced row echelon form gives , , and . With free parameter ,
The displayed nonzero vector therefore forms a basis for the null space:
This agrees with the rank-nullity theorem: there are four columns and matrix rank three.
Let and . The remaining columns are and , so these two vectors span the column space. They are linearly independent: the determinant of their first two coordinates is . Therefore is a basis of the range, and the matrix rank is two.
Equivalently the range is the plane : both displayed basis vectors satisfy this equation, and their two-dimensional span fills that plane.