第四章多项式及插值.pptVIP

二、 多项式计算 1. 多项式求根 求多项式 p(x)的根的函数是roots(P)，这里，P是 p(x)的系数向量，该函数返回方程 p(x)=0 的全部根 (含重根，复根)。 2. 多项式求值 求多项式 p(x)在某点或某些点的函数值的函数是 polyval(P, x)。若x为一数值，则求多项式在该点的 值；若x为向量或矩阵，则对向量或矩阵中的每个元 素求其多项式的值。 例1 已知一个多项式 (1)计算f (x)=0 的全部根。 (2)由方程f (x)=0的根构造一个多项式 g(x)，并与 f (x) 进行对比。 (3)计算f (5)、f (7.8)、f (9.6)、f (12.3)的值。 3. 多项式的四则运算 (1)多项式的加减法 例2 设有两个多项式，计算： (1)求f(x)+g(x)、f(x)-g(x)。 (2)求f(x)·g(x)、f(x)/g(x)。 4. 多项式的微分与积分 （1）对多项式求导数的函数是： p=polyder(P) 求多项式P的导函数 p=polyder(P,Q) 求P*Q的导函数 [p,q]=polyder(P,Q) 求P/Q的导函数，导函数的 分子存入p，分母存入q。 （2）对多项式的积分函数： d=poly_itg(c) d是多项式c积分后的系数，但 不包括积分常数 例3 求有理分式的导数。 为计算Lagrange插值多项式的一阶导数，可用polyder 函数将 p 的每一行转换为一阶导数的系数数组。 x = [1.1, 2.3, 3.9, 5.1]; y = [3.887, 4.276, 4.651, 2.117]; xi = [2.101 ,4.234]; np = length(x); p=shape_pw(x); s=0; for i=1:np s = s+polyder(p(i,:)).*y(i); end yi=polyval(s,xi) 结果为： yi＝ 0.6292 -1.4004 四、Chebyshev多项式 等距Lagrange插值的误差在中间部分最小，越靠近端点处误差越大，降低误差的一种方法是改变插值点的分布，即增大定义域中间部分的插值步长，同时减少定义域两端的插值步长，但是最优的插值点分布取决于插值多项式的目的。 若目的是近似一个函数，则最优插值点是Chebyshev多项式的零点，因为 （1）此时L(x)分布是最平坦的分布； （2）误差值不会像等距插值那样随不同插值区间 而变化。 Chebyshev多项式： 当 时，Chebyshev多项式的图形如下 幂级数形式Chebyshev多项式的系数可由函数 Cheby_ pw计算 function pn = Cheby_pw(n) pbb=[1]; if n==0, pn=pbb; break; end pb=[1 0]; if n==1, pn=pb; break; end for i=2:n; pn= 2*[pb,0] - [0, 0, pbb] ; pbb=pb; pb=pn; end 调用格式： p=Cheby_pw(n) n为多项式阶数，p是系数的行数组 Chebyshev多项式的根系数可由 roots命令计算 roots(Cheby_ pw(n)) 若令 ，则由 得 有 k 个根，都在[-1,1]，也可用下面公式计算： 若插值区间为[a,b]，可建立映射 下图比较了区间 上9个Chebyshev点和等距 点条件下L(x)的分布 五、Lobatto多项式 Lobatto点是Chebyshev多项式的零点加上 x ＝-1和 x＝1 对区间[a,b]，k +1个Lobatto点可由 k 阶Chebyshev 多项式得到 六、Legendre多项式 Legendre多项式： n 阶Legendre多项式的系数可由函数Legen_pw计算： function pn = Legen_pw(n) pbb=[1]; if n==0, pn=pbb; break; end pb=[1 0]; if n==1, pn=pb; bre