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

Diss / 10

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

Global Digital Signal Processing System Analysis

343

In accordance with the earlier discussion, a distribution of time intervals between the queuing requests served by the detector–selector (the detector–selector input) can coincide with distribution of the request queue time in the detector–selector. The parameter of request queue at the detector– selector output coincides with the parameter of the incoming request queue.

10.2.3  Analysis of MTI Characteristics

The MTI is considered as the one-channel queuing system with waiting time. The request queue served by the detector–selector with the parameter γin and the pdf of time intervals between the requests coinciding with the pdf of request queue time in the detector–selector given by (10.14) and

(10.15) comes in at the MTI input. The request queue time in MTI is constant and equal to τMTI = a. Depending on the relation between the constant constituent τ0 in (10.36) and the duration of the

request queuing time interval in MTI τMTI = a, the following cases are possible.

Case 1: a ≤ τ0. MTI can forever serve the request queue before the next request queue comes in, and the downtime interval is even possible. The average duration of downtime depends on the variance of the random component ξi of the time interval between the requests at the MTI input. If the pdf of time intervals between the requests coming in at the MTI input is subjected to the exponential pdf with the shift given by (10.14), then the average downtime is defined as

 

down = τDS a ≥ στDS .

(10.38)

t

If the pdf of time intervals between the requests coming in at the MTI input is subjected to the pdf given by (10.15), we have

 

down = τDS a

τDS

.

(10.39)

t

 

 

 

 

Case 2: τ0 < a ≤ γ in−1. There is a request queue at the MTI input. It is impossible to compute the request queue length or the waiting time using analytical methods because the input request queue is not simple. To define the mentioned characteristics there is a need to apply a simulation. The greatest difficulty is to simulate the request queue with the pdf given by (10.14) or (10.15). The curves of relative waiting time for request queue by MTI as a function of the loading factor and the shift coefficient α = τ0στDS1 as the parameter are shown in Figure 10.7. As we can see from Figure 10.7,

tMTIwait

 

 

 

a

 

 

 

 

 

5

 

 

 

 

 

 

 

1

2

 

 

 

1

α = 1

 

 

 

 

 

 

2

α = 2

 

 

 

 

4

 

 

3

α = 3

 

 

 

3

 

 

 

4

α = 4

 

 

 

 

 

 

 

 

 

 

5

α = 5

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

χMTI

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0.8

 

0.9

1.0

 

 

FIGURE 10.7  Relative waiting time versus the loading factor at different values of shift coefficient α.

344

Signal Processing in Radar Systems

the average waiting time at the MTI input under the fixed loading factor is a function of the shift coefficient α. If α ≥ 3, the input request queue degenerates into a regular request queue and the MTI request queue time is approximately equal to τ0. In this case, the request queue at the MTI input is absent and the input register is used as the buffer memory.

10.3  ANALYSIS OF “n n – 1” MTI SYSTEM

The “n – n – 1” MTI system is presented in Figure 10.4. For this system the request queuing by the detector–selector is started immediately after filling out the memory matrix. If, as before, we assume that the pdf of time to fill out the memory is given by (10.6), then the pdf of request queuing time by each channel of the system memory–detector–selector is a combination of the pdfs given by (10.6) and (10.14) or (10.6) and (10.15) and the average request queuing time is defined as

τmemory− DS = τmemory + τDS,

(10.40)

where

τmemory is the average time to fill out the memory

τDS is the average time for request queuing by the detector–selector given by (10.17) or (10.29)

Now, if the request queuing parameter is given (the input request queuing is considered as the simplest, as before), the number of channels of the system memory–detector–selector can be deter-

mined using the Erlang formula (10.18) at the given probability of failure Pfailure.

Let us discuss the request queuing pdf at the output of n-channel system memory–detector– selector. Action of the considered system memory–detector–selector on the incoming request queue can be presented as an expansion of the simplest request queue on the elementary request queues, the number of which is equal to the number of channels of the system memory–detector–selector. In a general case, these elementary request queues may not be the simplest. The output request queue of the system memory–detector–selector is a superposition of elementary request queues and can be considered as the simplest request queue with the parameter equal to the parameter of the incoming

or input request queue at low values of Pfailure based on Sevastyanov’s theorem [5]. The request queue time in MTI is constant and equal to a, as earlier. Therefore, the average waiting time for request

queue by MTI is defined as

 

wait

=

γ ina2

=

χMTIa

.

 

tMTI

 

 

2(1 − γ ina)

2(1 − χMTI )

 

 

 

 

 

The average number of request in queue is given by

Nwait = χMTI 3 2χMTI ,

2(1 − χMTI )

and the variance

