International Scholarly Research Notices

International Scholarly Research Notices / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 604074 |

Cherif Bendjaballah, "Statistical Properties and Applications of Correlated Length Intervals of Some Stationary Random Point Processes", International Scholarly Research Notices, vol. 2012, Article ID 604074, 7 pages, 2012.

Statistical Properties and Applications of Correlated Length Intervals of Some Stationary Random Point Processes

Academic Editor: J. Li
Received15 Nov 2011
Accepted18 Dec 2011
Published23 Jan 2012


The point process, a sequence of random univariate random variables derived from correlated bivariate random variables as modeled by Arnold and Strauss, has been examined. Statistical properties of the time intervals between the points as well as the probability distributions of the number of points registered in a finite interval have been analyzed specifically in function of the coefficient of correlation. The results have been applied to binary detection and to the transmission of information. Both the probability of error and the cut-off rate have been bounded. Several simulations have been generated to illustrate the theoretical results.

1. Introduction

It is known that the detection of an optical field at a low level of power is a sequence of events which is a set of distinct time instants {𝜗𝑗},0𝑗<, such that 𝜗𝑗+1𝜗𝑗,forall𝑗. These {𝜗𝑗}𝑠, which are the time instants of interaction between the photons and the detector device, for example, a photomultiplier, constitute a random point process (RPP) for which a positive instantaneous density 𝜆(𝜗𝑗) can be defined. We here consider only simple and homogeneous stationary processes.

Reciprocally, the existence of a physical optical field, given the knowledge of all of the properties of such an RPP, is not an easy problem to solve. For example, given an RPP of nonclassical properties [1] does not necessarily correspond to a nonclassical physical optical field, although its feasibility can be demonstrated [2]. This problem, mainly due to the quantum nature of the interaction between radiation and matter is not examined here.

Nevertheless, using a parameterized RPP, whose special values of such a parameter (denoted here 𝑐) correspond to a physical field, for example, a coherent or a thermal state, then it is reasonable to admit that intermediate values of 𝑐 corresponds to realistic optical fields. Despite these limitations, the properties of the RPPs studied here are important in statistical optics.

As already pointed out, there are two types of processing that can be utilized to characterize such RPP: the time interval distributions (TIDs) and the probability of number distributions (PNDs) [3]. If we choose to characterize the RPP by the TID, we can define the time interval between points 𝜃𝑗=𝜗𝑗+1𝜗𝑗, called residual time (or lifetime sometimes). Its probability distribution function (PDF) 𝑤(𝜃) will be called here a triggered PDF. When 𝜗𝑗=𝑡0 is arbitrary, that is to say not a point of the (RPP), the 𝑣(𝜃), PDF of the corresponding 𝜃𝑗 will be called a relaxed PDF.

If we now choose to characterize the RPP by the PND, we can as above, define the relaxed PND 𝑃(𝑛,𝑡) of the random variable (RV) 𝑁(𝑡0,𝑡0+𝑡) which is the number of instants occurring within the time interval [𝑡0,𝑡0+𝑡],𝑡0and𝑡𝜗𝑗,forall𝑗. As previously we may as well define 𝑄(𝑛,𝑡) a triggered PND of the RV 𝑁(𝜗𝑗,𝜗𝑗+𝑡),𝑡𝜗𝑗, the counting process being started by a point of the RPP.

The purpose of this paper is firstly to calculate the statistical properties of the RPPs, both in terms of the TIDs and the PNDs. Secondly we apply the results to the calculation of the performances, namely, the probability of error in binary detection and the cut-off rate in the binary transmission of the information, and their variations with respect to the coefficient of correlation 𝑐. The one-dimensional processes, such as the Poisson and the geometric PNDs, are finally utilized for establishing bounds to the performances.

In Section 2, the notations are defined and the basic equations are briefly recalled. Section 3 establishes the main results of the PDFs both those dealing with TIDs and PNDs depending on extreme values of 𝑐: low values and high values. Finally, Section 4 presents theoretical results and curves that explain and illustrate the results of the numerical simulation.

2. Basic Equations

