Ask your own question, for FREE!
Linear Algebra 21 Online
OpenStudy (konradzuse):

Given the linear system Ax = b with (a) Find all least squares solutions of the linear system Ax = b. (b) What is the normal system for this least squares problem? (c) Compute the least squares error to three decimal places.

OpenStudy (konradzuse):

\[A =\left[\begin{matrix}1 & 3 \\ 3 & 9\\ -2 & -6\end{matrix}\right]\] \[B = \left(\begin{matrix}1 \\ 1\\ 0\end{matrix}\right)\]

OpenStudy (konradzuse):

@UnkleRhaukus @Zarkon

OpenStudy (konradzuse):

Since my work was deleted I will explain that I got up to the end where I need to solve, I am confused exactly how to solve the matrix [x1x2] col matrix = [4...12] which was my col matrix answer... Flipping out trying to get this done by midnight.... :(

OpenStudy (konradzuse):

There was also a trick someone taught me that if I take the A * A^-1 I can use that to find my answer; however i got 0 for my determinant so........ Idk how that works, everything = 0? I don't think so...... Maybe?

OpenStudy (konradzuse):

nope, it's correct.....

OpenStudy (konradzuse):

Tk have you come to save me? :P

OpenStudy (tkhunny):

Ax = b Introduce the Transpose \(A^{T}Ax = A^{T}b\) Hope that \(A^{T}A\) is non-singular and calculate it's inverse: \(\left(A^{T}A\right)^{-1}\left(A^{T}A\right)x = \left(A^{T}A\right)^{-1}A^{T}b\) \(Ix = \left(A^{T}A\right)^{-1}A^{T}b\)

OpenStudy (konradzuse):

Oh were you the one who showed me that trick with the -1? When taking the determinant I get 0... What happens then? Is it all 0? lol

OpenStudy (tkhunny):

Determinant? It's an inverse.

OpenStudy (tkhunny):

Unfortunately, your \(\left(A^{T}A\right)^{-1}\) is singular. We should have seen this from the original data. Column 2 is rather obviuosly twice Column 1. Where does that leave us?

OpenStudy (konradzuse):

I thought taking the inverse is the same as taking the determinant?

OpenStudy (konradzuse):

I guess I konfused myself again :p.

OpenStudy (konradzuse):

1/ad-bc [d -b -c a]?

OpenStudy (konradzuse):

:P

OpenStudy (tkhunny):

No, the Determinat CAN be used as part of finding an Inverse, depending on how you find an Inverse. However, in this case, there is no Inverse. The determinant will not help you find it.

OpenStudy (konradzuse):

ic, so what happens since there is none?

OpenStudy (konradzuse):

I got 0 so that makes sense.

OpenStudy (tkhunny):

Well, just for practice, can you write the Normal Equations? I'll get you started: \(\sum x_{1} + 3\cdot\sum x_{2} = \sum y\) You write the other two.

OpenStudy (konradzuse):

3x1 +9x2= y

OpenStudy (konradzuse):

-2x1 - 6x2

OpenStudy (tkhunny):

Whoops! Using the model \(I = \left(a\cdot x_{1} + b\cdot x_{2} - y\right)^{2}\) Forget that first thing. You did exactly what I told you and that was incorrect. We're trying to FIND the coefficients, not just invent them! Find \(\dfrac{\partial I}{\partial a}\) and \(\dfrac{\partial I}{\partial b}\)

OpenStudy (konradzuse):

I'm soo konfused.......

OpenStudy (konradzuse):

I am able to solve everything up until the end where it asks for the final part, that's where IU'm confused. I'll show yu what my book explains.

OpenStudy (konradzuse):

OpenStudy (zarkon):

for your first question just solve the system (without inverses). you will get an infinite number of solutions. which is fine

OpenStudy (konradzuse):

yeah I figured that out 2/7-3*t is x1 and x2 is t.

OpenStudy (konradzuse):

part b we already solved the normal eq, now I'm working on the error part 3...

OpenStudy (konradzuse):

I'm not sure how to get the decimal answer. My matrix comes out as Vector(3, {(1) = 5/7, (2) = 1/7, (3) = 4/7})

OpenStudy (zarkon):

are you taking the error to be \[\|Ax-b\|\]

OpenStudy (konradzuse):

first it shows the error vector, then shows that.

OpenStudy (konradzuse):

is that asking for the norm?

OpenStudy (zarkon):

then just compute it...you have all the pieces yes

OpenStudy (konradzuse):

does that error play into it, or what exactly am I computing the norm with?

OpenStudy (konradzuse):

error vector*

OpenStudy (zarkon):

you are computing \(\|Ax-b\|\) where for any vector \(v=<a,b,c>\) \[\|v\|=\sqrt{a^2+b^2+c^2}\]

OpenStudy (konradzuse):

SO I guess we do use that error vector then lol.

OpenStudy (konradzuse):

(1/7)*sqrt(42)

OpenStudy (konradzuse):

0.9258

OpenStudy (zarkon):

seems fine

OpenStudy (konradzuse):

Thanks for the help :).

Can't find your answer? Make a FREE account and ask your own questions, OR help others and earn volunteer hours!

Join our real-time social learning platform and learn together with your friends!
Can't find your answer? Make a FREE account and ask your own questions, OR help others and earn volunteer hours!

Join our real-time social learning platform and learn together with your friends!