該樓層疑似違規(guī)已被系統(tǒng)折疊?隱藏此樓查看此樓

首先，我先描述一下我的問題，我得到了一個三維數(shù)組，數(shù)字代表了三維空間中的勢能，我想把等勢面畫出來，也就是把勢能相同的點(diǎn)連接起來形成面，問下有沒有大牛會的，之前我用的contour可以畫二維的，三維的真心不會了，我附上二維的程序還有球三維數(shù)組的程序，希望有好心人幫我，十分感謝！

function [U1,k]=cff(a)%二維等勢面，可以實(shí)現(xiàn)%

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函數(shù)實(shí)現(xiàn)不了

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

總結(jié)

以上是生活随笔為你收集整理的用matlab画三维凸起,求助大牛MATLAB画三维等势面的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。