| 查看: 1992 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
用DL_POLY模拟了硅酸盐熔体后，如何计算硅氧四面体Qn结构？
|
|||
| 模拟了CaO,Al2O3，SiO2的体系，然后想计算其中硅氧四面体类型分布，即计算与所有硅配位的氧原子是否桥氧， |
回复此楼» 收录本帖的淘帖专辑推荐
 动力学模拟 |
» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
-
是这周出结果还是下周出结果？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
filecode
已经有10人回复
6楼2016-04-07 14:57:13
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
