Processing math: 61%

数学物理学报, 2025, 45(3): 960-971

关于高斯乘积不等式的新结果 (I)

马丽,1,2, 陈蓬颖,1, 韩新方,1,2,*

1海南师范大学数学与统计学院 海口 571158

2海南师范大学数据科学与智慧教育教育部重点实验室 海口 571158

New Results On Gauss Product Inequalities (I)

Ma Li,1,2, Chen Pengying,1, Han Xinfang,1,2,*

1Department of Mathematics and Statistics, Hainan Normal University, Haikou 571158

2Key Laboratory of Data Science and Intelligence Education (Hainan Normal University), Ministry of Education, Haikou 571158

通讯作者: 韩新方,Email: xfanghan@163.com

收稿日期: 2024-05-7   修回日期: 2025-01-13  

基金资助: 海南省自然科学基金(122MS056)
海南省自然科学基金(124MS056)

Received: 2024-05-7   Revised: 2025-01-13  

Fund supported: Hainan Provincial Natural Science Foundation(122MS056)
Hainan Provincial Natural Science Foundation(124MS056)

作者简介 About authors

马丽,Email:malihnsd@163.com;

陈蓬颖,Email:1353525952@qq.com

摘要

(X1,X2,X3) 为中心化的高斯随机变量, 其协方差矩阵的对角线元素均为 1, 该文借助于超几何函数的性质及因式分解得到了

E[|X41X32X33|]E|X41|E|X32|E|X33|,

等号成立当且仅当 X1,X2,X3 相互独立. 从而补充了现有文献中三维高斯乘积不等式的结果.

关键词: 高斯乘积不等式; 正态分布; 超几何函数; 因式分解

Abstract

Let (X1,X2,X3) be a centered Gaussian random vector with D(Xi)=1, i=1,2,3. By means of the properties of hypergeometric function and factorization, we prove that

E[|X41X32X33|]E|X41|E|X32|E|X33|,

and the equal sign holds if and only if X1,X2,X3 are independent. This complements the results of the three dimensional Gauss product inequality in the existing literature.

Keywords: Gauss Product Inequality; normal Distribution; hypergeometric Function; factorization

