The chirp pulse compression process transforms a long duration frequency-coded pulse into a narrow pulse of greatly increased amplitude. It is a technique used in radar and sonar systems because it is a method whereby a narrow pulse with high peak power can be derived from a long duration pulse with low peak power. Furthermore, the process offers good range resolution because the half-power beam width of the compressed pulse is consistent with the system bandwidth.

The basics of the method for radar applications were developed in the late 1940s and early 1950s,^[1][2][3] but it was not until 1960, following declassification of the subject matter, that a detailed article on the topic appeared the public domain.^[4] Thereafter, the number of published articles grew quickly, as demonstrated by the comprehensive selection of papers to be found in a compilation by Barton.^[5]

Briefly, the basic pulse compression properties can be related as follows. For a chirp waveform that sweeps over a frequency range F1 to F2 in a time period T, the nominal bandwidth of the pulse is B, where B = F2 – F1, and the pulse has a time-bandwidth product of T×B. Following pulse compression, a narrow pulse of duration τ is obtained, where τ ≈ 1/B, together with a peak voltage amplification of √T×B.

The chirp compression process – outline

In order to compress a chirp pulse of duration T seconds, which sweeps linearly in frequency from F1 Hz to F2 Hz, a device with the characteristics of a dispersive delay line is required. This provides most delay for the frequency F1, the first to be generated, but with a delay which reduces linearly with frequency, to be T second less at the end frequency F2. Such a delay characteristic ensures that all frequency components of the chirp pass through the device, to arrive at the detector at the same time instant and so augment one another, to produce a narrow high amplitude pulse, as shown in the figure:

An expression describing the required delay characteristic is

${\displaystyle Y(f)=\exp[j\pi (f-f_{0})^{2}/k]}$

This has a phase component ψ(f), where

${\displaystyle \psi (f)=\pi .(f-f_{0})^{2}/k}$, and the instantaneous delay t_d is given by

${\displaystyle t_{d}=-{\frac {1}{2\pi }}.{\frac {d\psi }{df}}={\frac {1}{k}}.(f_{0}-f)}$

which has a linear slope with frequency, as required. In this expression, the delay characteristic has been normalized (for convenience), so as to give zero delay when the frequency f equals the carrier frequency f₀. Consequently, when the instantaneous frequency is (f₀ − B/2) or (f₀ + B/2), the required delay is +T/2 or −T/2, respectively, so k = B/T.

The required dispersive characteristic may be obtained from a lumped element delay network,^[6][7][8][9] a SAW device,^{[10][11][12][13][14]} or by means of digital signal processing ^[15][16][17]

An overview of pulse compression concepts

A general closed form solution for chirp compression

The characteristic of a single linear chirp pulse, of unit amplitude, can be described by

${\displaystyle s_{1}(t)=\operatorname {rect} \left({\frac {t}{T}}\right).e^{2\pi j(f_{0}t+kt^{2}/2)}}$

where rect(z) is defined by rect(z) = 1 if |z| < 1/2 and rect(z) = 0 if |z| > 1/2

The phase response φ(t) is given by

${\displaystyle \phi (t)=2\pi (f_{0}t+kt^{2}/2)}$

and the instantaneous frequency f_I is

${\displaystyle f_{i}={\frac {1}{2\pi }}.{\frac {d\phi }{dt}}=f_{0}+k.t}$

So, during the T second duration of the pulse, the frequency changes in a linear fashion from f₀ – kT/2 to f₀ + kT/2. With the net frequency sweep defined as B, where B = (F1- F2), then k = B/T, as stated previously.

The spectrum of this waveform may be found from its transform

${\displaystyle S_{1}(f)=\int _{-\infty }^{\infty }s_{1}(t).e^{-j2\pi .ft}.dt=\int _{-{\frac {T}{2}}}^{\frac {T}{2}}e^{j2\pi .[(f_{0}-f)t+{\frac {kt^{2}}{2}}]}.dt}$

which is an integral which has been evaluated in chirp spectrum.

