celiang1.m

数学建模在的插值问题的求解

clear
x=[129.0000 140.0000 108.5000 88.0000 185.5000 195.0000 105.5000 157.5000 107.5000 77.0000 81.0000 162.0000 117.5000 162.0000];
y=[7.5000 141.5000 28.0000 147.0000 22.5000 147.5000 85.5000 -6.5000 -81.0000 3.0000 56.5000 84.0000 -38.5000 -66.5000];
z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -4 -9 -9];
X=75:0.5:200;
Y=-70:0.5:150;
Z=griddata(x,y,z,X,Y','v4');
meshz(X,Y,Z),rotate3d
figure(2),contour(X,Y,Z,[-5,-5]);grid
[i1,j1]=find(Z<-5);
for k=1:length(i1)
    Z(i1(k),j1(k))=-5;
end
figure(3),meshc(X,Y,Z)