第六章matlab数据分析与多项式计算课件.ppt

n n 数据统计处理数据统计处理n n 数据插值数据插值n n 曲线拟和曲线拟和n n离散傅立叶变换离散傅立叶变换n n多项式计算多项式计算2024/9/616.1 数据统计处理数据统计处理6.1.1 最大值和最小值最大值和最小值 MATLAB提供的求数据序列的最大值和最提供的求数据序列的最大值和最小值的函数分别为小值的函数分别为max和和min，，两个函数的两个函数的调用格式和操作过程类似调用格式和操作过程类似 1. 求向量的最大值和最小值求向量的最大值和最小值 求一个向量求一个向量X X的最大值的函数有两种调用的最大值的函数有两种调用 格式格式, ,分别是分别是：： (1) y=max(X)：：返回向量返回向量X的最大值存入的最大值存入y，， 如果如果X中包含复数元素，则按模取最大值中包含复数元素，则按模取最大值2024/9/62(2) [y,I]=max(X)(2) [y,I]=max(X)：：：：返回向量返回向量返回向量返回向量X X X X的最大值存入的最大值存入的最大值存入的最大值存入y y，，，，最大最大最大最大值的序号存入值的序号存入值的序号存入值的序号存入I I，，，，如果如果如果如果X X中包含复数元素，则按模取中包含复数元素，则按模取中包含复数元素，则按模取中包含复数元素，则按模取最大值。

最大值üü求向量求向量求向量求向量X X的最小值的函数是的最小值的函数是的最小值的函数是的最小值的函数是min(X),min(X),用法和用法和用法和用法和max(X)max(X)完全相同完全相同完全相同完全相同例例6-1 求向量求向量x的最大值的最大值2024/9/632. 求矩阵的最大值和最小值求矩阵的最大值和最小值 求矩阵求矩阵A的最大值的函数有的最大值的函数有3种调用格式，种调用格式，分别是：分别是：(1)max(A)：：返回一个行向量，向量的第返回一个行向量，向量的第i个元个元素是矩阵素是矩阵A的第的第i列上的最大值列上的最大值2)[Y,U]=max(A)：：返回行向量返回行向量Y和和U，，Y向量向量记录记录A的每列的最大值的每列的最大值，，U向量记录每列最向量记录每列最大值的行号大值的行号2024/9/64(3) max(A,[ ],dim)：：dim取取1或或2dim取取1时时，，该函数和该函数和max(A)完全相同完全相同；；dim取取2时时，，该该函数返回一个列向量函数返回一个列向量,其第其第i个元素是个元素是A矩阵矩阵的第的第i行上的最大值行上的最大值 求最小值的函数是求最小值的函数是min，，其用法和其用法和max完全完全相同。

相同例例6-2 分别求分别求3×4矩阵矩阵x中各列和各行元素中的中各列和各行元素中的最大值，并求整个矩阵的最大值最大值，并求整个矩阵的最大值.2024/9/653. 两个向量或矩阵对应元素的比较两个向量或矩阵对应元素的比较 函数函数max和和min还能对两个同型的向量或矩还能对两个同型的向量或矩阵进行比较阵进行比较,调用格式为调用格式为：：(1)U=max(A,B)：：A,B是两个同型的向量或矩是两个同型的向量或矩 阵，结果阵，结果U是与是与A,B同型的向量或矩阵同型的向量或矩阵，，U的每个元素等于的每个元素等于A,B对应元素的较大者对应元素的较大者2)U=max(A,n)：：n是一个标量，结果是一个标量，结果U是与是与A同型的向量或矩阵同型的向量或矩阵,U的每个元素等于的每个元素等于A对应对应元素和元素和n中的较大者中的较大者 min函数的用法和函数的用法和max完全相同完全相同例例6-3 求两个求两个2×3矩阵矩阵x, y所有同一位置上的较所有同一位置上的较大元素构成的新矩阵大元素构成的新矩阵p2024/9/666.1.2 求和与求积求和与求积 数据序列求和与求积的函数是数据序列求和与求积的函数是sum和和prod，，其使用方法类似。

设其使用方法类似设X是一个向量是一个向量,A是一个是一个矩阵矩阵,函数的调用格式为：函数的调用格式为： sum(X)：：返回向量返回向量X各元素的和各元素的和 prod(X)：：返回向量返回向量X各元素的乘积各元素的乘积 sum(A)：：返回一个行向量返回一个行向量,其第其第i个元素是个元素是A的第的第i列的元素和列的元素和2024/9/67 prod(A)：：返回一个行向量返回一个行向量,其第其第i个元素是个元素是A的第的第i列的元素乘积列的元素乘积 sum(A,dim)：：当当dim为为1时时,该函数等同于该函数等同于sum(A)；；当当dim为为2时时,返回一个列向量，其返回一个列向量，其第第i个元素是个元素是A的第的第i行的各元素之和行的各元素之和 prod(A,dim)：：当当dim为为1时，该函数等同于时，该函数等同于prod(A)；；当当dim为为2时，返回一个列向量，时，返回一个列向量，其第其第i个元素是个元素是A的第的第i行的各元素乘积行的各元素乘积例例6-4 求矩阵求矩阵A的每行元素的乘积和全部元素的每行元素的乘积和全部元素的乘积2024/9/686.1.3 平均值和中值平均值和中值 求数据序列平均值的函数是求数据序列平均值的函数是求数据序列平均值的函数是求数据序列平均值的函数是meanmean，，，，求数据序列中值求数据序列中值求数据序列中值求数据序列中值的函数是的函数是的函数是的函数是medianmedian。

两个函数的调用格式为两个函数的调用格式为两个函数的调用格式为两个函数的调用格式为：：：： mean(X)mean(X)：：：：返回向量返回向量返回向量返回向量X X的算术平均值的算术平均值的算术平均值的算术平均值 median(X)median(X)：：：：返回向量返回向量返回向量返回向量X X的中值的中值的中值的中值 mean(A)mean(A)：：：：返回一个行向量返回一个行向量返回一个行向量返回一个行向量, ,其第其第其第其第i i个元素是个元素是个元素是个元素是A A的第的第的第的第i i列列列列的算术平均值的算术平均值的算术平均值的算术平均值 median(A)median(A)：：：：返回一个行向量返回一个行向量返回一个行向量返回一个行向量, ,其第其第其第其第i i个元素是个元素是个元素是个元素是A A的第的第的第的第i i列的中值列的中值列的中值列的中值 mean(A,dim)mean(A,dim)：：：：当当当当dimdim为为为为1 1时时时时, ,该函数等同于该函数等同于该函数等同于该函数等同于mean(A)mean(A)；；；；当当当当dimdim为为为为2 2时，返回一个列向量，其第时，返回一个列向量，其第时，返回一个列向量，其第时，返回一个列向量，其第i i个元素是个元素是个元素是个元素是A A的的的的第第第第i i行的算术平均值。

行的算术平均值行的算术平均值行的算术平均值 median(A,dim)median(A,dim)：：：：当当当当dimdim为为为为1 1时，该函数等同于时，该函数等同于时，该函数等同于时，该函数等同于median(A)median(A)；；；；当当当当dimdim为为为为2 2 2 2时，返回一个列向量，其第时，返回一个列向量，其第时，返回一个列向量，其第时，返回一个列向量，其第i i个元素是个元素是个元素是个元素是A A的第的第的第的第i i行的中值行的中值行的中值行的中值例例例例6-5 6-5 求矩阵求矩阵求矩阵求矩阵A A A A平均值和中值平均值和中值平均值和中值平均值和中值2024/9/696.1.4 累加和与累乘积累加和与累乘积(examp6_6.m) 在在在在MATLABMATLAB中，使用中，使用中，使用中，使用cumsumcumsum和和和和cumprodcumprod函数能方函数能方函数能方函数能方便地求得向量和矩阵元素的累加和与累乘积向量，便地求得向量和矩阵元素的累加和与累乘积向量，便地求得向量和矩阵元素的累加和与累乘积向量，便地求得向量和矩阵元素的累加和与累乘积向量，函数的调用格式为：函数的调用格式为：函数的调用格式为：函数的调用格式为： cumsum(X)cumsum(X)：：：：返回向量返回向量返回向量返回向量X X累加和向量累加和向量累加和向量累加和向量。

cumprod(X)cumprod(X)：：：：返回向量返回向量返回向量返回向量X X累乘积向量累乘积向量累乘积向量累乘积向量 cumsum(A)cumsum(A)：：：：返回一个矩阵，其第返回一个矩阵，其第返回一个矩阵，其第返回一个矩阵，其第i i列是列是列是列是A A的第的第的第的第i i列的列的列的列的累加和向量累加和向量累加和向量累加和向量 cumprod(A)cumprod(A)：：：：返回一个矩阵，其第返回一个矩阵，其第返回一个矩阵，其第返回一个矩阵，其第i i列是列是列是列是A A的第的第的第的第i i列列列列的累乘积向量的累乘积向量的累乘积向量的累乘积向量 cumsum(A,dim)cumsum(A,dim)：：：：当当当当dimdim为为为为1 1 1 1时，该函数等同于时，该函数等同于时，该函数等同于时，该函数等同于cumsum(A)cumsum(A)；；；；当当当当dimdim为为为为2 2 2 2时，返回一个矩阵，其第时，返回一个矩阵，其第时，返回一个矩阵，其第时，返回一个矩阵，其第i i行行行行是是是是A A的第的第的第的第i i行的累加和向量。

行的累加和向量行的累加和向量行的累加和向量 cumprod(A,dim)cumprod(A,dim)：：：：当当当当dimdim为为为为1 1 1 1时，该函数等同于时，该函数等同于时，该函数等同于时，该函数等同于cumprod(A)cumprod(A)；；；；当当当当dimdim为为为为2 2 2 2时，返回一个向量，其第时，返回一个向量，其第时，返回一个向量，其第时，返回一个向量，其第i i行是行是行是行是A A的第的第的第的第i i行的累乘积向量行的累乘积向量行的累乘积向量行的累乘积向量2024/9/6106.1.5 标准方差与相关系数标准方差与相关系数 1. 1. 求标准方差求标准方差求标准方差求标准方差 在在在在MATLABMATLAB中，提供了计算数据序列的标准方中，提供了计算数据序列的标准方中，提供了计算数据序列的标准方中，提供了计算数据序列的标准方差的函数差的函数差的函数差的函数stdstd对于向量对于向量对于向量对于向量X X，，，，std(X)std(X)返回一个标准方返回一个标准方返回一个标准方返回一个标准方差。

对于矩阵差对于矩阵差对于矩阵差对于矩阵A A，，，，std(A)std(A)返回一个行向量，它的各返回一个行向量，它的各返回一个行向量，它的各返回一个行向量，它的各个元素便是矩阵个元素便是矩阵个元素便是矩阵个元素便是矩阵A A各列或各行的标准方差各列或各行的标准方差各列或各行的标准方差各列或各行的标准方差stdstd函数函数函数函数的一般调用格式为的一般调用格式为的一般调用格式为的一般调用格式为：：：： Y=std(A,flag,dim)Y=std(A,flag,dim) 其中其中其中其中dimdim取取取取1 1或或或或2 2当dim=1dim=1时，求各列元素的标时，求各列元素的标时，求各列元素的标时，求各列元素的标准方差；当准方差；当准方差；当准方差；当dim=2dim=2时，则求各行元素的标准方差时，则求各行元素的标准方差时，则求各行元素的标准方差时，则求各行元素的标准方差flagflag取取取取0 0或或或或1 1，，，，当当当当flag=0flag=0时，按时，按时，按时，按σ1σ1所列公式计算标准所列公式计算标准所列公式计算标准所列公式计算标准方差，当方差，当方差，当方差，当flag=1flag=1时，按时，按时，按时，按σ2σ2所列公式计算标准方差。

所列公式计算标准方差所列公式计算标准方差所列公式计算标准方差缺省缺省缺省缺省flag=0flag=0，，，，dim=1dim=1例例6-7 对二维矩阵对二维矩阵A，，从不同维方向求出其标准从不同维方向求出其标准方差方差2024/9/6112. 相关系数相关系数 MATLAB提供了提供了corrcoef函数，可以求出数函数，可以求出数据的相关系数矩阵据的相关系数矩阵corrcoef函数的调用格函数的调用格式为：式为： corrcoef(X)：：返回从矩阵返回从矩阵X形成的一个相关形成的一个相关系数矩阵此相关系数矩阵的大小与矩阵系数矩阵此相关系数矩阵的大小与矩阵X一样它把矩阵一样它把矩阵X的每列作为一个变量，然的每列作为一个变量，然后求它们的相关系数后求它们的相关系数 corrcoef(X,Y)：：在这里在这里，，X,Y是向量，它们是向量，它们与与corrcoef([X,Y])的作用一样的作用一样例例6-8 生成满足正态分布的生成满足正态分布的3×3随机矩阵随机矩阵,然后求然后求各列元素的均值和标准方差，再求这各列元素的均值和标准方差，再求这3列随机列随机数据的相关系数矩阵数据的相关系数矩阵。

