The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » sci.math.* » sci.math

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Least-squares scaling
Replies: 4   Last Post: Nov 13, 2012 6:57 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]

Posts: 289
Registered: 5/23/11
Least-squares scaling
Posted: Nov 11, 2012 1:55 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply



R in R^{d times n}
P in R^{d times n}, and
S in R^{d times d}, S symmetric positive semi-definite.

The problem is to find a matrix S such that the squared Frobenius norm

E = |SP - R|^2

is minimized. Geometrically, find a scaling which best relates the
paired vector sets P and Q. The E can be rewritten as

E = tr((SP - R)^T (SP - R))
= tr(P^T S^2 P) - 2tr(P^T SR) + tr(R^T R)
= tr(S^2 PP^T) - 2tr(SRP^T) + tr(RR^T).

Taking the first variation of E, with symmetric variations, and setting
it to zero gives that

SPP^T + PP^T S = RP^T + PR^T

holds in the minimum point. One can rearrange this to

(SPP^T - RP^T)^T = -(SPP^T - RP^T),

which says that SPP^T - RP^T is skew-symmetric. But I have no idea how
to make use of this fact. Anyone?


Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.