| 查看: 1277 | 回复: 9 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
sgyin铁杆木虫 (正式写手)
|
[交流]
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 ] |
回复此楼» 猜你喜欢
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
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 ] |
6楼2008-01-02 10:30:29
y1ding
铁杆木虫 (著名写手)
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
2楼2008-01-02 09:02:28
y1ding
铁杆木虫 (著名写手)
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
3楼2008-01-02 09:02:42
y1ding
铁杆木虫 (著名写手)
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
4楼2008-01-02 09:02:57
30