10_函数、插值与曲线拟和.pdf

第1 0章 函数、插值和曲线拟合分析 M AT L A B可以处理有估值和没有估值的多项式，还有一些强大的数值分析命令，如求零 值和最小值。M AT L A B中还有数据集合的插值、曲线拟合的命令和函数，还提到了经典的贝 赛尔( B e s s e l )函数。 10.1MAT 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的逆卷积。[u v k]= 计算p ( x ) / q ( x )的部分展开式： r e s i d u e ( p , q ) 向量p和q分别是多项式p ( x )和q ( x )的系数。得到的余数在 向量u中，极点在列向量v中，商多项式在向量k中。 [p q]=residue(u,v,x)从同上的部分展开式u、v和x计算得到多项式p和q。