This paper presents new theoretical results comparing the distributions of Gaussian quadratic forms, advancing understanding in probability theory and statistical analysis.
Contribution
It introduces two theorems that provide novel insights into the distributional properties of Gaussian quadratic forms.
Findings
01
Derived new comparison theorems for Gaussian quadratic forms
02
Enhanced understanding of distributional behavior of quadratic forms
03
Potential applications in statistical hypothesis testing
Abstract
New results on comparison of distributions of Gaussian quadratic forms are presented
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsAnalytic Number Theory Research
Full text
Problems of Information Transmission,
v. 53, no. 3, pp. 3–15, 2017.
M. V. Burnashev
TWO THEOREMS ON DISTRIBUTION OF GAUSSIAN QUADRATIC
FORMS
111The research was carried out at the IITP RAS at the
expense of the Russian Foundation for Sciences
(project 14-50-00150).
New results on comparison of
distributions of Gaussian quadratic forms are presented.
Let ξ1,…,ξn – independent
N(0,1)-Gaussian random variables. For
a=(a1,…,an)∈R+n and
x≥0 consider the following probability
[TABLE]
We are interested in what
a,b∈R+n and x
the following inequality holds
[TABLE]
Below the vector a∈R+n
is called monotone, if
a1≥a2≥…≥an≥0.
As usually, a≥b means
ai≥bi, i=1,…,n.
1. Known comparison theorem.
In [1, Theorem 1] the following result was proved.
Let a,b∈R+n -
monotone vectors and the following condition for them is fulfilled
[TABLE]
Then for any x≥2i=1∑nbi
the inequality (1) holds.
Remark 1. In order (2) to be valid, we need, in
particular, imaxai≥imaxbi,
what is rather restrictive.
Inequality (1) is useful in problems of detection of
stochastic signals in Gaussian noise. But application of Theorem 1
from [1] in such problems is rather difficult because of
the restrictive assumption (2) and the requirement
x≥2i=1∑nbi (in stochastic signals
detection problems the inequality (1) is usually
required for x<i=1∑nbi).
In the paper the inequality (1) is proved under different
from (2) assumption. First, a simple similar
Proposition 1 is proved, and then it is strengthened using
additional arguments (Theorems 1 and 2).
Concerning applications of the inequality (1) it should
be mentioned that such results are helpful in problems of detection
of Gaussian stochastic signals in the background of independent
additive Gaussian noise [2–5].
Consider, for example, the problem of detection of Gaussian
stochastic signal vector s in the background of
independent additive Gaussian noise ξ. If the
vector σ of the vector s
intensities is known, then the logarithm of the corresponding
likelihood ratio in that problem is a Gaussian quadratic form,
similar to one considered above. Assume that we know only that
the vector σ belongs to the given set
E. Then natural question arises: is it possible to
replace the set E by a smaller set E0
without loss of detection quality (in particular, to replace
E by a single point σ0) ?
Such problem will be considered by author in the paper
[6].
Some results showing validity of the inequality (1)
can also be found in [7, 1].
2. The first result. For vectors
a,b∈R+n
introduce functions
[TABLE]
P r o p o s i t i o n 1. *Let
a,b∈R+n and
D(a,b)>1.
If x satisfies the condition *
[TABLE]
thenβ(x,a)≤β(x,b).
P r o o f. Consider first the case c=a,
d=b. We have
[TABLE]
where the ellipsoid
A(x,a)∈Rn has axes
{x/ai,i=1,…,n} and volume
V(A(x,a)), proportional to
(i=1∏nai)−1/2. In order to
compare probabilities β(x,a) and
β(x,b), consider the difference
(see (5))
[TABLE]
Changing variables
zi=ai/biyi, i=1,…,n, in the second
integral we have (D=D(a,b))
[TABLE]
Therefore Δ(x,a,b)≤0, if
1−De−xd/2≤0, i.e. if
[TABLE]
where W(a)=i=1∏nai,
W(b)=i=1∏nbi, from which
Proposition 1 for c=a,
d=b follows. If
c≤a and
d≥b, then
β(x,a)≤β(x,c) and
β(x,d)≤β(x,b) for any
x (see the statement 1 of Lemma below in section 3). Moreover,
β(x,c)≤β(x,d), if
x satisfies the condition (4). ▲
Remark 2. Auxiliary vectors c and
d in (4) allow sometimes to relax
the constraint (4) in comparison with
c=a and
b=d. One of further extensions
(Theorem 2) concerns the choice of vectors c and
d for given a and b.
Below Proposition 1 is strengthened increasing the right-hand
side of the condition (4) based on various
additional arguments.
3. Strengthening 1. We use the following auxiliary result.
L e m m a 1. 1) Let
a,b∈R+n and for
components with indices i,j we have
max{ai,aj}≥max{bi,bj}
and aiaj≥bibj. Let also ak≥bk
for all k=i, k=j. Then
β(x,a)≤β(x,b)
for any x.
If b≤a, then
β(x,a)≤β(x,b) for any
x.
Assume that it is possible to partition the set of indices
I={1,2,…,n} on k≥1 parts I1,…,Ik,
such that
I=j=1⋃kIj, Ii∩Ij=∅,
i=j, and the following conditions are fulfilled
[TABLE]
where
[TABLE]
Then β(x,a)≤β(x,b)
for any x.
P r o o f. 1) Assume first that n=2. Then the integration
region A(x,b)∈R2 for
β(x,b) has volume
V(A(x,b)), proportional to
(b1b2)−1/2, and V(A(x,b))≥V(A(x,a)). The random vector
ξ=(ξ1,ξ2) has the distribution density
p(y), proportional to
e−r2(y)/2,
r2(y)=y12+y22,
monotonically decreasing in r≥0. It can be checked that
[TABLE]
Therefore for given volume V (i.e. for given product
b1b2=T) the value β(x,b) attains its
maximum when b1=b2=T, and monotonically decreases
when b1 deviates from T, from which the inequality
(1) follows for any x.
If n>2, then the inequality (1) holds for any fixed
{ξk}, k=i, k=j and any x, from which
necessary assertion follows.
If b≤a, then
A(x,a)⊆A(x,b), and therefore
β(x,a)≤β(x,b) for any x.
That assertion follows from part 1). It is sufficient to
consider the case k=2, I1={1,…,m} and
I2={m+1,…,n} for some 1<m<n. Introduce
n-vector a0=(a0,1,…,a0,1,a0,2,…,a0,2), consisting of
m components a0,1 and n−m components a0,2. Then
repeatedly applying lemma’s part 1) it is possible to show that
β(x,a)≤β(x,a0) for any
x. Since b≤a0, then
β(x,a0)≤β(x,b) for any
x. ▲
We strengthen the Proposition 1. Setting for convenience
c=a and
d=b, assume that there exists
i≤n−1, such that f(i+1,a,b)<f(i,a,b). Then ai+1<ai
(since bi+1≤bi) and bi+1<bi. We build the
vector a(1) such that
β(x,a)≤β(x,a(1))
for all x and
β(x,a(1))≤β(x,b)
for x, satisfying, generally, a weaker than (4)
constraint. For that purpose we use the following procedure.
Decrease ai down to the value ai(1) and increase
respectively ai+1 up to the value ai+1(1), such that
the following three conditions are satisfied:
ai(1)≥ai+1(1);
aiai+1=ai(1)ai+1(1);
f(i+1,a(1),b)=f(i,a(1),b), where
a(1) - obtained that way from a
the new monotone vector (it differs from
a only in components with indices i and i+1).
Then we have ai+1(1)<ai(1),
f(i,a(1),b)<f(i,a,b) and
f(i+1,a,b)<f(i+1,a(1),b). We also have
β(x,a)≤β(x,a(1))
for any x (due to Lemma 1). Then after standard calculations
we have
[TABLE]
where
[TABLE]
Note that values f(i,a,b),f(i,a(1),b) may be negative.
Similarly, instead of the vector a we may change
the vector b (but in opposite direction), replacing
it by the vector b(1) such that
β(x,b)≤β(x,b(1))
for all x and
β(x,a)≤β(x,b(1))
for x, satisfying, generally, a weaker than (4)
constraint. For that purpose, increase bi up to the value
bi(1) and decrease respectively bi+1 down to the value
bi+1(1), such that the following three conditions
are satisfied:
bi(1)≥bi+1(1);
bibi+1=bi(1)bi+1(1);
f(i+1,a,b(1))=f(i,a,b(1)), where
b(1) - obtained that way from b
the new monotone vector (it differs from
b only in components with indices i and i+1).
Then we have bi+1(1)<bi(1),
f(i,a,b(1))<f(i,a,b) and
f(i+1,a,b)<f(i+1,a,b(1)). We also have
β(x,b(1))≤β(x,b)
for any x (due to Lemma). Then after standard calculations
we have
Compare values f(i,a(1),b) and
f(i,a,b(1)). It is possible to check
that f(i,a(1),b)<f(i,a,b(1)) (i.e. changing the vector
a gives better result), if
bibi+1>aiai+1. If
bibi+1<aiai+1, then
f(i,a(1),b)>f(i,a,b(1)), i.e. changing the vector
b gives better result.
After that we apply the procedure described to obtained vectors
a(1) or b(1) and so on.
Unfortunately, the author was not able to investigate the optimal
sequence of changes the vectors a or
b. For that reason we limit ourselves to the case
when only the vector a (or only the vector
b) is changing.
Note that if zi≥0 (i.e.
bibi+1≥aiai+1), then ai≤bi and
f(i,a,b)≤0, or
ai+1≤bi+1 and
f(i+1,a,b)≤0.
Denote by i0=i0(a,b) the
minimal of indices i, such that
f(i0,a,b)=imaxf(i,a,b), and by
i1=i1(a,b)≥i0
the minimal of indices i, such that
f(i1,a,b)=imaxf(i,a,b). If
i1(a,b)=n, then the method used
does not improve the condition (4) and then in
Proposition 1 we have
imaxf(i,a,b)=1/bn−1/an. Therefore we assume that
i1(a,b)≤n−1.
Change vector a.
First, for fixed monotone b we change the vector
a. Choose arbitrary
i≥i0(a,b), such that
f(i+1,a,b)<f(i,a,b), and apply to a
the procedure described. Then for obtained in that way from
a the new monotone vector a(1)
(it differs from a only in components with indices
i and i+1) we have ai+1(1)<ai(1),
f(i,a(1),b)<f(i,a,b) and
f(i+1,a,b)<f(i+1,a(1),b). We also have
β(x,a)≤β(x,a(1))
for any x (due to Lemma).
As an initial index i we may set, for example,
i=i1(a,b). Then apply the
procedure described to the received vector a(1)
(i.e. find a new index i1(a(1),b),
corresponding component ai1(1) and transform components
ai1(1) and ai1+1(1), such that three
conditions above are satisfied). It will give a new monotone vector
a(2). Then apply that procedure to the vector
a(2) and so on.
As a result, we may get a sequence of monotone vectors
am, m=1,2,…, converging to the
monotone vector a0. Let
k1, 1≤k1≤i0(a,b) -
minimal of indices i, which were used on all stages of getting
the sequence {am}. Then for the limiting
monotone vector a0 the following conditions
are satisfied:
[TABLE]
i.e. the function f(i,a0,b) is
constant for k1≤i≤n. Components
{ai0,i=1,…,k1−1} of the vector
a0 coincide with corresponding components
{ai,i=1,…,k1−1} of the initial vector
a (they will not participate in getting the vector
a0).
We find the value f(k1,a0,b).
Since
[TABLE]
then due to the last of conditions (7) values
{f(i,a0,b)} satisfy also
equations
[TABLE]
or, equivalently
[TABLE]
Define the value T(k,a,b) as the
unique root of the equation
[TABLE]
Then (see (8) and (7))
T(k1,a,b)=T(k1,a0,b)=f(k1,a0,b). Note that if
i=k1∑nln(ai/bi)>0,
then T(k1,a,b)>0, and if
i=k1∑nln(ai/bi)<0, then
T(k1,a,b)<0.
The index k1(a,b),
1≤k1≤i0(a,b) is
defined as follows
[TABLE]
Consider now first k1−1 coordinates ai0,
i=1,…,k1−1 of the vector a0. Again
we find for them corresponding values i0′,i1′,k1′
and replace coordinates ai0, i=k1′,…,k1−1 by
corresponding coordinates ai0′, i=k1′,…,k1−1,
such that the function f(i,a0,b)
is constant for k1′≤i≤k1−1. Continuing that
process, we get the monotone vector a(0),
such that
β(x,a)≤β(x,a(0))
for any x. Moreover, the function
f(i,a(0),b) is piecewise constant
and does not decrease in i.
Change vector b. Similarly for fixed
a we may sequentially change the vector
b (but in opposite direction), again using the index
i1(a,b) and getting the sequence of
monotone vectors {b(i)}, i=1,2,…,
such that
β(x,b)≥β(x,b(1))≥β(x,b(2))≥… for any x (because of
Lemma). As in the case of vector a, assume that
i1(a,b)≤n−1. Then
f(i1+1,a,b)<f(i1,a,b) and
bi1+1<bi1. Apply to b the procedure
described and get new monotone vector b(1)
(it differs from b only in components with indices
i and i+1). We have for it
bi1+1(1)<bi1(1),
f(i1,a,b(1))<f(i1,a,b) and
f(i1+1,a,b)<f(i1+1,a,b(1)). We also have
β(x,b(1))≤β(x,b)
for any x (because of Lemma).
Then apply the procedure described to the obtained vector
b(1) (i.e. find the new index
i1(a,b(1)) and corresponding
component bi1(1) and transform components
bi1(1) and bi1+1(1), such that three
conditions above are satisfied). It will give the vector
b(2). Then apply that procedure to the vector
b(2) and so on.
As a result, we may get a sequence of monotone vectors
bm, m=1,2,…, converging to the
monotone vector b0. Let
k2, 1≤k2≤i1(a,b) -
minimal of indices i1, used on all stages of getting
the sequence {bm}. Then for the limiting
monotone vector b0 the following conditions,
similar to (7), are satisfied:
[TABLE]
i.e. the function f(i,a,b(0)) is
constant for k2≤i≤n. Components
{bi(0),i=1,…,k2−1} of the vector
b(0) coincide with corresponding components
{bi,i=1,…,k2−1} of the initial vector
b (they will not participate in getting the vector
b0). Now if
β(x,a)≤β(x,b(0)),
then β(x,a)≤β(x,b).
Since
[TABLE]
then due to the last of conditions (11) values
{f(i,a,b(0))} satisfy also
equations
[TABLE]
or, equivalently
[TABLE]
Define the value D(k,a,b) as the
unique root of the equation
[TABLE]
Then (see (8) and (7))
D(k2,a,b)=D(k2,a,b(0))=f(k2,a,b(0)).
Note that if
i=k2∑nln(ai/bi)>0,
then D(k2,a,b)>0, and if
i=k2∑nln(ai/bi)<0, then
D(k2,a,b)<0.
The index k2(a,b),
1≤k2≤i0(a,b) can be
defined as follows
[TABLE]
In order to formulate the result obtained note that similarly
to Proposition 1 we may additionally introduce arbitrary vectors
c≤a and
d≥b, such that
D(c,d)>1. Then we have
T h e o r e m 1. Let
a,b∈R+n - monotone
decreasing vectors. If x satisfies condition
[TABLE]
(values T(k,c,d)andD(k,c,d) are defined in
(9) and (12), and values
k1=k1(c,d) and
k2=k2(c,d) in (10) and
(13))
thenβ(x,a)≤β(x,b).
Remark 3. If
[TABLE]
then k1(a,b)=k2(a,b)=1 and
[TABLE]
The value T(k)=T(k,a,b) can be
upper bounded, for example, as follows. Using Jensen’s inequality
Elnξ≤lnEξ,
from (9) we have
[TABLE]
Applying the inequality ln(1+x)≤x, we get
[TABLE]
For the value D(k,a,b) those estimates
work in opposite direction:
[TABLE]
C o r o l l a r y 1. If a1≥b1,
the inequality holds
[TABLE]
and any of the following conditions is valid
[TABLE]
[TABLE]
then
β(x,a)≤β(x,b) for
any x.
P r o o f. We set c=a and
b=d. Show that if a1≥b1,
and conditions (16) and (17) are satisfied,
then k1(a,b)=1. Indeed, due to
(10) for that purpose it is sufficient to have
[TABLE]
where T=T(1,a,b)≥0
(see (9)) – the unique root of the equation
[TABLE]
Since F(0)≥0 and F′(T)<0, then the equation
(19) has the unique root
T=T(1,a,b)≥0. Moreover,
if F[f(1,a,b)]≤0, then
f(1,a,b)≥T(1,a,b) (and then
k1(a,b)=1), which coincides with
the condition (17). We also have
[TABLE]
Hence b≤a(0), and therefore
β(x,a)≤β(x,a(0))≤β(x,b) for any x.
Similarly we can show that if a1≥b1, and also
conditions (16) and (18) are satisfied, then
k2(a,b)=1, and therefore
β(x,a)≤β(x,b(0))≤β(x,b) for any x. ▲
The following result is also valid.
P r o p o s i t i o n 2. Let
a,b∈R+n. If
β(x,a)≤β(x,b)
for any x>0, then
imaxai≥imaxbi.
P r o o f. It is sufficient to consider the case
imaxai=a1, imaxbi=b1,
where b1>0. Assume that a1<b1. Show that then the
inequality β(x,a)≤β(x,b)
can not hold for large x. Denoting
a1=(a1,…,a1) and
b1=(b1,0,…,0), note that
β(x,a1)≤β(x,a) and
β(x,b)≤β(x,b1).
Therefore it is sufficient to show that for large x
the inequality
β(x,a1)≤β(x,b1)
does not hold. We may limit ourselves to the case
b1=1, a1<1. Then the inequality,
equivalent to (1), takes the form
[TABLE]
The left side of that formula can be bounded using exponential
Chebychev inequality
[TABLE]
and therefore
[TABLE]
On the other hand, using the standard estimate
[TABLE]
we have
[TABLE]
from which it follows that the inequality (21) can not
hold for large x. ▲
C o r o l l a r y 2. Let
a,b∈R+n -
monotone vectors, such that
[TABLE]
*Then *
[TABLE]
4. Strengthening 2. For vectors
a,b∈R+n we use
the function f(i,a,b) from
(3). In each of vectors
a,b we change in some way
numeration of their components {i}, such that the function
f(i,a,b) become monotone increasing in
i and, in particular,
[TABLE]
Of course, then the values
β(x,a),β(x,b) will not be
changed. It is desirable to have the value
bn−1−an−1 as small as possible.
Then for c=a,
d=b the condition (4)
takes the form
[TABLE]
Note that if in (22)
i=1∑nln(ai/bi)>0, then d>0. Denote
[TABLE]
Then f(i,a,b)=f(t1,a,b)=bt1−1−at1−1, t1≤i≤n.
We try to increase the right-hand side of (22),
increasing bi, t1≤i≤n, but not changing
{ai}. Then the value d will decrease. For some
ε≥0 we set
[TABLE]
Then
[TABLE]
and (since lnz≥1−1/z)
[TABLE]
Therefore if
i=1∑t1−1lnbiai>0, then
decrease all bi−1, t1≤i≤n, on the value
ε, such that for the new vector
b′≥b (and then
β(x,b′)≤β(x,b) for any
x) we will have
[TABLE]
i.e. the function f(i,a,b′) becomes
constant for t1−1≤i≤n. For that purpose we set
[TABLE]
and then get
[TABLE]
Similarly we repeat that procedure for the obtained vector
b′. Denote
[TABLE]
Next, again if
i=1∑t2−1lnbiai>0, then
decrease all (bi′)−1, t2≤i≤n on the value
ε1, such that for the new vector
b′′≥b′ we have
[TABLE]
i.e. the function f(i,a,b′′) become
constant for t2−1≤i≤n. For that purpose we set
[TABLE]
Then we get
[TABLE]
Repeating that process, we get the following result. Denote
T h e o r e m 2. Let
a,b∈R+n and let
the function f(i,a,b) monotonically
increases in i.
If n1=n1(a,b)≥2
and x satisfies the condition
[TABLE]
where
[TABLE]
thenβ(x,a)≤β(x,b).
If n1(a,b)=1, then
β(x,a)≤β(x,b) for allx.
Explain only the statement 2 of that Theorem. If
n1(a,b)=1, then a1≥b1.
Since the function f(i,a,b)
monotonically increases in i, then ai≥bi for all
i, and therefore a≥b. Then from
the statement 2 of Lemma we get
β(x,a)≤β(x,b) for all x.
5. Examples. 1. Let b=1. If
i=1∏nai≥1, then for any x we have
(see the statement 3 of Lemma)
β(x,a)≤β(x,1), i.e.
[TABLE]
Compare results which give in that example [1, Theorem 1],
Proposition 1 and Theorems 1, 2 of the paper.
If additionally the condition (2) is also satisfied then
according to [1, Theorem 1] the inequality (27)
holds for x≥2n.
According to Proposition 1 (i.e. the condition (4)) for
c=a,
d=1, the inequality (27)
holds for x, satisfying the condition
[TABLE]
In order to apply Theorem 1 (i.e. the condition (14))
note that if a - monotone decreasing vector, then
the function f(i,a,1)=1−1/ai
also monotonically decreases in i, and therefore
k1(a,1)=1. Using formulas
(9) and (10) we get
[TABLE]
Therefore if D>1, then the condition (14)
has the form
[TABLE]
where
[TABLE]
The function f(z) monotonically increases from f(1)=1 up to
f(∞)=∞. In particular, if D≥1, then the
inequality (29) holds for x≤n. The condition
(29) may be much better than the condition
(28).
Theorem 2 (i.e. the condition (25)) may give result
better than (29), and worse than it as well
(depending on the vector a).
Next example gives the less obvious inequality, opposite to
(27).
Let a=1,
b1>1, b1>b2≥…≥bn and
b1b2<1. Theorem 1 from [1] is not applicable
here, since the condition (2) is not fulfilled.
The function f(i,1,b)=1/bi−1
monotonically increases in i and
imaxf(i,1,b)=1/bn−1.
Therefore Theorem 1 does not improve the estimate (4).
If i=1∑nlnbi<0, then from (4)
we get
[TABLE]
In order to apply Theorem 2 (i.e. the condition (25))
notice that n1=n1(a,b)=2,
d=1/b2−1. Therefore we get
The author thanks the reviewer for constructive critical remarks,
which helped to improve the paper material presentation.
REFERENCES
Список литературы
[1]
Bakirov N. K. Comparison Theorems for Distribution
Functions of Quadratic Forms of Gaussian Vectors //
Theory of Probability and Its Applications. 1995. V. 40. no. 2.
P. 404–412.
[2]
Wald A. Statistical Decision Functions. New York:
Wiley, 1950.
[3]
Lehmann E. L. Testing of Statistical Hypotheses.
New York: Wiley, 1959.
[4]
Burnashev M. V.
On the Minimax Detection of an Inaccurately Known Signal in a
White Gaussian Noise Background //
Theory of Probability and Its Applications. 1979. V. 24. no. 1.
P. 106–118.
[5]
Zhang W., Poor H.V. On Minimax Robust Detection of Stationary
Gaussian Signals in White Gaussian Noise // IEEE Trans. Inform.
Theory. 2011. V. 57. № 6. P. 3915–3924.
[6]
Burnashev M. V. On Detection of Gaussian Stochastic
Sequences // Problems of Information Transmission. 2017 (in print).
[7]
Ponomarenko L. S. On Estimating Distributions of
Normalized Quadratic Forms of Normally Distributed Random
Variables // Theory of Probability and Its
Applications. 1985. V. 30. no. 3. P. 545–549.
*Burnashev Marat Valievich
Kharkevich Institute for Information Transmission Problems,
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] Bakirov N. K. Comparison Theorems for Distribution Functions of Quadratic Forms of Gaussian Vectors // Theory of Probability and Its Applications. 1995. V. 40. no. 2. P. 404–412.
2[2] Wald A. Statistical Decision Functions. New York: Wiley, 1950.
3[3] Lehmann E. L. Testing of Statistical Hypotheses. New York: Wiley, 1959.
4[4] Burnashev M. V. On the Minimax Detection of an Inaccurately Known Signal in a White Gaussian Noise Background // Theory of Probability and Its Applications. 1979. V. 24. no. 1. P. 106–118.
5[5] Zhang W., Poor H.V. On Minimax Robust Detection of Stationary Gaussian Signals in White Gaussian Noise // IEEE Trans. Inform. Theory. 2011. V. 57. № 6. P. 3915–3924.
6[6] Burnashev M. V. On Detection of Gaussian Stochastic Sequences // Problems of Information Transmission. 2017 (in print).
7[7] Ponomarenko L. S. On Estimating Distributions of Normalized Quadratic Forms of Normally Distributed Random Variables // Theory of Probability and Its Applications. 1985. V. 30. no. 3. P. 545–549.