如何提取6 neasrest neigbors粒子吗

2视图(30天)

显示旧的评论

阿斯利瓦斯塔瓦 2021年5月11日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/827505-how-to-extract-6-neasrest-neigbors-for-particles

评论道: 阿斯利瓦斯塔瓦2021年6月3日

答:接受为副总经理

particle-distance.txt

大家好。

我有一个文本文件(请参见附加的文本文件)包含的粒子邻居(从1到64)。我想提取6最近的邻居的每个粒子在一个独立的文件。

下面是示例的输出文件。列1和图2显示粒子数,第三列显示的距离。

1 19 0.874004

1 41 1.130905

1 42 1.131306

1 45 0.841402

1 50 1.393870

1 52 0.961277

1 61 0.802248

1 63 1.327221

1 64 1.408026

2 16 0.819803

2 33 1.153725

2 37 0.954796

2 39 1.439849

2 42 1.419793

2 58 1.140030

2 60 0.912613

2 63年1.454144

3 6 0.843783

3 7 1.071621

3 9 1.009899

3 10 1.254684

3 13 1.108856

3 14 1.015914

3 20 0.885181

3 35 1.263609

3 45 1.356776

3 51 1.268250

4 5 1.197906

4四十1.238262

4 49 1.131004

4 53 1.207126

4 54 0.848770

4 57 1.278407

4 59 1.398415

4 62 0.782374

4 64 0.872153

我想要的是输出

1 19 0.874004

1 41 1.130905

1 42 1.131306

1 45 0.841402

1 50 1.393870

1 52 0.961277

2 16 0.819803

2 33 1.153725

2 37 0.954796

2 39 1.439849

2 42 1.419793

2 58 1.140030

3 6 0.843783

3 7 1.071621

3 9 1.009899

3 10 1.254684

3 13 1.108856

3 14 1.015914

4 5 1.197906

4四十1.238262

4 49 1.131004

4 53 1.207126

4 54 0.848770

4 57 1.278407

任何帮助将不胜感激。

提前谢谢,

阿

2的评论
显示隐藏 1年长的评论

阿斯利瓦斯塔瓦 2021年5月11日

下面是我的代码:

                                %代码6测定工业最近的邻居
                               
                                清晰的
                               
                                clc
                               
                                一个=负载(“particle-distance.txt”);
                               
                                b =一个;
                               
                                数= 0;
                               
                                natoms = 64;
                               
                                为我= 1:natoms
                               
                                为j = 1:长度(a)
                               
                                如果(j, 1) = =我
                               
                                数=计数+ 1;
                               
                                如果数> = 6
                               
                                str = [int2str ((j, 1)),' ',int2str ((j, 2)),' 'num2str ((j, 3)));
                               
                                disp (str)
                               
                                结束
                               
                                结束
                               
                                结束
                               
                                结束

登录置评。

在回答这个问题。

接受的答案

为副总经理 2021年5月11日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/827505-how-to-extract-6-neasrest-neigbors-for-particles answer_697375

编辑:为副总经理 2021年5月11日

这可能是改善,但这是我所做的:

                              A = [1 19 0.874004;
                             
                              1 41 1.130905;
                             
                              1 42 1.131306;
                             
                              1 45 0.841402;
                             
                              1 50 1.393870;
                             
                              1 52 0.961277;
                             
                              61 0.802248;
                             
                              63 1.327221;
                             
                              64 1.408026;
                             
                              2 16 0.819803;
                             
                              2 33 1.153725;
                             
                              2 37 0.954796;
                             
                              2 39 1.439849;
                             
                              2 42 1.419793;
                             
                              2 58 1.140030;
                             
                              2 60 0.912613;
                             
                              63 1.454144;
                             
                              3 6 0.843783;
                             
                              3 7 1.071621;
                             
                              3 9 1.009899;
                             
                              3 10 1.254684;
                             
                              3 13 1.108856;
                             
                              3 14 1.015914;
                             
                              3 20 0.885181;
                             
                              3 35 1.263609;
                             
                              3 45 1.356776;
                             
                              3 51 1.268250;
                             
                              4 5 1.197906;
                             
                              4个40 1.238262;
                             
                              4 49 1.131004;
                             
                              4 53 1.207126;
                             
                              4 54 0.848770;
                             
                              4 57 1.278407;
                             
                              4 59 1.398415;
                             
                              62 0.782374;
                             
                              4 64 0.872153)
                             
                              %按距离排序
                             
                              [~,深度贴图]=排序((:,3),“提升”);
                             
                              =(深度贴图:);
                             
                              npoints =元素个数(独特的((:1)));
                             
                              idx = [];
                             
                              为p = 1: npoints
                             
                              %这将找到最多6场比赛
                             
                              idx = [idx找到((:1)= = p, 6)];
                             
                              结束
                             
                              一个= (idx:)

