Least-squares line fitting: Difference between revisions
From JSXGraph Wiki
A WASSERMANN (talk | contribs) No edit summary |
A WASSERMANN (talk | contribs) No edit summary |
||
| Line 23: | Line 23: | ||
n = p.length; | n = p.length; | ||
for (i=0;i<n;i++) { | for (i=0;i<n;i++) { | ||
M.push([ | M.push([p[i].X(), p[i].Y(), 1.0]); | ||
y.push( | y.push(p[i].X()*p[i].X() + p[i].Y()*p[i].Y()); | ||
} | } | ||
| Line 39: | Line 39: | ||
// Finally, we can read from the solution vector z the coordinates [xm, ym] of the center | // Finally, we can read from the solution vector z the coordinates [xm, ym] of the center | ||
// and the radius r and draw the circle. | // and the radius r and draw the circle. | ||
var zm = z[0]*0.5; | //var zm = z[0]*0.5; | ||
var xm = z[ | var xm = z[0]*0.5; | ||
var ym = z[ | var ym = z[1]*0.5; | ||
var r = Math.sqrt(z[ | var r = Math.sqrt(z[2]+zm*zm+xm*xm+ym*ym); | ||
brd.create('circle',[ [zm,xm,ym], r]); | brd.create('circle',[ [zm,xm,ym], r]); | ||
Revision as of 11:50, 7 November 2010
This little JXSGraph application finds the line - described by homogeneous coordinates [a,b,c] - that minimizes
- [math]\displaystyle{ \sum_{i=1}^n (ax_i+by_i+cz_i)^2. }[/math]
Coming soon...
