空间几何体的直观图matlab

function B = transformation(A)
sqrt2=1.414;
n=size(A,1);
B=zeros(n,2);
B(:,1)=A*[-sqrt2/4,1,0]';
B(:,2)=A*[-sqrt2/4,0,1]';

bmin=min(min(B(:,1),B(:,2)));

if bmin < 0
    B = B - bmin;
end


p=max(max(B)-min(B));
k=300/p;
B=k*B;
B(:,2)=300-B(:,2);


B=B+10;

sprintf('%.0f ',B')

zimu = ["a","b","c","d", "e" , "f" "g" "h"];
calptal = ["A","B","C","D", "E" , "F" "G" "H"];

for m =1:n
fprintf('circle.move( name="%s"', zimu(m));
fprintf(' cx=%3.0f',B(m,1));
fprintf(' cy=%3.0f',B(m,2));
fprintf(' label="%s")
', calptal(m))
end

end
%A=[0,0,0; 0,4,0;-4,4,0;-4,0,0; 0,0,4; 0,4,4;-4,4,4;-4,0,4]
%
% 2020-2-29 A=[A;A';B;B';C;Cp;D;Dp];
% A=[18,0,0;18,23,0;0,23,0;0,0,0;18,0,23;18,23,23;0,23,23;0,0,23];
% B = transformation(A);

  320*320的画布；10边框宽度

>> A=[18,0,0;18,23,0;0,23,0;0,0,0;18,0,23;18,23,23;0,23,23;0,0,23];
>> B = transformation(A);

ans =

    '10 310 245 310 310 245 75 245 10 75 245 75 310 10 75 10 '

circle.move( name="a" cx= 10 cy=310 label="A")
circle.move( name="b" cx=245 cy=310 label="B")
circle.move( name="c" cx=310 cy=245 label="C")
circle.move( name="d" cx= 75 cy=245 label="D")
circle.move( name="e" cx= 10 cy= 75 label="E")
circle.move( name="f" cx=245 cy= 75 label="F")
circle.move( name="g" cx=310 cy= 10 label="G")
circle.move( name="h" cx= 75 cy= 10 label="H")