Literature DB >> 27504217

Calibration of Speed Enforcement Down-The-Road Radars.

John Jendzurski1, Nicholas G Paulter1.   

Abstract

We examine the measurement uncertainty associated with different methods of calibrating the ubiquitous down-the-road (DTR) radar used in speed enforcement. These calibration methods include the use of audio frequency sources, tuning forks, a fifth wheel attached to the rear of the vehicle with the radar unit, and the speedometer of the vehicle. We also provide an analysis showing the effect of calibration uncertainty on DTR-radar speed measurement uncertainty.

Entities:  

Keywords:  DTR radar; calibration; down the road radar; radar; traffic radar; uncertainty

Year:  2009        PMID: 27504217      PMCID: PMC4646567          DOI: 10.6028/jres.114.009

Source DB:  PubMed          Journal:  J Res Natl Inst Stand Technol        ISSN: 1044-677X


1. Introduction

Down-the-road (DTR) radar (radio detection and ranging) is the most common type of speed measurement device used in the United States and in the rest of the world for traffic speed enforcement. There are approximately 150 000 DTR radar units in use in the United States, and almost every legal jurisdiction in the United States accepts DTR radar for evidentiary purposes in speed enforcement cases in traffic court. Consequently, the performance of these devices must be sufficient to assure the courts’ confidence in using these devices. This assurance comes from 1) compliance to established minimum performance requirements, 2) scheduled calibration by testing labs to these requirements, and 3) routine calibration in the field by the radar operator. DTR radars are provided as either hand-held radar “guns,” dashboard-mounted units, or similar vehicle-mounted units. The DTR radar has three basic components: the transceiver/mixer head, the processing component, and the display/human interface component. DTR radars typically are found in one of two forms, one in which the display and transceiver head are combined into a common unit, and another form in which the two heads are separate but connected by an appropriate electrical interconnect. DTR radars operate within narrowly allocated (by the Federal Communications Commission) spectral bands in one of the following three broad frequency bands: X (nominally 8 GHz to 12 GHz in the United States), K (nominally 18 GHz to 26 GHz in the United States), and Ka (nominally 26 GHz to 40 GHz in the United States). The transceiver head contains the microwave source that illuminates the target vehicle and the mixer that is used to detect the Doppler-shifted microwave signal reflected from the target vehicle. As will be discussed in the next section, the signal output from the mixer circuit (which contains a diode) is at a frequency that is the difference between the Doppler-shifted signal and the incident (illuminating) signal. This difference frequency is much lower than the frequency of the incident or the Doppler-shifted signals, and is the signal upon which a speed measurement is based. The primary purpose of the processing part of the DTR unit is to compute the speed of a target (or patrol) vehicle from the difference frequency. The display/human interface provides a readout of the computed speed, operator controls, status indicators, etc. In this paper we are concerned with the uncertainty contributions to speed measurement using DTR radar that can be attributed to the DTR radar calibration process. This analysis does not consider uncertainties in speed measurement during actual use, that is, that are caused by operational issues with DTR radar usage.

2. Background

As just mentioned, DTR radar uses the Doppler frequency shift, Δf, in the frequency of the reflected microwave signal that is used to determine speed. The value of Δf is dependent on v, which is the radial component of the relative velocity, v, and f0, which is the frequency of the microwave source. Figure 1 shows the relationship between v and its radial and tangential components, v and v. It is worth noting that v = v only when θ = 0, that is when the target is moving directly towards the patrol vehicle. For all other values of θ, v < v.
Fig. 1

Radial and tangential components (v and v) of the relative velocity. This diagram is not drawn to scale; typical values for θ are on the order of a few degrees for normal DTR radar operation

