在 Polars 中高效计算指数移动平均线 (EMA) 并避免常见陷阱

在 polars 中高效计算指数移动平均线 (ema) 并避免常见陷阱

本教程详细介绍了如何在 Polars 中计算指数移动平均线 (EMA)。文章首先解释了 EMA 的基本概念和 Polars 中 `ewm_mean` 方法的使用。接着,重点阐述了在 Polars 中处理空值(`None` 与 `np.NaN`)的关键差异,并提供了一个经过优化的 `polars_ema` 函数,该函数能够正确地进行初始 SMA 计算并生成准确的 EMA 序列,从而避免了常见的 `NaN` 输出问题。

1. 理解指数移动平均线 (EMA) 与 Polars ewm_mean

指数移动平均线 (EMA) 是一种常用的技术分析指标,用于平滑价格数据,识别趋势。与简单移动平均线 (SMA) 不同,EMA 对近期数据赋予更高的权重,使其对市场变化更加敏感。

在 Polars 中,计算指数加权平均值的主要方法是 Series 上的 ewm_mean() 函数。该函数提供了灵活的参数来控制加权方式和空值处理,是进行金融数据分析的强大工具

2. EMA 计算中的初始 SMA 处理

在某些 EMA 的实现中,为了提供一个稳定的起始点,通常会在计算前 length 个数据点的 EMA 之前,先计算这 length 个点的简单移动平均 (SMA) 作为第一个有效的 EMA 值。这种处理方式旨在减少初期数据对 EMA 曲线的剧烈影响,使 EMA 曲线更加平滑地启动。

本教程所讨论的 EMA 实现也采用了这种策略:

首先计算前 length 个数据点的 SMA。将 SMA 值放置在序列的第 length – 1 个位置(索引从 0 开始),而前 length – 1 个位置则填充空值。随后对整个序列应用指数加权平均计算。

3. Polars 中空值处理的陷阱:np.NaN 与 None

在 Polars 中处理空值时,理解 numpy.NaN (简称 np.NaN) 和 Python 原生的 None 之间的区别至关重要。尽管两者都表示缺失数据,但 Polars 在内部处理和函数行为上可能有所不同。

np.NaN: 这是一个浮点数类型,表示“非数字”。当在 Polars Series 中使用 np.NaN 填充空值时,Polars 的某些函数(例如 ewm_mean)可能不会将其直接识别为可忽略的缺失值,而是作为一个特殊的浮点值进行处理。这可能导致在计算过程中,如果序列包含 np.NaN,ewm_mean 会在结果中产生一连串的 NaN,即使 min_periods 设置为 0 也无济于事。

None: Python 原生的 None 在 Polars 中通常被更明确地视为缺失值。当创建一个包含 None 的 Series 并指定浮点类型(例如 pl.Float64)时,Polars 能够正确地将这些 None 值识别为缺失数据,并在 ewm_mean 等函数中根据 ignore_nulls 和 min_periods 参数进行适当处理。

问题的根源:原始代码在使用 pl.Series([npNaN] * (length – 1)) 创建空值序列时,ewm_mean 无法正确处理这些 np.NaN,导致输出全为 NaN。

4. 在 Polars 中正确实现 EMA 计算

为了在 Polars 中正确计算 EMA,并妥善处理初始 SMA 填充的空值,我们需要使用 Python 原生的 None 来表示缺失数据。以下是修正后的 polars_ema 函数实现:

import polars as plimport numpy as np # 仅用于上下文说明,实际解决方案中避免使用 np.NaNdef polars_ema(close: pl.Series, length: int = 10, adjust: bool = False, sma: bool = True) -> pl.Series:    """    计算 Polars Series 的指数移动平均线 (EMA)。    该函数支持在 EMA 启动前,使用初始简单移动平均线 (SMA) 进行预热。    参数:        close (pl.Series): 输入的收盘价或待计算 EMA 的序列。        length (int): EMA 的周期长度。默认为 10。        adjust (bool): 是否根据权重调整除数。默认为 False (与 pandas 默认行为不同)。                       当 adjust=False 时,递归公式为:                       EMA_t = (value_t + (1 - alpha) * EMA_{t-1}) / (1 + (1 - alpha))                       当 adjust=True 时,EMA_t = alpha * value_t + (1 - alpha) * EMA_{t-1}                       其中 alpha = 2 / (span + 1),span 等同于 length。        sma (bool): 是否使用初始 SMA 来预热 EMA 序列。如果为 True,                    前 length-1 个值将为 None,第 length 个值为前 length 个数据的 SMA。                    默认为 True。    返回:        pl.Series: 包含 EMA 值的 Polars Series。    """    if close is None or not isinstance(close, pl.Series):        raise ValueError("输入 'close' 必须是一个 Polars Series。")    if length <= 0:        raise ValueError("长度 'length' 必须大于 0。")    # 如果需要初始 SMA 预热    if sma:        # 1. 计算前 'length' 个数据点的简单移动平均 (SMA)        # 使用 .head(length) 或 .slice(0, length) 获取前 'length' 个元素        initial_sma = close.head(length).mean()        # 2. 创建一个包含 'length - 1' 个 None 值的 Polars Series        # 关键修复: 使用 None 而不是 np.NaN,并明确指定 dtype 为 pl.Float64        nones_series = pl.Series([None] * (length - 1), dtype=pl.Float64)        # 3. 创建一个包含初始 SMA 值的 Series        initial_sma_series = pl.Series([initial_sma], dtype=pl.Float64)        # 4. 拼接 None 值和初始 SMA 值,形成序列的起始部分        initial_segment = nones_series.append(initial_sma_series)        # 5. 获取原始 'close' 序列中剩余的部分        rest_of_close = close.slice(length, close.len())        # 6. 将起始部分与剩余部分拼接,形成最终用于 EMA 计算的序列        # 这个序列的结构是:[None, ..., None, initial_sma, close_value_at_length, ...]        processed_close = initial_segment.append(rest_of_close)    else:        # 如果不使用 SMA 预热,直接使用原始序列        processed_close = close    # 7. 计算指数移动平均线 (EMA)    # span: 等同于 length,用于计算指数加权因子 alpha = 2 / (span + 1)    # adjust: 控制加权方式,False 更接近 pandas 的默认行为(但仍有细微差异,需注意)    # ignore_nulls=False: 表示在计算时考虑 None 值。    # min_periods=0: 允许在少于 span 个非 None 值的情况下也进行计算。    #                这对于我们序列开头有 None 的情况至关重要,确保从第一个有效值开始计算。    ema = processed_close.ewm_mean(        span=length,        adjust=adjust,        ignore_nulls=False,        min_periods=0    )    return ema

5. 示例与应用

让我们使用提供的示例数据来测试上述修正后的 polars_ema 函数。

import polars as pl# 示例数据sample_data = pl.Series('close', [    1.08086, 1.08069, 1.08077, 1.08077, 1.08052,    1.08055, 1.08068, 1.08073, 1.08077, 1.08073,    1.08068, 1.08062, 1.08052, 1.0806, 1.08063,    1.08064, 1.08063, 1.08053, 1.08067, 1.08058], dtype=pl.Float64)# 调用函数计算 EMA,周期为 10ema_result = polars_ema(sample_data, length=10)print("计算得到的 EMA 序列:")print(ema_result)

预期输出(部分):

