《MATLAB实验-6 MATLAB插值与拟合实验.docx》由会员分享,可在线阅读,更多相关《MATLAB实验-6 MATLAB插值与拟合实验.docx(6页珍藏版)》请在优知文库上搜索。
1、新乡学院数学与信息科学系试验报告试脸名称插值与拟合所属课程数学软件与试验试舱类型综合型试脸.专业信息与计算科学班级2011级1班学号姓名李欢丽指导老拜朱推生一、试验概述试触目的】学会用一维括值函数yi=interl(xo,yo,menthod,)求出函数在括值点处的函数值,和用二维函数plot()作图。用二维插值函数zi=interp2(x,y.z,x,y.,method*)求其在网格节点数据的插值,和用三维函数SUrfCQ作图.【试3原理】1.一维插值函数yi=interpl(xo,yo,x,menthod),yi被插值点处的函数值,xo.yo插值节点,x被插值点,nearest最邻近痛值,
2、linear线性插值,SPIine三次样条插值,CUbiC立方插值,缺省时,为分段段性插值.2,二维作图函数plot(x,y),x,y都是向量时,则以X中元素为横坐标,y中元素为纵坐标,且X,y长度相同。x,y都是矩阵,X列与y列结合,绘制多条平面曲线,且必同大小。3,用作网格节点数据的插值zi=interp2(x,y,z,x,y,method),zi为被插值点的函数值x.y,z点,x,y被插值点,Dielhod中nearest/为最邻近差值,linear为双线性插值,cubic为双三次势值,缺省时为双线性毋值。4,网格生成函数X,Ymeshgrid(x,y),x.y为给定的向量,X.Y是网格
3、划分后得到的网格矩阵绘制三维曲面图;三维曲面绘制函数SUrfC(X.y,z)【试酷环境】MaHabR201Ob二、试验内容HA1对函数/(.r)=r,%-5,5,分别用分段战性插做和三次样条插值作插值(其l+x中插值节点不少于20),并分别作出每种捕值方法的误差曲税.1 .分析问题通过一维插值函数yi=interpl(xo,yo,x,Inenthod)求出函数在插值点处的函数值,然后用二维函数pkt()作图。2 .问题求解functiony=IagrangMxO,y,x)%x插值节点,插值节点处的函数tfi,X要计糅函数值的点:n=length(x):%计算x的长度*=length(x:%计物
4、的长度fori=l:ms=O:z=x(i):fork=l:nP=1.0;forj=l:nif=kp=p*(z-x(j)(x(k)-x(j):。计基数endCndsp*y(k)*s;endy;s:%计算在x(i)处的函数依(拉格朗H)endXo=-5:0.5:5:y=l.(lx2);x=-5:0.1:5;y=l.(l+x.2):ylslagrange(xly.x);zl=yl-y:plotylk*)%拉格。Hi插值作图holdanplot(x,zlk+,)M(ftfflH卜的误差曲线y2=interpl(x.y,x);holdonplot(y2/r,用分段线性插脩作图z2=y2-y;Dkn(扁z2,r)%分段然性插依卜的误差曲践y3=interpl(x,y,x,spline,);plot(y3.111.,)%holdmiz3=y3-y:plotz3)%:次样条插tfl下的误差曲线holdoffaxis2):x=-5:0.1:5:yllagrange(x,y,x):zlyl、”1。1小1,|0)%拉格朗日下的误差曲她y”.+xJ2);必不行少2 .二维插值函数zi=interp2(x.y.z.x.y.method*),其中x.y为矩阵,或X为行向量,y为列向量。且X.y范围不行超过x,y3 ,y的转置可表示为y.4,输人一个矩阵的元素时,空格和分号的个数要做到不重不漏。