第4章Matlab其他函数库.pptVIP

第4章 Matlab的其他函数库;第1节 数据分析数据库4.1.1 基本的数据分析;4.1.2 用于场论的数据分析;4.1.3 用于随机数据分析的函数;4.1.4 用于相关分析和傅立叶分析的函数;第2节 矩阵的分解与变换（matfun函数库） 4.2.1 线性方程组的系数矩阵;4.2.2 矩阵的分解p68;常用的矩阵函数;第3节 多项式函数库（polyfun）;4.3.1 多项式简介;例：已知向量A=[1 –34 –80]，用此向量构造一多项式并显示结果。 A=[1 -34 -80]; PA=poly(A) PA = 1 113 2606 -2720 PAX=poly2sym(PA,X) PAX = X^3+113*X^2+2606*X-2720;4.3.2 多项式运算;多项式乘法 多项式乘法采用conv函数。 据上例求积： d=conv(a,b) d = 15 12 14 10 6 2 1 多项式除法 由deconv函数完成，结果包括商和余数两部分。 [div,rest]=deconv(d,a) div = 3 0 1 rest = 0 0 0 0 0 0 0;微分 提供了微分函数　polyder 例：求多项式p(x)=2x4-6x3+3x2+0x+7的微分。 p=[2 -6 3 0 7]; q=polyder(p) q = 8 -18 6 0 pd=poly2sym(q,x) pd = 8*x^3-18*x^2+6*x ;求根 多项式求根可用roots函数。 例：求多项式p(x)=2x4-6x3+3x2+0x+7的根。 p=[2 -6 3 0 7]; x=roots(p) x = 1.9322 + 0.4714i 1.9322 - 0.4714i -0.4322 + 0.8355i -0.4322 - 0.8355i;求值 采用polyval函数可以求出当多项式中的未知数为某个特定值时该多项式的值。 例：求上例中x=1时的值。 p=[2 -6 3 0 7]; polyval(p,1) ans = 6 ;常用的多项式函数;应用举例;例2：已知一元四次方程所对应的四个根为 -5.0000 4.0000 3.0000 3.0000 求这个方程所对应的表达式原型。 x=[-5 4 3 3]; p=poly(x); px=poly2sym(p,x) px = x^4-5*x^3-17*x^2+129*x-180 ;习题;习题6：求解一元六次方程 3x6+12x5+4x4+7x3+8x+1=0 的根。 习题7：求多项式 3x6+12x5+4x4+7x3+8x+1=0 被(x-3)(x3+5x)除后的结果。;4.3.3 多项式拟合;x=0:0.1:1;y=[-0.447,1.978,3.28,6.16,7.08,7.34,7.66,9.56,9.48,9.30,11.2]; a1=ployfit(x,y,1); xi=linspace(0,1); yi1=polyval(a1,xi);plot(x,y,o,xi,yi1,b),hold on a2=polyfit(x,y,2);yi2=polyval(a2,xi);plot(x,y,o,xi,yi2,m);hold on a3=polyfit(x,y,3);yi3=polyval(a3,xi);plot(x,y,o,xi,yi3,r);hold on a9=polyfit(x,y,9);yi9=polyval(a9,xi);plot(x,y,o,xi,yi9,c);hold on a10=polyfit(x,y,10);yi10=polyval(a10,xi);plot(x,y,o,xi,yi10,g);hold on;4.3.4 多项式插值 ;s p l i n e ( x , y , x x ) 等同于interp1(x, y, xx, ’s p l i n e’)，但是参数必须是向量。 s p l i n e ( x , y ) 返回三次样条插值向量的 p p形式，它是函数y=f ( x )的近似值。 p p是’piecewise polynomial’的缩写，得到的向量元素包含 计算的三次样条系数。这个命令可以被p p v a l函数使用。;4.3.5线性微分方程的解;第4节 函数功能和数值积分函数库（funfun）;定积分子程序 quad（’函数名’,初值x0,终值xf） S=quad(‘humps’,1,2