Since the DTR radar measures v, for the remainder of this paper, the subscript “rad” will not be explicitly shown, but it is implied that v is the radial velocity for all subsequent mentions and calculations. To compute v, we must first have a value for Δf. To compute Δf, we start with the basic formula equating the frequency, f0, of an electromagnetic wave propagating in air to the speed of light, c, in air, and the wavelength, λ, of that wave: The moving vehicle causes the observed wavelength to get shorter or longer, depending on whether the vehicle is moving toward the source or away from it. This change in λ, Δλ, is given by: where the factor of 2 arises because of the round trip propagation of the microwave signal from the patrol vehicle to the target vehicle and then back to the patrol vehicle. The corresponding change in frequency, Δf, is given by: where f is the Doppler frequency and the approximation in (3) is valid for v ≪ c, which, hopefully, is the case for ground-based vehicles. The Doppler shifted microwave signal, at the frequency f, that is reflected from the target vehicle is collected by the transceiver and mixed with the source signal, at frequency f0, in a diode, to yield the signal at Δf. (Signal components at other frequencies are also present, but their amplitudes are much smaller than the signal at Δf.) To understand this mixing, we start with the well-known diode equation (see Ref. [1], page 126), which is: where I is the saturation current, q is the electronic charge, k is Boltzmann’s constant, T is the temperature in Kelvin, n is a constant for a given diode type, and V is the voltage applied to the diode. In our case: where V is the amplitude of the local oscillator signal (LO), which in this case is the radar’s microwave source operating at f0, V is the amplitude of the reflected signal (at f), t is time, and θ is the phase difference between f0 and f. The effect of mixing in a diode is more readily seen using a Taylor series expansion of the exponential term in the diode equation: Limiting this to first four (i = 0, 1, 2, 3) terms to facilitate discussion, we get: where . It will be apparent shortly why the series was truncated. By substituting the right side of Eq. (5) for V in Eq. (7) and expanding and combining terms, it can be shown that mixing in the diode results in frequency components at f0, f, 2f0, 2f, (f0 + f), (f0 − f), 3f0, 3f, (2f0 + f), (2f + f0), (2f0 − f), (2f − f0), etc. The output of the mixer is filtered to eliminate all but the first-order difference frequency, (f0 − f), which is the Δf computed in (3). This Δf is in the audio frequency range and can range from several hundred hertz to 20 kilohertz, depending on the values of v and f0. Present DTR radar technology and design uses a digital signal processing (DSP) unit to analyze the audio frequency signal that is output from the filter. The accuracy of the estimate of Δf and, consequently, the accuracy of the estimate of target vehicle speed, is dependent on the stability and known value of the DSP clock. What is important to note here is that the computed value of Δf is dependent on the frequency, f0, of the microwave source, the speed of the vehicle, v, and the DSP clock frequency. There are also secondary effects that are dependent on v, the bandwidth of the audio filter, and the relative amplitude of the second-order and fourth-order difference frequencies.

3. Calibration Methods

There are four primary tools that are used in DTR radar calibration methods: a tuning fork, a speed simulator based on amplitude modulation (AM) of the reflected radar signal, the vehicle’s speedometer, and a fifth-wheel. A fifth wheel is a thin wheel and tire, similar in appearance to a unicycle without a seat, which is attached to the rear of the vehicle in which the radar is mounted (for a description of its usage, see Ref. [2]). Before discussing each of these tools and associated calibration methods, it is important to note the definition of a few terms typically used to describe measurement results. The first term is measurement error. Error is simply the difference between the measured value and the actual value of the measurand, which in this case is vehicle speed. An error is not an uncertainty. The second term is measurement accuracy, which describes how close a measured value is to the actual value. Measurement error affects measurement accuracy; however, a known measurement error can be removed from a measured value to improve measurement accuracy. The last term is measurement uncertainty. Measurement uncertainty characterizes the dispersion of the values that could reasonably be attributed to the measurand. Uncertainty includes the limitations of measurement instrumentation, environmental effects, measurement standards, etc. Errors are not part of an uncertainty. However, if an error is removed from a measurement, the uncertainty in the value of the error must be considered in establishing the uncertainty of a measurand. In the following sections, approximate values for uncertainties will be given to provide the reader an opportunity to compare different calibration methods. So that this comparison is uniform, we have arbitrarily chosen the target speed to be 96.6 km/h (60 mph). Furthermore, the uncertainties provided in all but Sec. 4 and as noted, are at the one standard deviation (1σ) level and represent the range in v corresponding to a confidence level of about 68.3 %. This means, for example, if v = 96.6 km/h (60 mph) and the uncertainty is 2 km/h (1.24 mph), we can say with 68.3 % confidence that the vehicle was traveling between 94.6 km/h (58.8 mph) and 98.6 km/h (61.3 mph).

