24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sgyin

铁杆木虫 (正式写手)

应助: 13 (小学生)
金币: 5519.4
散金: 100
红花: 8
帖子: 680
在线: 279.4小时
虫号: 405447
注册: 2007-06-17
性别: GG
专业: 材料物理化学

[交流] matlab编程问题

A是m*3的矩阵，是m个原子的坐标值，请问如何求的所有两两原子之间的距离R?(R应该有m*(m-1)/2个分量。)
我试着用for

for ii=1:m,
r(ii)=sqrt((A(1,：)-A(ii,：)).^2);
end
这样只得到m个值。
该如何把这个程序完成啊？
谢谢

[ Last edited by csfn on 2008-12-29 at 20:20 ]

回复此楼

» 猜你喜欢

有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复

&amp;lt;ahref=&amp;quot;http://www.TickerFactory.com/&amp;quot;&amp;gt;&amp;lt;imgborder=&amp;quot;0&amp;quot;src=&am

1楼 2008-01-01 23:31:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dnp

荣誉版主 (知名作家)

小木虫浪子

应助: 74 (初中生)
贵宾: 18.476
金币: 18799.9
散金: 2040
红花: 77
沙发: 6
帖子: 7304
在线: 593.8小时
虫号: 437184
注册: 2007-08-26
性别: GG
专业: 生物信息学
管辖: 有机交流

★ ★ ★ ★
spur(金币+2,VIP+0):感谢帮助！欢迎常来！
csfn(金币+2,VIP+0):thanks~

我想你的矩阵A中的坐标值有三个吧？分别代表三维的坐标，你可以这样算：
R=[];
for i=1:m
for j=i:m
R(i,j)=sqrt((A(i,: )-A(j,: )).^2);
R(j,i)=R(i,j);
end
end
这样应该得到的是一个m*m的对称矩阵，对角线上的元素为0。
如果m<1000的话，这样算很快，而且matlab占的内存也不大。