| 查看: 1955 | 回复: 8 | |||
[求助]
用DL_POLY模拟了硅酸盐熔体后，如何计算硅氧四面体Qn结构？
|
| 模拟了CaO,Al2O3，SiO2的体系，然后想计算其中硅氧四面体类型分布，即计算与所有硅配位的氧原子是否桥氧， |
回复此楼» 收录本帖的淘帖专辑推荐
 动力学模拟 |
» 猜你喜欢
filecode变化情况
已经有16人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
-
filecode
已经有19人回复
-
基金中了
已经有11人回复
-
听说今天filecode变了
已经有44人回复
-
固定端突然变了，今天
已经有8人回复
-
filecode
已经有11人回复
-
求各位大神看下
已经有33人回复
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
100