最小二乘法算多次项的解(求指教)

sintai

请问可以提供M文件吗?我刚学,遇到这个问题.
不知道变程序,大概要求拟合的精度到1.0E-10~15的精度.
数据有60000个点,是轮廓测定仪器的实际点坐标数据,需要求
R,K,A1,A2,A4,A6......,点的数据是*.TXT文件格式,也可以作其他格式.

用最小二乘法算多次项的解,函数的方程和上面的一样的:y=x^2/R/(1+sqrt(1-(1-K)*x^2/R^2))+A*x^4+B*x^6+C*x^8+D*x^10+E*x^12

附件是简单的点数据文件,谢谢.

[ 本帖最后由 bainhome 于 2007-1-26 01:17 编辑 ]

shamohu

这是针对你数据的答案：

均方差(RMSE): 8.08420524279771E-5
残差平方和(SSE): 2.48346622749176E-7
相关系数(R): 0.983099371625185
相关系数之平方(R^2): 0.966484374489834

r            227.342858670186
k            -14882.3692110109
a            -0.00596362353181927
b            0.00677374903347943
c            -0.00399395071700348
d            0.00111622672317863
e            -0.000119864513133811

sintai

非常感谢,可以告诉我怎么计算的吗?程序?谢谢!~

shamohu

是1stOpt的功劳！试试吧，本论坛有1stOpt的介绍。

sintai

1stOpt计算出来的结果于我的实际数据相差好远,如R,K等等完全不对,是否计算的方式有问题?