scipy.linalg.

invhilbert#

scipy.linalg.invhilbert(n, exact=False)[源代码]#

计算阶数为 n 的希尔伯特矩阵的逆矩阵。

希尔伯特矩阵的逆矩阵中的条目是整数。当 n 大于 14 时，逆矩阵中的某些条目将超出 64 位整数的上限。exact 参数提供了处理这些大整数的两种选择。

参数:

nint: 希尔伯特矩阵的阶数。
exactbool, 可选: 如果为 False，则返回的数组的数据类型为 np.float64，并且该数组是逆矩阵的近似值。如果为 True，则该数组是精确的整数逆矩阵数组。为了在 n > 14 时表示精确的逆矩阵，返回的数组是长整数的对象数组。对于 n <= 14，精确的逆矩阵将作为数据类型为 np.int64 的数组返回。

返回:

invh(n, n) ndarray: 如果 exact 为 False，则数组的数据类型为 np.float64。如果 exact 为 True，则数据类型为 np.int64（对于 n <= 14）或 object（对于 n > 14）。在后一种情况下，数组中的对象将是长整数。

另请参阅

hilbert: 创建希尔伯特矩阵。

说明

在 0.10.0 版本中添加。

示例

>>> from scipy.linalg import invhilbert
>>> invhilbert(4)
array([[   16.,  -120.,   240.,  -140.],
       [ -120.,  1200., -2700.,  1680.],
       [  240., -2700.,  6480., -4200.],
       [ -140.,  1680., -4200.,  2800.]])
>>> invhilbert(4, exact=True)
array([[   16,  -120,   240,  -140],
       [ -120,  1200, -2700,  1680],
       [  240, -2700,  6480, -4200],
       [ -140,  1680, -4200,  2800]], dtype=int64)
>>> invhilbert(16)[7,7]
4.2475099528537506e+19
>>> invhilbert(16, exact=True)[7,7]
42475099528537378560