如何拟合一个封闭的曲面

cichishia

请教大家：
已知若干个离散点的三维坐标，如何拟合成一个封闭的曲面。

  a. W1 b# M5 g# V程序：
%--------------------------------------------------------------------------
A=xlsread('PointData.xls');
6 \. [+ P* y% k% 下部分网格曲面
x=A(1:85,1);   % 根据dz的大小，相应调整行的范围--及确定好哪些点属于下部分
y=A(1:85,2);
z=A(1:85,3);
. F2 L  }: L/ u8 w* Uminx = min(x);
maxx = max(x);
miny = min(y);
1 _6 |$ s9 [* W3 a' Xmaxy = max(y);
xtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
ytemp=linspace(min(y),max(y),50);
[X,Y]=meshgrid(xtemp,ytemp);
# r% c1 A7 C$ b/ L& \* {+ UZ=griddata(x,y,z,X,Y,'cubic');
figure('numbertitle','off','name','网格图')
suRF(X,Y,Z)
shading faceted

%--------------------------------------------------------------------------
% 上部分网格曲面
. T& |8 y, V  U/ f7 v* chold on
; F& X6 L  L1 r) `( wx=A(61:404,1);
y=A(61:404,2);
8 Q9 t- s9 C/ Ez=A(61:404,3);
( {$ ~" @" o8 A9 z+ x% n5 M+ ~; Nminx = min(x);
maxx = 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);
Z=griddata(x,y,z,X,Y,'cubic');
surf(X,Y,Z)
) |; q  I, }+ ~! sshading faceted
' j! s9 p9 b; \9 F: G) W) t
figure('name','散点图')
) U+ D+ k9 b! E4 f9 M( n: o  rplot3(A(:,1),A(:,2),A(:,3),'*')
: c. P6 t5 q% v, C9 P

我自己用上述方法分别拟合上下两部分，但是图中有缺口，如何才能得到封闭的图形呢？

- ]1 A# T- _" ~2 D- s

IBB-EUT

帮你顶一下

shuddkk

给出数据，我可以试试。

小白的白