第7章-插值问题实验.pptVIP

第七章 一维、二维插值方法实验 一维插值方法 二维插值方法 l一维插值方法 1.引例1 已经测得在北纬32.3? 海洋不同深度处的温度如下表： 表7.1.1 根据这些数据，我们希望能合理地估计出其它深度（如500米、600米、1000米…）处的水温。 解决这个问题，可以通过构造一个与给定数据相适应的函数来解决，这是一个被称为插值的问题。 l一维数据插值的Matlab实现 l二维插值方法 二维规则插值 二维散乱插值 l二维数据插值的Matlab实现 （4）求解的Matlab程序如下： clear; x=[129,140,108.5,88,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]; subplot(2,2,1),plot(x,y,+),title(测量点分布图); z=[-4,-8,-6,-8,-6,-8,-8,-9,-9,-8,-8,-9,-4,-9]; [fnodes,minnq,rnw,rnq,ifail]=e01sef(x,y,z); nx=100; px=linspace(75,200,nx); ny=200; py=linspace(-50,150,ny); for i=1:ny for j=1:nx [pf(i,j),ifail]=e01sff(x,y,z,rnw,fnodes,px(j),py(i)); end end subplot(2,2,2),meshz(px,py,pf+5),title((75，200)x(-50，150) 范围内的海底地形图); [a,b]=find(pf=-5); amin=min(a); amax=max(a); bmin=min(b); bmax=max(b); xmin=75+((200-75)/100)*bmin xmax=75+((200-75)/100)*bmax ymin=-50+((150+50)/200)*amin ymax=-50+((150+50)/200)*amax for k=1:length(b) i0(k)=75+((200-75)/100)*b(k); end for k=1:length(a) j0(k)=-50+((150+50)/200)*a(k); end subplot(2,2,3),plot(i0,j0,+),title(水深不超过5米的危险区域的平面图); [i1,j1]=find(pf-5); for k=1:length(i1) pf(i1(k),j1(k))=-5; end subplot(2,2,4),meshc(px,py,pf),title(水深不超过5米的危险区域的海底地貌图); rotate3d X=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; Y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]; plot(X,Y,*); Z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4??-9]; %a=linspace(75,200,100); %b=linspace(-50,150,200); [x,y]=meshgrid(75:.5:200,-50:.5:150); z=griddata(X,Y,Z,x,y,cubic);  figure(2); meshz(x,y,z+5); figure(3); contour(x,y,z,[-5,5]);figure(4); meshz(x,y,z);rotate3d 同样可以作出海底地貌图，只是效果不是太好 * 深度x (m) 466 714 950 1422 1634 水温y (C?) 7.04 4.28 3.40 2.54 2.13 2．插值问题的基本提法 对于给定的函数表 x x0 x1 … xn y = f (x) y0 y1 … yn 其中f (x) 在区间 [a, b] 上连续，x0，x1，…，xn为 [a, b] 上n ? 1个互不相同的点，要求在一个性质优良、便于计算的函数类 {P(x)} 中，选出一个使 P(xi ) = yi，i = 0, 1, …, n (7.