PDF (522KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

马丽, 陈蓬颖, 韩新方. 关于高斯乘积不等式的新结果 (I)[J]. 数学物理学报, 2025, 45(3): 960-971

Ma Li, Chen Pengying, Han Xinfang. New Results On Gauss Product Inequalities (I)[J]. Acta Mathematica Scientia, 2025, 45(3): 960-971

1 引言

期望为零的高斯随机向量称为中心高斯随机向量. 关于中心高斯随机向量, 文献 [a] 提出了如下猜想: 对任意 Rn 值中心高斯随机向量 (X1,,Xn) 和正实数 α1,αn, 有

E[nj=1|Xj|αj]nj=1E|Xj|αj.
(1.1)

该猜想被称为高斯不等式猜想. 高斯不等式 (1.1) 蕴含着泛函分析中的实极化猜想, 且与统计学中的U猜想有关. 近年来许多学者关注于此猜想的证明.

对复数值中心高斯随机向量, 高斯乘积不等式已经被完全证明了对实数值中心高斯随机向量, 二维 (即 n=2) 高斯乘积不等式也已经被解决, 见文献[2,3,4]. 三维高斯乘积不等式只是部分地被证明了. 本文中 (X1,,Xn) 均为实数值的中心高斯随机向量, N 表示自然数, ρijXiXj 的相关系数. 对三维高斯随机向量, 文献[5,6,7] 得到了α1=2m1,α2=2m2,α3=2m3,m1,m2,m3N时,

(1.1) 式成立. 关于更高维的高斯乘积不等式, 文献 [8] 证明了对任意的 α1,,αn(1,0), 任意的 1kn1, 有

E[nj=1|Xj|αj]E[kj=1|Xj|αj]E[nj=k+1|Xj|αj].

文献 [9] 得到了对任意的 α1(1,0), 有

E[|X1|α1nj=2|Xj|2]ni=2(1ρ21i)E[|X1|α1]nj=2E[|Xj|2].

对任意的 α1(1,0) 和正偶数 α2,α3,α4, 有

E[|X1|α1|X2|α2|X3|α3|X4|α4]4i=2(1ρ21i)αi2E[|X1|α1]4j=2E[|Xj|αi].

综上, n=3,α1,α2,α3 均为正偶数时, (1.1) 式已经被证明了. α1,α2,α3 不全为偶数时,(1.1)式是否成立还尚未解决. 本文借鉴了文献 [10] 第六章降阶的方法, 运用超几何函数的性质得到了 α1=4,α2=3,α3=3 时, (1.1) 式是成立的.不同于文献 [10], 本文得到了含有四次方和三次方的不等式, 从而不能直接用文献[10] 的方法考虑判别式与零的关系. 本文通过换元把含有高次项的不等式转化为含有二次项不等式, 再通过基本不等式得到新的变元为负的.通过证明二次项的系数大于 0, 一次项的系数小于 0, 常数项小于 0, 可得(1.1) 式成立.

2 主要结果及其证明

引理 2.1n3,m1,,mn 为自然数, 若对任意的中心高斯随机向量 (Y1,,Yn)Yn=α1Y1++αn1Yn1,α1,,αn1 为常数, 有

E[n1j=1|Ymjj|Y2kn][n1j=1E|Ymjj|]E[Y2kn],0kmn,

则对任意中心高斯随机向量 (X1,,Xn),

E[n1j=1|Xmjj|X2mnn][n1j=1E|Xmjj|]E[X2mnn].

Z0=E[Xn|X1Xn1],Z1=XnZ0, 那么

X2mnn=(Z0+Z1)2mn=2mni=0Ci2mnZ2mni0Zi1,

Z1 独立于 X1Xn1.i 为奇数时,

E[Z2mni0Zi1|X1Xn1]=Z2mni0E[Zi1]=0,
E[X2mnn|X1Xn1]=mni=0C2i2mnZ2mn2i0E(Z2i1),

注意到 Z0=α1X1++αn1Xn1, 因此,

E[n1j=1|Xmjj|Z2mn2i0][n1j=1E|Xmjj|]E(Z2mn2i0),
E[n1j=1|Xmjj|X2mnn]=E[n1j=1|Xmjj|E(X2mnn|X1Xn1)]=E[n1j=1|Xmjj|(mni=0C2i2mnZ2mn2i0E(Z2i1))]=mni=0C2i2mnE(Z2i1)E[n1j=1|Xmjj|Z2mn2i0]mni=0C2i2mnE(Z2i1)E(Z2mn2i0)[n1j=1E|Xmjj|]=E[mni=0C2i2mnZ2i1Z2mn2i0][n1j=1E|Xmjj|]=E[(Z0+Z1)2mn][n1j=1E|Xmjj|]=[n1j=1E|Xmjj|]E(X2mnn).

推论 2.1 下面两个断言等价

(1) 设 n3,m1,,mn 为自然数, 对任意中心高斯随机向量(Y1,,Yn)Yn=α1Y1++αn1Yn1, 其中 α1,,αn1 为常数, 有

E[n1j=1|Ymjj|Y2kn][n1j=1E|Ymjj|]E(Y2kn),0kmn;

(2) 设 n3,m1,,mn 为自然数, 对任意的中心高斯随机向量(Y1,,Yn), 有

E[n1j=1|Ymjj|Y2mnn][n1j=1E|Ymjj|]E(Y2mnn),

等号成立当且仅当 Y1,,Yn 独立.

定理 2.1(X1,X2,X3) 为中心化的高斯随机变量, 其协方差矩阵的对角线元素均为 1, 则有

E[|X41X32X33|]E|X41|E|X32|E|X33|.
(2.1)

等号成立当且仅当 X1,X2,X3 相互独立.

当 {X1} 与 {X2X3} 独立时, 由二维高斯乘积不等式可得

E[|X41X32X33|]E|X41|E|X32X33|E|X41|E|X32|E|X33|.

当 {X1} 与 {X2,X3} 不独立时, 由推论 2.1 只需证明X1=a2X2+a3X3 时,

E[|X2k1X32X33|]>E|X2k1|E|X32|E|X33|,k=0,1,2.
(2.2)

即可. 由二维高斯乘积不等式可知 k=0 时, (2.2) 式成立. 由文献 [10] 第六章可知 k=1 时, (2.2) 式成立. 下面证明 k=2 时,(2.2) 式成立.a2=0 时, 由 X1=a3X3 及二维高斯乘积不等式可得 (2.2) 式成立. 当 a20 时, 令 Z1=X1a2,Z2=X2,Z3=a3a2X3,α=a3a2,则 Z1=Z2+αZ3, 且 (2.2) 式等价于

E[|Z41Z32Z33|]>E|Z41|E|Z32|E|Z33|,

因此, 为了简化符号, 不妨设 X1=X2+αX3, 只需证明

E[|X41X32X33|]>E|X41|E|X32|E|X33|
(2.3)

即可. 设

F(α,b;c;z)=n=0(α)n(b)n(c)nznn!

为超几何函数, (α)0=1, 当 n1 时, (α)n=α(α+1)(α+n1).x=E[X2X3], 由于 X2X2 同分布, X3X3 同分布, 不妨设 α<0, 0<x<1.由文献 [10,(2.1),(2.2) 式] 即

E[X2m22X2m33]=(2m21)!!(2m31)!![var(X2)]m2[var(X3)]m3F(m2,m3;12;x2),E[X2m2+12X2m3+13]=(2m2+1)!!(2m3+1)!![var(X2)]2m2+12[var(X3)]2m3+12F(m2,m3;23;x2),

可得

E|X41|E|X32|E|X33|=8πE[(X2+αX3)4]=8π[(α4+1)E(X42)+6α2E(X22X23)+(4α+4α3)E(X2X33)]=8π[3(α4+1)+6α2F(1,1;12;x2)+(4α+4α3)3!!xF(0,1;23;x2)]=8π[3(α4+1)+6α2(1+2x2)+3x(4α+4α3)],

由文献 [10,(5.2) 式] 即

E[|X2|y2|X3|y3+2]=(y3+1)2y2+y32Γ(y2+12)Γ(y3+12)πF(y321,y22;12;x2),E[|X2|y2+2|X3|y3]=(y2+1)2y2+y32Γ(y2+12)Γ(y3+12)πF(y221,y32;12;x2),E[|X2|y2X2|X3|y3X3]=x(y2+1)(y3+1)2y2+y32Γ(y2+12)Γ(y3+12)πF(y22,y32;32;x2),

可得

E[|X41X32X33|]=E[(X42+4αX32X3+6α2X22X23+4α3X2X33+α4X43)|X32X33|]=E[|X72X33|+4αX32X3|X32X33|+6α2|X52X53|+4α3X2X33|X32X33|+α4|X32X73|]=128Γ(3)Γ(2)πF(72,32;12;x2)+192xα8Γ(3)Γ(2)πF(52,32;32;x2)+48α28Γ(3)Γ(2)πF(52,52;12;x2)+192xα38Γ(3)Γ(2)πF(32,52;32;x2)+12α48Γ(3)Γ(2)πF(72,32;12;x2).

注意到 F(y22,y3+22;32;x2)=F(y3+22,y22;32;x2), 从而 (2.1) 式等价于

8F(72,32;12;x2)+128xαF(52,32;32;x2)+32α2F(52,52;12;x2)+128xα3F(52,32;32;x2)+8α4F(72,32;12;x2)(α4+1)+2α2(1+2x2)+x(4α+4α3),

上式两边同除以 α2, 可得

(α2+1α2)[8F(72,32;12;x2)1]+(α+1α)[128xF(52,32;32;x2)4x]+32F(52,52;12;x2)2(1+2x2)>0,
(2.4)

α+1α+2=b, 则 α2+1α2=b24b+2,b0, (2.4) 式等价于

b2[8F(72,32;12;x2)1]4b[8F(72,32;12;x2)132xF(52,32;32;x2)+x]+16F(72,32;12;x2)256xF(52,32;32;x2)+32F(52,52;12;x2)4+8x4x2>0.

ddxF(a,b;c;x)=abcF(a+1,b+1;c+1;x),
(2.5)
F(a,b;c;x)=(1x)cabF(ca,cb;c;x),
(2.6)

可得当 0x1 时,

ddx[8F(72,32;12;x2)1]=168xF(52,12;32;x2)=168x(1x2)92F(3,2;32;x2)0.

因此, 二次项系数

8F(72,32;12;x2)18F(72,32;12;0)1=70.

(1z)F(α,b;c;z)=F(α1,b;c;z)c1(cb)zF(α,b;c+1;z),
(2.7)

可得

(1x2)F(y221,y22;12;x2)+(y2+1)x2F(y221,y22;32;x2)=F(y222,y22;12;x2),(1x2)F(52,32;12;x2)+4x2F(52,32;32;x2)=F(72,32;12;x2),(1x2)F(32,52;12;x2)=F(52,52;12;x2)6x2F(32,52;32;x2).

F(b,x)=b2[8F(72,32;12;x2)1]4b[8F(72,32;12;x2)132xF(52,32;32;x2)+x]+16F(72,32;12;x2)256xF(52,32;32;x2)+32F(52,52;12;x2)4(1x)2=b2[8F(72,32;12;x2)1]4b[8(1x2)F(52,32;12;x2)+32x(x1)F(52,32;32;x2)1+x]+48(1x2)F(52,32;12;x2)+256(x2x)F(52,32;32;x2)4(1x)2=8b2[F(72,32;12;x2)18]+32b(1x)[4xF(52,32;32;x2)(1+x)F(52,32;12;x2)+18]+48(1x2)F(32,52;12;x2)256x(1x)F(52,32;32;x2)4(x1)2,

下面证明

8(1x2)F(52,32;12;x2)+32x(x1)F(52,32;32;x2)1+x=(1x)[8(1+x)F(52,32;12;x2)32xF(52,32;32;x2)1]0.

从而F(b,x)的第二项大于等于零.由 (2.6) 式和附录中 (3.2) 式可得

8(1+x)F(52,32;12;x2)32xF(52,32;32;x2)1=8(1+x)(1x2)92F(2,3;12;x2)32x(1x2)112F(4,3;32;x2)1=(1+x)(1x2)32(22.5+30x2)F(1,1;12;x2)(1+x)(1x2)32(14.5+3x2)(1x2)32(7.5+90x2+60x42x)F(1,1;12;x2)+(1x2)32(7.5+41x2+4x42x)1=(1x2)32F(1,1;12;x2)(22.57.52x135x6+30x2)1+(1x2)32(14.5+22.56x+6x3x2x3)=sin1x(22.5x7.52135x26+30x3)+1x2(22.57.52x135x6+30x2)1+1x2(14.5+22.56x+13.5x6+11.5x27x3+3x4+x5)=sin1x(22.5x7.52135x26+30x3)+1x2(820.25x+41.5x27x3+3x4+x5)1,

由于0<x<1, x<sin1x<x1x2, 所以只需证

x(22.5x7.52135x26+30x3)+1x2(820.25x+41.5x27x3+3x4+x5)1>0
(2.8)

x1x2(22.5x7.52135x26+30x3)+1x2(820.25x+41.5x27x3+3x4+x5)1>0.
(2.9)

即可.

先证 (2.8) 式. 由于

820.25x+41.5x27x3+3x4+x5=820.25x+34.5x2+7x2(1x)+3x4+x5=34.5[x220.25x34.5+(20.2534.5×2)2](20.25)234.5×4+8+7x2(1x)+3x4+x5=34.5(x20.2569)2+5.028533+7x2(1x)+3x4+x5>0

1x2>1x, 所以只需证

x(22.5x7.52135x26+30x3)+(1x)(820.25x+41.5x27x3+3x4+x5)1>0,

即可.

x(22.5x7.52135x26+30x3)+(1x)(820.25x+41.5x27x3+3x4+x5)1=732x+84.25x271x3+40x42x5x6=50(x0.32)2+1.88+37x2(x7174)2+737x2+x4(1x)(3+x)>0.

因此, (2.8) 式成立. 下面证 (2.9) 式.

x1x2(22.5x7.52135x26+30x3)+1x2(820.25x+41.5x27x3+3x4+x5)1=11x2(824x+56x29.25x38.5x4+8x53x6x7)1=11x2[3.2+30(x0.4)2+8.25x2+9.25x2(1x)+8.5x2(1x2)+4x5+3x5(1x)+x5(1x2)]111x21>0.

综上所述: 8(1+x)F(52,32;12;x2)32xF(52,32;32;x2)1>0.

下面证明 F(b,x) 的第三项大于等于零. 令

G(x)=48(1x2)F(32,52;12;x2)256x(1x)F(52,32;32;x2)4(x1)2=4(1x)[12(1+x)F(32,52;12;x2)(1x)64xF(52,32;32;x2)].
(2.10)

注意到 G(1)=0, 只需证明 0<x<1 时, G(x)<0. 由 (2.5) 式可得

G(x)8=12xF(32,52;12;x2)+90x(1x2)F(12,32;32;x2)(x1)32(12x)F(52,32;32;x2)160x2(1x)F(32,12;52;x2)=12x(1x2)92F(2,3;12;x2)+90x(1x2)92F(2,3;32;x2)(x1)32(12x)(1x2)112F(4,3;32;x2)160x2(1x)(1x2)92F(4,3;52;x2),

由附录 (3.3) 式和

F(1,1;32;x2)=(1x2)12F(12,12;32;x2)=(1x2)12sin1xx,(1x2)F(1,1;12;x2)=1+x2F(1,1;32;x2)=1+x2(1x2)12sin1xx,(1x2)32F(1,1;12;x2)=(1x2)12+xsin1x,

可得

G(x)8=12x(1x2)92F(2,3;12;x2)+90x(1x2)92F(2,3;32;x2)(x1)32(12x)(1x2)112F(4,3;32;x2)160x2(1x)(1x2)92F(4,3;52;x2)=3x(1x2)32[22.5+30x22F(1,1;12;x2)14.5+3x22]+22.5(1x2)52[1.5+6x22xF(1,1;12;x2)1.5+x22x](x1)(12x)(1x2)32[7.5+90x2+60x42x2F(1,1;12,x2)7.5+41x2+4x42x2]5(1x)(1x2)32[0.75+9x2+18x4x2F(1,1;12;x2)+0.757.5x22x4x2]=41.25x180x2+303.75x3240x4+75x52x2(1x2)32F(1,1;12;x2)(x1)+(1x2)3241.25x+116x2102.25x3+24x4+3.5x52x2=41.25x180x2+303.75x3240x4+75x52xsin1x+41.25x180x2+303.75x3240x4+75x52x2(1x2)12+(1x2)3241.25x+116x2102.25x3+24x4+3.5x52x2(x1)=41.25180x+303.75x2240x3+75x42sin1x+(1x2)1264+242.75x332x2+180.75x324x43.5x52(x1)=(x1)(75x3165x2+138.75x41.25)2sin1x+(1x2)12(x1)(64178.75x+153.25x227.5x33.5x4)2(x1),

因此,

\begin{eqnarray*} \frac{G'(x)}{8(1-x)}&=&\frac{-64+178.75x-153.25x^2+27.5x^3+3.5x^4}{2}(1-x^2)^{\frac{1}{2}}+1\\ &&+\frac{41.25-138.75x+165x^2-75x^3}{2}\sin^{-1}x. \end{eqnarray*}

由于 0<x<1, x\leq{\sin^{-1}x}\leq\frac{x}{\sqrt{1-x^2}}, 因此, 只需要证

\begin{matrix} &&\frac{(75x^3-165x^2+138.75x-41.25)x}{2}\nonumber\\ &&+(1-x^2)^{\frac{1}{2}}\frac{(64-178.75x+153.25x^2-27.5x^3-3.5x^4)}{2}-1>0 \end{matrix}
(2.11)

\begin{matrix} &&\frac{x}{\sqrt{1-x^2}}\frac{(75x^3-165x^2+138.75x-41.25)}{2}\nonumber\\ &&+(1-x^2)^{\frac{1}{2}}\frac{(64-178.75x+153.25x^2-27.5x^3-3.5x^4)}{2}-1>0. \end{matrix}
(2.12)

先证 (2.11) 式. 用 Mathematica 软件对64-178.75x+153.25x^2-27.5x^3-3.5x^4 进行因式分解, 可得

\begin{eqnarray*} && 64-178.75x+153.25x^2-27.5x^3-3.5x^4\\ &=&-3.5(-2.55432+x)(11.9059+x)[(x-0.747225)^2+0.04293179]>0, \end{eqnarray*}

所以 (2.11) 式左侧为

\begin{eqnarray*} &&\frac{75x^4-165x^3+138.75x^2-41.25x}{2}\\ &&+(1-x^2)^{\frac{1}{2}}\frac{(64-178.75x+153.25x^2-27.5x^3-3.5x^4)}{2}-1\\ &\geq&\frac{75x^4-165x^3+138.75x^2-41.25x}{2}\\ &&+(1-x^2)\frac{(64-178.75x+153.25x^2-27.5x^3-3.5x^4)}{2}-1\\ &=&\frac{62-220x+228x^2-13.75x^3-81.75x^4+27.5x^5+3.5x^6}{2}\\ &=&\frac{3.5(1.81006+x)(10.0682+x)(2.12587-2.67622x+x^2)(0.457239-1.34367x+x^2)}{2}\\ &=&\frac{3.5(1.81006+x)(10.0682+x)\big[(x-1.33811)^2+0.3353316279\big]\big[(x-0.67244)^2+0.0050634464\big]}{2}\\ &>&0, \end{eqnarray*}

其中倒数第二个等号用到了 Mathematica 因式分解. 下面证明 (2.12) 式.

\begin{eqnarray*} &&\frac{x}{\sqrt{1-x^2}}\frac{(75x^3-165x^2+138.75x-41.25)}{2}\\ &&+(1-x^2)^{\frac{1}{2}}\frac{(64-178.75x+153.25x^2-27.5x^3-3.5x^4)}{2}-1\\ &=&\frac{64-220x+228x^2-13.75x^3-81.75x^4+27.5x^5+3.5x^6}{2}(1-x^2)^{-\frac{1}{2}}-1\\ &=&\frac{64-220x+228x^2-13.75x^3-81.75x^4+27.5x^5+3.5x^6-2(1-x^2)^{\frac{1}{2}}}{2(1-x^2)^{\frac{1}{2}}}\\ &\geq&\frac{62-220x+228x^2-13.75x^3-81.75x^4+27.5x^5+3.5x^6}{2(1-x^2)^{\frac{1}{2}}}\\ &=&\frac{3.5(1.81006+x)(10.0682+x)(2.12587-2.67622x+x^2)(0.457239-1.34367x+x^2)}{2(1-x^2)^{\frac{1}{2}}}\\ &=&\frac{3.5(1.81006\!+\!x)(10.0682\!+\!x)[(x\!-\!1.33811)^2\!+\!0.3353316279][(x\!-\!0.67244)^2\!+\!0.0050634464]}{2(1-x^2)^{\frac{1}{2}}}\\ &>&0. \end{eqnarray*}

因此, G'(x)<0. 定理证明完毕.

3 附录

本节给出第二节中用到的超几何函数的详细计算过程.

由相邻函数 (2.7) 式和

\big[2\alpha-c+(b-\alpha)z\big]F(\alpha,b;c;z)=\alpha(1-z)F(\alpha+1,b;c;z)-(c-\alpha)F(\alpha-1,b;c;z),
(3.1)

可得

\begin{matrix} F(3,2;\frac{1}{2};x^2)&=&\frac{1.5+2x^2}{1-x^2}F(1,3;\frac{1}{2};x^2)-\frac{0.5}{1-x^2}F(0,3;\frac{1}{2};x^2)\nonumber\\ &=&\frac{1.5+2x^2}{1-x^2}\frac{15}{8(1-x^2)^2}F(1,1;\frac{1}{2};x^2)-\frac{1.5+2x^2}{1-x^2}\frac{7-2x^2}{8(1-x^2)^2}-\frac{0.5}{1-x^2}\nonumber\\ &=&\frac{22.5+30x^2}{8(1-x^2)^3}F(1,1;\frac{1}{2};x^2)-\frac{14.5+3x^2}{8(1-x^2)^3},\nonumber\\ F(4,3;\frac{3}{2},x^2)&=&\frac{(6-\frac{3}{2})F(3,3;\frac{3}{2};x^2)-\frac{3}{2}F(2,3;\frac{3}{2};x^2)}{3(1-x^2)}\nonumber\\ &=&\frac{3F(3,3;\frac{3}{2};x^2)-F(2,3;\frac{3}{2};x^2)}{2(1-x^2)}\nonumber\\ &=&\frac{13.5+108x^2+36x^4}{64x^2(1-x^2)^4}F(1,1;\frac{1}{2};x^2)-\frac{13.5-39x^2}{64x^2(1-x^2)^4}\nonumber\\ &-&\frac{1.5+6x^2}{16x^2(1-x^2)^3}F(1,1;\frac{1}{2};x^2)+\frac{1.5+x^2}{16x^2(1-x^2)^3}\nonumber\\ &=&\frac{7.5+90x^2+60x^4}{64x^2(1-x^2)^4}F(1,1;\frac{1}{2};x^2)-\frac{7.5+41x^2+4x^4}{64x^2(1-x^2)^4}. \end{matrix}
(3.2)

由 (2.7) 式和 (3.1) 式可得

\begin{matrix} F(3,3;\frac{1}{2};x^2)&=&\frac{3.5+x^2}{2(1-x^2)}F(2,3;\frac{1}{2};x^2)-\frac{1.5}{2(1-x^2)}F(1,3;\frac{1}{2};x^2)\nonumber\\ &=&\frac{3.5+x^2}{2(1-x^2)}\Big[\frac{1.5+2x^2}{1-x^2}F(1,3;\frac{1}{2};x^2)-\frac{0.5}{1-x^2}F(0,3;\frac{1}{2};x^2)\Big]\nonumber\\ &&-\frac{1.5}{2(1-x^2)} F(1,3;\frac{1}{2};x^2)\nonumber\\ &=&\frac{7.5+20x^2+4x^4}{4(1-x^2)^2}F(1,3;\frac{1}{2};x^2)-\frac{3.5+x^2}{4(1-x^2)^2}\nonumber\\ &=&\frac{7.5+20x^2+4x^4}{4(1-x^2)^2}\Big[\frac{15}{8(1-x^2)^2}F(1,1;\frac{1}{2};x^2)-\frac{7-2x^2}{8(1-x^2)^2}\Big]-\frac{3.5+x^2}{4(1-x^2)^2}\nonumber\\ &=&\frac{112.5+300x^2+60x^4}{32(1-x^2)^4}F(1,1;\frac{1}{2};x^2)-\frac{80.5+77x^2}{32(1-x^2)^4},\nonumber\\ F(3,1;\frac{1}{2};x^2)&=&\frac{15}{8(1-x^2)^2}F(1,1;\frac{1}{2};x^2)-\frac{7-2x^2}{8(1-x^2)^2},\nonumber\\ F(3,3;\frac{3}{2};x^2)&=&\frac{1}{5x^2}(1-x^2)F(3,3;\frac{1}{2};x^2)-\frac{1}{5x^2}F(2,3;\frac{1}{2};x^2)\nonumber\\ &=&\frac{22.5+60x^2+12x^4}{32x^2(1-x^2)^3}F(1,1;\frac{1}{2};x^2)-\frac{16.1+15.4x^2}{32x^2(1-x^2)^3}\nonumber\\ &&-\frac{18+24x^2}{32x^2(1-x^2)^3}F(1,1;\frac{1}{2};x^2)+\frac{11.6+2.4x^2}{32x^2(1-x^2)^3}\nonumber\\ &=&\frac{4.5+36x^2+12x^4}{32x^2(1-x^2)^3}F(1,1;\frac{1}{2};x^2)-\frac{4.5+13x^2}{32x^2(1-x^2)^3},\nonumber\\ F(2,2;\frac{1}{2};x^2) &=&\frac{1.5+x^2}{1-x^2}F(1,2;\frac{1}{2};x^2)-\frac{0.5}{1-x^2}F(0,2;\frac{1}{2};x^2)\nonumber\\ &=&\frac{1.5+x^2}{1-x^2}\Big[\frac{1.5}{1-x^2}F(1,1;\frac{1}{2};x^2)-\frac{0.5}{1-x^2}F(0,1;\frac{1}{2};x^2)\Big]-\frac{0.5}{1-x^2}\nonumber\\ &=&\frac{9+6x^2}{4(1-x^2)^2}F(1,1;\frac{1}{2};x^2)-\frac{5}{4(1-x^2)^2},\nonumber\\ F(3,2;\frac{3}{2};x^2) &=&\frac{1-x^2}{3x^2}F(3,2;\frac{1}{2};x^2)-\frac{1}{3x^2}F(2,2;\frac{1}{2};x^2)\nonumber\\ &=&\frac{1-x^2}{3x^2}\Big[\frac{22.5+30x^2}{8(1-x^2)^3}F(1,1;\frac{1}{2};x^2)-\frac{14.5+3x^2}{8(1-x^2)^3}\Big]\nonumber\\ &&-\frac{1}{3x^2}\Big[\frac{9+6x^2}{4(1-x^2)^2}F(1,1;\frac{1}{2};x^2)-\frac{5}{4(1-x^2)^2}\Big]\nonumber\\ &=&\frac{1.5+6x^2}{8x^2(1-x^2)^2}F(1,1;\frac{1}{2};x^2)-\frac{1.5+x^2}{8x^2(1-x^2)^2},\nonumber\\ F(4,3;\frac{5}{2};x^2)&=&(\frac{1}{x^2}-1)F(4,3;\frac{3}{2};x^2)-\frac{1}{x^2}F(3,3;\frac{3}{2};x^2)\nonumber\\ &=&(\frac{1}{x^2}-1)\Big[\frac{7.5+90x^2+60x^4}{64x^2(1-x^2)^4}F(1,1;\frac{1}{2},x^2)-\frac{7.5+41x^2+4x^4}{64x^2(1-x^2)^4}\Big]\nonumber\\ &&-\frac{1}{x^2}\Big[\frac{4.5+36x^2+12x^4}{32x^2(1-x^2)^3}F(1,1;\frac{1}{2};x^2)-\frac{4.5+13x^2}{32x^2(1-x^2)^3}\Big]\nonumber\\ &=&\frac{-0.75+9x^2+18x^4}{32x^4(1-x^2)^3}F(1,1;\frac{1}{2};x^2)+\frac{0.75-7.5x^2-2x^4}{32x^4(1-x^2)^3}. \end{matrix}
(3.3)

参考文献

Li W V, Wei A.

A Gaussian inequality for expected absolute products

J Theor Probab, 2012, 25: 92-99

Frenkel P E. Pfaffians,

Hafnians and product of real linear functionals

Math Res Lett, 2008, 15: 351-358

[本文引用: 1]

Karlin S, Rinott Y.

Total positivity properties of absolute value multinormal variables with applications to confidence interval estimates and related probabilistic inequalities

Ann Stat, 1981, 9: 1035-1049

[本文引用: 1]

Lan G L, Hu Z C, Sun W.

The three-dimensional gaussian product inequality

J Math Anal, 2020, 485(2): 123858

[本文引用: 1]

Russell O, Sun W.

Moment ratio inequality of bivariate gaussian distribution and three-dimensional gaussian product inequality

J Math Anal Appl, 2023, 527(1): 127410

[本文引用: 1]

Herry R, Malicet D, Poly G. A short proof of the strong three dimensional gaussian product inequality. arXiv:2211.07314v1

[本文引用: 1]

Wei A.

Representations of the absolute value function and applications in gaussian estimates

J Theor Probab, 2014, 27: 1059-1070

[本文引用: 1]

Russell O, Sun W.

An opposite gaussian product inequality

Statist Probab Lett, 2022, 191: 109656

[本文引用: 1]

Zhou Q Q, Zhao H, Hu Z C, Song R M.

Song new results on gaussian product inequalities

Journal of Mathematical Analysis and Applications, 2024, 531(2): 127907

[本文引用: 1]

Russell O. On the Gaussian Product Inequality. PhD thesis. Montreal: Concordia University, 2023. https://spectrum.library.concordia.ca/id/eprint/992526/

URL     [本文引用: 6]

Rainville E D. Special Functions. New York: The Macmillan Company, 1960

/