[yhcode]

线性最小二乘基本模型：y = Ax+b
对于你的问题： y = [y1 y2 ... yn]';
A = [1 1/(x1^2) 1/(x1^4);1 1/(x2^2) 1/(x2^4);.....;1 1/(xn^2) 1/(xn^4);]
b为零均值残差。 x = [a b c]'
b的得到需要对y=Ax进行关于x进行一阶泰勒展开，之后得到b和新的A，就可以迭代求解了