Monitor (synchronization)

In mathematics, particularly in linear algebra and functional analysis, the polar decomposition of a matrix or linear operator is a factorization analogous to the polar form of a nonzero complex number z as $z=re^{i\theta }\,$ where r is the absolute value of z (a positive real number), and $e^{i\theta }$ is called the complex sign of z.

Matrix polar decomposition

The polar decomposition of a square complex matrix A is a matrix decomposition of the form

$A=UP\,$

where U is a unitary matrix and P is a positive-semidefinite Hermitian matrix. Intuitively, the polar decomposition separates A into a component that stretches the space along a set of orthogonal axes, represented by P, and a rotation represented by U. The decomposition of the complex conjugate of $A$ is given by ${\overline {A}}={\overline {U}}$ ${\overline {P}}$.

This decomposition always exists; and so long as A is invertible, it is unique, with P positive-definite. Note that

$\det A=\det U\,\det P=re^{i\theta }$

gives the corresponding polar decomposition of the determinant of A, since $\det P=r=|\det A|$ and $\det U=e^{i\theta }$.

The matrix P is always unique, even if A is singular, and given by

$P={\sqrt {A^{*}A}}$

where A* denotes the conjugate transpose of A. This expression is meaningful since a positive-semidefinite Hermitian matrix has a unique positive-semidefinite square root. If A is invertible, then the matrix U is given by

$U=AP^{-1}.\,$

In terms of the singular value decomposition of A, A = W Σ V^*, one has

$P=V\Sigma V^{*}\,$

$U=WV^{*}\,$

confirming that P is positive-definite and U is unitary.

One can also decompose A in the form

$A=P'U\,$

Here U is the same as before and P′ is given by

$P'=UPU^{-1}={\sqrt {AA^{*}}}=W\Sigma W^{*}.$

This is known as the left polar decomposition, whereas the previous decomposition is known as the right polar decomposition. Left polar decomposition is also known as reverse polar decomposition.

The matrix A is normal if and only if P′ = P. Then UΣ = ΣU, and it is possible to diagonalise U with a unitary similarity matrix S that commutes with Σ, giving S U S* = Φ⁻¹, where Φ is a diagonal unitary matrix of phases e^iφ. Putting Q = V S^*, one can then re-write the polar decomposition as

$A=(Q\Phi Q^{*})(Q\Sigma Q^{*}),\,$

so A then thus also has a spectral decomposition

$A=Q\Lambda Q^{*}\,$

with complex eigenvalues such that ΛΛ^* = Σ² and a unitary matrix of complex eigenvectors Q.

Bounded operators on Hilbert space

The polar decomposition of any bounded linear operator A between complex Hilbert spaces is a canonical factorization as the product of a partial isometry and a non-negative operator.

The polar decomposition for matrices generalizes as follows: if A is a bounded linear operator then there is a unique factorization of A as a product A = UP where U is a partial isometry, P is a non-negative self-adjoint operator and the initial space of U is the closure of the range of P.

The operator U must be weakened to a partial isometry, rather than unitary, because of the following issues. If A is the one-sided shift on l²(N), then |A| = {A*A}^½ = I. So if A = U |A|, U must be A, which is not unitary.

The existence of a polar decomposition is a consequence of Douglas' lemma:

Lemma If A, B are bounded operators on a Hilbert space H, and A*A ≤ B*B, then there exists a contraction C such that A = CB. Furthermore, C is unique if Ker(B*) ⊂ Ker(C).

The operator C can be defined by C(Bh) := Ah for all h in H, extended by continuity to the closure of Ran(B), and by zero on the orthogonal complement to all of H. The lemma then follows since A*A ≤ B*B implies Ker(A) ⊂ Ker(B).

In particular. If A*A = B*B, then C is a partial isometry, which is unique if Ker(B*) ⊂ Ker(C). In general, for any bounded operator A,

$A^{*}A=(A^{*}A)^{\frac {1}{2}}(A^{*}A)^{\frac {1}{2}},$

where (A*A)^½ is the unique positive square root of A*A given by the usual functional calculus. So by the lemma, we have

$A=U(A^{*}A)^{\frac {1}{2}}$

for some partial isometry U, which is unique if Ker(A*) ⊂ Ker(U). Take P to be (A*A)^½ and one obtains the polar decomposition A = UP. Notice that an analogous argument can be used to show A = P'U' , where P' is positive and U' a partial isometry.

When H is finite dimensional, U can be extended to a unitary operator; this is not true in general (see example above). Alternatively, the polar decomposition can be shown using the operator version of singular value decomposition.

