在Python中使用三维系数数组在点(x,y)处评估二维Legendre级数
要在点x、y处评估二维Legendre级数,请使用Python NumPy中的polynomial.legendre.legval2d()方法。该方法返回由x和y的对应值组成的点处的二维Legendre级数的值。
第一个参数是x,y。在点(x,y)处评估二维级数,其中x和y必须具有相同的形状。如果x或y是列表或元组,则首先将其转换为ndarray,否则将其保留不变,如果它不是ndarray,则将其视为标量。
第二个参数是c。系数数组,其顺序为多重度i,j的项的系数包含在c[i,j]中。如果c的维数大于2,则剩余的指数枚举多个系数集。
步骤
首先导入所需的库 −
import numpy as np
from numpy.polynomial import legendre as L
创建一个三维系数数组 –
c = np.arange(24).reshape(2,2,6)
显示数组 –
print("我们的数组...\n",c)
检查维数 –
print("\n我们数组的维度...\n",c.ndim)
获取数据类型 –
print("\n我们数组对象的数据类型...\n",c.dtype)
获取形状 –
print("\n我们数组对象的形状...\n",c.shape)
要在点x、y处评估二维Legendre级数,请使用Python NumPy中的polynomial.legendre.legval2d()方法 –
print("\n结果...\n",L.legval2d([1,2],[1,2],c))
样例
import numpy as np
from numpy.polynomial import legendre as L
# Create a 3D array of coefficients
c = np.arange(24).reshape(2,2,6)
# Display the array
print("我们的数组...\n",c)
# Check the Dimensions
print("\n我们数组的维度...\n",c.ndim)
# Get the Datatype
print("\n我们数组对象的数据类型...\n",c.dtype)
# Get the Shape
print("\n我们数组对象的形状...\n",c.shape)
# To evaluate a 2D Legendre series at points x, y, use the polynomial.legendre.legval2d() method in Python NumPy
print("\n结果...\n",L.legval2d([1,2],[1,2],c))
输出
我们的数组...
[[[ 0 1 2 3 4 5]
[ 6 7 8 9 10 11]]
[[12 13 14 15 16 17]
[18 19 20 21 22 23]]]
我们数组的维度...
3
我们数组对象的数据类型...
int64
我们数组对象的形状...
(2, 2, 6)
结果...
[[ 36. 108.]
[ 40. 117.]
[ 44. 126.]
[ 48. 135.]
[ 52. 144.]
[ 56. 153.]]