| 查看: 2131 | 回复: 8 | |||
[求助]
用DL_POLY模拟了硅酸盐熔体后，如何计算硅氧四面体Qn结构？
|
| 模拟了CaO,Al2O3，SiO2的体系，然后想计算其中硅氧四面体类型分布，即计算与所有硅配位的氧原子是否桥氧， |
回复此楼» 收录本帖的淘帖专辑推荐
 动力学模拟 |
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有23人回复
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
15216105346
金虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 2062.4
- 散金: 2693
- 红花: 7
- 帖子: 1205
- 在线: 298小时
- 虫号: 3441372
- 注册: 2014-09-25
- 专业: 固体无机化学
5楼2016-04-03 13:47:49
6楼2016-04-07 14:57:13
15216105346
金虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 2062.4
- 散金: 2693
- 红花: 7
- 帖子: 1205
- 在线: 298小时
- 虫号: 3441372
- 注册: 2014-09-25
- 专业: 固体无机化学
7楼2016-04-07 17:59:05
8楼2016-04-27 21:37:47
9楼2016-04-27 21:45:14
