matlab第六章课件.ppt

第六章 数值计算 多项式运算 2. 多项式函数 2. 插值运算 2．求矩阵的最大值和最小值 求矩阵A的最大值的函数有3种调用格式，分别是： (1) max(A)：返回一个行向量，向量的第i个元素是矩阵A的第i列上的最大值。 * * 多项式计算 插值运算 数据分析 1. 表示方法 ★ MATLAB采用行向量表示多项式系数，多项式系数按降幂排列。 例如： f(x)=anxn+an-1xn-1+……+a0 可用行向量 p=[an an-1 …… a1 a0]表示 ★ 函数poly2str()将多项式系数向量转换为完整形式。 例如： a=[2 3 4 5 0]; poly2str(a,’x’) 多项式求值 y = polyval(p,x)，p为多项式系数行向量，x代入多项式的值（对每一个元素计算多项式的值）； Y = polyvalm(p,X)，把矩阵X代入多项式p中进行计算。 多项式乘法和除法 w = conv(u,v)，实现多项式乘法，返回结果多项式的系数行向量； [q,r] = deconv(u, v)，实现多项式除法。q和r仍是多项式系数向量。满足u=conv(q,v)+r 例如：计算多项式x4+8x3-10与多项式2x2-x+3的乘积和除法，并显示完整结果 a=[1 8 0 0 -10]; b=[2 -1 3]; c=conv(a,b); [q,r]=deconv(a,b); sc=poly2str(c,’x’) sq=poly2str(q,’x’) sr=poly2str(r,’x’) 对于次数相同的若干个多项式，可直接对多项式系数向量进行加、减的运算。如果多项式的次数不同，则应该把低次的多项式系数不足的高次项用零补足，使同式中的各多项式具有相同的次数。 多项式加法和减法 例如：计算多项式x4+8x3-10与多项式2x2-x+3的加法和减法 a=[1 8 0 0 -10]; b=[0 0 2 -1 3]; p1=a-b; p2=a+b; sp1=poly2str(p1,’x’) sp2=poly2str(p2,’x’) 多项式的微积分 k = polyder(p)，返回多项式p微分的系数向量； k = polyder(a,b)，返回多项式a*b乘积微分的系数向量； [q,d] = polyder(b,a)，返回多项式b/a微分的系数向量。 s=polyint(p,k)，返回多项式p不定积分的系数向量。（k为积分常数项，缺省为0） 例如：用函数polyint()求多项式x2+3x+4的不定积分 a=[1 3 4]; S1=polyint(a) s2=polyint(a,3) 多项式的根和由根创建多项式 r = roots(c)，返回多项式c的所有根r，r是向量，其长度等于根的个数。 p = poly(r)，输入r是多项式所有根，返回值为多项式的系数向量 例如：已知f(x)=x4+8x3-10，计算f(x)=0的全部根；利用这些根构造一个多项式g(x),并与f(x)进行对比 p=[1 8 0 0 -10]; x=roots(p) %求方程f(x)=0的根 g=poly(x) %求多项式g(x) 多项式部分分式展开 [r,p,k] = residue(b,a)，求多项式b/a的部分分式展开； r是分子部分，p是分母部分，k是 余式 [b,a] = residue(r,p,k)，从部分分式得到多项式向量 例如：求多项式的部分分式 p1=[1 2 0 0 0 0 1]; p2=[1 1 -8 -6 8 24]; [r,p,k]=residue(p1,p2) 多项式曲线拟合 曲线拟合的目的用一个较简单的函数去逼近一个复杂的或未知的函数；matlab曲线拟合的最优标准是最小二乘拟合法。 设测得n个离散数据点(xi,yi),今欲构造一个m(m≤n)次多项式p(x): 原理:使拟合多项式在各节点出的偏差p(xi)-yi的平方和达到最小。 p = polyfit(x,y,n)，采用n次多项式p来拟合数据x和y。 在工程测量和科学实验中，所得到的数据通常都是离散的。如果要得到这些离散点以外的其他点的数值，就需要根据这些已知数据进行估算，即插值。 例如，测量得n个点的数据为(x1,y1), (x2,y2),…, (