The spectrum of the compressed pulse can be found from

${\displaystyle S_{\text{out}}(f)=Y(f).S_{1}(f)}$

Where Y(f) is the spectrum of the compression filter.

The time domain waveform ${\displaystyle s_{\text{out}}(t)}$ of the compressed pulse can be found as the inverse transform of ${\displaystyle S_{\text{out}}(f)}$. (This procedure has been described in a paper by Chin and Cook.^[9][43])

Here it is more convenient to find ${\displaystyle s_{\text{out}}(t)}$ from the convolution of the two time domain responses, i.e.

${\displaystyle s_{\text{out}}(t)=s_{1}(t)^{*}y(t)}$

where the convolution of two arbitrary functions is defined by

${\displaystyle f(t)^{*}g(t)=\int _{-\infty }^{\infty }f(\tau ).g(t-\tau ).d\tau }$

However, in order to use this method, the impulse response of Y(f) is first needed. This is y(t), which is obtained from

${\displaystyle y(t)=\int _{-\infty }^{\infty }Y(f).e^{j2\pi .ft}.df=\int _{-\infty }^{\infty }e^{j\pi .(f-f_{0})^{2}/k}.e^{j2\pi .ft}.df}$

${\displaystyle {{\text{Put}}\quad f-f_{0}=u\quad {\text{so}}\quad f=u-f_{0}\quad {\text{and}}\quad df=du.\qquad {\text{Also when}}\quad f=\pm \infty \quad {\text{then}}\quad u=\pm \infty }}$

${\displaystyle y(t)=\int _{-\infty }^{\infty }e^{j\pi .u^{2}/k}.e^{j2\pi .(u+f_{0})t}.du=e^{j2\pi .f_{0}t}\int _{-\infty }^{\infty }e^{j\pi .u^{2}/k}.e^{j2\pi .ut}.du}$

A table of standard integrals^[44] gives the following transform ${\displaystyle \int _{-\infty }^{\infty }\exp(-\pi \beta .u^{2}).\exp(j2\pi .ut).du={\frac {1}{\sqrt {\beta }}}.\exp \left(-{\frac {\pi .t^{2}}{\beta }}\right)}$

Comparing the equations, they are equivalent if β = -j/k, so y(t) becomes

${\displaystyle y(t)={\sqrt {jk}}.e^{j2\pi .f_{0}t}.e^{-j\pi .t^{2}k}={\sqrt {\frac {jB}{T}}}.e^{j2\pi .f_{0}t}.e^{-j\pi .t^{2}.k}={\sqrt {\frac {jB}{T}}}.e^{j2\pi (f_{0}t-t^{2}k/2)}}$

[Note: the same transform can also be found in Fourier transforms, no. 206, but with α replacing πβ]

With y(t) determined, the output s_out(t) can be obtained from the convolution of s₁(t) and y(t), i.e.

${\displaystyle s_{\text{out}}(t)={\sqrt {\frac {jB}{T}}}.\int _{-T/2}^{T/2}e^{j2\pi (f_{0}\tau +k\tau ^{2}/2)}\cdot e^{j2\pi (f_{0}(t-\tau )-k(t-\tau )^{2}/2)}.d\tau }$

which can be simplified to

${\displaystyle s_{\text{out}}(t)={\sqrt {\frac {jB}{T}}}.e^{j2\pi (f_{0}t-kt^{2}/2)}.\int _{-T/2}^{T/2}e^{j2\pi kt\tau }.d\tau }$

now as ${\displaystyle \int _{}^{}e^{ax}.dx={\frac {e^{ax}}{a}}}$ then

${\displaystyle \int _{-T/2}^{T/2}e^{j2\pi .kt\tau }.d\tau ={\frac {1}{j2\pi kt}}.{\Bigl [}e^{j2\pi .kt\tau }{\Bigr ]}_{\text{-T/2}}^{\text{T/2}}={\frac {1}{j2\pi kt}}[e^{j2\pi kT/2}-e^{-j2\pi kT/2}]}$

