|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
大家可以试一下
A=[500 528.56 557.12 592.82 617.14 649.98 664.26 682.82 700 704.28];
B=[0.47712 0.60206 0.77815 1 1.30103 1.77815 2 2.30103 3 3.30103]; pp=spline(A,B); [breaks,coefs,npolys,ncoefs]=unmkpp(pp)
y=ppval(pp,500);
得到得系数,代入拟合后得到得三次多项式中f(x)=a0*x^3+a1*x^2+a2*x+a3。
按道理f(500)应该等于y,但是就是不相等,不知道出现什么问题?难道是漏洞? 取其他值也不相等!
coefs =
-0.0000 0.0001 0.0019 0.4771
-0.0000 0.0000 0.0061 0.6021
0.0000 -0.0001 0.0055 0.7782
-0.0000 0.0002 0.0096 1.0000
0.0000 0.0000 0.0138 1.3010
-0.0000 0.0001 0.0158 1.7781
0.0000 -0.0002 0.0141 2.0000
0.0000 0.0007 0.0236 2.3010
0.0000 0.0016 0.0631 3.0000
f(500)= -0.0000*500^3+0.0001*500^2+ 0.0019 *500 +0.4771=26.4271并不等于y=ppval(pp,500)=0.47712
[ 本帖最后由 eight 于 2008-1-8 15:09 编辑 ] |
|