第二课堂插值与拟合4.14x.pptVIP

在某海域测得一些点(x, y)处的水深 z(单位：英尺)由下表给出, 水深数据是在低潮时测得的. 船的吃水深度为 5 英尺, 问在矩形区域 (75,200)?(?50,150) 里的哪些地方船要避免进入. 例4 水道测量数据(AMCM 86A题) 水道水深测量数据(单位：英尺) 9 4 9 8 8 9 9 z ?38.5 84.0 ?66.5 56.5 3.0 ?81.0 ?6.5 y 117.5 162.0 162.0 81.0 77.0 107.5 157.5 x 8 8 6 8 6 8 4 z 85.5 137.5 22.5 147.0 28.0 141.5 7.5 y 105.5 195.0 185.5 88.0 108.5 140.0 129.0 x 1. 输入插值点数据 2. 在矩形区域(75,200)?(?50,150)做二维插值 3. 做海底地形图 4. 作出水深小于5的海域范围, 即z=5的等高线. 程序见hd.m clc clear all x=[129 140 108.5 88.5 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; y=[7.5 141.5 28 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -38.5]; z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9]; cx=75:1:200; cy=-50:1:150; cz=griddata(x,y,z,cx,cy,v4); mesh(cx,cy,cz) xlabel(X),ylabel(Y),zlabel(Z) figure(2) contour(cx,cy,cz,[-5 -5],r); grid hold on plot(x,y,*) xlabel(X),ylabel(Y) (75，200)?(?50，150) 范围内的海底地形图 水深不超过5米的危险区域的平面图 其调用格式与interp2很类似, 这里就不再重述了,具体调用格式读者可利用help来得. §2.3 高维插值 高维插值：三维插值及三维以上的插值. 命令： interp3(三维插值)、interpn(n维插值) griddata3(三维插值)、griddatan(n维插值) 网格数据 散乱数据 ndgrid(n维数据网格) §3 拟合问题的Matlab实现 3.1 一维数据拟合 1. 多项式拟合 2. 一般曲线拟合 如， 如， 调用格式：a= polyfit(x, y, n)， 求出已知数据x,y的n阶拟合多项式 的系数a 调用格式：a= lsqcurvefit(‘fun’, a0, x,y) 做各种类型的曲线拟合, 如非线性函数. 3.1 一维数据拟合 1. 多项式拟合 调用格式：a= polyfit(x, y, n)， 2. 一般曲线拟合 如， 如， 求出已知数据x,y的n阶拟合多项式 的系数a 调用格式：a= lsqcurvefit(fun, a0, x,y) 或 a= lsqnonlin(fun, a0, x,y) 例5 对下面一组数据进行2次多项式拟合，数据见下表： 8.60 7.00 4.80 3.81 2.45 1.75 y 3.0 2.5 2.0 1.5 1.0 0.5 x Matlab程序： clc clear all 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]; a=polyfit(x,y,2) y1=polyval(a,x); RR=1-sum((y-y1).^2)/sum((y-mean(y)).^2) 计算多项式 在x处的值 拟合效果分析 程序见nihe1.m a= polyfit(x, y, n) 程序 RR=1-sum((y-y1).^2)/sum((y-mean(y)).^2) 曲线拟合效果分析——可决系数 y1为预测值 2次多项式的拟合效果 例6 给定下列数据： 0.72 1.50 1.06 0.86 0.84 0.95 y 0.3 -0.2 0 0.15 0.2 0.1 x 考虑非线性拟合, 用函数 . Matlab程序： clc clear all fun=inline(a(1)*exp(a(2)*x),a,x); x=[0.1 0.2 0.15 0 -0.2 0.3]; y=[0.95 0.84 0.86 1.06 1.50 0.72]; a=lsqcurvefit(fun,[0 0],