2

 

 

 

χMTI

 

1

 

1

 

χMTI

 

 

 

 

 

 

 

 

 

DNwait = σNwait

= χMTI 1

+

 

 

 

 

+ χMTI

 

+

 

 

 

,

1 − χMTI

2

3

 

 

 

 

 

 

 

4(1 − χMTI )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

where χMTI is the MTI loading factor.

(10.41)

(10.42)

(10.43)

Global Digital Signal Processing System Analysis

345

Knowing Nwait and σ2Nwait, we can define the buffer memory capacity at the MTI input under the given acceptable probability to loss the request and assuming, for example, that the pdf of the

number ofrequests in queue is the normal Gaussian. This approximation makes sense only at high values of Nwait. In this case, the probability of losing the request at the MTI input is defied as [6]

 

 

 

 

 

 

 

 

 

loss

 

1

 

(Nwait Nwait )2

 

PMTI

=

 

exp −

 

 

 

dNwait ,

(10.44)

2

2

 

 

 

 

QBM

2πσNwait

 

Nwait

 

 

 

 

 

 

 

 

 

 

 

where QBM is the required buffer memory. There is a need to note that the acceptable probability of losing the request at the MTI input PMTIloss must be, at least, on the order less than the acceptable probability of losing the request at the system input. Only in this case the main requirement to the system operation will be satisfied, namely, all queue requests that pass the first phase of queue must be served by MTI.

10.4  ANALYSIS OF “n m – 1” MTI SYSTEM

The considered system is the three-phase queuing system with losses at the input. The first phase is the n-channel queuing system with losses. Definition of the required number of channels of this system is carried out according to the procedure discussed earlier. The second phase is the m-channel queuing system with waiting. Analysis of QoS factors of this system is the subject of the present section. The input request queue is the simplest one with the parameter γin. We assume that the splitter operates as a counter with respect to the base m sending to the ith channel, i = 1, 2,…, m, the requests with numbers i, i + m, i + 2m. This allocation method of requests is called the cyclic way. The simplest request queue thinned in m − 1 times, i.e., the Erlang request queue of the (m − 1)th order, with the following pdf

pm−1(t) =

γ in int)m−1 exp(−γ int)

(10.45)

(m − 1)!

 

 

comes in at the input of each channel of the one-channel queuing system. The condition of the m-channel queuing system stationary mode is the following:

m

 

χi = χm = γ in τDS < m.

(10.46)

i=1

At this time we assume that the average request queue time τDS is the same for all channels. Thus, in the considered case, the one-channel queuing system with the Erlang incoming request queue possessing the pdf equal to γinm−1 and the request queue time subjected to the pdf given by (10.14)

and (10.15) with the parameters τ and σ2 is investigated. The result of investigation must be the

DS τDS

statistical characteristics of the request queue waiting time, namely, the average time twait and the

variance of this time σ2twait.

As noted earlier, the analytical investigation of queuing systems with the incoming request queue different from the simplest one is a complicated problem. A general approach to solving this problem

is discussed in Ref. [5]. However, an implementation of this general approach to specific conditions of

2

our problem is difficult, not obvious, and does not lead to final formulae for twait and σtwait. The explicit

solutions for twait / τDS obtained by numerical simulation are shown in Figure 10.8, where the solid lines represent the exponential with shift pdf of request queue time and the dashed lines represent the truncated and shifted normal Gaussian pdf. As we can see from Figure 10.8, with an increase in m the

346

Signal Processing in Radar Systems

twait τDS

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

1

α = 0.5

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

α = 1.0

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

3

α = 2.0

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

4

α = 3.0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

1

2

3

4

5

6

 

 

 

 

 

FIGURE 10.8  Average waiting time versus the number of queuing system channels.

relative average waiting time decreases exponentially. The shift coefficient α affects the average waiting time by an analogous way as for previously discussed systems. At the given parameters γin, m, α, we can compute τDS based on the curves presented in Figure 10.8. We are able to define the average waiting time of request queue and, taking into consideration the average waiting time of request queue, we are able to define the required number of channels n for the first-phase queuing system.

Dependences of the required speed of operation for each mth channel of the second-phase queuing system for several values ofthe input request queue parameter γin as a function of the average number of reduced operations N required to process a single request under the loading factor of each mth channel χi = 0.9m are shown in Figure 10.9 for the truncated and shifted normal Gaussian pdf given by (10.15). The required speed of operation is computed by the following formula:

 

 

 

γ in

 

 

Vef =

N

.

(10.47)

 

 

 

 

0.9m

 

