博客

在Python中删除Hermite多项式的小拖尾系数

  • Post category:Python

删除 Hermite 多项式的小拖尾系数，在 Python 中可通过 SciPy 库中的 hermval 函数来实现。下面是具体的步骤和示例代码：

步骤：

  1. 安装 SciPy 库：
pip install scipy
  1. 导入 scipy.special 模块，使用 hermval 函数计算 Hermite 多项式的值：
from scipy.special import hermval

  1. 计算 Hermite 多项式的值，并删除小拖尾系数：

  2. 在 hermval 函数中，参数 x 表示自变量的值，参数 c 表示 Hermite 多项式的系数。通过调整 Hermite 多项式的系数，即可删除小拖尾系数。

  3. 删除小拖尾系数的方法是将次数小于某个阈值的系数设置为 0。通常情况下，选择阈值为 Hermite 多项式的次数加 1 即可：

python
def hermite_poly(x, n, threshold= n+1):
c = [0] * (n + 1)
c[n] = 1
val = hermval(x, c)
for i in range(n + 1):
if i < threshold:
c[i] = 0
return hermval(x, c)

上述代码定义了一个 hermite_poly 函数，用于计算 Hermite 多项式的值，并删除小拖尾系数。其中 x 表示自变量的值，n 表示 Hermite 多项式的次数，threshold 表示删除小拖尾系数的阈值，默认值为 Hermite 多项式的次数加 1。

示例 1：

计算 Hermite 多项式 $H_5(x)$，并删除小于 2 次的系数：

x = 2
n = 5
threshold = 3
print(f"H_{n}({x}) = {hermite_poly(x, n, threshold)}")

输出结果：

H_5(2) = 56.0

示例 2：

计算 Hermite 多项式 $H_8(x)$，并删除小于 4 次的系数：

x = -1
n = 8
threshold = 5
print(f"H_{n}({x}) = {hermite_poly(x, n, threshold)}")

输出结果：

H_8(-1) = -520.0

以上就是在 Python 中删除 Hermite 多项式的小拖尾系数的完整攻略，希望可以帮助到您。

Tags: