第三章_MATLAB数值计算.pptVIP

例 9 用一个8次多项式在区间[0,3]内逼近函数sin(x). 程序：x=0:0.1:3; y=sin(x); p=polyfit(x,y,8); y1=polyval(p,x); plot(x,y1,’k*’,x,y,’k-’) *曲线是p(x),实线是sin(x)。 数值插值 插值是在已知数据之间计算估计值的过程，在信号处理和图形分析中，插值运算具有广泛的应用。 插值运算分为：一维插值、二维 插值和多维插值。其中一维插值是最 常用的插值运算。 一维插值 格式：yi=interp1(x,y,xi,method) 其中, x为自变量的取值范围； y为函数值矢量； xi为插值点的自变量矢量； yi为自变量xi所对应的插值; method为插值方法选项: 1) method=‘nearest’ 则为邻近点插值； 2) method=‘linear’,则为线性插值; 该方法是一维插值函数interp1的缺省法。 3) method=‘spline’,则为三次样条插值; 4) method=‘pchip’或‘cubic’， 则为立方插值。共4种插值方法。 例10 一维插值函数4种插值方法的比较 程序如下： x0=0:10; %x为函数y的自变量取值范围,即:0,1,2,3,4,5,6,7,8,9,10 y0=sin(x0) % y为函数值矢量 xi=0:0.25:10; % xi为插值点的自变量矢量 % 下面将插值方法定义为单元数组 a={‘nearest’,‘linear’,‘spline’,‘cubic’} ； %即分别为邻近点插值法,线性插值,三次样条插值,立方插值 % 下面将X轴标识定义为单元数组 b={‘(a) method=nearest’, ‘(b) method=linear’,‘(c) method=spline’, ‘(d) method=cubic’}; %注意此语句只能占一个程序行 for i=1:4 yi=interp1(x0,y0,xi,a{i}); subplot(2,2,i) %注意subplot(m,n,i) 是在一个图形窗口中建立 m*n 个子窗口,当i=2,是指第2个窗口. plot(x0,y0,‘ro’,xi,yi,‘b’) % 在一个图形窗口绘制多幅图形 xlabel(b(i)) % 将x轴上标注由b(i)所定义的内容 grid %绘制二维网格图 end 从曲线可见，三次样条插值方法的插值效果最好，即三次样条插值的平滑性最好。 但是不能认为什么情况下都是这样的，插值方法的好坏取决于被插值函数，没有一种对所有函数都是最好的插值方法。 二维插值 格式: ZI=interp2(X,Y,Z,XI,YI,method) 其中,X和Y为自变量组成的数组； Z为自变量X,Y的二维函数数组; XI和YI为插值点的自变量数组； ZI为自变量XI、YI所对应的插值; method为插值方法选项： 1) method=‘nearest’, 则为邻近点插值； 2) method=‘linear’, 则为双线性插值； 该方法是二维插值函数interp2的缺省方法。 3) method=‘spline’， 则为三次样条插值； 4) method=‘pchip’或‘cubic’ 则为二重立方插值。 共4种二维插值方法。 例11 4种二维插值方法比较 程序如下： [x,y,z]=peaks(7) %生成双峰函数z(x和y为自变量), 同时,并产生x、y和z三个都是7×7阶的方阵。 mesh(x,y,z) % 绘制双峰函数三维网格曲面图,类似二维作图函数plot(x,y)的作用. title(‘(1) mesh(x,y,z)---绘制双峰函数网格图) 第三章 MATLAB数值计算 3.1.1 特殊矩阵 1．通用的特殊矩阵常用的产生通用特殊矩阵的函数有：zeros：产生全0矩阵(零矩阵)