${\displaystyle ={\frac {1}{\pi kt}}.{\frac {e^{j\pi Bt}-e^{-j\pi Bt}}{2j}}={\frac {1}{\pi kt}}.\sin(\pi Bt)}$

and finally

${\displaystyle s_{\text{out}}(t)={\frac {1}{\pi kt}}\cdot {\sqrt {\frac {jB}{T}}}\cdot \sin(\pi Bt).\exp[j2\pi (f_{0}t-kt^{2}/2)]}$

${\displaystyle ={\sqrt {T.B}}\cdot {\frac {\sin(\pi Bt)}{\pi Bt}}\cdot {\sqrt {j}}\cdot \exp[j2\pi (f_{0}t-kt^{2}/2)]}$

So, for a unit amplitude linear chirp, with a pulse duration T seconds and a frequency sweep B Hz (i.e. with a ‘time-bandwidth product’ T.B), pulse compression gives a waveform with a magnitude given by

${\displaystyle \left|{\text{Compressed Output}}\right|\approx {\sqrt {T.B}}\times {\frac {\sin(\pi Bt)}{(\pi Bt)}}}$

which has the form of the familiar sinc function. The collapsed pulse width τ, is of the order 1/B (with τ measured at the −4 dB points). Consequently, a pulse width reduction has occurred given by the ratio T/τ where ${\displaystyle {\frac {T}{\tau }}\approx T\times B}$

Also there is a signal amplification of ${\displaystyle {\sqrt {T\times B}}}$

The main parameters are displayed in the figures below. The product T.B gives the compression ratio of the system and it equates, approximately, to the improvement in signal to noise ratio of the main lobe of the compressed pulse relative to that of the original chirp.

Properties of linear chirps

Pulse degradation caused by Fresnel ripples

In the closed form solution, just presented, the compressed waveform has the standard sinc function response, because a rectangular shape was assumed for the amplitude of the pulse spectrum. In practice, the spectrum of a linear chirp has a rectangular profile only when the time-bandwidth product of a pulse is large, i.e. when T×B exceeds 100, say. When the product is small, the spectral profile of the chirp pulse is seriously degraded by Fresnel ripples as shown in chirp spectrum and so is that of the matched filter. To investigate fully the consequences of these ripples it is advisable to consider each situation individually, either by evaluating convolution integrals, or more conveniently, by means of FFTs.

Some examples are shown below, for TB = 1000, 250, 100 and 25. They are dB plots which have all been normalised to have their pulse peaks set at 0 dB.

As can be seen, at high values of TB, the plots match the sinc characteristic closely, but at low values, significant differences can be seen. As already stated, these degradations in the waveforms, at low TB values, are because the spectral characteristics are no longer truly rectangular. In all cases, the close-in sidelobe levels are consistently high, at about −13.5 dB relative to the main lobe.

These range sidelobes are an unwelcome presence on the compressed pulse, because they will obscure signals of lower amplitude, which may also be present.

Reducing sidelobes by weighting functions

As the sinc-like characteristics of a compressed pulse are due to the near-rectangular profile of its spectrum, then by modifying that characteristic into a bell shape, for example, it is possible to reduce the sidelobe levels considerably. Previous work for antenna arrays and digital signal processing has already addressed this same problem. So, for example, in the case of antennas, the spatial sidelobes on the beam pattern are improved by applying a weighting function to the array elements,^[45] and in the case of digital signal processing, window functions are used to reduce the amplitude of unwanted sidelobes^[18] on the sampled functions.

In an example of the process, the spectrum of a chirp pulse, with a time-bandwidth product of 250, is shown on the left, and is approximately rectangular in profile. Below this plot, also on the left, is shown the waveform after compression of the chirp by its matched filter and is similar to the sinc function, as expected. The top plot, on the right, is the spectrum after Hamming weighting. (This was achieved by applying a root-Hamming characteristic to both the chirp spectrum and to the compressor spectrum.) The compressed pulse corresponding to this spectrum, shown in the lower plots on the right, has much lower sidelobe levels.

