如何拟合一个封闭的曲面

cichishia

请教大家：
已知若干个离散点的三维坐标，如何拟合成一个封闭的曲面。
- M, K" N# S3 p& g
程序：
%--------------------------------------------------------------------------
: d9 I8 M( R( BA=xlsread('PointData.xls');
2 t8 Y0 h7 y9 {" B5 O6 J+ V% 下部分网格曲面
x=A(1:85,1);   % 根据dz的大小，相应调整行的范围--及确定好哪些点属于下部分
) x, \7 e# i7 k1 U. ty=A(1:85,2);
0 N: T# m/ H$ a4 Jz=A(1:85,3);
minx = min(x);
maxx = max(x);
miny = min(y);
: K# b7 L! o; K( `# T& u1 Amaxy = max(y);
xtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
/ d! ]. j2 B/ X! ]9 }! G/ }, Mytemp=linspace(min(y),max(y),50);
2 l. t  N5 W. d6 Y$ q3 W; J* U[X,Y]=meshgrid(xtemp,ytemp);
1 }4 i8 j: o$ E3 a6 gZ=griddata(x,y,z,X,Y,'cubic');
3 v; B  M& {% K& U# |' Nfigure('numbertitle','off','name','网格图')
suRF(X,Y,Z)
shading faceted
. r' n* q% J; E" z9 |9 F
, T, ~1 W/ X& P, N%--------------------------------------------------------------------------
1 X8 O# P) U7 O4 f- h4 \! E, g% 上部分网格曲面
hold on
4 I8 i( n9 m+ _x=A(61:404,1);
y=A(61:404,2);
9 z. U  T* {& e* ^, r. U$ Iz=A(61:404,3);
minx = min(x);
maxx = max(x);
miny = min(y);
& R$ G" D" u0 c- \6 T6 `5 omaxy = max(y);
xtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
ytemp=linspace(min(y),max(y),50);
[X,Y]=meshgrid(xtemp,ytemp);
Z=griddata(x,y,z,X,Y,'cubic');
surf(X,Y,Z)
shading faceted
) F6 }8 W2 v' G7 S
- P5 ^% L4 ~# ?( r8 G, I+ Z% Qfigure('name','散点图')
plot3(A(:,1),A(:,2),A(:,3),'*')
# F$ [* i* N3 T; T: a7 n4 }
1 u1 T7 q7 n' f
我自己用上述方法分别拟合上下两部分，但是图中有缺口，如何才能得到封闭的图形呢？
* `5 }) A& `" o# P/ ?6 p4 t  S

IBB-EUT

帮你顶一下

shuddkk

给出数据，我可以试试。

小白的白