Hermite插值条件下的三点曲线拟合程序设计.docVIP

本文已经被《湖南工程学院学报》录用 Hermite插值条件下的三点曲线拟合程序设计 李庆福、张光远、颜宁生 1）北京服装学院基础教学部 2）北京服装学院商学院 北京100029 摘要:利用插值法和最小二乘法解决Hermite插值条件下的三点曲线拟合问题，并通过刀具磨损问题实例检验了Hermite插值条件下的三点曲线拟合效果，最后对Hermite插值条件下的三点曲线拟合问题作了程序设计。 关键词: 最小二乘法；Hermite插值条件； 曲线拟合 中图分类号 O241.5 文献码 A 1 引言 在文献中提出了一种将插值与拟合相结合的方法,文献在此基础上讨论了一个插 值点和三个拟合点的问题,得到了拟合曲线,其中为插值点 (可以将实际问题中最重要的一个数据作为插值点), 为最小二乘拟合曲线的待定的参 数.本文在文献中关于一个插值点的基础上，增加了在这个点上的导数要求,我们将要求一 条次数尽可能低的曲线经过,且在处满足导数条件 ， 对另外三点,j=1,2,3，进行曲线拟合的问题称为Hermite插值条件下的三点曲线拟合， 不难看出，在Hermite插值条件下的三点曲线拟合模型是： 其中 2 五次曲线拟合模型中参数的确定 令 将看成自变量和的一个二元函数，那么问题就可以归结为求函数在哪些点处取得最小值。 解方程组： 即 设 则 3 VC辅助算法 程序使用VC编制界面，通过在VC下调用Matcom进行绘图。将两者结合采用MATCOM转化法， 将Windows中的系统文件v4501v.dll与ago4501.dll以及MATCOM中的文件matlib.h与 v4501.lib加入VC++工程目录，然后在程序中包含文件matlib.h和链接文件v4501.lib。主要 代码如下： //matcom接口初始化 initM(MATCOM_VERSION); //设置绘图句柄 Mm h = (BR(0.0),100,0.0,100); CWnd *p1=NULL; p1=(CWnd *)GetDlgItem(IDC_PIC); Mm plothandle=winaxes(p1-m_hWnd); //设置坐标 Mm pos=(BR(0.0),0.0,300,300); set(plothandle,TM(RealPosition),pos); //根据上文公式计算拟合方程系数a和b（其中m_pt[i][1]、m_pt[i][2]分别为数据点的x、y坐标，count为数据点数目） for(int i = 1; i (count-1); i++) { M = m_pt[i][1]*(m_pt[i][1]-m_pt[0][1])*(m_pt[i][1]-m_pt[datacount-1][1])*(m_pt[i][1]-m_pt[0][1])*(m_pt[i][1]-m_pt[datacount-1][1]); N = (m_pt[i][1]-m_pt[0][1])*(m_pt[i][1]-m_pt[datacount-1][1])*(m_pt[i][1]-m_pt[0][1])*(m_pt[i][1]-m_pt[datacount-1][1]); L = m_pt[i][2]-m_pt[0][2]*(1-2*(m_pt[i][1]-m_pt[0][1])/(m_pt[0][1]-m_pt[datacount-1][1])) *((m_pt[i][1]-m_pt[datacount-1][1])*(m_pt[i][1]-m_pt[datacount-1][1]))/((m_pt[0][1]-m_pt[datacount-1][1])*(m_pt[0][1]-m_pt[datacount-1][1])) -m_y0*(m_pt[i][1]-m_pt[0][1]) *((m_pt[i][1]-m_pt[datacount-1][1])*(m_pt[i][1]-m_pt[datacount-1][1]))/((m_pt[0][1]-m_pt[datacount-1][1])*(m_pt[0][1]-m_pt[datacount-1][1])) -m_pt[datacount-1][2]*(1-2*(m_pt[i][1]-m_pt[datacount-1][1])/(m_pt[datacount-1][