Shape: (20,)Series: 'close' [f64][    null    null    null    null    null    null    null    null    null    1.080697    // 前10个值的SMA    1.0806954545454546    1.0806817363636363    1.0806605115684364    1.0806489640096297    1.0806436978260244    1.0806429345851109    1.0806414919332726    1.080628493399945    1.0806387673272277    1.0806298914500954]

可以看到,

以上就是在 Polars 中高效计算指数移动平均线 (EMA) 并避免常见陷阱的详细内容,更多请关注创想鸟其它相关文章!

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 chuangxiangniao@163.com 举报,一经查实,本站将立刻删除。
发布者:程序猿,转转请注明出处:https://www.chuangxiangniao.com/p/1379582.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
Python环境管理:解决Pip更新时的权限问题 (WinError 5)
上一篇 2025年12月14日 20:47:15
在Polars中高效计算指数移动平均线(EMA)及其初始化策略
下一篇 2025年12月14日 20:47:23

相关推荐

  • composer require-dev和require有什么不同_Composer Require与Require-Dev区别解析

    require用于声明项目运行必需的依赖,如框架、数据库组件和第三方SDK,这些包会随项目部署到生产环境;2. require-dev用于声明仅在开发和测试阶段需要的工具,如PHPUnit、PHPStan、Faker等,不会默认部署到生产环境;3. 安装时composer install根据环境决定…

    2026年5月10日
    1000
  • Matplotlib 地图中多类型图例的创建与优化

    Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化

    本教程旨在解决matplotlib地图可视化中,如何在一个图例中同时展示颜色块(如区域分类)和自定义标记(如特定兴趣点)的问题。文章详细介绍了当传统`patch`对象无法正确显示标记时,如何利用`matplotlib.lines.line2d`创建标记图例句柄,并将其与颜色块图例句柄合并,从而生成一…

    2026年5月10日 用户投稿
    300
  • Golang JSON序列化:控制敏感字段暴露的最佳实践

    本教程探讨golang中如何高效控制结构体字段在json序列化时的可见性。当需要将包含敏感信息的结构体数组转换为json响应时,通过利用`encoding/json`包提供的结构体标签,特别是`json:”-“`,可以轻松实现对特定字段的忽略,从而避免敏感数据泄露,确保api…

    2026年5月10日
    000
  • 利用海象运算符简化条件赋值:Python教程与最佳实践

    本文旨在探讨Python中海象运算符(:=)在条件赋值场景下的应用。通过对比传统if/else语句与海象运算符,以及条件表达式,分析海象运算符在简化代码、提高可读性方面的优势与局限性。并通过具体示例,展示如何在列表推导式等场景下合理使用海象运算符,同时强调其潜在的复杂性及替代方案,帮助开发者更好地掌…

    2026年5月10日
    100
  • Debian syslog性能优化技巧有哪些

    提升Debian系统syslog (通常基于rsyslog)性能,关键在于精简配置和高效处理日志。以下策略能有效优化日志管理,提升系统整体性能: 精简配置,高效加载: 在rsyslog配置文件中,仅加载必要的输入、输出和解析模块。 使用全局指令设置日志级别和格式,避免不必要的处理。 自定义模板: 创…

    2026年5月10日
    000
  • 比特币新手教程 比特币交易平台有哪些

    比特币是一种去中心化的数字货币,基于区块链技术实现点对点交易,具有匿名性、有限发行和不可篡改等特点;新手可通过交易所购买,P2P交易获得比特币,常用平台包括Binance、OKX和Huobi;交易流程包括注册账户、实名认证、绑定支付方式、充值法币并下单购买,可选择市价单或限价单;比特币存储方式有交易…

    2026年5月10日
    000
  • c++中的SFINAE技术是什么_c++模板编程中的SFINAE原理与应用

    SFINAE 是“替换失败不是错误”的原则,指模板实例化时若参数替换导致错误,只要存在其他合法候选,编译器不报错而是继续重载决议。它用于条件启用模板、类型检测等场景,如通过 decltype 或 enable_if 控制函数重载,实现类型特征判断。尽管 C++20 引入 Concepts 简化了部分…

    2026年5月10日
    000
  • Go语言mgo查询构建:深入理解bson.M与日期范围查询的正确实践

    本文旨在解决go语言mgo库中构建复杂查询时,特别是涉及嵌套`bson.m`和日期范围筛选的常见错误。我们将深入剖析`bson.m`的类型特性,解释为何直接索引`interface{}`会导致“invalid operation”错误,并提供一种推荐的、结构清晰的代码重构方案,以确保查询条件能够正确…

    2026年5月10日
    100
  • RichHandler与Rich Progress集成:解决显示冲突的教程

    在使用rich库的`richhandler`进行日志输出并同时使用`progress`组件时,可能会遇到显示错乱或溢出问题。这通常是由于为`richhandler`和`progress`分别创建了独立的`console`实例导致的。解决方案是确保日志处理器和进度条组件共享同一个`console`实例…

    2026年5月10日
    000
  • Golang goroutine与channel调试技巧

    使用go run -race检测数据竞争,结合runtime.NumGoroutine监控协程数量,通过pprof分析阻塞调用栈,利用select超时避免永久阻塞,有效排查goroutine泄漏、死锁和数据竞争问题。 Go语言的goroutine和channel是并发编程的核心,但它们也带来了调试上…

    2026年5月10日
    000
  • 使用 Jupyter Notebook 进行探索性数据分析

    Jupyter Notebook通过单元格实现代码与Markdown结合,支持数据导入(pandas)、清洗(fillna)、探索(matplotlib/seaborn可视化)、统计分析(describe/corr)和特征工程,便于记录与分享分析过程。 Jupyter Notebook 是进行探索性…

    2026年5月10日
    000
  • 《魔兽世界》将于6月11日开启国服回归技术测试

    《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试

    《%ign%ignore_a_1%re_a_1%》官方宣布,将于6月11日开启国服回归技术测试,时间为7天,并称可以在6月内正式开服,玩家们可以访问官网下载战网客户端并预下载“巫妖王之怒”客户端,技术测试详情见下图。 WordAi WordAI是一个AI驱动的内容重写平台 53 查看详情 以上就是《…

    2026年5月10日 用户投稿
    200
  • php常量怎么用_PHP常量(define/const)定义与使用方法

    PHP中可通过define函数和const关键字定义常量,用于存储不可变值。define适用于全局作用域,支持动态名称和条件定义,如define(‘SITE_NAME’, ‘MyWebsite’);const在编译时生效,语法简洁但限制多,只能在类或全…

    2026年5月10日
    000
  • 如何在HTML中插入表单元素_HTML表单控件与输入类型使用指南

    HTML表单通过标签构建,包含action和method属性定义数据提交目标与方式,常用input类型如text、password、email等适配不同输入需求,配合label、required、placeholder提升可用性,结合textarea、select、button等控件实现完整交互,是…

    2026年5月10日
    100
  • 网站标题关键词更新后,搜索引擎为何仍显示旧标题?

    网站标题更新后,搜索引擎为何显示旧标题? 网站SEO优化中,站长常修改网站标题关键词,期望搜索结果显示自定义标题。然而,即使更新标签、meta keywords、meta description和结构化数据中的name属性后,搜索结果仍显示旧标题,这令人费解。本文将对此进行解释。 问题:站长修改了网…

    2026年5月10日
    300
  • 创建指定大小并填充特定数据的Golang文件教程

    本文将介绍如何使用Golang创建一个指定大小的文件,并用特定数据填充它。我们将使用 `os` 包提供的函数来创建和截断文件,从而实现快速生成大文件的目的。示例代码展示了如何创建一个10MB的文件,并将其填充为全零数据。掌握这些方法,可以方便地在例如日志系统或磁盘队列等场景中,预先创建测试文件或初始…

    2026年5月10日
    000
  • 深入理解 Express.js 中 next() 参数的作用与中间件机制

    本文深入探讨 express.js 中间件函数中的 `next()` 参数。它负责将控制权传递给请求-响应周期中的下一个中间件或路由处理程序。文章将详细解释 `next()` 的工作原理、中间件的注册与执行顺序,以及不正确使用 `next()` 可能导致请求挂起的风险,并通过代码示例和实际应用场景,…

    2026年5月10日
    000
  • Python命令怎样使用profile分析脚本性能 Python命令性能分析的基础教程

    使用Python的cProfile模块分析脚本性能最直接的方式是通过命令行执行python -m cProfile your_script.py,它会输出每个函数的调用次数、总耗时、累积耗时等关键指标,帮助定位性能瓶颈;为进一步分析,可将结果保存为文件python -m cProfile -o ou…

    2026年5月10日
    000
  • 使用 WebCodecs VideoDecoder 实现精确逐帧回退

    本文档旨在解决在使用 WebCodecs VideoDecoder 进行视频解码时,实现精确逐帧回退的问题。通过比较帧的时间戳与目标帧的时间戳,可以避免渲染中间帧,从而提高用户体验。本文将提供详细的解决方案和示例代码,帮助开发者实现精确的视频帧控制。 在使用 WebCodecs VideoDecod…

    2026年5月10日
    000
  • 如何插入查询结果数据_SQL插入Select查询结果方法

    如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法

    使用INSERT INTO…SELECT语句可高效插入数据,通过NOT EXISTS、LEFT JOIN、MERGE语句或唯一约束避免重复;表结构不一致时可通过别名、类型转换、默认值或计算字段处理;结合存储过程可提升可维护性,支持参数化与动态SQL。 将查询结果数据插入到另一个表中,可以…

    2026年5月10日 用户投稿
    300

发表回复

登录后才能评论
关注微信