| 查看: 1018 | 回复: 9 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
sgyin铁杆木虫 (正式写手)
|
[交流]
【求助】Matlab如何编写，请指教
|
||
|
文件名：ZnO.xyz(文本可打开)，存储10个原子坐标。 内容如下： ============ 10 Zn -0.00133 1.53941 0.00000 Zn -1.33383 0.76855 2.47350 Zn 1.33117 0.76855 2.47350 Zn -0.00133 3.07651 2.47350 Zn 0.00000 1.53710 0.00000 O -1.33117 0.76855 2.47350 O 1.33383 0.76855 2.47350 O 0.00133 3.07651 2.47350 O 0.00133 1.53941 0.00000 O -1.33250 0.77086 2.47350 ============= 已知： f_Zn=1673; f_O =2432; 请问如何求：（1）r_nm （两两原子之间的距离，n、m=Zn或O）； （2）对n、m求和{ f_n*f_m*(1+r_nm) } [ Last edited by csfn on 2008-12-29 at 20:18 ] |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
dnp
荣誉版主 (知名作家)
小木虫浪子
- 应助: 74 (初中生)
- 贵宾: 18.476
- 金币: 18827.9
- 散金: 2040
- 红花: 77
- 沙发: 6
- 帖子: 7319
- 在线: 598.4小时
- 虫号: 437184
- 注册: 2007-08-26
- 性别: GG
- 专业: 生物信息学
- 管辖: 有机交流
|
不知道你的那个f_Zn是不是已知的，如果是已知的话，你在编程的时候可以当作输入参数； 要是你直接在命令窗口计算，可以手动输入； 如果需要计算的话，当你将数据导入的时候可以利用导入的数据进行计算再结合你的计算公式进行计算，这个也很容易实现。 说实话，我看不懂你的第二个算式是怎么回事，求和？后面的算式好像不像是一个求和算法。 你在用fgetl读取数据之后，你可以进行操作的，不要把他们都放到一个字符串里面，也就是说，把前面代表原子信息的放到一个字符串向量(char格式)或者是单元数组(cell)里面，而剩余的坐标数据就存为一个三列的矩阵，这样，你就可以随意地调用坐标数据来进行你的计算了。这里我想建议你用一个命令，就是判断空格的命令：isspace (详细看help)。这样你读取数据的时候就方便多了。 祝楼主好运！ |
10楼2008-01-18 09:23:13
davin
木虫 (著名写手)
一只来自东北的大虫
- 应助: 6 (幼儿园)
- 贵宾: 0.549
- 金币: 4293.9
- 散金: 676
- 红花: 59
- 帖子: 1341
- 在线: 718小时
- 虫号: 420443
- 注册: 2007-07-10
- 性别: GG
- 专业: 传热传质学
,首先，你最好将Zn和O删除。
;2楼2008-01-14 09:03:06
sgyin
铁杆木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5519.4
- 散金: 100
- 红花: 8
- 帖子: 680
- 在线: 279.4小时
- 虫号: 405447
- 注册: 2007-06-17
- 性别: GG
- 专业: 材料物理化学
3楼2008-01-14 11:30:45
|
rn_m= 0 2.9134 2.9134 2.9122 0.0026655 2.9122 2.9146 2.9122 0.00266 2.9122 2.9134 0 2.665 2.665 2.9134 0.00266 2.6677 2.6663 2.9146 0.0026655 2.9134 2.665 0 2.665 2.9122 2.6623 0.00266 2.6637 2.9122 2.6637 2.9122 2.665 2.665 0 2.9134 2.6637 2.6663 0.00266 2.9122 2.6623 0.0026655 2.9134 2.9122 2.9134 0 2.9122 2.9134 2.9134 0.0026655 2.9122 2.9122 0.00266 2.6623 2.6637 2.9122 0 2.665 2.665 2.9134 0.0026655 2.9146 2.6677 0.00266 2.6663 2.9134 2.665 0 2.665 2.9134 2.6663 2.9122 2.6663 2.6637 0.00266 2.9134 2.665 2.665 0 2.9122 2.6637 0.00266 2.9146 2.9122 2.9122 0.0026655 2.9134 2.9134 2.9122 0 2.9134 2.9122 0.0026655 2.6637 2.6623 2.9122 0.0026655 2.6663 2.6637 2.9134 0 |
4楼2008-01-14 15:17:29
100