2026a

# cauchy


柯西矩阵

函数库: TyMath

# 语法

A = cauchy([T],x,y)
A = cauchy([T],x)

# 说明

A = cauchy([T],x,y) 返回一个 n×n 矩阵,其中包含的项满足 A[i,j] = 1/(x[i]+y[j])。参数 x 和 y 是长度为 n 的向量,T 为可选返回类型参数。如果您为 x 和 y 传入标量,它们将解释为向量 1:x 和 1:y。示例


A = cauchy([T], x) 返回与上述 y = x 的情况相同的结果。换言之,A[i,j] = 1/(x[i]+x[j])。示例

# 示例

柯西矩阵的生成及其性质

对给定的恒正、递增、不含有共有元素的 X 和 Y 生成对应的 cauchy 矩阵。

using TyMath
X = [1,3,5,6,7]
Y = [2,6,7,8,10]
C = TestArrays.cauchy(X,Y)
C = 5×5 Matrix{Float64}:

 0.333333  0.142857   0.125      0.111111   0.0909091
 0.2       0.111111   0.1        0.0909091  0.0769231
 0.142857  0.0909091  0.0833333  0.0769231  0.0666667
 0.125     0.0833333  0.0769231  0.0714286  0.0625
 0.111111  0.0769231  0.0714286  0.0666667  0.0588235

由于 X 和 Y 含有共有元素, det(C) 为零。

det(C)
ans = 2.357483478207392e-17

由于 X 恒正、递增,使用 X 生成的 cauchy 矩阵为对称正定阵。

C = TestArrays.cauchy(X)
eigenC = eigen(C).values
C = 5×5 Matrix{Float64}:

 0.5       0.25      0.166667   0.142857   0.125
 0.25      0.166667  0.125      0.111111   0.1
 0.166667  0.125     0.1        0.0909091  0.0833333
 0.142857  0.111111  0.0909091  0.0833333  0.0769231
 0.125     0.1       0.0833333  0.0769231  0.0714286

eigenC = 5-element Vector{Float64}:

 1.460969227992038e-7
 5.045242222152628e-5
 0.0037801834568022743
 0.10550546622236223
 0.8120923232302627

对于任何非负数 r(这里取为 1.5),C.^r 为对称半正定阵。

r = 1.5
Cr = C.^r
eigenCr = eigen(Cr).values
eigenCr = 5-element Vector{Float64}:

 1.012860912605185e-7
 3.1917855616886315e-5
 0.0024246625640654954
 0.05713628615052903
 0.43677093100689623

# 输入参数

x - 输入向量
标量 | 向量

输入向量,指定为向量或标量。对于标量输入,将会被解释为 1:x。

y - 输入向量
标量 | 向量

输入向量,指定为向量或标量。对于标量输入,将会被解释为 1:y。

T - 输出类型
类型

输出类型,指定为类型(Type)。该参数为可选参数,如果不指定,会依据输入的默认值。如果可能,函数将会输出 T 指定类型的矩阵。

# 输出参数

A - 输出矩阵
矩阵

输出矩阵,以矩形形式返回。

# 参考

[1] R. Bhatia, Infinitely divisible matrices, Amer. Math. Monthly, (2005), to appear.

[2] D. E. Knuth, The Art of Computer Programming, Volume 1, Fundamental Algorithms, third edition, Addison-Wesley, Reading,Massachusetts, 1997.

[3] E. E. Tyrtyshnikov, Cauchy-Toeplitz matrices and some applications, Linear Algebra and Appl., 149 (1991), pp. 1-18.

[4] O. Taussky and M. Marcus, Eigenvalues of finite matrices, in Survey of Numerical Analysis, J. Todd, ed., McGraw-Hill, New York, 1962, pp. 279-313.

# 另请参阅

eigen