给了

                              一个=
                             
                              1 61 0.80225
                             
                              1 45 0.8414
                             
                              1 19 0.874
                             
                              1 52 0.96128
                             
                              1 41 1.1309
                             
                              1 42 1.1313
                             
                              2 16 0.8198
                             
                              2 60 0.91261
                             
                              2 37 0.9548
                             
                              2 58 1.14
                             
                              2 33 1.1537
                             
                              2 42 1.4198
                             
                              3 6 0.84378
                             
                              3 20 0.88518
                             
                              3 9 1.0099
                             
                              3 14 1.0159
                             
                              3 7 1.0716
                             
                              3 13 1.1089
                             
                              4 62 0.78237
                             
                              4 54 0.84877
                             
                              4 64 0.87215
                             
                              4 49 1.131
                             
                              4 5 1.1979
                             
                              4 53 1.2071

12个评论
显示隐藏 11年长的评论

里克 2021年5月12日

我不假装自己是个专家,但我知道,使用我建议这样的编码风格是Mathworks工程师更容易优化比特定于JIT代码,利用行为在一个特定的版本,或CPU、或操作系统,等等。所以通过使用下面这样的代码,我可以合理地预期性能改进的版本比我。

我不知道谁最严格的人将会在这儿。我不认为我是。我认为这个论坛是一个挑战来解决难题,以及帮助和教育人们。尤其是最后一部分在这种情况下是很重要的。如果你教好的编码实践当没关系,人们更有可能去应用它。既然你似乎在坚持做,我只是想确保你正在一个深思熟虑的选择。

为副总经理 2021年5月12日

编辑:为副总经理 2021年5月12日

我想我从来没有真正想过它。我可能是优化自己的非常利基案例(旧版本在一个特定的操作系统硬件),不知道是典型的性能。我知道我的很多优化MIMT只有有用的旧版本(R2009b)和99.9%的用户可能是无用的。

关于我的目的在这里,我才开始回答问题时通过时间等待我自己的答案。现在只是变成一个懒惰的习惯。我承认大部分只是娱乐和游戏对我来说,虽然这并不是全部。我一直惊讶很多事情我从未有理由甚至尝试之前,它是很高兴帮助人们(当事实证明我可以)。我希望我的意见我自己的习惯将改变和塑造他们的前进。至少我写(略)更好的格式化代码,因为现在人们阅读。我怀疑是我珍惜最后的坏习惯遭受侵蚀。

里克 2021年5月12日

自从我发现,R13 (v6.5)是非常快的在某些任务(由于1月),我确定我的代码可以运行。我已经确定我的代码运行在八度所以我将能够分享我写了以法律的方式。

因为我已经覆盖范围广泛,它不是很多额外的努力来检查我的代码是如何工作的不同版本。随意查看我FEX提交的文档,如果你想看一看,困扰了到目前为止,)

有些人在这个论坛我之前一直在使用Matlab完成小学教育。如果你愿意,你可以学到很多。

关于别人阅读你的代码:你应该对待你未来的自己是一个不同的人。你可能不会理解你写的注释代码2年前。你至少应该写评论这。

阿斯利瓦斯塔瓦 2021年5月19日

6 - nn -框架- 1. - zip

你好,

