用多项式拟合出的函数来计算f(x)值怎么就不对呢

seeyou8

这样两列数据,多项式拟合一个函数y = 44195x4 - 132962x3 + 149979x2 - 75172x + 14127,然后我取一个数据计算f(0.728803),值怎么不和原来的值接近呢.拟合挺好的啊.，这是什么原因啊？

0.728803	0.59077
0.748929	0.628106
0.774194	0.692193
0.786858	0.745445
0.789025	0.758609
0.789691	0.763105
0.791325	0.775284
0.792127	0.782015
0.792904	0.789121
0.796175	0.828962
0.797693	0.856583
0.797903	0.861061
0.798057	0.864457
0.803538	1.001588
0.804494	1.021224
0.804628	1.023867
0.806469	1.057236
0.809097	1.097809

[ 本帖最后由 eight 于 2008-5-14 15:31 编辑 ]

wwbeyondww

怎么不接近啊?蛮好啊.
[a,S]=polyfit(x,y,4);
[ye,delta]=polyval(a,0.728803,S);
ye

ye =

   0.58990327400716

seeyou8

您好,我拟合出的函数y = 44195x4 - 132962x3 + 149979x2 - 75172x + 14127
能不能直接用这个函数计算值,比如
y=44195*0.728803^4-132962*0.728803^3+149979*0.728803^2-75172*0.728803+14127

结果就是1.3718了.

seeyou8

不好意思,我知道错误了,大概是拟合出来的系数小数位不够,所以直接用公式计算值会有误差.,精度不够.

wwbeyondww

估计是你取的位数不够.

a =

  1.0e+005 *

   0.44198651960080  -1.32972123361196   1.49990439386474  -0.75178293998984   0.14127617613639
b=0.728803;
a(1)*b^4+a(2)*b^3+a(3)*b^2+a(4)*b+a(5)

ans =

   0.58990327401079