China-pub.com 下载 第10章 函数、插值和曲线拟合分析 MATLAB可以处理有估值和没有估值的多项式，还有一些强大的数值分析命令，如求零 值和最小值。MATLAB中还有数据集合的插值、曲线拟合的命令和函数，还提到了经典的贝 赛尔(Bessel)函数。 10.1 MATLAB中的多项式 MATLAB将阶为n的多项式px)存储在长度为n+I的行向量p中。元素为多项式的系数，并 按x的幂降序排列，表示为： p=(an an-1...al ao) 代表的多项式为： p(x)=anx"an-1x-I+...+aix+ao 令A是一个稀疏矩阵，向量p和q的长度分别为n+1和m+1,它们分别表示次数为n和m的多 项式。MATLAB中处理多项式的命令如下： 命令集99 多项式 polyval(p,x) 计算多项式P。如果x是一个标量，则计算出多项式在x点 的值：如果x是一个向量或者一个矩阵，则计算出多项式 在x中所有元素上的值。 [y,err]= 计算向量x的多项式p的值。同上，计算结果在y中，同时 polyval(p,x,E) 还根据polyfit命令给出的矩阵E返回一个误差估计向量 err。见help polyvaf和nelp polyfit参见10.4节。 polyvalm(p,A) 直接对矩阵A进行多项式计算。不是象上个命令一样对每个 元素进行多项式计算，而是计算p(A)=pA+BA-++pl。 poly(A) 计算矩阵A的特征多项式向量。 poly(x) 给出一个长度为+1的向量，其中的元素是次数为n的多项 式的系数。这个多项式的根是长度为n的向量x中元素。 compan (p) 计算带有系数P的多项式的友矩阵A,这个矩阵的特征多 项式为p。 roots(p) 计算特征多项式即的根，是一个长度为m的向量，也就是方程 0的解。表达式oly(roots(p)p为真。结果可以是复数。 conv(p,q) 计算多项式p和q的乘积，也可以认为是p和q的卷积。 [k,r]=deconv (p,q) 计算多项式p除q。k是商多项式，r是残数多项式。这个 计算等价于p和q的逆卷积。下载 第1 0章 函数、插值和曲线拟合分析 M AT L A B可以处理有估值和没有估值的多项式，还有一些强大的数值分析命令，如求零 值和最小值。M AT L A B中还有数据集合的插值、曲线拟合的命令和函数，还提到了经典的贝 赛尔( B e s s e l )函数。 10.1 MAT L A B中的多项式 M AT L A B将阶为n的多项式p ( x )存储在长度为n+ 1的行向量p中。元素为多项式的系数，并 按x的幂降序排列，表示为： 代表的多项式为： 令A是一个稀疏矩阵，向量 p和q的长度分别为n+ 1和m+ 1，它们分别表示次数为 n和m的多 项式。M AT L A B中处理多项式的命令如下： 命令集9 9 多项式 p o l y v a l ( p , x ) 计算多项式 p。如果x是一个标量，则计算出多项式在 x点 的值；如果 x是一个向量或者一个矩阵，则计算出多项式 在x中所有元素上的值。 [ y , e r r ] = 计算向量x的多项式p的值。同上，计算结果在y中，同时 p o l y v a l ( p , x , E ) 还根据p o l y f i t命令给出的矩阵E返回一个误差估计向量 e r r。见help polyval和help polyfit，参见1 0 . 4节。 p o l y v a l m ( p , A ) 直接对矩阵A进行多项式计算。不是象上个命令一样对每个 元素进行多项式计算，而是计算p(A) =p1An+p2An－1+…+pn + 1I。 p o l y ( A ) 计算矩阵A的特征多项式向量。 p o l y ( x ) 给出一个长度为n+ 1的向量，其中的元素是次数为 n的多项 式的系数。这个多项式的根是长度为 n的向量x中元素。 c o m p a n ( p ) 计算带有系数 p的多项式的友矩阵 A，这个矩阵的特征多 项式为p。 r o o t s ( p ) 计算特征多项式p的根，是一个长度为n的向量，也就是方程 p(x) = 0的解。表达式p o l y ( r o o t s ( p ) ) = p为真。结果可以是复数。 c o n v ( p , q ) 计算多项式p和q的乘积，也可以认为是p和q的卷积。 [ k , r ] = d e c o n v ( p , q ) 计算多项式 p除q。k是商多项式， r是残数多项式。这个 计算等价于p和q的逆卷积