实验32用待定系数法求插值多项式 求通过指定的几个点的多项式方程 以下仅对三次和四次多项式的情况进行讨论 我们知道对给定的点D=(x1y1),E=(x2,y2)和F=(x3y3),求出拟合它们的函数 p(x)=ax2+bx+c的系数是非常麻烦的。求解这个问题需要解下列的线性方程组 axl+bxl+cxl =y1 ax2 + b x2 +Cx2 =y2 ax3+bx3+Cx3=y3 此问题徒手来做也是很容易,但是如果我们给定四个点或者给定六个点来拟合三次 或五次多项式便没有那么容易了.我们使用 Mathcad的求解模块来解决这个问题 定义数据点向量,x代表点的横坐标,y代表点的纵坐标: 4 这个点集看上去可以用三次多项式拟合,我们用函数c(x)=ax2+bx2+cx+d 建立方程组和求解模块为求解模块定义初始值: d:=1 Given d a()3+b()2+cx+d= (*3)+b(*3) + cx,+d=ya x 3 ( ) 3 × b x 3 ( ) 2 + × c x 3 + × + d y 3 = a x 2 ( ) 3 × b x 2 ( ) 2 + × c x 2 + × + d y 2 = a x 1 ( ) 3 × b x 1 ( ) 2 + × c x 1 + × + d y 1 = a x 0 ( ) 3 × b x 0 ( ) 2 + × c x 0 + × + d y 0 = Given a := 1 b := 1 c := 1 d := 1 这个点集看上去可以用三次多项式拟合, 我们用函数 c(x) a x 3 × b x 2 = + × + c×x + d 建立方程组和求解模块. 为求解模块定义初始值: 0 1 2 3 4 5 5 10 y x y i := 0.. 3 1 5 2 9 æ ç ç ç ç è ö ÷ ÷ ÷ ÷ ø x := 1 2 3 4 æ ç ç ç ç è ö ÷ ÷ ÷ ÷ ø := 定义数据点向量, x代表点的横坐标, y代表点的纵坐标: 此问题徒手来做也是很容易, 但是如果我们给定四个点或者给定六个点来拟合三次 或五次多项式便没有那么容易了. 我们使用Mathcad的求解模块来解决这个问题. a×x3 + b×x3 + c×x3 =y3 a×x2 + b×x2 + c×x2 =y2 a×x1 + b×x1 + c×x1 =y1 以下仅对三次和四次多项式的情况进行讨论. 我们知道, 对给定的点D=(x1,y1), E=(x2,y2) 和 F=(x3,y3), 求出拟合它们的函数 p(x) a x 2 = × + b×x + c 的系数是非常麻烦的。求解这个问题需要解下列的线性方程组: 求通过指定的几个点的多项式方程 实验32 用待定系数法求插值多项式 1