EDA365电子论坛网

标题: 如何拟合一个封闭的曲面 [打印本页]

---

作者: cichishia    时间: 2021-1-27 13:22
标题: 如何拟合一个封闭的曲面
请教大家：
9 I; |, Z) Y; F# M8 ^: C* X+ W已知若干个离散点的三维坐标，如何拟合成一个封闭的曲面。
- c. [, I+ d. O  s4 t
程序：
%--------------------------------------------------------------------------
A=xlsread('PointData.xls');
' F) L$ B9 ?8 D% 下部分网格曲面
x=A(1:85,1);   % 根据dz的大小，相应调整行的范围--及确定好哪些点属于下部分
y=A(1:85,2);
8 H/ G# O( h. U. @4 zz=A(1:85,3);
minx = min(x);
3 @% N6 U. L$ \8 r0 `maxx = max(x);
miny = min(y);
- M1 G" B3 u& g* C+ i, Y4 n$ bmaxy = max(y);
xtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
ytemp=linspace(min(y),max(y),50);
* H' @3 k2 m! D4 n8 F; e  x[X,Y]=meshgrid(xtemp,ytemp);
Z=griddata(x,y,z,X,Y,'cubic');
figure('numbertitle','off','name','网格图')
surf(X,Y,Z)
shading faceted
* |: n! B' _6 s7 K
%--------------------------------------------------------------------------
% 上部分网格曲面
hold on
' f* q" H/ }8 N# v$ }" fx=A(61:404,1);
' _) a) y, C! @: H1 Jy=A(61:404,2);
1 D+ b% L9 Y+ lz=A(61:404,3);
minx = min(x);
maxx = max(x);
! B: l: v2 p2 w3 Yminy = min(y);
maxy = max(y);
$ E$ l+ L( y* ?# P& F! y; Xxtemp=linspace(min(x),max(x),50);%这个是选择的数值样本大小，越大曲线越平滑%
3 g* P0 U( B. cytemp=linspace(min(y),max(y),50);
$ {" k0 K# A- S, N$ O$ w[X,Y]=meshgrid(xtemp,ytemp);
7 K8 r0 p0 M: L5 NZ=griddata(x,y,z,X,Y,'cubic');
: K; t$ O$ k1 O  vsurf(X,Y,Z)
5 d# V4 O, J4 y+ V0 S7 Ushading faceted
& L1 I0 [$ n0 o9 d" j
! Q( z+ Z0 k  r* v" a( ~$ ofigure('name','散点图')
plot3(A(:,1),A(:,2),A(:,3),'*')

! R% Q2 s/ O/ l  F
1 g# |+ y8 M* s" }, e4 j/ j我自己用上述方法分别拟合上下两部分，但是图中有缺口，如何才能得到封闭的图形呢？

) Q# p% E3 R0 Y

---

作者: IBB-EUT    时间: 2021-1-27 13:32
帮你顶一下

---

作者: shuddkk    时间: 2021-1-27 16:59
给出数据，我可以试试。

---

作者: 小白的白    时间: 2021-1-27 17:50

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2