根据矩阵数字特征的定义,结合矩阵Lie积的定义,利用矩阵迹的性质,矩阵数字特征的性质,以及密度矩阵的性质,对由Hermite矩阵组成的实线性空间Hn上的一些数字特征进行了数值计算并且给出了一个例子,其次对WY相关系数进行改造使其成为矩阵空间上的半内积并且利用矩阵范数、矩阵迹及矩阵数字特征的性质得出了改造后的WY相关系数与原WY相关系数、斜信息以及协方差的关系,获得了原WY相关系数和斜信息的一些性质,并给出改造后的WY相关系数的一个应用。这些结果将在量子信息论中具有一定的应用。 According to the definition of matrix numerical characteristics, combined with the definitions of Lie product of matrix, using the properties of matrix traces, matrix numerical characteristics, and density matrix, numerical calculations were performed on some numerical characteristics on real linear space Hncomposed of Hermite matrices, and an example was given, secondly, the WY corre-lation coefficient was modified to become a semi inner product in the matrix space, and the rela-tionship between the modified WY correlation coefficient and the original WY correlation coefficient and skew information, as well as some properties of covariance, the original WY correlation coefficient and skew information, were obtained using the properties of matrix norm, matrix trace, and matrix numerical characteristics, and an application of the modified WY correlation coefficient was also presented. These results will have certain applications in quantum information theory.
根据矩阵数字特征的定义,结合矩阵Lie积的定义,利用矩阵迹的性质,矩阵数字特征的性质,以及密度矩阵的性质,对由Hermite矩阵组成的实线性空间Hn上的一些数字特征进行了数值计算并且给出了一个例子,其次对WY相关系数进行改造使其成为矩阵空间上的半内积并且利用矩阵范数、矩阵迹及矩阵数字特征的性质得出了改造后的WY相关系数与原WY相关系数、斜信息以及协方差的关系,获得了原WY相关系数和斜信息的一些性质,并给出改造后的WY相关系数的一个应用。这些结果将在量子信息论中具有一定的应用。
密度矩阵,矩阵期望,矩阵方差,矩阵协方差,WY相关系数,矩阵斜信息
Ouyang Yu, Fangguo Ren*
School of Mathematics and Statistics, Shaanxi Normal University, Xi’an Shaanxi
Received: Dec. 25th, 2023; accepted: Jan. 5th, 2024; published: Feb. 23rd, 2024
According to the definition of matrix numerical characteristics, combined with the definitions of Lie product of matrix, using the properties of matrix traces, matrix numerical characteristics, and density matrix, numerical calculations were performed on some numerical characteristics on real linear space Hncomposed of Hermite matrices, and an example was given, secondly, the WY correlation coefficient was modified to become a semi inner product in the matrix space, and the relationship between the modified WY correlation coefficient and the original WY correlation coefficient and skew information, as well as some properties of covariance, the original WY correlation coefficient and skew information, were obtained using the properties of matrix norm, matrix trace, and matrix numerical characteristics, and an application of the modified WY correlation coefficient was also presented. These results will have certain applications in quantum information theory.
Keywords:Density Matrix, Matrix Expectation, Matrix Variance, Matrix Covariance, WY Correlation Coefficient, Matrix Skew Information
Copyright © 2024 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
近年来量子信息的研究与应用受到国内外的广泛关注,而概率与矩阵是研究量子信息的重要的不可缺少的工具,如量子信息论中主要的研究对象(态、量子测量及量子运算)与概率及矩阵密不可分 [
为了简述方便,我们对文中符号进行说明:Mn表示所有n阶复矩阵组成的复线性空间,Hn表示所有n阶Hermite矩阵组成的实线性空间,In表示n阶单位矩阵, A * 和 T r ( A ) 分别表示矩阵A的共轭转置和迹;不特别说明,矩阵空间Mn上的默认标准内积为 ( X , Y ) = T r ( X * Y ) ,其中 X , Y ∈ M n 。其它未作说明的符号、概念及术语参见文献 [
定义1 [
定义2 [
E P ( X ) ≡ T r ( P X ) ;
V P ( X ) ≡ T r ( P ( X − E P ( X ) I n ) ( X − E P ( X ) I n ) * ) ;
C o v P ( X , Y ) ≡ T r ( P ( X − E P ( X ) I n ) ( Y − E P ( Y ) I n ) * ) ;
C o r r P ( X , Y ) ≡ T r ( P X Y * − P X P Y * ) ;
I P ( X ) ≡ T r ( P X X * − P X P X * ) .
分别称为X关于P的期望、X关于P的方差、X与Y关于P的协方差、WY相关系数及斜信息。
定义3 [
定义4 在矩阵空间Mn上定义二元函数 C o r r P G ( ⋅ , ⋅ ) ,其中
C o r r P G ( X , Y ) = 1 2 T r ( P X Y * + P Y * X ) − T r ( P X P Y * ) , ∀ X , Y ∈ M n .
定义5 ∀ A ∈ M n ,矩阵空间Mn上的Frobenius范数 ‖ A ‖ F 定义为 ‖ A ‖ F = [ T r ( A * A ) ] 1 2
引理1设 A , B ∈ M n 是Hermite矩阵。则
(1) T r ( A B ) 是实数;
(2) 如果 A , B 是半正定矩阵,则 T r ( A B ) 是非负实数,且 T r ( A B ) = 0 当且仅当 A B = 0 。
证明 (1) 由于A是Hermite矩阵,则存在酉矩阵U,使得 U * A U = Λ = d i a g ( λ 1 , ⋯ , λ n ) ,其中 λ i ∈ R 。再由B是Hermite矩阵知, U * B U = ( b i j ) n 的对角元 b i i ∈ R ,则 T r ( A B ) = T r ( U * A U U * B U ) = T r ( Λ U * B U ) = ∑ i = 1 n λ i b i i ∈ R 。
(2) T r ( A B ) = T r ( A A B B ) = T r ( A B B A ) = T r ( ( A B ) ( A B ) * ) = ‖ A B ‖ 2 2 ,则 T r ( A B ) ≥ 0 ,且 T r ( A B ) = 0 ,当且仅当 A B = 0 ,当且仅当 A B = 0 。
根据矩阵迹的性质及引理1易得:
引理2 给定半正定矩阵 P ∈ M n , M n 上的二元函数 f i ( i = 1 , ⋯ , 6 ) 分别为 M n 上的半内积,定义如下: ∀ X , Y ∈ M n
(1) f 1 ( X , Y ) = T r ( P X Y * ) ;
(2) f 2 ( X , Y ) = T r ( P Y * X ) ;
(3) f 3 ( X , Y ) = T r ( P X ) T r ( P Y * ) ;
(4) f 4 ( X , Y ) = T r ( P X P Y * ) ;
(5) f 5 ( X , Y ) = T r ( P X Y * + P X P Y * ) ;
(6) f 6 ( X , Y ) = T r ( P Y * X + P X P Y * ) 。
引理3 [
(1) | f ( x , y ) | 2 ≤ f ( x , x ) f ( y , y ) , ∀ x , y ∈ V ;
(2) 由f诱导的V上的一元实函数 f ( x , x ) , ∀ x ∈ V 是V上的一个半范数。
引理4 设 P ∈ M n 是密度矩阵, X , Y , Z ∈ M n , λ , μ ∈ C , U ∈ M n 是酉矩阵,则
(1) C o v P ( X , Y ) = T r ( P X Y * ) − T r ( P X ) T r ( P Y ) ¯ ;
(2) C o v P ( X , Y ) ¯ = C o v P ( Y , X ) ;
(3) C o v P ( λ X + μ Y , Z ) = λ C o v P ( X , Z ) + μ C o v P ( Y , Z ) ;
(4) V P ( X ) ≥ 0 ;
(5) C o r r P ( X , Y ) ¯ = C o r r P ( Y , X ) ;
(6) C o r r P ( λ X + μ Y , Z ) = λ C o r r P ( X , Z ) + μ C o r r P ( Y , Z ) ;
(7) (i) C o v P ( X , Y ) = C o v U * P U ( U * X U , U * Y U ) ;
(ii) C o r r P ( X , Y ) = C o r r U * P U ( U * X U , U * Y U ) ;
(8) C o v P ( ⋅ , ⋅ ) 是矩阵空间 M n 上的一个半内积;
(9) V P ( ⋅ ) 是矩阵空间 M n 上的一个半范数;
(10) V U * P U ( U * X U ) = V P ( X ) , I U * P U ( U * X U ) = I P ( X ) 。
证明 由定义2且通过直接计算可知(2),(3),(5),(6),(7)成立,下证(1),(4),(8),(9),(10)。
(1) C o v P ( X , Y ) = T r ( P ( X Y * − E P ( Y ) ¯ X − E P ( X ) Y * + E P ( X ) E P ( Y ) ¯ I n ) ) = T r ( P X Y * ) − E P ( Y ) ¯ T r ( P X ) − E P ( X ) T r ( P Y * ) + E P ( X ) E P ( Y ) ¯ T r ( P ) = T r ( P X Y * ) − E P ( Y ) ¯ T r ( P X ) − E P ( X ) T r ( Y P ) ¯ + E P ( X ) E P ( Y ) ¯ = T r ( P X Y * ) − E P ( X ) E P ( Y ) ¯ = T r ( P X Y * ) − T r ( P X ) T r ( P Y ) ¯
(4) 由 V P ( ⋅ ) 的定义及(1)知, V P ( X ) = T r ( P X X * ) − | T r ( P X ) | 2 。再由柯西不等式知, | T r ( P X ) | 2 = | T r ( ( P X ) P ) | 2 ≤ T r ( ( P X ) ( P X ) * ) T r ( P P ) = T r ( P X X * P ) T r ( P ) = T r ( P X X * ) ,所以 V P ( X ) ≥ 0 。
(8) 由(2),(3),(4)可知,(8)成立。
(9) 由(8)及引理3可知,(9)成立。
(10) 由(7)可知,(10)成立。
引理5 设 P ∈ M n 是密度矩阵, A , B ∈ H n 。则 I m C o v P ( A , B ) = 1 2 i T r ( P [ A , B ] ) 。
证明 由引理4 (1)可得
2 i I m C o v P ( A , B ) = C o v P ( A , B ) − C o v P ( A , B ) ¯ = T r ( P ( A B − B A ) ) = T r ( P [ A , B ] ) ,所以 I m C o v P ( A , B ) = 1 2 i T r ( P [ A , B ] ) 。
讨论在实线性空间Hn上数字特征的数值计算。由引理4 (7),(10)可知,下面讨论的数字特征具有酉相似不变性,再由于任何一个Hermite矩阵都酉相似于一个对角矩阵,因此将密度矩阵取为非负对角矩阵。下面结合矩阵Lie积的定义,利用迹的循环性(即 T r ( X Y ) = T r ( Y X ) ,其中X是 m × n 矩阵,Y是 n × m 矩阵),密度矩阵P的性质(即 T r ( P ) = 1 ),矩阵数字特征的性质(引理4),对由Hermite矩阵组成的实线性空间Hn上的一些数字特征进行数值计算。
定理1 设 P = d i a g ( λ 1 , ⋯ , λ n ) ∈ M n 是密度矩阵, A = ( a i j ) , B = ( b i j ) ∈ H n 。则
(1) E P ( A ) = ∑ i = 1 n λ i a i i ,且 min { a i i } ≤ E P ( A ) ≤ max { a i i } ;
(2) V P ( A ) = ∑ i = 1 n ∑ j = 1 n λ i | a i j | 2 − ( ∑ i = 1 n λ i a i i ) 2 ;
(3) C o v P ( A , B ) = ∑ i = 1 n ∑ j = 1 n λ i a i j b i j ¯ − ( ∑ i = 1 n λ i a i i ) ( ∑ i = 1 n λ i b i i ) ;
(4) T r ( P [ A , B ] ) = 2 i ∑ 1 ≤ k < j ≤ n n ( λ k − λ j ) I m ( a k j b k j ¯ ) ;
(5) I P ( A ) = ∑ 1 ≤ k < j ≤ n n ( λ i + λ j − 2 λ i λ j ) | a i j | 2 ;
(6) C o r r P ( A , B ) = ∑ 1 ≤ k < j ≤ n n ( ( λ k + λ j − 2 λ k λ j ) R e ( a k j b k j ¯ ) + i ( λ k − λ j ) I m ( a k j b k j ¯ ) ) 。
证明 (1) 由定义易知, E P ( A ) = ∑ i = 1 n λ i a i i 。由引理1可知 E P ( A ) 是实数,又由于
T r P = 1 ,所以 ∑ i = 1 n λ i = 1 ,故 E P ( A ) = ∑ i = 1 n λ i a i i 是 a i i 的一个凸组合 i = 1 , 2 , ⋯ , n ,因此
min { a i i } ≤ E P ( A ) ≤ max { a i i } 。
(2) V P ( A ) = T r ( P A 2 ) − ( T r ( P A ) ) 2 = ∑ i = 1 n ∑ j = 1 n λ i | a i j | 2 − ( ∑ i = 1 n λ i a i i ) 2 。
(3) 由于 T r ( P A B ) = ∑ i = 1 n ∑ j = 1 n λ i a i j b i j ¯ ,则 C o v P ( A , B ) = T r ( P A B ) − T r ( P A ) T r ( P B ) = ∑ i = 1 n ∑ j = 1 n λ i a i j b i j ¯ − ( ∑ i = 1 n λ i a i i ) ( ∑ i = 1 n λ i b i i ) 。
(4) T r ( P [ A , B ] ) = T r ( P A B − P B A ) = ∑ k = 1 n ∑ j = 1 n λ k a k j b k j ¯ − ∑ k = 1 n ∑ j = 1 n λ j b j k a j k ¯ = ∑ 1 ≤ k < j ≤ n n ( λ k − λ j ) ( a k j b k j ¯ − b k j a k j ¯ ) = ∑ 1 ≤ k < j ≤ n n ( ( λ k − λ j ) ( R e ( a k j b k j ¯ ) + i I m ( a k j b k j ¯ ) − R e ( a k j ¯ b k j ) − i I m ( a k j ¯ b k j ) ) ) = 2 i ∑ 1 ≤ k < j ≤ n n ( λ k − λ j ) I m ( a k j b k j ¯ ) .
(5) I P ( A ) = T r ( P A 2 − P A P A ) = ∑ i = 1 n ∑ j = 1 n λ i | a i j | 2 − ∑ i = 1 n ∑ j = 1 n λ i λ j | a i j | 2 = ∑ 1 ≤ k < j ≤ n n ( λ i + λ j − 2 λ i λ j ) | a i j | 2 .
(6) C o r r P ( A , B ) = T r ( P A B − P A P B ) = ∑ k = 1 n ∑ j = 1 n λ k a k j b k j ¯ − ∑ k = 1 n ∑ j = 1 n λ k λ j a k j b k j ¯ = ∑ k ≠ j n ( λ k − λ k λ j ) a k j b k j ¯ = ∑ 1 ≤ k < j ≤ n n ( ( λ k − λ k λ j ) a k j b k j ¯ + ( λ j − λ k λ j ) a j k b k j ¯ ) = ∑ 1 ≤ k < j ≤ n n ( ( λ k − λ k λ j ) ( R e ( a k j b k j ¯ ) + i I m ( a k j b k j ¯ ) ) + ( λ j − λ k λ j ) ( R e ( a k j b k j ¯ ) + i I m ( a k j ¯ b k j ) ) ) = ∑ 1 ≤ k < j ≤ n n ( ( λ k + λ j − 2 λ k λ j ) R e ( a k j b k j ¯ ) + i ( λ k − λ j ) I m ( a k j b k j ¯ ) )
推论1 设 P = d i a g ( λ 1 , λ 2 ) ∈ M 2 是密度矩阵, A = ( a i j ) , B = ( b i j ) ∈ H 2 , X = ( x i j ) ∈ M 2 。则
(1) | T r ( P [ A , B ] ) | 2 = 4 ( λ 1 − λ 2 ) 2 | I m ( a 12 b 12 ¯ ) | 2 ;
(2) I P ( A ) = ( λ 1 + λ 2 − 2 λ 1 λ 2 ) | a 12 | 2 = ( 1 − 2 λ 1 λ 2 ) | a 12 | 2 ;
(3) C o r r P ( A , B ) = ( λ 1 + λ 2 − 2 λ 1 λ 2 ) R e ( a 12 b 12 ¯ ) + i ( λ 1 − λ 2 ) I m ( a 12 b 12 ¯ ) ;
(4) I P ( X ) = T r ( P X X * − P X P X * ) = ∑ i = 1 2 ∑ j = 1 2 λ i | x i j | 2 − ∑ i = 1 2 ∑ j = 1 2 λ i λ j | x i j | 2 。
定理1对Hn上的一些数字特征进行了计算,并且给出了具体的计算结果。由于Mn中的任何一个矩阵都可以写成一个Hermite矩阵与一个反Hermite矩阵的和,所以研究Hn上的数字特征对Mn上的数字特征的研究有着很大的价值,故定理1中对Hn上数字特征的计算结果将有助于对Mn上数字特征的研究。
下面给出一个计算H2上数字特征的例子。
例1 设 P = ( 1 3 0 0 2 3 ) , A = ( 0 1 1 0 ) , B = ( 0 − i i 0 ) 。则
(1) 由 E P ( A ) = ∑ i = 1 2 λ i a i i 知, E P ( A ) = 0 , E P ( B ) = 0 ;
(2) 由 T r ( P A B ) = ∑ i = 1 2 ∑ j = 1 2 λ i a i j b i j ¯ 知, C o v P ( A , B ) = i 1 3 − i 2 3 = − i 1 3 ;
(3) 由 V P ( A ) = ∑ i = 1 n ∑ j = 1 n λ i | a i j | 2 − ( ∑ i = 1 n λ i a i i ) 2 知, V P ( A ) = 1 , V P ( B ) = 1 ;
(4) 由 C o r r P ( A , B ) = ( λ 1 + λ 2 − 2 λ 1 λ 2 ) R e ( a 12 b 12 ¯ ) + i ( λ 1 − λ 2 ) I m ( a 12 b 12 ¯ ) 知, C o r r P ( A , B ) = i ( 1 3 − 2 3 ) I m ( i ) = − i 1 3 。由 I P ( A ) = ( λ 1 + λ 2 − 2 λ 1 λ 2 ) | a 12 | 2 知, I P ( A ) = 1 − 2 2 3 = I P ( B ) ,由 | T r ( P [ A , B ] ) | 2 = 4 ( λ 1 − λ 2 ) 2 | I m ( a 12 b 12 ¯ ) | 2 知, | T r ( P [ A , B ] ) | 2 = 4 9 ,故 V P ( A ) V P ( B ) = 1 ≥ 1 9 = | T r ( P [ A , B ] ) | 2 4 ,即 V P ( A ) V P ( B ) ≥ | T r ( P [ A , B ] ) | 2 4 ,这便是量子信息中Heisenberger不确定性关系的常见形式。
由例1可以看出矩阵数字特征与量子信息中的不确定性关系密切相关。
例2 设 P = ( 1 3 0 0 2 3 ) , X = ( 1 + i 2 + 4 i 0 2 + 3 i ) 。则
T r ( P X X * ) = T r ( ( 1 3 0 0 2 3 ) ( 1 + i 2 + 4 i 0 2 + 3 i ) ( 1 − i 0 2 − 4 i 2 − 3 i ) ) = 16
T r ( P X P X * ) = T r ( ( 1 3 0 0 2 3 ) ( 1 + i 2 + 4 i 0 2 + 3 i ) ( 1 3 0 0 2 3 ) ( 1 − i 0 2 − 4 i 2 − 3 i ) ) = 28 + 20 2 3
C o r r P ( X , X ) = I P ( X ) = T r ( P X X * − P X P X * ) = 16 − 28 + 20 2 3 = 20 ( 1 − 2 ) 3 < 0 。
由例2可以看到存在密度矩阵P,使得矩阵空间 M 2 上的WY相关系数 C o r r P ( ⋅ , ⋅ ) 不满足半内积的非负性,所以存在密度矩阵P使得 C o r r P ( ⋅ , ⋅ ) 不是矩阵空间 M 2 上的半内积。对任意的 X , Y ∈ M n 以及密度矩阵P,现在对WY相关系数 C o r r P ( X , Y ) ≡ T r ( P X Y * − P X P Y * ) 进行改造,将其变为 C o r r P G ( X , Y ) ≡ 1 2 T r ( P X Y * + P Y * X ) − T r ( P X P Y * ) 。
下面证明 C o r r P G ( ⋅ , ⋅ ) 是矩阵空间 M n 上的半内积并利用矩阵迹的性质(即 T r ( X Y ) = T r ( Y X ) ,其中X是 m × n 矩阵,Y是 n × m 矩阵)、矩阵范数的性质(即 ‖ A ‖ F ≥ 0 , ∀ A ∈ M n 且 ‖ A ‖ F = 0 当且仅当 A = 0 )、以及矩阵数字特征的性质(引理4),得出 C o r r P G ( ⋅ , ⋅ ) 以及 C o r r P ( ⋅ , ⋅ ) 、协方差、斜信息的一些性质并给出 C o r r P G ( ⋅ , ⋅ ) 的一个应用。
定理2设 P ∈ M n 是密度矩阵, A , B ∈ H n , X , Y ∈ M n 。则
(1) C o r r P G ( ⋅ , ⋅ ) 是矩阵空间 M n 上的一个半内积;
(2) 若X是正规矩阵,则 C o r r P G ( X , X ) = C o r r P ( X , X ) = I P ( X ) ≥ 0 ;
(3) 如果P是一个纯态,则 C o r r P ( X , Y ) = C o v P ( X , Y ) ;
(4) (i) C o r r P G ( A , B ) = R e C o r r P ( A , B ) = 1 4 ( I P ( A + B ) − I P ( A − B ) ) ;
(ii) I m C o r r P ( A , B ) = 1 2 i T r ( P [ A , B ] ) = I m C o v P ( A , B ) ;
(iii) 2 I P ( A ) = ‖ [ P , A ] ‖ F 2 ≥ 0 ,即当 A ∈ H n 时, I P ( A ) ≥ 0 ;
(iv) I P ( A ) = 0 当且仅当 A ∈ C ( P ) = { X ∈ H n | P X = X P } , d i m { A ∈ H n | I P ( A ) = 0 } ≥ n .
证明 (1) C o r r P G ( ⋅ , ⋅ ) 是矩阵空间 M n 上的一个半内积;根据前面引理2、引理3和引理4的论述,只需证明非负性即可。由柯西不等式及算术-几何不等式知,
T r ( P X P X * ) = | T r ( P X ( X P ) * ) | ≤ T r ( P X X * P ) T r ( X P P X * ) = T r ( P X X * ) T r ( P X * X ) ≤ 1 2 ( T r ( P X X * ) + T r ( P X * X ) ) 。
再由 C o r r P G ( I n , I n ) = 1 2 T r ( P + P ) − T r ( P ) = 0 知, C o r r P G ( ⋅ , ⋅ ) 是矩阵空间 M n 上的一个半内积。
(2) 若X是正规矩阵,则由(1)可知 C o r r P G ( X , X ) = C o r r P ( X , X ) = I P ( X ) ≥ 0 。
(3) 如果P是一个纯态,由 r ( P ) = 1 及 T r ( P ) = 1 知,P只有唯一的非零特征值1,再由P是半正定矩阵知,存在单位向量 α ∈ C n ,使得 P = α α * ,于是P是一个正交投影矩阵,有 P = P = P 2 ,则
T r ( P X P Y * ) = T r ( α α * X α α * Y * ) = T r ( ( α * X α ) ( α * Y * α ) ) = T r ( ( α * X α ) ( α * Y α ) * ) = ( α * X α ) ( α * Y α ) ¯ = T r ( α α * X ) T r ( α α * Y ) ¯ = T r ( P X ) T r ( P Y ) ¯ ,
所以 C o r r P ( X , Y ) = C o v P ( X , Y ) 。
(4) (i) 由 T r ( P ( A B + B A ) ) = 1 2 T r ( P ( ( A + B ) 2 − ( A − B ) 2 ) ) ,
C o r r P ( A , B ) + C o r r P ( A , B ) ¯ = T r ( P A B ) + T r ( P B A ) − T r ( P A P B ) − T r ( P B P A ) = T r ( P A B ) + T r ( P B A ) − 2 T r ( P A P B )
以及
T r ( P A P B ) = T r ( P ( A + B ) + ( A − B ) 2 P ( A + B ) − ( A − B ) 2 ) = 1 4 T r ( P ( A + B ) P ( A + B ) − P ( A + B ) P ( A − B ) + P ( A − B ) P ( A + B ) − P ( A − B ) P ( A − B ) ) = 1 4 T r ( P ( A + B ) P ( A + B ) − P ( A − B ) P ( A − B ) )
可知,
R e C o r r P ( A , B ) = 1 2 ( C o r r P ( A , B ) + C o r r P ( A , B ) ¯ ) = 1 2 T r ( P ( A B + B A ) ) − T r ( P A P B ) = 1 4 T r ( P ( A + B ) 2 − P ( A + B ) P ( A + B ) ) − 1 4 T r ( P ( A − B ) 2 − P ( A − B ) P ( A − B ) ) = 1 4 ( I P ( A + B ) − I P ( A − B ) )
(ii) 由 C o r r P ( A , B ) − C o r r P ( A , B ) ¯ = T r ( P A B − P A P B ) − T r ( P B A − P B P A ) = T r ( P A B ) − T r ( P B A ) + T r ( P B P A ) − T r ( P A P B ) = T r ( P A B ) − T r ( P B A )
及引理5可知
I m C o r r p ( A , B ) = 1 2 i ( C o r r P ( A , B ) − C o r r P ( A , B ) ¯ ) = 1 2 i T r ( P A B − P B A ) = 1 2 i T r ( P [ A , B ] ) = I m C o v P ( A , B ) .
(iii) 2 I P ( A ) = ‖ [ P , A ] ‖ F 2 ,因为 I P ( A ) = T r ( P A 2 − P A P A ) = T r ( P A A P − P A P A ) = − T r ( P A ( P A − A P ) ) = − 1 2 T r ( [ P A − A P ] 2 ) = 1 2 ‖ [ P , A ] ‖ F 2 ,所以 I P ( A ) ≥ 0 。
(iv) 由 2 I P ( A ) = ‖ [ P , A ] ‖ F 2 知, I P ( A ) = 0 当且仅当 ‖ [ P , A ] ‖ F 2 = 0 ,当且仅当 [ P , A ] = 0 ,进而有 P A = P ( P A ) = P ( A P ) = ( P A ) P = ( A P ) P = A P ;
如果 P A = A P ,于是由 P 是P的多项式可知, P A = A P ,故 I P ( A ) = 0 。由于 A ∈ C ( P ) ,P是半正定矩阵,则存在酉矩阵U,使得 U * P U = ( λ 1 I k 1 ⋱ λ m I k m ) ,其中 λ i ( i = 1 , ⋯ , m ) 非负互不相同且 ∑ i = 1 m k i = n ,于是 C ( U * P U ) = { ( X k 1 ⋱ X k m ) | X k i ∈ H k i , i = 1 , ⋯ , m } ,再由 H n 作为实线性空间是 n 2 维的知, d i m C ( U * P U ) = ∑ i = 1 m k i 2 ,故由 I P ( ⋅ ) 具有酉相似不变性及酉相似保持可交换性知, d i m C ( P ) = d i m C ( U * P U ) ≥ n 。
下面给出改造后的WY相关系数 C o r r P G ( ⋅ , ⋅ ) 的一个应用。
推论2 设 P ∈ M n 是密度矩阵, A , B ∈ H n ,则
1 16 ( I P ( A + B ) − I P ( A − B ) ) 2 ≤ I P ( A ) I P ( B ) ≤ 1 16 ( I P ( A + B ) + I P ( A − B ) ) 2
证明 因为 A , B ∈ H n ,所以A和B是正规矩阵,故由定理2 (2)知, C o r r P G ( A , A ) = I P ( A ) 且 C o r r P G ( B , B ) = I P ( B ) ,由定理2 (4)的(i)和(iii)知, C o r r P G ( A , B ) = 1 4 ( I P ( A + B ) − I P ( A − B ) ) 且 I P ( A + B ) 和 I P ( A − B ) 是实数,又由于 C o r r P G ( ⋅ , ⋅ ) 是矩阵空间 M n 上的一个半内积,则由引理3 (1)知,
C o r r P G ( A , B ) ≤ | C o r r P G ( A , B ) | ≤ C o r r P G ( A , A ) C o r r P G ( B , B ) = I P ( A ) I P ( B ) 。即
1 16 ( I P ( A + B ) − I P ( A − B ) ) 2 ≤ I P ( A ) I P ( B ) 。
另一方面,由于 I P ( ⋅ ) 是由矩阵空间 M n 上的一个半内积 C o r r P G ( ⋅ , ⋅ ) 诱导出来的半范数,满足平行四边形法则,于是 2 I P ( A ) I P ( B ) ≤ I P ( A ) + I P ( B ) = 1 2 ( I P ( A + B ) + I P ( A − B ) ) ,即有 I P ( A ) I P ( B ) ≤ 1 16 ( I P ( A + B ) + I P ( A − B ) ) 2 ,所以综上,有 1 16 ( I P ( A + B ) − I P ( A − B ) ) 2 ≤ I P ( A ) I P ( B ) ≤ 1 16 ( I P ( A + B ) + I P ( A − B ) ) 2
注意:不一定有 I P ( A ) I P ( B ) ≥ | T r ( P [ A , B ] ) | 2 4 成立。
例3 设 P = ( 1 5 0 0 4 5 ) , A = ( 0 1 1 0 ) , B = ( 0 − i i 0 ) 。 I P ( A ) I P ( B ) ≥ | T r ( P [ A , B ] ) | 2 4 不成立。
由 I P ( A ) = ( 1 − 2 λ 1 λ 2 ) | a 12 | 2 知, I P ( A ) = 1 − 4 5 = 1 5 = I P ( B ) 。由 | T r ( P [ A , B ] ) | 2 = 4 ( λ 1 − λ 2 ) 2 | I m ( a 12 b 12 ¯ ) | 2 知, 1 4 | T r ( P [ A , B ] ) | 2 = 9 25 ,则 I P ( A ) I P ( B ) − | T r ( P [ A , B ] ) | 2 4 = 1 25 − 9 25 = − 8 25 < 0 ,即有 I P ( A ) I P ( B ) < | T r ( P [ A , B ] ) | 2 4 。
例4 对于给定的密度矩阵P,不一定有 C o r r P G ( X , Y ) = C o r r P ( X , Y ) , X , Y ∈ M n 。
当P为纯态时,设 P = ( 1 0 0 0 ) , X = ( 0 1 1 0 ) , Y = ( 0 − i i 0 ) ,可得 X Y * = ( i 0 0 − i ) , Y * X = ( − i 0 0 i ) , T r ( P X Y * ) = i , T r ( P Y * X ) = − i ,于是 C o r r P G ( X , Y ) = 1 2 T r ( P X Y * + P Y * X ) − T r ( P X P Y * ) = − T r ( P X P Y * ) , C o r r P ( X , Y ) = T r ( P X Y * ) − T r ( P X P Y * ) = i − T r ( P X P Y * ) ,所以 C o r r P G ( X , Y ) ≠ C o r r P ( X , Y ) 。
本文首先对Hermite矩阵组成的实线性空间Hn上的一些数字特征进行了数值计算,并且给出了一个具体的例子,通过例子体现了矩阵数字特征与量子信息中的不确定性关系密切相关。其次对原有的WY相关系数进行改造,使其成为矩阵空间Mn上的内积并且通过矩阵范数、矩阵迹以及矩阵数字特征的性质,得出了改造后的WY相关系数与原WY相关系数和斜信息之间的关系以及矩阵协方差,原WY相关系数和斜信息的一些性质,并且给出了改造后的WY相关系数的一个应用以及斜信息不一定满足Heisenberger不确定性关系的一个例子。所获得的结果充实并完善了矩阵数字特征的性质,为矩阵在量子信息中的研究和应用提供了理论资源。
国家自然科学基金资助项目(批准号:11471200)。
禹鸥洋,任芳国. 关于矩阵数字特征的一些注记Some Notes on the Numerical Characteristics of Matrix[J]. 理论数学, 2024, 14(02): 482-491. https://doi.org/10.12677/PM.2024.142047
https://doi.org/10.1017/9781316848142
https://doi.org/10.1017/CBO9780511810817
https://doi.org/10.1073/pnas.49.6.910
https://doi.org/10.1103/PhysRevLett.91.180403
https://doi.org/10.1109/TIT.2004.831853
https://doi.org/10.1103/PhysRevA.72.042110
https://doi.org/10.1109/TIT.2005.858971
https://doi.org/10.1103/PhysRevA.73.022324
https://doi.org/10.1103/PhysRevA.79.052106
https://doi.org/10.1016/j.jmaa.2009.02.043
https://doi.org/10.1016/j.jmaa.2010.02.037
https://doi.org/10.1016/j.jmaa.2009.09.060
https://doi.org/10.1007/s11128-015-0949-5
https://doi.org/10.1007/s11128-016-1274-3
https://doi.org/10.1016/j.physleta.2021.127222