Analysis of (10.47) and dependences in Figure 10.9 show that at the fixed parameters of the input request queue and the already-given number of operations for a single realization of digital signal processing algorithm, the required speed of operation of the second-phase queuing system decreases proportionally to the number of channels in this system.

A set of requests processed by the second-phase queuing system form the incoming request queue for the third-phase queuing system that is a one-channel queuing system with waiting. According to the limiting theorem for summarized flows, the incoming request queue at the third-phase queuing system input is very close to the simplest one with the parameter γin. The request queue time in the third-phase queuing system is constant, as earlier. Therefore, a computation of required characteristics of the three-phase queuing system and buffer memory capacity is carried out by the procedure discussed in Section 10.3.

Global Digital Signal Processing System Analysis

347

Vef , operations/s

105

1 2 3

1

2

3

104

m= 3 m= 5

1 γin = 50

2 γin = 40

3 γin = 30

103

0

1

2

3

4

5

6

7

8

 

 

 

 

N × 10–3

 

 

 

FIGURE 10.9  Required speed of operation for each mth channel of the second-phase queuing system for several values γin versus the average number of operations under the loading factor χ i = 0.9m; truncated and shifted normal Gaussian pdf.

10.5  COMPARATIVE ANALYSIS OF TARGET TRACKING SYSTEMS

The total QoS factors of discussed versions of the target tracking systems are the following:

The required number of input memory channels under the given probability of failure

Pfailure

The required effective speed of operation

The average time required to process the request by the queuing system

Based on the discussions provided in the previous sections, we are able to compare the considered target tracking systems under the following initial data:

• The request queue at the target tracking system input is the simplest with γin = 40 s−1.

• The acceptable probability of failure at the target tracking system input Pfailure = 10−3.

• The loading factor of the detector–selector and MTI is χDS (MTI) = 0.9.

• The pdf of request queue time in the detector–selector is given by (10.15) with α = 2.

• The average number of reduced short operations required to process a single request by the

detector–selector system is NDS = 1500.

• The number of reduced short operations required to process a single request by the MTI

is NMTI = 2500.

Scanning by azimuth with the uniform velocity V0 = 250° s−1.

Minimal size of the target tracking gate by azimuth is Δβmingate = 2 .

Minimal time to fill out the memory matrix (the constant constituent of the request queue

min

min

.

time in memory) τmemory =

βgate /V0

= 0.008 s

348 Signal Processing in Radar Systems

If we assume in (10.6) α = τminmemory τmemory = 2, then the average request queue time in the memory is defined as

τmemory = τmemorymin +

τ

= 0.011 s.

(10.48)

 

 

 

To compare the target tracking system of various types we define the following:

Number of channels (matrices) of the input memory Nmemory

Average waiting time at the detector–selector input tDSwait

Average request queue time in the detector–selector τDS

Required speed of operation of the detector–selector VDSef

Number of detector–selectors NDS

Average waiting time at the MTI input tMTIwait

Request queue time in MTI τMTI

Number of the buffer memory cells at the MTI input NBM (MTI)

Required speed of MTI operation VMTIef

Average request delay by queuing system τΣ

Computational results are presented in Table 10.1, which shows that the system “n – 1 – 1” ensures the minimal average time of signal processing by one target equal to τΣ = 0.16s at the speed of operation of MTI equal to VMTIef = 150 × 103 operations per second and the speed of operation of detector–selector equal to

VDSef = 66.7 × 103 operations per second, i.e., providing a refreshment on nine targets in the course of each scanning. The same average time of signal processing is provided by the system “n – n – 1” when there are six detector–selectors with the speed of operation equal to VDSef = 66.7 × 103 operations per second but at the lesser speed of MTI operation. Obviously, a realization of this version is expensive in comparison with the system “n – 1 – 1.” Other versions of the system can be realized using the detector–selectors with the low speed of operation that can be a clincher in their favor [7]. The system “n – m – 1” can also be realized using the detector–selectors with a low speed of operation. In doing so, the number of detector–selectors is not high (two or three), which is an advantage of this system in comparison with the second type of system. The disadvantage is an increase in the request queue time.

Selection of definite version of the system is defined by computational resources and acceptable time of signal processing. Evidently, the most effective system from considered ones is the system “n – m – 1” with two or three detector–selectors.

TABLE 10.1

Comparative Analysis of Target Tracking Systems

 

 

 

Detector–Selector

 

 

 

 

 

MTI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

VDSef , 1000

 

 

 

 

MTwaitI

 

VMTIef , 1000

 

 

 

 

 

 

 

 

 

 

t

 

 

System

Nmemory

tDS s

 

 

DSwait s

Operations/s

 

NDS

τMTI s

 

s

NBM (MTI)

Operations/s

τΣ s

 

t

 

