Processing math: 57%

数学物理学报, 2021, 41(3): 770-782 doi:

论文

Rn中线性子空间束与凸体相交的几何概率

赵江甫,

Geometric Probability of Subspaces Intersecting with a Convex Body in Rn

Zhao Jiangfu,

收稿日期: 2019-09-18  

基金资助: 福建省中青年教师教育科研项目基金(科技类).  JT180585
福建江夏学院科研培育人才项目基金.  JXZ2019016

Received: 2019-09-18  

Fund supported: the Educational Research Project Fund of Young and Middle-Aged Teachers of Fujian Province.  JT180585
the Project Fund for Scientific Research and Cultivation of Talents of Fujian Jiangxia University.  JXZ2019016

作者简介 About authors

赵江甫,E-mail:2833811196@qq.com , E-mail:2833811196@qq.com

Abstract

The probability that three independent random subspaces in Rn intersecting a convex body K have their common point intersecting K is found by using of the mean curvature integral of convex sets. Then we focus on the particular case of hyperplanes. On the base of this, we state the geometric probability of hyperplanes that intersect a ball, a cube or a right parallelepiped having an intersection inside the same ball, the cube or the right parallelepiped respectively. Finally, the monotonicity, convergence, and size relationship of the geometric probabilistic sequence are discussed.

Keywords: Mean curvature integral ; Geometric probability ; Buffon needle throwing ; Elementary symmetric function ; Hyperplanes ; Subspaces

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

本文引用格式

赵江甫. Rn中线性子空间束与凸体相交的几何概率. 数学物理学报[J], 2021, 41(3): 770-782 doi:

Zhao Jiangfu. Geometric Probability of Subspaces Intersecting with a Convex Body in Rn. Acta Mathematica Scientia[J], 2021, 41(3): 770-782 doi:

1 引言

积分几何起源于几何概率, 其发展也始终与几何概率紧密相连. 最早利用积分学的方法解决的几何概率问题是经典的蒲丰投针问题: 设在平面上有一组平行线, 其行距都等于d, 把一根长度为l(l<d)的针随机地投上去, 则这根针和一条直线相交的概率为2lπd. 蒲丰投针问题不仅是最早的一个几何概率问题, 也是一个最具有代表性、影响最大的几何概率问题. 其优美性和吸引力历时两个多世纪经久不衰. 特别是积分几何的出现, 使人们得以从全新的角度对这类问题予以洞察. 近年来, 经典的蒲丰投针问题得到了较为广泛的推广与研究, 参见文献[1-9]. 但这些推广仅限于2维和3维空间.文献[10-11] 将几何概率问题的研究推广到了n维空间, 其中文献[11] 主要针对线性空间偶的情形.对于有三个以上的线性子空间组成的线性空间束, 仅仅给出了三维空间中的一个特例: 在三维空间中, 三个和凸体K相交的平面的公共点在K内的概率是

p=π4VM3,
(1.1)

其中, VK的体积, MK的边界曲面K的平均曲率积分. 特别地, 当K为半径为R的球体B时, V=43πR3,M=Bk1+k22ds=4πR, 概率为π248.

本文的主要的工作是将(1.1)式推广到Rn中, 即研究Rn中的线性子空间束与凸体相交的几何概率.

2 预备知识

定义2.1   设KRn中的点集, 若对于任意的x,yK都有λx+(1+λ)yK(0λ1), 则称KRn中的凸集.

定义2.2   具有非空内点的紧凸集称为凸体. 凸体K的边界K称为凸曲面.设Cn表示Rn中所有凸集的集合, Kn表示Rn中所有凸体的集合. 对于KKn, V(K)表示它的体积.

定义2.3   设LriRnri维线性子空间, 则称(Lri,Lrj)Rn中的线性子空间偶, 称(Lri,Lr2,,Lrk)为线性子空间束, 称(LriLr2Lrk)为线性子空间束的交, 其中k3.

定义2.4   如果LriK,(i=1,2,,k), 则称线性子空间束(Lri,Lr2,,Lrk)K相交.特别地, 若r1=r2=rk=n1, 则称(Lri,Lr2,,Lrk)为超平面束.

定义2.5  [12-13]nN,xRn+={(x1,x2,,xn)|xi[0,+),1in}, 则称

ek(x)=ek(x1,x2,,xn)=1i1<i2<<iknxi1xi2xik
(2.1)

k阶初等对称函数, 并补充规定e0(x)=e0(x1,x2,,xn)=1.

1阶初等对称函数e1(x1,x2,,xn)=ni=1xi; 第2阶初等对称函数e2(x1,x2,,xn)=n1i<jnxixj.

