如何拟合一个封闭的曲面

cichishia

请教大家：
! W- B5 T, b8 j* v( F  f已知若干个离散点的三维坐标，如何拟合成一个封闭的曲面。

程序：
/ h! R! h( g2 ?. d8 {; V9 N%--------------------------------------------------------------------------
A=xlsread('PointData.xls');
9 j6 L! F8 s* q% 下部分网格曲面
x=A(1:85,1);   % 根据dz的大小，相应调整行的范围--及确定好哪些点属于下部分
& c3 @5 f: l) ^' Z$ Q( m! A9 ]y=A(1:85,2);
# G1 B) B+ J' {, \& F, ~z=A(1:85,3);
minx = min(x);
maxx = max(x);
miny = min(y);
; H! \: s& W/ H% `4 C* h8 Kmaxy = max(y);
! U1 k- F$ K5 G: p4 _  @xtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
ytemp=linspace(min(y),max(y),50);
[X,Y]=meshgrid(xtemp,ytemp);
# i2 F3 O0 o" E# a" H, ^9 B2 i: W, IZ=griddata(x,y,z,X,Y,'cubic');
figure('numbertitle','off','name','网格图')
8 w/ g; O6 `; z3 m( z* I8 jsuRF(X,Y,Z)
shading faceted
: N3 @! P1 l' `1 p/ D9 w* z5 k% Q6 v
  B1 T- {7 E: o" T1 F0 a/ D%--------------------------------------------------------------------------
% 上部分网格曲面
/ ?0 D* c5 ^; ^hold on
x=A(61:404,1);
5 T1 B4 p, g' p& J1 Z' x; `y=A(61:404,2);
+ h1 I" F8 ~7 R% [( c' A% Mz=A(61:404,3);
minx = min(x);
maxx = max(x);
miny = min(y);
2 M9 D" y& `) B- W7 ]  Imaxy = max(y);
xtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
* q3 q% V( a; W1 }ytemp=linspace(min(y),max(y),50);
# u: I$ O4 R9 u/ l: V, [& c) i1 R6 v! ][X,Y]=meshgrid(xtemp,ytemp);
3 t; j2 Q$ s' i5 TZ=griddata(x,y,z,X,Y,'cubic');
+ {, N$ j. j$ S: b) @3 r0 y- d" ksurf(X,Y,Z)
; @& I$ h- i; e) Oshading faceted
2 r% g/ t0 R$ r. N; T
figure('name','散点图')
plot3(A(:,1),A(:,2),A(:,3),'*')
3 z6 F2 w+ L# l/ k& H
4 ]1 d* ^; A$ }- c# N" e5 D. J8 d
我自己用上述方法分别拟合上下两部分，但是图中有缺口，如何才能得到封闭的图形呢？

IBB-EUT

帮你顶一下

shuddkk

给出数据，我可以试试。

小白的白