Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Diss / 10

.pdf
Скачиваний:
143
Добавлен:
27.03.2016
Размер:
18.05 Mб
Скачать

Estimation of Stochastic Process Variance

473

stochastic process χ(t) and the low-frequency stochastic process θ(t) with zero mathematical expectations. In doing so, we assume that the stochastic process θ(t) possesses the unit variance, that is,θ2(t) = 1, and the correlation functions can be presented in the following form:

χ(t1)χ(t2 ) = Varn exp{−α|τ }cos ω0τ;

(13.177)

θ(t1)θ(t2 ) = exp{−η|τ }, τ = t2 t1.

(13.178)

Moreover, we assume that αT 1, α η, α γ. In this case, the variation in the variance estimate or the variance of the signal at the ideal integrator output can be presented in the following form:

Var {Vars*} =

(1+ q)

2

 

2

 

+ 8 α Varβ

γ T − 1+ exp{−γT}

 

 

 

 

 

 

Vars

1

 

 

 

 

αT

 

 

γT

 

 

 

 

 

 

 

 

γ

 

 

 

 

 

 

 

 

2

 

T − 1+ exp{−2ηT}

+ 16Varβ

(2η + γ )T − 1+ exp{−(2η + γ )T}

 

+ Varn

 

 

 

 

2

 

2

 

 

 

 

 

.

 

 

 

 

T

 

 

2

T

2

 

 

 

 

 

 

η

 

 

 

 

(2η + γ )

 

 

(13.179)

Comparing (13.179) with (13.168), we see that owing to the low-frequency fluctuations of amplifier intrinsic noise the variance of the variance estimate increases as a result of the value yielded by the second term of (13.179).

13.5.2  Method of Comparison

This method is based on a comparison of constant signal components formed at the output of twochannel amplifier by the investigated stochastic process and the generator signal coming in at the amplifier input. Moreover, the generator signal is calibrated by power or variance. The wideband Gaussian stochastic process with known power spectral density or the deterministic harmonic signal can be employed as the calibrated generator signal. Flowchart of stochastic process variance measurer employing a comparison of signals at the low-pass filter output is presented in Figure 13.9.

Assuming that the receiver channels are identical and independent, the signal at the squarer inputs can be presented in the following form:

x (t)

1

x2 (t)

= [1 + υ1(t)] [s(t) + n1(t)];

(13.180)

= [1 + υ2 (t)] [s0 (t) + n2 (t)],

s(t)

 

x

1

(t)

 

x2(t)

 

y1(t)

 

Vars0

 

 

Amplifier

 

 

 

Quadrator

1

 

Low-pass

 

 

 

 

 

 

 

 

 

filter

+

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x22(t)

 

 

 

 

 

 

