ClosedForm Bounds to the Rice and Incomplete Toronto Functions and Incomplete LipschitzHankel Integrals
Abstract
This article provides novel analytical results for the Rice function, the incomplete Toronto function and the incomplete LipschitzHankel Integrals. Firstly, upper and lower bounds are derived for the Rice function, . Secondly, explicit expressions are derived for the incomplete Toronto function, , and the incomplete LipschitzHankel Integrals of the modified Bessel function of the first kind, , for the case that is an odd multiple of and . By exploiting these expressions, tight upper and lower bounds are subsequently proposed for both function and integrals. Importantly, all new representations are expressed in closedform whilst the proposed bounds are shown to be rather tight. Based on these features, it is evident that the offered results can be utilized effectively in analytical studies related to wireless communications. Indicative applications include, among others, the performance evaluation of digital communications over fading channels and the informationtheoretic analysis of multipleinput multipleoutput systems.
Closedform representations, Rice function, Incomplete Toronto function, Incomplete LipschitzHankel Integrals, Marcum Qfunction, lower and upper bounds
1 Introduction
It is widely accepted that special functions constitute invaluable mathematical tools in the majority of fields in natural sciences and engineering. In the area of telecommunications, their utilization in various studies studies often allows the derivation of analytic expressions for important performance measures such as for example error probability and channel capacity. Furthermore, it has been shown that the computational realization special functions is not generally laborious since the most of them are included as builtin functions in popular scientific mathematical packages such as , and . To this effect, both the algebraic representation and computation of any derived analytical expressions have been significantly simplified.
Among others, the Rice function, the incomplete Toronto function and the incomplete LipschitzHankel integrals (ILHIs) appear in analytical solutions of numerous problems in telecommunications. They were all proposed several decades ago and they are denoted as , and , respectively.
The Rice function was firstly proposed by S. O. Rice in [1] and has been largely exploited in the study of zero crossings, in the analysis of angle modulation systems, in radar pulse detection and in error rate analysis of differentially encoded systems, [2][5]. Its definition is typically given in integral form which involves an exponential term and a modified Bessel function of the first kind and zeroth order. Alternative representations include two infinite series and one closedform expression  reported in [4] and [5], respectively. The former, are expressed in terms of the modified Struve function and the modified Bessel function of the first kind whereas the latter is expressed in terms of the Marcum Qfunction, , [6][10].
In the same context, the incomplete Toronto function constitutes a special case of the Toronto function, which was initially proposed by Hatley in [11]. It also includes as a special case the Marcum Qfunction and has been used in studies related to statistics, signal detection and estimation, radar systems and error probability analysis, [12][14]. Its definition is also given in integral form while alternative representations include two infinite series, which were proposed in [15].
Finally, the ILHIs are a class of incomplete cylindrical functions that have been largely encountered in analytical solutions of numerous problems in electromagnetics, [16][17] and the references therein. In communication theory, they have been sufficiently utilized in investigations associated with the error rate analysis of MIMO systems under imperfect channel state information (CSI) that employ adaptive modulation, transmit beamforming and maximal ratio combining (MRC), [18].
However, in spite of the evident importance of the and functions and the integrals, it is noted that they are all neither tabulated, nor included as builtin functions in the aforementioned popular software packages. As a consequence, they are rather inconvenient to handle both analytically and computationally.
Motivated by these issues, this work is devoted in deriving novel representations for these functions and integrals. Specifically, upper and lower bounds to the function and explicit expressions and upper and lower bounds to the function and integrals, are derived for the case that and . Notably, the offered results are expressed in closedform and have a tractable algebraic form. As a result, they can be meaningfully utilized in various analytical studies associated to wireless communications such as the performance evaluation of digital communications over fading channels and the informationtheoretic analysis of MIMO systems.
2 Definitions and Existing Representations
2.1 The Rice function
The Rice function is defined as [4][5],
(1) 
where is the modified Bessel function of the first kind and zero order, [6][7]. An equivalent integral representation to was given in [4], namely,
(2) 
along with the following alternative series representations,
(3) 
and
(4) 
The notations and denote the modified Struve function and the gamma function, respectively [6][7]. Recalling [4], equation converges relatively quickly for the case that is large and is small. On the contrary, equation converges quickly when is small and is large. Nevertheless, this way of computation of function is rather inefficient due to the following three facts: a) two relationships are required; b) the above series are relatively unstable due to their infinite form; c) the function is not builtin in widely used mathematical software packages.
An adequate way of resolving this issue was reported in [5]. There, the Rice function is related to the Marcum Qfunction of the first order, by the following relationships,
(5) 
and
(6) 
where, and .
2.2 The Incomplete Toronto Function
The incomplete Toronto function is defined as,
(7) 
Importantly, for the special case that , it can be equivalently expressed in terms of the Marcum Qfunction by the following relationship,
(8) 
Alternative representations to the function, in the form of infinite series, were reported in [15], namely,
(9) 
and
(10) 
where the notations and denote the Pochhammer symbol and the lower incomplete gamma function, respectively [6][7]. Also,
and
Although equation is exact, its algebraic representation is rather inconvenient both analytically and numerically. Equation is significantly more tractable than , yet its infinite form eventually raises convergence and truncation issues.
2.3 The Incomplete LipschitzHankel Integrals
The general ILHI is defined as,
(11) 
where ,,, may be complex [16][17]. The notation denotes one of the cylindrical functions
(12) 
where the set of coefficients and can be obtained recursively, [18]. As aforementioned, the above relationship was found useful in analytical investigations related to error rate of MIMO systems under imperfect channel state information (CSI).
3 Performance Bounds to the Rice function
This section is devoted to the derivation of upper and lower bounds for the Rice function. To this end, it is a critical to express function alternatively.
3.1 An Alternative Representation to the Function
Lemma 1. For and , the following relationship holds,
(13) 
Proof. By integrating once equation by part, it follows that,
(14) 
According to the basic principles of integration, the first integral yields straightforwardly,
(15) 
Subsequently, based on [6][7], the derivative of the modified Bessel function of the first kind and order is rewritten as,
(16) 
which for reduces to
(17) 
Consequently, by setting , substituting and into and noticing that , equation is deduced.
3.2 An Upper Bound to the Function
Theorem 1. For and , the following inequality holds,
(18) 
Proof. It is recalled here that the modified Bessel function of the first kind is strictly decreasing with respect to its order, . Therefore, for it immediately follows that . Thus, for , one obtains the following inequality,
(19) 
Notably, for the special case that is an odd multiple of , i.e. , a closedform representation for the function is given in , namely,
(20) 
Therefore, for , eq. reduces to,
(21) 
Evidently, the proof of the theorem is subject to evaluation of the integral in . To this end, by substituting into , one obtains,
(22) 
which has the following closedform solution,
(23) 
where,
is the error function, [6][7]. Finally, by substituting into , the proof is completed.
Remark. The authors in [19] derived closedform bounds to the Marcum Qfunction. Thus, by making the necessary change of variables and make the according substitution in equations and/or , an alternative expression to can be deduced. However, this expression is significantly less compact and convenient than both analytically and numerically.
3.3 A Lower Bound to the Function
Theorem 2. For and , the following inequality holds,
(24) 
where is the one dimensional Gaussian Qfunction [6],
Proof. According to the aforementioned monotonicity property of the function it follows that . Therefore, by making the necessary change of variables and substituting in , one obtains,
(25) 
Importantly, a similar inequality may be also deduced by exploiting equations and . To this end, it is firstly recalled that the Marcum Qfunction is strictly increasing with respect to its order . Based on this, it follows that
(26) 
Subsequently, by substituting into , the following inequality is deduced,
(27) 
Of note, the authors in [19] show that,
(28) 
Therefore, by substituting into , eq. is obtained and thus the proof is completed.
Remark. A lower bound to could be theoretically derived by following the same methodology as in Theorem . Nevertheless, this approach ultimately renders a representation that is both complex and divergent.
3.4 Numerical Results
The behaviour and tightness of the derived bounds is illustrated in figures . In more details, figure depicts the bounds in and for along with results obtained from numerical integrations for comparative purposes. Evidently, the upper bound is tighter than the lower bound for small values of . However, its tightness degrades as increases while the lower bound becomes much tighter. According to figure , this appears to be also the case for a relatively small and different values of .
Interestingly enough, for large values of  typically  the lower bound in is so tight that it eventually becomes a highly accurate approximation to . This is clearly illustrated in figure where for and different , the plotted curve provides an excellent match to the corresponding theoretical results. This is also justified by the overall small absolute relative error, , which is over almost the whole range of values of .
4 A Closedform Representation and Bounds to the Incomplete Toronto Function
As mentioned in section , the is neither expressed in terms of other elementary and/or special functions, nor is a builtin function in popular mathematical software packages. Motivated by this, we derive a closedform representation for the case that is an odd multiple of . Subsequently, we exploit this representation to propose novel closedform upper and lower bounds.
4.1 A Closedform Solution to the Function
Theorem 3. For , and , the following closed form relationship holds,
(29) 
where and , denote the upper and lower incomplete gamma functions, respectively.
Proof. By setting in , the function is rewritten as,
(30) 
By substituting into and utilizing the basic identity , it follows that,
(31) 
where . Evidently, a closedform solution to the above expression is subject to evaluation of the two involved integrals. To this end, with the aid of , the above expression is rewritten equivalently as follows,
(32) 
Finally, with the aid of , eq. is deduced and thus, the proof is completed.
4.2 Upper and lower bounds to the function
With the aid of Theorem , explicit bounds to the incomplete Toronto function may be straightforwardly deduced.
Corollary 1. For and , , the following inequality holds,
(33) 
where is given in closedform in eq..
Proof. The incomplete Toronto function is strictly decreasing with respect to . Therefore, for an arbitrary real positive value , it follows that . Hence, for the case that and , the lower bound in is straightforwardly deduced in closedform.
Corollary 2. For and , , the following inequality holds,
(34) 
where is given in closedform in eq..
Proof. The proof follows immediately from Theorem and Corollary .
4.3 Numerical Results
The validity of the closedform expression and the behaviour and tightness of the derived bounds are shown in figures and with respect to . In the former, eq. is depicted for and along with results obtained from numerical integrations for and , and . In the latter, eq. is depicted for and along with numerical results for , and . It is evident that is in exact agreement with the numerical results while the overall tightness of the derived bounds is quite adequate over the whole range of values of .
5 A closedform representation and performance bounds to the incomplete LipschitzHankel integrals
Likewise the function, the ILHIs are neither tabulated, nor builtin in popular mathematical software packages. However, their algebraic form renders possible the derivation of a closedform expression for the case that is an odd multiple of . By exploiting this result, upper and lower bounds are deduced based on the same criteria as in the previous section.
5.1 A ClosedForm Solution to the Integrals
Theorem 4. For , and , the following closedform relationship holds,
(35) 
where .
Proof. By substituting the function with its closedform representation in , equation can be equivalently rewritten as,
(36) 
Notably, the integrals in clearly belong to the family of gamma special functions. Therefore, after some basic algebraic manipulation and with the aid of , eq. is deduced.
Remark. The present analysis was limited in the consideration of only the function in . However, similar results may be also derived analogously for the Bessel functions , , and the Hankel functions , .
5.2 Upper and Lower Bounds to the Integrals
Corollary 3. For , and , the following inequality holds
(37) 
where the is given in in closedform.
Proof. The integrals are strictly decreasing with respect to . Therefore, for an arbitrary real positive value , it follows that . Thus, for the case that and , one obtains the the closedform lower bound in .
Corollary 4. For , and , the following inequality holds
(38) 
where is given in in closedform.
Proof. The proof follows immediately from Theorem and Corollary .
5.3 Numerical Results
The validity and behaviour of the offered results are explicitly illustrated in figure . Specifically, one can observe the exactness of and the rather tight performance of the performance bounds in all areas of values of their arguments.
6 Conclusion
In this work, explicit representations and performance bounds for the Rice function, the incomplete Toronto function and the incomplete LipschitzHankel integrals of the modified Bessel function of the first kind were derived. The offered results are novel and are all expressed in closedform. This property is sufficiently advantageous since it renders them suitable for efficient utilization in various analytical applications in the wide area of digital communications.
Footnotes
 Only the function is considered in the present analysis.
