没啥好说的,上代码
A=importdata('E:\各种数据集\220417\y001.txt');
B=importdata('E:\各种数据集\220417\001yy.txt');
% 重新定义或者不定义都行
Ax=A(:, 1);
Ay=A(:, 2);
Az=A(:, 3);
Bx=B(:, 1);
By=B(:, 2);
Bz=B(:, 3);
view(3) % 设置默认三维视图
% convhulln函数就是算凸包的,得到的是凸包的坐标们
f = convhulln(A);
% patch函数,将坐标点连成面,形成包络体
patch('vertices',A,'faces',f,'facecolor','r')
axis equal
f1 = convhulln(B);
patch('vertices',B,'faces',f1,'facecolor','r')
axis equal
% 这里开始是算包络体的体积的,k是凸包点坐标,就跟上面的f一样,是个矩阵,av就是该包络体的体积
[k1,av1] = convhull(Ax,Ay,Az);
[k2,av2] = convhull(Bx,By,Bz);
av1
av2
我画了两组数据的,所以,效果图:
每天学一点点,就能少菜一点点,加油( •̀ ω •́ )y