定义2.6  [9, 11, 14-16]ΣRn中的C2类超曲面, k1,k2,,knΣn1个主曲率函数, 则Σ的第r个平均曲率积分可用这n1个主曲率的初等对称函数定义为

Mr(Σ)=1(n1  r)Σer(k1,k2,,kn)dσ,r=1,2,,n1,
(2.2)

其中, dσ表示Σ的面积元.另外, 补充规定, M0(Σ)=F (即Σ的面积).

ΣRn中半径为的R的球面, 其主曲率为k1=k2=1R, 则Σ3个平均曲率积分分别为

M0(Σ)=F=4πR2;

M1(Σ)=1(21)Σ2Rdσ=1R4πR2=4πR;

M2(Σ)=1(22)Σ1R2dσ=1R24πR2=4π.

注2.1   R3中曲面的第1个平均曲率积分即通常微分几何中所指的平均曲率积分M, 即M1(Σ)=M, 因此文中出现的M1(Σ)将以M替代.

KCn, 若KC2, 则平均曲率积分与均质积分有如下关系[9, 11, 14-18], 即Cauchy公式

Mr(K)=nWr+1(K),r=0,1,,n1.
(2.3)

注2.2   均质积分Wr(K)对于任意KCn都有定义, 而平均曲率积分Mr(K)则要求KC2, 因此可以借助Cauchy公式, 对任意的KCn给出平均曲率积分的定义.当KC2时, 此定义与前述定义一致.

注2.3   为简单起见, 把凸体K的边界曲面K的平均曲率积分简称为K的平均曲率积分.

下面给出几种特殊凸体的平均曲率积分.

当凸体K是棱长为ai (i=1,2,,n)的长方体C时, 其平均曲率积分为

Mr(C)=nOr(r+1)( nr+1)enr1(a1,a2,,an).
(2.4)

当凸体K是棱长为a的正方体Ce时, 其平均曲率积分为

Mr(Ce)=nOranr1(r+1).
(2.5)

当凸体K是半径为R的球体B时, 其平均曲率积分为

Mr(B)=On1Rnr1.
(2.6)

Lr[O]Rn中过定点Or维平面, 则Lr[O]的总测度, 即Grassmann流形Gr,nr的体积为[9, 11]

m(Gr,nr)=Gr,nrdLr[O]=On1On2OnrOr1Or2O0,
(2.7)

式中Oi=2πi+12Γ(i+12)i维单位球面的表面积, 且有(i1)Oi=2πOi2.

Lnr[O]Rn中过定点O且垂直于Lr[O]nr维平面, 则由密度公式dLr[O]=dLnr[O]

m(Gnr,r)=Gnr,rdLnr[O]=Gr,nrdLr[O]=On1On2OnrOr1Or2O0=m(Gr,nr).
(2.8)

3 主要结论

本章主要介绍本文所得出的结论, 共6个定理, 3个推论.其中, 定理3.1和定理3.2给出线性子空间束与凸体相交时的相关测度; 定理3.3是本文的核心结论, 它给出线性子空间束的交与固定凸体相交的几何概率; 取定理3.3中的线性子空间束为特殊的子空间束-超平面束, 便得到定理3.4, 即超平面束与凸体相交的几何概率; 推论3.1–3.3分别给出当定理3.4中的凸体分别为正方体、长方体、球体时的几何概率的具体结果; 定理3.5讨论了超平面束与球体相交时所得的几何概率序列的单调性; 定理3.6描述了推论3.1–3.3所得的三个结果的大小关系.

定理 3.1  设KKn, (Lp,Lq,Ls)Rn中与K相交的线性子空间束, 且维数满足p+qn, 则此线性子空间束的交与K相交的测度m1

(1) 当p+qn+1时,

1) 当p+q+s2n+1时,

m1=4π2(p1)!(q1)!(s1)!Op1Oq1Os1O3npqs+1((n1)!)2(p+q+s2n1)!Onp+1Onq+1Ons+1Op+q+s2n1TnpTnqTnsMp+q+s2n1(K);
(3.1)

2) 当p+q+s=2n时,

m1=8π3(p1)!(q1)!(s1)!Op1Oq1Os1On1((n1)!)2Onp+1Onq+1Ons+1TnpTnqTnsV(K);
(3.2)

3) 当p+q+s<2n时,

m1=0;
(3.3)

(2) 当p+q=n时,

m1=2πp!q!O2n1(n1)!Ons+1TnpTnqTnsV(K),
(3.4)

其中Tni=On2On3Oni1(ni)Oi1Oi2O0.

定理 3.2   设KKn, (Lp,Lq,Ls)Rn中的线性子空间束, 且维数满足p+qn, 则此线性子空间束与K相交的测度m2

m2=TnpTnqTnsMp1(K)Mq1(K)Ms1(K).
(3.5)

