我需要“压缩”代表信号的python数组的大小。信号看起来像下面的例子。在
signal = [
[0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1], #time values
[1,1,1,2,3,4,4,4,4,2,1,1] #function values
]
压缩后,信号应该如下面的代码所示。在
^{pr2}$你看,如果有一些区域具有常量值,那么只存储该区域的第一个和最后一个值。在
我写了下面的算法来实现这一点。在
signal_compressed = [[],[]]
old_value = None
for index, value in enumerate(signal[1]):
if value != old_value:
if index > 0:
if signal_compressed[0][-1] != signal[0][index - 1]:
signal_compressed[0].append(signal[0][index - 1])
signal_compressed[1].append(signal[1][index - 1])
signal_compressed[0].append(signal[0][index])
signal_compressed[1].append(value)
old_value = value
if signal_compressed[0][-1] < signal[0][-1]:
signal_compressed[0].append(signal[0][-1])
signal_compressed[1].append(signal[1][-1])
这个算法运行良好。对于有很多恒定段的信号,他工作得相当快。但是如果我试图压缩没有恒定段的信号(例如正弦信号或噪声信号),算法的工作速度非常慢。在
如何加速算法并保留功能?在
这里有一种使用发电机的方法:
我认为转置
^{pr2}$signal
,这样存储会更自然:这样,这两个
zip(*seq)
调用就不需要了,整个处理过程可以动态完成。在最后,如果这对于大的输入来说仍然太慢,那么使用NumPy可能是值得的。下面是一个这样的解决方案的概要:
您可以使用
itertools.groupby()
:相关问题 更多 >
编程相关推荐