Example 8.  Use the linear least squares method to find the polynomial curve fit of degree = 3  for the points [Graphics:Images/LeastSqPolyMod_gr_120.gif].

Solution 8.

Enter the n points.

[Graphics:../Images/LeastSqPolyMod_gr_121.gif]


[Graphics:../Images/LeastSqPolyMod_gr_122.gif]

[Graphics:../Images/LeastSqPolyMod_gr_123.gif]

[Graphics:../Images/LeastSqPolyMod_gr_124.gif]

[Graphics:../Images/LeastSqPolyMod_gr_125.gif]

Enter the degree of the polynomial.

[Graphics:../Images/LeastSqPolyMod_gr_126.gif]
[Graphics:../Images/LeastSqPolyMod_gr_127.gif]
[Graphics:../Images/LeastSqPolyMod_gr_128.gif]

Construct the matrices  [Graphics:../Images/LeastSqPolyMod_gr_129.gif], and the vector  [Graphics:../Images/LeastSqPolyMod_gr_130.gif]

[Graphics:../Images/LeastSqPolyMod_gr_131.gif]



[Graphics:../Images/LeastSqPolyMod_gr_132.gif]
[Graphics:../Images/LeastSqPolyMod_gr_133.gif]
[Graphics:../Images/LeastSqPolyMod_gr_134.gif]
[Graphics:../Images/LeastSqPolyMod_gr_135.gif]
[Graphics:../Images/LeastSqPolyMod_gr_136.gif]
[Graphics:../Images/LeastSqPolyMod_gr_137.gif]

Solve the linear system  [Graphics:../Images/LeastSqPolyMod_gr_138.gif].

[Graphics:../Images/LeastSqPolyMod_gr_139.gif]


[Graphics:../Images/LeastSqPolyMod_gr_140.gif]
[Graphics:../Images/LeastSqPolyMod_gr_141.gif]
[Graphics:../Images/LeastSqPolyMod_gr_142.gif]
[Graphics:../Images/LeastSqPolyMod_gr_143.gif]

Construct the polynomial.

[Graphics:../Images/LeastSqPolyMod_gr_144.gif]


[Graphics:../Images/LeastSqPolyMod_gr_145.gif]

We are done.

Aside.  We can compare this with the result obtained by using Mathematica's procedure "Fit" to find the polynomial.

[Graphics:../Images/LeastSqPolyMod_gr_146.gif]

[Graphics:../Images/LeastSqPolyMod_gr_147.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004