References
 S. O. Rice, Statistical properties of a sine wave plus random noise, Bell Syst. Tech. J., 1948, 27, pp. 109157
 J. H. Roberts, Angle Modulation, Stevenage, England: Peregrinus, 1977
 R. F. Pawula, S. O. Rice and J. H. Roberts, Distribution of the phase angle between two vectors perturbed by Gaussian noise, IEEE Trans. Commun. vol. COM30, pp. 18281841, Aug. 1982
 B. T. Tan, T. T. Tjhung, C. H. Teo and P. Y. Leong, Series representations for Rice’s function, IEEE Trans. Commun. vol. COM32, No. 12, Dec. 1984
 R. F. Pawula, Relations between the Rice function and the Marcum Qfunction with applications to error rate calculations, Elect. Lett. vol. 31, No. 24, pp. 20782080, Nov. 1995
 I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, ed. New York: Academic, 2007.
 M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables., New York: Dover, 1974.
 J. G. Proakis, Digital Communications, 3rd ed. New York: McGraw  Hill, 1995
 J. I. Marcum, A statistical theory of target detection by pulsed radar, IRE Trans. Inf. Theory, 1960, IT6, pp. 59267
 M. K. Simon and M.S. Alouni, Digital Communication over Fading Channels, New York: Wiley, 2005
 A. H. Heatley, A short table of the Toronto functions, Trans. Roy. Soc. (Canada), vol. 37, sec. III. 1943
 R. A. Fisher, The general sampling distribution of the multiple correlation coefficient, Proc. Roy. Soc. (London), Dec. 1928
 J. I. Marcum, A statistical theory of target detection by pulsed radar, IRE Trans. on Inf. Theory, vol. IT6, pp. 59267, April 1960
 P. Swerling, Probability of detection for fluctuating targets, IRE Trans. on Inf. Theory, vol. IT6, pp. 269  308, April 1960
 H. Sagon, Numerical calculation of the incomplete Toronto function, Proceedings of the IEEE, vol. 54, Issue 8, pp. 1095  1095, Aug. 1966
 M. M. Agrest and M. Z. Maksimov, Theory of incomplete cylindrical functions and their applications, New York: SpringerVerlag, 1971
 S. L. Dvorak, Applications for incomplete LipschitzHankel integrals in electromagnetics, IEEE Antennas Prop. Mag. vol. 36, no. 6, pp. 2632, Dec. 1994
 J. F. Paris, E. MartosNaya, U. FernandezPlazaola and J. LopezFernandez Analysis of Adaptive MIMO transmit beamforming under channel prediction errors based on incomplete LipschitzHankel integrals, IEEE Trans. Veh. Tech., vol. 58, no. 6, July 2009
 V. M. Kapinas, S. K. Mihos and G. K. Karagiannidis, On the Monotonicity of the Generalized Marcum and Nuttall QFunctions, IEEE Trans. Inf. Theory, vol. 55, no. 8, pp. 37013710, Aug. 2009
 A. R. Miller, Incomplete LipschitzHankel integrals of Bessel functions, J. Math. Anal. Appl. vol. 140, pp. 476484, 1989
 A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev, Integrals and Series, 3rd ed. New York: Gordon and Breach Science, 1992, vol. 1, Elementary Functions.
 G. N. Watson, A Treatise on the theory of Bessel functions, 2nd ed. London, UK, Cambridge Univ. Press, 1944
 P. C. Sofotasios, On Special Functions and Composite Statistical Distributions and Their Applications in Digital Communications over Fading Channels, Ph.D Dissertation, University of Leeds, UK, 2010

A. J. Goldsmith,
Wireless Communications, 1st ed. New York, Cambridge univ. Press, 2005.
Figure 2: Behaviour of the bounds in equations and for
Figure 3: Behaviour of the lower bound in for
Figure 4: Behaviour of the exact solution and the performance bounds to for and different
Figure 5: Behaviour of the exact solution and the performance bounds to for and different
Figure 6: Behaviour of the exact solution and the performance bounds to for different values of and