24小时热门版块排行榜

>论坛更新日志 (4226)
>虫友互识 (436)
>文献求助 (400)
>导师招生 (389)
>论文投稿 (135)
>硕博家园 (128)
>考博 (113)
>休闲灌水 (110)
>博后之家 (100)
>招聘信息布告栏 (98)
>催化 (59)
>公派出国 (56)
>基金申请 (51)
>考研 (51)
>论文道贺祈福 (46)
>教师之家 (45)

返回列表

[交流] 给一个画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;

回复此楼

» 猜你喜欢

2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有4人回复
磺酰氟产物，毕不了业了！已经有6人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有9人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
论文终于录用啦！满足毕业条件了已经有25人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

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

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的角度的数据，可以建立一下。。。