fpoly.m

该算法是bimoment非线性归一化的改进

function  [a,b,c] = fpoly(X,Y)
syms a b c;
x0 = X(1); x1 = X(2); x2 = X(3);
y0 = Y(1); y1 = Y(2); y2 = Y(3);
f0 = a * x0.^2 + b * x0 + c -y0;
f1 = a * x1.^2 + b * x1 + c -y1;
f2 = a * x2.^2 + b * x2 + c -y2;
[a,b,c] = solve(f0,f1,f2,'a','b','c');
a=double(a);
b=double(b);
c=double(c);