第四章的应用.pptVIP

第四章的应用第1页，共26页，星期日，2025年，2月5日4.1多项式的运算多项式的生成和表达1、多项式的表达在MATLAB环境下多项式是用向量的形式表达的。向量最右边的元素表示多项式的0阶，向左数依次表示多项式的第1阶、第2阶、第3阶……。例如多项式5x4+3x2+2x+1表示为[50321]。第2页，共26页，星期日，2025年，2月5日2、多项式的生成语法：poly(MA)说明：1、若MA为方阵，则生成的多项式P为方阵MA的特征多项式。2、若MA为向量、则向量和多项式满足这样一种关系：MA＝[r1r2r3…rn]生成的多项式为：(x-r1)(x-r2)(x-r3)…(x-rn)=a0xn+a1xn-1+a2xn-2+…+an即生成的向量为[a0a1a2…an]3、当然也可以用直接输入的方式生成多项式。第3页，共26页，星期日，2025年，2月5日3.多项式的加减乘除加减法多项式的加减法与向量的加减法完全相同。乘除法A.c=conv(a，b)乘积运算[qr]=deconv(b,a)除法运算(即有b=conv(a,q)+r)1、a、b和c,q,r均为以向量形式表示的多项式。2、q表示除运算的商，r表示除运算的余数。第4页，共26页，星期日，2025年，2月5日例：求多项式F(x)=x2+5x和G(x)=2x+1的乘积M(x)。f=[150];g=[21];m=conv(f,g)m=21150即M(x)=2x3+11x2+5x+0第5页，共26页，星期日，2025年，2月5日已知：s(x)=5x4+3x3+2x2+x+1,a(x)=x+1求s(x)被a(x)除的结果s=[53211];a=[11];[q,r]=deconv(s,a);qq=5-24-3rr=00004即有q(x)=5x3-2x2+4x-3r(x)=4第6页，共26页，星期日，2025年，2月5日3.多项式求导dp=polyder(p)多项式p(x)的导数dp(x)dp=polyder(p,q)多项式p(x)*q(x)的导数[b,a]=polyder(p,q)多项式p(x)/q(x)的导数b(x)/a(x)第7页，共26页，星期日，2025年，2月5日例：求多项式:p(x)=x3+2x2+x+1的一阶和二阶导数。p=[1211];dp=polyder(p)dp2=polyder(dp)计算结果dp=341%3x2+4x+1dp2=64%6x+4第8页，共26页，星期日，2025年，2月5日例：求的一阶导数[b,a]=polyder([151],[35])b=31022a=93025polyder([151],[35])ans=94028第9页，共26页，星期日，2025年，2月5日4.多项式的求根r＝roots(p)p为多项式的向量表示式r为多项式的根，以列向量表示第10页，共26页，星期日，2025年，2月5日例：求方程x3-3x2-3x+1=0的解。p=[1-3-31];r=roots(p)r=3.7321-1.00000.2679也即有(x-3.7321)(x-(-1))(x-0.2679)=x3-3x2-3x+1第11页，共26页，星期日，2025年，2月5日注意：poly和roots功能刚好相反p=[12345];%x4+2x3+3x2+4x+5r=roots(p)ans=0.2878+1.4161i0.2878-1.4161i-1.2878+0.8579i-1.2878-0.8579ipoly(r)ans=1.00002.00003.0000