在Python中沿特定轴差分一个具有多维系数的Chebyshev级数
要对Chebyshev级数进行差分,请使用Python Numpy中的polynomial.chebder()方法。该方法返回导数的Chebyshev级数。返回沿轴差分m次的Chebyshev级数系数c。在每次迭代中,结果都乘以scl。参数c是一个沿着每个轴从低到高次数的系数数组,例如,[1,2,3]表示系列1T_0 + 2T_1 + 3T_2,而 [[1,2],[1,2]] 表示 1T_0(x)T_0(y) + 1T_1(x)T_0(y) + 2T_0(x)T_1(y) + 2T_1(x)*T_1(y) 如果轴=0是x,轴=1是y。
第1个参数是c,Chebyshev系数的数组。如果c是多维的,则不同的轴对应于不同的变量,每个轴上的度数由相应的索引给出。第2个参数是m,需要取的导数的数量,必须是非负数。(默认值:1)。第3个参数是scl,即每个导数都乘以scl。最终结果是乘以scl**m。这是用于线性变量的变化。(默认值:1)。第4个参数是轴,即沿轴差分的轴。(默认值:0)。
步骤
首先,导入所需的库:
import numpy as np
from numpy.polynomial import chebyshev as C
创建一个Chebyshev级数系数的多维数组:
c = np.arange(4).reshape(2,2)
显示系数数组:
print("我们的系数数组...\n",c)
检查维度:
print("\n我们的Array的维度...\n",c.ndim)
获取数据类型:
print("\n我们Array对象的数据类型...\n",c.dtype)
获取形状:
print("\n我们Array对象的形状...\n",c.shape)
要对Chebyshev级数进行差分,请使用Python Numpy中的polynomial.chebder()方法:
print("\n结果...\n",C.chebder(c, axis = 1))
示例
import numpy as np
from numpy.polynomial import chebyshev as C
# 创建一个Chebyshev级数系数的多维数组
c = np.arange(4).reshape(2,2)
# 显示系数数组
print("我们的系数数组...\n",c)
# 检查维度
print("\n我们的Array的维度...\n",c.ndim)
# 获取数据类型
print("\n我们Array对象的数据类型...\n",c.dtype)
# 获取形状
print("\n我们Array对象的形状...\n",c.shape)
# 要对Chebyshev级数进行差分,请使用Python Numpy中的polynomial.chebder()方法。
# 该方法返回导数的Chebyshev级数。
print("\n结果...\n",C.chebder(c, axis = 1))
输出
我们的系数数组...
[[0 1]
[2 3]]
我们的Array的维度...
2
我们Array对象的数据类型...
int64
我们Array对象的形状...
(2, 2)
结果...
[[1.]
[3.]]