z(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Amplifier

x2(t)

Quadrator

Low-pass

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s0(t)

 

 

 

 

 

 

filter

y2(t)

 

 

Var*s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

FIGURE 13.9  Flowchart of the stochastic process variance measurer.

474

Signal Processing in Radar Systems

where

υ1(t) and υ2(t) are the realizations of Gaussian stationary stochastic processes describing the random variations β1(t) and β2(t) relative to the amplification coefficient value at the first and second channels of signal amplifier

n1(t) and n2(t) are the realizations of signal amplifier noise at the first and second channels s(t) and s0(t) are the investigated stochastic process and the reference signal, respectively

Let the reference signal be the stationary Gaussian stochastic process with the variance Vars0. The signal at the subtractor output takes the following form:

 

1

T

2

2

 

 

 

 

 

 

 

z(t) =

T x1

(t) x2

(t) dt.

(13.181)

 

 

0

 

 

 

 

The variance estimate is determined in the following form:

Var*(t) = z(t) + Var

.

(13.182)

s

s0

 

 

Define the mathematical expectation and the variance of the variance estimate under assumptions made earlier

 

2

(t)

=

2

(t)

= Varβ ;

β1

β2

 

 

 

=

 

 

(13.183)

 

2

(t)

2

(t)

= Var .

 

n1

 

n2

n

The mathematical expectation of variance estimate can be presented in the following form:

Vars* = z + Vars0 = (1 + Varβ )(Vars − Vars0 ) + Vars0 .

(13.184)

As we can see from (13.184), the bias of variance estimate can be presented in the following form:

b{Vars*} = Varβ (Vars − Vars0 ).

(13.185)

If the variance of the reference stochastic process is controlled then the variance measurement process can be reduced to zero signal aspects at the output indicator. This procedure is called the method with zero instant. Naturally, the variance estimate is unbiased.

The variance of the variance estimate can be presented in the following form:

Var {Vars*} =

 

4 Varβ

T

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

1

 

 

rβ (τ)[2 + rβ (τ)Varβ ]

 

 

 

T

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×{(Vars + Varn)2 + (Vars0

+ Varn )2 + 2{[Rs (τ) + Rn (τ)]2 + [Rs0 (τ) + Rn (τ)]2}}dτ

 

 

 

 

 

2

4

T

 

 

τ

2

2

 

 

+ (1 + Varβ )

 

 

 

 

 

 

1

 

{[Rs (τ) + Rn (τ)]

+ [Rs0 (τ) + Rn (τ)] }dτ.

(13.186)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

T

 

 

 

0

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

 

 

 

 

 

 

 

476

Signal Processing in Radar Systems

As applied to the exponential normalized correlation function given by (13.165), we can transform (13.191) into the following form:

Var {Vars*} = 2Vars2 (1 + 2q + 2q2 ) T − 1 + exp{−2αT} .

(13.192)

 

(2αT )2

 

At αT 1, we have

 

 

 

Var {Vars*}

Vars2 (1 + 2q + 2q2 )

.

(13.193)

 

 

αT

 

Comparing (13.193) with (13.189) at Varβ = 0, we can see that in the case of weal signals, that is, the signal-to-noise ratio is small or q 1, with the use of the harmonic signal the variance of the variance estimate is twice as less as when the wideband stochastic process is used. Otherwise, at the high signal- to-noise ratio, q > 1, under the use of the harmonic signal the variance of the variance estimate can be higher compared to the case when we use the stochastic process. This phenomenon is explained by the presence of uncompensated components of the output signal formed by the terms of high order under expansion of the amplifier noise in series. In addition, the considered procedure possesses the errors caused by the nonidentity of channels and presence of statistical dependence between amplifier channels.

13.5.3  Correlation Method of Variance Measurement

While using the correlation method to measure the stochastic process variance, the stochastic process comes in at the inputs of two channels of amplifier. The intrinsic amplifier channel noise samples are independent of each other. Because of this, their mutual correlation functions are zero and the mutual correlation function of the investigated stochastic process is not equal to zero and the coincidence instants are equal to the variance of the investigated stochastic process. A flowchart of the variance measurement made using the correlation method is presented in Figure 13.10.

Let the independent channels of amplifier operate at the same frequency. Then the stochastic processes forming at the amplifier outputs can come in at the mixer inputs directly. We assume, as before, that the integrator is the ideal, that is, h(t) = T−1. The signal at the integrator output defines the variance estimate

 

 

 

 

 

Vars* =

1

T

x1(t)x2 (t)dt,

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

where

 

 

 

 

 

 

 

 

 

 

 

x1(t) = [1 + υ1(t)] [s(t) + n1(t)];

 

 

 

(t) = [1 + υ2 (t)] [s(t) + n2

 

 

 

 

 

x2

(t)]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1(t)

 

 

 

 

s(t)

 

 

Amplifier

 

 

Var*s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Integrator

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Amplifier

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.194)

(13.195)

FIGURE 13.10  Correlation method of variance measurement.

Estimation of Stochastic Process Variance

477

are the realizations of stochastic processes at the outputs of the first and second channels, respectively. The mathematical expectation of estimate can be presented in the following form:

Vars* =

1

T

s2 (t) dt = Vars ;

(13.196)

T

0

 

 

that is, the variance estimate is unbiased under the correlation method of variance measurement and when the identical channels are independent.