Although the sidelobe level has been much reduced, there are some undesirable consequences of the weighting process. Firstly, there is an overall loss of gain, with the peak amplitude of the main lobe reduced by about 5.4 dB and, secondly, the half-power beam width of the main lobe has increased by nearly 50%. In, say, a radar system, these effects would cause loss of range and reduced range resolution, respectively.

In general, the more the sidelobe levels are lowered, the broader the main lobe will become. However, the various windowing functions do perform differently from one other, with some giving main lobes which are unnecessarily broad for the sidelobe levels achieved. The most efficient function is the Dolph–Chebyshev window (see window functions) as this gives the narrowest pulse at a given sidelobe level.^[18] A selection of the better performing windowing functions are shown in the graph of Beamwidth × Bandwidth as sidelobe level.

The lowest full line on the graph is for Dolph-Chebyshev weighting which, as already mentioned, sets the narrowest lobe possible for a given sidelobe level. So, from this plot, if a sidelobe level of −40 dB is desired, the graph shows that the smallest achievable half-power beamwidth × bandwidth is 1.2. Thus a chirp sweeping over a 20 MHz frequency band will have a compressed pulsewidth of 60 nanoseconds (at least).

As can be seen from the diagram, Taylor weighting is particularly efficient, and Hamming and the three and four term Blackman-Harris functions all give good results, too. Although the cos^N functions perform badly, they have been included because they are amenable to mathematical manipulation and were studied in some detail in early work.^[23][46]

Far-out sidelobes on compressed pulses

The example of a chirp with TB = 250 and Hamming weighting, given earlier, illustrates the benefits of weighting but is not representative of a normal situation because the results there were achieved by applying weighting equally to both the signal chirp and to its compressor. However, in a typical radar system the chirp pulse is usually transmitted by an amplifier operating in or near to compression, to maximise transmitter efficiency. In such a case, amplitude modulation of the chirp waveform or its spectrum is not possible, so the full window characteristic has to be incorporated into the compressor response. Unfortunately, this arrangement has undesirable consequences for the far-out sidelobes of the compressed pulse, especially when the time-bandwidth of the chirp is small.

Consider first the compressed pulse when TB = 250, which is shown in the left-hand figure below. For this result no weighting has been applied to the transmit pulse, but the full Hamming weighting has been applied in the compressor. As can be seen, the close-in sidelobe levels are consistent with the Hamming weighting (−42 dB), but further out, the sidelobe levels rise to a peak value of −45 dB at +/-T/2 each side of the main lobe. In the right-hand figure, where TB =25, the problems with the far-out sidelobes are much more serious. Here these sidelobes now rise to −25 dB at +/-T/2.

As a guide, the far-out sidelobe levels are given by

${\displaystyle FarSL(dB)\approx -20\times \log _{\text{10}}(T.B)}$

Slight variations on this equation are given in the literature,^[47][48][49] but they only differ by a few dBs. Best results seem to be obtained when the window function is applied in the time domain to the compressor waveform (as amplitude modulation) rather than in the frequency domain to its spectrum.^[50]

Reducing far-out sidelobes

As the far-out sidelobes are a consequence of the Fresnel ripple on the compressed pulse spectrum, any technique that reduces this ripple will also reduce the level of the sidelobes. In fact, there are several ways of achieving this reduction,^[51] as shown below. Several of the methods are featured in chirp spectrum.

Introducing finite duration rise and fall times

A chirp with slow rise and fall times has reduced ripple on its spectrum (see chirp spectrum), so will result in lower time sidelobes on the compressed pulse. As an example, consider first the figure shows the compressed spectrum of a linear chirp, which has fast rise and fall times, with T×B = 100 and where Blackman-Harris weighting has been applied. The waveform corresponding to this spectrum has time sidelobes rising to about −40 dB, as predicted.

