Diss / 10
.pdf
Estimation of Stochastic Process Variance |
475 |
We can simplify (13.186) applied to zero measurement procedure. In doing so, owing to identity of two channels of amplifier we can think that Rs(τ) = Rs0(τ). Taking into consideration the condition that Varβ 1 and neglecting the terms with double frequency 2ω0, we can write
Var {Vars*} = (1 |
+ q) |
2 4Vars2 |
T |
1 |
− |
τ |
[r |
2 |
(τ) + 4rβ (τ)Varβ + 4r |
2 |
(τ)rβ (τ)Varβ ]dτ. |
(13.187) |
||
|
|
|
|
|
|
|
||||||||
|
T |
|
|
|
||||||||||
|
|
|
∫ |
|
|
T |
|
|
|
|
|
|
||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
By analogy with the compensation method, (13.187) allows us to obtain the time interval of observation corresponding to the sensitivity threshold. To satisfy this condition, the following equality needs to be satisfied:
Var {Vars*} |
= 1. |
(13.188) |
|
Vars2 |
|||
|
|
Comparing (13.187) with (13.162), we can see that under the considered procedure the variance of the variance estimate is twice as high as the variance of the variance estimate at the compensation method. This increase in the variance of the variance estimate is explained by the presence of the second channel, which results in an increase in the total variance of the output signal. We need to note that although there is an increase in the variance of the variance estimate for the considered procedure, the present method is a better choice compared to the compensation method if the variance of intrinsic noise of amplifier varies after compensation procedure.
Let the normalized correlation functions, as before, be described by (13.165). Then (13.187) based on the following conditions αT 1 and α γ can be presented in the following form:
2 |
|
|
|
γ T − 1+ exp{−γT} |
|
|
|
Var {Vars*} = (1+ q)2 2Vars |
1 |
+ 8 |
α Varβ |
. |
(13.189) |
||
γT |
|||||||
αT |
|
|
γ |
|
|
In other words, the variance of the variance estimate in the case of the considered procedure of measurement is twice as high as the variance of the variance estimate given by (13.168) obtained using compensation method. This is true also when the random variations of the amplification coefficient are absent. In doing so, the time interval of observation corresponding to the sensitivity threshold increases twofold compared to the time interval of observation of the investigated stochastic process while using the compensation method to measure the variance.
Using the deterministic harmonic signal
s0 (t) = A0 cos(ω0t + ϕ0 ) |
(13.190) |
as a reference signal, applied to the zero measurement procedure, that, Vars0 = 0.5A02, and when the random variations of the amplification coefficient are absent, the variance of the variance estimate can be presented in the following form:
2 |
|
2 |
|
1 |
T |
|
|
τ |
2 |
|
|
|
Var {Vars*} = 2 Vars |
(1 + 2q + 2q |
|
) × |
|
|
1 |
− |
|
r |
|
(τ)dτ. |
(13.191) |
|
T |
|
|
|||||||||
|
|
|
|
∫ |
|
|
T |
|
|
|
||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
478 |
Signal Processing in Radar Systems |
As we can see from (13.200), it is not difficult to define the time interval of observation corresponding to the sensitivity threshold in the case of the correlation method of measurement of the stochastic process variance.
Comparing (13.198) through (13.200) and (13.162), (13.166), and (13.169), we can see that the sensitivity of the correlation method of the stochastic process variance measurement is higher compared to the compensation method sensitivity. This difference is caused by the compensation of noise components with high order while using the correlation method and by the compensation of errors caused by random variations of the amplification coefficients. However, when the channels are not identical and there is a statistical relationship between the intrinsic receiver noise and random variations of the amplification coefficients, there is an estimate bias and the variance of estimate also increases, which is undesirable.
The correlation method used for variance measurement leads to additional errors caused by the difference in the performance levels of the used mixers and the ideal mixers. As a rule, to multiply two processes, we use the following operation:
(a + b)2 − (a − b)2 = 4ab. |
(13.201) |
In other words, a multiplication of two stochastic processes is reduced to quadratic transformation of sum and difference of multiplied stochastic processes and subtraction of quadratic forms. The highest error is caused by quadratic operations.
Consider briefly an effect of spurious coupling between the amplifier channels on the accuracy of measurement of the stochastic process variance. Denote the mutual correlation functions of the amplification coefficients and intrinsic noise of amplifier at the coinciding instants as
β |
(t)β |
2 |
(t) = R |
, |
|
|
1 |
|
β12 |
|
|
|
|
|
|
|
(13.202) |
n1 |
(t)n2 (t) = Rn12 . |
||||
As we can see from (13.195), the bias of variance estimate is defined as
b{Var*} = Varn (1 + Rβ12 ) + Rβ12 Vars. |
(13.203) |
Naturally, the variance of estimate increases due to the relationship between the amplifier channels. However, the spurious coupling, as a rule, is very weak and can be neglected in practice.
13.5.4 Modulation Method of Variance Measurement
While using the modulation method to measure the stochastic process variance, the received realization x(t) of stochastic process is modulated at the amplifier input by the deterministic signal u(t) with audio frequency. Modulation, as a rule, is carried out by periodical connection of the investigated stochastic process to the amplifier input (see Figure 13.11) or amplifier input to the investigated stochastic processes and reference process (see Figure 13.12). After amplification and quadratic transformation, the stochastic processes come in at the mixer input, the second input of which is used by the deterministic signal u1(t) of the same frequency of the signal u(t). The mixer output process comes in at the low-pass filter input that make smoothing or averaging of stochastic fluctuations. As a result, we obtain the variance estimate of the investigated stochastic process.
Estimation of Stochastic Process Variance |
|
|
|
|
|
|
|
479 |
|||||||||
|
s(t) |
|
x(t) |
|
y1(t) |
|
|
y2(t) |
|
|
|
||||||
|
Amplifier |
Quadrator |
|
|
Low-pass |
|
Var*s |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
filter |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
u(t) |
|
u1(t) |
|
|
|
|
|
|
|
|
|
|||
|
|
|
Switching |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
FIGURE 13.11 Modulation method of variance measurement.
s(t) |
|
x(t) |
|
|
y1(t) |
|
|
y2(t) |
|
|
|
|
|
|
|
||||
|
|
Amplifier |
Quadrator |
|
|
Low-pass |
|||||||||||||
|
|
|
|
||||||||||||||||
s0(t) |
|
|
|
|
|
|
|
|
|
|
filter |
||||||||
u1(t) |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
Vars0 |
|
|
|
|
|
||||||||
|
u(t) |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
Switching |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Var*
s
FIGURE 13.12 Modulation method of variance measurement using the reference process.
Define the characteristics of the modulation method of the stochastic process variance measurement assuming that the quadratic transformer has no inertia. The ideal integrator with the integration time T serves as the low-pass filter. Under these assumptions and conditions we can write
|
1 at kT0 ≤ t ≤ kT0 + 0.5T0 ; |
|
||
|
|
|
|
|
|
|
at |
kT0 + 0.5T0 < t < (k +1)T0 ; |
(13.204) |
u(t) = 0 |
||||
|
|
= 0, ±1,…; |
|
|
|
k |
|
||
|
|
|
|
|
|
1 |
at kT0 ≤ t ≤ kT0 + 0.5T0; |
|
|
|
|
|
|
|
u1 |
|
|
kT0 + 0.5T0 < t < (k + 1)T0; |
(13.205) |
(t) = −1 at |
||||
|
|
= 0, ± 1,…. |
|
|
|
k |
|
||
|
|
|
|
|
As we can see from (13.204), u2(t) = u(t) and the average value in time of the function u(t) can be presented as u(t) = 0.5. Analogously we have
[1 − u(t)]2 = 1 − u(t). |
(13.206) |
Note that due to the orthogonality of the functions u(t) and 1 − u(t), their product is equal to zero. At first, we consider the modulation method of the stochastic process variance measurement pre-
sented in Figure 13.11. The total realization of stochastic process coming in at the quadrator input can be presented in the following form:
x(t) = [1 + υ(t)]{[s(t) + n(t)]u(t) + [1 − u(t)]n(t)}. |
(13.207) |
482 |
Signal Processing in Radar Systems |
When the random variations of the amplification coefficients are absent, based on (13.217), we can write
Var0 {Vars*} = |
4Var |
2 |
∞ |
|
(1+ 2q + 2q2 )∫r2 (τ)dτ. |
|
|||
|
s |
(13.219) |
||
T |
|
|||
|
|
|
0 |
|
Comparing (13.219) and (13.162) at Varβ = 0 while using the compensation method to measure variance, we can see that the relative value of the variance of the variance estimate defining the sensitivity of the modulation method is twice as high compared to the relative variation in the variance estimate obtained by the compensation method at the same conditions of measurement. Physically this phenomenon is caused by a twofold decrease in the time interval of observation of the investigated stochastic process owing to switching.
Now, consider the modulation method used for variance measurement based on a comparison between the variance Vars of the observed stochastic process ζ(t) within the limits of amplifier bandwidth and the variance Vars of the reference calibrated stochastic process s0(t) in accordance with the block diagram shown in0 Figure 13.12. In this case, the signal at the amplifier input can be presented in the following form:
x(t) = [1 + υ(t)]{[s(t) + n(t)]u(t) + [s0 (t) + n(t)] [1 − u(t)]}. |
(13.220) |
||||
The signal at the ideal integrator output can be presented in the following form: |
|
||||
z(t) = |
2 |
∫T [1 + υ(t)]2{[s(t) + n(t)]2 u(t) − [s0 (t) + n(t)]2 [1 − u(t)]}dt. |
(13.221) |
||
T |
|||||
0 |
|
|
|
||
The variance estimate is given by |
|
|
|
||
|
|
Var* = z + Var |
. |
(13.222) |
|
|
|
s |
s0 |
|
|
Averaging by realizations, we have |
|
|
|
||
|
|
Vars* = (1 + Varβ )(Vars − Vars0 ) + Vars0; |
(13.223) |
||
that is, the bias of variance estimate can be presented in the following form:
b{Vars*} = Varβ (Vars − Vars0 ). |
(13.224) |
When the random variations of the amplification coefficients are absent, the process at the modulation measurer output can be calibrated in values of difference between the variance of investigated stochastic process and the variance of reference stochastic process. This measurement method is called the modulation method of variance measurement with direct reading. When the variances