2024/9/6126.1.6 排序排序 MATLAB中对向量中对向量X X排序函数是排序函数是sort(X)，，函数返回一个对函数返回一个对X中的元素按升序排列的新中的元素按升序排列的新向量向量 sort函数也可以对矩阵函数也可以对矩阵A的各列或各行重新的各列或各行重新排序，其调用格式为：排序，其调用格式为： [Y,I]=sort(A,dim) 其中其中dim指明对指明对A的列还是行进行排序若的列还是行进行排序若dim=1，，则按列排；若则按列排；若dim=2，，则按行排则按行排Y是排序后的矩阵，而是排序后的矩阵，而I记录记录Y中的元素在中的元素在A中位置中位置例例6-9 对二维矩阵做各种排序对二维矩阵做各种排序2024/9/6136.2 数据插值数据插值6.2.1 一维数据插值一维数据插值n ny1=interp1(x,y,x1)返回在插值向量返回在插值向量xi处的函处的函数向量数向量y1，它是根据向量，它是根据向量x和和y插值而来若插值而来若y是矩阵，则对是矩阵，则对y每一列进行插值，如每一列进行插值，如xi中元素中元素不在不在x内，返回内，返回NaNn ny1=interp1(y,x1)省略省略x，表示，表示x=1:N，此时，此时N为向量为向量y的长度或为矩阵的长度或为矩阵y的行数。

的行数n ny1=interp1(x,y,x1,’method’)表示用表示用method指定的插值方法进行插值指定的插值方法进行插值2024/9/614Method可取如下的值：可取如下的值：n n‘linear’线性插值线性插值n n‘nearest’最近插值最近插值n n‘spline’三次样条插值三次样条插值n n‘cubic’三次插值三次插值n nMethod默认值为线性插值，上述插值要求默认值为线性插值，上述插值要求向量向量x单调2024/9/615例：例：>> x=[1 2 4 6 8 9 10 13 15 16];>> y=[5 7 8 10 13 14 15 17 19 20];>> x1=[1.2 2.1 3];>> y1=interp1(x,y,x1)y1 = 5.4000 7.0500 7.50002024/9/616>> x=[1 2 4 6 8 9 10 13 15 16];>> y=[5 7 8 10 13 14 15 17 19 20];>> x1=[1.2 2.1 3];>> y1=interp1(x,y,x1,'linear')y1 = 5.4000 7.0500 7.50002024/9/617>> x=[1 2 4 6 8 9 10 13 15 16];>> y=[5 7 8 10 13 14 15 17 19 20];>> x1=[1.2 2.1 3];>> y1=interp1(x,y,x1,'nearest')y1 = 5 7 82024/9/618>> x=[1 2 4 6 8 9 10 13 15 16];>> y=[5 7 8 10 13 14 15 17 19 20];>> x1=[1.2 2.1 3];>> y1=interp1(x,y,x1,'spline')y1 = 5.5529 7.1110 7.67472024/9/619>> x=[1 2 4 6 8 9 10 13 15 16];>> y=[5 7 8 10 13 14 15 17 19 20];>> x1=[1.2 2.1 3];>> y1=spline(x,y,x1)y1 = 5.5529 7.1110 7.67472024/9/620>> x=[1 2 4 6 8 9 10 13 15 16];>> y=[5 7 8 10 13 14 15 17 19 20];>> x1=[1.2 2.1 3];>> y1=interp1(x,y,x1,'cubic')y1 = 5.5006 7.0814 7.54762024/9/621例例6-10 用不同的插值方法计算。

用不同的插值方法计算 MATLAB中有一个专门的中有一个专门的3次样条插值函数次样条插值函数y1=spline(x,y,x1)，，其功能及使用方法与函数其功能及使用方法与函数y1=interp1(x,y,x1,‘spline’)完全相同完全相同例例6-11 某观测站测得某日某观测站测得某日6:006:00时至时至18:0018:00时之时之间每隔间每隔2 2小时的室内外温度小时的室内外温度(℃)(℃)，用，用3 3次样条次样条插值分别求得该日室内外插值分别求得该日室内外6:306:30至至17:3017:30时之时之间每隔间每隔2 2小时各点的近似温度小时各点的近似温度(℃)(℃) 2024/9/6226.2.2 6.2.2 二维数据插值二维数据插值二维数据插值二维数据插值n n z1=interp2(x,y,z,x1,y1)返回在插值向量返回在插值向量x1、、y1处的函数值向量，它是根据向量处的函数值向量，它是根据向量x、、y与与z插值而来，如果插值而来，如果x1、、y1有元素不在有元素不在x、、y范围内，则返回范围内，则返回NaNn nz1=interp2(z,x1,y1)省略省略x、、y，表示，表示x=1:N，，y=1:M, [M,N] = size(z).。

