[信息与通信]第5章 MATLAB 的数学运算.pptVIP

第5章 MATLAB 的数学运算 多项式与插值 多项式在数学中有着极为重要的作用，同时多项式的运算也是工程和应用中经常遇到的问题。MATLAB 提供了一些专门用于处理多项式的函数，用户可以应用这些函数对多项式进行操作。MATLAB 中对多项式的操作包括多项式求根、多项式的四则运算及多项式的微积分。 多项式的表示 在 MATLAB 中多项式用一个行向量表示，向量中的元素为该多项式的系数，按照降序排列。如多项式 可以表示为向量 p=[9 7 4 3]。用户可以创建向量的方式创建多项式，再将其显示为多项式， 多项式的四则运算 由于多项式是利用向量来表示，多项式的四则运算可以转化为向量的运算。 多项式的加减为对应项系数的加减，因此可以通过向量的加减来实现。但是在向量的加减中两个向量需要有相同的长度，因此在进行多项式加减时，需要将短的向量前面补 0。 多项式的乘法实际上是多项式系数向量之间的卷积运算，可以通过 MATLAB 中的卷积函数 conv 来完成。 多项式的除法为乘法的逆运算，可以通过反卷积函数 deconv 来实现。 多项式的其他运算 除多项式的四则运算外，MATLAB 还提供了多项式的一些其他运算。这些运算及其对应的函数如表所示。 多项式的运算函数（1/2） roots 函数和 poly 函数 这两个函数为功能互逆的两个函数。roots 函数用于求解多项式的根。该函数的输入参数为多项式的系数组成的行向量，返回值为由多项式的根组成的列向量。poly 函数用于生成根为制定数值的多项式。 polyval 函数 polyval 函数用于多项式求值。对于给定的多项式，利用该函数可以计算该多项式在任意点的值。 多项式的运算函数（2/2） polyder 函数 函数 polyder 用于多项式求导。该函数可以用于求解一个多项式的导数、两个多项式乘积的导数和两个多项式商的导数。该函数的用法为： q = polyder(p) 该命令计算多项式 p 的导数。 c = polyder(a,b) 该命令实现多项式 a、b 的积的导数。 [q,d] = polyder(a,b) 该命令实现多项式a、b 的商的导数，q/d 为最后的结果。 多项式运算 MATLAB语言把多项式表达成一个行向量，该向量中的元素是按多项式降幂排列的。 f(x)=anxn+an-1xn-1+……+a0可用行向量 p=[an an-1 …… a1 +a0]表示。 1、poly 产生特征多项式系数向量。 (1)特征多项式一定是n+1维的； (2)特征多项式第一个元素一定是1。 例: a=[1 2 3;4 5 6;7 8 0]; p=poly(a) p = 1.00? -6.00? -72.00? -27.00 p是多项式p(x)=x3-6x2-72x-27的MATLAB描述方法，我们可用： p1=poly2str(p,x) — 函数文件，显示数学多项式的形式p1 =x^3 - 6 x^2 - 72 x - 27 多项式拟合 曲线拟合是工程中经常要用到的技术之一。MATLAB 提供了曲线拟合工具箱满足用户要求，另外，还提供了多项式拟合函数。函数 polyfit 给出在最小二乘意义下最佳拟合系数。该函数的调用格式为： p = polyfit(x,y,n) 其中x、y分别为待拟合数据的 x 坐标和 y 坐标，n 用于指定返回多项式的次数。 数据插值 多项式插值是指根据给定的有限个样本点，产生另外的估计点以达到数据更为平滑的效果。该技巧在信号处理与图像处理上应用广泛。 插值的定义——是对某些集合给定的数据点之间函数的估值方法。 当不能很快地求出所需中间点的函数时，插值是一个非常有价值的工具。 MATLAB提供了一维、二维、 三次样条等许多插值选择。 所用指令有一维的interp1、二维的interp2、三维的interp3。这些指令分别有不同的方法（method），设计者可以根据需要选择适当的方法，以满足系统属性的要求。Help polyfun可以得到更详细的内容。 一维插值 一维插值在曲线拟合和数据分析中具有重要的地位。在 MATLAB 中，一维插值由函数 interp1 实现。该函数的调用格式为 yi = interp1(x,y,xi,method) x、y：采用数据的 x 坐标和 y 坐标 xi ：待插值的位置 method：采用的插值方法 该语句返回函数在点 xi 处的插值结果。该语句中的参数 method 可以选择的内容如表所示。 数据插值 在MATLAB中，实现一维插值的函数是interp1，其调用格式为： y=interp1(xs,ys,x,’method’) 在有限样本点向量xs与ys中，插