资源描述:
《第四章字符串数组、元胞数组和构架数组》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、3.8多项式的表达方式及其操作1、poly:Polynomialwithspecifiedroots0例如p=poly(A);r二roots(p),p二poly(r)Poly:Characteristicpolynomialofmatrixo例如1:p二poly(A)p=poly(A,v);要注意A是否是symbolic・2、conv:Convolutionandpolynomialmultiplicationo例女II:W二conv(u,v),只能两两相乘。但可以相互套用。deconv:Deconvolutionandpolynomialdivisi
2、on。[q,r]=dcconv(v,u)o其中q为商多项式,r为余多项式。3>pascal:Pascalmatrixo4、polyval:Polynomialevaluation。例女口:y=polyval(p,x)returnsthevalueofapolynomialofdegreenevaluatedatx.polyvalm:Matrixpolynomialevaluation。例女13:Y二polyvalm(p,X)evaluatesapolynomialinamatrixsense。5、residue:Convertbetweenpartia
3、lfractionexpansionand。polynomialcoefficientSo即两多项式相除而没有商,故把它们进行部分分式展开。3.9标准数组牛成函数和数组操作函数Kmatlab的随机数算法是通过一个“种子”递推出来的,种子在程序初始的时候都一样。RAND('state',0)resetsthegeneratortoitsinitialstate.重新恢复初始状态,刚运行matlab时候的状态RAND('state',J),forintegerJ,resetsthegeneratortoitsJ-thstate・安排种子产生笫j个状态RAX
4、D(,state,,sum(100*clock))resetsittoadifferentstateeachtime.根据时间产生不同的随机数matlab的rand函数生的是伪随机数,即由种了递推出來的,相同的种了,生成相同的随机数.matlab刚运行起来时,种子都为初始值,因此毎次第一次执行rand得到的随机数都是相同的.remd('state',sum(lOO^clock));A=rand(5,5);rand('state',sum(100*clock));B=rand(5,5);A和B是相同•所以建议再增加一个随机变量,变成:rand('stat
5、e',sum(100*clock)*rand(1));2、diag:Diagonalmatricesanddiagonalsofmatrix0例如:diag(D)%取。阵的对角元。al=diag(A,l)*取人阵“第一上对角线”的元素。diag(diag(X))isadiagonalmatrixo即此命令能够取矩阵对角线上的元素,也可以冇一组向量生成对角线矩阵。3、repmat(D,1,3)%在水平方向“铺放”三个D阵4、reshape:Reshapearray□例如:B二reshape(A,m,n)o即重新排列矩阵A,但必须保证矩阵A得个数等于5、f
6、lipud:矩阵上下对称交换Fliplr:矩阵左右对称交换Rot90:矩阵逆吋针旋转90读6、C二reshape(1:4,2,2)%利用重组操作产生(2x2)矩阵,结果事:c=13247、kron:Kroneckertensorproducto演示Kronecker乘法不具备"口J交换规律”。即kron(A,B)与kron(B,A)是不同的。3.10数组构作技法综合1、数组的赋值扩展法,没有被赋值而扩展的部分为0。2、多次寻访扩展法,例如:AA二A(:,[1:6,1:6])%相当于指令repmat(A,1,2)3、合成扩展法,即将两个矩阵进行行数或列数
7、扩展。如:r二[A;B]or[A,B]。4、triu:Uppertriangularpartofmatrixo例如I:U=triu(X)returnstheuppertriangularpartofX.U=triu(X,k)(注意:不包描主对角线元素)Tril:Lowertriangularpartofmatrixtriu(A,l)+tril(A,-1)%利用操作函数,使主对角元索为全05、[A(l:2,end:-l:1);B(1,:)]%注意灵活合成以及end的应用。6、“对列(或行)同加一个数”三种的操作方法:如:b二[1,2,3];A二resha
8、pe(1:9,3,3);C1二A-b([l,1,1],:);orCl=A-repmat(b,3