Linear algebra least squares solution

[SpiffyEh]

Homework Statement

Suppose you have a set S of three points in R^2,
S1 = {(1, 12), (2, 15), (3, 16)}
S2 = {(1, 12), (1, 15), (3, 16)}
S3 = {(1, 12), (2, 15), (2, 15)}

which you seek to interpolate with the quadratic polynomial p(t) = a_0 + a_1t + a_2t^2.
Problem: Least-Squares Approximation. Find the least-squares solution for S3 which we previously found to have no solution.

Homework Equations

The Attempt at a Solution

I understand the idea of least squares solutions I think, if I have a matrix and a solution but I'm confused about this one. I'm just not seeing where to start this. If someone could please get me going that would be great. Thanks

 

[hgfalling]

Recall that in least squares analysis, our goal is the set of coefficients that you have identified as a_i.

So we could write down the equations of this regression as:

a_0 + a_1 t_1 + a_2 t_1^2 = y_1
a_0 + a_1 t_2 + a_2 t_2^2 = y_2
a_0 + a_1 t_3 + a_2 t_3^2 = y_3

This is the same as the matrix equation:

<br /> \left(\begin{array}{ccc}<br /> 1 &amp; t_1 &amp; t_1^2 \\ 1 &amp; t_2 &amp; t_2^2 \\ 1 &amp; t_3 &amp; t_3^2<br /> \end{array}\right)<br /> \left(\begin{array}{c}<br /> a_0\\a_1\\a_2<br /> \end{array}\right) <br /> = <br /> \left(\begin{array}{c}<br /> y_1\\y_2\\y_3<br /> \end{array}\right) <br />

Now you can use your normal least squares technique to solve this for \vec a.

 

[SpiffyEh]

I don't see how to set s3 up in this way though, that's what I'm really having issues with.

 

[hgfalling]

Well, the a vector is going to be an unknown. You need to plug in for t and for y.

Your data has (independent, dependent) data, right? So plug in the data to the equation there.

 

[SpiffyEh]

umm... since S3 = {(1, 12), (2, 15), (2, 15)}

\left(\begin{array}{ccc}1 &amp; 1 &amp; 12 \\ 1 &amp; 2 &amp; 12 \\ 1 &amp; 2 &amp; 15\end{array}\right)\left(\begin{array}{c}a_0\\a_1\\a_2\end{array}\right) = \left(\begin{array}{c}y_1\\y_2\\y_3\end{array}\right)

Is that right? Then what do I do for y? I don't think I'm getting it

 

[Dickfore]

Recall that in least squares analysis, our goal is the set of coefficients that you have identified as a_i.

So we could write down the equations of this regression as:

a_0 + a_1 t_1 + a_2 t_1^2 = y_1
a_0 + a_1 t_2 + a_2 t_2^2 = y_2
a_0 + a_1 t_3 + a_2 t_3^2 = y_3

This is the same as the matrix equation:

<br /> \left(\begin{array}{ccc}<br /> 1 &amp; t_1 &amp; t_1^2 \\ 1 &amp; t_2 &amp; t_2^2 \\ 1 &amp; t_3 &amp; t_3^2<br /> \end{array}\right)<br /> \left(\begin{array}{c}<br /> a_0\\a_1\\a_2<br /> \end{array}\right) <br /> = <br /> \left(\begin{array}{c}<br /> y_1\\y_2\\y_3<br /> \end{array}\right) <br />

Now you can use your normal least squares technique to solve this for \vec a.

This is not least squares method. There is no more than one parabola passing through three points, so there is no fitting involved.

 

[Mark44]

This is not least squares method. There is no more than one parabola passing through three points, so there is no fitting involved.

But S3 has only two distinct points.

This problem is not stated very clearly, IMO.

Suppose you have a set S of three points in R^2,
S1 = {(1, 12), (2, 15), (3, 16)}
S2 = {(1, 12), (1, 15), (3, 16)}
S3 = {(1, 12), (2, 15), (2, 15)}

I don't see any set S. I see three sets, S1, S2, and S3, the first two of which have three points, and the last of which has two points.

 

[SpiffyEh]

I'm told to use s3 which is why I'm so confused about how to set up a matrix for it and a vector b in Ax = b. I understand where to go from there.

 

[Mark44]

Can you give us the exact wording of this problem? As you have posted it, there seems to be some incorrect and/or missing information.

 

[SpiffyEh]

I attached a screen shot of the exact problem. I don't think I missed any information.

 

[Mark44]

No, it looks like you copied the information faithfully, but whoever wrote the problem wasn't very clear. It says that S is a set with three points, then gives three sets, S1, S2, and S3, each with (sort of) three points.

I would be inclined to ask the instructor for clarification.

 

[Dickfore]

The point is that for every set S_{i}, (i=1,2,3) there are three points in the plane and you can draw at most one interpolating parabola y = a_{0} + a_{1} \, x + a_{2} \, x^{2}, i.e. one that passes through each of the three points according to the system:

<br /> \left(\begin{array}{ccc}<br /> 1 &amp; x_{1} &amp; x^{2}_{1} \\<br /> <br /> 1 &amp; x_{2} &amp; x^{2}_{2} \\<br /> <br /> 1 &amp; x_{3} &amp; x^{2}_{3}<br /> \end{array}\right) \cdot \left(\begin{array}{c}<br /> a_{0} \\ a_{1} \\ a_{2}<br /> \end{array}\right) = \left(\begin{array}{c}<br /> y_{1} &amp; y_{2} &amp; y_{3}<br /> \end{array}\right)<br />

However, this is not least squares fitting. This is quadratic interpolation. If the determinant of the system is zero, there is no solution for this linear system. I think the problem says that for the set S₃ this is the case.

In fitting, you don't require that the fitting curve passes through all of the experimental points (it can even not pass through any of them). You simply require that the coeffiicients are chosen such that the sum of the squares of the errors:

<br /> S(a_{0}, a_{1}, a_{2}) = \sum_{i = 1}^{3}{(a_{0} + a_{1} \, x_{i} + a_{2} \, x^{2}_{i} - y_{i})^{2}}<br />

be minimal. For this, you have the conditions:

<br /> \frac{\partial S}{\partial a_{0}} = \frac{\partial S}{\partial a_{1}} = \frac{\partial S}{\partial a_{2}} = 0<br />

By evaluating the derivatives and simplifying, you should obtain another linear system for the fitting coefficients. Try to find what the system looks like (in general) and post it here so that we can verify that your work is correct.

 

[SpiffyEh]

I'm confused. I don't know what to do with s3. If you could explain in more detail that would really help. I also emailed my professor but he hasn't gotten back to me and this is a practice problem and my exam is tomorrow so I'm hoping I can understand it before then.

 

[SpiffyEh]

ok, I found out from the professor that
A = \left(\begin{array}{ccc}1 &amp; 1 &amp; 1 \\1 &amp; 1 &amp; 1 \\1 &amp; 3 &amp; 9\end{array}\right)
and b = \left(\begin{array}{ccc}12 \\ 15 \\16\end{array}\right)

So, I went through and did A^TA and A^Tb and augmented the matrix and row reduced to solve for x
I got

x =\left(\begin{array}{ccc}391/60 \\ 5/4\\0\end{array}\right) + x_3 \left(\begin{array}{ccc}3\\ -4\\1\end{array}\right)

Can someone please check my work on this? I'm worried about the 391/60, but I went through the work and I don't see an error.