a=importdata('cost.txt');
x=a(:,1);
y=a(:,2);
plot(x,y,'o')
z=a(:,3);
h=-Z;
xi=1:1:4;yi=0:1:3;
[X,Y]=meshgrid(xi,yi);
H=griddata(x,y,h,X,Y,'v4');
mesh(X,Y,H);
view(-60,30);
a=importdata('shuju.txt');
x=a(:,1);
y=a(:,2);
z=a(:,3);
xx=meshgrid(x1);
yy=meshgrid(y1);
zz=meshgrid(z1);
mesh(xx,yy,zz)
x=cost(:,1);
y=cost(:,2);
plot(x,y,'o')
z=cost(:,3);
xi=0:5:100;yi=0:0.5:5;
[X,Y]=meshgrid(xi,yi);
Z=griddata(x,y,z,X,Y,'v4');
mesh(X,Y,Z);
view(-60,30);
x=[140.0 88.0 195.0 105.5 157.5 107.5 81.0 162.0 162.0];
y=[141.5 147.0 137.5 85.5 60.5 121.0 56.5 116.5 84.0 ];
plot(x,y,'o')
Z=[14 16 14 14 15 15 16 15 13 ];
h=-Z;
xi=70:3:200;yi=0:3:150;
[X,Y]=meshgrid(xi,yi);
H=griddata(x,y,h,X,Y,'v4');
mesh(X,Y,H);
view(-60,30);
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];
z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9];
cx=75:0.5:200;
cy=-70:0.5:150;
cz=griddata(x,y,z,cx,cy','cubic');
meshz(cx,cy,cz)
x=[1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4];
y=[0 1 2 3 0 1 2 3 0 1 2 3 0 1 2 3];
z=[70 60 30 0 80 75 60 30 85 80 75 70 90 85 80 75];
cx=1:1:4;
cy=0:1:4;
cz=griddata(x,y,z,cx,cy','cubic');
meshz(cx,cy,cz)
xyz=[40 2 1.4
40 5 2.5
40 7 1.4
40 9 0.9
50 2 1
50 5 3.5
50 7 0.7
50 9 1.8
60 2 0.4
60 5 0.5
60 7 4.9
60 9 4.5
70 2 0.2
70 5 1
70 7 3.5
70 8 5.6 ];
tri = delaunay(xyz(:,1), xyz(:,2));
trimesh(tri, xyz(:,1), xyz(:,2),xyz(:,3));
shading interp