n nzi=interp2(x,y,z,xi,yi,’method’)表示用表示用method指定的插值方法进行插值指定的插值方法进行插值Method取值同上，要求取值同上，要求x与与y都单调且具有都单调且具有相同格式相同格式2024/9/623 曲线拟合涉及到两个基本问题：什么是曲线拟合涉及到两个基本问题：什么是最佳拟合？用什么样的曲线进行拟合可以最佳拟合？用什么样的曲线进行拟合可以用许多方法定义最佳拟合，而且存在无穷数用许多方法定义最佳拟合，而且存在无穷数目的曲线目的曲线 当最佳拟合定义为数据最小误差平方和，当最佳拟合定义为数据最小误差平方和，所用的曲线限定为多项式时，拟合曲线就相所用的曲线限定为多项式时，拟合曲线就相对简单 数学上称为多项式的最小二乘曲线拟合数学上称为多项式的最小二乘曲线拟合6.3 曲线拟合曲线拟合2024/9/624 在在在在MATLABMATLAB中，用中，用中，用中，用polyfitpolyfit函数来求得最小二乘函数来求得最小二乘函数来求得最小二乘函数来求得最小二乘拟合多项式的系数，再用拟合多项式的系数，再用拟合多项式的系数，再用拟合多项式的系数，再用polyvalpolyval函数按所得的多函数按所得的多函数按所得的多函数按所得的多项式计算所给出的点上的函数近似值。

项式计算所给出的点上的函数近似值项式计算所给出的点上的函数近似值项式计算所给出的点上的函数近似值 polyfitpolyfit函数的调用格式为：函数的调用格式为：函数的调用格式为：函数的调用格式为： [P,S]=polyfit(X,Y,m)[P,S]=polyfit(X,Y,m) 函数根据采样点函数根据采样点函数根据采样点函数根据采样点X X和采样点函数值和采样点函数值和采样点函数值和采样点函数值Y Y，，，，产生一产生一产生一产生一个个个个mm次多项式次多项式次多项式次多项式P P及其在采样点的误差向量及其在采样点的误差向量及其在采样点的误差向量及其在采样点的误差向量S S其中其中其中其中X,YX,Y是两个是两个是两个是两个等长等长等长等长的向量的向量的向量的向量，，，，P P是一个长度为是一个长度为是一个长度为是一个长度为m+1m+1的向的向的向的向量量量量，，，，P P的元素为多项式系数的元素为多项式系数的元素为多项式系数的元素为多项式系数例例例例6-146-14求如下给定数据的拟合曲线，求如下给定数据的拟合曲线，x=[0.5,1.0,1.5,2.0,2.5,3.0]，， y=[1.75,2.45,3.81,4.80,7.00,8.60]。

2024/9/6256.5 多项式计算多项式计算2024/9/6266.5.1 多项式的四则运算多项式的四则运算 1. 多项式的加减运算多项式的加减运算 对于次数相同的若干个多项式，可直接对对于次数相同的若干个多项式，可直接对多项式系数向量进行加、减的运算如果多多项式系数向量进行加、减的运算如果多项式的次数不同，则应该把低次的多项式系项式的次数不同，则应该把低次的多项式系数不足的高次项用零补足，使同式中的各多数不足的高次项用零补足，使同式中的各多项式具有相同的次数项式具有相同的次数 2024/9/6272. 多项式乘法运算多项式乘法运算 函数函数conv(P1,P2)用于求多项式用于求多项式P1和和P2的的乘积这里，乘积这里，P1、、P2是两个多项式系数向是两个多项式系数向量例例6-16 求多项式求多项式x4+8x3-10与多项式与多项式2x2-x+3的的乘积2024/9/6283. 多项式除法多项式除法 函数函数[Q,r]=deconv(P1,P2)用于对多项式用于对多项式P1和和P2作除法运算作除法运算其中其中Q返回多项式返回多项式P1除除以以P2的商式的商式，，r返回返回P1除以除以P2的余式。