By property of the continuous functional calculus, |A| is in the C*-algebra generated by A. A similar but weaker statement holds for the partial isometry: U is in the von Neumann algebra generated by A. If A is invertible, the polar part U will be in the C*-algebra as well.

Unbounded operators

If A is a closed, densely defined unbounded operator between complex Hilbert spaces then it still has a (unique) polar decomposition

$A=U|A|\,$

where |A| is a (possibly unbounded) non-negative self adjoint operator with the same domain as A, and U is a partial isometry vanishing on the orthogonal complement of the range Ran(|A|).

The proof uses the same lemma as above, which goes through for unbounded operators in general. If Dom(A*A) = Dom(B*B) and A*Ah = B*Bh for all h ∈ Dom(A*A), then there exists a partial isometry U such that A = UB. U is unique if Ran(B)^⊥⊂ Ker(U). The operator A being closed and densely defined ensures that the operator A*A is self-adjoint (with dense domain) and therefore allows one to define (A*A)^½. Applying the lemma gives polar decomposition.

If an unbounded operator A is affiliated to a von Neumann algebra M, and A = UP is its polar decomposition, then U is in M and so is the spectral projection of P, 1_B(P), for any Borel set B in [0, ∞).

Quaternion polar decomposition

The polar decomposition of quaternions H depends on the sphere $\lbrace xi+yj+zk\in H:x^{2}+y^{2}+z^{2}=1\rbrace$ of square roots of minus one. Given any r on this sphere, and an angle –π < a ≤ π, the versor $e^{ar}=\cos(a)+r\ \sin(a)$ is on the 3-sphere of H. For a = 0 and a = π, the versor is 1 or −1 regardless of which r is selected. The norm t of a quaternion q is the Euclidean distance from the origin to q. When a quaternion is not just a real number, then there is a unique polar decomposition $q=te^{ar}\!.$

Alternative planar decompositions

In the Cartesian plane, alternative planar ring decompositions arise as follows:

• If x ≠ 0, z = x ( 1 + (y/x) ε) is a polar decomposition of a dual number z = x + y ε, where ε² = 0, i.e. ε is nilpotent. In this polar decomposition, the unit circle has been replaced by the line x = 1, the polar angle by the slope y/x, and the radius x is negative in the left half-plane.

• If x² ≠ y², then the unit hyperbola x² − y^{2 = 1 and its conjugate x2 − y2 = −1 can be used to form a polar decomposition based on the branch of the unit hyperbola through (1,0). This branch is parametrized by the hyperbolic angle a and is written}

$\cosh(a)+j\ \sinh(a)=\exp(aj)=e^{aj}$

where j ² = +1 and the arithmetic of split-complex numbers is used. The branch through (−1,0) is traced by −e^{a j}. Since the operation of multiplying by j reflects a point across the line y = x, the second hyperbola has branches traced by je^{a j} or −je^{a j}. Therefore a point in one of the quadrants has a polar decomposition in one of the forms:

$re^{aj},-re^{aj},rje^{aj},-rje^{aj},\quad r>0$

The set { 1, −1, j, −j } has products that make it isomorphic to the Klein four-group. Evidently polar decomposition in this case involves an element from that group.

Numerical determination of the matrix polar decomposition

To compute an approximation of the polar decomposition A=UP, usually the unitary factor U is approximated. The iteration is based on Heron's method for the square root of 1 and computes, starting from $U_{0}=A$, the sequence

$U_{k+1}={\frac {1}{2}}\left(U_{k}+(U_{k}^{*})^{-1}\right)$, k=0,1,2,...

The combination of inversion and Hermite conjugation is chosen so that in the singular value decomposition, the unitary factors remain the same and the iteration reduces to Heron's method on the singular values.

This basic iteration may be refined to speed up the process:

• Every step or in regular intervals, the range of the singular values of $U_{k}$ is estimated and then the matrix is rescaled to $\gamma _{k}U_{k}$ to center the singular values around 1. The scaling factor $\gamma _{k}$ is computed using matrix norms of the matrix and its inverse. Examples of such scale estimates are:

$\gamma _{k}={\sqrt[{4\;}]{\frac {\|U_{k}^{-1}\|_{1}\,\|U_{k}^{-1}\|_{\infty }}{\|U_{k}\|_{1}\,\|U_{k}\|_{\infty }}}}$

using the row-sum and column-sum matrix norms or

$\gamma _{k}={\sqrt {\frac {\|U_{k}^{-1}\|_{F}}{\|U_{k}\|_{F}}}}$

