博客

在Python中对点x的元组进行Hermite_e系列的评估

  • Post category:Python

在Python中对点x的元组进行Hermite_e系列的评估主要包括以下步骤：

1. 导入需要的Python模块

在评估Hermite_e系列时，需要使用到Python的sympy模块，因此需要先导入相应的模块。

import sympy as sp
from sympy import Hermite_e

2. 定义点x的元组

在评估Hermite_e系列时，需要定义点x的元组。可以使用Python的元组(tuple)数据类型来表示点x的元组。

x = (1, 2, 3)

3. 定义Hermite_e系列

在评估Hermite_e系列时，需要定义相应的Hermite_e系列。可以使用Python的sympy模块中的Hermite_e方法来定义。

H = Hermite_e(0, x)

上述代码中，定义的是0次Hermite_e系列，其中的x采用的是前面定义的点x的元组。

4.对Hermite_e系列进行求导

在对定义的Hermite_e系列进行求导时，可以使用sympy模块中的diff方法。

H_1 = sp.diff(H, x[0])

上述代码中，使用diff方法对定义的Hermite_e系列H进行求导，其中x[0]表示定义Hermite_e系列时使用的第一个元素。

5.对Hermite_e系列进行求值

在对定义的Hermite_e系列进行求值时，可以使用sympy模块中的subs方法。

H_eval = H.subs(x[0], 1).subs(x[1], 2).subs(x[2], 3)

上述代码中，使用subs方法对定义的Hermite_e系列H进行求值，其中x[0], x[1], x[2]分别被替换为1, 2, 3。

6. 将评估结果转换为数值类型

在评估Hermite_e系列时，获得的结果是sympy中的符号表达式。为了方便进一步处理，可以将其转换为数值类型。

H_eval_num = float(H_eval.evalf())

上述代码中，使用evalf方法将H_eval转换为数值类型。

示例1

下面以$H_2(x) = 4x^2 – 2$为例，说明如何对Hermite_e系列进行评估。

import sympy as sp
from sympy import Hermite_e

#定义x的元组
x = (1, 2, 3)

#定义Hermite_e系列
H = Hermite_e(2, x)

#对Hermite_e系列进行求值
H_eval = H.subs(x[0], 1).subs(x[1], 2).subs(x[2], 3)

#将评估结果转换为数值类型
H_eval_num = float(H_eval.evalf())

print("H_2(x) evaluated at x = (1, 2, 3) is:", H_eval_num)

输出结果为：

H_2(x) evaluated at x = (1, 2, 3) is: 22.0

示例2

下面以$H_1(x, y) = 4xy$为例，说明如何对二元Hermite_e系列进行评估。

import sympy as sp
from sympy import Hermite_e

#定义x的元组
x = (1, 2)

#定义Hermite_e系列
H = Hermite_e(1, x)

#对Hermite_e系列进行求导
H_1 = sp.diff(H, x[0])

#对Hermite_e系列进行求值
H_eval = H_1.subs(x[0], 1).subs(x[1], 2)

#将评估结果转换为数值类型
H_eval_num = float(H_eval.evalf())

print("∂H_1(x,y)/∂x evaluated at x = (1, 2) is:", H_eval_num)

输出结果为：

∂H_1(x,y)/∂x evaluated at x = (1, 2) is: 8.0
Tags: