我有一个类似的文件:
0.000000
0.000051
0.000076
0.000102
0.000152
0.000178
0.000203
0.000229
0.000254
0.000279
0.000305
0.000356
0.000381
0.000406
我想找到两个数字及其最近的距离。 我正在使用shell脚本。当时我有:
points.dat cat | cut-f1 | sort | uniq
在下一种情况下,我想考虑两列。那么x轴的列和y轴的列也是如此。 这就像最近的一对积分问题。 我也想帮忙。
感谢您的帮助!
答案 0 :(得分:2)
这是评论无法回答。
检查算法时使用此输出。
Python代码:
import itertools
data = [
0.000000, #data[0]
0.000051, #data[1]
0.000076, #data[2]
0.000102, #data[3]
0.000152, #data[4]
0.000178, #data[5]
0.000203, #data[6]
0.000229, #data[7]
0.000254, #data[8]
0.000279, #data[9]
0.000305, #data[10]
0.000356, #data[11]
0.000381, #data[12]
0.000406 #data[13]
]
for j in sorted((abs(i[0][0]-i[1][0]), (i[0][1],i[1][1])) for i in itertools.combinations( zip(data,range(len(data))), 2)):
print str(j[1]).ljust(8), '->',j[0]#indices -> distance
输出:
(7, 8) -> 2.5e-05
(1, 2) -> 2.5e-05
(5, 6) -> 2.5e-05
(8, 9) -> 2.5e-05
(11, 12) -> 2.5e-05
(12, 13) -> 2.5e-05
(4, 5) -> 2.6e-05
(9, 10) -> 2.6e-05
(2, 3) -> 2.6e-05
(6, 7) -> 2.6e-05
(7, 9) -> 5e-05
(3, 4) -> 5e-05
(11, 13) -> 5e-05
(4, 6) -> 5.1e-05
(6, 8) -> 5.1e-05
(8, 10) -> 5.1e-05
(10, 11) -> 5.1e-05
(0, 1) -> 5.1e-05
(1, 3) -> 5.1e-05
(5, 7) -> 5.1e-05
(7, 10) -> 7.6e-05
(3, 5) -> 7.6e-05
(0, 2) -> 7.6e-05
(2, 4) -> 7.6e-05
(5, 8) -> 7.6e-05
(6, 9) -> 7.6e-05
(10, 12) -> 7.6e-05
(9, 11) -> 7.7e-05
(4, 7) -> 7.7e-05
(3, 6) -> 0.000101
(1, 4) -> 0.000101
(5, 9) -> 0.000101
(10, 13) -> 0.000101
(2, 5) -> 0.000102
(4, 8) -> 0.000102
(6, 10) -> 0.000102
(8, 11) -> 0.000102
(9, 12) -> 0.000102
(0, 3) -> 0.000102
(7, 11) -> 0.000127
(1, 5) -> 0.000127
(2, 6) -> 0.000127
(3, 7) -> 0.000127
(4, 9) -> 0.000127
(5, 10) -> 0.000127
(8, 12) -> 0.000127
(9, 13) -> 0.000127
(7, 12) -> 0.000152
(0, 4) -> 0.000152
(1, 6) -> 0.000152
(3, 8) -> 0.000152
(8, 13) -> 0.000152
(4, 10) -> 0.000153
(6, 11) -> 0.000153
(2, 7) -> 0.000153
(7, 13) -> 0.000177
(3, 9) -> 0.000177
(0, 5) -> 0.000178
(2, 8) -> 0.000178
(5, 11) -> 0.000178
(6, 12) -> 0.000178
(1, 7) -> 0.000178
(0, 6) -> 0.000203
(1, 8) -> 0.000203
(2, 9) -> 0.000203
(3, 10) -> 0.000203
(5, 12) -> 0.000203
(6, 13) -> 0.000203
(4, 11) -> 0.000204
(1, 9) -> 0.000228
(5, 13) -> 0.000228
(2, 10) -> 0.000229
(4, 12) -> 0.000229
(0, 7) -> 0.000229
(0, 8) -> 0.000254
(1, 10) -> 0.000254
(3, 11) -> 0.000254
(4, 13) -> 0.000254
(0, 9) -> 0.000279
(3, 12) -> 0.000279
(2, 11) -> 0.00028
(3, 13) -> 0.000304
(0, 10) -> 0.000305
(1, 11) -> 0.000305
(2, 12) -> 0.000305
(1, 12) -> 0.00033
(2, 13) -> 0.00033
(1, 13) -> 0.000355
(0, 11) -> 0.000356
(0, 12) -> 0.000381
(0, 13) -> 0.000406
所以最小的集是(data [7],data [8])。
如果您在计算之前对数据进行排序,则可以从输出中轻松找到这样的最小距离。
#This is pseudocode.
x[n] = data[n] - data[n-1]
min(x)
答案 1 :(得分:1)
根据您当前的输入文件,这将计算最近的距离:
sort -u points.dat \
| xargs -n2 -J{} echo {} "- _1 * p" \
| dc \
| sort -u \
| sed '1!d;s/^\./0./'
您应该发布最终输入文件以获得更准确的答案。
修改强> 嗯,这将是一个更准确的距离:
sort -n -u points.dat \
| sed '1!p' \
| paste - - \
| head -n-1 \
| xargs -n2 -I{} echo {} "- _1 * p" \
| dc \
| sort -n -u \
| sed '1!d;s/^\./0./'