Chirp

From Wikipedia, the free encyclopedia
Jump to: navigation, search
For other uses, see Chirp (disambiguation).

A chirp is a signal in which the frequency increases ('up-chirp') or decreases ('down-chirp') with time. In some sources, the term chirp is used interchangeably with sweep signal.[1] It has also been called quadratic-phase signal.[2] It is commonly used in sonar and radar, but has other applications, such as in spread spectrum communications. In spread spectrum usage, SAW devices such as RACs are often used to generate and demodulate the chirped signals. In optics, ultrashort laser pulses also exhibit chirp, which, in optical transmission systems interacts with the dispersion properties of the materials, increasing or decreasing total pulse dispersion as the signal propagates. The name is a reference to chirping in analogy to the sound made by some birds, see bird vocalization.

Types of chirp[edit]

Linear chirp[edit]

A linear chirp waveform; a sinusoidal wave that increases in frequency linearly over time
Spectrogram of a Linear Chirp. The Spectrogram plot demonstrates the linear rate of change in frequency as a function of time, in this case from 0 to 7 kHz repeating every 2.3 seconds. The intensity of the plot is proportional to the energy content in the signal at the indicated frequency and time.
Sound example for linear chirp (5 repetitions).

Problems playing this file? See media help.

In a linear chirp, the instantaneous frequency f(t) varies linearly with time:

f(t) = f_0 + k t

where f_0 is the starting frequency (at time t = 0), and k is the rate of frequency increase or chirp rate.

k = \frac{f_1-f_0}{t_1}

where f_1 is the final frequency and  f_0 is the starting frequency.

The corresponding time-domain function for the phase of any oscillating signal is the integral of the frequency function, since one expects the phase to grow like \phi(t+\Delta t)\simeq\phi(t)+2\pi f(t)\,\Delta t, i.e., that the derivative of the phase is the angular frequency \phi'(t)=2\pi\,f(t).

For the linear chirp, this results in:

\begin{align}
\phi(t) &= \phi_0 + 2\pi\int_0^t f(\tau)\, d\tau\\
& = \phi_0 + 2\pi\int_0^t (f_0 + k \tau)\, d\tau\\
& = \phi_0 + 2\pi \left(f_0 t + \frac{k}{2} t^2 \right),
\end{align}

where \phi_0 is the initial phase (at time t = 0).

The corresponding time-domain function for a sinusoidal linear chirp is the sine of the phase in radians:

x(t) = \sin\left[\phi_0 + 2\pi \left(f_0 t + \frac{k}{2} t^2 \right) \right]

In the frequency domain, the instantaneous frequency described by the equation f(t) = f_0 + k t is accompanied by additional frequencies (harmonics) which exist as a fundamental consequence of frequency modulation.[citation needed] These harmonics are quantifiably described through the use of Bessel functions. However, with the aid of a frequency vs. time profile spectrogram one can readily see that the linear chirp has spectral components at harmonics of the fundamental chirp.

Exponential chirp[edit]

An exponential chirp waveform; a sinusoidal wave that increases in frequency exponentially over time
Spectrogram of an exponential chirp. The exponential rate of change of frequency is shown as a function of time, in this case from nearly 0 up to 8 kHz repeating every second. Also visible in this Spectrogram is a frequency fallback to 6 kHz after peaking, likely an artifact of the specific method employed to generate the waveform.
Sound example for exponential chirp (5 repetitions).

Problems playing this file? See media help.

In a geometric chirp, also called an exponential chirp, the frequency of the signal varies with a geometric relationship over time. In other words, if two points in the waveform are chosen, t_1 and t_2, and the time interval between them t_2 - t_1 is kept constant, the frequency ratio f(t_2)/f(t_1) will also be constant.

In an exponential chirp, the frequency of the signal varies exponentially as a function of time:

f(t) = f_0 k^t

where f_0 is the starting frequency (at t = 0), and k is the rate of exponential increase in frequency. Unlike the linear chirp, which has a constant chirp rate, an exponential chirp has an exponentially increasing chirp rate.

The corresponding time-domain function for the phase of an exponential chirp is the integral of the frequency:

\begin{align}
\phi(t)
& = \phi_0 + 2\pi \int_0^t f(\tau)\, d\tau \\
& = \phi_0 + 2\pi f_0 \int_0^t k^{\tau} d\tau \\
& = \phi_0 + 2\pi f_0 \left( \frac{k^t - 1}{\ln(k)} \right)
\end{align}

Where \phi_0 is the initial phase (at t = 0).

The corresponding time-domain function for a sinusoidal exponential chirp is the sine of the phase in radians:

x(t) = \sin\left[\phi_0 + 2\pi f_0 \left( \frac{k^t - 1}{\ln(k)} \right)\right]

