简单移动平均线(SMA)的功能



i具有以下功能,用于计算Python中的SMA:

import numpy as np
def calcSma(data, smaPeriod):
    sma = []
    count = 0
    for i in xrange(data.size):
        if data[i] is None:
            sma.append(None)
        else:
            count += 1
            if count < smaPeriod:
                sma.append(None)
            else:
                sma.append(np.mean(data[i-smaPeriod+1:i+1]))
    return np.array(sma)

此功能有效,但我发现它几乎没有Pythonic。我不喜欢我在做的索引和计数,也不喜欢我必须附加到列表的方式,然后在返回之前将其变成一个numpy阵列。

我必须处理所有这些的原因,是因为我想以与输入数组相同的大小返回一个数组。这使得稍后在一般层面上绘制和处理变得更加容易。我可以轻松地做这样的事情:

sma = calcSma(data=data, smaPeriod=20)
sma2 = calcSma(data=sma, smaPeriod=10)
plt.plot(data)
plt.plot(sma)
plt.plot(sma2)
plt.show()

那么,关于如何更漂亮,更重要的Pythonic的任何想法?

Pythonic,我希望

import numpy as np

def calcSma(data, smaPeriod):
    j = next(i for i, x in enumerate(data) if x is not None)
    our_range = range(len(data))[j + smaPeriod - 1:]
    empty_list = [None] * (j + smaPeriod - 1)
    sub_result = [np.mean(data[i - smaPeriod + 1: i + 1]) for i in our_range]
    return np.array(empty_list + sub_result)

这是移动平均值的另一个实现,仅使用标准python库:

from collections import deque
import itertools
def moving_average(iterable, n=3):
    # http://en.wikipedia.org/wiki/Moving_average
    it = iter(iterable) 
    # create an iterable object from input argument
    d = deque(itertools.islice(it, n-1))  
    # create deque object by slicing iterable
    d.appendleft(0)
    s = sum(d)
    for elem in it:
        s += elem - d.popleft()
        d.append(elem)
        yield s / n
# example on how to use it
for i in  moving_average([40, 30, 50, 46, 39, 44]):
    print(i)
# 40.0
# 42.0
# 45.0
# 43.0

最新更新