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 » Software » comp.soft-sys.matlab

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

Topic: PROOF? XS'*XS = eye ==> (XS\X0)' = X0'*XS
Replies: 1   Last Post: Mar 13, 2013 1:49 AM

Advanced Search

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

Posts: 9,822
Registered: 7/26/08
Re: PROOF? XS'*XS = eye ==> (XS\X0)' = X0'*XS
Posted: Mar 13, 2013 1:49 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

"Greg Heath" <> wrote in message <khoe76$bsd$>...
> help plsregress
> contains the line
> XL = (XS\X0)' = X0'*XS

If XS'*XS = eye(n), XS is full column-rank (the columns of XS is orthonormal)

(XS \ X0) is a not an overdetermined system, and there gives a unique solution Y of the least square solution:

Y = argmin |XS*Y - X0|^2

where |.| is the L2 norm. The Euler-Lagrange condition is:

XS'*XS*Y-XS'*X0 = 0.

Y = XS'*X0. (since XS'*XS = eye)

That means
Y = XS\X0 = XS'*X0

Transpose that you get your identity.


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.