四MATLAB数值计算功能9.docVIP

第四章 MATLAB 的数值计算功能 Chapter 4: Numerical computation of MATLAB 一、多项式(Polynomial)` 1 二. 线性代数(Linear Algebra) 8 三． 逆矩阵及行列式(Revers and determinant of matrix) 12 四． 矩阵分解(Matrix decomposition) 15 五． 数据分析(Data Analyaia) 20 六． 插值: (Interpolation) 27 七． 数字信号处理初步(Introduction to Signal Processing ) 30 一、多项式(Polynomial)` 1．多项式的表达与创建(Expression and Creating of polynomial) (1) 多项式的表达（expression of polynomial）_ Matlab用行矢量表达多项式系数(Coefficient)，各元素按变量的降幂顺序排列，如多项式为： P(x)=a0xn+a1xn-1+a2xn-2…an-1x+an 则其系数矢量(Vector of coefficient)为：P=[a0 a1 … an-1 an] 如将根矢量(Vector of root)表示为： ar=[ ar1 ar2 … arn] 则根矢量与系数矢量之间关系为： （x-ar1）(x- ar2) … (x- arn)= a0xn+a1xn-1+a2xn-2…an-1x+an (2)多项式的创建（polynomial creating） a)系数矢量的直接输入法 利用poly2sym函数直接输入多项式的系数矢量，就可方便的建立符号形式的多项式。 例：创建多项式x3-4x2+3x+2 poly2sym([1 -4 3 2]) ans = x^3-4*x^2+3*x+2 POLY Convert roots to polynomial. POLY(A), when A is an N by N matrix, is a row vector with N+1 elements which are the coefficients of the characteristic polynomial, DET(lambda*EYE(SIZE(A)) - A) . POLY(V), when V is a vector, is a vector whose elements are the coefficients of the polynomial whose roots are the elements of V . For vectors, ROOTS and POLY are inverse functions of each other, up to ordering, scaling, and roundoff error. b) 由根矢量创建多项式 通过调用函数 p=poly(ar)产生多项式的系数矢量, 再利用poly2sym函数就可方便的建立符号形式的多项式。 注：（1）根矢量元素为n ，则多项式系数矢量元素为n+1； （2）函数poly2sym(pa) 把多项式系数矢量表达成符号形式的多项式，缺省情况下自变量符号为x，可以指定自变量。 （3）使用简单绘图函数ezplot可以直接绘制符号形式多项式的曲线。 例 1：由根矢量创建多项式。将多项式(x-6)(x-3)(x-8)表示为系数形式 a=[6 3 8] %根矢量 pa=poly(a) %求系数矢量 ppa=poly2sym(pa) %以符号形式表示原多项式 ezplot(ppa,[-50,50]) pa = 1 -17 90 -144 ppa = x^3-17*x^2+90*x-144 注：含复数根的根矢量所创建的多项式要注意： （1）要形成实系数多项式，根矢量中的复数根必须共轭成对； （2）含复数根的根矢量所创建的多项式系数矢量中，可能带有很小的虚部，此时可采用取实部的命令(real)把虚部滤掉。 进行多项式的求根运算时，有两种方法，一是直接调用求根函数roots，poly (POLY Convert roots to polynomial)和 roots 互为逆函数。另一种是先把多项式转化为伴随矩阵，然后再求其特征值，该特征值即是多项式的根。 POLY2SYM Polynomial coefficient vector to symbolic polynomi