由定理3.1、3.2可得出定理3.3.

定理 3.3   设KKn, (Lp,Lq,Ls)Rn中与K相交的线性子空间束, 且维数满足p+qn, 则此线性子空间束的交与K相交的概率pn(K)

(1) 当p+qn+1时,

1) 当p+q+s2n+1时,

pn(K)=4π2(p1)!(q1)!(s1)!Op1Oq1Os1O3npqs+1((n1)!)2(p+q+s2n1)!Onp+1Onq+1Ons+1Op+q+s2n1Mp+q+s2n1(K)Mp1(K)Mq1(K)Ms1(K);
(3.6)

2) 当p+q+s=2n时,

pn(K)=8π3(p1)!(q1)!(s1)!Op1Oq1Os1On1V(K)((n1)!)2Onp+1Onq+1Ons+1Mp1(K)Mq1(K)Ms1(K);
(3.7)

3) 当p+q+s<2n时,

pn(K)=0;
(3.8)

(2) 当p+q=n时,

pn(K)=2πp!q!O2n1V(K)(n1)!Ons+1Mp1(K)Mq1(K)Ms1(K).
(3.9)

特别地, 当p=q=s=2,n=3时, 根据(3.7)式可得, 在三维空间中, 三个和K相交的平面的公共点落在K内的概率为

p3(K)=8π3(21)!(21)!(21)!O21O21O21O31V(K)((31)!)2O32+1O32+1O32+1M21(K)M21(K)M21(K)=8π3O31O2V(K)4O32M31(K)=π4V(K)M3,

这与(1.1)式吻合.

在定理3.3中, 取p=q=s=n1, 可得出定理3.4.

定理 3.4   设KKn, (Ln1,Gn1,Hn1)Rn中与K相交的超平面束, 则此超平面束的交与K相交的概率pn(K)

