| 查看: 872 | 回复: 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 ] |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
dnp
荣誉版主 (知名作家)
小木虫浪子
- 应助: 74 (初中生)
- 贵宾: 18.476
- 金币: 18799.9
- 散金: 2040
- 红花: 77
- 沙发: 6
- 帖子: 7304
- 在线: 593.8小时
- 虫号: 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