第五讲MATLAB在高等代数中的应用之二.ppt

* * 第五讲 MATLAB在高等代数中的应用之二 一、多项式的建立与表示 任意多项式都可以用一个行向量来表示，即n+1维向量a表示多项式y(x)=a(1)*x^n+a(2)*x^(n-1)++a(n)*x+a(n+1),可经看出n+1维向量表示阶数为n的多项式，且任意一个向量就可以作为多项式。 p=[1 -6 11 -6] poly2sym(p,x) Poly(A)A是一个n×n的矩阵，此函数返回矩阵Ａ特征多项式p，p是n+1维向量，特征多项式的根就是矩阵Ａ的特征值 高等代数Ｐ３０７习题１９－３ A=[1 1 1 1;1 1 -1 -1;1 -1 1 -1;1 -1 -1 1] ｐ＝poly(A) roots(p)　　　　%求多项式的根 eig(A)　　　　　％求Ａ的特征值 1、poly函数 P=poly(r）r为向量此函数返回以向量中的元素为根的多项式 高等代数p307 19-1 r=[7 -2],p=poly(r) polyvalm(p,A)　％求Ａ关于多项式的值 ２、多项式的符号表示 Poly2sym(p)函数将向量表示为多项式表示为符号表示的多项式默认的字符为‘x’ poly2sym(p,v) 表示以字符‘v’的多项式。 也可以sym v Poly2sym(p,v) 二、多项式的基本运算 １、多项式的加减运算 function f=sumpoly(p1,p2) n=abs(length(p1)-length(p2)); if length(p1)length(p2) p12=[zeros(1,n),p1]; p22=p2; elseif length(p1)length(p2) p22=[zeros(1,n),p2]; p12=p1; Else p12=p1;p22=p2; end f=p12+p22; p1=[1 2 3] p2=[1 1 1 1] sumpoly(p1,p2) 2、多项式的乘除法运算 conv(p1,p2) 此函数求多项式p1与p2的乘积，即求向量p1与p2的卷积。如果m=length(p1),n=length(p2),则w=conv(p1,p2)的长度为m+n-1,并且满足w(k)=a(1)b(k)+a(2)b(k+1)++a(k)b(2k-1) 高等代数p14例 f=[1 3 -1 -4 -3] g=[ 3 10 2 -3] [py,r]=deconv(f,g) 3、多项式的根 roots(x)此函数返回多项式的根组成的向量，同时也是友元阵的特征向量对系数为实数的多项式的根，则如果其根出现复数，则复数必是成对出现的。 p=[2 -1 0 2 -3] roots(p) B=compan(p) eig(B) 4、多项式的数组和矩阵运算 y=polyval(p,x)计算多项式在x处的值，x可以为矩阵或向量，此时函数计算多项式在x的每个元素处的值。 polyval(p,1) polyval(p,0) x=[1 2] polyval(p,x) b=compan(p) polyvalm(p,b) 5、多项式的部分分式展开（留数计算） 函数residue()在部分分式展开和多项式系数之间进行转换。 命令格式及功能 [r,p,k]=resicue(b,a)找出两个多项式b(x)和a(x)之比b(x)/a(x)的留数，极点和直项向量，分别用r,p,k表示。 [b,a]=residue(r,p,k)将部分分式展开的形式还原为两个多项式b(x)和a(x)相除的形式。 如果多项式a(x)不含重根，则两个多项式之比可以写成b(x)/a(x)=r1/(x-p1)+r2/(x-p2)+ … +rn/(x-pn)+k(x),其中p1,p2,p3, …,pn称为极点，r1,r2, …,rn称为留数，k(x)称为直项。留数和极点的个数n满足 n=length(a)-1=length(r)=length(p)。 如果b的次数小于a的次数，则直项的系数向量的系数为空；否则，，它们之间满足length(k)=length(b)-length(a)+1. 如果a(x)含有m重p(j)==p(j+m-1),则展开的这m项应该写成r(j)/r-pj)+r(j+1)/(x-pj)2+…+r(j+m-1)/(s-pj)m a=[2 4 6 8 1] ,b=[1 3 9 2] [r,p,k]=residue(b,a) %将多项式除部分分式展开 [b1,a1]=residue(r,p,k) %将部分分式还原 例：吉米多维奇习题集1869题 a=[1 -5 6 0] b=[1 0 0 1] format rat [r,p,k]=residue(b,a) r =28/3