pn(K)={π(n2)(n3)(On2)3Mn4(K)6(n1)2On4(Mn2(K))3,n4,π4V(K)M3,n=3.
(3.10)

K为一些特殊凸体时, 将对应的平均曲率积分代入(3.10)式, 可得到推论3.1–3.3.

推论 3.1   设CeKn中棱长为a的正方体, (Ln1,Gn1,Hn1)Rn中与Ce相交的超平面束, 则此超平面束的交与Ce相交的概率pn(Ce)

pn(Ce)={π(n1)(n2)6n2,n4,π27,n=3.
(3.11)

推论 3.2   设CKn中棱长为ai,(i=1,2,,n)的长方体, (Ln1,Gn1,Hn1)Rn中与C相交的超平面束, 则此超平面束的交与C相交的概率pn(C)

pn(C)={πe3(a1,a2,,an)(e1(a1,a2,,an))3,n4,πa1a2a327(a1+a2+a3)3,n=3.
(3.12)

推论 3.3   设BKn中半径为R的球体, (Ln1,Gn1,Hn1)Rn中与B相交的超平面束, 则此超平面束的交与B相交的概率pn(B)

pn(B)={π(n2)Γ2(n2)3(n1)2Γ2(n12),n4,π248,n=3.
(3.13)

推论3.1–3.3所得的三个几何概率pn(Ce)pn(C)pn(B)具有以下性质:

定理 3.5   pn(Ce)pn(B)均关于维数n严格单调递增且收敛于常数π6.

定理 3.6   当维数n固定时, pn(Ce)pn(C)pn(B)具有以下大小关系:

pn(C)pn(Ce)<pn(B).
(3.14)

4 定理的证明

4.1 定理3.1的证明

定理3.1的证明需要用到引理4.1–4.5.

引理 4.1[9, 11]  设KKn,Lp,LqRn中与K相交的两个线性子空间, 且维数满足p+qn, 则LpLqK相交的测度为

(1) 当p+qn+1时,

m=2π(p1)!(q1)!Op1Oq1O2npq+1(n1)!(p+qn1)!Onp+1Onq+1Op+qn1TnpTnqMp+qn1(K);
(4.1)

(2) 当p+q=n时,

m=p!q!On1(n1)!TnpTnqV(K).
(4.2)

引理 4.2[9, 11]   设Krr维平面LrRn中的凸体, 并假定其边界Kr二次可微, 当把Kr视为Lr中的凸曲面时, 以M(r)q(Kr)表示其平均曲率积分.当把Kr视为Rn(r<n)中的平坦凸体时, Kr的平均曲率积分则以M(n)q(Kr)记之, 其中, q=0,1,,r1. 二者关系如下:

(1) 当qnr时,

M(n)q(Kr)=( r1qn+r)Oq(n1 q)Oqn+rM(r)qn+r(Kr);
(4.3)

(2) 当q=nr1时,

M(n)nr1(Kr)=Onr1( n1nr1)Vr(Kr),
(4.4)

其中Vr(Kr)Krr维体积;

(3) 当q<nr1时,

M(n)q(Kr)=0.
(4.5)

引理 4.3[9, 11]   平均曲率积分M(r)q((LrK))在集{Lr:LrK}上的积分为

LrKM(r)q((LrK))dLr=On2On3OnrOnqOr2Or3O0OrqMq(K).
(4.6)

引理 4.4   设LrRnr维平面, σr(LrK)表示LrKr维体积, 则有积分公式

LrKσr(LrK)dLr=On1On2OnrOr1Or2O0V(K).
(4.7)

   设p=LrLnr[O],dσnrLnr[O]p点的体积元, 则有密度公式

dLr=dσnrdLnr[O].

对上式两边积分得

LrKσr(LrK)dLr=LrKσr(LrK)dσnrdLnr[O]=Gnr,r(LrKσr(LrK)dσnr)dLnr[O]=V(K)Gnr,rdLnr[O].

将(2.8)式代入上式, 引理4.4得证.

下面证明定理3.1.

   LpLqLsK相交的测度为

m1=(LpLqLsK)dLpdLqdLs=LsK((LpLq(LsK)dLpdLq)dLs.
(4.8)

下面分两种情况计算内层积分.

(1) 当p+qn+1时, 由引理4.1

(LpLq(LsK)dLpdLq=2π(p1)!(q1)!Op1Oq1O2npq+1(n1)!(p+qn1)!Onp+1Onq+1Op+qn1TnpTnqMp+qn1((LsK)).
(4.9)

下面分三种情况计算积分LsKM(s)p+q+s2n1(Ks)dLs.

1) 当p+q+s2n+1时, 即p+qn1ns, 由(4.3)式得

Mp+qn1((LsK))=M(n)p+qn1(Ks)=(s1p+q+s2n1)Op+qn1(n1p+qn1)Op+q+s2n1M(s)p+q+s2n1(Ks).
(4.10)

由引理4.3

LsKM(s)p+q+s2n1(Ks)dLs=On2On3OnsO3npqs+1Os2Os3O0O2n+1pqMp+q+s2n1(K)=2πOs1O3npqs+1Ons+1O2n+1pqTnsMp+q+s2n1(K).
(4.11)

将(4.9)–(4.11)式代入(4.8)式, (3.1)式得证.

2) 当p+q+s=2n时, 即p+qn1=ns1, 由(4.4)式得

Mp+qn1((LsK))=M(n)p+qn1(Ks)=Ons1( n1ns1)Vs(Ks)=Ons1( n1ns1)σs(LsK).
(4.12)

由引理4.4

LsKσs(LsK)dLs=On1On2OnsOs1Os2O0V(K)=2πOn1Ons+1TnsV(K).
(4.13)

将(4.9)、(4.12)和(4.13)式代入(4.8)式, (3.3)式得证.

3) 当p+q+s<2n时, 即p+qn1<ns1, 由(4.5)式得

Mp+qn1((LsK))=0;
(4.14)

将(4.9)和(4.14)式代入(4.8)式, (3.2)式得证.

(2) 当p+q=n时, 由引理4.1

LpLq(LsK)dLpdLq=p!q!On1(n1)!TnpTnqV(LsK),
(4.15)

其中V(LsK)表示LsKs维体积, 这与σs(LsK)表示的意义相同, 因此有

V(LsK)=σs(LsK).
(4.16)

将(4.13)、(4.15)和(4.16)式代入(4.8)式, (3.4)式得证.

定理3.1得证.

4.2 定理3.2的证明

定理3.2的证明需要用到引理4.5.

引理 4.5[9, 11]   设KKn, LrRnr维平面, 则LrK相交的测度为

mr{Lr:LrK}=LrKdLr=On2On3Onr1(nr)Or1Or2O0Mr1(K).
(4.17)

下面证明定理3.2.

   由引理4.5得

m2=m{(Lp,Lq,Ls):LpK,LqK,LsK}=LpK,LqK,LsKdLpdLqdLs=On2On3Onp1(np)Op1Op2O0Mp1(K)On2On3Onq1(nq)Oq1Oq2O0Mq1(K)On2On3Ons1(ns)Os1Os2O0Ms1(K)=TnpTnqTnsMp1(K)Mq1(K)Ms1(K).

证毕.

4.3 定理3.5的证明

   pn(Ce)显然收敛于π6, 因此定理3.5只需证明以下三个结论:

(1) pn(Ce)关于维数n严格单调递增;

(2) pn(B)关于维数n严格单调递增;

(3) limnpn(B)=π6.

(1) 证明 p_n(C_e) 严格单调递增.

n = 3 时, p_{4}(C_e) = \frac{\pi}{16}>\frac{\pi}{27} = p_3(C_e) ;

n\geq4 时, 因为

\begin{eqnarray*} p_{n+1}(C_e)-p_n(C_e)& = & \frac{\pi n(n-1)}{6(n+1)^2}-\frac{\pi (n-1)(n-2)}{6n^2}\\ & = &\frac{\pi(n-1)}{6}\left(\frac{n^3-(n-2)(n+1)^2}{n^2(n+1)^2}\right)\\ & = &\frac{\pi (n-1)(3n+2)}{6n^2(n+1)^2})>0, \end{eqnarray*}

所以 p_{n+1}(C_e)>p_n(C_e) .

综上, 对任意 n\geq3 , 都有 p_n(C_e) 关于维数 n 严格单调递增.

(2) 证明 p_n(B) 严格单调递增.

n = 3 时, p_3(B) = \frac{\pi^2}{48}<\frac{8}{27} = p_4(B);

n\geq4 时, 因为

\frac{n^4}{(2n+1)^2}-\frac{n^5}{4(n-1)(n+1)^2} = \frac{-n^4(5n+4)}{4(n-1)(n+1)^2(2n+1)^2}<0,

所以

\frac{n^2}{2n+1}<\sqrt{\frac{n^5}{4(n-1)(n+1)^2}}.

由不等式[19-20]

\frac{\Gamma^2(\frac{n+1}{2})}{\Gamma^2(\frac{n}{2})}\leq\frac{n^2}{2n+1}

\frac{\Gamma^2(\frac{n+1}{2})}{\Gamma^2(\frac{n}{2})}<\sqrt{\frac{n^5}{4(n-1)(n+1)^2}},

整理得

\frac{(n-1)\Gamma^2(\frac{n+1}{2})}{n^2\Gamma^2(\frac{n}{2})}<\sqrt{\frac{n(n-1)}{4(n+1)^2}},

两边取对数得

f_n = \ln{\frac{(n-1)\Gamma^2(\frac{n+1}{2})}{n^2\Gamma^2(\frac{n}{2})}} = 2\ln{\Gamma(\frac{n+1}{2})}-2\ln{\Gamma(\frac{n}{2})}+\ln{\frac{n-1}{n^2}}<\ln{\sqrt{\frac{n(n-1)}{4(n+1)^2}}},

所以

\begin{eqnarray*} f_{n+1}-f_n& = & 2\ln{\Gamma(\frac{n+2}{2})}-2\ln{\Gamma(\frac{n+1}{2})}+\ln{\frac{n}{(n+1)^2}}-f_n\\ & = &2\ln{\frac{n}{2}\Gamma(\frac{n}{2})}-2\ln{\Gamma(\frac{n+1}{2})}+\ln{\frac{n}{(n+1)^2}}-f_n\\ & = &\ln{\frac{n(n-1)}{4(n+1)^2}}-2f_n\\ & = &2\left(\ln{\sqrt{\frac{n(n-1)}{4(n+1)^2}}-f_n}\right)>0. \end{eqnarray*}

这样就证明了序列 f_n 是关于维数 n 严格单调递增的. 又因为 p_n(B) = \frac{\pi}{3}e^{f_{n-1}} , 所以当 n\geq4 时, p_n(B) 关于维数 n 严格单调递增.

综上, 对于任意的 n\geq3 , 都有 p_n(B) 关于维数 n 严格单调递增.

(3) 证明 \lim \limits_{n\rightarrow \infty}p_n(B) = \frac{\pi}{6} .

这里需要用到著名的Wallis's公式[19]

\begin{eqnarray*} \lim\limits_{n\rightarrow \infty}\frac{(2m)!!^2}{m(2m-1)!!^2} = \pi. \end{eqnarray*}

(1) 当 n = 2m 时, 由Wallis's公式可得偶子列的极限为

\begin{eqnarray*} \lim\limits_{m\rightarrow \infty}p_{2m}(B)& = & \lim\limits_{m\rightarrow \infty}\frac{\pi (2m-2)\Gamma^2(m)}{3(2m-1)^2\Gamma^2(\frac{2m-1}{2})}\\ & = & \lim\limits_{m\rightarrow \infty}\frac{\pi (2m-2)(m-1)!^2}{3(2m-1)^2\pi \left(\frac{(2m-3)!!}{2^{m-1}}\right)^2}\\ & = & \lim\limits_{m\rightarrow \infty}\frac{(2m-2)(2m-2)!!^2}{3(2m-1)^2(2m-3)!!^2}\\ & = & \lim\limits_{m\rightarrow \infty}\frac{(2m)!!^2}{m(2m-1)!!^2}\cdot\frac{m-1}{6m} = \frac{\pi}{6}; \end{eqnarray*}

(2) 当 n = 2m+1 时, 由Wallis's公式可得奇子列的极限为

\begin{eqnarray*} \lim\limits_{m\rightarrow \infty}p_{2m+1}(B)& = & \lim\limits_{m\rightarrow \infty}\frac{\pi (2m-1)\Gamma^2(\frac{2m+1}{2})}{3(2m)^2\Gamma^2(m)}\\ & = & \lim\limits_{m\rightarrow \infty}\frac{\pi^2 (2m-1)\left(\frac{(2m-1)!!}{2^{m}}\right)^2}{12m^2(m-1)!^2}\\ & = & \lim\limits_{m\rightarrow \infty}\frac{\pi^2 (2m-1)(2m-1)!!^2}{48m^2(2m-2)!!^2}\\ & = & \lim\limits_{m\rightarrow \infty}\frac{m(2m-1)!!^2}{(2m)!!^2}\cdot\frac{\pi^(2m-1)}{12m} = \frac{\pi}{6}. \end{eqnarray*}

综上可知, \lim \limits_{n\rightarrow \infty}p_n(B) = \frac{\pi}{6} . 定理3.5得证.

4.4 定理3.6的证明

定理3.6描述的是三种特殊凸体: 长方体、正方体、球体所对应的几何概率的大小关系, 本节将分两步进行证明.首先根据均值不等式, 讨论长方体与正方体的情形, 然后根据文献[19] 中的两个不等式讨论正方体与球体的情形.

   (1) 首先证明 p_n(C)\leq p_n(C_e) .

(i) 当 n = 3 时, 由 a_1, a_2, a_3>0 可得

(a_1+a_2+a_3)^3\geq27a_1a_2a_3,

所以

p_3(C) = \frac{\pi a_1a_2a_3}{27(a_1+a_2+a_3)^3}\leq\frac{\pi a_1a_2a_3}{27\times27a_1a_2a_3} = \frac{\pi}{729}<\frac{\pi}{27} = p_3(C_e),

{\rm{ (当且仅当}} a_1 = a_2 = a_3 {\rm{时等号成立);}}

(ii) 当 n\geq4 时, 令 A = \sum\limits_{i = 1}^na^3_i, B = \sum\limits_{i\neq j}^na_ia^2_j, C = \sum\limits_{1\leq i<j\leq n}^na_ia_ja_k , 由 a_i, a_j, a_k>0

\begin{eqnarray*} \left\{\begin{array}{ll} a_ia_j^2+a_ia^2_k\geq2a_ia_ja_k , \\ a_ja_i^2+a_ja^2_k\geq2a_ia_ja_k , \\ a_ka_j^2+a_ka^2_i\geq2a_ia_ja_k , \end{array}\right. \end{eqnarray*}

将以上三个不等式相加并两边作和, 得

(n-2)\sum\limits_{i\neq j}^na_ia^2_j\geq6\sum\limits_{1\leq i<j\leq n}^na_ia_ja_k,

(n-2)B\geq6C, 所以

\begin{equation} B\geq\frac{6C}{n-2} \; \; {\rm{(当且仅当}}| a_1 = a_2 = \cdots = a_n {\rm{时等号成立)}}. \end{equation}
(4.18)

又因为

\begin{eqnarray*} (a_i+a_j+a_k)^3& = & a^3_i+a^3_j+a^3_k+3a_ia^2_j+3a_ia^2_k+3a_ja^2_k+3a_ja^2_i+3a_ka^2_i+3a_ka^2_j + 6a_ia_ja_k\\ &\geq&27a_ia_ja_k, \end{eqnarray*}

所以

a^3_i+a^3_j+a^3_k+3a_ia^2_j+3a_ia^2_k+3a_ja^2_k+3a_ja^2_i+3a_ka^2_i+3a_ka^2_j\geq21a_ia_ja_k.

将以上不等式两边作和可得

\left(^{n-1}_{\ \, 2}\right)\sum\limits_{i = 1}^na^3_i+3\left(^{n-2}_{\ \ 1}\right)\sum\limits_{i\neq j}^na_ia^2_j\geq21\sum\limits_{1\leq i<j\leq n}^na_ia_ja_k,

\begin{equation} \frac{(n-1)(n-2)}{2}A+3(n-2)B\geq21C. \end{equation}
(4.19)

由(4.18)–(4.19)式得

\begin{eqnarray*} &&\frac{(n-1)(n-2)}{2}(A+3B+6C) \\ & = &\frac{(n-1)(n-2)}{2}A+3(n-2)B+\frac{3(n-2)(n-3)}{2}B+3(n-1)(n-2)C\\ &\geq & 21C+\frac{3(n-2)(n-3)}{2}\cdot\frac{6C}{n-2}+3(n-1)(n-2)C\\ & = & 3n^2C, \end{eqnarray*}

所以

A+3B+6C\geq\frac{6n^2C}{(n-1)(n-2)},

所以

\begin{eqnarray*} p_n(C)& = & \frac{\pi e_3(a_1, a_2, \cdots, a_n)}{(e_1(a_1, a_2, \cdots, a_n))^3} = \frac{\pi\sum\limits_{1\leq i<j\leq n}^na_ia_ja_k }{\Big(\sum\limits_{i = 1}^na_i\Big)^3} = \frac{\pi C}{A+3B+6C}\\ & \leq& \frac{\pi C}{\frac{6n^2C}{(n-1)(n-2)}} = \frac{\pi (n-1)(n-2)}{6n^2} = p_n(C_e), \end{eqnarray*}

上式子中当且仅当 a_1 = a_2 = \cdots = a_n 时等号成立.

综上, 对于任意的 n\geq3 都有 p_n(C)\leq p_n(C_e) .

(2) 证明p_n(C_e)< p_n(B). 这里需要用到以下两个不等式[19]

\begin{equation} \frac{(2m)!!^2}{(2m-1)!!^2}>\frac{\pi (4m+1)}{4}, m = 1, 2, \cdots, \end{equation}
(4.20)

\begin{equation} \frac{(2m-1)!!^2}{(2m)!!^2}>\frac{4m+3}{\pi(2m+1)^2}, m = 1, 2, \cdots. \end{equation}
(4.21)

1) 当 n = 3 时, p_3(C_e) = \frac{\pi}{27}<\frac{\pi^2}{48} = p_3(B) ;

2) 当 n\geq4 时, 令

\begin{eqnarray*} h_n = \frac{p_n(B)}{p_n(C)} = \frac{\pi (n-2)\Gamma^2(\frac{n}{2})}{3(n-1)^2\Gamma^2(\frac{n-1}{2})}\cdot\frac{6n^2}{\pi(n-1)(n-2)} = \frac{2n^2\Gamma^2\frac{n}{2}}{(n-1)^3\Gamma^2(\frac{n-1}{2})}, \end{eqnarray*}

(i) 当 n = 2m 时, 由(4.20)式得

\begin{eqnarray*} h_{2m} = \frac{2(2m)!!^2}{(2m-1)(2m-1)!!^2\pi}>\frac{2}{\pi(2m-1)}\cdot\frac{\pi(4m+1)}{4} = \frac{4m+1}{2(2m-1)}>1; \end{eqnarray*}

(ii) 当 n = 2m+1 时, 由(4.21)式得

\begin{eqnarray*} h_{2m+1} = \frac{2(2m+1)^2(2m+1)!!^2\pi}{2m(2m)!!^2}>\frac{\pi(2m+1)^2}{4(m-1)}\cdot\frac{4m+3}{\pi(2m+1)^2} = \frac{4m+3}{4(m-1)}>1. \end{eqnarray*}

所以, 对于任意的 n\geq4 都有 h_n>1 , 即 p_n(C)<p_n(C_e) .

综上可知, 对于任意的 n\geq3 都有 p_n(C)< p_n(C_e) ,

因此, 对于任意的 n\geq3 都有 p_n(C)\leq p_n(C)< p_n(C_e) .定理3.6得证.

5 结语

本文虽然仅讨论由三个线性子空间组成的线性子空间束与凸体相交的几何概率问题, 但却给出了一种求解此类问题的方法.比如, 可以进一步讨论由 h 个线性子空间组成的线性子空间束与凸体相交的几何概率问题.另外, 本文仅讨论了超平面束与特殊凸体相交时所得的几何概率的单调性、收敛性以及大小关系, 任意维的线性子空间束与一些特殊凸体相交的几何概率的性质并未讨论.这些问题有待于进一步研究.

参考文献

Xie F F , Li D Y .

On generalized Buffon needle problem for lattices

Acta Mathematica Scientia, 2011, 31B (1): 303- 308

URL     [本文引用: 1]

谢凤繁, 胡适耕.

E3中有界网格的Buffon问题

数学物理学报, 2011, 31A (2): 410- 414

URL    

Xie F F , Hu S G .

On Buffon problem for bounded lattices in E3

Acta Math Sci, 2011, 31A (2): 410- 414

URL    

李德宜, 蒋慧峰, 熊革.

凸体形状匹配的参考点

数学杂志, 2005, 25 (3): 336- 340

DOI:10.3969/j.issn.0255-7797.2005.03.020     

Li D Y , Jiang H F , Xiong G .

The reference point of matching ahapes

J of Math, 2005, 25 (3): 336- 340

DOI:10.3969/j.issn.0255-7797.2005.03.020     

李寿贵, 韩汇芳, 杨佩佩.

平行四边形内点与边界点的平均距离

武汉科技大学学报, 2011, 34 (5): 376- 380

DOI:10.3969/j.issn.1674-3644.2011.05.014     

Li S G , Han H F , YANG P P .

The average distance from the interior points to the boundary points of a parallelogram

Journal of Wuhan University of Science and Technology, 2011, 34 (5): 376- 380

DOI:10.3969/j.issn.1674-3644.2011.05.014     

李德宜, 杨佩佩, 李亭.

矩形的弦长分布

武汉科技大学学报, 2011, 34 (5): 381- 383

DOI:10.3969/j.issn.1674-3644.2011.05.015     

Li D Y , Yang P P , Li T .

Chord length distribution of rectangles

Journal of Wuhan University of Science and Technology, 2011, 34 (5): 381- 383

DOI:10.3969/j.issn.1674-3644.2011.05.015     

黄朝霞.

蒲丰投针问题研究

集美大学学报(自然科学版), 2005, 10 (4): 381- 384

DOI:10.3969/j.issn.1007-7405.2005.04.021     

Huang Z X .

Study on Buffon throwing problem

Journal of Jimei University (Natural Science), 2005, 10 (4): 381- 384

DOI:10.3969/j.issn.1007-7405.2005.04.021     

王媛媛, 李德宜, 李雪婵, .

正三棱柱的Buffon投针问题

数学杂志, 2013, 33 (5): 887- 890

URL    

Wang Y Y , Li D Y , Li X C , et al.

Buffon needle problem in regula tri-prism

J of Math, 2013, 33 (5): 887- 890

URL    

邹明田, 李寿贵, 陈莉莉.

一类特殊网格的几何概率

数学杂志, 2014, 34 (2): 374- 378

URL    

Zou M T , Li S G , Chen L L .

A special class of the geometric probability

J of Math, 2014, 34 (2): 374- 378

URL    

Ren D L . Topics in Integral Geometry. Singapore: World Scientific, 1994

[本文引用: 8]

Zhang G Y .

Dual kinematic formulas

Transactions of the American Mathematical Society, 1999, 351 (3): 985- 995

DOI:10.1090/S0002-9947-99-02053-X      [本文引用: 1]

Santalo L A . Integral Geometry and Geometric Probability. London: Addidion-Wesley, 1976

[本文引用: 9]

Merca M .

New convolutions for complete and elementary symmetric functions

Integral Transforms and Special Functions, 2016, 27 (12): 965- 973

DOI:10.1080/10652469.2016.1233405      [本文引用: 1]

Macdonald I G . Symmetric Functions and Hall Polynomials. 2nd ed Oxford: Clarendon Press, 1995

[本文引用: 1]

曾春娜, 姜德烁.

关于凸集平均曲率积分的注记

重庆师范大学学报(自然科学版), 2013, 30 (5): 62- 65

URL     [本文引用: 2]

Zeng C N , Jiang D S .

Some notes on mean curvature integral of convex sets

Journal of Chongqing Normal University (Natural Science), 2013, 30 (5): 62- 65

URL     [本文引用: 2]

曾春娜, 柏仕坤.

关于凸集的平均曲率积分的不等式(Ⅱ)

重庆师范大学学报(自然科学版), 2017, 34 (6): 57- 60

URL    

Zeng C N , Bo S K .

Some inequalities on mean curvature integral of convex sets (Ⅱ)

Journal of Chongqing Normal University (Natural Science), 2017, 34 (6): 57- 60

URL    

Zeng C N , Ma L , Xia Y W .

On mean curvature integrals of the outer parallel body of the projection of a convex body

Journal of Inequalities and Applications, 2014, Article number, 415

[本文引用: 1]

周家足.

\mathbb{S}.4中的Willmore泛函与包含问题

中国科学, 2007, 37 (2): 249- 256

URL    

Zhou J Z .

Willmore functional and inclusion problems in \mathbb{S}.4

Science in China, 2007, 37 (2): 249- 256

URL    

Zhang G Y, Zhou J Z. Containment Measures in Integral Geometry//Grinberg E, Li S G, Zhang G Y, Zhou J Z. Integral Geometry and Convexity. Singapore: World Scientific, 2005: 153-168

[本文引用: 1]

John G .

On Wallis formula

The American Mathematical Monthly, 1956, 63 (9): 643- 645

URL     [本文引用: 4]

Guo B N , Q F .

On the increasing monotonicity of sequence originating from computation of the probability of intersecting between a plane couple and a convex body

Turkish Journal of Analysis and Number Theory, 2015, 3 (1): 21- 23

URL     [本文引用: 1]

/