We need to define the Laplace transform of 𝑓(𝑡)𝜌(𝑠)=𝑓(𝑡)0e𝑠𝑡𝑓(𝑡)d𝑡(1) and the functions𝜌1(𝑠)=𝜌(𝑠),𝜌𝑠(1𝜌(𝑠))2𝜌(𝑠)=(𝑠)(1+𝜌(𝑠)).𝑠(1𝜌(𝑠))(2) The PDF of the number of events registered between 0 and 𝑡, denoted by PND, reads𝑃𝑡(𝑛;𝑡)=Prob𝑡𝑇=𝑛,𝑡𝑛+1𝑡=𝖤𝑛<𝑡<𝑡𝑛+1,(3) where is the indicator function (𝑡𝑇)=1(0), in case the event occurs or not within 𝑡𝑇. Therefore𝑃(𝑛;𝑡)=1𝑃(𝑛;𝑠)=11𝜌(𝑠)𝑠𝜌(𝑠)𝑛,(4) where the symbol 1 denotes the inverse Laplace transform. From (4), it is easily seen that (see e.g., [4, 5])𝖤[𝑁](𝑡)=1𝜌1𝖤𝑁(𝑠),2(𝑡)=1𝜌2(𝑠).(5) The symbol 𝖤 denotes the mathematical expectation. In the following 𝑃(𝑛;𝑡) will be called the relaxed PDF. Another interesting PDF that can be derived from (4),𝑄(𝑛;𝑡)=𝑛+1𝖤[𝑁]𝑃(𝑛+1;𝑡),(6) will be called the triggered PDF. In fact, it is known [3, 6] that these PDFs can be expressed using the moments of the time-integrated density𝒥(𝑡)=𝑡0𝑡+𝑡0𝜆(𝜗)d𝜗,(7)𝒥(𝑡)=𝜆𝑡 for the Poisson RPP. Therefore1𝑃(𝑛;𝑡)=𝖤𝒥𝑡𝑛!0+𝑡𝑛e𝒥(𝑡0+𝑡),1𝑄(𝑛;𝑡)=𝖤[𝑁]1𝖤𝒥𝑡𝑛!0+𝑡𝑛+1e𝒥(𝑡0+𝑡).(8) For all RPPs that we deal with here, all the PDFs depend only on 𝑡 (not on 𝑡0) due to the stationary property. When 𝑡 is a parameter, 𝑃(𝑛;𝑡) and 𝑄(𝑛;𝑡) are simply denoted 𝑃(𝑛) and 𝑄(𝑛).

In terms of TIDs, we just recall the basic formulas of 𝑤(𝑡) and 𝑣(𝑡) (see [6])1𝑤(𝑡)=𝖤[𝒥]𝖤𝒥𝑡0𝒥𝑡0e+𝑡𝒥(𝑡0+𝑡),𝒥𝑡𝑣(𝑡)=𝖤0e+𝑡𝒥(𝑡0+𝑡).(9) Notice that 𝑔(0)𝑤(𝑡=0)/𝑣(𝑡=0) the ratio of the values at 𝑡=0 is related to 𝜎2𝒥, the variance of 𝒥, such that 𝑔(0)1=𝜎2𝒥/𝖤[𝒥]2, where 𝜎2𝒥 is the variance of the time-integrated density.

3. Model for Correlated Variables

