Radar Systems Analysis and Design Using MATLAB 3ed. Edition Bassem R. Mahafza
Radar Systems Analysis and Design Using MATLAB 3ed. Edition Bassem R. Mahafza
com
https://ebookname.com/product/radar-systems-analysis-and-
design-using-matlab-3ed-edition-bassem-r-mahafza/
OR CLICK BUTTON
DOWLOAD NOW
https://ebookname.com/product/radar-systems-analysis-and-design-
using-matlab-second-edition-bassem-r-mahafza/
https://ebookname.com/product/radar-signal-analysis-and-
processing-using-matlab-1st-edition-bassem-r-mahafza/
https://ebookname.com/product/matlab-simulations-for-radar-
systems-design-1st-edition-elsherbeni/
https://ebookname.com/product/object-oriented-systems-analysis-
and-design-using-uml-4th-edition-simon-bennett/
Numerical Analysis Using MATLAB and Spreadsheets 2nd
Edition Steven T. Karris
https://ebookname.com/product/numerical-analysis-using-matlab-
and-spreadsheets-2nd-edition-steven-t-karris/
https://ebookname.com/product/vibration-simulation-using-matlab-
and-ansys-1st-edition-michael-r-hatch/
https://ebookname.com/product/numerical-analysis-using-matlab-
and-spreadsheets-2nd-ed-edition-steven-t-karris/
https://ebookname.com/product/micromechatronics-modeling-
analysis-and-design-with-matlab-2nd-ed-edition-giurgiutiu/
https://ebookname.com/product/linear-feedback-control-analysis-
and-design-with-matlab-advances-in-design-and-control-1st-
edition-dingyu-xue/
RADAR SYSTEMS
ANALYSIS AND
DESIGN USING
®
MATLAB
THIRD EDITION
BASSEM R. MAHAFZA
deciBel Research Inc.
Huntsville, Alabama, USA
MATLAB® is a trademark of The MathWorks, Inc. and is used with permission. The MathWorks does not warrant the
accuracy of the text or exercises in this book. This book’s use or discussion of MATLAB® software or related products
does not constitute endorsement or sponsorship by The MathWorks of a particular pedagogical approach or particular
use of the MATLAB® software.
CRC Press
Taylor & Francis Group
6000 Broken Sound Parkway NW, Suite 300
Boca Raton, FL 33487-2742
© 2013 by Taylor & Francis Group, LLC
CRC Press is an imprint of Taylor & Francis Group, an Informa business
This book contains information obtained from authentic and highly regarded sources. Reasonable efforts have been
made to publish reliable data and information, but the author and publisher cannot assume responsibility for the valid-
ity of all materials or the consequences of their use. The authors and publishers have attempted to trace the copyright
holders of all material reproduced in this publication and apologize to copyright holders if permission to publish in this
form has not been obtained. If any copyright material has not been acknowledged please write and let us know so we may
rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced, transmitted, or uti-
lized in any form by any electronic, mechanical, or other means, now known or hereafter invented, including photocopy-
ing, microfilming, and recording, or in any information storage or retrieval system, without written permission from the
publishers.
For permission to photocopy or use material electronically from this work, please access www.copyright.com (http://
www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC), 222 Rosewood Drive, Danvers, MA 01923,
978-750-8400. CCC is a not-for-profit organization that provides licenses and registration for a variety of users. For
organizations that have been granted a photocopy license by the CCC, a separate system of payment has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are used only for
identification and explanation without intent to infringe.
Visit the Taylor & Francis Web site at
http://www.taylorandfrancis.com
and the CRC Press Web site at
http://www.crcpress.com
Table of Contents
Preface
PART I: Radar Principles
Chapter 1: Definitions and Nomenclature, 3
1.1. Radar Systems Classifications and Bands, 3
1.2. Pulsed and Continuous Wave (CW) Radars, 8
1.3. Range, 9
1.4. Range Resolution, 12
1.5. Doppler Frequency, 14
1.6. Coherence, 20
1.7. Decibel Arithmetic, 21
Problems, 23
Appendix 1-A: Chapter 1 MATLAB® Code Listings, 25
Function “pulse_train.m” Listing, 25
Function “range_resolution.m” Listing, 25
Function “doppler_freq.m” Listing, 26
vii
viii Table of Contents
Chapter 10: Moving Target Indicator (MTI) and Pulse Doppler Radars, 361
10.1. Clutter Power Spectrum Density, 361
Table of Contents xiii
Bibliography, 735
Index, 743
Preface
In the year 2000 the first edition of Radar Systems Analysis and Design Using MATLAB®1
was published. It was developed and organized based on my years of teaching graduate level
courses on radar systems analysis and design including advanced topics in radar signal pro-
cessing. At the time, the primary motivation behind the book was to introduce a college-suit-
able comprehensive textbook that provides hands-on experience with MATLAB companion
software. This book very quickly turned into a bestseller, which prompted the publication of its
second edition in the year 2005. The second edition continued in the same vein as its predeces-
sor. It was updated, expanded, and reorganized to include advances in the field and to be more
logical in sequence. New topics were introduced in the body of the text, and much of the MAT-
LAB code was updated and improved upon to reflect the advancements of the latest MATLAB
release.
Since the publication of the first edition, Radar Systems Analysis and Design Using MAT-
LAB filled a void in the market by presenting a comprehensive and self-contained text on radar
systems analysis and design. It was the first book on the market to provide companion MAT-
LAB software to support the theoretical and mathematical discussion found within the pages
of the text. These features were also supported with a detailed solutions manual of all end-of-
chapter problems. This book quickly became the standard adopted by many books published
on the subject; none of which, however, matched the clear presentation nor the transparency
offered by this author, particularly when considering the end-of-chapter solutions manual and
the complete and comprehensive set of MATLAB code, which was made available to all of the
book audience without any restrictions. Users of this book were not only able to reproduce all
plots found in the text, but they also had the ability to change the code by inputting their own
parameters so that they could generate their own specific plots and outputs that met their own
unique academic interest.
In addition to my academic tenure and experience in teaching the subject at the collegiate
level, I have also taught numerous industry courses and conducted many seminars on the sub-
ject of radar systems. Based on this teaching experience, the following conclusion has become
very evident to me: The need and the demand for a comprehensive textbook / reference book
focused on all aspects of radar systems design and analysis remain very strong. Add to this the
1. All MATLAB® functions and programs provided in this book were developed using MATLAB R2011a version
7.12.0.635 with the Signal Processing Toolbox, on a PC with Windows XP Professional operating system. MAT-
LAB® is a registered trademark of the The MathWorks, Inc. For product information, please contact: The Math-
Works, Inc., 3 Apple Hill Drive, Natick, MA 01760-2098 USA. Web: www.mathworks.com.
xix
xx Preface
fact that many college professors have adopted this book as the primary textbook for their
courses on radar systems. Therefore, my desire to write this third edition was turned into real-
ity and has materialized into this product.
It is my view that the third edition of Radar Systems Analysis and Design Using MATLAB is
warranted for the following reasons: (1) bring the text to a more modern status to reflect the
current state of the art; (2) incorporate into the new edition much of the feedback this author
has received from professors using this book as a text and from other practicing engineers; (3)
introduce several new topics that have not found much treatment by other authors, and even
when they did, it was not on a level comparable to the comprehensive and exhaustive approach
adopted by this author in the first two editions; (4) add many new end-of chapter problems; (5)
restructure the presentation to be more convenient for users to adopt the text for either three
graduate-level courses, or one senior-level and two graduate-level courses; and (6) take advan-
tage of the new features offered by the latest MATLAB releases.
Note that all MATLAB code provided in this book was designed as an academic standalone
tool and is not adequate for other purposes. The code was written in a way to assist the reader
in gaining better understanding of the theory. The code was not developed, nor is it intended to
be used as part of an open loop or a closed loop simulation of any kind. The MATLAB code
found in this textbook can be downloaded from this book’s web-page on the CRC Press web-
site. Simply use your favorite web browser, go to www.crcpress.com, and search for keyword
“Mahafza” to locate this book’s web page.
Just like the first and second editions, this third edition provides easy-to-follow mathemati-
cal derivations of all equations and formulas present within the book, resulting in a user
friendly coverage suitable for advanced as well as introductory level college courses. This
third edition provides comprehensive up-to-date coverage of radar systems design and analysis
issues. Users of this book will need only one book instead of several, to gain essential under-
standing of radar design, analysis, and signal processing. This edition contains numerous
graphical plots and supporting artwork. The MATLAB code companion of this edition will
help users evaluate the trade-offs between different radar parameters.
This book is composed of 18 chapters and is divided into 5 parts: Part I, Radar Principles,
Part II, Radar Signals and Signal Processing, Part III, Special Radar Considerations, Part IV,
Radar Detection, and Part V, Radar Special Topics. Part I comprises Chapters 1 and 2. Chapter
1, Definitions and Nomenclature, presents the basic radar definitions and establishes much of
the nomenclature used throughout the text. In Chapter 2, Basic Pulsed and Continuous Wave
(CW) Radar Operations, the radar equation is derived for both pulsed and CW radars, while
other related material such as radar losses and noise are also discussed in details. The radar
equation in the presence of electronic counter measures (ECM) is derived, as well as the
bistatic radar equation.
Part II comprises Chapters 3 through 7. The main thrust of this part of the book is radar sig-
nals or waveforms and radar signal processing. Chapter 3, Linear Systems and Complex Signal
Representation, contains a top-level discussion of elements of signal theory that are relevant to
radar design and radar signal processing. It is assumed that the reader has sufficient and ade-
quate background in signals and systems as well as in the Fourier transform and its associated
properties. Lowpass and bandpass signals are discussed in the context of radar applications.
Continuous as well as discrete systems are analyzed, and the sampling theorem is presented.
Chapter 4, The Matched Filter Radar Receiver, is focused on the matched filter. It presents
the unique characteristic of the matched filter and develops a general formula for the output of
Preface xxi
the matched filter that is valid for any waveform. Chapter 5, Ambiguity Function - Analog
Waveforms, and Chapter 6, Ambiguity Function - Discrete Coded Waveforms, analyze the out-
put of the matched filter in the context of the ambiguity function. In Chapter 5 the most com-
mon analog radar waveforms are analyzed; this includes the single unmodulated pulse, Linear
Frequency Modulation (LFM) pulse, unmodulated pulse train, LFM pulse train, stepped fre-
quency waveforms, and nonlinear FM waveforms. Chapter 6 is concerned with discrete coded
waveforms. In this chapter, unmodulated pulse-train codes are analyzed as well as binary
codes, polyphase codes, and frequency codes. Chapter 7, Pulse Compression, contains details
of radar signal processing using pulse compression. The correlation processor and stretch pro-
cessor are presented. High range resolution processing using stepped frequency waveforms is
also analyzed.
Part III comprises three chapters. Chapter 8, Radar Wave propagation, extends the free
space analysis presented in the earlier chapters to include the effect of the atmosphere on radar
performance. Topics such as refraction, diffraction, atmospheric attenuation, surface reflection,
and multipath are discussed in a fair amount of detail. The subject of radar clutter is in Chapter
9, Radar Clutter. Area clutter as well as volume clutter are defined and the radar equation is re-
derived to reflect the importance of clutter, where in this case, the signal to interference ratio
becomes more critical than the signal to noise ratio. A step-by-step mathematical derivation of
clutter RCS is presented, and the statistical models for the clutter backscatter coefficient is also
presented. Chapter 10, Moving Target Indicator (MTI) and Pulse Doppler Radars, discusses
how delay line cancelers can be used to mitigate the impact of clutter within the radar signal
processor. PRF staggering is analyzed in the context of blind speeds and in the context of
resolving range and Doppler ambiguities. Finally, pulsed Doppler radars are briefly analyzed.
In Part IV, radar detection is discussed and analyzed. The material presented in this part of
the book requires a strong background in random variables and random processes. Therefore,
Chapter 11, Random Variables and Random Processes, presents a review of the subject, and is
written in such a way that it only highlights the major points of the subject. Users of this book
are advised to use this chapter as a means for a quick top-level review of random variables and
random processes. Instructors using this book as a text may assign Chapter 11 as a reading
assignment to their students. Single pulse detection with known and unknown signal parame-
ters is in Chapter 12, Single Pulse Detection. Chapter 13, Detection of Fluctuating Targets,
extends the analysis of Chapter 12 to include target fluctuation where the Swerling target mod-
els are discussed. Detailed discussion of coherent and noncoherent integration in the context of
a square law detector is in this chapter. An overview of CFAR, cumulative probability of detec-
tion, and M-out-of-N detection are also discussed.
Part V of this book addresses a few specialized topics in radar systems. In Chapter 14, Radar
Cross Section (RCS), the RCS dependency on aspect angle, frequency, and polarization are dis-
cussed. A target scattering matrix is developed. RCS formulas for many simple objects are pre-
sented. Complex object RCS is discussed, and RCS prediction methods are introduced.
Chapter 15, Phased Array Antennas, starts by developing the general array formulation. Linear
arrays and several planar array configurations such as rectangular, circular, rectangular with
circular boundaries, and concentric circular arrays are discussed. Beam steering with and with-
out using a finite number of bits is analyzed. Scan loss is also presented. A concept of a multi-
ple input multiple output radar system developed by this author is discussed and analyzed. In
Chapter 16, Adaptive Signal Processing, the concept behind conventual and adaptive beam-
forming is discussed. Adaptive signal processing using the least mean square algorithm is ana-
lyzed. Adaptive linear arrays and complex weights computation in the context of the least
xxii Preface
mean square algorithm are presented. Finally, this chapter discusses, space time adaptive pro-
cessing.
Chapter 17, Target Tracking, discusses target tracking radar systems. The first part of this
chapter covers the subject of single target tracking. Topics such as sequential lobing, conical
scan, monopulse, and range tracking are discussed in detail. The second part of this chapter
introduces multiple target tracking techniques. Fixed gain tracking filters such as the DE and
the DEJ filters are presented in detail. The concept of the Kalman filter is introduced. Special
cases of the Kalman filter are analyzed in depth and a MATLAB-based simulation of the
Kalamn filter is developed. The last chapter of this book is Chapter 18, Tactical Synthetic
Aperture Radars. The topics of this chapter include: SAR signal processing, SAR design con-
siderations, and the SAR radar equation. Arrays operated in sequential mode are discussed in
this chapter.
This book is written primarily as a graduate-level textbook, although parts of it can be used
as a senior level course on radar systems. A companion solutions manual has been developed
for use by professors that adopt this book as a text. This solutions manual is available through
the publisher. Based on my own teaching experience, the following breakdown can be utilized
by professors using this book as a text:
1. Option I: Chapters 1-4 (with omission of certain advanced sections) can be used as a
senior-level course. Chapters 5-10 and the omitted sections in the previous course can be
used as a first graduate level course. Finally, Chapters 11-18 can be used as a second
advanced graduate-level course.
2. Option II: Chapters 1-4 can be used as an introductory graduate-level course. Chapters 5 10
can be used as a second graduate-level course, while Chapters 11-18 can be used as an
advanced graduate course on the subject.
Bassem R. Mahafza
Huntsville, Alabama
United States of America
November, 2012
Part I
Radar Principles
Chapter 1:
Definitions and Nomenclature
Radar Systems Classification and Bands
Pulsed and Continuous Wave (CW) Radars
Range
Range Resolution
Doppler Frequency
Coherence
Decibel Arithmetic
Problems
Appendix 1-A: Chapter 1 MATLAB Code Listings
Chapter 2:
Basic Pulsed and Continuous Wave (CW) Radar Operations
The Radar Range Equation
Low PRF Radar Equation
High PRF Radar Equation
Surveillance Radar Equation
Radar Equation with Jamming
Range Reduction Factor
Bistatic Radar Equation
Radar Losses
Noise Figure
Continuous Wave (CW) Radars
1
2 Radar Systems Analysis and Design Using MATLAB®
Chapter 1
Definitions and
Nomenclature
This chapter presents some basic radar definitions and establishes much of the nomencla-
ture used throughout this text. The word radar is an abbreviation for radio detection and rang-
ing. In most cases, radar systems use modulated waveforms and directive antennas to transmit
electromagnetic energy into a specific volume in space to search for targets. Objects (targets)
within a search volume will reflect portions of the incident energy (radar returns or echoes) in
the direction of the radar. These echoes are then processed by the radar receiver to extract tar-
get information such as range, velocity, angular position, and other target identifying charac-
teristics.
3
4 Radar Systems Analysis and Design Using MATLAB®
Mainly military
radars
High Frequency (HF) and Very High Frequency (VHF) Radars (A- and B-Bands):
These radar bands below 300MHz represent the frontier of radio technology at the time during
the World War II. However, in the modern radar era, these frequencies bands are used for early
warning radars. These radars utilize the electromagnetic waves’ reflection off the ionosphere to
detect targets beyond the horizon, and so they are called Over-the-Horizon Radars (OTHR).
Some examples include the United States (U.S.) Navy Relocatable over-the-horizon Radar
(ROTHR) shown in Fig. 1.2, and the Russian Woodpecker radar shown in Fig. 1.3. By using
these low HF and VHF frequency bands, one can use high-power transmitters. At these fre-
quencies, the electromagnetic wave atmospheric attenuation is small and can be overcome by
using high-power transmitters. Radar angular measurement accuracies are limited in these
bands because lower frequencies require antennas with significant physical size, thus limiting
Radar Systems Classifications and Bands 5
the radar’s angle accuracy and angle resolution. Other communication and broadcasting ser-
vices typically use these frequency bands. Therefore, the available bandwidth for military
radar systems is limited and highly contested throughout the world. Low-frequency systems
can be used for Foliage Penetration (FoPen) applications, as well as in Ground Penetrating
(GPen) applications.
Figure 1.3. Russian Woodpecker OTHR radar. Photograph obtained via the
Internet (http://passingstrangeness.wordpress.com/2010/04/23/the-russian-woodpecker/).
6 Radar Systems Analysis and Design Using MATLAB®
Ultra High Frequency (UHF) Radars (C-Band): UHF bands are used for very long range
Early Warning Radars (EWR). Some examples include the Ballistic Missile Early Warning
System (BMEWS) search-and-track monopulse radar that operates at 245MHz (see Fig. 1.4),
the Perimeter and Acquisition Radar (PAR), which is a very long range multifunction phased
array radar; and the early warning PAVE PAWS multifunction UHF phased array radar. This
frequency band is also used for the detection and tracking of satellites and ballistic missiles
over a long range. In recent years, ultra wideband (UWB) radar applications use all frequencies
in the A- to C-Bands. UWB radars can be used in GPen applications as well as in see-through-
the-wall applications.
L-Band Radars (D-Band): Radars in the L-band are primarily ground-based and ship-
based systems that are used in long range military and air traffic control search operations for
up to 250 (~500Km) nautical miles. Therefore, due to earth curvature their maximum achiev-
able range is limited when detecting low-altitude targets which can disappear very quickly
below the horizon. The Air Traffic Management (ATM) long-range surveillance radars like the
Air Route Surveillance Radar (ARSR), work in this frequency band. These radar systems are
relatively large and demand sizable footprints. Historically, the designator L-Band was
adopted since the “L” represent with large antenna or long range radars.
S-Band Radars (E- and F-Bands): Most ground- and ship-based medium range radars
operate in the S-band. For example, the Airport Surveillance Radar (ASR) used for air traffic
control, and the ship-based U.S. Navy AEGIS (Fig. 1.5) multifunction phased array are S-band
radars, and the Airborne Warning and Control System (AWACS) shown in Fig. 1.6. The atmo-
spheric attenuation in this band is higher than in the D-Band, and they are also more suscepti-
ble to weather conditions. Radar in this band usually need considerably high transmitting
power as compared to the lower-frequency radars in order to achieve maximum detection
range. Even with the considerable weather susceptibility, the National Weather Service Next
Generation Doppler Weather Radar (NEXRAD) uses an S-band radar, because it can see
Radar Systems Classifications and Bands 7
beyond a severe storm. Special Airport Surveillance Radars (ASR) used at some civilian air-
ports are also in this band where they can detect aircrafts for up to 60 nautical miles. The des-
ignator S-Band (contrary to L-Band) was adopted since the “S” represents the smaller antennas
or shorter range radars.
Figure 1.5. U. S. Navy AEGIS. Photograph obtained via the Internet (http://mostlymis-
siledefense.com/2012/08/03/ballistic-missile-defense-the-aegis-spy-1-radar-august-3-2012/).
Figure 1.6. U. S. Air Force AWACS. Photograph obtained via the Internet (http://
www.globalsecurity.org/military/systems/aircraft/e-3-pics.htm).
8 Radar Systems Analysis and Design Using MATLAB®
C-Band Radar (G-Band): Many of the mobile military battlefield surveillance, missile-
control and ground surveillance radar systems operate in this band. Most weather radar sys-
tems are also C-band radars. Medium range search and fire control military radars and metric
instrumentation radars are C-band systems. In this band, the size of the antenna allows for
achieving excellent angular accuracies and resolution. Performance of systems operating in
this band suffer severely from bad weather conditions and to counter that, they often employ
antenna feeds with circular polarization.
X- and Ku-Band Radars (I- and J-Bands): In the X-band frequency range (8 to 12GHz)
the relationship between the wave length and size of the antenna is considerably better than in
lower-frequency bands. Radar systems that require fine target detection capabilities and yet
cannot tolerate the atmospheric attenuation of higher-frequency bands are typically X-Band.
The X- and Ku-bands are relatively popular radar frequency bands for military applications
like airborne radars, since the small antenna size provides good performance. Missile guidance
systems use the Ku-Band (I- and J-Bands) because of the convenient antenna size where
weight is a limiting requirement. Space borne or airborne imaging radars used in Synthetic
Aperture Radar (SAR) for military electronic intelligence and civil geographic mapping typi-
cally use these frequency bands. Finally, these frequency bands are also widely used in mari-
time civil and military navigation radars.
K- and Ka- Band Radars (J- and K-Bands): These high-frequency bands suffer severe
weather and atmospheric attenuation. Therefore, radars utilizing these frequency bands are
limited to short range applications, such as police traffic radars, short range terrain avoidance,
and terrain following radars. Alternatively, the achievable angular accuracies and range resolu-
tion are superior to other bands. In ATM applications these radars are often called Surface
Movement Radar (SMR) or Airport Surface Detection Equipment (ASDE) radars.
Millimeter Wave (MMW) Radars (V- and W-Bands): Radars operating in this frequency
band also suffer from severe high atmospheric attenuation. Radar applications are limited to
very short range of up to a tens of meters. In the W-Band maximum attenuation occurs at about
75GHz and at about 96GHz. Both of these frequencies are used in practice primarily in auto-
motive industry where very small radars (~ 75-76GHz) are used for parking assistants, blind
spot and brake assists. Some radar systems operating at 96 to 98GHz are used as laboratory
experimental or prototype systems.
Pulsed radars use a train of pulsed waveforms (mainly with modulation). In this category,
radar systems can be classified on the basis of the Pulse Repetition Frequency (PRF), as low
PRF, medium PRF, and high PRF radars. Low PRF radars are primarily used for ranging where
target velocity (Doppler shift) is not of interest. High PRF radars are mainly used to measure
target velocity. Continuous wave as well as pulsed radars can measure both target range and
radial velocity by utilizing different modulation schemes. The design, operation, and analysis
of CW and pulsed radar systems are found in subsequent chapters of this book.
1.3. Range
Figure 1.7 shows a simplified pulsed radar block diagram. The time control box generates
the synchronization timing signals required throughout the system. A modulated signal is gen-
erated and sent to the antenna by the modulator/transmitter block. Switching the antenna
between the transmitting and receiving modes is controlled by the duplexer. The duplexer
allows one antenna to be used to both transmit and receive. During transmission it directs the
radar electromagnetic energy toward the antenna. Alternatively, on reception, it directs the
received radar echoes to the receiver. The receiver amplifies the radar returns and prepares
them for signal processing. Extraction of target information is performed by the signal proces-
sor block. The target’s range, R , is computed by measuring the time delay, 't ; it takes a pulse
to travel the two-way path between the radar and the target. Since electromagnetic waves travel
8
at the speed of light, c = 3 u 10 m e s , then
where R is in meters and 't is in seconds. The factor of 1 e 2 is used to account for the two-
way time delay.
In general, a pulsed radar transmits and receives a train of pulses, as illustrated by Fig. 1.8.
The Inter Pulse Period (IPP) is T , and the pulse width is W . The IPP is often referred to as the
Pulse Repetition Interval (PRI). The inverse of the PRI is the PRF, which is denoted by f r ,
During each PRI the radar radiates energy only for W seconds and listens for target returns for
the rest of the PRI. The radar transmitting duty cycle (factor) d t is defined as the ratio
d t = W e T . The radar average transmitted power is
R
Transmitter/ Duplexer
Time Modulator
Control
Signal
Signal
processor
processor Receiver
transmitted pulses
IPP
pulse 1 W pulse 2 pulse 3 time
P av = P t u d t Eq. (1.3)
where P t denotes the radar peak transmitted power. The pulse energy is
E p = P t W = P av T = P av e f r . Eq. (1.4)
The range corresponding to the two-way time delay T is known as the radar unambiguous
range, R u . Consider the case shown in Fig. 1.9. Echo 1 represents the radar return from a target
at range R 1 = c't e 2 due to pulse 1. Echo 2 could be interpreted as the return from the same
target due to pulse 2, or it may be the return from a faraway target at range R 2 due to pulse 1
again. In this case,
T + 't -
R 2 = c't
-------- or R 2 = c---------------------- . Eq. (1.5)
2 2
Clearly, range ambiguity is associated with echo 2. Therefore, once a pulse is transmitted the
radar must wait a sufficient length of time so that returns from targets at maximum range are
back before the next pulse is emitted. It follows that the maximum unambiguous range must
correspond to half of the PRI,
R u = cT e 2 = c e 2f r . Eq. (1.6)
t = 0 t = 1 e fr
PRI
transmitted pulses W pulse 1 pulse 2 time or range
't
echo1 echo 2 time or range
received pulses
c't
R 1 = --------
2 't
Ru
R2
Example:
A certain airborne pulsed radar has peak power P t = 10KW , and uses two PRFs,
f r1 = 10KHz and f r2 = 30KHz . What are the required pulse widths for each PRF so that the
average transmitted power is constant and is equal to 1500watts ? Compute the pulse energy
in each case.
Solution:
Since P av is constant, then both PRFs have the same duty cycle. More precisely,
1500
d t = -------------------3- = 0.15 .
10 u 10
The pulse repetition intervals are
1
T 1 = -------------------3- = 0.1ms
10 u 10
1
T 2 = -------------------3- = 0.0333ms .
30 u 10
It follows that
W 1 = 0.15 u T 1 = 15Ps
W 2 = 0.15 u T 2 = 5Ps
3 –6
E p1 = P t W 1 = 10 u 10 u 15 u 10 = 0.15Joules
3 –6
E p2 = P 2 W 2 = 10 u 10 u 5 u 10 = 0.05Joules .
R max – R min
M = ---------------------------
-. Eq. (1.7)
'R
Targets separated by at least 'R will be completely resolved in range, as illustrated in Fig.
1.10. Targets within the same range bin can be resolved in cross range (azimuth) utilizing sig-
nal processing techniques.
Consider two targets located at ranges R 1 and R 2 , corresponding to time delays t 1 and t 2 ,
respectively. Denote the difference between those two ranges as 'R :
t2 – t1 Gt
'R = R 2 – R 1 = c ------------------ = c ---- . Eq. (1.8)
2 2
Cluster 2
'R 'R
cross range
range
Cluster 1
R min Cluster 3
R max
Now, try to answer the following question: What is the minimum time, Gt , such that target 1
at R 1 and target 2 at R 2 will appear completely resolved in range (different range bins)? In
other words, what is the minimum 'R ?
First, assume that the two targets are separated by cW e 4 , W is the pulse width. In this case,
when the pulse trailing edge strikes target 2, the leading edge would have traveled backward a
distance cW , and the returned pulse would be composed of returns from both targets (i.e., unre-
solved return), as shown in Fig. 1.11a. However, if the two targets are at least cW e 2 apart, then
as the pulse trailing edge strikes the first target, the leading edge will start to return from target
2, and two distinct returned pulses will be produced, as illustrated by Fig. 1.11b. Thus, 'R
should be greater or equal to cW e 2 . And since the radar bandwidth B is equal to 1 e W , then
cW c
'R = ----- = ------- . Eq. (1.9)
2 2B
Range Resolution 13
R1 R2
incident pulse cW
cW
-----
4
reflected pulse return return
tgt1 tgt2
tgt1 tgt2
3
--- cW
2
shaded area has returns from both targets
(a) R1 R2
cW
-----
2
reflected pulses return return
tgt1 tgt2
cW cW
tgt1 tgt2
(b)
Figure 1.11. (a) Two unresolved targets. (b) Two resolved targets.
In general, radar users and designers alike seek to minimize 'R in order to enhance the
radar performance. As suggested by Eq. (1.9), in order to achieve fine range resolution one
must minimize the pulse width. However, this will reduce the average transmitted power and
increase the operating bandwidth. Achieving fine range resolution while maintaining adequate
average transmitted power can be accomplished by using pulse compression techniques.
Example:
A radar system has an unambiguous range of 100Km, and a bandwidth 0.5MHz. Compute the
required PRF, PRI, 'R , and W .
Solution:
8
c- 3 u 10
PRF = -------- = ----------------5- = 1500Hz
2R u 2 u 10
1 - = -----------
PRI = ---------- 1 - = 0.6667ms
PRF 1500
Using the function “range_resolution” yields
8
c 3 u 10
'R = ------- = ------------------------------6- = 300m
2B 2 u 0.5 u 10
u 300
W = 2'R
----------- = 2-----------------
- = 2Ps .
c 3 u 10
8
14 Radar Systems Analysis and Design Using MATLAB®
O Oc
closing target
O ! Oc
radar
O Oc
opening target
O Oc
radar
incident
reflected
Consider a pulse of width W (seconds) incident on a target that is moving toward the radar at
velocity v , as shown in Fig. 1.13. Define d as the distance (in meters) that the target moves
into the pulse during the interval 't ,
where 't is equal to the time between the pulse leading edge striking the target and the trailing
edge striking the target. Since the pulse is moving at the speed of light and the trailing edge has
moved distance cW – d , then
cW = c't + v't Eq. (1.11)
v+c
Wc = ----------- W . Eq. (1.15)
c–v
To derive an expression for Doppler frequency, consider the illustration shown in Fig. 1.14.
It takes 't seconds for the leading edge of pulse 2 to travel a distance c e f r – d to strike the
target. Over the same time interval, the leading edge of pulse 1 travels the same distance c't .
More precisely,
trailing leading
incident pulse edge
edge v
at time t = t 0 L = cW
s = c't
s = c't
d = v't
at time t = t 0 + 't L' = cW'
reflected pulse
leading trailing
edge edge
cW c e fr
v
incident pulse 2 pulse 1
TE LE TE LE
cW' cW d
s – d = c e f r' cW' 2d
c
--- – d = c't . Eq. (1.17)
fr
Solving for 't yields
cef
't = ----------r- Eq. (1.18)
c+v
cv e f
d = -------------r . Eq. (1.19)
c+v
The reflected pulse spacing is now s – d and the new PRF is f r c , where
c cv e f
s – d = ----- = c't – -------------r Eq. (1.20)
fr c c+v
It follows that the new PRF is related to the original PRF by
c+v
f r c = ----------- f r . Eq. (1.21)
c–v
However, since the number of cycles does not change, the frequency of the reflected signal
will go up by the same factor. Denoting the new frequency by f 0 c , it follows that
c+v
f 0 c = ----------- f 0 Eq. (1.22)
c–v
Doppler Frequency 17
where f 0 is the carrier frequency of the incident signal. The Doppler frequency f d is defined as
the difference f 0 c – f 0 . More precisely,
c+v 2v
f d = f 0 c – f 0 = ----------- f 0 – f 0 = ----------- f 0 , Eq. (1.23)
c–v c–v
but since v « c and c = Of 0 , then
2v 2v
f d | ------ f 0 = ------ . Eq. (1.24)
c O
Eq. (1.24) indicates that the Doppler shift is proportional to the target velocity, and, thus, one
can extract f d from range rate and vice versa.
The result in Eq. (1.24) can also be derived using the following approach: Fig. 1.15 shows a
closing target with velocity v . Let R 0 refer to the range at time t 0 (time reference); then the
range to the target at any time t is
R t = R0 –v t – t0 . Eq. (1.25)
xr t = x t – \ t Eq. (1.26)
x r t = x § § 1 + 2v
------· t – \ 0· Eq. (1.28)
©© c¹ ¹
where the constant phase \ 0 is
2R 2v
\ 0 = --------0- + ------ t 0 . Eq. (1.29)
c c
Define the compression or scaling factor J by
J = 1 + 2v e c Eq. (1.30)
Note that for a receding target the scaling factor becomes J = 1 – 2v e c . Utilizing Eq.
(1.30), one can rewrite Eq. (1.28) as
x r t = x Jt – \ 0 . Eq. (1.31)
Eq. (1.31) represents a time-compressed version of the return signal from a stationary target
( v = 0 ). Hence, based on the scaling property of the Fourier transform, the spectrum of the
received signal will be expanded in frequency to a factor of J .
Consider the special case when
where Z 0 is the radar center frequency in radians per second. The received signal x r t is then
given by
18 Radar Systems Analysis and Design Using MATLAB®
R0
1 Z Z
X r Z = ----- § Y § ---- – Z 0· + Y § ---- + Z 0· · , Eq. (1.34)
2J © © J ¹ ©J ¹¹
where for simplicity the effects of the constant phase \ 0 have been ignored in Eq. (1.34).
Therefore, the bandpass spectrum of the received signal is now centered at JZ 0 instead of Z 0 .
The difference between the two values corresponds to the amount of Doppler shift incurred
due to the target motion,
Z d = Z 0 – JZ 0 f d = f 0 – Jf 0 . Eq. (1.35)
Z d and f d are the Doppler frequency in radians per second and in Hz, respectively. Substitut-
ing the value of J in Eq. (1.35) yields
2v
f d = ------ f 0 = 2v
------ , Eq. (1.36)
c O
which is the same as Eq. (1.24). It can be shown that for a receding target, the Doppler shift is
f d = – 2v e O . This is illustrated in Fig. 1.16.
In both Eq. (1.36) and Eq. (1.24) the target radial velocity with respect to the radar is equal to
v , but this is not always the case. In fact, the amount of Doppler frequency depends on the tar-
get velocity component in the direction of the radar (radial velocity). Fig. 1.17 shows three tar-
gets all having velocity v : target 1 has zero Doppler shift; target 2 has maximum Doppler
frequency as defined in Eq. (1.36). The amount of Doppler frequency of target 3 is
f d = 2v cos T e O , where v cos T is the radial velocity, and T is the total angle between the
radar line of sight and the target. Thus, a more general expression for f d that accounts for the
total angle between the radar and the target is
2v
f d = ------ cos T Eq. (1.37)
O
and for an opening target
–2 v
f d = --------- cos T Eq. (1.38)
O
where cos T = cos T e cos T a . The angles T e and T a are, respectively, the elevation and azi-
muth angles; see Fig. 1.18.
Doppler Frequency 19
amplitude
amplitude
fd fd
f0 frequency f0 frequency
closing target receding target
v
v v
T
v
Ta Te
Figure 1.18. Radial velocity is proportional to the azimuth and elevation angles.
20 Radar Systems Analysis and Design Using MATLAB®
Example:
Compute the Doppler frequency measured by the radar shown in the figure below.
O = 0.03m
vtarget = 175 m/sec
line of sight
target
Solution:
The relative radial velocity between the radar and the target is v radar + v t arg et . Thus, using Eq.
(1.36), we get
250 + 175
f d = 2 ----------------------------- = 28.3KHz .
0.03
Similarly, if the target were opening the Doppler frequency is
250 – 175
f d = 2 ------------------------ = 5KHz .
0.03
1.6. Coherence
A radar is said to be coherent if the phase of any two transmitted pulses is consistent, i.e.,
there is a continuity in the signal phase from one pulse to the next, as illustrated in Fig. 1.19a.
One can view coherence as the radar’s ability to maintain an integer multiple of wavelengths
between the equiphase wavefront from the end of one pulse to the equiphase wavefront at the
beginning of the next pulse, as illustrated by Fig. 1.19b. Coherency can be achieved by using a
STAble Local Oscillator (STALO). A radar is said to be coherent-on-receive or quasi-coherent
if it stores in its memory a record of the phases of all transmitted pulses. In this case, the
receiver phase reference is normally the phase of the most recent transmitted pulse.
Decibel Arithmetic 21
(a)
(b) distance
Figure 1.19. (a) Phase continuity between consecutive pulses. (b) Maintaining an integer
multiple of wavelengths between the equiphase wavefronts of any two successive pulses
guarantees coherency.
Coherence also refers to the radar’s ability to accurately measure (extract) the received sig-
nal phase. Since Doppler represents a frequency shift in the received signal, then only coherent
or coherent-on-receive radars can extract Doppler information. This is because the instanta-
neous frequency of a signal is proportional to the time derivative of the signal phase. More pre-
cisely,
1 d
f i = ------ \ t Eq. (1.39)
2S d t
where f i is the instantaneous frequency, and \ t is the signal phase.
For example, consider the following signal:
where the scaling factor J is defined in Eq. (1.30), and \ 0 is a constant phase. It follows that
the instantaneous frequency of x t is
f i = Jf 0 Eq. (1.41)
2v 2v
f i = f 0 § 1 + ------· = f 0 + ------ Eq. (1.42)
© c¹ O
where the relation c = Of is utilized. Note that the second term of the most right-hand side of
Eq. (1.42) is a Doppler shift.
type. The unit dB is named after Alexander Graham Bell, who originated the unit as a measure
of power attenuation in telephone lines. By Bell’s definition, a unit of Bell gain is
P
log § -----0· Eq. (1.43)
© Pi ¹
where the logarithm operation is base 10, P 0 is the output power of a standard telephone line
(almost one mile long), and P i is the input power to the line. If voltage (or current) ratios are
used instead of the power ratio, then a unit Bell gain is defined as
V 2 I 2
log § -----0· or log § ---0-· . Eq. (1.44)
© Vi ¹ © Ii ¹
–1
A decibel, dB, is 1 e 10 of a Bell (the prefix “deci” means 10 ). It follows that a dB is
defined as
P V 2 I 2
10 log § -----0· = 10 log § -----0· = 10 log § ---0-· . Eq. (1.45)
© Pi ¹ © Vi ¹ © Ii ¹
The decibel nomenclature is widely used by radar designers and users for several reasons,
and perhaps, the most important one is that representing radar-related physical quantities using
dBs drastically reduces the dynamic range that a designer or a user has to use. For example, an
–9
incoming radar signal may be as weak as 1 u 10 V , which can be expressed in dBs as
–9
10 log 1 u 10 = – 90dB . Alternatively, a target may be located at range R = 1000Km ,
which can be expressed in dBs as 60dB . Another advantage of using dB in radar design and
analysis is to facilitate the arithmetic associated with calculating the different radar parameters.
This is true since multiplication in base-10 arithmetic translates into addition in dB-arithmetic,
and division translates into subtraction. For example,
250 u 0.0001 = > 10 log 250 + 10 log 0.0001 – 10 log 455 @dB = – 42.6dB
------------------------------- . Eq. (1.47)
455
In general,
AuB
10 log § -------------· = 10 log A + 10 log B – 10 log C Eq. (1.48)
© C ¹
q
10 log A = q u 10 log A . Eq. (1.49)
Other dB ratios that are often used in radar analysis include the dBsm (dB, squared meters).
This definition is very important when referring to target Radar Cross Section (RCS), whose
2
units are in squared meters. More precisely, a target whose RCS is V m can be expressed in
2 2
dBsm as 10 log V m . For example, a 10m target is often referred to as a 10dBsm target,
2
and a target with RCS 0.01m is equivalent to a – 20dBsm .
Problems 23
Finally, the units dBm (dB, milliwatt) and dBW (dB, Watt) are power ratios of dBs with ref-
erence to one milliwatt and one Watt, respectively.
P
dBm = 10 log § -------------· Eq. (1.50)
© 1mW¹
P
dBW = 10 log § --------· Eq. (1.51)
© 1W¹
To find dBm from dBW, add 30dB, and to find dBW from dBm, subtract 30dB. Other common
dB units include dBz and dBi. dBz is used to measure weather radar reflectivity representing
6 –3
the amount of returned power received by the radar referenced to mm m . The unit dBi (dB,
isotropic) represents the forward gain of an antenna compared to an ideal isotropic antenna that
emits energy equally in all directions.
Problems
1.1. (a) Calculate the maximum unambiguous range for a pulsed radar with PRF of 200Hz
and 750Hz . (b) What are the corresponding PRIs?
1.2. For the same radar in Problem 1.1, assume a duty cycle of 30% and peak power of
5KW . Compute the average power and the amount of radiated energy during the first 20ms .
1.3. A certain pulsed radar uses pulse width W = 1Ps . Compute the corresponding range
resolution.
1.4. An X-band radar uses PRF of 3KHz . Compute the unambiguous range and the
required bandwidth so that the range resolution is 30m . What is the duty cycle?
1.5. Compute the Doppler shift associated with a closing target with velocity 100, 200, and
350 meters per second. In each case, compute the time dilation factor. Assume that O = 0.3m .
1.6. Compute the round-trip delays, minimum PRIs, and corresponding PRFs for targets
located 30Km, 80Km, and 150Km away from the radar.
1.7. Assume an S-band radar, what are the Doppler frequencies for the following target
range rates: 50m/s; 200m/s; and 250m/s.
1.8. Repeat the previous problem for an X-Band radar (9.5GHz).
1.9. A certain L-band radar has center frequency 1.5GHz , and PRF f r = 10KHz . What is
the maximum Doppler shift that can be measured by this radar?
1.10. Starting with a modified version of Eq. (1.25), derive an expression for the Doppler
shift associated with a receding target.
1.11. In reference to Fig. 1.18, compute the Doppler frequency for v = 150m e s ,
T a = 30q , and T e = 15q . Assume that O = 0.1m .
1.12. A pulsed radar system has a range resolution of 30cm . Assuming sinusoid pulses at
45KHz , determine the pulse width and the corresponding bandwidth.
1.13. (a) Develop an expression for the minimum PRF of a pulsed radar. (b) Compute f rmin
for a closing target whose velocity is 400m e s . (c) What is the unambiguous range? Assume
that O = 0.2m .
24 Radar Systems Analysis and Design Using MATLAB®
1.14. A certain radar is tasked with detecting and tracking the moon. Assume that the aver-
8
age distance to the moon is 3.844 u 10 m , and its average radar cross section is
11 2
6.64 u 10 m . (a) Compute the delay to the moon. (b) What is required PRF so the range to
the moon is unambiguous. (c) What is the moon’s radar cross section in dBsm.
1.15. An L-band pulsed radar is designed to have an unambiguous range of 100Km and
range resolution 'R d 100m . The maximum resolvable Doppler frequency corresponds to
v t arg et d 350m e sec . Compute the maximum required pulse width, the PRF, and the average
transmitted power if P t = 500W .
1.16. A certain target has the following characteristics: its range away from the radar given
in its corresponding x- y- and z- components is ^ 25Km 32Km 12Km ` . The target velocity
vector is v z = v y = 0 , and v x = – 250m e s . Compute the composite target range and range
rate. If the radar’s operating frequency is 9GHz, what is the corresponding Doppler frequency.
Appendix 1-A: Chapter 1 MATLAB Code Listings 25
Chapter 2
Watts
P D = Peak
transmitted power --------------
-.
------------------------------------------------------------------- 2
Eq. (2.1)
area of a sphere m
2
The power density, in Watts e m , at range R away from the radar (assuming a lossless propa-
gation medium) is
2
P D = P t e 4SR Eq. (2.2)
2
where P t is the peak transmitted power and 4SR is the surface area of a sphere of radius R .
Radar systems utilize directional antennas in order to increase the power density in a certain
direction. Directional antennas are usually characterized by the antenna gain G and the
antenna effective aperture A e . They are related by
2
G = 4SA e e O Eq. (2.3)
where O is the radar operating wavelength. The relationship between the antenna’s effective
aperture A e and the physical aperture A is
A e = UA Eq. (2.4)
0dUd1
where U is referred to as the aperture efficiency, and good antennas require U o 1 . In this
book, unless otherwise noted, A and A e are used interchangeably to refer to the antenna’s
aperture, and will assume that antennas have the same gain in the transmitting and receiving
modes. In practice, U | 0.7 is widely accepted.
The gain is also related to the antenna’s azimuth and elevation antenna beamwidths by
4S
G = K ----------- Eq. (2.5)
Te Ta
27
28 Radar Systems Analysis and Design Using MATLAB®
where K d 1 and depends on the physical aperture shape, and the angles T e and T a are,
respectively, the antenna’s elevation and azimuth beamwidths in radians. An excellent com-
monly used approximation of Eq. (2.5) is
26000
G | --------------- Eq. (2.6)
Te Ta
where in this case the azimuth and elevation beamwidths are given in degrees.
The power density at a distance R away from a radar using a directive antenna of gain G is
then given by
Pt G
P D = ------------
-
2
Eq. (2.7)
4SR
When the radar radiated energy impinges upon a target, the induced surface currents on that
target radiate electromagnetic energy in all directions. The amount of the radiated energy is
proportional to the target size, orientation, physical shape, and material, which are all lumped
together in one target-specific parameter called the Radar Cross Section (RCS) denoted sym-
bolically by the Greek letter V .
The radar cross section is defined as the ratio of the power reflected back to the radar to the
power density incident on the target,
P 2
V = ------r m Eq. (2.8)
PD
where P r is the power reflected from the target. Thus, the total power delivered to the radar
signal processor by its antenna is
P t GV
P Dr = -------------------
- Ae . Eq. (2.9)
2 2
4SR
Substituting the value of A e from Eq. (2.3) into Eq. (2.9) yields
2 2
Pt G O V
P Dr = ---------------------
- Eq. (2.10)
3 4
4S R
Let S min denote the minimum detectable signal power by the radar. It follows that the
maximum radar range R max is
2 2 1e4
§ Pt G O V ·
R max = ¨ -----------------------
3
-¸ . Eq. (2.11)
© 4S S min¹
Eq. (2.11) suggests that in order to double the radar maximum range, one must increase the
peak transmitted power P t sixteen times; or equivalently, one must increase the effective aper-
ture four times.
In practical situations the returned signals received by the radar will be corrupted with noise,
which introduces unwanted voltages at all radar frequencies. Noise is random in nature and can
be characterized by its Power Spectral Density (PSD) function. The noise power N is a
function of the radar operating bandwidth, B . More precisely,
The Radar Range Equation 29
N i = kT s B Eq. (2.13)
– 23
where k = 1.38 u 10 Joule e degree Kelvin is Boltzmann’s constant, and T s is the total
effective system noise temperature in degrees Kelvin. It is always desirable that the minimum
detectable signal ( S min ) be greater than the noise power. The fidelity of a radar receiver is nor-
mally described by a figure of merit referred to as the noise figure, F . The noise figure is
defined as
SNR Si e Ni
F = ------------------i = --------------
- Eq. (2.14)
SNR o So e No
where SNR i and SNR o are, respectively, the Signal to Noise Ratios (SNR) at the input
and output of the receiver. The input signal power is S i , and the input noise power immedi-
ately at the antenna terminal is N i . The values S o and N o are, respectively, the output signal
and noise powers.
The receiver effective noise temperature excluding the antenna is
Te = To F – 1 Eq. (2.15)
where T 0 = 290K and F is the receiver noise figure. It follows that the total effective system
noise temperature T s is given by
Ts = Te + Ta = T0 F – 1 + Ta = To F – T o + Ta Eq. (2.16)
In many radar applications it is desirable to set the antenna temperature T a to T 0 and thus,
Eq. (2.16) is reduced to
Ts = To F . Eq. (2.17)
Using Eq. (2.17) in Eq. (2.13) and substituting the result into Eq. (2.14) yields
The radar detection threshold is set equal to the minimum output SNR, SNR o min .
Substituting Eq. (2.19) in Eq. (2.11) gives
2 2 1e4
§ Pt G O V ·
R max = ¨ ------------------------------------------------------
3
¸ Eq. (2.20)
© 4S kT o BF SNR omin¹
or equivalently,
2 2
Pt G O V
SNR o min = ------------------------------------------
3
. Eq. (2.21)
4
4S kT o BFR max
Another random document with
no related content on Scribd:
Juliet started, coloured, and looked round in some
alarm, for Flossie had not thought it necessary to lower her
voice.
"Of course you are. That is what I told him. You are a
darling, Juliet. There is no one like you. You are not to be
tried by ordinary rules."
Juliet felt impatient of this flow of adulation. She looked
back, and saw her mother coming alone.
"I must go back and meet mother," she said; "I cannot
let her walk home alone. Good-bye, Flossie; I shall see you
again, I hope, when we return from North Devon."
That night, after her uncle had gone away, she lingered
alone in the little garden attached to the house in which she
and her mother were lodging. It was growing dusk. Already
stars were appearing in the clear sky above her head. A
light breeze rustled the trees. Behind her lay the vast,
mysterious moor. In front, far down beneath the trees, out
of sight, but making its presence known by the low, distant
moaning of its waves, was the sea.
"Yes," the voice of her better self made answer, "it did
matter. It must be better to take God's way, even if it
seemed steep and hard, for it would lead upward."
CHAPTER IX
GRATIFIED DESIRES
"Poor child! You are hardly used," said her uncle, with a
merry twinkle in his eyes; "but now let me hear more about
this concert. How would it do if I were to take you?"
Juliet did not again turn her pretty head towards the
seat where young Chalcombe sat, but she was aware,
without seeming to be so, that his opera-glasses were often
directed to the spot where she sat. Her uncle observed it
also, and felt enraged with the fellow for his impudence.
When the concert was over, Juliet and her uncle took a
walk in the grounds. Juliet spoke with rapture of the
concert, and the delight with which she had listened to
Adelina Patti.
Not only was Ralph Tracy his own medical man, he liked
to be his own lawyer as well. It vexed him to think that the
simple will which he could have drawn up himself would
hardly meet the necessities of the case. To secure the
property to Juliet, and protect it from unscrupulous hands,
it would be necessary to make careful provisions. He hardly
knew if he were equal to framing them himself, but ere he
called in the aid of a solicitor, he would have a try at it, with
the help of sundry large books crammed full of legal
information which was often very hard to digest.
"You need not save money for that, Juliet. I will pay for
the lessons. Why did you not speak to me about it before?"
"Oh, uncle, you cannot think how happy you make me!
You don't know how I have longed for lessons!"
"He may be all that," said Juliet; "but I do not like good
young men."
CHAPTER X
A PERILOUS PATH
She should have been very happy, one might think. But
human happiness does not consist of anything external, and
it will not come even with the realisation of all one's desires.
And Juliet would have said that she was far indeed from
such an attainment. No acquisition satisfied her. She was
ever reaching after something beyond. Her uncle's
indulgence had the effect of making her more of a spoilt
child than ever. Restless, petulant, and perverse, she was
constantly working herself into a fever over something or
other. She continually destroyed the peace of the household
by her irritability and impatience.
"I asked him how you were getting on," said Algernon,
responding to her look, "and he said you were doing well."
But when Juliet was questioned about it, she could not
say that her master had given her much encouragement.
"And yet I do feel encouraged somehow," she added, with a
sunny smile.
Juliet could not reply. She had felt vexed that he had
waited for her, and she wanted to tell him that he must not
do so again. But she could not say that or anything now.
CHAPTER XI
HIS LAST MESSAGE