0

0

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

DDD

DDD

发布时间:2025-11-12 13:46:01

|

779人浏览过

|

来源于php中文网

原创

在 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 实现也采用了这种策略:

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

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 也无济于事。

    Sora
    Sora

    Sora是OpenAI发布的一种文生视频AI大模型,可以根据文本指令创建现实和富有想象力的场景。

    下载
  • 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 pl
import numpy as np # 仅用于上下文说明,实际解决方案中避免使用 np.NaN

def 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,周期为 10
ema_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
]

可以看到,

热门AI工具

更多
DeepSeek
DeepSeek

幻方量化公司旗下的开源大模型平台

豆包大模型
豆包大模型

字节跳动自主研发的一系列大型语言模型

WorkBuddy
WorkBuddy

腾讯云推出的AI原生桌面智能体工作台

腾讯元宝
腾讯元宝

腾讯混元平台推出的AI助手

文心一言
文心一言

文心一言是百度开发的AI聊天机器人,通过对话可以生成各种形式的内容。

讯飞写作
讯飞写作

基于讯飞星火大模型的AI写作工具,可以快速生成新闻稿件、品宣文案、工作总结、心得体会等各种文文稿

即梦AI
即梦AI

一站式AI创作平台,免费AI图片和视频生成。

ChatGPT
ChatGPT

最最强大的AI聊天机器人程序,ChatGPT不单是聊天机器人,还能进行撰写邮件、视频脚本、文案、翻译、代码等任务。

相关专题

更多
length函数用法
length函数用法

length函数用于返回指定字符串的字符数或字节数。可以用于计算字符串的长度,以便在查询和处理字符串数据时进行操作和判断。 需要注意的是length函数计算的是字符串的字符数,而不是字节数。对于多字节字符集,一个字符可能由多个字节组成。因此,length函数在计算字符串长度时会将多字节字符作为一个字符来计算。更多关于length函数的用法,大家可以阅读本专题下面的文章。

954

2023.09.19

C# ASP.NET Core微服务架构与API网关实践
C# ASP.NET Core微服务架构与API网关实践

本专题围绕 C# 在现代后端架构中的微服务实践展开,系统讲解基于 ASP.NET Core 构建可扩展服务体系的核心方法。内容涵盖服务拆分策略、RESTful API 设计、服务间通信、API 网关统一入口管理以及服务治理机制。通过真实项目案例,帮助开发者掌握构建高可用微服务系统的关键技术,提高系统的可扩展性与维护效率。

76

2026.03.11

Go高并发任务调度与Goroutine池化实践
Go高并发任务调度与Goroutine池化实践

本专题围绕 Go 语言在高并发任务处理场景中的实践展开,系统讲解 Goroutine 调度模型、Channel 通信机制以及并发控制策略。内容包括任务队列设计、Goroutine 池化管理、资源限制控制以及并发任务的性能优化方法。通过实际案例演示,帮助开发者构建稳定高效的 Go 并发任务处理系统,提高系统在高负载环境下的处理能力与稳定性。

38

2026.03.10

Kotlin Android模块化架构与组件化开发实践
Kotlin Android模块化架构与组件化开发实践

本专题围绕 Kotlin 在 Android 应用开发中的架构实践展开,重点讲解模块化设计与组件化开发的实现思路。内容包括项目模块拆分策略、公共组件封装、依赖管理优化、路由通信机制以及大型项目的工程化管理方法。通过真实项目案例分析,帮助开发者构建结构清晰、易扩展且维护成本低的 Android 应用架构体系,提升团队协作效率与项目迭代速度。

83

2026.03.09

JavaScript浏览器渲染机制与前端性能优化实践
JavaScript浏览器渲染机制与前端性能优化实践

本专题围绕 JavaScript 在浏览器中的执行与渲染机制展开,系统讲解 DOM 构建、CSSOM 解析、重排与重绘原理,以及关键渲染路径优化方法。内容涵盖事件循环机制、异步任务调度、资源加载优化、代码拆分与懒加载等性能优化策略。通过真实前端项目案例,帮助开发者理解浏览器底层工作原理,并掌握提升网页加载速度与交互体验的实用技巧。

97

2026.03.06

Rust内存安全机制与所有权模型深度实践
Rust内存安全机制与所有权模型深度实践

本专题围绕 Rust 语言核心特性展开,深入讲解所有权机制、借用规则、生命周期管理以及智能指针等关键概念。通过系统级开发案例,分析内存安全保障原理与零成本抽象优势,并结合并发场景讲解 Send 与 Sync 特性实现机制。帮助开发者真正理解 Rust 的设计哲学,掌握在高性能与安全性并重场景中的工程实践能力。

223

2026.03.05

PHP高性能API设计与Laravel服务架构实践
PHP高性能API设计与Laravel服务架构实践

本专题围绕 PHP 在现代 Web 后端开发中的高性能实践展开,重点讲解基于 Laravel 框架构建可扩展 API 服务的核心方法。内容涵盖路由与中间件机制、服务容器与依赖注入、接口版本管理、缓存策略设计以及队列异步处理方案。同时结合高并发场景,深入分析性能瓶颈定位与优化思路,帮助开发者构建稳定、高效、易维护的 PHP 后端服务体系。

458

2026.03.04

AI安装教程大全
AI安装教程大全

2026最全AI工具安装教程专题:包含各版本AI绘图、AI视频、智能办公软件的本地化部署手册。全篇零基础友好,附带最新模型下载地址、一键安装脚本及常见报错修复方案。每日更新,收藏这一篇就够了,让AI安装不再报错!

169

2026.03.04

Swift iOS架构设计与MVVM模式实战
Swift iOS架构设计与MVVM模式实战

本专题聚焦 Swift 在 iOS 应用架构设计中的实践,系统讲解 MVVM 模式的核心思想、数据绑定机制、模块拆分策略以及组件化开发方法。内容涵盖网络层封装、状态管理、依赖注入与性能优化技巧。通过完整项目案例,帮助开发者构建结构清晰、可维护性强的 iOS 应用架构体系。

246

2026.03.03

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
最新Python教程 从入门到精通
最新Python教程 从入门到精通

共4课时 | 22.5万人学习

Django 教程
Django 教程

共28课时 | 4.9万人学习

SciPy 教程
SciPy 教程

共10课时 | 1.9万人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号