As was the case for the Linear Chirp, the instantaneous frequency of the Exponential Chirp consists of the fundamental frequency f(t) = f_0 k^t accompanied by additional harmonics.[citation needed]

Generation of a chirp signal[edit]

A chirp signal can be generated with analog circuitry via a VCO, and a linearly or exponentially ramping control voltage. It can also be generated digitally by a DSP and DAC, using a Direct digital synthesizer (DDS) and by varying the step in the numerically controlled oscillator. It can also be generated by a YIG oscillator.

Relation to an impulse signal[edit]

Chirp and impulse signals and their (selected) spectral components.

A chirp signal shares the same spectral content with an impulse signal. However, unlike in the impulse signal, spectral components of the chirp signal have different phases.[3] Dispersion of a signal propagation medium may result in unintentional conversion of impulse signals into chirps. On the other hand, many practical applications, such as chirped pulse amplifiers or echolocation systems,[4] use chirp signals instead of impulses because of their inherently lower PAPR.

Uses and occurrences[edit]

Chirp modulation[edit]

Chirp modulation, or linear frequency modulation for digital communication was patented by Sidney Darlington in 1954 with significant later work performed by Winkler in 1962. This type of modulation employs sinusoidal waveforms whose instantaneous frequency increases or decreases linearly over time. These waveforms are commonly referred to as linear chirps or simply chirps.

Hence the rate at which their frequency changes is called the chirp rate. In binary chirp modulation, binary data is transmitted by mapping the bits into chirps of opposite chirp rates. For instance, over one bit period "1" is assigned a chirp with positive rate a and "0" a chirp with negative rate −a. Chirps have been heavily used in radar applications and as a result advanced sources for transmission and matched filters for reception of linear chirps are available.

(a) In image processing, direct periodicity seldom occurs, but, rather, periodicity-in-perspective is encountered. (b) Repeating structures like the alternating dark space inside the windows, and light space of the white concrete, "chirp" (increase in frequency) towards the right. (c) Thus the best fit chirp for image processing is often a projective chirp.

Chirplet transform[edit]

Main article: Chirplet transform

Another kind of chirp is the projective chirp, of the form:

g = f\left[\frac{a \cdot x + b}{c \cdot x + 1}\right],

having the three parameters a (scale), b (translation), and c (chirpiness). The projective chirp is ideally suited to image processing, and forms the basis for the projective chirplet transform.[5]

Key chirp[edit]

A change in frequency of Morse code from the desired frequency, due to poor stability in the RF Oscillator is known as chirp,[6] and in the RST code is given an appended letter 'C'.

The spectrum of a chirp pulse[edit]

A general expression for an oscillatory waveform, centered on frequency ω0 is

s(t)= a(t).exp[j(\omega_0.t + \theta(t))]

where a(t) and θ0(t) give the amplitude and phase variations of the waveform s, with time.
The frequency spectrum of this waveform is obtained by calculating the Fourier Transform of s(t), i.e.