n – 1 – 1

14

0.0225

0.107

66.7

1

0.0167

0

1

150

0.16

n – n – 1

6

0.0225

0

66.7

6

0.0225

0.1

15

110

0.16

 

8

0.045

0

33.3

8

0.0225

0.1

15

110

0.18

 

15

0.15

0

10

15

0.0225

0.1

15

110

0.28

n – m – 1

14

0.045

0.12

33.3

2

0.0225

0.1

15

110

0.3

 

14

0.0685

0.123

22.2

3

0.0225

0.1

15

110

0.325

 

13

0.09

0.112

16.7

4

0.0225

0.1

15

110

0.335

 

13

0.113

0.107

13

5

0.0225

0.1

15

110

0.355

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Global Digital Signal Processing System Analysis

349

10.6  SUMMARY AND DISCUSSION

Detection and tracking of all targets in the all-round surveillance radar coverage with accuracy sufficient to reproduce and estimate data of the air situation are carried out by the so-called rough channel of global digital signal processing system covering the whole scanned area. The rough channel consists of the following subsystems: the binary signal quantization, the specific microprocessor network for target return signal preprocessing, and the microprocessor network for digital signal reprocessing and control. Employment of binary quantization and simplified versions of digital signal processing algorithms allows us to implement the microprocessor networks and sets in this channel of global digital signal processing system of the all-round surveillance radar with the uniform antenna rotation.

The detection and selection problem of a single target within the limits of the gate is assigned in the following way based on information contained in the likelihood ratio maxima. First, we take the hypothesis that there is only a single target within the limits of the gate. In the considered case, the event, when several targets are within the limits of the gate, is possible but it is improbable. Using the relief of the two-dimensional likelihood ratio surface with M peaks of different heights, there is a need to define the maximum (the peak) formed by the target return signal and, if a “yes,” there is a need to define the coordinates of this two-dimensional likelihood ratio surface maximum. The amplitudes of peaks Zl(l = 1, 2,…, M) and their coordinates ξl and ηl with respect to the gate center are used as the input parameters, based on which the decision is made. If the hypothesis about the statistical independence of the two-dimensional likelihood ratio surface peak amplitudes is true and the coordinates of the two-dimensional likelihood ratio surface maxima are known within the limits of range where the target return signal is present, the optimal detection–selection problem of the target return signal pips within the limits of gate is solved in two steps.

All considered and discussed versions of target tracking by several MTI systems are threephase queuing systems. Each phase is represented by one or several devices of queuing system connected in series. Generally, the request queuing time for each device is the random variable with pdf available to investigate and define. It is assumed that there is a simple request queue at first-phase input of the queuing system. The request queue coming in at the first device input is immediately served if, at least, one channel of the queuing system is free from service; otherwise, the request is rejected and flushed. The requests carried out by the first device are processed sequentially at the next phases; in other words, the request loss at each next phase or stage is inadmissible. In line with this, the memory device with the purpose storing requests in line should be provided before the second and third phases of the queuing system.

Difficulties under analysis of multiphase queuing systems are the following. At all cases of practical importance, the output stream of phase takes a more complex form in comparison with the incoming request queue. In some cases, the output request queue can be approximated by the simplest incoming stream with the same parameters. Then we can use analytical procedures and methods of the queuing theory to analyze the next phase or stage. If this approximation is impossible, then the only method to investigate the stream is the simulation. A rational combination of analytical and simulation methods and procedures allows us to solve the problem of the three-phase signal processing system analysis using the MTI system for any design and construction version.

REFERENCES

1.Lyons, R.G. 2004. Understanding Digital Signal Processing. 2nd edn. Upper Saddle River, NJ: Prentice Hall, Inc.

2.Harris, F.J. 2004. Multirate Signal Processing for Communications Systems. Upper Saddle River, NJ: Prentice Hall, Inc.

3.Barton, D.K. 2005. Modern Radar System Analysis. Norwood, MA: Artech. House, Inc.

350

Signal Processing in Radar Systems

4.Skolnik, M.I. 2008. Radar Handbook. 3rd edn. New York: McGraw-Hill, Inc.

5.Gnedenko, V.V. and I.N. Kovalenko. 1966. Introduction to Queueing Theory. Moscow, Russia: Nauka.

6.Skolnik, M.I. 2001. Introduction to Radar Systems. 3rd edn. New York: McGraw-Hill, Inc.

7.Hall, T.M. and W.W. Shrader. 2007. Statistics of clutter residue in MTI radars with IF limiting, in IEEE Radar Conference, April 2007, Boston, MA, pp. 01–06.

Part III

Stochastic Processes

Measuring in Radar Systems

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