![](https://iknow-base.cdn.bcebos.com/lxb/notice.png)
1个回答
展开全部
绘制立方网格matlab代码如下
m=3;n=3;h=3;
%定义数据的规模
Hx=1;Hy=1;Hz=1;
H=(h-1)*m*n;
CN=m*n*h;
N=m*n;
[x,y,z]=meshgrid(0:Hx:(n-1)*Hx,0:Hy:(m-1)*Hy,0:Hz:(h-1)*Hz);
Cube=reshape(1:N*h,m,n,h);
MN=2*m*n-m-n;
X=zeros(CN,1);Y=X;Z=X;
for i=1:m
for j=1:n
for k=1:h
X(Cube(i,j,k))=x(i,j,k);
Y(Cube(i,j,k))=y(i,j,k);
Z(Cube(i,j,k))=z(i,j,k);
end
end
end
linkN=0;
Img=sqrt(-1);
AA=zeros(H+h*MN,1);
for k=1:h-1
for j=1:n
for i=1:m
linkN=linkN+1;
AA(linkN)=Cube(i,j,k)+Img*Cube(i,j,k+1);
end
end
end
for k=1:h
for j=1:n-1
for i=1:m
linkN=linkN+1;
AA(linkN)=Cube(i,j,k)+Img*Cube(i,j+1,k);
end
end
for j=1:n
for i=1:m-1
linkN=linkN+1;
AA(linkN)=Cube(i,j,k)+Img*Cube(i+1,j,k);
end
end
end
figure;
plot3(X,Y,Z,'r.')
%--plot bond-----
axis equal
P1=real(AA);P2=imag(AA);
line([X(P1)';X(P2)'],[Y(P1)';Y(P2)'],[Z(P1)';Z(P2)'],'color','b')
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询