24小时热门版块排行榜

返回列表

[交流] 给一个画C60图的MATLAB程序，但是数据文件不知道

CODE:

function C60

N=60;R=3.517;

load C60.mat theta phi;

x=R*sin(theta).*cos(phi);

y=R*sin(theta).*sin(phi);

z=R*cos(theta);

X=repmat(x,N,1);Y=repmat(y,N,1);

Z=repmat(z,N,1);

D=(X-X').^2+(Y-Y')^2+(Z-Z')^2;

M=sqrt(D);

bond=(M<2.1);

M=sort(M,1);

M=M.*(M<=2.1);

CNT=(0.01
M=M.*(M>1.8);

cnt=(1.8
CNT=CNT+cnt;

corN=sum(CNT,1);

total=sum(corN,2)/N;

drawAtoms(x,y,z,N,bond);

function drawAtoms(x,y,z,N,bond)

figure;

plot3(x,y,z,'.g','MarkerSize',36)

title('Fullerene Nucleation');

xlabel('X');ylabel('Y');zlabel('Z');

grid on;axis equal;

for i=1:N

    for j=i+1:N

        xtmp=[x(i) x(j)];ytmp=[y(i) y(j)];ztmp=[z(i) z(j)];

        if(bond(i,j)~=0)h=line(xtmp,ytmp,ztmp);

            Set(h,'color','b');end;

    end;

end;

回复此楼

» 猜你喜欢

体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有17人回复
面上可以超过30页吧？已经有7人回复
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？已经有5人回复
“人文社科而论，许多学术研究还没有达到民国时期的水平” 已经有6人回复
版面费该交吗已经有13人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有10人回复
什么是人一生最重要的？已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab如何读取一个文件夹下的所有图片已经有9人回复
帮忙看看这个平面图用matlab怎么画已经有19人回复
MATLAB 画图小问题已经有11人回复
如何用matlab处理txt文件里的数据已经有7人回复
求助：谁会用matlab画图？已经有13人回复
用Matlab软件如何画出三维荧光光谱图已经有11人回复
【求助】matlab 中，几个m文件调用已经有14人回复
【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧已经有105人回复
Matlab 画图后怎么获得图上的数据点已经有6人回复
【求助】请教一个matlab程序，谢谢帮忙！！已经有15人回复
【求助】求助帮助编一个读取TEXT文件画图的程序已经有8人回复
【求助】matlab中如何将程序的结果输出到一个文件中已经有11人回复
用matlab画的图保存什么格式插入word里比较清楚？已经有45人回复
【求助】MATLAB多个mat中数据画图问题已经有7人回复
【求助】柱形统计图matlab怎么画？已经有4人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-11-08 17:07:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

1楼: Originally posted by zyj8119 at 2011-11-08 17:07:22:
[code]function C60
N=60;R=3.517;
load C60.mat theta phi;
x=R*sin(theta).*cos(phi);
y=R*sin(theta).*sin(phi);
z=R*cos(theta);
X=repmat(x,N,1);Y=repmat(y,N,1);
Z=repmat(z,N,1);
D=(X-X').^2+(Y ...

如果谁知道C60的角度的数据，可以建立一下。。。