[理学]第七章 matlab在数值分析中的应用.pptVIP

多项式运算 （62~68页） 多项式的表示和创建 在Matlab中，任意的多项式都是用一个行向量表示，将多项式的系数按降幂排列存放在行向量中。多项式p(x)=a0xn+a1xn-1+…+an-1x+an对应的行向量为： P=[a0 a1 … an]。 1. 直接创建多项式 例： p(x)=2x4+5x3-8x+3?[2 5 0 -8 3] 2.通过根创建多项式poly 其调用格式如下： P=poly(A) 如A为向量，则P是多项式的系数行向量，该多项式的根A。 注意： 如A为n×n的矩阵，则P是一个长度为n+1的行向量，它是A的特征多项式的系数； 例： A=[6 -8 6;1 0 0;0 2 0]; p=poly(A), roots(p), D=eig(A) 多项式求根 r=roots(p)：求多项式p的根，r是根组成的向量。 例： p=[2 3 3 0 6]; roots(p) 多项式的基本运算 1. 多项式的乘法conv 多项式的乘法本质上是多项式系数向量的卷积运算。卷积的定义为： k c(k)= ? a(i)b(k+1-i) i=1 在Matlab，完成此功能的命令是conv： C= conv(a,b)。 例： 求a=x3+3x2+2x1+1和b=4x3+3x2+9x1+10的乘积。 a=[1 3 2 1]; b=[4 3 9 10]; c=conv(a,b) 2. 多项式的除法dconv 多项式的除法的运算实质就是多项式系数的解卷积运算。向量a对向量c进行解卷积得到的商向量q和余向量r，满足： k c(k)-r(k)= ? a(i)q(k+1-i) i=1 在Matlab，完成此功能的命令是dconv： [q,r]= dconv(c,a) 其中： c=conv(a,q)+r 例: a=[1 3 2 1]; c=[4 15 26 47 51 29 10]; [q,r]=deconv(c,a) 3. 多项式的加减 对多项式所对应的向量利用数组加减运算规则执行加减运算。 注：多项式所对应的向量如果大小不相等，必须用0补齐。 例： 求a=x3+3x2+2x1+1和b=4x3+3x2+9x1+10的和与差。 a=[1 3 2 1]; b=[4 3 9 10]; c=a+b, c1=a-b, 4. 多项式的求值polyval 在Matlab中，多项式求值有两种方式：按数组规则和按矩阵规则计算多项式值。 y=polyval(p,x)：按数组规则运算。用来计算多项式p在自变量x处的值，可以是数、向量和矩阵。如果x是向量或矩阵，则该命令将对x的每个元素计算p的值，结果大小和x的大小相同 例： p=[2 3 3 0 6]; x=[1 -1]; y=polyval(p,x) 5. 多项式的部分分式展开 没有多重根： B(s) R(1) R(2) R(n) ---- = -------- + -------- + ... + -------- + K(s) A(s) s - P(1) s - P(2) s - P(n) 有多重根时，则包含有： R(j) R(j+1) R(j+m-1) -------- + ------------ + ... + ------------ s - P(j) (s - P(j))^2 (s - P(j))^m 在Matlab中完成此任务的指令