【求助】偶次多项式的曲线拟合

dirac · 发表于 2006-4-8 09:59

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

对读取的面型轮廓进行曲线拟合，已知的是该轮廓曲线是偶次多项式的形式，最高阶次为6，也就是说可以表达为：f(x)=a*x^6+b*x^4+c*x^2，顶点在原点，故没有常数项。 但是读取的轮廓与实际轮廓肯定不完全相等，而是大小成比例；而且所读取的坐标并非完全对称的。也就是说读取到的曲线多项式应该为： F(x)=A*(E*x+F)^6+B*(E*x+F)^4+C*(E*x+F)^2+D MATLAB的多项式拟合后的表达式是从高阶到低阶的依次展开式，而要想得到我所需要的上述形式的话，可以利用各阶次系数相等来得到一组非线性方程。但是非线性方程组的求解需要对其赋初值x0，因为x0的选择对结果是有很大影响的。所以，请问，在这种情况下，如何选择x0比较合适呢？ 或者，还有别的办法可以得到我所需要的表达式吗？ 谢谢！

happy · 发表于 2006-4-8 14:56

可以取x1=x^2 这样你的拟合方程可以变成3次方程，求解就方便多了

dirac · 发表于 2006-4-8 19:19

对啊，我怎么没想到。谢谢啦！

dirac · 发表于 2006-4-8 21:17

再请教，取x1＝x^2以后，拟合的曲线好像非常失真啊，与数据非常不符。怎么回事？ 帮忙看看： xdata=[linspace(168,298,14) linspace(318,418,6) 518 818 linspace(888,988,6) linspace(998,1158,17)]; ydata=[749 677 605 547 494 444 404 364 330 300 272 250 228 210 178 155 139 127 118 112 102 102 108 112 119 128 140 156 168 179 194 210 228 250 274 303 333 367 407 447 499 552 612 680 749]; polytool(xdata.^2,ydata,3) 这样得到的曲线有问题。为什么？

happy · 发表于 2006-4-9 14:46

你这个数据最好能够采用分段拟合

dirac · 发表于 2006-4-9 18:44

我的数据就是直接从一条6阶曲线读出来的啊 曲线方程：y=3*x^6+8*x^4+11*x^2 不需要分段拟合吧

账号		自动登录	找回密码
密码			我要加入

【求助】偶次多项式的曲线拟合

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（dirac）【求助】偶次多项式的曲线拟合

回复：（happy）回复：（dirac）【求助】偶次多项式的...

回复：（dirac）【求助】偶次多项式的曲线拟合

回复：（happy）回复：（dirac）【求助】偶次多项式的...