The variance of the output signal or the variance of the variance estimate of the input stochastic process in the case when the input stochastic process is Gaussian takes the following form:

Var {Vars*} =

2Vars2

 

T 1 −

τ

 

Rβ

(τ) + Rβ

 

(τ) + Rβ

(τ)Rβ

 

(τ) dτ

 

 

 

 

 

 

2

2

 

 

 

 

T

 

 

 

T

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

2 T

1

τ

1+ Rβ (τ) + Rβ

 

(τ) + Rβ (τ)Rβ

 

(τ)

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

T

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

2Rs2

(τ) + Rn

(τ)Rs (τ) + Rn (τ)Rs (τ) + Rn (τ)Rn

(τ) dτ,

(13.197)

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

1

 

 

2

 

 

where Rβ1(τ) and Rβ2(τ) are the correlation functions of random components of the amplification coefficients of the first and second channels of amplifier.

We can simplify (13.197) taking into consideration that the correlation interval of random components of amplification coefficients is longer compared to the correlation interval of the intrinsic noise n1(t) and n2(t) and the investigated stochastic process ς(t). At the same time, we can introduce the function Rβ(0) = Varβ instead of the correlation function Rβ(τ) in (13.197). Then, taking into consideration the condition accepted before, that is, Varβ < 1 or Varβ2 1, identity of channels, and (13.153) through (13.155), we obtain

Var {Vars } =

2Vars2

(1+ q + 0.5q

2

 

T

τ

2

(τ)dτ +

4Vars2Varβ

 

 

)

1

 

r

 

 

T

 

 

 

T

 

 

 

 

 

T

 

 

 

 

 

 

 

0

 

 

 

 

 

 

T

1

τ

r

(τ)dτ. (13.198)

 

 

 

 

 

β

 

 

T

 

 

0

 

 

 

 

 

Formula (13.198) allows us to define the time interval of observation T corresponding to the sensitivity threshold at the condition given by (13.188). Let the normalized correlation functions r(τ) and rβ(τ) be defined by the exponents given by (13.165) and substituting into (13.198), we obtain

2

 

2

 

T − 1+ exp{−2αT}

2

 

γ T − 1+ exp{−γT}

 

Var {Vars*} = 2Vars

(1+ q + 0.5q

 

)

(2αT )2

+ 4Vars

Varβ

 

.

 

T )2

(13.199)

When the random variations of the amplification coefficients are absent and the condition T τcor is satisfied, the variance of the variance estimate can be presented in the following form:

Var0

{Vars*} = Vars2 1+ q + 0.5q2 .

(13.200)

 

αT

 

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)

480

 

 

Signal Processing in Radar Systems

The process at the quadrator output takes the following form:

 

 

y1(t) = x2 (t) = [1 + υ(t)]2{[s(t) + n(t)]2 u(t) + [1 − u(t)]n2 (t)}.

(13.208)

The process at the mixer output takes the following form:

 

 

y2 (t) = [1 + υ(t)]2{[s(t) + n(t)]2 u(t) − [1 − u(t)]n2 (t)}.

(13.209)

The signal at the ideal integrator output (the variance estimate) takes the following form:

 

 

Vars* =

2

T [1 + υ(t)]2{[s(t) + n(t)]2 u(t) − [1 − u(t)]n2 (t)}dt.

(13.210)

 

T

 

0

 

The factor 2 before the integral is not important, in principle, but it allows us, as we can see later, to obtain the unbiased variance estimate and the variance of the variance estimate that are convenient to compare with the variance of the variance estimates obtained by other procedures and methods.

By averaging the random estimate Vars*, we obtain

 

Vars* = Vars (1 + Varβ ).

(13.211)

As we can see from (13.211), when the random variations of the amplification coefficients are absent Varβ = 0, the variance estimate of stochastic process using the modulation method of variance measurement is unbiased as is the case with the compensation method. To determine the variance of the variance estimate we need to present the function u(t) by the Fourier series expansion

 

 

 

 

 

 

 

u(t) =

1

+