3.1 Tuning Fork

The tuning forks used to calibrate DTR radar are usually machined from a solid aluminum blank to provide an acoustic resonant frequency that gives a specific speed indication for a given f0. To calibrate the DTR radar, the tuning fork is struck and then placed in front of the transceiver. The vibrations of the tuning fork modulate the phase of a radar signal reflected from it. This modulation simulates the difference frequency that is determined by the radar processing unit from a Doppler shifted return from a moving vehicle. However, the tuning fork does not actually simulate a moving vehicle. This difference will be clarified with the following discussion. The tuning fork imposes a phase modulation on the incident microwave signal. The frequency, f1, of the modulation is described by: where A is the amplitude of the modulation, f is the designed resonant frequency of the tuning fork, and t is time. The signal reflected from the tuning fork is detected in the diode mixer and subsequently filtered to yield only the difference frequency, similar to what happens for the signal returned from the moving target vehicle (see Sec. 2). For the tuning fork modulation, this difference frequency, Δf, is given by: Note, this frequency is sinusoidally varying with time so that the actual signal, V, out of the filter is described by: To understand how this modulation results in a signal inside the DTR radar from which the modulation frequency (and hence the speed simulated by the tuning fork) is determined, V, must be expanded. The expansion of V is: where the J are Bessel functions of the first kind. This equation shows that the frequency components will be at the odd-numbered frequency harmonics of f. Since DTR radar units are capable of measuring speeds of up to 322 km/h (200 mph), and we may want to measure speeds down to 16 km/h (10 mph), we maintain the first few Bessel functions for this discussion, which will now be determined. The Bessel functions of the first kind can be further expanded using: Because of the two factorials in the denominator and because |x| ≤ 1, only J1 needs to be considered (J5 ≈ 1/50 J1 or less and J3 ≈ 1/10 J1 or less). Therefore, the dominant signal component for which the DSP will compute the simulated Doppler frequency shift is at f. The tuning fork is machined so that f is equal to Δf, for a particular pair of v and f0. The value of f is dependent on the material of construction and the length and cross sectional area of the tines. Consequently, a given tuning fork should only be used to calibrate a DTR radar operating at a specified frequency and for a given speed. The main advantages of the tuning fork calibrations are that the tuning forks are inexpensive to make, are small and portable, and require no power. They are also easy to use in the field, and this allows the officer to test the operation of the DTR radar often. However, because the material most often used is aluminum, f will depend on temperature, which is now discussed. An extensive test over the entire range of tuning fork operating temperatures was not performed, but f for several tuning forks was measured over a smaller range of temperatures to demonstrate the relationship between temperature and f. The tuning forks were placed in a temperature-controlled chamber for several hours before each measurement. Typically tuning forks are stored in the interior of a vehicle and are subject to the ambient temperature of the vehicle and possibly direct sunlight. The tuning forks were all actual forks used to calibrate DTR radar. Figure 2 shows the results of the measurements for one particular tuning fork (referred to here as TF1) designed to simulate a vehicle moving at 56.3 km/h (35 mph) for a K-band DTR radar.
Fig. 2

Plot of frequency (f) vs. temperature for a tuning fork designed for 56.3 km/h (35 mph) at K-band. The Type A standard uncertainty is approximately 0.025 Hz.