After linear rise and fall times are introduced, using the amplitude template shown, the ripple on the spectrum is much reduced and the time sidelobes are significantly lower, as shown.

The procedure is most effective when both the signal chirp and the compressor chirp have modified rise-times, when sidelobe levels can be lowered by 15 – 20 dB. However, it is not always possible to apply amplitude modulation in the transmitter, so there is less improvement when only the compressor waveform is modified. Even so, a sidelobe reduction of about 6 dB can still be achieved.

The precise manner in which the rise and fall times are made less severe is not very critical, so the technique of adding cosine tapers to the compressed pulse spectrum (as with the Tukey^[18] weighting function) gives a similar improvement – of several dBs.^[21]

Improvements achieved by the method are tolerant of Doppler shifts.

Introducing 'tweaking' of the phase characteristic

An alternative form of waveform 'tweaking' is where frequency modulation distortion is applied to the chirps, instead of amplitude modulation.^[23][52][53] The two types of distortion are functionally similar when the distortion levels are low. As with amplitude modulation, best results are obtained when both the expander and compressor waveforms are modified.

For best results Cook and Paolillo recommend δf = 0.75×B and δ = 1/B.

As an example, a pulse considered earlier, with T×B = 100 and Blackman–Harris weighting, is modified with phase tweaking and the results are shown. There is reduced ripple on the compressed pulse spectrum and the far-out sidelobes have been reduced.

The improvements are maintained even when Doppler frequency shifts are present on the signals. In a more recent paper ^[54] slightly different parameters have been suggested, namely δ = 0.86/B and δf = 0.73×B.

