第十三讲 多项式的计算.ppt

第十三讲 多项式运算 作 业 1．已知某班的5名学生的三门课成绩列表如下： 学生序号 1 2 3 4 5 高等数学 78 89 64 73 68 外语 83 77 80 78 70 MATLAB语言 82 91 78 82 68 试写出有关命令，先分别找出三门课的最高分及其学生序号；然后找出三门课总分的最高分及其学生序号。 2．针对上小题的成绩表，求出其三门课总分存入数组ZF，再利用sort命令对之按升序排序，同时把相应的学生序号存入数组XH。 chengji=[78 89 64 73 68;...] 83 77 80 78 70;... 82 91 78 82 68] [singlesc,xuhao]=max(chengji) [zf]=sum(chengji) [maxzf,nmb1]=max(zf) [zfsort,xh]=sort(zf) [zfsort,xh]=sort(zf,’descend’) 4.3 多项式运算 多项式的一般形式表达为： poly2sym(A) %给出A多项式的表达式。 例：A=[1 2 3] f=poly2sym(A) f = x^2+2*x+3 3、poly(x) %产生以x的元素为根的多项式 例：x=[1,-1] A=poly(x) f=poly2sym(A) 二. 多项式求根 命令格式：x=roots(A)。这里A为多项式的系数A(1),A(2),…,A(N),A(N+1)；解得的根赋值给数组X，即X(1),X(2), …,X(N)。 【例】试用roots函数求多项式x4+8x3-10=0的根 这是一个4次多项式，它的五个系数依次为：1,8,0,0,-10。下面先产生多项式系数的向量A，然后求根： A=[1 8 0 0 -10] A = 1 8 0 0 -10 x=roots(A) 三. 求多项式的值 polyval函数用来求代数多项式的值，调用的命令格式为： Y=polyval(A,x) 本命令将polyval函数返回的多项式的值赋值给Y。若x为一数值，则Y也为一数值；若x为向量或矩阵，则对向量或矩阵中的每个元素求其多项式的值。 【例】求多项式x4+8x3-10=0的值，分别取x=1和下面的矩阵的2×3个元素为自变量，计算该多项式的值。 A=[1 8 0 0 -10]; % 例7.8的4次多项式系数 x=1; % 取自变量为一数值 y1=polyval(A,x) y1 = -1 x=[-1 0 1;2 3 4] % 给出一个矩阵x y1=polyval(A,x) y1 = -17 -10 -1 70 287 758 四、多项式的四则运算 (1)多项式加、减 对于次数相同的若干个多项式，可直接对多项式系数向量进行加、减的运算。如果多项式的次数不同，则应该把低次的多项式系数不足的高次项用零补足，使同式中的各多项式具有相同的次数。 a=[2,4,6,8] b=[3,6,9] d=a+[0,b] %补零 (2)多项式乘法 若A、B是由多项式系数组成的向量，则CONV函数将返回这两个多项式的乘积。调用它的命令格式为： C=conv(A,B) 命令的结果C为一个向量，由它构成一个多项式。 【例】求多项式x4+8x3-10与多项式2x2-x+3的乘积。 A=[1 8 0 0 -10]; B=[2 -1 3] B = 2 -1 3 C=conv(A,B) C = 2 15 -5 24 -20 10 -30 本例的运行结果是求得一个6次多项式 2x6+15x5-5x4+24x3-20x2+10x-30 (3)多项式除法 当A、B是由多项式系数组成的向量时，DECONV函数用来对两个多项式作除法运算。调用的命令格式为： [Q,r]=deconv(A,B) 本命令的结果：多项式A除以多项式B获商多项式赋予Q（也为多项式系数向量）；获余项多项式赋予r（其系数向量的长度与被除多项式相同，通常高次项的系数为0）。 DECONV是CONV的逆函数，即有A=conv(B,Q)+r。 【例】试用多项式x4+8x3-10与多项式2x2-