A linear least squares approximation was used to obtain a formula to fit the data shown in Fig. 2 for TF1 as well as data (not shown) for the other tuning forks. The result is the following formula, where S is the slope, in Hz/°C, of the fit to the f vs. T data and f,0 is the intercept of the fit at T = 0 °C. The values of S and f,0 are unique for each tuning fork and, for the tuning forks examined here, these values are given in Table 1. Equation (13) can be used to approximate the variation in speed measurement for a given tuning fork. For example, for TF1 and a K-band radar unit transmitting at 24.05 GHz, a temperature variation from −12.2 °C to 71.1 °C (10 °F to 160 °F) would cause a variation in Δf that corresponds to a target vehicle speed ranging from 55.8 km/h to 57.2 km/h (34.7 mph to 35.5 mph).
Table 1

Parameters for the Tuning Forks used in this Study

Tuning forkFrequency bandDesign speedkm/h (mph)STF(Hz/°C)fTF,0(Hz)
TF1K56.3 (35)−0.6882535.8
TF2  Ka32.2 (20)−0.5522028.3
TF3X30 (18.6)−0.152  589.4
It should be noted that f is not dependent on f0, and so the signal analyzed by the DSP is not dependent on f0. Consequently, the computed tuning-fork speed, v, of the tuning-fork simulated target vehicle speed is: which means that the tuning fork must be used at the radar frequency, f0, for which it was designed otherwise the value of v will be incorrect. Similarly, drift in f0 will cause errors in v. To minimize the uncertainty in the value of f and subsequently v, f should be measured using a frequency counter or spectrum analyzer. In this case the standard uncertainty in f ≈ 10−5 f.

3.1.1 Uncertainty in Tuning-Fork-Computed Target Vehicle Speed,

The standard uncertainty in the speed computed using a tuning fork, , is computed and presented using standard methods [2,3] to be: where is the standard uncertainty in f0, which is about 10−5 f0 for typical microwave sources under normal operating conditions, and is the standard uncertainty in f, which can be expanded using (13) to give: where is the standard deviation in the measurement of f, which is about 10−5 f, is the variation in the value of in the fit to the f vs. T data, σS is the variation (residuals) in the fit to the f vs. T data, and σ is the uncertainty in T (0.5 °C for this measurement). For typical values of the contributing uncertainty parameters, or about 0.3 km/h (0.2 mph) for v = 96.6 km/h (60 mph).

3.2 Speed Simulator Based on Amplitude Modulation of the Reflected Signal

DTR radars may be calibrated in a laboratory using a moving target simulator. The particular simulator discussed here is a small benchtop anechoic chamber comprising a radar transceiver that receives the radar signal (at f0) and subsequently retransmits this signal to the radar under test after this signal has been amplitude modulated. To the best of our knowledge, this is the only simulator design currently used. An audio function generator provides the modulation via a PiN voltagetunable diode attenuator. As with the tuning fork, we will provide an expression describing the effect of this modulation on the signal. The amplitude modulated signal is: where A is the modulation amplitude, f is the modulation frequency, and φ is the modulation phase. Expanding this gives: Equation (18) describes the signal that will be mixed in the diode with local oscillator (LO) signal from the radar. As we did for the tuning fork (which is a case describing phase modulation), we use the approximation to the diode equation response to find that the frequencies of the signal output from the diode mixer and filter are at ± f. Therefore, to calibrate a DTR radar, f has to equal the Δf. This calibration process is more versatile than using a tuning fork because f can be adjusted to almost any desired frequency. It should be noted that f is not dependent on f0, and so the signal analyzed by the DSP is not dependent on f0. The speed estimate, v, based on the simulator, is given by:

3.2.1 Uncertainty in Simulator-Computed Target Vehicle Speed,

The uncertainty in the speed computed using a amplitude modulating moving target simulator, , is computed from (19) to be: where is the uncertainty in the frequency of the modulation, which is equal to the clock frequency uncertainty of the source, so that . Consequently or about 0.0014 km/h (0.0008 mph) for v = 96.6 km/h (60 mph).