我已经成功地使用您的代码生成的六个邻居。无论在我的系统中,有一些粒子没有确切的六个neigbors。所以我想eleminate那些台词文本文件。

的例子,

1 9 1.484157

1 16 0.828089

1 21 1.481878

2 14 1.084812

2 19 0.212855

2 23日0.745430

2 27日0.709192

2 29日0.485757

2 36 1.351770

3 8 1.485509

3 10 0.789934

3 12 0.996034

3 22 1.273936

3 28日1.240151

3 31日1.209610

4 1 0.445649

4 5 0.491061

4 6 0.611696

4 9 1.460465

4 16 1.081862

4 24 1.222578

5 1 0.868254

5 4 0.491061

5 6 0.366953

5 9 1.269292

5 16 1.212136

5 18 1.255205

6 1 1.052711

6 4 0.611696

6 5 0.366953

6 18 1.467878

6 20 1.142525

6 34 0.708171

7 8 0.398548

7 10 1.086699

7 14 1.288270

7 32 0.816957

7 37 0.881323

8 3 1.485509

8 7 0.398548

8 10 1.228414

8 32 0.781834

8 37 0.906441

9 1 1.484157

9 4 1.460465

9 5 1.269292

9 16 0.938138

9 18 0.803028

9 24 1.018272

粒子数7,8(第一列)5的邻居。

我想要有一个数据文件包含6个邻居。如果你能帮助,那将是非常好的。

提前谢谢,

阿

为副总经理 2021年6月3日

编辑:为副总经理 2021年6月3日

很抱歉。通知本网站往往出现的顺序,所以有时候我完全错过后续问题如果他们不出现,直到队列已经过去几个小时的时间戳。

例如:

                                   A = [1 19 0.874004;
                                  
                                   1 41 1.130905;
                                  
                                   1 50 1.393870;
                                  
                                   63 1.327221;
                                  
                                   64 1.408026;
                                  
                                   2 16 0.819803;
                                  
                                   2 33 1.153725;
                                  
                                   2 37 0.954796;
                                  
                                   2 39 1.439849;
                                  
                                   2 42 1.419793;
                                  
                                   2 58 1.140030;
                                  
                                   2 60 0.912613;
                                  
                                   63 1.454144;
                                  
                                   3 6 0.843783;
                                  
                                   3 7 1.071621;
                                  
                                   3 13 1.108856;
                                  
                                   3 35 1.263609;
                                  
                                   3 45 1.356776;
                                  
                                   4 5 1.197906;
                                  
                                   4个40 1.238262;
                                  
                                   4 49 1.131004;
                                  
                                   4 53 1.207126;
                                  
                                   4 54 0.848770;
                                  
                                   4 57 1.278407;
                                  
                                   4 59 1.398415;
                                  
                                   62 0.782374;
                                  
                                   4 64 0.872153)
                                  
                                   n = 6;%最小数量的邻居
                                  
                                   [~,深度贴图]=排序((:,3),“提升”);
                                  
                                   =(深度贴图:);
                                  
                                   npoints =元素个数(独特的((:1)));
                                  
                                   idx = [];
                                  
                                   为p = 1: npoints
                                  
                                   %这最多能找到n匹配
                                  
                                   %如果有少于n邻居,跳过它
                                  
                                   thisblock =找到((:1)= = p, n);
                                  
                                   如果元素个数(thisblock) > = n
                                  
                                   idx = [idx thisblock];
                                  
                                   结束
                                  
                                   结束
                                  
                                   一个= (idx:)

阿斯利瓦斯塔瓦 2021年6月3日

上面的代码工作。

非常感谢你的帮助。

登录置评。

类别

游戏只是为了好玩

找到更多的在只是为了好玩在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

如何提取6 neasrest neigbors粒子吗

2的评论
显示隐藏 1年长的评论

接受的答案

12个评论
显示隐藏 11年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

如何提取6 neasrest neigbors粒子吗

2的评论 显示隐藏 1年长的评论

接受的答案

12个评论 显示隐藏 11年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

2的评论
显示隐藏 1年长的评论

12个评论
显示隐藏 11年长的评论