尧图网站开发 
SymPy 一元幂级数环完全指南Power Series Ring 的精度语义、运算体系与 FLINT 后端【免费下载链接】sympyA computer algebra system written in pure Python项目地址: https://gitcode.com/GitHub_Trending/sy/sympy本文基于 SymPy 仓库 doc/src/modules/polys/series.rst 及其对应实现系统讲解一元幂级数环Univariate Power Series Ring模块的完整用法如何通过power_series_ring创建环、理解其独特的精确多项式自动保真、超精度自动截断语义、在整数环ZZ与有理数域QQ上执行加减乘除、求逆、复合、微积分及log/exp/三角级数展开并深入双后端架构纯 Python 实现与 FLINT 高性能后端的自动选择原理。读完本文你将能够像使用普通多项式一样精确且高效地处理带截断误差的幂级数计算。一、幂级数环概念与设计动机幂级数Power Series是形如a_0 a_1*x a_2*x**2 ...的无穷级数。在实际计算中无法表示无穷项因此 SymPy 的幂级数环模块引入精度precision概念对于一个精度为prec的环级数只保留次数不超过prec - 1的项任何会产生次数 prec项的操作都会被截断截断结果用大 O 记号O(x**prec)显式标注见 sympy/polys/series/base.py 中series_pprint的实现截断时输出形如poly O(sym**prec)。这一设计解决了数值与符号计算中的经典矛盾精确性当参与运算的都是有限多项式且结果次数仍小于环精度时返回无任何截断的精确多项式源码中精度记录为None完备性一旦运算结果触及精度边界则自动截断并携带O(x**prec)保证每次运算结果仍是环内元素。模块目前的适用范围是一元univariate幂级数底层域仅支持整数环ZZ与有理数域QQ两种源码 sympy/polys/series/ring.py 对K.is_ZZ/K.is_QQ分支处理其他域抛出ValueError: Unsupported ground domain。二、快速上手创建环并做算术原文档给出的最小示例是创建精度为 8 的整数幂级数环 from sympy.polys.series import power_series_ring from sympy import ZZ, QQ R, x power_series_ring(x, ZZ, 8) f 1 2*x 3*x**2 g 4 x工厂函数返回一个二元组第一个元素是环实例第二个元素是环的生成元generator即变量x。此后即可直接用 Python 运算符操作级数 f g 5 3*x 3*x**2 f * g 4 9*x 14*x**2 3*x**3加法与乘法都符合多项式直觉——因为f、g及运算结果次数均低于 8系统直接返回精确结果__repr__输出的字符串由环的pretty方法生成见 sympy/polys/series/ring.py。三、核心语义何时精确、何时截断这是本模块最重要的设计。原文档用两组对照示例说明场景一结果次数低于环精度返回精确多项式 p 2 - 3*x q 7*x 6*x**2 x**3 p * q 14*x - 9*x**2 - 16*x**3 - 3*x**4p*q的最高次数为 4小于精度 8因此结果是数学上完全精确的多项式不带任何O(x**8)。这在源码 sympy/polys/series/ringpython.py 的_useries辅助函数中实现当series_prec is None且deg ring_prec时直接返回(coeffs, None)——注意第二个元素None正是精确多项式的标记乘法_useries_mul中若两个操作数都精确且次数和小于环精度同样走dup_mul并返回None精度ringpython.py。场景二结果超出精度边界自动截断并标注 O r 1 2*x 3*x**2 4*x**3 5*x**4 6*x**5 7*x**6 8*x**7 s x x**2 r * s x 3*x**2 5*x**3 7*x**4 9*x**5 11*x**6 13*x**7 O(x**8)r本身已是满精度8 项与s相乘会产生高达 9 次的项系统自动截断到精度 8并追加O(x**8)。对应源码路径_useries_mul中经_unify_prec统一精度后调用dup_series_mul(coeffs1, coeffs2, min_prec, dom)完成带截断的级数乘法ringpython.py。精度合并规则可从_unify_prec源码推断两个精确级数None相遇时取环精度一方精确一方有精度时取有精度一方的值双方都有精度时取较小者并统一截断系数列表。测试 sympy/polys/series/tests/test_ring.py 验证了R.series_prec(R([1, 2, 3, 4, 5, 6, 7])) 6即默认环精度 6 下输入 7 个系数会被截断为 6 项。四、工厂函数power_series_ring详解函数签名含类型重载位于 sympy/polys/series/ring.pydef power_series_ring(symbol: str, K: Domain, prec: int 6) - tuple[PowerSeriesRingRing | PowerSeriesRingField, PowerSeriesElement]参数类型说明symbolstr也可传Expr级数变量名如x内部通过Symbol(symbol)转为 SymPy 符号KDomain底域仅支持ZZ与QQ传其他域抛出ValueErrorprecint默认6环的默认精度必须非负底层__init__对prec 0抛ValueError: Power series precision must be non-negative返回的二元组(ring, gen)其中ring是PowerSeriesRingRing底域为环时或PowerSeriesRingField底域为域时gen是生成元PowerSeriesElement类型。选择逻辑非常直接K.is_ZZ→PowerSeriesRingRingK.is_QQ→PowerSeriesRingField。值得注意PowerSeriesRingRing/PowerSeriesRingField的构造函数同样接受(domain, symbolx, prec6)可以直接实例化例如文档用例PowerSeriesRingRing(ZZ, x, 5)。但官方推荐始终经由工厂函数因为它内部通过 sympy/polys/series/tring.py 的_power_series_ring完成了后端的透明选择。五、两种环类Ring 与 Field 的能力差异原文档将高层 API 分为两个类PowerSeriesRingRing整数环ZZ之上除基本四则运算外支持square平方、compose复合、inverse乘法逆、reversion复合逆/反演、differentiate求导、truncate截断、from_expr/to_expr/from_list/from_int/from_ground等转换方法以及order_term()生成O(x**prec)。PowerSeriesRingField有理数域QQ之上继承 Ring 的全部能力额外提供integrate积分、sqrt平方根以及一整套超越函数log、log1p、exp、expm1、atan、atanh、asin、asinh、tan、tanh、sin、sinh、cos、coshsympy/polys/series/ring.py。原因很直观ZZ上的系数无法自由做除法与积分1/2不属于整数而QQ是域所有非零系数都可逆因此微分、积分与解析函数级数展开成为可能。5.1 环级操作示例来自PowerSeriesRingRing源码 docstring from sympy import ZZ from sympy.polys.series import PowerSeriesRingRing R PowerSeriesRingRing(ZZ, x, 5) x R.gen R.square(1 x x**2) 1 2*x 3*x**2 2*x**3 x**4 R.inverse(1 x) 1 - x x**2 - x**3 x**4 O(x**5) R.compose(3 x 14*x**2, x**2 x**3) 3 x**2 x**3 14*x**4 O(x**5) R.reversion(x 4*x**2 8*x**3) x - 4*x**2 24*x**3 - 160*x**4 O(x**5) R.differentiate(10 x x**2 x**3 x**4 R.order_term()) 1 2*x 3*x**2 4*x**3 O(x**4)注意各操作的前置条件源码中均有显式校验compose(s, t)t的常数项必须为 0否则ValueError: Series composition requires the constant term of the second series to be zeroinverse(s)s的常数项必须是底域中的单位NotReversible: Series inverse requires the constant term to be a unitreversion(s)常数项必须为 0、一次项系数必须是单位square(s)即multiply(s, s)。5.2 域级操作示例来自PowerSeriesRingField源码 docstring from sympy import QQ from sympy.polys.series import PowerSeriesRingField R PowerSeriesRingField(QQ, x) x R.gen R.sqrt(1 2*x 4*x**2 8*x**3) 1 x 3/2*x**2 5/2*x**3 - 29/8*x**4 - 1/8*x**5 O(x**6) R.integrate(1 x x**2 x**3 x**4 x**5 R.order_term()) x 1/2*x**2 1/3*x**3 1/4*x**4 1/5*x**5 O(x**6)六、超越函数级数log、exp 与三角展开原文档专门演示了有理数域上常见的数学函数展开这是幂级数环最实用的场景之一 R, x power_series_ring(x, QQ, 8) R.log(1 x) x - 1/2*x**2 1/3*x**3 - 1/4*x**4 1/5*x**5 - 1/6*x**6 1/7*x**7 O(x**8) R.tan(x) x 1/3*x**3 2/15*x**5 17/315*x**7 O(x**8)所有超越函数都在 sympy/polys/series/ringpython.py 中以_useries_*前缀实现其输入前置条件非常统一见各函数源码的ValueError分支函数输入要求实现策略从源码可读log(s)常数项必须为 1项数 20 时用log(f)f/f再积分否则直接复合标准级数ln(1u)log1p(s)常数项必须为 0先加 1 再调logexp(s)常数项必须为 0项数 20 时用 Newton 迭代校正exp (s - log(exp))*exp否则复合e^u级数expm1(s)常数项必须为 0先exp再减 1sin/cos常数项必须为 0项数大时借助tan半角公式否则复合标准麦克劳林级数tan/tanh常数项必须为 0Newton 迭代y_{n1} y_n (s − atan(y_n))·(1y_n²)asin/atan等常数项必须为 0利用导数关系如atan(f)f/(1f²)积分或复合级数以原文档的R.log(1 x)为例log要求常数项为 11x恰好满足输入仅 2 项走复合标准级数分支——先构造系数[0, 1, -1/2, 1/3, ...]_useries_log中的c列表再与(1x)-1 x复合得到x - x²/2 x³/3 - ...ringpython.py。而tan则通过_useries_tan的 Newton 迭代逐步逼近到目标精度ringpython.py。这些函数同样支持级数作为自变量例如工厂函数 docstring 中的示例 from sympy import QQ from sympy.polys.series import power_series_ring R, x power_series_ring(x, QQ) R.sin(x x**2) x x**2 - 1/6*x**3 - 1/2*x**4 - 59/120*x**5 O(x**6) R.log(1 7*x**2) 7*x**2 - 49/2*x**4 O(x**6)七、元素类PowerSeriesElement运算符重载与转换PowerSeriesElementsympy/polys/series/ring.py是环中元素的封装类继承DomainElement实现了完整的 Python 运算符重载加减乘除-*/含与int或底域元素混合运算的__radd__/__rsub__/__rmul__/__rtruediv__幂**支持整数次幂n0返回环的onen1返回自身否则走_pow_int一元、-比较基于表示比较equal_repr不同环的元素相加/相乘会抛ValueError: Cannot add/multiply power series with different rings.。元素还提供以下实用方法/属性 s.as_expr() # 转回 SymPy 表达式含 O(x**prec) s.removeO() # 去掉大 O 记号返回纯多项式部分 s.prec # 该级数的精度int 或 None s.constant_coefficient() # 常数项 s.is_ground # 是否退化为常数底域元素环侧配套的转换方法包括from_expr(expr)将 SymPy 表达式嵌入环若表达式含Order项则按其精度截断见 ring.py 的_rebuild_expr、from_list、to_list、to_dense、to_expr。ring_new即R(...)调用协议可接受PowerSeriesElement、Expr、int或底域元素四种输入。八、双后端架构纯 Python 与 FLINT 的自动切换原文档明确说明系统含两个后端纯 Python 实现与基于FLINT库的高性能实现前者是回退方案后者在可用时被自动选用。这一机制在 sympy/polys/series/tring.py 中实现模块根据sympy.external.gmpy.GROUND_TYPES判断当前 SymPy 使用的整数/有理数底层类型gmpy或flinttring.py若GROUND_TYPES flint导入 sympy/polys/series/ringflint.py 中的FlintPowerSeriesRingZZ/FlintPowerSeriesRingQQ底层级数类型为ZZSeries fmpz_series | fmpz_poly、QQSeries fmpq_series | fmpq_poly否则将 FLINT 类置为None回退到 sympy/polys/series/ringpython.py 的PythonPowerSeriesRingZZ/PythonPowerSeriesRingQQ底层级数类型为USeries tuple[list[Er], int | None]系数列表 精度标记。PowerSeriesRingZZ(prec)/PowerSeriesRingQQ(prec)两个工厂函数执行具体选择tring.py再由_power_series_ring(K, prec)根据底域分派到 ZZ 或 QQ 版本最终由高层power_series_ring包装为带生成元的结果。整个链路对用户透明。两个后端都实现了 sympy/polys/series/base.py 中PowerSeriesRingProto/PowerSeriesRingFieldProto协议定义的统一接口add/subtract/multiply/divide/pow_int/square/truncate/compose/inverse/reversion/differentiate域级再加integrate/sqrt/log/exp/tan/sin/cos等因此上层PowerSeriesElement无需关心具体后端。关于 FLINT 实现的两个工程细节ringflint.py_get_series_precision通过级数对象的prec属性获取精度对 python-flint ≤ 0.7.1 的旧版本退而解析repr字符串中的prec字段代码注释注明该兼容分支可在 python-flint 0.8 成为最低版本后移除_global_cap上下文管理器临时调整ctx.capFLINT 全局级数截断上限退出时恢复保证并发/嵌套环境下的精度隔离。此外 ringflint.py 的类声明带doctest_depends_on(ground_types[flint])装饰器表明其文档测试仅在 FLINT 后端下运行这从侧面印证了双后端的运行时条件依赖。九、与多项式域框架的集成作为 SymPy Domain 使用除独立使用的PowerSeriesRingRing之外模块还提供了领域domain层的集成类。原文档Domain for Univariate Power Series Rings一节指向 sympy/polys/domains/powerseriesring.pySeriesRingProto定义领域层接口的 Protocolring、dtype、domain、symbol、prec等属性以及from_expr/to_expr/domain_new/ring_new等方法PowerSeriesRing继承Ring与CompositeDomain的领域类其构造函数内部直接调用power_series_ring(str(symbol), domain, prec)powerseriesring.py并提供new、of_type、zero、one、prec、is_unit、canonical_unit、to_sympy、from_sympy等标准域方法使幂级数环可以嵌入 SymPy 的Domain体系例如供其他多项式域工具使用。__repr__形如ZZ[[x], 6]即底域 双括号 变量 精度的数学记号。类源码注释明确标注该领域类仍处于实验阶段experimental stage后续应补充与现有领域框架的完整集成使用时应注意这一前提。十、类型协议与测试验证10.1 协议层base.pysympy/polys/series/base.py 定义了三个核心协议/类型PowerSeriesRingProto[TSeries, Er]环级操作协议ZZ 后端PowerSeriesRingFieldProto[TSeries, Ef]域级操作协议继承前者并追加integrate/sqrt/log/exp及全部三角函数base.pyTSeries级数类型的 TypeVar而 sympy/polys/series/tring.py 中的TSeriesElement是仅供静态类型检查使用的底层级数占位类型。此外series_pprint(series, prec, symx, ascendingTrue)是统一的字符串渲染函数prec0时直接输出O(x**0)系数为空或全零时视prec输出O(x**prec)或0否则输出poly O(x**prec)base.py。这就是原文档示例中O(x**8)打印格式的来源。10.2 测试覆盖测试集中在 sympy/polys/series/tests/test_ring.py另有 test_hypothesis.py 用于基于假设的性质测试。从测试命名可见覆盖范围test_equal精度边界下的相等性、test_basicsone/zero/gen/prec/精度截断如R.series_prec(R([1,2,3,4,5,6,7])) 6、test_positive/test_negative正负号与精度收束、test_add/test_int_subtract加减法、精度统一、以及后续大量针对乘法、幂、逆、复合、反演、导数和超越函数的用例。测试通过参数化夹具ring_int/ring_rational同时跑纯 Python 与 FLINT 两种后端保证行为一致。十一、实战小结与注意事项回顾全文使用 SymPy 幂级数环的关键要点用工厂函数创建R, x power_series_ring(x, ZZ|QQ, prec)默认精度 6prec必须非负底域决定能力ZZ环只支持代数运算与求导、逆、复合、反演QQ域额外支持积分、开方与全部超越函数级数精度语义要牢记次数低于prec的结果是精确多项式打印无O触及prec自动截断并带O(x**prec)混合精度按取较小精度合并前置条件要满足inverse要求常数项为单位、reversion要求常数项为 0 且一次项为单位、compose要求第二参数常数项为 0、log要求常数项为 1、exp/sin/cos等要求常数项为 0违反时抛出带明确信息的异常后端自动选择安装 FLINT 支持GROUND_TYPES flint时自动启用高性能后端否则回退纯 Python接口完全一致转换与集成from_expr/to_expr/as_expr/removeO实现与 SymPy 表达式的双向互转PowerSeriesRing领域类可融入Domain框架但后者处于实验阶段。深入阅读建议本指南依据 doc/src/modules/polys/series.rst 撰写实现与验证请继续查看 sympy/polys/series/ring.py、sympy/polys/series/ringpython.py、sympy/polys/series/ringflint.py、sympy/polys/series/tring.py、sympy/polys/series/base.py、sympy/polys/domains/powerseriesring.py 与测试文件 sympy/polys/series/tests/test_ring.py。【免费下载链接】sympyA computer algebra system written in pure Python项目地址: https://gitcode.com/GitHub_Trending/sy/sympy创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考 
SymPy 一元幂级数环完全指南Power Series Ring 的精度语义、运算体系与 FLINT 后端【免费下载链接】sympyA computer algebra system written in pure Python项目地址: https://gitcode.com/GitHub_Trending/sy/sympy本文基于 SymPy 仓库 doc/src/modules/polys/series.rst 及其对应实现系统讲解一元幂级数环Univariate Power Series Ring模块的完整用法如何通过power_series_ring创建环、理解其独特的精确多项式自动保真、超精度自动截断语义、在整数环ZZ与有理数域QQ上执行加减乘除、求逆、复合、微积分及log/exp/三角级数展开并深入双后端架构纯 Python 实现与 FLINT 高性能后端的自动选择原理。读完本文你将能够像使用普通多项式一样精确且高效地处理带截断误差的幂级数计算。一、幂级数环概念与设计动机幂级数Power Series是形如a_0 a_1*x a_2*x**2 ...的无穷级数。在实际计算中无法表示无穷项因此 SymPy 的幂级数环模块引入精度precision概念对于一个精度为prec的环级数只保留次数不超过prec - 1的项任何会产生次数 prec项的操作都会被截断截断结果用大 O 记号O(x**prec)显式标注见 sympy/polys/series/base.py 中series_pprint的实现截断时输出形如poly O(sym**prec)。这一设计解决了数值与符号计算中的经典矛盾精确性当参与运算的都是有限多项式且结果次数仍小于环精度时返回无任何截断的精确多项式源码中精度记录为None完备性一旦运算结果触及精度边界则自动截断并携带O(x**prec)保证每次运算结果仍是环内元素。模块目前的适用范围是一元univariate幂级数底层域仅支持整数环ZZ与有理数域QQ两种源码 sympy/polys/series/ring.py 对K.is_ZZ/K.is_QQ分支处理其他域抛出ValueError: Unsupported ground domain。二、快速上手创建环并做算术原文档给出的最小示例是创建精度为 8 的整数幂级数环 from sympy.polys.series import power_series_ring from sympy import ZZ, QQ R, x power_series_ring(x, ZZ, 8) f 1 2*x 3*x**2 g 4 x工厂函数返回一个二元组第一个元素是环实例第二个元素是环的生成元generator即变量x。此后即可直接用 Python 运算符操作级数 f g 5 3*x 3*x**2 f * g 4 9*x 14*x**2 3*x**3加法与乘法都符合多项式直觉——因为f、g及运算结果次数均低于 8系统直接返回精确结果__repr__输出的字符串由环的pretty方法生成见 sympy/polys/series/ring.py。三、核心语义何时精确、何时截断这是本模块最重要的设计。原文档用两组对照示例说明场景一结果次数低于环精度返回精确多项式 p 2 - 3*x q 7*x 6*x**2 x**3 p * q 14*x - 9*x**2 - 16*x**3 - 3*x**4p*q的最高次数为 4小于精度 8因此结果是数学上完全精确的多项式不带任何O(x**8)。这在源码 sympy/polys/series/ringpython.py 的_useries辅助函数中实现当series_prec is None且deg ring_prec时直接返回(coeffs, None)——注意第二个元素None正是精确多项式的标记乘法_useries_mul中若两个操作数都精确且次数和小于环精度同样走dup_mul并返回None精度ringpython.py。场景二结果超出精度边界自动截断并标注 O r 1 2*x 3*x**2 4*x**3 5*x**4 6*x**5 7*x**6 8*x**7 s x x**2 r * s x 3*x**2 5*x**3 7*x**4 9*x**5 11*x**6 13*x**7 O(x**8)r本身已是满精度8 项与s相乘会产生高达 9 次的项系统自动截断到精度 8并追加O(x**8)。对应源码路径_useries_mul中经_unify_prec统一精度后调用dup_series_mul(coeffs1, coeffs2, min_prec, dom)完成带截断的级数乘法ringpython.py。精度合并规则可从_unify_prec源码推断两个精确级数None相遇时取环精度一方精确一方有精度时取有精度一方的值双方都有精度时取较小者并统一截断系数列表。测试 sympy/polys/series/tests/test_ring.py 验证了R.series_prec(R([1, 2, 3, 4, 5, 6, 7])) 6即默认环精度 6 下输入 7 个系数会被截断为 6 项。四、工厂函数power_series_ring详解函数签名含类型重载位于 sympy/polys/series/ring.pydef power_series_ring(symbol: str, K: Domain, prec: int 6) - tuple[PowerSeriesRingRing | PowerSeriesRingField, PowerSeriesElement]参数类型说明symbolstr也可传Expr级数变量名如x内部通过Symbol(symbol)转为 SymPy 符号KDomain底域仅支持ZZ与QQ传其他域抛出ValueErrorprecint默认6环的默认精度必须非负底层__init__对prec 0抛ValueError: Power series precision must be non-negative返回的二元组(ring, gen)其中ring是PowerSeriesRingRing底域为环时或PowerSeriesRingField底域为域时gen是生成元PowerSeriesElement类型。选择逻辑非常直接K.is_ZZ→PowerSeriesRingRingK.is_QQ→PowerSeriesRingField。值得注意PowerSeriesRingRing/PowerSeriesRingField的构造函数同样接受(domain, symbolx, prec6)可以直接实例化例如文档用例PowerSeriesRingRing(ZZ, x, 5)。但官方推荐始终经由工厂函数因为它内部通过 sympy/polys/series/tring.py 的_power_series_ring完成了后端的透明选择。五、两种环类Ring 与 Field 的能力差异原文档将高层 API 分为两个类PowerSeriesRingRing整数环ZZ之上除基本四则运算外支持square平方、compose复合、inverse乘法逆、reversion复合逆/反演、differentiate求导、truncate截断、from_expr/to_expr/from_list/from_int/from_ground等转换方法以及order_term()生成O(x**prec)。PowerSeriesRingField有理数域QQ之上继承 Ring 的全部能力额外提供integrate积分、sqrt平方根以及一整套超越函数log、log1p、exp、expm1、atan、atanh、asin、asinh、tan、tanh、sin、sinh、cos、coshsympy/polys/series/ring.py。原因很直观ZZ上的系数无法自由做除法与积分1/2不属于整数而QQ是域所有非零系数都可逆因此微分、积分与解析函数级数展开成为可能。5.1 环级操作示例来自PowerSeriesRingRing源码 docstring from sympy import ZZ from sympy.polys.series import PowerSeriesRingRing R PowerSeriesRingRing(ZZ, x, 5) x R.gen R.square(1 x x**2) 1 2*x 3*x**2 2*x**3 x**4 R.inverse(1 x) 1 - x x**2 - x**3 x**4 O(x**5) R.compose(3 x 14*x**2, x**2 x**3) 3 x**2 x**3 14*x**4 O(x**5) R.reversion(x 4*x**2 8*x**3) x - 4*x**2 24*x**3 - 160*x**4 O(x**5) R.differentiate(10 x x**2 x**3 x**4 R.order_term()) 1 2*x 3*x**2 4*x**3 O(x**4)注意各操作的前置条件源码中均有显式校验compose(s, t)t的常数项必须为 0否则ValueError: Series composition requires the constant term of the second series to be zeroinverse(s)s的常数项必须是底域中的单位NotReversible: Series inverse requires the constant term to be a unitreversion(s)常数项必须为 0、一次项系数必须是单位square(s)即multiply(s, s)。5.2 域级操作示例来自PowerSeriesRingField源码 docstring from sympy import QQ from sympy.polys.series import PowerSeriesRingField R PowerSeriesRingField(QQ, x) x R.gen R.sqrt(1 2*x 4*x**2 8*x**3) 1 x 3/2*x**2 5/2*x**3 - 29/8*x**4 - 1/8*x**5 O(x**6) R.integrate(1 x x**2 x**3 x**4 x**5 R.order_term()) x 1/2*x**2 1/3*x**3 1/4*x**4 1/5*x**5 O(x**6)六、超越函数级数log、exp 与三角展开原文档专门演示了有理数域上常见的数学函数展开这是幂级数环最实用的场景之一 R, x power_series_ring(x, QQ, 8) R.log(1 x) x - 1/2*x**2 1/3*x**3 - 1/4*x**4 1/5*x**5 - 1/6*x**6 1/7*x**7 O(x**8) R.tan(x) x 1/3*x**3 2/15*x**5 17/315*x**7 O(x**8)所有超越函数都在 sympy/polys/series/ringpython.py 中以_useries_*前缀实现其输入前置条件非常统一见各函数源码的ValueError分支函数输入要求实现策略从源码可读log(s)常数项必须为 1项数 20 时用log(f)f/f再积分否则直接复合标准级数ln(1u)log1p(s)常数项必须为 0先加 1 再调logexp(s)常数项必须为 0项数 20 时用 Newton 迭代校正exp (s - log(exp))*exp否则复合e^u级数expm1(s)常数项必须为 0先exp再减 1sin/cos常数项必须为 0项数大时借助tan半角公式否则复合标准麦克劳林级数tan/tanh常数项必须为 0Newton 迭代y_{n1} y_n (s − atan(y_n))·(1y_n²)asin/atan等常数项必须为 0利用导数关系如atan(f)f/(1f²)积分或复合级数以原文档的R.log(1 x)为例log要求常数项为 11x恰好满足输入仅 2 项走复合标准级数分支——先构造系数[0, 1, -1/2, 1/3, ...]_useries_log中的c列表再与(1x)-1 x复合得到x - x²/2 x³/3 - ...ringpython.py。而tan则通过_useries_tan的 Newton 迭代逐步逼近到目标精度ringpython.py。这些函数同样支持级数作为自变量例如工厂函数 docstring 中的示例 from sympy import QQ from sympy.polys.series import power_series_ring R, x power_series_ring(x, QQ) R.sin(x x**2) x x**2 - 1/6*x**3 - 1/2*x**4 - 59/120*x**5 O(x**6) R.log(1 7*x**2) 7*x**2 - 49/2*x**4 O(x**6)七、元素类PowerSeriesElement运算符重载与转换PowerSeriesElementsympy/polys/series/ring.py是环中元素的封装类继承DomainElement实现了完整的 Python 运算符重载加减乘除-*/含与int或底域元素混合运算的__radd__/__rsub__/__rmul__/__rtruediv__幂**支持整数次幂n0返回环的onen1返回自身否则走_pow_int一元、-比较基于表示比较equal_repr不同环的元素相加/相乘会抛ValueError: Cannot add/multiply power series with different rings.。元素还提供以下实用方法/属性 s.as_expr() # 转回 SymPy 表达式含 O(x**prec) s.removeO() # 去掉大 O 记号返回纯多项式部分 s.prec # 该级数的精度int 或 None s.constant_coefficient() # 常数项 s.is_ground # 是否退化为常数底域元素环侧配套的转换方法包括from_expr(expr)将 SymPy 表达式嵌入环若表达式含Order项则按其精度截断见 ring.py 的_rebuild_expr、from_list、to_list、to_dense、to_expr。ring_new即R(...)调用协议可接受PowerSeriesElement、Expr、int或底域元素四种输入。八、双后端架构纯 Python 与 FLINT 的自动切换原文档明确说明系统含两个后端纯 Python 实现与基于FLINT库的高性能实现前者是回退方案后者在可用时被自动选用。这一机制在 sympy/polys/series/tring.py 中实现模块根据sympy.external.gmpy.GROUND_TYPES判断当前 SymPy 使用的整数/有理数底层类型gmpy或flinttring.py若GROUND_TYPES flint导入 sympy/polys/series/ringflint.py 中的FlintPowerSeriesRingZZ/FlintPowerSeriesRingQQ底层级数类型为ZZSeries fmpz_series | fmpz_poly、QQSeries fmpq_series | fmpq_poly否则将 FLINT 类置为None回退到 sympy/polys/series/ringpython.py 的PythonPowerSeriesRingZZ/PythonPowerSeriesRingQQ底层级数类型为USeries tuple[list[Er], int | None]系数列表 精度标记。PowerSeriesRingZZ(prec)/PowerSeriesRingQQ(prec)两个工厂函数执行具体选择tring.py再由_power_series_ring(K, prec)根据底域分派到 ZZ 或 QQ 版本最终由高层power_series_ring包装为带生成元的结果。整个链路对用户透明。两个后端都实现了 sympy/polys/series/base.py 中PowerSeriesRingProto/PowerSeriesRingFieldProto协议定义的统一接口add/subtract/multiply/divide/pow_int/square/truncate/compose/inverse/reversion/differentiate域级再加integrate/sqrt/log/exp/tan/sin/cos等因此上层PowerSeriesElement无需关心具体后端。关于 FLINT 实现的两个工程细节ringflint.py_get_series_precision通过级数对象的prec属性获取精度对 python-flint ≤ 0.7.1 的旧版本退而解析repr字符串中的prec字段代码注释注明该兼容分支可在 python-flint 0.8 成为最低版本后移除_global_cap上下文管理器临时调整ctx.capFLINT 全局级数截断上限退出时恢复保证并发/嵌套环境下的精度隔离。此外 ringflint.py 的类声明带doctest_depends_on(ground_types[flint])装饰器表明其文档测试仅在 FLINT 后端下运行这从侧面印证了双后端的运行时条件依赖。九、与多项式域框架的集成作为 SymPy Domain 使用除独立使用的PowerSeriesRingRing之外模块还提供了领域domain层的集成类。原文档Domain for Univariate Power Series Rings一节指向 sympy/polys/domains/powerseriesring.pySeriesRingProto定义领域层接口的 Protocolring、dtype、domain、symbol、prec等属性以及from_expr/to_expr/domain_new/ring_new等方法PowerSeriesRing继承Ring与CompositeDomain的领域类其构造函数内部直接调用power_series_ring(str(symbol), domain, prec)powerseriesring.py并提供new、of_type、zero、one、prec、is_unit、canonical_unit、to_sympy、from_sympy等标准域方法使幂级数环可以嵌入 SymPy 的Domain体系例如供其他多项式域工具使用。__repr__形如ZZ[[x], 6]即底域 双括号 变量 精度的数学记号。类源码注释明确标注该领域类仍处于实验阶段experimental stage后续应补充与现有领域框架的完整集成使用时应注意这一前提。十、类型协议与测试验证10.1 协议层base.pysympy/polys/series/base.py 定义了三个核心协议/类型PowerSeriesRingProto[TSeries, Er]环级操作协议ZZ 后端PowerSeriesRingFieldProto[TSeries, Ef]域级操作协议继承前者并追加integrate/sqrt/log/exp及全部三角函数base.pyTSeries级数类型的 TypeVar而 sympy/polys/series/tring.py 中的TSeriesElement是仅供静态类型检查使用的底层级数占位类型。此外series_pprint(series, prec, symx, ascendingTrue)是统一的字符串渲染函数prec0时直接输出O(x**0)系数为空或全零时视prec输出O(x**prec)或0否则输出poly O(x**prec)base.py。这就是原文档示例中O(x**8)打印格式的来源。10.2 测试覆盖测试集中在 sympy/polys/series/tests/test_ring.py另有 test_hypothesis.py 用于基于假设的性质测试。从测试命名可见覆盖范围test_equal精度边界下的相等性、test_basicsone/zero/gen/prec/精度截断如R.series_prec(R([1,2,3,4,5,6,7])) 6、test_positive/test_negative正负号与精度收束、test_add/test_int_subtract加减法、精度统一、以及后续大量针对乘法、幂、逆、复合、反演、导数和超越函数的用例。测试通过参数化夹具ring_int/ring_rational同时跑纯 Python 与 FLINT 两种后端保证行为一致。十一、实战小结与注意事项回顾全文使用 SymPy 幂级数环的关键要点用工厂函数创建R, x power_series_ring(x, ZZ|QQ, prec)默认精度 6prec必须非负底域决定能力ZZ环只支持代数运算与求导、逆、复合、反演QQ域额外支持积分、开方与全部超越函数级数精度语义要牢记次数低于prec的结果是精确多项式打印无O触及prec自动截断并带O(x**prec)混合精度按取较小精度合并前置条件要满足inverse要求常数项为单位、reversion要求常数项为 0 且一次项为单位、compose要求第二参数常数项为 0、log要求常数项为 1、exp/sin/cos等要求常数项为 0违反时抛出带明确信息的异常后端自动选择安装 FLINT 支持GROUND_TYPES flint时自动启用高性能后端否则回退纯 Python接口完全一致转换与集成from_expr/to_expr/as_expr/removeO实现与 SymPy 表达式的双向互转PowerSeriesRing领域类可融入Domain框架但后者处于实验阶段。深入阅读建议本指南依据 doc/src/modules/polys/series.rst 撰写实现与验证请继续查看 sympy/polys/series/ring.py、sympy/polys/series/ringpython.py、sympy/polys/series/ringflint.py、sympy/polys/series/tring.py、sympy/polys/series/base.py、sympy/polys/domains/powerseriesring.py 与测试文件 sympy/polys/series/tests/test_ring.py。【免费下载链接】sympyA computer algebra system written in pure Python项目地址: https://gitcode.com/GitHub_Trending/sy/sympy创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考