3.3 Vehicle Speedometer

The speedometer is a device that counts the number of times an object of a given radius rotates for a given time. Since the late 1980s, the most common implementation in passenger vehicles is an electronic speedometer [4]. In this design, an electronic sensor measures the number of rotations, N, of the transmission’s output shaft during an interval of time, t. Accordingly, the speed obtained from the speedometer, v, is computed from a measure of the number, N, of rotations; the effective radius of the tire, r; the known differential gear ratio, g; and the time, t. N and t can be determined relatively more accurately than r because the object is typically a pneumatic tire, and the effective radius can change due to the variety of influences that follow.

3.3.1 Uncertainty in Speedometer-Computed Target Vehicle Speed,

The uncertainty in the speed computed using the vehicle’s speedometer, , using (21), is given by: where is the uncertainty in the determination of the time period required for N revolutions, is the uncertainty in the number of revolutions counted, , is the uncertainty in the differential’s input/output gear ratio, and is the uncertainty in the effective radius of the tire. Typically, is fixed by the uncertainties in the clock frequency, so we can approximate by: . The uncertainty in N will be less than N will be less than 1 revolution, and typical values for The value of is, based on typical manufacture, to be about 0.005 , and in typical automotive applications 2.5 ≤ ≤ 3.5. An approximate value for is 1.65 cm, and the Appendix contains a discussion and derivation of this value. Using these values for the uncertainties and g = 3, t = 0.5 s, and r = 0.326 m (typical values) for this example gives or 4.9 km/h (3.1 mph) for v = 96.6 km/h (60 mph).

3.4 Fifth Wheel

As previously mentioned, speedometers determine speed by directly measuring the rotation of the transmission output shaft, and then applying a conversion factor to give the vehicle’s speed. The conversion factor manifests itself differently in mechanical and electronic speedometers. In mechanical speedometers it can be adjusted by changing the gearing on the speedometer cable, the strength of the permanent magnet or the stiffness of the restraining spring. In electronic speedometers, the conversion is handled by a DSP which converts pulses from the detector to a signal for display. In either case, accurately setting the conversion factor requires knowledge of the vehicle speed. The fifth wheel is towed behind the vehicle and measures the vehicle’s speed, v5, which may be used to set the vehicle speedometer’s conversion factor The fifth wheel is itself a calibrated tool. Calibration is performed by timing several passes at constant speed in a straight line across a calibrated distance, typically 0.5 miles. The ASTM calibration document describes the calibration technique in detail [1].

3.4.1. Uncertainty in Speed Obtained From the Fifth Wheel,

Fifth wheels are immune to many of the factors that affect speedometer uncertainty on passenger vehicles. Fifth wheels are always free rolling, operate within a narrow range of normal load (vertical load), have tires which vary little with speed and are not intended for severe cornering. This reduces the number of previously mentioned parameters that can introduce errors in a speedometer measurement. Nevertheless some of the tire effects that are described in the previous section should be accounted for in the fifth wheel. The uncertainty in v5 can be approximated by: where is the uncertainty in the fifth wheel calibration process, which is about and the remaining uncertainties are same as those described in Sec. 3.3, but as applied to the fifth wheel. The uncertainty in r for the fifth wheel is similar to that shown for the speedometer (see Eq. A8) but with some terms missing, as mentioned above. This uncertainty, is given by: Using the typical values for uncertainty contributions for a fifth wheel, we approximate 1.1 km/h (0.68 mph) for v = 96.6 km/h (60 mph).

4. Uncertainty in Vehicle Speed, u