Also, Kowatsch and Stocker^[21] has reported improved results by applying a cubic distortion function (whereas Cook and Paolillo's technique may be termed ‘square-law modulation distortion’). This new characteristic is also tolerant of Doppler frequency shifts.

Reciprocal ripple correction

The spectral response of the matched filter has a magnitude which is the mirror image of that of the expanded pulse, when a chirp spectrum has symmetry about its centre frequency, so the Fresnell ripples on the spectrum are augmented by the compression process. What is needed to reduce the ripples is a compression filter whose spectrum has the inverse (reciprocal) ripple to that of the expander.^[23] As this will no longer be a matched filter, there will be increased mismatch loss. In his early work^[6][9][23] Cook did not recommend attempting such a procedure because the filters needed were considered too difficult to make. However, with the advent of SAW technology it became possible to achieve the required characteristics.^{[11][12][22][33]} More recently, digital techniques with mathematically derived look-up tables have provided a convenient way of introducing reciprocal ripple correction.^[16]

The spectrum of the compressed pulse is the product of the spectra of the expander and compressor filters, as given earlier. Now, instead of C(ω), a new output spectrum C'(ω) is defined, which has no Fresnell ripples but which defines a desired sidelobe structure (such as that defined by a Hamming window). The compression filter that will achieve this requirement is determined by the equation

${\displaystyle K(\omega )={\frac {C'(\omega )}{H(\omega )}}}$

where H(ω) is the spectrum of the signal, C’(ω) is the target spectrum for the compressed pulse and has the low sidelobes of the chosen weighting function and K(ω) is the spectrum of the compression filter which has the reciprocal ripple properties. Close-in sidelobes are automatically dealt with in the process.

As an example of the procedure consider a linear chirp with T×B =100. The left-hand figures shows (one half of) the spectrum of the chirp, and the right-hand figure shows the waveform after compression. As expected, close-in sidelobes start at −13.5 dB.

In the next figure, Blackman-Harris weighting has been applied to the compressed pulse spectrum. Although the close-in sidelobes have been reduced, the far-out sidelobes remain high with a predicted level of, approximately −20×log₁₀(100) = -40 dB, as predicted for a time-bandwidth product of 100. With lower time-bandwidth products, these sidelobes will be even higher.

Next, a compression filter that provides reciprocal-ripple correction has been used. As can be seen, a ripple-free spectrum has been achieved resulting in a waveform that is free from high level far-out sidelobes.

However, this procedure has a problem. Although the process has found a compressor spectrum that leads to low sidelobes on the compressed pulse, no account was taken of the waveform this spectrum might have. When an inverse Fourier transform is carried out on this spectrum, in order to determine the characteristics of its waveform, it is found that the waveform is of extremely long duration, typically exceeding 10T. Even assuming the waveform is no longer than 10T, it means that the total time needed to process one chirp will be at least 11T, in total, a length of time unacceptable in most circumstances.

In order to achieve a practical solution Judd^[22] proposed that the total length of the compression pulse be truncated to 2T, whereas Butler^[11] suggested 1.6T and 1.3T. Extensions as low as 10% have also been used^[55]

Unfortunately, when the new compressor waveform is truncated, then far-out sidelobes reappear once more. The next figures shows what happens to the compressed pulse when the compressor is set at 2T duration and then at 1.1T duration. New far-out sidelobes have appeared with amplitudes that make them clearly visible. These sidelobes are often referred to as "gating sidelobes".^[54] They can be irritatingly high but, fortunately, even if the compressor is set have just 10% extension, the sidelobes are still at a level than that achieved without correction.

Any Doppler frequency shift on the received signals will degrade the ripple cancelation process^[11][21] and this is discussed in more detail next.

Non-linear chirps

To ensure that a compressed pulse has low time sidelobes, its spectrum should be approximately bell-shaped. With linear chirp pulses this can be achieved by applying a window function either in the time domain or in the frequency domain, i.e. by amplitude modulating the chirp waveforms or by applying weighting to the compressed pulse spectra. In either case there is a mismatch loss of 1½dB, or more.

An alternative way to obtain the required spectral shape is to use a non-linear frequency sweep in the chirp. In this case, to achieve the required spectral shape, the frequency sweep changes very rapidly at band edges and more slowly around band centre. Consider, as an example, the frequency versus time plot that achieves the Blackman-Harris windowing profile. When T×B =100, the spectrum of the compressed pulse and the compressed waveform are as shown.

The required non-linear characteristic can be derived using the method of stationary phase.^[24][57] As this technique does not take account of the Fresnel ripples, these have to be dealt with in additional ways, as was the case with linear chirps.

In order to achieve the required spectral shape for low time sidelobes, linear chirps require amplitude weighting and consequently incur a mismatch loss. Non-linear chirps, however, have the advantage that by achieving the spectral shaping directly, close-in sidelobe levels can be made low with negligible mismatch loss (typically less than 0.1 dB). Another benefit is that the far out sidelobes, due to Fresnel ripples on the spectrum, tend to be lower than for a linear chirp with the same T×B product (4 to 5 dB lower with large T×B).

However, for chirps where the T×B product is low, the far-out sidelobe levels of the compressed pulse can still be disappointingly high, because of high amplitude Fresnel ripples on the spectrum. As with linear chirps, results can be improved by means of reciprocal ripple correction but, as previously, truncation of the compression waveform results in the appearance of gating sidelobes.

An example of reciprocal ripple and truncation is shown below. The left hand figure shows the spectrum of a non-linear chirp, with a time bandwidth product of 40, aiming to have a Blackman-Harris profile. The right-hand figure shows the compressed pulse for this spectrum,

The next figures shows the spectrum after r-r compensation, but with truncation of the compression waveform to 1.1T, and the final compressed waveform.

Doppler tolerance of non-linear chirps

A major disadvantage of non-linear chirps is their sensitivity to Doppler frequency shifts. Even modest values of Doppler will result in broadening of the main pulse, raising of the sidelobe levels, increase in mismatch loss and the appearance of new spurious sidelobes.

An example of a non-linear chirp pulse and the effects of Doppler are shown. The non-linear characteristic is chosen to achieve −50 dB sidelobes using Taylor weighting. The first figure shows the compressed pulse for a non-linear chirp, with bandwidth 10 MHz, pulse duration 10usec, so T×B = 100, and with no Doppler shift. The next two figures shows the pulse degradation cause by 10 kHz and 100 kHz Doppler, respectively. In addition to the waveform degradation, the mismatch loss increases to 0.5 dB. The final figure shows the effect of 100 kHz Doppler on a linear chirp which has had amplitude weighting applied to give the same spectral shape as that of the non-linear chirp. The greater tolerance to Doppler is clearly seen.

Cook,^[23] using paired-echo distortion methods,^[41] estimated that in order to keep sidelobe levels below −30 dB, the maximum allowed Doppler frequency is given by

${\displaystyle f_{d}\leq {\frac {0.06}{T}}}$

so, for a 10μs pulse, the maximum Doppler frequency that can be tolerated is 6 kHz. However, more recent work suggests that this is unduly pessimistic.^[33] In addition, as the new sidelobes, when at a low level, are very narrow. Consequently, it may be possible to ignore them initially, as they may not be resolvable by the receiver's D to A.

Using a combination of non-linear and linear characteristics to improve Doppler tolerance

A way of reducing the susceptibility of non-linear chirps to Doppler is to use a ‘hybrid’ scheme, where part of the spectral shaping is achieved by a non-linear sweep, but with additional spectral shaping achieved by amplitude weighting.^[11][12] Such a scheme will have greater mismatch loss than a true non-linear scheme, so the advantage of greater Doppler tolerance has to be weighed against the disadvantage of the increased mismatch loss.

In the two examples below, the chirps have a non-linear sweep characteristic which gives a spectrum with Taylor weighting which, used alone, will achieve a sidelobe level of −20 dB on its compressed pulses. To achieve lower level sidelobes, this spectral shape is augmented by amplitude weighting so that the final target sidelobe level for the compressed pulses is −50 dB. Comparing the results for Doppler shifts of 10 kHz and 100 kHz with those shown earlier it is seen that the new spurious sidelobes, caused by the Doppler, are seen to be 6 dB lower than before. However, the mismatch loss has increased from 0.1 dB to 0.6 dB, but this is still better than the 1.6 dB figure for linear chirps.

Signal-to-noise ratio improvements by pulse compression

The amplitude of random noise is not changed by the compression process, so the signal to noise ratios of received chirp signals are increased in the process. In the case of a high power search radars, this extends the range performance of the system, while for stealth systems the property will permit lower transmitter powers to be used.

As an illustration, a possible received noise sequence is shown, which contains a low amplitude chirp signal obscured within it. After processing by the compressor, the compressed pulse is clearly visible above the noise floor.

When pulse compression is carried out in digital signal processing, after the incoming signals are digitised by A/D converters, it is important that level of the noise floor is correctly set. The noise floor at the A/D must be high enough to ensure that the noise is adequately characterised. If the noise level is too low, Nyquist will not be satisfied, and any embedded chirp will not be recovered correctly. On the other hand, setting the noise level unnecessarily high will reduce the dynamic range capability of the system.

For systems using digital processing, it is important to carry out the chirp compression in the digital domain, after the A/D converters. If the compression process is carried out in the analogue domain before digitization (by a SAW filter, for example), the resulting high-amplitude pulses will place excessive demands on the dynamic range of the A/D converters.^[17]

Pre-correction of system characteristics

The transmitter and receiver subsystems of a radar are not distortion free. In consequence system performance is often less than optimum. In particular, the time sidelobe levels of the compressed pulses are found to be disappointingly high.

Some of the characteristics which degrade performance are:

• Gain slope, or non-linear phase slope, across the system passband.
• Amplitude and phase ripple across the passband (which may be caused by mismatches on interconnecting cables^[58] as well as by imperfections in amplifiers).

Delay modulation by the transmitter (if power supply regulation is poor).

In addition, filters employed in the frequency conversion processes of the transmitter and receiver all contribute to gain and phase variations across the system passband, especially near to band edges. In particular, major contributors to overall non-linear phase characteristics are the low-pass filters preceding the A/D converters, which are usually sharp-cut filters chosen to ensure maximum bandwidth while minimizing aliased noise. The transient response characteristics of these filters contribute another (unwanted) source of time sidelobes.

Fortunately it is possible to compensate for several system properties, provided they are stable and can be characterized adequately when a system is first assembled. This is not difficult to implement in radars using digital look-up tables, since these tables can be easily amended to include compensation data. Phase pre-corrections can be included in the expander tables and phase and amplitude corrections can be included in the compressor tables, as required.

So, for example, the earlier equation, defining the compressor characteristic to minimize spectral ripple, could be expanded to include additional terms to correct for known amplitude and phase impairments, thus:

${\displaystyle K(\omega )={\frac {C'(\omega )}{\Phi (\omega ).A(\omega ).H(\omega )}}}$

where, as before, H(ω) is the initial chirp spectrum and C'(ω) is the target spectrum, such as a Taylor window, but now additional terms have been included, namely, Φ((ω) ) and A(ω) which are the phase and amplitude characteristics that require compensation.

A compressor chirp waveform that includes phase correction data will have additional ripple components present at each end of the waveform (pre-shoots and after-shoots). Any truncation procedure should not remove these new features.

In addition, it is easy to time shift the compressed pulses by ±t₀, by multiplying the compressor spectrum by the unity amplitude vector, i.e.

${\displaystyle \cos(\omega .t_{0})\pm j.\sin(\omega .t_{0})=\exp(\pm j\omega .t_{0})}$.

Time shift can be useful to position the main lobes of compressed pulses at a standard location, regardless of chirp pulse length. However, care has to be taken with the overlap and save or overlap and discard algorithm, should time shift be used, to ensure only valid waveform sequences are retained.

There has been a growth in interest in adaptive filters for pulse compression, made possible by the availability of small fast computers, and some relevant articles are mentioned in the next section. These techniques will also compensate for hardware deficiencies, as part of their optimization procedure^[59]

More recent work on chirp compression techniques – some examples

The growth in digital processing and methods had a significant influence in the field of chirp pulse compression. An introduction to these techniques is provided in a chapter of the Radar Handbook (3rd ed.), edited by Skolnik.^[17]

The main aims of most investigations into pulse compression has been to obtain narrow main lobes, with low sidelobe levels, a tolerance to Doppler frequency shifts and to incur low system losses. The availability of computers has led to a growth in numerical processing and much interest in adaptive networks and optimization methods, to achieve these aims. For example, see the comparison of the various techniques made by Damtie and Lehtinen^[60] and, also, various articles by Blunt and Gerlach on these topics.^{[61][62][63][64]} A number of other contributors in the field has included Zrnic et al.^[65] Li et al.^[49] and Scholnik.^[59]

A number of other works, with a variety of approaches to pulse compression, are listed below:

• New methods of generating non linear chirp waveforms and of improving their Doppler tolerance has been investigated by Doerry^[66][67]
• Further studies of Hyperbolic chirps have been carried out by Kiss,^[68] Readhead,^[69] Nagajyothi and Rajarajeswari^[70] and Yang and Sarkar.^[71]
• Convolution windows have been investigated by Sahoo and Panda who show that they can result in very low sidelobes yet be Doppler tolerant, but may suffer from some pulse broadening.^[72] Wen and his co-workers have also discussed convolution windows.^[73][74]
• Some new window functions have been proposed by Samad^[75] and Sinha and Ferreira,^[76] which claim improved performance over the familiar functions.
• Several techniques to lower the sidelobe levels of the compressed pulses for non-linear FM chirps are compared by Varshney and Thomas.^[77]
• In a paper by Vizitui,^[78] sidelobe reduction is considered where phase pre-distortion is applied to non-linear FM chirps, rather than to linear chirps. Lower sidelobes and some improvement in Doppler tolerance is claimed.

There have been extensive investigations of phase modulation for pulse compression schemes, such as biphase (binary phase-shift keying) and polyphase coding methods, but this work is not considered here.

See also

• Spread spectrum

References