2

sin(2k − 1)Ωt

,

(13.212)

2

π

 

 

k =1

2k − 1

 

 

 

 

 

 

 

 

 

where

 

 

 

 

 

 

 

 

 

 

Ω =

 

(13.213)

 

 

 

 

 

T0

 

 

is the switching function frequency.

Ignoring oscillating terms and taking into consideration that Varβ 1 and the correlation interval of random component of the amplification coefficients is much more the correlation ­interval of the investigated stochastic process, the variance of the variance estimate, can be presented in the following form:

Var {Vars } =

32

1

 

T

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

π T

k =1

(2k − 1)

 

0

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

2

 

 

× { Rs (τ) + 2Rs (τ)Rn (τ) + 4Rn

(τ)

+ 2Rβ (τ)(Vars

+ 4VarsVarn + 4Varn )}

 

 

 

 

 

 

 

 

 

4

 

T

 

 

 

τ

2

 

2

 

 

 

× cos[(2k − 1)Ωτ]dτ +

 

 

 

 

1

 

 

 

(τ) + 2Rs (τ)Rn (τ) + 2Rβ (τ)Vars

τ

 

T

 

 

 

 

Rs

d

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

(13.214)

Estimation of Stochastic Process Variance

481

Taking into consideration (13.153) through (13.155), (13.214) can be written in the following form:

Var {Vars } =

16Vars2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π2T

 

 

(2k − 1)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

 

T

 

τ

 

+ 2q +

 

2

 

2

(τ) + 4Varβrβ (τ)(1 + 4q + 4q

2

 

 

 

1

 

 

 

4q

 

)r

 

 

 

 

 

 

(1

 

 

 

)

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2Vars2 T

 

τ

 

 

 

 

 

 

2

 

 

 

 

+

 

 

 

 

 

 

 

1 −

 

 

(1

+ 2q)r

 

(τ) + 4Varβrβ (τ) dτ.

 

 

 

 

T

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos[(2k − 1)Ωτ]dτ

(13.215)

In practice, the correlation interval of the investigated stochastic process is much less than the modulation period T0 and, consequently, the time interval T of observation, since in real applications, the effective spectrum bandwidth of the investigated stochastic process is more than 104–105 Hz and at the same time, the modulation frequency fmod = Ω × (2π)−1 is several hundred hertz. In this case, we can assume that the functions cos[(2k − 1)Ωτ] are not variable functions within the limits of the correlation interval, and we can think that this function can be approximated by unit, that is, cos[(2k − 1)Ωτ] ≈ 1. This statement is true for the components of variance of the variance estimate caused by stochastic character of variation of the processes s(t) and n(t). Taking into consideration that

π2

 

 

(2k − 1)−2 =

,

(13.216)

8

k=1

 

 

 

 

 

we obtain

Var {Vars*} =

4Var2

 

 

 

 

 

 

 

 

 

 

8Var2Var

 

 

(1 + 2q + 2q

 

)r

 

(τ)dτ +

 

 

 

T

 

s

 

2

2

 

 

T

β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

τ

 

 

8(1 + 4q + 4q

2

)

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

 

 

 

 

1 −

 

 

rβ (τ)dτ +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

2

 

 

 

(2k − 1)

2

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

T

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

 

 

1

 

 

rβ (τ)cos[(2k − 1)Ωτ]dτ dτ .

 

 

(13.217)

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

The approximation cos[(2k − 1)Ωτ] ≈ 1 is true within the limits of the correlation interval of investigated stochastic process owing to fast convergence of the series (13.216) to its limit. If, in addition to the aforementioned condition, the correlation interval of the amplification coefficients is much less than the time interval of observation, that is, τβ 1 but τβ > T0, then

Var {Vars*}

4Var

2

 

2

2

 

8Var2Var

 

s

(1+ 2q + 2q

)r

(τ)dτ +

s

β

rβ (τ)dτ.

(13.218)

T

 

 

 

T

 

 

 

 

 

 

0

 

 

 

 

0

 

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

Соседние файлы в папке Diss