python - 在python中计算指数移动平均线

标签 python signal-processing average digital-filter

我有一个日期范围和每个日期的测量值。我想计算每个日期的指数移动平均值。有人知道怎么做吗?

我是 python 新手。标准python库中似乎没有内置平均值,这让我觉得有点奇怪。也许我没找对地方。

那么,给定以下代码,我如何计算日历日期的 IQ 点的移动加权平均值?

from datetime import date
days = [date(2008,1,1), date(2008,1,2), date(2008,1,7)]
IQ = [110, 105, 90]

(可能有更好的方法来构建数据,任何建议将不胜感激)

最佳答案

编辑: 看来mov_average_expw()函数来自 scikits.timeseries.lib.moving_funcs来自 SciKits 的子模块(补充 SciPy 的附加工具包)更适合您问题的措辞。


计算 exponential smoothing具有平滑因子 alpha 的数据(在维基百科中它是 (1 - alpha)):

>>> alpha = 0.5
>>> assert 0 < alpha <= 1.0
>>> av = sum(alpha**n.days * iq 
...     for n, iq in map(lambda (day, iq), today=max(days): (today-day, iq), 
...         sorted(zip(days, IQ), key=lambda p: p[0], reverse=True)))
95.0

上面的不太好看,我们稍微重构一下:

from collections import namedtuple
from operator    import itemgetter

def smooth(iq_data, alpha=1, today=None):
    """Perform exponential smoothing with factor `alpha`.

    Time period is a day.
    Each time period the value of `iq` drops `alpha` times.
    The most recent data is the most valuable one.
    """
    assert 0 < alpha <= 1

    if alpha == 1: # no smoothing
        return sum(map(itemgetter(1), iq_data))

    if today is None:
        today = max(map(itemgetter(0), iq_data))

    return sum(alpha**((today - date).days) * iq for date, iq in iq_data)

IQData = namedtuple("IQData", "date iq")

if __name__ == "__main__":
    from datetime import date

    days = [date(2008,1,1), date(2008,1,2), date(2008,1,7)]
    IQ = [110, 105, 90]
    iqdata = list(map(IQData, days, IQ))
    print("\n".join(map(str, iqdata)))

    print(smooth(iqdata, alpha=0.5))

例子:

$ python26 smooth.py
IQData(date=datetime.date(2008, 1, 1), iq=110)
IQData(date=datetime.date(2008, 1, 2), iq=105)
IQData(date=datetime.date(2008, 1, 7), iq=90)
95.0

关于python - 在python中计算指数移动平均线,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/488670/

相关文章:

python - Django 模型 : how to return a default value in case of a non-existing foreign-key relationship?

windows - 用于 VC++ 的 DSP SDK

算法:寻找具有不完美值的模式

python每秒平均比特率

python - 如何快速不返回加盐密文,而是返回未加盐版本的密文?

python - 基于同一模型中其他字段的django字段验证器

r - 如何使用R中的seewave库从wave文件中获取最小和最大基频

R 取一个变量在另一个数值变量的区间内的平均值

python - 从具有相同名称的脚本导入已安装的包引发 "AttributeError: module has no attribute"或 ImportError 或 NameError

matlab - 如何在MATLAB中实现阶次分析