Log-Quadratic Bounds for the Gaussian Q-function
We present bounds of quadratic form for the logarithm of the Gaussian Q-function. We also show an analytical method for deriving log-quadratic approximations of the Q-function and give an approximation with absolute error less than $10^{-3}$.
\section{Introduction}
Approximations and bounds for the Gaussian Q-function have been studied extensively, motivated by the fact that the Q-function is not an elementary function \cite{mills26, gordon41, birnbaum42, tate53, boyd59, borj79, chiani03, deabreu09, chang11}\footnote{The Q-function is the tail probability for a standard normal random variable: \begin{eqnarray}
Q(x) \triangleq \frac{1}{\sqrt{2\pi}} \int_x^\infty e^{-\frac{u^2}{2}} \mathrm{d} u.
\end{eqnarray}}. There has been specific effort in deriving bounds that remain tractable under multiplication and exponentiation, especially those with a single exponential term. The classical example is the Chernoff bound \cite{chernoff52},
\begin{equation}
Q(x) \le \frac{1}{2}e^{-\frac{x^2}{2}}, \quad x \ge 0.
\end{equation}
In recent work, C{\^o}t{\'e} et al. \cite{cote12} showed strong Chernoff-type lower bounds for the Q-function of the form 
\begin{equation}
Q(x) \ge \alpha e^{-\beta x^2}, \quad x \in \mathbb{R}.
\end{equation}
Chernoff-type bounds are sufficient for many applications but may be loose for small values of $x$.  L{\'o}pez-Ben{\'\i}tez and Casadevall \cite{lopez11} considered Q-function approximations of the log-quadratic form
\begin{eqnarray}
Q(x) &\approx& e^{-ax^2-bx-c}, \quad x \ge 0,
\end{eqnarray}
and determined the parameters $a,b,c$ using numerical fitting techniques for different ranges of $x$. We consider bounds and approximations of this form from an analytical perspective. The bounds are stated by the following theorem.

\begin{thm}
\begin{eqnarray}
Q(x) &\ge& \frac{1}{2} e^{-\frac{1}{2}x^2 - \sqrt{\frac{2}{\pi}}x}, \quad x \ge 0, \\
Q(x) &\le& \frac{1}{2} e^{-\frac{1}{\pi}x^2 - \sqrt{\frac{2}{\pi}}x}, \quad x \ge 0.
\end{eqnarray}
\end{thm}
\noindent We derive the following approximation, which has absolute error less than $10^{-3}$.
\begin{eqnarray}
Q(x) &\approx& \frac{1}{2} e^{-0.374 x^2 - 0.777 x}, \quad x \ge 0.
\end{eqnarray}


\section{Bounds and Approximations}
We show the proof of Theorem 1 and then the derivation of the approximation. ~\\ ~\\
\textit{Proof of Theorem 1.} 
Let the standard normal density be given by
\begin{equation}
\phi(x) \triangleq \frac{e^{-\frac{x^2}{2}}}{\sqrt{2\pi}}. 
\end{equation}
Define
\begin{eqnarray}
L(x) &\triangleq& \log(Q(x)), \\
\ell(x) &\triangleq& \frac{\partial L(x)}{\partial x} = -\frac{\phi(x)}{Q(x)}, \\
\ell'(x) &\triangleq& \frac{\partial^2 L(x)}{\partial^2 x} =\frac{\phi(x)Q(x)x-\phi(x)^2}{Q(x)^2} , \\
L_l(x) &\triangleq& -\frac{1}{2}x^2-\sqrt{\frac{2}{\pi}}x-\log(2), \\
\ell_l(x) &\triangleq& \frac{\partial L_l(x)}{\partial x} = -x-\sqrt{\frac{2}{\pi}}, \\
\ell_l'(x) &\triangleq& \frac{\partial^2 L_l(x)}{\partial x^2} = -1.
\end{eqnarray}
For the lower bound, we show that $L_l(x) \le L(x)$ for $x\ge0$. It is clear that $L_l(0) = L(0)$, so it is sufficient to show that $\ell_l(x) \le \ell(x)$ for $x \ge 0$. Again, we have $\ell_l(0) = \ell(0) = -\sqrt{\frac{2}{\pi}}$, so we are left to show that $\ell_l'(x) \le \ell'(x)$, or equivalently $\ell'(x) \ge -1$, for $x \ge 0$. The desired property 
\begin{eqnarray}
\ell'(x) &=& \frac{\phi(x)Q(x)x-\phi(x)^2}{Q(x)^2} \ge -1
\end{eqnarray}
is equivalent to the statement
\begin{eqnarray}
Q(x)^2+Q(x)\phi(x)x-\phi(x)^2 \ge 0.
\end{eqnarray}
The above inequality follows from a simple bound on Mill's ratio (e.g. \cite{sampford53}).

The upper bound is found with similar reasoning. It is sufficient to show that $\ell'(x) \le -\frac{2}{\pi}$, which is equivalent to the property
\begin{eqnarray}
\frac{2}{\pi}Q(x)^2+Q(x)\phi(x)x-\phi(x)^2 \le 0.
\end{eqnarray}
An upper bound referred to in \cite{borj79} and \cite{dasgupta10} (as the Mitrinovi{\'c} inequality) is
\begin{eqnarray}
Q(x) \le \frac{2\phi(x)}{x+\sqrt{x^2+\frac{8}{\pi}}}.
\end{eqnarray}
The substitution yields
\begin{eqnarray}
\frac{2}{\pi}Q(x)^2+Q(x)\phi(x)x-\phi(x)^2 &\le& \frac{\left(\frac{8}{\pi}\right)\phi(x)^2}{\left(x+\sqrt{x^2+\frac{8}{\pi}}\right)^2} + \frac{2x\phi(x)^2}{x+\sqrt{x^2+\frac{8}{\pi}}} - \phi(x)^2 = 0.
\end{eqnarray}\qed
 ~\\
 \clearpage
\noindent\textit{Q-function approximation}.
The approximation is of the form
\begin{eqnarray}
Q(x) &\approx& e^{-ax^2-bx-c} ~\triangleq~ \hat{Q}(x).
\end{eqnarray}
for $x \ge 0$.
We use the function
\begin{eqnarray}
\hat{L}(x) &\triangleq&  \log(\hat{Q}(x)) ~=~ -ax^2-bx-c
\end{eqnarray}
and select parameters to approximately minimize the mean squared error between derivatives of $L(x)$ and $\hat{L}(x)$ for $0 \le x \le t$. Define
\begin{eqnarray}
\hat{\ell}(x,a,b) &\triangleq& \frac{\partial \hat{L}(x)}{\partial x} = -2ax-b, \\
\hat{\ell}'(x,a) &\triangleq& \frac{\partial^2 \hat{L}(x)}{\partial^2 x} = -2a.
\end{eqnarray}
The mean squared value for the difference of second derivatives of the two terms is proportional to
\begin{eqnarray}
\int_0^t (\ell'(x)-\hat{\ell}'(x,a))^2 \mathrm{d} x.
\end{eqnarray}
Differentiating with respect to $a$ and equating to zero gives an optimal value of $a$ for a given $t$:
\begin{eqnarray}
\int_0^t \ell'(x) \mathrm{d} x &=& \int_0^t \hat{\ell'}(x,a) \mathrm{d} x, \\
\int_0^t\ell'(x) \mathrm{d} x &=& \frac{\phi(0)}{Q(0)} - \frac{\phi(t)}{Q(t)} ~=~ \sqrt{\frac{2}{\pi}} - \frac{\phi(t)}{Q(t)}, \\
\int_0^t \hat{\ell'}(x,a) \mathrm{d} x &=& -2at, \\
a^*(t) &=& -\frac{1}{2t}\left (\sqrt{\frac{2}{\pi}} - \frac{\phi(t)}{Q(t)} \right ).
\end{eqnarray}
The mean squared value for the difference of first derivatives is proportional to
\begin{eqnarray}
\int_0^t (\ell(x)-\hat{\ell}(x,a,b))^2 \mathrm{d} x.
\end{eqnarray}
Differentiating with respect to $b$, equating to zero, and substituting $a = a^*(t)$ gives the optimal value of $b$ as a function of $t$,
\begin{eqnarray}
\int_0^t \ell(x) \mathrm{d} x &=& \int_0^t \hat{\ell}(x,a,b) \mathrm{d} x, \\
\int_0^t \ell(x)\mathrm{d} x &=& \log(Q(t)) - \log(Q(0)) ~=~ \log \left(2Q(t) \right), \\
\int_0^t \hat{\ell}(x,a,b)\mathrm{d} x &=& -at^2-bt, \\
b^*(t) &=& -a^*(t) t - \frac{1}{t} \log(2Q(t)), \\
b^*(t)&=& \frac{1}{2}\left (\sqrt{\frac{2}{\pi}} - \frac{\phi(t)}{Q(t)} \right ) - \frac{1}{t} \log(2Q(t)).
\end{eqnarray}
We simply use $c = \log(2)$ to match $Q(0) = \frac{1}{2}$. This gives the general approximation for a selected $t$
\begin{eqnarray}
Q(x) \approx \frac{1}{2} \exp \left ({\frac{1}{2t}\left (\sqrt{\frac{2}{\pi}} - \frac{\phi(t)}{Q(t)} \right )x^2 + \left ( \frac{1}{t} \log(2Q(t)) - \frac{1}{2}\left (\sqrt{\frac{2}{\pi}} - \frac{\phi(t)}{Q(t)} \right ) \right )x} \right),\quad x \ge 0.
\end{eqnarray}
For a given $t$, an initial evaluation $Q(t)$ is required.

By making the parameter $t$ large, the approximation becomes increasingly more accurate for large values of $x$ and less accurate for small values of $x$. Selecting $t = 1.295$ approximately minimizes the maximum absolute error $|Q(x) - \hat{Q}(x)|$ for all $x \ge 0$; the corresponding error is $9.485\times 10^{-4}$. This gives the approximation
\begin{equation}
Q(x) \approx \frac{1}{2} e^{-0.374 x^2 - 0.777 x}, \quad x \ge 0.
\end{equation}

\qed


\bibliographystyle{jota}

2013-04-10T02:02:19.000
1304
1304.2488
en
https://arxiv.org/abs/1304.2488
We present bounds of quadratic form for the logarithm of the Gaussian Q-function. We also show an analytical method for deriving log-quadratic approximations of the Q-function and give an approximation with absolute error less than $10^{-3}$.
Probability (math.PR)
Log-Quadratic Bounds for the Gaussian Q-function
Qwen/Qwen-72B
1. YES
2. YES