[求助]如何建立三角网？

[yanli111]

用matlab建立三角网，有已知点的三维坐标：程序如下x1=10338.864
x2=10376.258
x3=x1+5.8
x4=x1+11.24
x5=x1+29.45
x6=x1+47.66
x7=x1+67.66
x8=x1+88.98
x9=x1+106.78
x10=x1+124.16
x11=x1+129.05
x12=x1+141.61
x13=x1+145.59

y1=20806.339
y2=20949.963
y3=5.8*(y2-y1)/(x2-x1)+y1
y4=11.24*(y2-y1)/(x2-x1)+y1
y5=29.45*(y2-y1)/(x2-x1)+y1
y6=47.66*(y2-y1)/(x2-x1)+y1
y7=67.66*(y2-y1)/(x2-x1)+y1
y8=88.98*(y2-y1)/(x2-x1)+y1
y9=106.78*(y2-y1)/(x2-x1)+y1
y10=124.16*(y2-y1)/(x2-x1)+y1
y11=129.05*(y2-y1)/(x2-x1)+y1
y12=141.61*(y2-y1)/(x2-x1)+y1
y13=145.59*(y2-y1)/(x2-x1)+y1
%x=[x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13]
%y=[y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13]
z=[403.85,399.75,399.34,398.73,398.98,399.34,399.52,399.17,398.77,399.05,399.32,402.16,405.13]

%xlabel('Longitude');ylabel('Latitude');
%grid on
%tri=delaunay(y,z);

%h = trisurf(tri,y,x,z);
%grid on
%xlabel('Longitude');ylabel('Latitude');zlabel('Depth in Feet')
x1q=10240.736
x2q=10276.579
x3q=x1q+4.2
x4q=x1q+24.44
x5q=x1q+41.89
x6q=x1q+51.06
x7q=x1q+58.40
x8q=x1q+66.19
x9q=x1q+76.89
x10q=x1q+88.46
x11q=x1q+100.13
x12q=x1q+110.80
x13q=x1q+120.40
x14q=x1q+124.02
x15q=x1q+128.33
x16q=x1q+131.41


y1q=20825.512
y2q=20955.067
y3q=4.2*(y2q-y1)/(x2q-x1q)+y1q
y4q=24.44*(y2q-y1q)/(x2q-x1q)+y1q
y5q=41.89*(y2q-y1q)/(x2q-x1q)+y1q
y6q=51.06*(y2q-y1q)/(x2q-x1q)+y1q
y7q=58.40*(y2q-y1q)/(x2q-x1q)+y1q
y8q=66.19*(y2q-y1q)/(x2q-x1q)+y1q
y9q=76.89*(y2q-y1q)/(x2q-x1q)+y1q
y10q=88.46*(y2q-y1q)/(x2q-x1q)+y1q
y11q=100.13*(y2q-y1q)/(x2q-x1q)+y1q
y12q=110.80*(y2q-y1q)/(x2q-x1q)+y1q
y13q=120.40*(y2q-y1q)/(x2q-x1q)+y1q
y14q=124.02*(y2q-y1q)/(x2q-x1q)+y1q
y15q=128.33*(y2q-y1q)/(x2q-x1q)+y1q
y16q=131.41*(y2q-y1q)/(x2q-x1q)+y1q

%xq=[x1q,x2q,x3q,x4q,x5q,x6q,x7q,x8q,x9q,x10q,x11q,x12q,x13q]
%yq=[y1q,y2q,y3q,y4q,y5q,y6q,y7q,y8q,y9q,y10q,y11q,y12q,y13q]
%zq=[404.35,399.62,399.44,398.93,398.87,399.01,399.16,399.17,399.15,398.95,398.87,399.15,399.23]
x=[x1,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x2,x1q,x3q,x4q,x5q,x6q,x7q,x8q,x9q,x10q,x11q,x12q,x13q,x14q,x15q,x16q,x2q]
y=[y1,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y2,y1q,y3q,y4q,y5q,y6q,y7q,y8q,y9q,y10q,y11q,y12q,y13q,y14q,y15q,y16q,y2q]
z=[403.85,399.75,399.34,398.73,398.98,399.34,399.52,399.17,398.77,399.05,399.32,402.16,405.13,
404.35,399.62,399.44,398.93,398.87,399.01,399.16,399.17,399.15,398.95,398.87,399.15,399.23,399.24,402.18,402.30,405.19]


xlabel('Longitude');ylabel('Latitude');
tri=delaunay(x,y);
hold on
triplot(tri,x,y),hold off
figure
h = trimesh(tri,x,y,z);
xlabel('Longitude');ylabel('Latitude');zlabel('Depth in Feet')
为什么运行时出错呢？谁能帮我看看这个程序？并且进行改正？谢谢！