Diss / 10
.pdf
Estimate of Stochastic Process Frequency-Time Parameters |
577 |
of the trigger forming at the output normalized by amplitude and shape function ητ(t), the counter of spikes, the integrator, and the divisor determining the estimate of spike duration average τ*. The threshold M is given by external source.
To define the statistical characteristics of estimates τ* and θ* we assume that the condition T τcor |
||||||||
|
|
|
|
|
|
|
|
– |
is satisfied. In this case, we can approximately assume N ≈ NT. Then |
||||||||
|
|
|
|
1 |
T |
|
||
|
|
|
|
|
||||
|
|
|
|
|
|
|||
|
|
|
NT ∫ητ (t)dt, |
|||||
τ |
|
= |
|
|
||||
|
|
|
|
|
|
|
0 |
(15.315) |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
T |
|
||
θ |
= |
|
|
η (t)dt. |
||||
|
|
|
|
|||||
|
|
|
|
NT ∫ |
θ |
|||
|
|
|
|
|
|
|
0 |
|
The mathematical expectation of estimates can be presented in the following form:
|
|
|
1 |
T |
|
||
|
|
|
|
||||
|
|
|
|
|
|||
|
|
NT ∫ ητ (t) dt, |
|||||
τ |
|
= |
|
||||
|
|
|
|
|
|
0 |
(15.316) |
|
|
|
|
|
|
|
|
|
|
|
1 |
T |
|
||
θ |
= |
|
η (t) dt. |
||||
|
|
|
|
||||
|
|
|
NT ∫ |
θ |
|||
|
|
|
|
|
|
0 |
|
Taking into consideration (15.291) and (15.292), we see that τ = τ and θ = θ. In other words, the estimates of the spike duration average τ* and the average of interval θ* between the spikes are unbiased as a first approximation.
Determine the estimate variance of the average spike duration of stochastic process at level M:
Var{τ } = (τ − τ )2 = |
1 |
|
∫T ∫T ητ (t1)ητ (t2 ) dt1dt2 − τ 2. |
(15.317) |
|||
[ |
|
|
]2 |
||||
NT |
|||||||
|
|
|
|
0 |
0 |
|
|
The mathematical expectation |
|
|
|
|
|||
ητ (t1)ητ (t2 ) = τ (t1 − t2 ) |
(15.318) |
||||||
is defined by (15.298). By analogy with (15.299), we can define the estimate variance of the average spike duration τ*
Var{τ } = |
1 |
|
2 |
T |
− |
τ |
τ (t)dt − τ |
2 |
|
(15.319) |
|||
|
|
|
|
1 |
|
|
|
. |
|||||
|
|
2 |
|
|
|
||||||||
N |
|
|
|||||||||||
|
|
|
T ∫ |
|
T |
|
|
|
|
||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
The estimate variance of the average of interval θ* between the spikes can be presented in the following form:
Var{θ } = |
1 |
|
2 |
T |
|
τ |
|
|
2 |
|
|
|||
|
− |
θ (t)dt − θ |
|
|
||||||||||
|
|
|
|
1 |
|
|
|
. |
(15.320) |
|||||
|
|
2 |
|
|
|
|||||||||
N |
|
|
||||||||||||
|
|
|
T ∫ |
|
T |
|
|
|
|
|
||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
580 Signal Processing in Radar Systems
As we discussed before, the estimates of variances are unbiased and the variances of variance estimates are defined by (13.62) where we need to use the corresponding correlation function of the investigated stochastic process and its first derivative instead σ2 (τ).
Going forward, we assume that the condition T τcor is satisfied. In this case, the errors Varx
and |
Varx will be small compared to Varx and Varx. To define the bias of the mean-square frequency |
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
· |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
· |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
under previous assumptions, we can think that |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
Varx + |
|
|
|
|
|
|
|
|
|
|
1+ |
Varx |
|
|
|
1+ 1 |
|
|
|
|
|
|
|
Varx2 |
1 − 1 |
|
|
|
Varx2 |
. |
|||||||||||||||
|
|
= |
1 |
|
|
|
Varx |
= |
|
Varx |
|
≈ |
|
|
Varx |
− 1 |
|
|
|
Varx |
+ 3 |
||||||||||||||||||||||||||||||
|
f |
f |
|
|
|
|
|
f |
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|||||||||||||||||||||||||
|
|
|
|
2π |
|
Varx + |
|
|
Varx |
|
|
|
|
|
|
|
Varx |
|
|
|
|
|
|
2 |
|
Varx |
8 |
|
|
|
|
2 |
|
2 |
Varx |
8 |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
1+ |
|
|
|
|
|
|
|
|
Varx |
|
Varx |
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(15.330) |
|
Now, we are able to obtain the relative bias of estimate |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
f |
|
f |
− f |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
= |
|
≈ − |
1 |
|
|
Varx |
|
− 3 |
Varx |
|
+ 2 |
|
Varx |
Varx |
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
(15.331) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8 |
|
2 |
|
2 |
|
|
|
Varx Varx |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
f |
|
|
|
|
f |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
|
|
|
Varx |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
where under the condition T τcor we have |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
T |
dR(τ) |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
Varx = |
|
|
|
∫ |
|
|
dτ. |
|
|
|
|
|
|
(15.332) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
0 |
|
dτ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
As a result, the relative bias can be presented in the following form:
|
|
|
|
|
|
|
∞ |
|
′′(τ) |
|
2 |
∞ |
|
∞ |
′2 |
(τ) |
|
|
|
|
f |
|
1 |
|
|
|
2 |
|
|
|
|
||||||||||
|
|
|
|
|
= |
|
−∫ |
|
|
|
dτ + 3∫ |
|
(τ)dτ + 2∫ |
|
|
dτ |
, |
(15.333) |
||
|
|
|
|
|
|
′′(0) |
|
|
′′(0) |
|||||||||||
|
f |
|
2T |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
0 |
|
|
|
|
|
where (τ), ′(τ), and ″(τ) are the normalized correlation function of the investigated stochastic process and its first and second derivatives. As applied to the exponential normalized correlation function of the investigated stochastic process given by (15.305), we can write
|
|
|
|
5 2π |
|
||
f |
= |
(15.334) |
|||||
|
|
|
|
|
32αT , |
||
|
|
|
|
|
|||
|
f |
|
|
|
|||
that is, it means that the bias of estimate is inversely proportional to the observation time interval T. Define a dispersion of the mean-square frequency estimate
|
|
|
|
|
|
|
|
|
|
+ |
Varx |
1 + |
Varx |
|
|
|
||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
|
Varx |
|
|
|
||
D{f } = ( f − f )2 = f 2 |
|
+1 − 2 |
|
|
(15.335) |
|||||||||||||||
|
|
|
|
|
|
. |
||||||||||||||
|
|
|
|
|
|
|
|
|
1 |
+ |
Varx |
|
1 + |
Varx |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
Varx |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Estimate of Stochastic Process Frequency-Time Parameters |
581 |
Using two-dimensional Taylor expansion in series for the first and third terms in (15.335) about the points
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx = 0 |
|
and |
|
|
Varx |
= 0 |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
|
|
|
|
Varx |
|
|
|
|
|
|
|
|
|
|
||||||||||
and limiting by terms of the second order, we obtain |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
1 |
|
|
|
|
1 |
2 |
|
|||||||
|
D{f |
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
|
|
|
Varx |
|
||||||||||||||||||||||||||||
|
} |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
+ 1 − |
|
|
Varx |
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
≈ |
1 |
+ |
|
|
|
|
|
1 − |
|
|
|
|
|
|
|
|
+ |
|
|
|
2 1 |
+ |
|
|
|
|
− |
|
|
|
|||||||||||||
|
|
|
|
2 |
|
|
|
|
Varx |
|
|
Varx |
|
|
|
2 |
|
2 |
Varx |
|
8 |
2 |
||||||||||||||||||||||||||||
|
|
f |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
|
|
|
|
|
Varx |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
1 − 1 |
Varx |
+ |
3 |
|
|
|
Varx2 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
(15.336) |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
Varx |
|
8 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Varx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Under averaging to a first approximation, we obtain |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
} |
|
|
∞ |
|
′′(τ) |
|
2 |
∞ |
|
|
|
|
|
∞ |
′2 (τ) |
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
D{ |
f |
|
1 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
≈ |
|
∫ |
|
|
|
|
|
|
|
|
|
dτ + ∫ |
|
(τ)dτ + 2∫ |
|
|
|
dτ . |
|
(15.337) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
′′(0) |
|
′′(0) |
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
f |
|
|
|
T |
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
||||
As applied to the exponential normalized correlation function given by (15.305), we obtain |
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D{ |
|
|
} 3 2π 0.47 |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
|
|
|
|
|
|
|
|
|
(15.338) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
≈ 16αT |
≈ αT . |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Comparing (15.338) with the relative variance of estimate of the average number of Gaussian stochastic process spikes for the same normalized correlation function given by (15.308), we see that according to (15.326), the average mean-square frequency estimate or the average number of spikes leads to the bias of estimate and decrease in estimate dispersion approximately in 2.8 times.
15.8 SUMMARY AND DISCUSSION
Methods of the correlation function estimate can be classified on three groups subject to a principle of realization of delay and other elements of correlators: analog, digital, and analog-to- digital. In turn, the analog measurement procedures can be divided based on the methods using a representation of the investigated stochastic process both as the continuous process and as the sampled process. As a rule, physical delays are used by analog methods with continuous representation of the investigated stochastic process. Under discretization of investigated stochastic process in time, the physical delay line can be changed by corresponding circuits. Under the use of digital procedures to measure the correlation function estimate, the stochastic process is sampled in time and transformed into binary number by analog-to-digital conversion. Further operations associated with the signal delay, multiplication, and integration are carried out by the shift registers, summator, and so on.
We can see that the maximum value of variance of the correlation function estimate corresponds to the case τ = 0 and is equal to the variance of the stochastic process variance estimate given by (13.61) and (13.62). Minimum value of variance of the correlation function estimate corresponds to the case when τ τcor and is equal to one-half of the variance of the stochastic process variance estimate.
582 |
Signal Processing in Radar Systems |
The correlation function of stationary stochastic process can be presented in the form of expansion in series with respect to earlier-given normalized orthogonal functions. The variance of correlation function estimate increases with an increase in the number of terms of expansion in series v. Because of this, we must take into consideration this fact choosing the number of terms under expansion in series.
In some practical problems, the correlation function of stochastic process can be measured accurately with some parameters defining a character of its behavior. In this case, the measurement of correlation function can be reduced to measurement or estimation of unknown parameters of correlation function. Because of this, we consider the optimal estimate of correlation function arbitrary parameter assuming that the investigated stochastic process ξ(t) is the stationary Gaussian stochastic process observed within the limits of time interval [0, T] in the background of Gaussian stationary noise ζ(t) with known correlation function.
The optimal estimate of stochastic process correlation function can be found in the form of estimations of the elements Rij of the correlation matrix R or elements Cij of the inverse matrix C. In the case of Gaussian stationary stochastic process with the multidimensional pdf given by (12.169), the solution of likelihood ratio equation allows us to obtain the estimates Cij and, consequently, the estimates of elements Rij of the correlation matrix R.
Based on (15.272), we can design the flowchart of spectral density measurer shown in Figure 15.13. The spectral density value at the fixed frequency coincides accurately within the constant factor with the variance of stochastic process at the filter output with known bandwidth. Operation principles of the spectral density measurer are evident from Figure 15.13. To define the spectral density for all possible values of frequencies, we need to design the multichannel spectrum analyzer and the central frequency of narrowband filter must be changed discretely or continuously. As a rule, we need to carry out a shift by the spectrum frequency of the investigated stochastic process using, for example, the linear law of frequency transformation instead of filter tuning by frequency. The structure of such measurer is depicted in Figure 15.14. The sawtooth generator controls the operation of measurer changing a frequency of heterodyne.
In many applications, we need to know the statistical parameters of stochastic process spike (see Figure 15.15a): the spike mean or the average number of down-up cross sections of some horizontal level M within the limits of the observation time interval [0, T], the average duration of spike, and the average interval between spikes. To measure these parameters of spikes, the stochastic process realization x(t) is transformed by the nonlinear transformer (threshold circuitry) into the pulse sequence normalized by the amplitude ητ with duration τi (Figure 15.15b) or normalized by the amplitude ηθ with duration θi (Figure– 15.15c).
The mean-square frequency f given by (15.6) is widely used as a parameter characterizing the spectral density of the stochastic process. The value f– defines a dispersion of component of the stochastic process spectral density relative to zero frequency. As applied to low-frequency stationary stochastic processes, the mean-square frequency f– characterizes the effective bandwidth of spectral density.
REFERENCES
1.Lunge, F. 1963. Correlation Electronics. Leningrad, Russia: Nauka.
2.Ball, G.A. 1968. Instrumental Correlation Analysis of Stochastic Processes. Moscow, Russia: Energy.
3.Kay, S.M. 1993. Fundamentals of Statistical Signal Processing: Estimation Theory. Upper Saddle River, NJ: Prentice Hall, Inc.
4.Lampard, D.G. 1955. New method of determining correlation functions of stationary time series.
Proceedings of the IEE, C-102(1): 343.
5.Kay, S.M. 1998. Fundamentals of Statistical Signal Processing: Detection Theory. Upper Saddle River, NJ: Prentice Hall, Inc.
6.Gribanov, Yu.I. and V.L. Malkov. 1978. Selected Estimates of Spectral Characteristics of Stationary Stochastic Processes. Moscow, Russia: Energy.


Var{