S(\omega) =\int_{-\infty}^{\infty}s(t).exp(-j\omega t).dt =
\int_{-\infty}^{\infty}a(t).exp[(-j(\omega_0 t + \theta(t)].exp(-j\omega t).dt
so
S(\omega)=\int_{-\infty}^{\infty}a(t).exp[j \left \{(\omega_0 - \omega).t + \theta(t) \right \}].dt

In a few special cases, the integral can be solved to give an analytical expression, but often the characteristics of a(t) and θ(t) are such that the integral can only be evaluated by an approximation algorithm or by numerical integration.

The linear chirp[edit]

In the special case where a(t) is constrained to be a flat topped pulse with its instantaneous frequency varying as a linear function of time, then an analytical solution is possible.
For convenience, the pulse is considered to have unit amplitude and be of duration T, with the amplitude and phase defined over the time interval -T/2 to +T/2. The total frequency sweep is ΔF, varying in a linear manner from -ΔF/2 to +ΔF/2 in the defined time interval.
When the frequency is a linear function of time, the phase is a quadratic function, and s(t) can be written

s(t)= 1.exp[j(\omega_0.t + \frac{\Delta \Omega}{2.T}.t^2)]\qquad where \quad\Delta \Omega=2\pi.\Delta F \qquad and \quad\frac{-T}{2} < t < \frac{T}{2}

The spectrum of this linear FM signal is
S(\omega)=\int_{-T/2}^{T/2}exp \left [j(\omega_0.t+\frac{\Delta\Omega}{2.T}.t^2) \right ].exp(-j\omega.t).dt
=\int_{-T/2}^{T/2}exp \left[ j \left\{(\omega_0-\omega).t+\frac{\Delta\Omega}{2.T}.t^2 \right \} \right ].dt

By recourse to the Fresnel integrals C(X) and S(X)[7]:35[8]:300, defined by
 C(X)=\int_{0}^{X}cos \frac{\pi.y^2}{2}.dy \qquad and \qquad S(X)=\int_{0}^{X}sin \frac{\pi.y^2}{2}.dy

the expression can be evaluated[9][10][11][12]:138[13] to give

S(\omega)=  \sqrt{\left(\frac{T}{2.\Delta F}\right )}.exp \left[-j.(\frac{\pi.n^2}{4}.T.\Delta F)\right]
.[C(X_1)+j.S(X_1)+C(X_2)+j.S(X_2)]

where X_1 and X_2 are given by X_1=\frac{\sqrt{T.\Delta F}}{\sqrt{2}}.(1+n) and X_2=\frac{\sqrt{T.\Delta F}}{\sqrt{2}}.(1-n)

and where n is a frequency scaling term, given by n=2.\frac{(\omega-\omega_0)}{\Delta \Omega}

The spectrum amplitude of this chirp is
|S(\omega)|= \sqrt{\left( \frac{T}{2.\Delta F}\right)}. \left [ [C(X_1)+C(X_2)]^2 +[S(X_1)+S(X_2)]^2 \right ]^ \frac{1}{2}

The product T.ΔF, which occurs frequently in the equations, is often referred to as the time-bandwidth product of the chirp.
Tables of the Fresnel integrals have been published[7]:32–35[8]:321–322, together with mathematical routines with which to compute the integrals manually or by means of a computer program. In addition, a number of mathematical software programs, such as Mathcad, MATLAB and Mathematica have built-in routines to evaluate the integrals, either as standard functions or in extension packages.

Some plots of the power spectrum |S(ω)|2 as a function of frequency are shown, for time-bandwidth products of 25, 100, 250 and 1000. When the product is small, the Fresnel ripples are very much in evidence, but the spectrum does tend to a more rectangular profile for larger values.

Spectra of Linear Chirps TB=25,100.png


Spectra of Linear Chirps TB=250,1000.png

Non-linear chirps[edit]

When the Frequency-Time characteristic is non-linear, the Fourier integral is difficult to evaluate. In such cases, it is possible to resort to an approximation method such as the stationary phase approximation, or to use numerical methods.

The stationary phase method[edit]

Often (as in radar applications) a(t) is a slowly varying function of time and the phase θ(t) is oscillatory and varies rapidly, over the range of integration. With such waveforms, the stationary phase approximation can be used to investigate the spectrum[12]:34.[14][15][16] The method relies on the fact that the major contributions to the Fourier integral come from the region where the rate of change of phase is minimal, i.e. when

\frac{d}{dt}[(\omega_0-\omega)t + \theta(t)]=0]\qquad or\qquad (\omega - \omega_0) -\theta'(t)=0

Unless θ(t) is a constant, the point in time ts at which the phase is stationary will vary according to the instantaneous frequency ωs.
Expressing the difference between (ωs-ω0).t and θ(t) as a Taylor series about the time ts, but discarding all but the first three terms (of which the second term is zero, here), the Fourier integral can be written, approximately, as

S(\omega)\approxeq a(t_s)\int_{t_s-\delta}^{t_s+\delta}exp \left [ -j \left \{ (\omega_s-\omega_0).t-\theta(t_s)-\frac{\theta''(t_s)}{2}.(t-t_s)^2\right \} \right ] .dt

In this equation ts represents a constant time point, so terms depending on ts alone may be taken outside the integral. The expression simplifies to[12]:39[16]
S(\omega)\approxeq \sqrt{2\pi}.\frac{a(t_s)}{\sqrt {|\theta''(t)|}}.exp \left [j \left \{(\omega_0-\omega_s)t+\theta(t_s)+\frac{\pi}{4} \right \} \right ]
so
|S(\omega_t)|^2 \approxeq 2\pi.\frac{a^2(t)}{|\theta''(t)|}

where ωt is used to indicate the dependence of the frequency variable on t.
This is a very useful expression linking, as it does, the spectrum profile to the amplitude and phase characteristics of the chirp.

To carry out the inverse process, i.e. to find the time domain function s(t) given frequency domain data, the inverse Fourier transform is derived.

s(t)=\frac{1}{2\pi}\int_{-\infty}^{\infty}|S(\omega)|.exp[j(\Phi(\omega)+\omega.t)].d\omega

where Φ(x) is the phase function of the spectrum. The stationary phase points for this integrand are located at

\Phi'(\omega)=-t

and the corollary relationship, equivalent to that derived for the spectrum, can be obtained by the stationary phase method, and is

a^2(t_\omega)\approxeq \frac{1}{2\pi}.\frac{|S(\omega)|^2}{|\Phi''(\omega)|}

In effect, stationary phase analysis gives the following (approximate) Fourier pair relationships[12]:43:
a(t).exp[j\theta(t)] \approxeq \frac{1}{2\pi}.\int_{-\infty}^{\infty}S(\omega)|.exp[j \{\Phi(\omega)+\omega.t \}].d\omega
and
|S(\omega)|.exp[j\Phi(\omega)] \approxeq \int_{-\infty}^{\infty}a(t).exp[-j\{\omega t-\theta(t) \}]dt

Consequently, approximate expressions for a(t) and θ(t) can be obtained when the spectrum, including its phase function Φ(ω) is given and, similarly, approximate expressions for |S(ω| and Φ(ω) can be obtained when the signal characteristics are given. Several examples of the procedure are given in the literature[12]:43[14][16]

Although the relationships are only approximate, their accuracy improves as the time-bandwidth product increases. In cases where the signal envelope and spectrum modulus are defined by smoothly varying Gaussian function then a T.ΔF product as low as 15 will give acceptable results, but if both a(t) and |S(ω)| are defined by rectangular functions, then the product T.ΔF needs to be much greater, typically over 100[12]:49.

Examples of the stationary phase method[edit]

Typically, in the radar case, a(t) is a constant over the duration of the signal and, for convenience, is assumed here to be unity. So the phase and amplitude characteristics, in the frequency domain, are related by

\Phi''(\omega)=\pm \frac{1}{2\pi}.|S(\omega)|^2

There are two solutions for Φ(ω), which are complex conjugates of each other. The two filters with these characteristics can be used as the transmitter and receiver filters of a radar system and are interchangeable.
The group delay characteristic D(ω), (where D(ω)=-dΦ/dω), is

D(\omega)=-\Phi'(\omega)=-\int_{0}^{\infty}\Phi''(\omega).d\omega + K
so
D(\omega)=\pm \frac{1}{2\pi}.\int_{0}^{\infty}|S(\omega)|^2.d\omega+K

So in the case of a rectangular time envelope, the dispersive delay characteristic is given by the integral of the square of the envelope.[16] If the positive sign is taken, then the group delay increases with increasing frequency and vice versa. The result is only approximate, but is more accurate for large values of the time bandwidth product.
Consider, as an example, the case of a spectrum that is uniform over the range -ωmax/2 to ωmax/2, then

|S(\omega)|^2=A \qquad for \qquad |\omega|< \frac{\omega_{max}}{2}
so
D(\omega)=\frac{1}{2\pi}.\int A.d\omega+K = \frac{A}{2\pi}.\omega + K

Put D(-ωmax/2) = 0 and D(ωmax/2 = T, where T is the pulse duration, then K = T/2 and A = (2πT)/ωmax
so, finally
D(\omega)= T. \left[ \frac{1}{2} +\frac{\omega}{\omega_{max}} \right ]

As expected, a flat topped frequency spectrum corresponds to a linear frequency sweep.

The linear chirp is just one special case which, in any case, can be calculated more precisely by the methods of the earlier section. The particular usefulness of the stationary phase method lies in its ability to provide results when the frequency sweep is non linear. In such cases the spectral response can be shaped to meet some desired design criteria, for example, low side-lobes when a chirp is compressed. One such family of spectral functions that has been studied[12]:51 is given by

|S(\omega)|^2 = A_n.cos^n \left (\frac{\pi \omega}{\omega_{max}} \right ) \qquad where \qquad |\omega|< \frac{\omega_{max}}{2} \qquad and\ n\ is\ an\ integer

It is possible to find the group delay characteristics of these functions in a similar manner to that carried out above and the results for n = 1 to 4 have been calculated[12]:51.
Although these cosine functions are amenable to mathematical manipulation, they are rarely chosen to define the spectral characteristics of a chirp, in practice, because when compressed they give broad main pulses with high side-lobe levels. A better characteristic (among many)[17] is the Hamming function, given by

|S(\omega)|^2 = A. \left [ 0.54+0.46.cos \left ( \frac{2\pi \omega}{\omega_{max}} \right ) \right ] =  A. \left [ 0.08 + 0.92.cos^2 \left ( \frac{\pi \omega}{\omega_{max}} \right ) \right ]

Spec. having Hamming Profile.png

A plot of this characteristic is shown, plotted over the range -ωmax/2 to ωmax/2.

Applying the equations given above, the group delay characteristic which achieves this spectral shape can be obtained. It is

D_H(\omega)= T. \left [ \frac{1}{2} + \frac{\omega}{\omega_{max}} + \frac{1.7037}{4\pi}.sin \left ( \frac{2\pi \omega}{\omega_{max}} \right ) \right ]

Now because the principle of stationary phase shows that there is a direct relationship between elapsed time and the instantaneous signal delay then, for the Hamming window, t/T can be related to ω/ωmax by

\frac{t}{T}= \frac{1}{2} + \frac{\omega}{\omega_{max}} + \frac{1.7037}{4\pi}.sin \left ( \frac{2\pi \omega}{\omega_{max}} \right )

This characteristic which is time as a function of frequency is shown here. Inverting the plot gives the more usual (and more useful) plot of frequency as a function of time, which is also shown.

Double Graph for Frequency and Time.png


Other spectral shapes can be investigated in the same way and the results, although approximate, are surprisingly accurate, especially when the time bandwidth product of the pulse is high.

The stationary phase method does not predict or deal with Fresnell ripples, so it is unable to offer any means by which these ripples can be minimized. As an example, the figure below shows a chirp spectrum with T.ΔF =250 obtained for a non-linear chirp aiming to match the Hamming window, using the methods described above. The figure shows that the spectral profile matches the Hamming characteristic quite well, but Fresnell ripples, not predicted by the method, are very much in evidence.

NonLinear Chirp with Hamming Profile, TB=250.png

Numerical methods[edit]

Sampling[edit]

Whenever a Fourier integral cannot be evaluated by analytical means, an approximate solution is usually possible by numerical analysis. Such a procedure requires the function to be sampled, usually at equi-spaced intervals in time.
One consequence of sampling is that the resultant spectrum is periodic in the frequency domain. In addition to the (desired) baseband spectrum, additional versions of the spectrum occur, centered on multiples of the sampling frequency. To ensure that there is no overlapping of frequency data (i.e. no aliasing) the Nyquist sampling theorem must be satisfied. In practice, a sampling rate substantially higher than that dictated by the sampling theorem is advisable[18]:11

Spectrum of a sampled Signal - the Fourier transform of a discrete time signal[edit]

A straightforward way to approximate an integral, such as a Fourier integral, is to use the standard 'rectangle rule' for numerical integration. The method assumes the signal value taken at a sample instant remains constant for one sampling interval, until the next sample is taken. This procedure is sometimes referred to as a 'box-car generator', or a zero order sample and hold[19]:114[20]:34. If the time interval between samples is W, then sn = s(nW), and the desired integral is obtained, approximately, by summing the rectangular areas.
The result so obtained is the convolution of a rectangular pulse with step size W with the impulses located at the sampling instants with weights equal to the sample values[18]:12. In consequence, the spectrum of interest will have superimposed upon it the frequency response of the sample and hold[19]:135[20]:36, and the spectrum of the sampled singnal Ss is given by[18]:12:

Ss(\omega))= \frac{W(sin \omega W/2)}{\omega W/2}. \left [\sum_{n=-\infty}^{\infty} s_n.exp(-jn \omega W) \right ]

The first part of the expression, i.e. the 'sin(x)/x' part, is the frequency response of the sample and hold. Its amplitude decreases with frequency and it falls to 63% of its peak value at half the sampling frequency and it is zero at multiples of that frequency (since fs =1/W).
The second term in the equation is called the Fourier transform of the discrete signal sn[18]:12.[21] It is a continuous function over all ω and involves an infinite number of summations. In practice the summation process can be truncated to a finite number of samples, N, possibly because the waveform is periodic or zero outside the range of samples. Furthermore, because the same spectrum is endlessly repeated, it is possible to confine interest to spectral data within the range -ωs/2 to +ωs/2.

As an example, an exponential chirp (with its top frequency well below the Nyquist limit) is sampled at 256 points, as shown.

Plot of Exponential Chirp, N=256.png

The sampled spectrum, Ss(ω) of this waveform, calculated using the equation given above, is shown. To simplify the plot, only the results at positive frequencies have been displayed. The influence of the frequency spectrum of the zero order hold circuit is clearly seen in the diagram.

Spectrum of Exponential Chirp, N=256.png

The baseband portion of the spectrum is shown in more detail in the next figure and the response shows a distinct slope, being significantly lower at the higher frequencies.

Spectrum of Exponential Chirp, N=256 (detail).png


Although the characteristic of the zero order hold has a small influence on this result, the slope is mainly due to the properties of the chirp. The waveform sweeps relatively quickly over the high frequencies and spends more time sweeping the low frequencies, consequently there is less energy content at the high frequencies with more at the lower ones. (A linear chirp, on the other hand, has a nominally flat spectrum because its frequencies are swept at the same rate, as shown in some earlier plots).

The discrete Fourier transform[edit]

If we limit interest in the output spectrum to a finite number of discrete data points (= N), at frequencies ωm given by
\omega_m= \frac{2\pi m}{NW} \qquad for \qquad m=0,1,2,...N-1

then the formula for calculating the discrete Fourier transform is

 Ss_m \quad = Ss \left ( j \frac{2\pi m}{NW} \right ) \quad = \quad \sum_{n = 0}^{N-1} s_n.exp(-j \left( \frac{2\pi mn}{N} \right )

The calculations can be carried out by means of a straightforward computer algorithm[18]:21, but this is not very efficient in computer usage. Consequently, more efficient algorithms have been developed, especially Fast Fourier Transforms (FFT). Computer programs which implement the FFT are widely available in the literature[18]:54[21]:119,412[22] and in proprietary CAD programs such as Mathcad, MATLAB, and Mathematica.
In the following example a linear chirp with time bandwidth product of 25 is sampled at 128 points (i.e. N = 128). In the figure samples of the real part of the waveform are shown. The FFT process assumes the waveform is cyclic, so these 128 data points can be considered to be part of an endlessly repeating sequence.

Linear chirp with TB=25 and N=128

By calculating the N-point FFT of this data, the discrete spectrum of the sequence is obtained. The magnitude of this spectrum is shown in the attached figure. The data is cyclic, so the zero frequency point is at n = 0 and also at n = 64 (i.e. both points are the same frequency). The point n = 64 corresponds to +fs/2 (and also to -fs/2).

Spectrum of Linear Chirp, TB=25, N=128

To display the spectrum in more detail (but not necessarily with more resolution[23]), the time sequence can be extended by zero padding[21]:80–85.[24][25] For example, extending the 128 point time sequence with zeros to give N = 4096 results in that part of the spectrum originally presented in 16 samples, now being presented in 512 samples, as shown.

Spectrum of Linear Chirp, TB=25, N=4096, (Detail).png

The spectral spread of a chirp[edit]

There is very little spectral content beyond the sweep frequency range of a chirp pulse and this is especially true for waveforms where the time-bandwidth product is large. The full line on the graph of the adjacent figure shows results for linear chirps. It shows, for example, that only about 2% of the total power resides at frequencies outside the sweep range ΔF when the time-bandwidth is 100, and it is less than 1/2% when T.ΔF is 500.
In the case of a non-linear chirp, or a linear chirp shaped by amplitude weighting, the fraction of power outside ΔF is even lower, as is shown on the graph, where the dashed line is for spectra with Hamming profiles.
This low spectral spread is particularly significant when baseband signals are to be digitized since it permits a sampling frequency to be chosen which is only slightly higher than twice the maximum frequency excursion of the chirp.

Fraction of Power Outside Sweep Range.png

Reduction of spectral ripple[edit]

The Fresnel ripples on a chirp spectrum are very obtrusive, especially when time-bandwidth products are low (under 50, say) and their presence leads to high time sidelobe levels when chirps are subject to pulse compression as in radar and sonar systems. They arise because of the sudden discontinuities in the chirp waveform at the commencement and termination of the pulse. Fortunately, there are a number of procedures that can be applied to reduce the ripple levels, but they are not all equally effective. Furthermore, some of the methods require amplitude shaping, or amplitude modulation, of the chirp pulse and this makes those methods unsuitable when, for example, the chirp pulses are to be transmitted by a power amplifier operating in a near-limiting condition. For such systems only the methods using frequency (or phase) pre-distortion are appropriate.

Introducing rise and fall times of finite duration[edit]

If the transitions at the start and end of the chirp are made less sudden (or more 'rounded'), then a reduction in ripple amplitude is achieved[12]:213.[26][27] The durations of the two transition regions need only be a small fraction of the pulse duration, and suggested values are between 2/ΔF and 3/ΔF [26] but, as expected, when the time-bandwidth product of the pulse is small, longer transition periods are needed. The actual profiles of these rise and fall regions of a pulse do not seem to be critical and may be provided, for example, by band limiting filters in analogue implementations and a linear slope in digital ones.
Two examples show the spectra of linear chirps with finite rise-times. The first is for a chirp with time-bandwidth of 250, where the rise and fall times are 4% of the total pulse duration and the second is for a chirp with time-bandwidth of 25, where the rise and fall times are 10% of the total. These two spectra show a marked reduction in ripple amplitude compared to the spectra of unmodified linear chirps shown earlier.

Chirps with Finite Rise and Fall Times.png

Applying phase or frequency distortion to the chirp pulse[edit]

A analogous technique can be applied to the frequency characteristic of the chirp waveform by adding linear FM distortion segments (quadratic phase modulation distortion) to the frequency characteristic of the chirp, as shown. The method is effective because amplitude and phase distortions having functional similarity can produce similar effects when the distortion factors are small.[26][28]

Freq Plot of Chirp with Freq Predistortion.png

Suggested values for these distortion regions, to give good results are:
\Delta f_p=0.75 \Delta F \qquad and \qquad \delta = 1/\Delta F

Later work[29] proposed slightly different values, namely:
\Delta f_p=0.73 \Delta F \qquad and \qquad \delta = 0.86/\Delta F
but the outcome can doubtless be improved by optimizing values for each particular situation.
Two plots show the effects of frequency pre-correction and can be compared to the results in the earlier sections.

Spectra with Freq Predistortion TB=250,25.png

The ripple reduction achieved by frequency pre-correction, although significant, is seen to be less successful than that achieved by the amplitude modulation methods of the previous section. However, it has been suggested[27] that by implementing cubic (rather than quadratic) phase pre-correction, comparable results can be achieved.

Deriving a waveform from a target frequency spectrum[edit]

This method uses an inverse Fourier transform in order to derive a waveform which has a spectrum with the phase characteristic of a chosen chirp but a new amplitude profile which is rectangular and ripple free. The method is very effective but, unfortunately, the waveform that is so derived has a semi-infinite time duration. If, for convenience, the newly derived waveform is truncated to a practical length, then some ripple is reintroduced onto the spectrum.
As an example, a linear chirp waveform with a time bandwidth of 25 is shown together with its spectrum magnitude (shown by a full line) which, as demonstrated earlier, has a large ripple component. It is possible to find, by means of an inverse FFT, a chirp waveform which, in the frequency domain, has the same phase characteristic as before, but with the rectangular magnitude characteristic shown by the dashed line on the plot. The chirp waveform resulting from this process has a very long time duration, but when it is truncated to say, a length 2T, then the spectrum acquires some ripple once more, as shown.

Chirp Waveform TB=25 and Target Spectrum.png
Truncated Chirp showing Wfm and Spectrum, TB=25.png

Applying window functions[edit]

There are many applications in which a spectrum with a rectangular magnitude profile is not ideal. For example, when a chirp waveform is compressed by means of its matched filter, then the resultant waveform approximates to the sinc function and, consequently, has annoyingly high sidelobes. Often, to improve the characteristics of the pulse and lower the sidelobe levels, its spectrum is modified, typically to a bell-shaped profile. Similar problems arise in digital signal processing where the spectral shaping is provided by a window function, a process sometimes called apodization. In the case of an antenna array, similar profiling by "weighting functions" is used to reduce the spatial sidelobes of the radiation pattern.
Although spectral shaping of a chirp could be applied in the frequency domain, better results are obtained if the shaping is carried out in the time domain.[30][31]
Examples of this process are shown for linear chirps with time-bandwidth products of 250 and 25. They have been shaped by a 3-term Blackman-Harris window[17] given by
 |U(\omega)|^2 = 0.42323-0.49755cos \left(\frac{2 \pi \omega}{\omega_{max}} \right)+0.07922cos \left(\frac{4 \pi \omega}{\omega_{max}} \right)
The spectra, now bell-shaped, are seen to be free of ripples.

Chirp and Spec TB=250 BH wgt.png
Chirp and Spec TB=25 BH wgt.png

Non-linear chirps can be devised that have a bell shaped spectrum, such as the Blackman-Harris window just discussed, and consequently will exhibit reduced ripple compared to the linear chirp. By means of the stationary phase method described earlier, an approximate relationship between time and frequency can be obtained and is:
 \frac{t}{T}= \frac{1}{2} + \frac{\omega}{\omega_{max}}+0.1871 sin \left(\frac{2 \pi \omega}{\omega_{max}} \right)+0.014895 sin \left(\frac{4 \pi \omega}{\omega_{max}} \right)

Frequency Plot for Non-linear Chirp with B-H wgt.png

Rearranging the equation, a plot of frequency against time can be plotted, as shown.

As examples, plots of the spectral magnitudes of non-linear chirps with spectral profiles of Blackman-Harris windows and with time-bandwidth products of 250 and 25 are shown below. As can be seen, there is some ripple reduction, but the disappointing performance can be attributed to the fact that these chirps, although they have reduced energy content in their outer frequency regions, they still have amplitude profiles with fast rise and fall times.

Non-linear Frequency Sweep for Chirp with B-H Wgt.png

See also[edit]

References[edit]

  1. ^ Weisstein, Eric W. "Sweep Signal." From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/SweepSignal.html
  2. ^ http://books.google.com.br/books?id=QuIHjnXQqM8C&lpg=PA703&dq=chirp%20quadratic%20phase&pg=PA703#v=onepage&q&f=false
  3. ^ http://setiathome.berkeley.edu/ap_chirp.php
  4. ^ http://www.dspguide.com/ch11/6.htm
  5. ^ Mann, Steve and Haykin, Simon; The Chirplet Transform: A Generalization of Gabor's Logon Transform; Vision Interface '91.[1]
  6. ^ The Beginner's Handbook of Amateur Radio By Clay Laster
  7. ^ a b Jahnke E. and Emde F., "Tables of functions", Dover Publications N.Y. 1945
  8. ^ a b Abramowitz M. and Stegun I.A.,"Handbook of Mathematical Functions", Nat. Bur. Standards 1964, reprinted by Dover Publications N.Y. 1965 (9th ed.1972)
  9. ^ Klauder J.R., Price A.C., Darlington S. and Albersheim W.J., The Theory and Design of Chirp Radars", The Bell system Technical Journal, Vol.39, July 1960 (pp.745-809)
  10. ^ Chin J.E. and Cook C.E., The Mathematics of Pulse compression", Sperry Eng. Review, Vol.12, Oct 1959. (pp.11-16)
  11. ^ Cook C.E., Pulse Compression - Key to More Efficient Radar Transmission",Proc.IRE, March 1960 (p.312)
  12. ^ a b c d e f g h i Cook C.E. and Bernfeld M., "Radar Signals - An Introduction to Theory and Application", Academic Press 1967,1987, reprinted by Artech House 1993.
  13. ^ Varhney L.R. and Thomas D.,"Sidelobe Reduction for Matched Filter range Processing", IEEE Radar Conference 2003
  14. ^ a b Fowle E.N., "The design of FM pulse compression signals", IEEE Trans. IT-10, 1964, (pp.61-67)
  15. ^ Key E.L., Fowle E.N., Haggarty R.D.., "A method of pulse compression employing nonlinear frequency modulation", M.I.T. Lincon Lab., Lexington, Mass., Tech. Rep. 207, 1959.
  16. ^ a b c d Key E.L., Fowle E.N., Haggarty R.D.., "A Method of Designing Signals of Large Time-Bandwidth Product", Proc. IRE Int. Con. Rec. Pt.4, March 1961 (pp.146-154)
  17. ^ a b Harris F.J., "On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform", Proc. IEEE Vol.66, Jan 1978
  18. ^ a b c d e f Burrus C.S. and Parks T.W., "DFT/FFT and Convolution Algorithms", Wiley & Sons, Interscience 1985.
  19. ^ a b Tou J.T., "Digital and Sampled-data Control Systems", McGraw-Hill N.Y. 1959
  20. ^ a b Ragazzini J.R. and Franklin G.F., "Sampled-Data Control Systems", McGraw-Hill N.Y. 1958
  21. ^ a b c Stearns S.D. and Hush D.R., "Digital signal Analysis", Prentice-Hall, 1990 (p.61)
  22. ^ Harris F.J. "Fast Fourier Transforms", San Diego State Univ. Cal. 1984
  23. ^ Anon, "Zero Padding does not buy Spectral Resolution" National Instruments 2006, http://www.ni.comwhite-paper/4880/en/
  24. ^ Harris F.J. "Signal Processing with Ones and Zeros and the FFT", San Diego State Univ., Cal. 1984
  25. ^ Lyons R., "How to Interpolate in the Time Domain by Zero-Padding in the Frequency Domain", http://www.dspguru.com/dsp/how-to-interpolate-in-time-domain-by-zero-padding-inthe-frequency-domain
  26. ^ a b c Cook C.E. & Paolillo J., "A Pulse Compression Predistortion Function for Efficient Sidelobe Reduction in a High-Power Radar", Proc. IEEE Vol.52, April 1964 (pp.377-384)
  27. ^ a b Kowatsch M. and Stocker H.R., "Effect of Fresnel ripples on sidelobe suppression in low time-bandwidth product linear FM pulse compression", IEE Proc. Vol. 129, Pf. F, No.1 Feb 1982
  28. ^ Wheeler H.A., "The Interpretation of Amplitude and Phase Distortion in Terms of Paired Echoes", Proc. IRE, June 1939
  29. ^ Solal M., "High Performance SAW Delay Lines for Low Time Bandwidth Using Periodically Sampled Transducers", Ultrasonics Symposium, IEEE, Nov. 1988.
  30. ^ Judd G.W., "Technique for Realizing Low Time Sidelobe Levels in Small Compression Ratio Chirp Waveforms", Proc. IEEE Ultrasonics Symposium, 1973, pp.478-483
  31. ^ McCue J.J.G., "A Note on the Hamming Weighting of Linear-FM Pulses", Proc. IEEE, Vol. 67, No. 11, Nov 1979, pp.1575-1577.

External links[edit]