首先先感謝Rasin神人幫忙。我用google找了一下，原來版上以前也有神人分享
類似問題，回饋給大家: (harry901提供)
我附上的code畫圖前半段都不變，畫圖部分:
tri = delaunay(vecx,vecy);
trisurf(tri,vecx,vecy,vecz);
或者
trimesh(tri,vecx,vecy,vecz);
唯一要注意的是，NaN要先處理掉 => 用find
此為三角剖分法。