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的回帖

sgyin

铁杆木虫 (正式写手)

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

谢谢y1ding，没看懂r(ii_ii)=sqrt((A(jj,：)-A(ii,：)).^2);里面r(ii_ii)是什么意思啊？运行也有错误：
??? Undefined function or variable 'ii_ii'.

Error in ==> D:\Program Files\MATLAB6.5\bin\win32\temp.m
On line 23  ==>    r(ii_ii)=sqrt((A(jj,：)-A(ii,：)).^2);

今天我也试了一下，
for iii=1:m,
for ii=1+iii:m,
      rr((iii-1)*m-iii*(iii-1)/2+ii-iii,：)=(A(iii,：)-A(ii,：)).^2;
      r((iii-1)*m-iii*(iii-1)/2+ii-iii)=(rr((iii-1)*m-iii*(iii-1)/2+ii-iii,1)^2+rr((iii-1)*m-iii*(iii-1)/2+ii-iii,2)^2+rr((iii-1)*m-iii*(iii-1)/2+ii-iii,3)^2)^(.5);
   end
end
这样是可以的，可是rr这个矩阵比较大，浪费内存啊

[ Last edited by sgyin on 2008-1-2 at 10:39 ]