Matlab第六节多项式运算.pptVIP

数学实验 多项式及相关运算 主要内容 多项式基本运算 多项式插值拟合 Matlab 多项式基本运算 多项式构造 Matlab 中多项式的表示方法 多项式构造 多项式四则运算 多项式加减运算 多项式四则运算 多项式乘法运算： 多项式的求导 有理多项式运算 部分分式展开与结合： 多项式的值 多项式的值 多项式的零点 多项式插值、拟合 多项式拟合 最小二乘拟合：polyfit 拟合举例 【例】 某种合金中的主要成分为A,B两种金属，经过试验发现：这两种金属成分之和x与合金的膨胀系数y有如下关系，建立描述这种关系的数学表达式. 多项式插值 ‘nearest’(最邻近插值):将插值点x处的函数值y置为与之最接近的已知点xi处的值，即为yi。 ‘linear’(线性插值):用直线连接相邻的已知数据点（ xi ，yi），插值点x处的函数值y由该直线确定。 ‘spline’(三次样条插值):用三次函数连接相邻的已知数据点（ xi ，yi），插值函数在已知点处连续并且它的一阶导数和二阶导数也连续。 ‘cubic’(立方插值):用于实现分段光滑hermite插值，并保留数据的单调性和形状。 多项式插值 实验内容 1、求多项式x4+8x3-10的根，并根据已经求出的根建立起该多项式，并与原多项式进行对比。 2、分别求f(x)=x3-10x+2在A=[1 2 3;4 5 6;2 3 5] 处的代数多项式和矩阵多项式的值. 3、课本：P130：3,4,5；8；12 * * 在 Matlab 中，n 次多项式是用一个长度为 n+1的向量来表示，缺少的幂次项系数为 0。 在 Matlab中表示为向量： 注：系数中的零不能省！ 例： 多项式的符号形式：poly2sym（p,’x’） poly2sym([2,-1,0,3]) poly2sym([2,-1,0,3],’x’) 特征多项式或指定根的多项式：poly（A） 若A是方阵，返回A的特征多项式；若A是向量，返回以A中元素为根的首1多项式。 A=[0,1,1;1,0,1;1,1,1]; poly(A), x=[1 2 0]; poly(x) Matlab 没有提供专门进行多项式加减运算的函数，事实上，多项式的加减就是其所对应的系数向量的加减运算 对于次数相同的多项式，可以直接对其系数向量进行加减运算； 如果两个多项式次数不同，则应该把低次多项式中系数不足的高次项用 0 补足，然后进行加减运算。 例： k = conv(p,q) 例：计算多项式 和 的乘积 p=[2,-1,0,3]; q=[2,1]; k=conv(p,q); 多项式除法运算： [k,r] = deconv(p,q) 其中 k 返回的是多项式 p 除以 q 的商，r 是余式。 [k,r]=deconv(p,q) p=conv(q,k)+r == 多项式的求导: polyder k=polyder(p) : 多项式 p 的导数； k=polyder(p,q): p*q 的导数； [k,d]=polyder(p,q): 有理分式p/q 的导数，k 是分子，d 是分母 k1=polyder([2,-1,0,3]); k2=polyder([2,-1,0,3],[2,1]); [k2,d]=polyder([2,-1,0,3],[2,1]); 例：已知 ， ， 求 [r,p,k] = residue(a,b) a,b为分子，分母多项式系数向量，r,p,k分别是留数，极点，直项 例：将 进行部分分式展开 a=[3 2 5 4 6];b=[1 3 4 2 7 2]; [r,p,k]=residue(a,b) [b,a] = residue(r,p,k)将部分分式展开的形式还原为两个多项式b(x)和a(x)相除的形式 计算多项式在给定点的值 代数多项式求值 y = polyval(p,x): 计算多项式 p 在 x 点的值 注：若 x 是向量或矩阵，则采用数组运算 (点运算)！ p=[2,-1,0,3]; x=2; y=polyval(p,x) x=[-1, 2;-2,1]; y=polyval(p,x) 例：已知 ，分别取 x=2 和一个 2?2 矩阵， 求 p(x) 在 x 处的值 矩阵多项式求值 Y=polyvalm(p,X) 采用的是普通矩阵运算 X 必须是方阵 例：已知 ，则 polyvalm(p,A) = 2*A*A*A