The uncertainty in v as determined by a DTR radar will be primarily affected by the following: calibration method, f0, and clock frequency, as mentioned previously. Although other parameters and effects may contribute to speed measurement uncertainty, those uncertainty contributions are not observable given the limited resolution of DTR radar measurements. To compute u, we need v, which we get by rearranging (3): The uncertainty can be written as: where u is the uncertainty associated with one of the methods used to calibrate the DTR radar, that is, , , , or . The uncertainty in the DSP clock frequency, f, will affect the uncertainty in Δf. The value of Δf in current DTR radars is computed from the periodic sampling of the difference frequency output from the diode mixer. The effect of f on Δf can be effectively understood by using: where r is a constant that relates f0 and f for a given DTR radar. Therefore, where and for typical clocks, with f ≈ 10 MHz, and Δf ≤ 20 kHz, r ≤ 2 × 10−3. To be conservative, we set the uncertainty in r, u, to be which is about 10−5r. Putting all this together gives a value of u of about 0.3 Hz. Table 2 contains uncertainty values of parameters that contribute to speed measurement uncertainty. These are standard uncertainty values, computed for a target vehicle speed of 96.6 km/h (60 mph) and a confidence interval of 68.3 % (1σ). Only the calibration uncertainty corresponding to the calibration method used to calibrate the DTR radar should be used to compute speed measurement uncertainty.
Table 2

Uncertainty components and their values for a 68.3 % confidence interval (1σ) for a target vehicle speed of 96.6 km/h (60 mph)

ComponentValue
uvtf 3.1×103vtf
uvAM 1.4×105vAM
uvsm 2.587×103vsm2+1.165×103
uv5th 1.236×104v5th2+4.075×103m2s2
uΔf0.3 Hz
uf010–5f0
Table 3 lists the uncertainty in speed measurement (see Eq. 27) at several confidence intervals for each calibration method. Figure 3 shows a plot of these uncertainties for a 99.73 % confidence level (3σ) as a function of speed.
Table 3

Uncertainty in v for the various calibration methods and for various confidence intervals for a target vehicle speed of 96.6 km/h (60 mph)

Confidence interval
1σ(68.3 %)km/h[mph]2σ(95.5 %)km/h[mph]3σ(99.73 %)km/h[mph]4σ(99.99 %)km/h[mph]5σ(99.99994 %)km/h[mph]
Calibration method
Vehicle Speedometer4.9[3.1]9.8[6.1]15[9.2]20[12]25[15]
Fifth Wheel1.1[0.68]2.2[1.4]3.3[2.0]4.4[2.7]5.5[3.4]
Tuning Fork0.30[0.19]0.60[0.37]0.90[0.56]1.2[0.74]1.5[0.93]
Simulator2.2 × 10−03[1.3 × 10−03]4.3 × 10−03[2.7 × 10−03]6.5 × 10−03[4.0 × 10−03]8.6 × 10−03[5.3 × 10−03]1.1 × 10−02[6.7 × 10−03]
Fig. 3

Uncertainty for each calibration technique for a 99.73 % confidence interval (3σ).

5. Conclusions

We have examined the uncertainty in vehicle speed measurements provided by the ubiquitous down-the-road (DTR) radar used by law enforcement. We included in this analysis the effects of the calibration (field performance verification) method, which uses one of the following tools: tuning fork, moving vehicle simulator, vehicle speedometer, or fifth wheel. This analysis does not consider uncertainties in speed measurement that are caused by operational issues with DTR radar usage. The uncertainties in speed measurement due to calibration vary significantly depending on the method used to calibrate the DTR radar. The uncertainty in vehicle speed is the greatest using the vehicle’s speedometer as the calibration reference and the least when using the laboratory simulator. Furthermore, analysis provides information on the confidence to which a speed measurement can be given. The resulting equations from the uncertainty analysis are tabulated in Table 2. Table 3 and Fig. 3 are examples of how the resulting equations may be used to interpret the uncertainty for specific confidence intervals and speed values.
  1 in total

1.  A Vehicular Mobile Standard Instrument for Field Verification of Traffic Speed Meters Based on Dual-Antenna Doppler Radar Sensor.

Authors:  Lei Du; Qiao Sun; Changqing Cai; Jie Bai; Zhe Fan; Yue Zhang
Journal:  Sensors (Basel)       Date:  2018-04-05       Impact factor: 3.576

  1 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.