| 查看: 1737 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
用DL_POLY模拟了硅酸盐熔体后，如何计算硅氧四面体Qn结构？
|
|||
| 模拟了CaO,Al2O3，SiO2的体系，然后想计算其中硅氧四面体类型分布，即计算与所有硅配位的氧原子是否桥氧， |
回复此楼» 收录本帖的淘帖专辑推荐
 动力学模拟 |
» 猜你喜欢
垃圾破二本职称评审标准
已经有19人回复
-
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
9楼2016-04-27 21:45:14
2楼2015-09-14 15:12:46
15216105346
金虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 2062.4
- 散金: 2693
- 红花: 7
- 帖子: 1205
- 在线: 298小时
- 虫号: 3441372
- 注册: 2014-09-25
- 专业: 固体无机化学
|
这是统计程序，matlab的！ 供你参考一下！ clear Ob=0;Onb=0;Of=0; A=xlsread('CS1'); [m,n]=size(A); Si=0;Ca=0;O=0; for i=2:m if A(i,4)==1 Si=Si+1; elseif A(i,4)==2 Ca=Ca+1; else O=O+1; end end for i=Si+Ca+2:Si+Ca+O+1 %O的个数 for j=2:Si+Ca+1 %Si和Ca的个数 B(i-Si-Ca-1,j-1)=sqrt((A(i,1)-A(j,1))^2+(A(i,2)-A(j,2))^2+(A(i,3)-A(j,3))^2); end [P,Q]=sort(B(i-Si-Ca-1,1:Si+Ca)); if Q(1)<Si+2&&Q(2)<Si+2 Ob=Ob+1; end if Q(1)<Si+2&&Q(2)>Si+1 Onb=Onb+1; end if Q(1)>Si+1&&Q(2)<Si+2 Onb=Onb+1; end if Q(1)>Si+1&&Q(2)>Si+1 Of=Of+1; end end Ob=Ob/O*100 Onb=Onb/O*100 Of=Of/O*100 发自小木虫IOS客户端 |
3楼2015-09-14 23:37:39
4楼2016-04-03 13:23:13