第3章MATLAB数值运算.pptVIP

主讲：刘忠伟 3.1 多项式 matlab语言把多项式表达成一个行向量，该向量中的元素是按多项式降幂排列的。 如： f(x)=anxn+an-1xn-1+……+a0 可用行向量 p=[an an-1 …… a1 a0]表示 3.1 多项式 多项式的四则运算包括多项式的加减运算及乘法、除法运算。 3.1 多项式 例:a(x)=x2+2x+3; b(x)=4x2+5x+6; c (x)= = (x2+2x+3)(4x2+5x+6) 》》a=[1 2 3];b=[4 5 6]; 》》c=conv(a,b) c = 4.00 13.00 28.00 27.00 18.00 c(x)= 4 x4 + 13 x3 + 28 x2 + 27 x + 18 3.1 多项式 函数[Q,r]=deconv(P1,P2)用于对多项式P1和P2作除法运算。其中Q返回多项式P1除以P2的商式，r返回P1除以P2的余式。这里，Q和r仍是多项式系数向量。 deconv是conv的逆函数，即有P1=conv(P2,Q)+r。 3.1 多项式 MATLAB提供了两种求多项式值的函数：polyval与polyvalm，它们的输入参数均为多项式系数向量P和自变量x。两者的区别在于前者是代数多项式求值，而后者是矩阵多项式求值。 3.1 多项式 （1）代数多项式求值 polyval函数用来求代数多项式的值，其调用格式为： Y=polyval(P,x) 若x为一数值，则求多项式在该点的值；若x为向量或矩阵，则对向量或矩阵中的每个元素求其多项式的值。 3.1 多项式 （2）矩阵多项式求值 polyvalm函数用来求矩阵多项式的值，其调用格式与polyval相同，但含义不同。polyvalm函数要求x为方阵，它以方阵为自变量求多项式的值。 3.1 多项式 n次多项式具有n个根，当然这些根可能是实根，也可能含有若干对共轭复根。MATLAB提供的roots函数用于求多项式的全部根，其调用格式为： x=roots(P) 其中P为多项式的系数向量，求得的根赋给向量x，即x(1),x(2),…,x(n)分别代表多项式的n个根。 3.1 多项式 例： 求多项式x4+8x3-10的根。 命令如下： A=[1,8,0,0,-10]; x=roots(A) 若已知多项式的全部根，则可以用poly函数建立起该多项式，其调用格式为： P=poly(x) 若x为具有n个元素的向量，则poly(x)建立以x为其根的多项式，且将该多项式的系数赋给向量P。 3.1 多项式 例： 构造多项式x4+3x3-15x3-2x+9的根。 命令如下： A=[1,3,-15,-2,9]; poly2sym(A) 3.2 插值和拟合 一、插值 插值的定义——是对某些集合给定的数据点之间函数的估值方法。 当不能很快地求出所需中间点的函数时，插值是一个非常有价值的工具。 Matlab提供了一维、二维、 三次样条等许多插值选择。 多项式插值函数：interp1 3.2 插值和拟合 二、拟合 拟合函数：polyfit 调用方法： P=polyfit(x,y,n) [p,s]=polyfit(x,y,n) 说明：x,y为已知的数据组，n为要拟合的多项式的阶次，向量p为返回的要拟合的多项式的系数，向量s为调用函数polyval获得的错误预估计值。 3.5 稀疏矩阵 矩阵的存储方式： 1.全元素(Full) 存储—— 完全矩阵 2.稀疏(Sparse)存储——稀疏矩阵 稀疏矩阵存在的必要性： 对大多数元素数值为0的矩阵，若采用满阵方式表示，则0元素将占用相当的存储空间。 稀疏矩阵的特点： 只存储“非零元素”值(按列)和“非零元素”的位置 3.5 稀疏矩阵 对大多数元素数值为0的矩阵，创建稀疏矩阵 稀疏矩阵建立指令 sparse 1、 B = sparse (A) 例. A = [2 0 0 0; 0 0 0 1; 0 4 0 0] A = 2 0 0 0 0 0 0 1 0 4 0 0 B = sparse (A) B = (1,1) 2 (3,2) 4 (2,4) 1 3.5 稀疏矩阵 2、 S = sparse (i,j,s,m,n) “三元组”表示法 i — 非零元素所在的行序