第3章_MATLAB数据分析与多项式计算.pptVIP

第3章 MATLAB多项式计算与数据分析 3.1 多项式计算 3.2 线性代数 3.3 数据统计处理 3.4 数据插值 3.5 曲线拟合 3.6 离散傅立叶变换 3.1 多项式计算 3.1.1 多项式的四则运算 1．多项式的加减运算 2．多项式乘法运算 函数conv(P1,P2)用于求多项式P1和P2的乘积。这里，P1、P2是两个多项式系数向量。 例3-1 求多项式x4+8x3-10与多项式2x2-x+3的乘积。 3．多项式除法 函数[Q,r]=deconv(P1,P2)用于对多项式P1和P2作除法运算。其中Q返回多项式P1除以P2的商式，r返回P1除以P2的余式。这里，Q和r仍是多项式系数向量。 deconv是conv的逆函数，即有P1=conv(P2,Q)+r。 例3-2 求多项式x4+8x3-10除以多项式2x2-x+3的结果。 3.1.2 多项式的导函数 对多项式求导数的函数是： p=polyder(P)：求多项式P的导函数 p=polyder(P,Q)：求P·Q的导函数 [p,q]=polyder(P,Q)：求P/Q的导函数，导函数的分子存入p，分母存入q。 上述函数中，参数P,Q是多项式的向量表示，结果p,q也是多项式的向量表示。 例3-3 求有理分式的导数。 命令如下： P=[1]; Q=[1,0,5]; [p,q]=polyder(P,Q) 3.1.3 多项式的求值 MATLAB提供了两种求多项式值的函数：polyval与polyvalm，它们的输入参数均为多项式系数向量P和自变量x。两者的区别在于前者是代数多项式求值，而后者是矩阵多项式求值。 1．代数多项式求值 polyval函数用来求代数多项式的值，其调用格式为： Y=polyval(P,x) 若x为一数值，则求多项式在该点的值；若x为向量或矩阵，则对向量或矩阵中的每个元素求其多项式的值。 例3-4 已知多项式x4+8x3-10，分别取x=1.2和一个2×3矩阵为自变量计算该多项式的值。 2．矩阵多项式求值 polyvalm函数用来求矩阵多项式的值，其调用格式与polyval相同，但含义不同。polyvalm函数要求x为方阵，它以方阵为自变量求多项式的值。设A为方阵，P代表多项式x3-5x2+8，那么polyvalm(P,A)的含义是： A*A*A-5*A*A+8*eye(size(A)) 而polyval(P,A)的含义是： A.*A.*A-5*A.*A+8*ones(size(A)) 例3-5 仍以多项式x4+8x3-10为例，取一个2×2矩阵为自变量分别用polyval和polyvalm计算该多项式的值。 3.1.4 多项式求根 n次多项式具有n个根，当然这些根可能是实根，也可能含有若干对共轭复根。MATLAB提供的roots函数用于求多项式的全部根，其调用格式为： x=roots(P) 其中P为多项式的系数向量，求得的根赋给向量x，即x(1),x(2),…,x(n)分别代表多项式的n个根。 例3-6 求多项式x4+8x3-10的根。 命令如下： A=[1,8,0,0,-10]; x=roots(A) 若已知多项式的全部根，则可以用poly函数建立起该多项式，其调用格式为： P=poly(x) 若x为具有n个元素的向量，则poly(x)建立以x为其根的多项式，且将该多项式的系数赋给向量P。 例3-7 已知 f(x) (1) 计算f(x)=0 的全部根。 (2) 由方程f(x)=0的根构造一个多项式g(x)，并与f(x)进行对比。 命令如下： P=[3,0,4,-5,-7.2,5]; X=roots(P) %求方程f(x)=0的根 G=poly(X) %求多项式g(x) 3.2代数方程组求解 matlab中有两种除运算左除和右除。 对于方程ax+b，a 为an×m矩阵，有三种情 况： ? 当n=m时，此方程成为“恰定”方程 ? 当nm时，此方程成为“超定”方程 ? 当nm时，此方程成为“欠定”方程 matlab定义的除运算可以很方便地解上 述三种方程 3.2.1恰定方程组的解 方程ax+b(a为非奇异) x=a-1 b 矩阵逆 两种解: x=inv(a)?b — 采用求逆运算解方程 x=a\b — 采用左除运算解方程 方程ax=b a=[1 2;2 3];b=[8;13]; ?x=inv(a)*b ? x=a\b x = x = 2.00