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

首先，我先描述一下我的问题，我得到了一个三维数组，数字代表了三维空间中的势能，我想把等势面画出来，也就是把势能相同的点连接起来形成面，问下有没有大牛会的，之前我用的contour可以画二维的，三维的真心不会了，我附上二维的程序还有球三维数组的程序，希望有好心人帮我，十分感谢！

function [U1,k]=cff(a)%二维等势面，可以实现%

x=5;y=5;

U1=[0 100 100 100 0

0 75 75 75 0

0 50 50 50 0

0 25 25 25 0

0 0 0 0 0];

U2=U1;maxt=1;k=0;

while(~(maxt<1e-3||k>100000))

k=k+1;

maxt=0;

for i=2:y-1,

for j=2:x-1,

U2(i,j)=U1(i,j)+a*(U1(i,j+1)+U1(i+1,j)+U2(i-1,j)+U2(i,j-1)-4*U1(i,j))/4;

t=abs(U2(i,j)-U1(i,j));

if(t>maxt) maxt=t;

end

end

end

U1=U2;

contour(U1);

grid on

end

end

function cff1(a)%三维的contour函数实现不了

x=5;y=5;z=5;

U1=[0 100 100 100 0

0 75 75 75 0

0 50 50 50 0

0 25 25 25 0

0 0 0 0 0];

C1=repmat(U1,[1,1,5]);

C2=C1;maxt=1;k=0;

while(~(maxt<1e-3||k>100000))

k=k+1;

maxt=0;

for i=2:y-1,

for j=2:x-1,

for m=2:z-1,

C2(i,j,m)=C1(i,j,m)+a*(C1(i,j+1,m)+C1(i+1,j,m)+C2(i-1,j,m)+C2(i,j-1,m)+C1(i,j,m+1)+C2(i,j,m-1)-6*C1(i,j,m))/6;

t=abs(C2(i,j,m)-C1(i,j,m));

if(t>maxt) maxt=t;

end

end

end

end

C1=C2

contour(C1);

grid on

end

end