最小二乘法拟合曲线C++编程

[lcymd]

i 0 1 2 3 4
xi 1.00 1.25 1.50 1.75 2.00
yi 5.10 5.79 6.53 7.45 8.46
-
yi 1.629 1.756 1.876 2.008 2.135

-
lnyi= yi y=ae^(bx)
急 求 C++ 程序 谢谢各位了!!

[ToBeTough]

不懂

[muyejingfeng100]

这个应该不是很难吧……只要懂曲线按按拟合就可以了。。

[muyejingfeng100]

如果将上面的拟合成y=ae^bx。只要将两边取对数lny=bx+lna再将lny=y',lna=a'就可以化为y'=bx+a'。然后由
xi 1.00 1.25 1.50 1.75 2.00
y' 1.629 1.756 1.876 2.008 2.135
利用p1=1,p2=x作基，权为1。
则有(p1,p1)=(1+1+1+1+1)=5 (p1,p2)=(p2,p1)=(x0+x1+x2+x3+x4)=(1.00+1.25+1.50+1.75+2.00)=7.5
(p2,p2)=(x0^2+x1^2+x2^2+x3^2+x4^2)=(1.00^2+1.25^2+1.50^2+1.75^2+2.00^2)=11.875
(y',p1)=(y'0+y'1+y'2+y'3+y'4)=(1.629+1.756+1.876+2.008+2.135)=9.404
(y',p2)=(y'0*x0+y'1*x1+y'2*x2+y'3+x3+y'4*x4)=(1.629*1.00+1.756*1.25+1.876*1.50+2.008*1.75+2.135*2.00)=14.422
则有：
5a'+7.5b=9.404
7.5a'+11.875b=14.422
根据上面方程就可以解出a',b。最后回代就可以了。
你根据上面的思路编程就可以了。。

[zhaohongbo83]

你按照最小二乘法的公式，写一个函数就可以了
然后自己调用。