using the Frobenius norm. Including the scale factor, the iteration is now

$U_{k+1}={\frac {1}{2}}\left(\gamma _{k}U_{k}+{\frac {1}{\gamma _{k}}}(U_{k}^{*})^{-1}\right)$, k=0,1,2,...

• The QR decomposition can be used in a preparation step to reduce a singular matrix A to a smaller regular matrix, and inside every step to speed up the computation of the inverse.

• Heron' method for computing roots of $x^{2}-1=0$ can be replaced by higher order methods, for instance based on Halley's method of third order, resulting in

$U_{k+1}=U_{k}\left(I+3\,U_{k}^{*}U_{k}\right)^{-1}\left(3\,I+U_{k}^{*}U_{k}\right)$, k=0,1,2,...

This iteration can again be combined with rescaling. This particular formula has the benefit that it also applicable to singular or rectangular matrices A.

See also

• Cartan decomposition

Literature

• Conway, J.B.: A Course in Functional Analysis. Graduate Texts in Mathematics. New York: Springer 1990
• Douglas, R.G.: On Majorization, Factorization, and Range Inclusion of Operators on Hilbert Space. Proc. Amer. Math. Soc. 17, 413-415 (1966)
• Many property agents need to declare for the PIC grant in Singapore. However, not all of them know find out how to do the correct process for getting this PIC scheme from the IRAS. There are a number of steps that you need to do before your software can be approved.
  
  Naturally, you will have to pay a safety deposit and that is usually one month rent for annually of the settlement. That is the place your good religion deposit will likely be taken into account and will kind part or all of your security deposit. Anticipate to have a proportionate amount deducted out of your deposit if something is discovered to be damaged if you move out. It's best to you'll want to test the inventory drawn up by the owner, which can detail all objects in the property and their condition. If you happen to fail to notice any harm not already mentioned within the inventory before transferring in, you danger having to pay for it yourself.
  
  In case you are in search of an actual estate or Singapore property agent on-line, you simply should belief your intuition. It's because you do not know which agent is nice and which agent will not be. Carry out research on several brokers by looking out the internet. As soon as if you end up positive that a selected agent is dependable and reliable, you can choose to utilize his partnerise in finding you a home in Singapore. Most of the time, a property agent is taken into account to be good if he or she locations the contact data on his website. This may mean that the agent does not mind you calling them and asking them any questions relating to new properties in singapore in Singapore. After chatting with them you too can see them in their office after taking an appointment.
  
  Have handed an trade examination i.e Widespread Examination for House Brokers (CEHA) or Actual Property Agency (REA) examination, or equal; Exclusive brokers are extra keen to share listing information thus making certain the widest doable coverage inside the real estate community via Multiple Listings and Networking. Accepting a severe provide is simpler since your agent is totally conscious of all advertising activity related with your property. This reduces your having to check with a number of agents for some other offers. Price control is easily achieved. Paint work in good restore-discuss with your Property Marketing consultant if main works are still to be done. Softening in residential property prices proceed, led by 2.8 per cent decline within the index for Remainder of Central Region
  
  Once you place down the one per cent choice price to carry down a non-public property, it's important to accept its situation as it is whenever you move in – faulty air-con, choked rest room and all. Get round this by asking your agent to incorporate a ultimate inspection clause within the possibility-to-buy letter. HDB flat patrons routinely take pleasure in this security net. "There's a ultimate inspection of the property two days before the completion of all HDB transactions. If the air-con is defective, you can request the seller to repair it," says Kelvin.
  
  15.6.1 As the agent is an intermediary, generally, as soon as the principal and third party are introduced right into a contractual relationship, the agent drops out of the image, subject to any problems with remuneration or indemnification that he could have against the principal, and extra exceptionally, against the third occasion. Generally, agents are entitled to be indemnified for all liabilities reasonably incurred within the execution of the brokers´ authority.
  
  To achieve the very best outcomes, you must be always updated on market situations, including past transaction information and reliable projections. You could review and examine comparable homes that are currently available in the market, especially these which have been sold or not bought up to now six months. You'll be able to see a pattern of such report by clicking here It's essential to defend yourself in opposition to unscrupulous patrons. They are often very skilled in using highly unethical and manipulative techniques to try and lure you into a lure. That you must also protect your self, your loved ones, and personal belongings as you'll be serving many strangers in your home. Sign a listing itemizing of all of the objects provided by the proprietor, together with their situation. HSR Prime Recruiter 2010