matlab 画包络面,求大神指点绘制空间内散点图的包络面，，，散点程序如下

红太狼 2023-01-16 03:25 92阅读 0赞

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

E=\[\];a=0;b=0;c=0;

for Zp=-50:2:50

for Xp=-200:2:200

for Yp=-200:2:200

P1\_1=\[cos(b)\*cos(c) -cos(b)\*sin(c) sin(b) Xp;sin(a)\*sin(b)\*cos(c)+cos(a)\*sin(c) cos(a)\*cos(c)-sin(a)\*sin(b)\*sin(c) -sin(a)\*cos(b) Yp;sin(a)\*sin(c)-cos(a)\*sin(b)\*cos(c) cos(a)\*sin(b)\*sin(c)+sin(a)\*cos(c) cos(a)\*cos(b) Zp;0 0 0 1\]\*\[42.01;50.25;0;1\]+\[0;0; -208.897;0\];

P2\_1=\[cos(b)\*cos(c) -cos(b)\*sin(c) sin(b) Xp;sin(a)\*sin(b)\*cos(c)+cos(a)\*sin(c) cos(a)\*cos(c)-sin(a)\*sin(b)\*sin(c) -sin(a)\*cos(b) Yp;sin(a)\*sin(c)-cos(a)\*sin(b)\*cos(c) cos(a)\*sin(b)\*sin(c)+sin(a)\*cos(c) cos(a)\*cos(b) Zp;0 0 0 1\]\*\[-42.01;50.25;0;1\]+\[0;0; -208.897;0\];

P3\_1=\[cos(b)\*cos(c) -cos(b)\*sin(c) sin(b) Xp;sin(a)\*sin(b)\*cos(c)+cos(a)\*sin(c) cos(a)\*cos(c)-sin(a)\*sin(b)\*sin(c) -sin(a)\*cos(b) Yp;sin(a)\*sin(c)-cos(a)\*sin(b)\*cos(c) cos(a)\*sin(b)\*sin(c)+sin(a)\*cos(c) cos(a)\*cos(b) Zp;0 0 0 1\]\*\[-64.52;11.26;0;1\]+\[0;0; -208.897;0\];

P4\_1=\[cos(b)\*cos(c) -cos(b)\*sin(c) sin(b) Xp;sin(a)\*sin(b)\*cos(c)+cos(a)\*sin(c) cos(a)\*cos(c)-sin(a)\*sin(b)\*sin(c) -sin(a)\*cos(b) Yp;sin(a)\*sin(c)-cos(a)\*sin(b)\*cos(c) cos(a)\*sin(b)\*sin(c)+sin(a)\*cos(c) cos(a)\*cos(b) Zp;0 0 0 1\]\*\[-22.51;-61.51;0;1\]+\[0;0; -208.897;0\];

P5\_1=\[cos(b)\*cos(c) -cos(b)\*sin(c) sin(b) Xp;sin(a)\*sin(b)\*cos(c)+cos(a)\*sin(c) cos(a)\*cos(c)-sin(a)\*sin(b)\*sin(c) -sin(a)\*cos(b) Yp;sin(a)\*sin(c)-cos(a)\*sin(b)\*cos(c) cos(a)\*sin(b)\*sin(c)+sin(a)\*cos(c) cos(a)\*cos(b) Zp;0 0 0 1\]\*\[22.51;-61.51;0;1\]+\[0;0; -208.897;0\];

P6\_1=\[cos(b)\*cos(c) -cos(b)\*sin(c) sin(b) Xp;sin(a)\*sin(b)\*cos(c)+cos(a)\*sin(c) cos(a)\*cos(c)-sin(a)\*sin(b)\*sin(c) -sin(a)\*cos(b) Yp;sin(a)\*sin(c)-cos(a)\*sin(b)\*cos(c) cos(a)\*sin(b)\*sin(c)+sin(a)\*cos(c) cos(a)\*cos(b) Zp;0 0 0 1\]\*\[64.52;11.26;0;1\]+\[0;0; -208.897;0\];

ZB1\_1=P1\_1(3)+sqrt(223^2-(P1\_1(1)-27.15)^2-(P1\_1(2)-126.87)^2);

ZB2\_1=P2\_1(3)+sqrt(223^2-(P2\_1(1)+27.15)^2-(P2\_1(2)-126.87)^2);

ZB3\_1=P3\_1(3)+sqrt(223^2-(P3\_1(1)+123.45)^2-(P3\_1(2)+39.92)^2);

ZB4\_1=P4\_1(3)+sqrt(223^2-(P4\_1(1)+96.29)^2-(P4\_1(2)+86.95)^2);

ZB5\_1=P5\_1(3)+sqrt(223^2-(P5\_1(1)-96.29)^2-(P5\_1(2)+86.95)^2);

ZB6\_1=P6\_1(3)+sqrt(223^2-(P6\_1(1)-123.45)^2-(P6\_1(2)+39.92)^2);

if ((ZB1\_1<=50)&&(ZB2\_1<=50)&&(ZB3\_1<=50)&&(ZB4\_1<=50)&&(ZB5\_1<=50)&&(ZB6\_1<=50)&&(ZB1\_1>=(-50))&&(ZB2\_1>=(-50))&&(ZB3\_1>=(-50))&&(ZB4\_1>=(-50))&&(ZB5\_1>=(-50))&&(ZB6\_1>=(-50)))

P0\_1=\[cos(b)\*cos(c) -cos(b)\*sin(c) sin(b) Xp;sin(a)\*sin(b)\*cos(c)+cos(a)\*sin(c) cos(a)\*cos(c)-sin(a)\*sin(b)\*sin(c) -sin(a)\*cos(b) Yp;sin(a)\*sin(c)-cos(a)\*sin(b)\*cos(c) cos(a)\*sin(b)\*sin(c)+sin(a)\*cos(c) cos(a)\*cos(b) Zp;0 0 0 1\]\*\[0;0;0;1\]+\[0;0; -208.897;0\];

E=\[E,P0\_1\];

end

Xp0\_1=E(1,:);

Yp0\_1=E(2,:);

Zp0\_1=E(3,:);

plot3(Xp0\_1,Yp0\_1,Zp0\_1,'.');

xlabel('x');ylabel('y');zlabel('z');

grid on