如何拟合一个封闭的曲面

cichishia

请教大家：
已知若干个离散点的三维坐标，如何拟合成一个封闭的曲面。
- b9 i: N) w5 t+ H$ z6 A  k7 V
  R7 u( ^$ x$ g$ Z7 u程序：
%--------------------------------------------------------------------------
A=xlsread('PointData.xls');
5 a/ w! R( c, K% _% 下部分网格曲面
x=A(1:85,1);   % 根据dz的大小，相应调整行的范围--及确定好哪些点属于下部分
y=A(1:85,2);
5 k7 y: z: @+ P  ?, T4 K- fz=A(1:85,3);
2 u5 _- |6 c! i! H' [- Y3 Y! j1 rminx = min(x);
* }( V, b" ~* T# w& `8 W3 lmaxx = max(x);
miny = min(y);
maxy = max(y);
xtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
ytemp=linspace(min(y),max(y),50);
[X,Y]=meshgrid(xtemp,ytemp);
0 ~$ \& X! J5 W) \Z=griddata(x,y,z,X,Y,'cubic');
6 o5 B7 a: W# t5 Sfigure('numbertitle','off','name','网格图')
3 o; u, O. E; W4 M0 e8 p9 f% Q' x4 TsuRF(X,Y,Z)
# ~' p' l; `. r% K" g: k3 U& I7 }shading faceted
- I" Z2 J7 [$ U! d4 z, N- b
%--------------------------------------------------------------------------
% 上部分网格曲面
& z$ B6 [. H+ g9 ~6 e: Qhold on
x=A(61:404,1);
# e9 x2 ^9 x% I% E" y- w: Gy=A(61:404,2);
: \' }, m8 K7 O4 `z=A(61:404,3);
minx = min(x);
. v+ m" h  w# \2 W! W5 ymaxx = max(x);
: k8 Y1 v- a+ U5 P. T/ tminy = min(y);
maxy = max(y);
7 [% P$ {, p( E* ]8 M/ oxtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
3 h& E! k& B+ I4 }* j; c% \ytemp=linspace(min(y),max(y),50);
: Y1 J1 f# f' L" }! Y: p. R7 _[X,Y]=meshgrid(xtemp,ytemp);
Z=griddata(x,y,z,X,Y,'cubic');
$ v" f2 y5 r% V4 a3 q' Q- @surf(X,Y,Z)
shading faceted

figure('name','散点图')
plot3(A(:,1),A(:,2),A(:,3),'*')

% L4 F; t: l5 {. M: C& R" W/ s
; ], h+ D* |2 a$ `( Z  F3 c我自己用上述方法分别拟合上下两部分，但是图中有缺口，如何才能得到封闭的图形呢？

IBB-EUT

帮你顶一下

shuddkk

给出数据，我可以试试。

小白的白