polynomials - What is the Numerically Stable Way to Compute the Coefficients of a Quadratic Function from Three Points -

-

i compute coefficients a0, a1, , a2 of quadradic function (polynomial of degree 2) given 3 points (x0, y0), (x1, y1), , (x2, y2), yi = a0 + a1*xi + a2*xi*xi?

i have tried following 2 formulas, not impressed precision of output

final double x0mx1, x0mx2, x1mx2, t0, t1, t2; double a2, a1, a0;  x0mx1 = (x0 - x1); x0mx2 = (x0 - x2); x1mx2 = (x1 - x2);   // method 1 t0 = (y0 / (x0mx1 * x0mx2)); t1 = (y1 / (-x0mx1 * x1mx2)); t2 = (y2 / (x0mx2 * x1mx2));  a2 = (t0 + t1 + t2); a1 = -((t0 * (x1 + x2)) + (t1 * (x0 + x2)) + (t2 * (x0 + x1))); a0 = (t0 * x1 * x2) + (t1 * x0 * x2) + (t2 * x0 * x1);  // method 2 a2 = ((((y1 - y0) * (x0mx2)) + ((y2 - y0) * ((-x0mx1)))) /                 (((x0mx2) * ((x1 * x1) - (x0 * x0)))                 + (((-x0mx1)) * ((x2 * x2) - (x0 * x0))))); a1 = (((y1 - y0) - (a2 * ((x1 * x1) - (x0 * x0)))) / ((-x0mx1))); a0 = y0 - (a2 * x0 * x0) - (a1 * x0);

the results sort of fit, i.e., seem within +/- 1e-5 * max{ |a0'|, |a1'|, |a2'| } window of real solution a0', a1', , a2'.

is there better, more numerically stable way compute coefficients?

i using java, btw, although think not matter.

cheers, thomas.


Comments

Post a Comment

Popular posts from this blog

java - Static nested class instance -

-
this question has answer here: java inner class , static nested class 23 answers i came across in java, know static nested class why create instance of 'static'. isn't whole idea of 'static' use without instance? understand concept of inner classes, why (and frankly how possible to) create 'instance' of 'static' member @ ? why this: 1- outerclass.staticnestedclass nestedobject = new outerclass.staticnestedclass(); 2- nestedobject.anestedclassmethod(); and not this: 1- outerclass outerinstance=new outerclass(); 2- outerinstance.staticnestedclass.anestedclassmethod(); use on inner classes, keyword static indicates can access inner class without instance of outer class. for example: public class outer { public static class inner { /* code here. */ } } with construction, can create instance of inn
Read more

c# - Bluetooth LE CanUpdate Characteristic property -

-
i building xamarin.forms cross platform mobile app, uses monkey.robotics bluetoth low energy functionality. connecting mbed based implimentation of custom gatt service . in xamarin c#, triggers characteristic. canupdate property in monkey.robotics? this standard example c# based on: if (characteristic.canupdate) { characteristic.valueupdated += (s, e) => { debug.writeline("characteristic.valueupdated"); device.begininvokeonmainthread( () => { updatedisplay(characteristic); }); isbusy = false; // spin until first result received }; isbusy = true; characteristic.startupdates(); } this has been working, since changed own custom gatt service connecting to, canupdate property false. property , how triggered? me debugging gatt service code. thanks monkey.robotics open sour
Read more

JavaScript - Replace variable from string in all occurrences -

-
ok, know if have character '-' , want remove in places in string javascript, ... someword = someword.replace(/-/g, ''); but, when applying array of characters, s not working ... const badchars = ('\/:*?"<>|').split(''); let filename = title.replace(/ /g, '-').tolocalelowercase(); (let item = 0; item < badchars.length; item++) { // below not work global '/ /g' filename = filename.replace(/badchars[item]/g, ''); } any ideas? /badchars[item]/g looks badchars , literally, followed i , t , e , or m . if you're trying use character badchars[item] , you'll need use regexp constructor, , you'll need escape regex-specific characters. escaping regular expression has been well-covered . using that: filename = filename.replace(new regexp(regexp.quote(badchars[item]), 'g'), ''); but , don't want that. want character class: let filename = title.replac
Read more