Among the several models proposed to deal with correlated variables, we consider the Arnold-Strauss model [7, 8]𝑓(𝑡,𝑢)=𝐾e(𝑎𝑡+𝑏𝑢+𝑐𝑡𝑢),(10) where 𝑎,𝑏,𝑐>0. The time intervals 𝑡,𝑢 are of course positive RVs in [0,[. The constant of normalization is𝐾=𝑐e𝑎𝑏/𝑐𝐾Ei(1,𝑎𝑏/𝑐)𝑐=𝑎𝑏+𝑐2𝐾𝑎𝑏for𝑐1,=𝑐1𝛾log(𝑎𝑏)for𝑐1,(11) where Ei(𝑚,𝑥) are the exponential integral functions [9] for 𝑚=1,2,, and 𝑥+. They verify the approximate expressions given in (A.1a)–(A.2a) in Appendix A.1. The constant 𝛾=0.5772157 is the Euler constant.

The marginal distribution of 𝑡 is deduced from (10) and is given by𝑓(𝑡)=0e𝑓(𝑡,𝑢)d𝑢=𝐾𝑎𝑡.𝑏+𝑐𝑡(12) In the following, we report some calculations that can be obtained under closed forms depending on the value of the parameter 𝑐. Let us first consider the case where 𝑐1.

3.1. 𝑐1

In what follows, most of the calculations are done up to 𝑂(𝑐3).

3.1.1. TID

It can be shown that (see, e.g., [6])𝑤e(𝑡)=𝐾𝑎𝑡,𝑏+𝑐𝑡𝑣(𝑡)𝜆𝑡𝑤(𝜃)d𝜃.(13) Given 𝖤[𝑡], the first moment of 𝑡, the average value of the density of the process 𝜆=1/𝖤[𝑡] is given by 𝑐𝜆𝑎+𝑏𝑐22𝑎𝑏2.(14) Equation (13) yields 𝜈(0)=𝜆 and d𝜈/d𝑡=𝜆𝑤(𝑡) which prove that 𝜈(𝑡) is a monotonic continuously decreasing function of 𝑡. More precisely, we can show that𝑤(𝑡)𝑎3𝑏2𝑏𝑐𝑡+𝑐2𝑡2𝑎2𝑏2𝑎𝑏𝑐+2𝑐2e𝑎𝑡,𝑎(15)𝜈(𝑡)𝑎2𝑏2𝑎𝑏𝑐+2𝑐2𝑎𝑐(𝑎𝑏2𝑐)𝑡+𝑎2𝑐2𝑡2𝑎2𝑏22𝑎𝑏𝑐+6𝑐2e𝑎𝑡.(16) From (15), we deduce𝜎2𝑡1𝑎212𝑐+𝑎𝑏9𝑐2𝑎2𝑏2.(17) On the other hand, we can see that𝑡𝐾(0,𝑐1)=𝑤(0)𝑣(0)=𝑏𝑐𝜆2𝑎𝑏2,(18) which positivity is a characteristic of classical processes [10].

3.1.2. PND

For a simple approximation, at a first order of 𝑐, we may use the following:1𝑓(𝑡)𝑎𝑏𝑎e𝑎𝑡𝑐𝑎2𝑏2𝑎2𝑡e𝑎𝑡,𝑃(𝑛,𝑡)(𝑎𝑡)𝑛e𝑛!𝑎𝑡𝑐𝑎2𝑏2×(𝑎𝑡)2𝑛+1(2𝑛)!2(𝑎𝑡)2𝑛+1+(2𝑛+1)!(𝑎𝑡)2𝑛1e(2𝑛1)!𝑎𝑡,(19) where 𝑃(0,𝑡)=(1+𝑎𝑡/2)e𝑎𝑡 (see [1] pages 789–792).

When only a few values of 𝑃(𝑛,𝑡) or only the moments are needed, it is however, better, to use (3) and (5). Based on their expansion as series in 𝑠, we obtain, at the second order of 𝑐,𝖤[𝑁]𝑎2𝑎𝑎𝜔𝑡+22(𝑎𝜔)2𝜔2𝑡2,𝖤𝑁2𝑎2𝑎𝑎𝜔𝑡+22𝑎2+𝜔2(𝑎𝜔)2𝜔2𝑡2,𝜎2𝑛1𝑎+𝜔+2+1𝑎𝜔𝑎2,(𝑡=1),(20) where we denote 𝜔=𝑐/𝑏.

Seeking the exact expressions of the PDFs of the number seems difficult to obtain. However, as just seen, approximations of closed expressions are simple. Thus, using (4), it is easy to calculate approximate expressions of the PDFs of the number for 𝑛=0,1,2 and 𝑡=1, as given in (A.3) in Appendix A.2. Therefore, we can prove that𝑛(0,𝑐1)=𝑃(0)𝑃(1)e𝑎1𝑎2𝑎𝑏=𝑐,𝑎=1,𝑏=11𝑒2𝑐𝑐3𝑒2,(21) which is negative forall𝑐1. We will numerically illustrate this property by simulations. When we choose the triggered processing, we obtain𝑘𝑛1(0,𝑐1)=𝑄(0)𝑄(1)𝑒1𝑐𝑐4𝑒2,(22) which is here again negative forall𝑐.

The case 𝑐1 (but finite) is perhaps more interesting although closed forms of the moments are difficult to obtain.

3.2. 𝑐1

The calculations are now done up to 𝑂(1/𝑐3).

3.2.1. TID

We can show that the normalized 𝑤(𝑡) is given by𝑤(𝑡)𝑐e𝑎𝑏/𝑐eEi(1,𝑎𝑏/𝑐)𝑎𝑡𝑐𝑡+𝑏𝑡0e𝑎𝑡𝑏𝑡(𝜒log𝑐)1𝑐1𝑎+𝑡,(23) and the unnormalized PDF, evaluated up to 𝑂(1/𝑐2) and 𝑂(𝑡),𝑣(𝑡)𝜆Ei(1,𝑎(𝑡+𝑏/𝑐)),Ei(1,𝑎+𝑏/𝑐)0<𝑡1𝜆Ei(1,𝑎𝑡)𝑏log(𝑐/𝑎𝑏)𝛾𝜆𝑐1elog(𝑎𝑏/𝑐)+𝛾𝑎𝑡𝑡𝜆𝑎𝑏𝑐Ei(1,𝑎𝑡)(log(𝑎𝑏/𝑐)+𝛾)2,𝜆𝛾+log(𝑎𝑡),log(𝑎𝑏/𝑐)𝛾(24) leading to𝑡𝐾(0,𝑐1)=𝑤(0)𝑣(0)=𝑏𝑐𝜆𝑏(1𝜒)𝑎log𝑐,(25) where we denote 𝜒=𝛾+log(𝑎𝑏) and assume that 𝜆𝑎(log𝑐𝜒). It is seen that 𝑡(0,𝑐1) given by (25) is now positive and increases more slowly than 𝑡(0,𝑐1)𝑐2 given by (18). From (23), we can also deduce that𝖤[𝑡]1𝑎1,𝜎log𝑐𝜒2𝑡1𝑎21=𝖤[𝑡]log𝑐𝜒𝑎.(26) Both moments tend to 0 when 𝑐.

3.2.2. PND

The approximate expressions given in (A.4) in Appendix A.2  lead to𝑛(0,𝑐1)=𝑃(0)𝑃(1)1+𝛾log𝑐𝑐.(27) We can calculate the first two moments of the number for the specific case 𝑎=1 and 𝑏=1, 𝖤[𝑁]𝜎log(𝑐),(28a)2𝑛[𝑁]𝖤2.(28b)As for the previous case, it can be shown that 𝑘𝑛(𝑐)0,forall𝑐.

4. Simulation and Results

We have used the algorithm recently described [11] for several values of 𝑐.

In Figure 1, results of the simulated data of the TIDs are plotted. The 𝑤(𝑡), the triggered PDF follows very well (23). The 𝑣(𝑡), the relaxed one has been well fitted by the expression𝑣(𝑡)=𝜈0𝜆Ei(1,𝑡)Ei(1,1+1/𝑐)0<𝑡1𝜈1𝜆log𝑡,log(1/𝑐)(29) for 𝑎=1,𝑏=1,𝑐=30 and where 𝜈0=0.07 and 𝜈1=0.2.

By the way, it is interesting to remark that the PDF of 𝑆=𝑗=0𝜃𝑗,1, which is the addition of several correlated identically distributed positive random time interval, deviates from the Gaussian profile. We may conclude that the application of the central limit theorem requires the addition of a very high number of correlated random variables.

On the other hand, for high values of 𝑐, the theoretical approximations of the PNDs are given by1𝑃(𝑛)=𝒩𝑏1e𝑎1𝑛+𝑏2e𝑎2𝑎𝑛2,𝑛!(30) where𝒩=𝑏2+𝑏11e𝑎1,𝖤[𝑁]=1(31a)𝒩𝑎2𝑏2+𝑏1e𝑎1(e𝑎1)12,(31b) leading to𝑄(𝑛)=𝑛+1𝖤[𝑁]𝑃(𝑛+1).(32) The values of the parameters that correctly fit the simulated results of the PDFs, as seen in Figure 2, are 𝑏1=0.2,𝑎1=1,𝑏2=0.62,𝑎2=5.2,and𝖤[𝑁]=3.64.

The variations of the moments with respect to 𝑐 are derived from the expressions recalled in the Section 2. We obtained𝖤[𝑁]11+𝛾+2𝜎log𝑐,(33a)2𝑛11+𝛾+4(1+log𝑐)2.(33b)Regarding the reduced moments 𝑛 and 𝑘𝑛, an excellent fit of the simulated results of Figure 3 is carried out with 𝛼1=4,𝛽1=0.17 and 𝛼2=3,𝛽1=0.24𝑛1𝛾2+12𝛼𝛾log1𝑐𝛼1𝑐𝛽1,𝑘(34a)𝑛1𝛾2+12𝛼𝛾log2𝑐𝛼2𝑐𝛽2,(34b)leading to 𝑛𝑐2/2log𝑐 and 𝑘𝑛𝑐4/24(log𝑐)3 which are only in qualitative agreement.

As an application of these results to communications, we consider a system of communication processed with a direct threshold detector. The decision device operates such that𝑛𝑛𝑠(𝑐)H0𝑃𝑐=0(𝑛)=e𝑎𝑎𝑛,𝑛!(35a)𝑛>𝑛𝑠(𝑐)H1𝑃𝑐0(𝑛),(35b)based on the binary hypotheses: H0 (no correlation) and H1 (correlation with the parameter 𝑐). The threshold is obtained from the likelihood ratio Λ(𝑛),𝑃Λ(𝑛)=𝑐0(𝑛)𝑃𝑐=0(𝑛)>1𝑛𝑠(𝑐).(36) We can also utilize the TIDs as a useful tool of processing [12]. The threshold would be based on the likelihood ratio Θ(𝑡)=𝑤𝑐(𝑡)/𝜆e𝜆𝑡>1𝑡𝑠(𝑐) and the decision would operate as follows𝑡𝑡𝑠(𝑐)H1𝑤𝑐e(𝑡)=𝐾𝑎𝑡,𝑏+𝑐𝑡𝑡𝑡𝑠(𝑐)H0𝑎e𝑎𝑡,(37) where 𝑡𝑠(𝑐)=(1/c)(𝐾/a𝑏), 𝐾 being given by (11).

Here, we focus on the method based on PNDs because it is generally more efficient. To simplify the calculations, the decision is not randomized [13].

Now, the probability of error in detection when processing with the relaxed PND is given by𝑃err=121𝑛𝑠𝑛=0𝑃𝑐=0(𝑛)+𝑛s𝑛=0𝑃𝑐0.(𝑛)(38) Similarly, we have for the processing with the triggered PND𝑄err=121𝑛𝑠𝑛=0𝑃𝑐=0(𝑛)+𝑛𝑠𝑛=0𝑄𝑐0.(𝑛)(39) It is first seen that 𝑄err𝑃err, the inequality which is demonstrated for a special case in the Appendix B. Furthermore, both probabilities of error in detection decrease with 𝖤[𝑁], then with 𝑐 because 𝖤[𝑁]log𝑐 as it is calculated in (28a), (33a) and is shown in Figure 3. This is within the range of 𝑐0𝑐100. For instance, with the help of (30)–(31b) and (38)-(39) calculated for 𝑛𝑠=1 and 𝑐=100, we obtained 𝑃err0.28 and 𝑄err0.16, the values which are in excellent agreement with the simulated results of Figure 4.

On the other hand, the bounds to 𝑄err and 𝑃err can easily be calculated. In fact, denoting𝖤[𝑁]𝜈=[𝑁]1+𝖤,(40) where 𝖤[𝑁]=𝑎+𝜉. Now, taking into account (35a), we have up to 𝑂(𝜉2),𝑃1(𝑛)=(1𝜈)𝜈𝑛,𝑃(1)err11𝑒𝜈2211𝑒𝜉8,(41) for 0𝜈1. Therefore 𝑃(1)err is the approximation of the exact 𝑃err for 𝑛𝑠=1 plotted as the curve quoted “1” in Figure 4 where 𝑎=1. Similarly, taking into account (35a), we have𝑃1(𝑛)=e𝖤[𝑁]𝖤[𝑁]𝑛,𝑃𝑛!(2)err121𝑒+[𝑁]1+𝖤2e𝖤[𝑁]121𝑒𝜉.(42) Again, 𝑃(2)err is the approximation of the exact 𝑃err for 𝑛𝑠=1 plotted as the the curve quoted “2” in Figure 4 for 𝑎=1. Finally, we have𝑃(2)err𝑄err𝑃err𝑃(1)err.(43) Let us conclude this analysis with a brief comment on information. We will concentrate on the cut-off rate which is known as a useful criterion for evaluating the performances of a channel. Thus, for a binary noiseless channel, when the transmission of messages “0” and “1” is done via the probabilities 𝑃𝑐=0 and 𝑃𝑐0 respectively (the probabilities a priori are taken equal to 1/2) the cut-off rate expresses as2𝑅=log1+𝑛=0𝑃𝑐=0(𝑛)𝑃𝑐0,(𝑛)(44) which is generally interpreted as a lower bound to the channel capacity [14]. The 𝑃𝑐=0(𝑛) is the Poisson PND of parameter 𝑎=1 and 𝑃𝑐0(𝑛) will be, as above, either the relaxed or the triggered PND yielding the cut-off rates 𝑅𝑝 and 𝑅𝑞, respectively. Because exact closed expressions seem difficult to attain, the bounds are very useful and can easily be established. The first bound is obtained using geometric PNDs𝑃0(𝑛)=1𝜈0𝜈𝑛0𝑅,(45a)12=log1+1𝜈0(1𝜈)/1𝜈0𝜈0𝜉2𝜉23213𝜉4+𝜉22,(45b)where 𝜈 is given by (40) and 𝜈0=𝑎/(1+𝑎). This is a lower bound to 𝑅𝑝.

The second bound is obtained using the Poisson PNDs𝑅22=log1+e1/2(𝖤[𝑁]𝑎)20𝜉2𝜉2𝜉1612+9𝜉2,32(46) which is an upperbound to 𝑅𝑝.

In Figure 5, the results of the simulation 𝑅𝑝 and 𝑅𝑞 and the bounds (curves quoted “1” and “2”) given by (45b) and (46) are plotted versus 𝖤[𝑁]. Here again, it is seen that the processing with the triggered PND performs much better that the processing with the relaxed PND𝑅1𝑅𝑝𝑅2𝑅𝑞,(47) in the limit of 0𝑐100. However, because the cut-off rate is a monotonic increasing function with respect to 𝖤[𝑁], the inequalities (47) may be extrapolated to for all 𝑐.

In conclusion, the binary performances, as summarized by the inequalities (43) and (47), show that the processing with the triggered PND is the preferable mode of operation. Both performances, in detection and information transmission, are improved with the coefficient of correlation.


A. Approximate Expressions

A.1. The Exponential Integral Functions

The exponential integral functions can be expressed, up to 𝑂(𝑥3), 𝑥Ei(1,𝑥)𝛾log𝑥+𝑥24,𝑥(A.1a)Ei(2,𝑥)1+(𝛾+log𝑥1)𝑥22,1(A.1b)Ei(3,𝑥)21𝑥21𝛾+23log𝑥4𝑥2,(A.1c)where 𝛾 is the Euler constant. For 𝑥eEi(1,𝑥)𝑥𝑥11𝑥,e(A.2a)Ei(2,𝑥)𝑥𝑥21𝑥e(A.2b)Ei(3,𝑥)𝑥𝑥31𝑥.(A.2c)

A.2. Probability of Number Distribution

Depending on the values of 𝑐, we have the following approximations of the first values of 𝑃(𝑛):𝑃(0)e𝑎11𝜔𝑎𝜔2,𝑃(1)e𝑎1𝑎+(1𝑎)𝜔+𝑎+𝑎632𝜔2,𝑃(2)e𝑎𝑎22𝑎+𝑎12𝑎𝜔+28+32𝜔2𝑎2,(A.3) for 𝑐1 and where 𝜔=𝑐/𝑏. Similarly, it can be shown that𝑃(0)1+𝛾log𝑐𝑐,𝑃(1)log𝑐𝛾𝑐,𝛾𝑃(2)22(1+𝛾)log𝑐𝑐2,(A.4) for 𝑐1.

B. Inequality between Probabilities of Error in Binary Detection

To demonstrate that the triggered processing yields better performance than the relaxed processing, Δ𝑒=𝑄err𝑃err0, we begin with (8)1𝑃(𝑛;𝑡)=𝖤𝑛!𝒥(𝑡)𝑛e𝒥(𝑡),1𝑄(𝑛;𝑡)=𝖤[𝑁]1𝖤𝑛!𝒥(𝑡)𝑛+1e𝒥(𝑡),(B.1) and show that for𝑃err=121𝑛𝑠𝑛=0𝑃𝑐=0(𝑛)+𝑛𝑠𝑛=0𝑃𝑐0,𝑄(𝑛)err=121𝑛𝑠𝑛=0𝑃𝑐=0(𝑛)+𝑛𝑠𝑛=0𝑄𝑐0,(𝑛)(B.2) we haveΔ𝑒=𝑄err𝑃err,=12𝑛𝑠𝑛=0𝖤e𝒥𝒥𝑛+1𝖤[𝒥]𝒥𝑛=12𝖤𝑛𝑠𝑛=0e𝒥𝒥𝑛+1𝖤[𝒥]𝒥𝑛12𝖤𝑛𝑠𝑛=0𝒥𝑛+1𝖤[𝒥]𝒥𝑛=12𝑛𝑠𝑛=0𝖤𝒥𝑛+1𝖤[𝒥]𝒥𝖤𝑛.(B.3) where we used e𝒥1,𝒥0.

For 𝑛𝑠=1, noticing that 𝖤[𝒥2]𝖤[𝒥]2, we have Δ𝑒0.

For higher values of 𝑛𝑠, this method seems not useful because it requires to prove that 𝖤[𝒥𝑛+1]/𝖤[𝒥]𝖤[𝒥𝑛],forall𝑛, which is not so easy although the inequality is true for several types of density distributions of interest in statistical optics.


Laboratoire des Signaux et Systèmes is a joint laboratory (UMR 8506) of CNRS. and École Supérieure d'Électricité is and associated with the Université Paris-Orsay, France.


  1. C. Bendjaballah, “Nonclassical second-order photon statistics,” Journal of Physics B, vol. 39, no. 4, pp. 783–803, 2006. View at: Publisher Site | Google Scholar
  2. C. Bendjaballah, “Nonclassical photon statistics,” Physical Review A, vol. 73, no. 5, Article ID 053816, 2006. View at: Publisher Site | Google Scholar | MathSciNet
  3. C. Bendjaballah, Introduction to Photon Communication, Springer, Heidelberg, Germany, 1995.
  4. J. A. McFadden, “On the lengths of intervals in a stationary point process,” Journal of the Royal Statistical Society B, vol. 24, pp. 364–382, 1962. View at: Google Scholar
  5. C. Godrèche and J. M. Luck, “Statistics of the occupation time of renewal processes,” Journal of Statistical Physics, vol. 104, no. 3-4, pp. 489–524, 2001. View at: Publisher Site | Google Scholar | MathSciNet
  6. C. Bendjaballah, “Statistical properties of several models of fractional random point processes,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, pp. 3224–3236, 2011. View at: Google Scholar
  7. B. C. Arnold and D. Strauss, “Bivariate distributions with exponential conditionals,” Journal of the American Statistical Association, vol. 83, pp. 522–527, 1988. View at: Google Scholar
  8. S. Nadarajah and S. Kotz, “Reliability for some bivariate exponential distributions,” Mathematical Problems in Engineering, vol. 2006, pp. 1–14, 2006. View at: Publisher Site | Google Scholar
  9. I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, New York, NY, USA, 6th edition, 2000.
  10. B. Picinbono and C. Bendjaballah, “Characterization of nonclassical optical fields by photodetection statistics,” Physical Review A, vol. 71, no. 1, Article ID 013812, 12 pages, 2005. View at: Publisher Site | Google Scholar
  11. Y. Yu, “Efficient simulation of a bivariate exponential conditionals distribution,” Computational Statistics and Data Analysis, vol. 52, no. 5, pp. 2273–2276, 2008. View at: Publisher Site | Google Scholar
  12. C. Bendjaballah and K. Hassan, “Probability of detecting a coherent optical signal in thermal noise,” Journal of the Optical Society of America A, vol. 73, pp. 1840–1843, 1983. View at: Google Scholar
  13. C. W. Helstrom, Quantum Detection and Estimation Theory, Academic Press, New York, NY, USA, 1976.
  14. J. L. Massey, “Capacity, cut-off rate and coding for a direct-detection optical channel,” IEEE Transactions on Communications, vol. 29, no. 11, pp. 1615–1621, 1981. View at: Google Scholar

Copyright © 2012 Cherif Bendjaballah. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.