第四章+MATLAB的数值计算功能.docVIP

第四章 MATLAB 的数值计算功能 Chapter 4: Numerical computation of MATLAB 数值计算是MATLAB最基本、最重要的功能，是MATLAB最具代表性的特点。MATLAB在数值计算过程中以数组和矩阵为基础。数组是MATLAB运算中的重要数据组织形式。前面章节对数组、矩阵的特征及其创建与基本运算规则等相关知识已作了较详尽的介绍，本章重点介绍常用的数值计算方法。 一、多项式(Polynomial)` 多项式在众多学科的计算中具有重要的作用，许多方程和定理都是多项式的形式。MATLAB提供了标准多项式运算的函数，如多项式的求根、求值和微分，还提供了一些用于更高级运算的函数，如曲线拟合和多项式展开等。 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函数直接输入多项式的系数矢量，就可方便的建立符号形式的多项式。 例1：创建多项式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) 由根矢量创建多项式 已知根矢量ar, 通过调用函数 p=poly(ar)产生多项式的系数矢量, 再利用poly2sym函数就可方便的建立符号形式的多项式。 例2：由根矢量创建多项式。将多项式(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）根矢量元素为n ，则多项式系数矢量元素为n+1； （2）函数poly2sym(pa) 把多项式系数矢量表达成符号形式的多项式，缺省情况下自变量符号为x，可以指定自变量。 （3）使用简单绘图函数ezplot可以直接绘制符号形式多项式的曲线。 例 3： 由给定复数根矢量求多项式系数矢量。 r=[-0.5 -0.3+0.4i -0.3-0.4i]; p=poly(r) pr=real(p) ppr=poly2sym(pr) p = 1.0000 1.1000 0.5500 0.1250 pr = 1.0000 1.1000 0.5500 0.1250 ppr = x^3+11/10*x^2+11/20*x+1/8 说明：含复数根的根矢量所创建的多项式要注意： （1）