第04章 MATLAB数值运算基础.pptVIP

第 页，共48页 多项式的加减运算 多项式积分 第四章 习题 第4章 MATLAB数值运算基础 在Matlab中用它的系数向量来表示 ： 例： a(x)=3x3+x2+2x+1 a＝[3,1,2,1] 多项式 一、多项式的表示 例： a(x)=3x3+2x+1 a＝[3,0,2,1] 若x某幂次的系数为零，则这个零必须列入系数向量中 例 a=[3,0,2,1]; PAX=poly2sym(a, x) PAX = 3*x^3+2*x+1 格式： poly2sym(a, x) 功能： 系数数组a转换为符号多项式 。 多项式的生成 格式： poly(A) 功能： 生成以A的分量为根，首项系数为1的多项式。 例：已知向量A=[1 –34 –80]，用此向量构造一多项式。 A=[1 -34 -80]; PA=poly(A) PA = 1 113 2606 -2720 (x-1)(x+34)(x+80) =x3+113x2+2606x-2720 poly2sym(PA, x) ans = x^3+113*x^2+2606*x-2720 a=[5 4 3 2 1];b=[3 0 1]; c=a+[0 0 b] 例：求多项式a(x)=5x4+4x3+3x2+2x+1和b(x)=3x2+1的和。 如果阶次不同，低阶的多项式必须用零添补至高阶多项式的阶次。 c = 5 4 6 2 2 多项式的运算 例：求多项式a(x)=5x2+2x+1和b(x)=3x2+1的和。 a=[5 2 1];b=[3 0 1]; c=a+b c = 8 2 2 8x2+2x+2 格式： conv(a,b) 功能： 计算多项式a(x)和b(x)的积。 多项式乘法——conv函数 例：求多项式a(x)=5x4+4x3+3x2+2x+1和b(x)=3x2+1的积。 a=[5 4 3 2 1];b=[3 0 1]; c=conv(a,b) c = 15 12 14 10 6 2 1 15x6+12x5+14x4+10x3+6x2+2x+1 格式： [div,rest]=deconv(a,b) 功能： 计算多项式a(x)和b(x)的商和余式。 多项式除法——deconv函数 例：求多项式a(x)= 15x6+12x5+14x4+10x3+6x2+2x+1和 b(x)=3x2+1的商及余式。 a=[15,12,14,10,6,2,1];b=[3,0,1]; c=deconv(a,b) 商 c = 5.0000 4.0000 3.0000 2.0000 1.0000 [c,d]=deconv(a,b) c = 5.0000 4.0000 3.0000 2.0000 1.0000 d = 0 0 0 0 0 0 余式 15x6+12x5+14x4+10x3+6x2+2x+1 =(5x4+4x3+3x2+2x+1)(3x2+1) 格式： polyder(f) 功能： 求以向量f的分量为系数的多项式的微分。 多项式微分——polyder函数 例：求多项式f(x)= 15x6+12x5+14x4+10x3+6x2+2x+1的微分。 f=[15,12,14,10,6,2,1]; c=polyder(f) c = 90 60 56 30 12 2 poly2sym(c,x) ans = 90*x^5+60*x^4+56*x^3+30*x^2+12*x+2 例 求多项式f(x)=90x5+60x4+56x3+30x2+12x+2的不定积分。 f =[90 60 56 30 12 2]； [f./(length(f):-1:1),0] ans = 15 12 14 10 6 2 0 算法： 设多项式f(x)的系数向量为f，则多项式f(x)的一个积分是[f./(length(f):-1:1),k]，其中k为常数。 格式： roots(f) 功能： 求