这里，的余式这里，Q和和r仍是多项式系数向量仍是多项式系数向量 deconv是是conv的逆函数，即有的逆函数，即有P1=conv(P2,Q)+r例例6-17 求多项式求多项式x4+8x3-10除以多项式除以多项式2x2-x+3的结果2024/9/6296.5.2 多项式的导数多项式的导数 对多项式求导数的函数是对多项式求导数的函数是：： p=polyder(P)：：求多项式求多项式P的导数的导数 p=polyder(P,Q)：：求求P·Q的导数的导数 [p,q]=polyder(P,Q)：：求求P/Q的导数，导数的的导数，导数的分子存入分子存入p，，分母存入分母存入q 上述函数中，参数上述函数中，参数P,Q是多项式的向量表示是多项式的向量表示,结果结果p,q也是多项式的向量表示也是多项式的向量表示例例6-18 求有理分式的导数求有理分式的导数2024/9/6306.5.3 多项式的求值多项式的求值 MATLAB提供了两种求多项式值的函数提供了两种求多项式值的函数：：polyval与与polyvalm，，它们的输入参数均为它们的输入参数均为多项式系数向量多项式系数向量P和自变量和自变量x。

两者的区别两者的区别在于前者是代数多项式求值，而后者是矩在于前者是代数多项式求值，而后者是矩阵多项式求值阵多项式求值2024/9/6311. 代数多项式求值代数多项式求值 polyval函数用来求代数多项式的值，其调函数用来求代数多项式的值，其调用格式为：用格式为： Y=polyval(P,x) 若若x为一数值，则求多项式在该点的值；若为一数值，则求多项式在该点的值；若x为向量或矩阵，则对向量或矩阵中的每个为向量或矩阵，则对向量或矩阵中的每个元素求其多项式的值元素求其多项式的值例例6-19 已知多项式已知多项式x4+8x3-10，，分别取分别取x=1.2和和一个一个2×3矩阵为自变量计算该多项式的值矩阵为自变量计算该多项式的值2024/9/632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))例例6-20 仍以多项式仍以多项式x4+8x3-10为例，取一个为例，取一个2×2矩阵为自变量分别用矩阵为自变量分别用polyval和和polyvalm计算计算该多项式的值。

该多项式的值2024/9/6336.5.4 多项式求根多项式求根 n次多项式具有次多项式具有n个根，当然这些根可能是个根，当然这些根可能是实根，也可能含有若干对共轭复根实根，也可能含有若干对共轭复根MATLAB提供的提供的roots函数用于求多项式的函数用于求多项式的全部根，其调用格式为：全部根，其调用格式为： x=roots(P) 其中其中P为多项式的系数向量，求得的根赋给为多项式的系数向量，求得的根赋给向量向量x，即，即x(1),x(2),…,x(n)分别代表多项式分别代表多项式的的n个根个根例例6-21 求多项式求多项式x4+8x3-10的根的根2024/9/634 若已知多项式的全部根，则可以用若已知多项式的全部根，则可以用poly函数函数建立起该多项式，其调用格式为：建立起该多项式，其调用格式为： P=poly(x) 若若x为具有为具有n个元素的向量，则个元素的向量，则poly(x)建立以建立以x为为其根的多项式，且将该多项式的系数赋其根的多项式，且将该多项式的系数赋给向量给向量P例例6-22 已知已知 f(x) (1) 计算计算f(x)=0 的全部根的全部根。

(2) 由方程由方程f(x)=0的根构造一个多项式的根构造一个多项式g(x)，，并与并与f(x)进行对比进行对比 2024/9/635多项式拟合多项式拟合 多项式拟合用多项式拟合用Polyfit(x,y,n)来实现，来实现，n是拟是拟合多项式的阶次合多项式的阶次例例6-23>> x=linspace(0,2*pi,100);>> y=cos(x);>> t=polyfit(x,y,6);>> y1= polyval(t,x);>> plot(x,y,'go',x,y1,'b--')2024/9/6362024/9/637>> x=linspace(0,2*pi,100);>> y=cos(x);>> t=polyfit(x,y,3);>> y1=polyval(t,x);>> plot(x,y,'go',x,y1,'b--')2024/9/6382024/9/639。