Distortion Mitigation in Cognitive Radio Receivers

DistortionMitigationin
CognitiveRadioReceivers
DISTORTION MITIGATION IN COGNITIVE
RADIO RECEIVERS
Dlovan Hoshiar Mahrof
Members of the dissertation committee:
Chairman:
Prof. Dr.
P.M.G. Apers
University of Twente
Prof. Dr.
P.M.G. Apers
University of Twente
B. Nauta
University of Twente
E.A.M. Klumperink
University of Twente
F.E. van Vliet
Henrik SjÖland
C.H. Slump
P.G.M. Baltus
TNO / University of Twente
Lund University
University of Twente
Eindhoven University of Technology
Secretary:
Promotor:
Prof. Dr. ir.
Assistant Promotor:
Dr. ing.
Members:
Prof. Dr. ir.
Prof. Dr.
Prof. Dr. ir.
Prof. Dr. ir.
This research is supported by the Dutch Technology Foundation
STW, which is part of the Netherlands Organisation for Scientific
Research (NWO) and partly funded by the Ministry of Economic
Affairs (08081).
CTIT Ph.D. Thesis Series No. 15-348
Centre for Telematics and Information Technology
P.O. Box 217, 7500 AE
Enschede, The Netherlands.
Title:
Distortion Mitigation in Cognitive Radio Receivers
ISBN: 978-90-365-3846-6
ISSN: 1381-3617 (CTIT Ph.D. Thesis Series No. 15-348 )
DOI: 10.3990/1.9789036538466
http://dx.doi.org/10.3990/1.9789036538466
Copyright © 2015 by Dlovan Hoshiar Mahrof, Enschede, The Netherlands
DISTORTION MITIGATION IN COGNITIVE RADIO
RECEIVERS
DISSERTATION
to obtain
the degree of doctor at the University of Twente,
on the authority of the rector magnificus,
prof.dr. H. Brinksma,
on account of the decision of the graduation committee,
to be publicly defended
on Friday the 20th March 2015 at 12:45
by
DLOVAN HOSHIAR MAHROF
born on the 26th of December 1973
in Al-Sulaimaniya, Iraq
This dissertation has been approved by:
Promotor: Prof. Dr. ir. Bram Nauta
Co-promotor: Dr. ing. Eric A.M. Klumperink
To Almas,
Varto and Zhareen,
Abstract
The exponential increase of wireless communication increasingly leads to spectrum congestion.
Attempts are being made to increase RF spectrum utilization efficiency by introducing Cognitive Radio
(CR) concept. A CR tries to intelligently solve the congestion problem via Dynamic Spectrum Access
(DSA), i.e. determine which frequencies are temporarily and locally free, and exploit this free spectrum.
Especially in the TV broadcasting bands below 1 GHz, such CR possibilities are being explored. Ideally, a
CR receiver should be able to operate directly adjacent to the primary service users, e.g. Digital TV
channels, which use high power levels and often leave the adjacent channels unused. Under such
conditions no or very low up-front filtering of the interferer is possible. Consequently, a CR receiver
must tolerate the existence of strong interferers, i.e. have a very high linearity front-end.
This thesis examines CR receiver linearity requirements and explores techniques that mitigate distortion.
The thesis starts by analyzing the linearity requirements for DSA. As the level and the spectral location of
the interferers can change instantaneously per location, it is relevant to monitor the spectrum and find
suitable opportunities for communication. The analysis is applied to a channelized spectrum in which a
number of interferers exist, and a CR tries to exploit any free channel. It is derived how the CR linearity
requirement depends not just on the power levels of the interferers but also on their spectral locations
around the desired CR frequency/channel. It is shown that the linearity requirement can be relaxed by
tens of dBs levels of 3rd order InterModulation product (IM3). The analysis also exploits the prediction of
the distortions in different channels for DSA. This prediction algorithm is denoted here as DPrA (i.e.
Distortion Prediction Algorithm). It processes the spectrum sensing information about the power level
and the spectral locations of the interferers to derive the linearity requirements for each potential CR
channel. Based on this information, a CR can choose the most suitable channel compatible with its
linearity capability.
A receiver with better linearity can work under worse interference conditions, and hence maximizing
linearity of a CR receiver is important. To increase the linearity of a CR receiver, CMOS receiver frontends with high linearity are explored. Receivers that exploit linear V-I conversion at RF, followed by
passive down-mixing and an OpAmp-based Transimpedance Amplifier at baseband, show high linearity
potential. However, it is shown that due to nonlinearity and finite gain in the OpAmp, the virtual ground
I
is imperfect, resulting in distortion currents. The concept of a negative conductance is proposed to
cancel such distortion currents. Through a simple intuitive analysis, the basic operation of the technique
is explained. By mathematical analysis the optimum negative conductance value is derived and related
to feedback theory. It is shown that a finite value of the negative conductance is needed to make the
feedback loopgain theoretically approach infinity, which is practically not easily possible by increasing
the gain of the OpAmp block. The technique is applied to linearize an RF receiver and a prototype is
implemented in 65nm technology. Measurement results show an increase of the in-band IIP3 from 9dBm
to >20dBm, and IIP2 from 51 to 61dBm, at the cost of an increase of the noise figure from 6 to 7.5dB and
<10% power penalty. The chip achieves a Spurious-Free Dynamic Range of 85dB in 1MHz.
II
Samenvatting
De exponentiële groei van draadloze communicatie leidt steeds meer tot congestie van het spectrum.
Pogingen zijn ondernomen om het gebruik van Radio Frequentie (RF) spectrum efficiënter te maken
door de invoering van Cognitieve Radio (CR) concept. Een CR probeert de congestie-probleem via
Dynamic Spectrum Access (DSA) op te lossen. DSA bepaalt welke frequenties tijdelijk en lokaal vrij zijn
om die toe te wijzen aan CR. De CR operatie mogelijkheden wordt vooral onderzocht in de TV uitzending
banden onder 1 GHz. In ideale geval moet een CR ontvanger kunnen ontvangen direct in de kanalen
naast de primaire servicegebruikers, zoals digitale TV zenders, die hoog vermogen sturen en vaak
verlaten de aangrenzende kanalen ongebruikt. Onder dergelijke omstandigheden zal er geen of zeer
lage pre-filteren van digitale TV signalen plaatsvinden. Bijgevolg moet een CR ontvanger het bestaan van
sterke interfereren tolereren. Dat wil zeggen dat de CR ontvanger een zeer hoge lineariteit moet
hebben.
Deze thesis onderzoekt de benodigde eisen aan de CR ontvanger lineariteit en verkent de technieken,
die de gevolge vervorming/distorsie verzachten.
Het proefschrift begint met het analyseren van de lineariteit voor DSA. Het vermogen niveau en de
spectrale locatie van de interferentie signalen verandert voortdurend per elke locatie, daarom is het
relevant om het spectrum voortdurend te meten en de geschikte vrij frequenties te bepalen voor CR
communicatie. De analyse wordt toegepast op een gekanaliseerd spectrum waarin een aantal kanalen
zijn bezet met interferentie signalen. Een CR probeert om van die vrij kanalen te benutten. Het is
afgeleid hoe afhankelijk de CR lineariteit is van niet alleen de vermogen van de interferentie signalen
maar ook van hun spectrale locaties rond de gewenste CR frequentie/kanaal. Het is aangetoond dat de
lineariteit met tientallen dBs niveaus van 3e orde intermodulatie product (IM3) kan worden versoepeld.
De analyse exploiteert het gebruik van de voorspelling van de verstoringen in verschillende kanalen voor
DSA. Deze voorspelling algoritme wordt ontwikkelt en aangeduid als DPrA (d.w.z. Distortion Prediction
Algorithm) hier. Het verwerkt de spectrum informatie over het vermogen en de spectrale locatie van de
interferentie signalen voor het afleiden van de lineariteit eisen voor elke potentiële CR kanaal. Op basis
van deze informatie, kan een CR het meest geschikt kanaal voor zijn lineariteit kiezen.
Een ontvanger met betere lineariteit kan werken onder slechtere omstandigheden van de interferentie
signalen, en daarom het maximaliseren van CR ontvanger lineariteit is belangrijk. Om de lineariteit van
III
een CR ontvanger te vergroten worden CMOS ontvanger met hoge lineariteit onderzocht. Ontvangers
die gebruikmaken van lineaire spanning naar stroom conversie bij RF, gevolgd door het passieve
frequentie translatie en een OpAmp gebaseerde transimpedantie versterker op lage frequentie, toont
hoge potentiële lineariteit. Het is aangetoond dat als gevolg van niet lineariteit en beperkte versterking
factor in de OpAmp, de virtueel grond onvolmaakt is waarmee distorsie stromingen resulteert. Het
concept van een negatieve weerstand wordt voorgesteld om die distorsie stromingen te annuleren.
Door middel van een eenvoudige intuïtieve analyse wordt de basiswerking van het techniek uitgelegd.
Door wiskundige analyse is de optimale waarde van de negatieve weerstand afgeleid en gerelateerd aan
de terugkoppeling theorie (Regel Techniek). Het is hier bewezen dat een eindige waarde van het
negatieve weerstand is nodig om de versterking loop van de terugkoppeling oneindig te maken
(theoretisch gezien), die praktisch door het verhogen van de OpAmp versterking niet gemakkelijk is. De
techniek wordt toegepast om de lineariteit van een RF IC chip ontvanger te verhogen. Een prototype is
geïmplementeerd in 65nm technologie. Meetresultaten tonen een toename van de IIP3 in-band van
9dBm naar > 20dBm, en IIP2 van 51 naar 61dBm, ten koste van een ruis toename van 6 naar 7.5dB en
een extra 10% vermogen gebruik. De chip bereikt een dynamic range (d.w.z. Spurious-Free Dynamic
Range) van 85dB in 1MHz.
IV
List of Main Abbreviations
RF: Radio Frequency
FCC: Federal Communications Commission
DTV and DVB-T: Digital TV broadcasting and Digital Video Broadcasting Terrestrial
DSA: Dynamic Spectrum Access
SR/SDR/CR: Software Radio/Software Defined Radio/Cognitive Radio
ADC/DAC: Analog-to-Digital Convertor/Digital-to-Analog Convertor
OFDM: Orthogonal Frequency Division Multiplexing
QPSK: Quadrature Phase Shift Keyed
SAW: Surface Acoustic Wave
NF: Noise Figure of the receiver
IIP3: 3rd order Input Intercept Point
IIP2: 2nd order Input Intercept Point
IM3: 3rd order InterModulation product
IM2: 2nd order InterModulation product
XM3: cross-modulation product
DPrA: Distortion Prediction Algorithm
LPF/ BPF: Low Pass Filter/ Band Pass Filter
X/ ⊗ : mathematical operation of Multiplication operation/ Convolution operation
OpAmp: Operation Amplifier
VGND: Virtual Ground of OpAmp
V
Contents
Abstract ................................................................................................................................................... I
Samenvatting ......................................................................................................................................... III
List of Main Abbreviations ....................................................................................................................... V
Chapter 1. Introduction ........................................................................................................................... 1
1.1 Cognitive Radio Motivation ............................................................................................................ 2
1.2 CR receiver hardware imperfections .............................................................................................. 5
1.3 Research Scope summary and questions ........................................................................................ 6
1.4 Outline of the Thesis ...................................................................................................................... 6
Chapter 2. Cognitive Radio Receiver Linearity Requirements [18] ............................................................ 9
2.1 Problem Definition....................................................................................................................... 10
2.2 Behavior Model for Nonlinearity .................................................................................................. 13
2.3 Analysis of the Spectral Locations of Distortion Products ............................................................. 13
2.4 Estimation of Linearity Requirements .......................................................................................... 20
2.5 Measurement Results .................................................................................................................. 24
2.6 Conclusion ................................................................................................................................... 26
Chapter 3. Distortion Prediction Algorithm and Frequency Selection ..................................................... 27
3.1 Distortion Prediction Algorithm (DPrA) ........................................................................................ 28
3.2 DPrA Application Examples .......................................................................................................... 43
3.3 Conclusions.................................................................................................................................. 48
Chapter 4. Analysis of a Very Linear Front-End [41] ............................................................................... 49
4.1 Motivation and Introduction ........................................................................................................ 50
4.2 Linearization Concept Analysis ..................................................................................................... 52
4.3 Mathematical Analysis ................................................................................................................. 56
4.4 Feedback Control Analysis ........................................................................................................... 58
VI
4.5 Stability Analysis .......................................................................................................................... 59
4.6 Conclusion ................................................................................................................................... 61
Chapter 5. Chip Design: UTCRxVLi [41]................................................................................................... 63
5.1 Receiver Implementation ............................................................................................................. 64
5.2 Noise Figure................................................................................................................................. 68
5.3 Measurement Results .................................................................................................................. 69
5.4 Benchmarking .............................................................................................................................. 74
5.5 Conclusions.................................................................................................................................. 74
Chapter 6. Conclusions and Recommendations ..................................................................................... 75
6.1 Summary and Conclusions ........................................................................................................... 76
6.2 Original Contributions .................................................................................................................. 80
6.3 Recommendations for Future Work ............................................................................................. 80
Appendix A.
Digital TV broadcasting (DTV) ......................................................................................... 83
Appendix B.
Wideband Correction Factor (∆Wideband) Effect on CR Linearity Requirement .............. 84
Appendix C.
DPrA Maple Program Implementation ............................................................................ 85
Appendix D.
Nonlinear Transimpedance Amplifier Analysis ................................................................ 87
Appendix E.
Latch-Up Analysis ........................................................................................................... 90
References ............................................................................................................................................ 93
List of Publications................................................................................................................................. 99
Acknowledgments ............................................................................................................................... 101
Biography ............................................................................................................................................ 103
VII
VIII
Chapter 1. Introduction
This chapter introduces the motivations and the main research questions regarding the importance of
distortion mitigation in Cognitive Radio. Section 1.1 begins with Cognitive Radio motivation to solve the
Frequency spectral utilization efficiency via Dynamic Spectrum Access. Afterwards in section 1.2, the
Cognitive Radio receiver hardware imperfections are discussed. Then the research scope summary is
presented in section 1.3 with the main research question. The chapter ends with the outline of the
thesis in section 1.4.
1
1.1 Cognitive Radio Motivation
The wireless communication revolution has led to a tremendous amount of applications and services in
our everyday life. Examples are the early success of GSM cellular telephony, DECT telephones at home,
Bluetooth technology for direct data exchange between devices, Wi-Fi internet access for local area high
data rate radio operation, GPS navigation, and currently the rapid expansion of wireless video and cloud
service access. To each application, government bodies allocate a specific radio frequency range in the
spectrum, called frequency band. Each band is derived by specific standards. Figure 1.1 shows the
frequency allocation made by the National Telecommunications and Information Administration (NTIA
[1]) in the USA. The Radio Frequency (RF) spectrum between 3kHz – 300GHz is divided in band(s) per
application, assigned through administrative licensing.
Figure 1.1: Spectrum allocation between 3kHz and 300GHz in the USA (NTIA [1])
Given the crowdedness of the spectrum plot of Figure 1.1 alongside with the observation that the
demand for wireless communication is growing exponentially, shortage of spectrum is foreseen, and in
this regard the chairman of the Federal Communications Commission (FCC) remarked in 2010 [2]: “Our
data shows there is a looming crisis. We may not run out of spectrum tomorrow or next month, but it’s
coming and we need to do something now”. In IEEE Spectrum [3], this problem is named “The great
spectrum famine”.
Different independent measurement campaigns (e.g. [4], [5]) found that most frequency bands were not
occupied continuously in time, i.e. spectral utilization efficiency is often poor. Although cellular network
bands are used intensively in most parts of the world, this is not true for many other frequency bands.
Measurement shows that less than 20% of the RF spectrum is actively used at any given time and place
[6]. Figure 1.2 presents an example of a spectrum measurement [4] for a part of the Digital TV (DTV)
2
broadcasting (Appendix A): channels 52 – 69, where each channel has 6MHz bandwidth, at the
Republican National Convention, New York City, USA.
Figure 1.2: Spectrum occupancy measurements [4], Location: Republican National Convention, New York City, USA
In an effort to improve the spectral utilization efficiency, FCC allows since 2008 [7] unlicensed (i.e.
secondary) users to opportunistically use TV broadcast channels that are licensed for TV broadcasters
(i.e. primary users), but are temporally locally unused.
Such free channels or bands are referred to as “white
spaces” [7] or “spectrum holes” [8]. Free channel can
for instance be utilized as shown in Figure 1.3,
referred to as “Dynamic Spectrum Access (DSA)” [9],
given the dynamic nature of this spectrum use. More
efficient use of such DTV frequency bands will be a
Figure 1.3: Opportunistic spectrum access [9]
key target of this thesis.
The DSA concept of operation poses several requirements on the radio hardware:
1. Flexibility and reconfigurability of the radio hardware to allow for dynamic switching between
different white spots across the frequency spectrum.
2. Supporting different standards, specified for different unoccupied primary users spots.
3. Sensing spectrum to detect the white spot instantaneously at any location.
Software (Defined) Radio towards Cognitive Radio:
Software Radio (SR) may achieve the ultimate level of flexibility, reconfigurability and can support all
different standards. A general schematic of a SR transceiver is shown in Figure 1.4. The SR involves
implementing all of the radio functionality in the software of the digital back-end (i.e. Digital Signal
3
Processing DSP) as described by Mitola [10]. The propagating Radio Frequency (RF) signal, with
frequency fC, is analog. This means that both the transmitted and received RF signals, via the antenna,
are also analog. Therefore a Digital-to-Analog Convertor
(DAC) in the transmitter and an Analog-to-Digital Convertor
(ADC) in the receiver are required. The focus of this thesis is
on the hardware requirements of the receiver side. The
transmitter part was the scope of another partner Ph D
project [11] inside the same research program at the
University of Twente. A true SR receiver is hard to
implement, because of the tough requirements on the ADC
resolution vs power consumption for mobile wireless
Figure 1.4: Software Radio receiver front-end
applications [12].
A solution that is more feasible is a radio receiver where part of the flexibility and configurability is
achieved by flexible analogue hardware instead of by software as shown in Figure 1.5, which is called
the Software Defined Radio (SDR) receiver.
The RF band filter, the channel filter and the
clock frequency to the mixer, which provide
the frequency down shift (or conversion)
functionality,
are
all
flexible
and
reconfigurable. The bandwidth of the RF
Amplifier can be wide to cover all the
Figure 1.5: Software Defined Radio receiver front-end
required bands, on condition that the RF filter
has enough out of band interferer rejection, so that the interferer will not distort the desired signal.
Extending the SDR and combining it with the Spectrum Sensing provides a radio that is aware of its
frequency environment and autonomously adjusts its communication parameters to fully realize DSA
(see Figure 1.3). This concept is called a Cognitive Radio (CR) [13]. The concept of CR was first proposed
by Joseph Mitola in a seminar at the Royal Institute of Technology in Stockholm (KTH) in 1998 and
published in an article by Mitola and Gerald Q. Maguire, Jr. in 1999 [14] .
In summary, the CR concept promises to increase the spectrum utilization efficiency via DSA. The
hardware physical implementation of the CR is based on SDR receiver plus Spectrum Sensing device. In
the next section, the SDR/CR receiver hardware imperfections especially the heavy linearity
4
requirements will be explained. This thesis aims to mitigate the linearity requirements challenges by
benefiting from the presence of spectrum sensing device, which is assumed to be on board.
1.2 CR receiver hardware imperfections
Imperfections of practical receivers decrease the quality of the received desired signal, via distortion and
noise [15].
The DTV interferers will often lead to distortion in practical receivers with a limited degree of up-front
filtering and linearity, which will be explained in Chapter 2. In traditional narrowband RF receivers, out
of band strong interferers can greatly suppressed by a Surface Acoustic Wave (SAW) RF band filter (a
fixed center frequency version of the RF band filter shown in Figure 1.5) after the antenna. This will
greatly reduce the distortion level of the received signal at the ADC input, i.e. the required out of band
linearity of the receiver will be also reduced. The SAW RF filter with a fixed center frequency does not
help to immune the CR receiver from the DTV interferers that actually exists in-band, hence the in-band
linearity becomes a critical specification for CR receiver. The ultimate solution is a flexible and
reconfigurable RF channel filter with high rejection (i.e. very high order filtering) of the interferers from
the adjacent channel(s). Note that guard channels are not desired to maximize the spectral utilization
efficiency. Such RF channel filters are very difficult to be implemented, especially on-chip. Different
implementations of flexible and reconfigurable RF filters are presented as example [16] and [17].
However their application as flexible and reconfigurable RF channel filter in the DTV band is bounded by
their limited adjacent channel(s) interference rejection ratio. As a consequence of the poor rejection
performance of RF channel filter, the in-band linearity of the receiver must be increased to tolerate the
received DTV interferers that actually exists in-band.
Any receiver contributes noise. The noise contribution, which is quantified by the Noise Figure (NF),
increases the noise floor. Hence the sensitivity level of the receiver that must be above the noise floor
by the minimum SNR requirement will decrease. Reducing the radio range (i.e. cell size) and/or
increasing the transmitted signal power helps always in reducing the NF requirement of the receiver,
note that this is not the case with linearity requirement, hence distortion becomes more important and
noise less, especially for the near future map of communication, which develops towards small cells.
5
The following section summarizes the research scope and presents the main research questions of the
thesis.
1.3 Research Scope summary and questions
In summary, a CR requires flexibility in receive frequency and such flexibility comes at the cost of less
up-front filtering, hence distortion products will be increased, especially because the in-band linearity is
normally limited. The high distortion products generated across the unoccupied DTV channels will finally
limit the spectral utilization efficiency. Hence distortion mitigation for CR is important regarding the
main idea of increasing the spectrum usage via DSA in the DTV spectrum band. It will be assumed that a
spectrum sensing device is available on-board, and options will be explored to benefit from the
presence of such a device. The main research questions are as follows:
1. How to analyze the distortion level (i.e. linearity requirement) across the unoccupied DTV band
for Dynamic Spectrum Access CR radio receiver?
2. How to predict the distortion level in Dynamic Spectrum Access, given Spectrum Sensing data?
3. How to increase the linearity of CR receiver hardware, given less up-front filtering?
Finally, the CR concept of operation demands both, (RF) analog and digital signal processing that
continuously interacts with each other. This makes CMOS technology to be the chosen candidate in this
thesis for possible monolithic integration between those two domains. In the following section, the
outline of the thesis will be discussed.
1.4 Outline of the Thesis
In this section, the outline of the thesis is presented in addressing the main research questions:
Chapter 2 investigates the first question:
Due to the existence of in-band strong DTV interferers and the nonlinearity of a radio receiver, different
types of distortion products will be generated inside a CR-receiver, namely intermodulation products,
cross-modulation products and self-interference products as analyzed in chapter 2. Using analysis based
on a memory-less nonlinearity model, it will be shown that CR receiver is always limited by crossmodulation and self-interference products in any white spots. While the level and the spectral location
6
of the intermodulation products across the white spots depends on the level and the spectral location of
the DTV signals. The cross-modulation and self-interference products are typically much weaker than
the intermodulation products. Thus it makes sense to monitor the level and spectral location of
interferes via spectrum sensing device, and classify the white spots into two types, namely
“intermodulation spots” and “intermodulation-free spots”. For the latter, the linearity requirement of a
CR receiver is more relaxed compared to the intermodulation spots. The analysis will be verified by
measurements.
Chapter 3 investigates the second question:
Based on the analysis in chapter 2, chapter 3 describes the development of a Distortion Prediction
Algorithm (DPrA) to quantify the level and the spectral location of the distortion products across the
white spots for Dynamic Spectrum Access. Then, according to the DPrA output, the CR receiver linearity
requirement can be quantified for different white spots as shown in chapter 3. Hence, the CR receiver
autonomously can select the most suitable white spots (i.e. unoccupied channel), taking into account
knowledge of linearity performance of its own CR receiver. Note, CR receiver with better linearity
performance can exploit more of the available white spots, hence higher spectral utilization efficiency
(i.e. the original motivation behind DSA). Therefore, increasing the CR receiver linearity is an important
target in this research.
Chapter 4 and 5 investigate the third question:
Chapter 4 focuses on increasing the linearity of CR receiver hardware. High linearity CMOS radio
receivers often exploit linear V-I conversion at RF, followed by passive down-mixing and an OpAmpbased Transimpedance Amplifier at baseband. Due to nonlinearity and finite gain in the OpAmp, virtual
ground is imperfect, inducing distortion currents. A negative conductance concept will be introduced
that allows for cancelling such distortion currents. Through a simple intuitive analysis, the basic
operation of the technique will be explained. By mathematical analysis the optimum negative
conductance value is derived and related to feedback theory. In- and out-of-band linearity, stability are
also analyzed.
The technique is applied to linearize the RF receiver in chapter 5, and a prototype is implemented in
65nm technology. Measurements show an increase of in-band linearity, quantified by IIP3 and IIP2, with
IIP3 from 9dBm to >20dBm, and IIP2 from 51 to 61dBm, at the cost of increasing the noise figure from 6
7
to 7.5dB and <10% power penalty. In 1MHz bandwidth, a Spurious-Free Dynamic Range of 85dB is
achieved at <27mA up to 2GHz for 1.2V supply voltage.
Chapter 6 summarizes the main research work:
Finally, Chapter 6 summarizes the main conclusion, makes a list of the main contributions of this work,
and provides recommendations for future work.
8
Chapter 2. Cognitive Radio Receiver Linearity Requirements [18]
This chapter investigates the first research questions “How to analyze the distortion level (i.e. linearity
requirement) across the unoccupied DTV band for Dynamic Spectrum Access CR radio receiver?”
The DSA operation inside the DTV band means that the power level and spectral location of the DTV
signals around the CR signal may instantaneously change; hence the traditional linearity analysis [15]
cannot be directly applied. The analysis of this chapter decomposes the complex view of the distortion
products into basic roots/types. Based on that, the linearity requirements for CR receiver are derived.
The content of this chapter is structured as follows. Section 2.1 introduces the problem. Then section
2.2 describes the receiver model and nonlinearity model. The spectral location analysis of the 3rd order
distortion products is done in Section 2.3. The effect of this spectral location of the distortion on the
linearity requirement is analyzed in Section 2.4. Section 2.5 verifies our theoretical analysis and presents
results from practical measurements. Finally, conclusions will be drawn in Section 2.6.
This chapter is a literal copy of the paper published in [18] with corrections related to Table 2.2 .
9
2.1 Problem Definition
Cognitive Radio is a new emerging radio communication paradigm [14] aiming at improving the
utilization efficiency of the scarce spectral resources. It senses unused “white spots” in the radio
spectrum that is licensed to a primary user and adapts its communication strategy to use these parts
while minimizing interference to the primary service.
This chapter studies requirements on the 3rd order linearity of a CR receiver for wireless applications
operating in the DTV bands [7]. Cognitive Radio requires high programmability for radio transmission
and reception because the position of the white spots changes dynamically with time (i.e. dynamic
spectrum). Traditional radio receivers are narrowband and are typically highly dedicated to a specific RFband, especially due to the fixed RF-filter (see Figure 2.1).
Figure 2.1: Traditional narrow-band RF receiver with one RF “Front-End”
The high out-of-band rejection of the RF-filter, usually a high-linearity highly selective Surface Acoustic
Wave (SAW) filter, suppresses the interference of out-of-band interferers. Traditional radio standards
define in detail how the radio-band is used, e.g. by specifying a multiple access method (TDMA, FDMA,
CDMA) and blocker profiles. Consequently in-band signals are “under control”, while RF-filtering reduces
the out-of-band interference. In such narrowband systems, typically third order intermodulation
products (IM3) produced by nonlinearities in the receiver dominate as is shown in Figure 2.2. Both the
signal and the IM3 product will be converted to the baseband and their ratio must be higher than the
minimum Signal to Distortion Ratio that is required to demodulate the transmitted information with
acceptable bit error rate.
10
Figure 2.2: Corruption of the desired signal due to intermodulation products caused by two interference signals
(Conventional narrow-band Radio)
For CR, the distinction between “out-of-band” and “in-band” signals more or less vanishes, as bands are
no longer reserved for one or a few radio standards. Also, the power of radio signals is not restricted to
a low level, as in ISM bands, in which more freedom is allowed provided low power levels are used
(typically 10-100mW). If we want to maximize spectral efficiency in the TV-bands, we like to use white
spots between DTV-signals and preferably even in the proximity of DTV signals. In such cases, very steep
high-order RF-filters would be needed, which are difficult to implement even with high-Q SAW filters. If
the frequency-distance between DTV signals is small, narrowband RF-filters would be needed. To cover
a significant part of the TV bands, many filters would be needed. From these observations we conclude
that a broadband RF receiver with relaxed RF-filter requirements is highly desired. In this chapter we
look for ways to benefit from spectral sensing which is available in a CR anyway. We will look at a
broadband receiver, but still limit the bandwidth to below an octave so that second order distortion
products are largely suppressed by the RF-filter and the 3rd order distortion is the main remaining
problem.
A broadband RF receiver implies the reception of undesired interference as shown in Figure 2.3. The
example of a DTV spectrum in Figure 2.3 consists of 20 numbered frequency spots: 5 spots contain
interference (i.e. DTV signals at spots 1, 4, 11, 16 and 19) and one spot contains the desired signal (i.e.
CR signal at spot 8).
11
Figure 2.3: An DTV instantaneous spectrum scenario at a wideband RF bandpass filter output with 20 frequency spots
In [19], the use of spectral sensing to relax linearity requirements of a receiver is addressed. The paper
assumes a dense spectrum with a large number of interferers. According to the law of large numbers,
the distribution of the intermodulation products due to the many interferers becomes even across the
spectrum and can be seen as a distortion induced (noise) floor. However, a different picture arises when
the spectrum is dominated by a few high power interference signals in the band of the RF bandpass
filter. According to our analysis and measurements, different spectral distributions of the interference
produce quite different distributions of the distortion products across the spots of the DTV spectrum.
Consequently, the linearity requirements of the CR receiver will be different for each white spot. A
prediction of the spectral location of the distortions provides the CR with the ability to choose an
appropriate white spot. The required information about the level and the spectral location of the
interference signals can be obtained from spectrum sensing, which is normally available anyhow in a CR.
Our analysis is based on a CR receiver with direct conversion architecture. The theory can be extended
to other (non-zero) IF architectures, where the image problem must be included as an extra contribution
of distortion to the desired signal, but this is outside of the scope of this chapter.
12
2.2 Behavior Model for Nonlinearity
It is common to model the nonlinearity of an RF receiver assuming a memory-less weakly nonlinear
model (e.g. [15] and [20]) with a truncated 3rd order Taylor series around the DC operating point:
s out = k 1s in + k 2 s in2 + k 3s 3in
Equation 2.1
where k1 specifies the gain, while k2 and k3 characterize the second and third order nonlinearity of the
circuit with input signal sin and output sout (s can be both a voltage or current). The sin contains all the
input signals at the output of the RF bandpass filter.
Equation 2.1 contains three terms. The first term with k1 is the amplified version of sin, the linear term.
The second and the third term will be called the 2nd order and 3rd order nonlinear terms, respectively.
Those nonlinear terms generate different distortion products. Some of those distortion products fall into
the band of interest. By using an octave RF bandpass filter [19] and RF components with differential
structure, the effect of 2nd order nonlinear term can often be handled. However, the effect of the 3rd
order nonlinear terms cannot be ignored because their distortion products always fall inside the desired
DTV spectrum.
2.3 Analysis of the Spectral Locations of Distortion Products
A CR scans the spectrum for white spots, which have a bandwidth in the order of 6-8 MHz in the DTV
spectrum. Although the DTV signals and the CR signal are wideband (i.e. the CR occupies the whole spot)
signals, The analysis starts representing them as tones with a power equal to the integral of the power
of the wideband signal. Afterwards, the analysis will be extended to deal with the spectral shape of the
wideband signals.
Referring to Figure 2.4, sin in Equation 2.1 contains three DTV signals around ω1, ω6 and ω19. To not lose
the general view, CR signal will be introduced later after specifying the spectral location of the distortion
products of the DTV signals.
13
Figure 2.4: An instantaneous scenario of the DTV spectrum dominated by three interference signals in the RF passband
The sin can be described by the following equation:


s in = Re 2 ∑ A TVn e jωn t 
 n =1,6,19

Equation 2.2
where ATV is the corresponding RMS value of the DTV (single tone) signals. Substituting sin in the 3rd
order nonlinear term of Equation 2.1 gives the following distortion products:
 3

3
3
2
2
k 3sin = Re 2  k 3A TV1 + 3k 3A TV1 A TV6 + 3k 3A TV1 A TV19 e jω1 t 

 2

 3

3
2
2
+ Re 2  k 3A TV6 + 3k 3A TV6 A TV1 + 3k 3A TV6 A TV19 e jω6 t 

 2

 3

2
+ Re 2  k 3A TV1 A TV6 e jω11 t 

 2

{ (
)
+ Re 2 3k 3A TV1 A TV6 A TV19 e jω14 t
Equation 2.3
}
 3

3
2
2
+ Re 2  k 3A TV19 + 3k 3A TV19 A TV1 + 3k 3A TV19 A TV6 e jω19 t 
2




Figure 2.5 provides a visual way to understand the spectral location of the distortion products that are
double underlined in Equation 2.3. These products are important because they fill the white spots,
where a CR may operate. They represent two types of intermodulation products:
1. The intermodulation products of two interferers (observe spot 11), which is the traditional IM3
that is already mentioned in literature [15].
2. The intermodulation products of three interferences (spot 14), which is 6 dB higher than the
traditional 2-tone IM3 for interferes with equal power. This type of distortion will also be called
IM3 in this paper.
14
Figure 2.5: Spectral location of the intermodulation products for the case of Figure 2.4
For simplicity, we will keep the power of the DTV signals equal in this paper. Of course, using Equation
2.3 case with different power can be analyzed easily.
Repeating the same analysis for another DTV spectrum scenario with three large interference signals at
different locations, gives the distortion distribution that is depicted in Figure 2.6.
Figure 2.6: Spectral location of the intermodulation products for another case
Comparing Figure 2.5 with Figure 2.6, we recognize that the first figure presents 15 white spots without
intermodulation distortion (i.e. IM3-free spots) and 2 white spots with IM3 (i.e. IM3-spots), while the
second figure presents 13 IM3-free spots and 4 IM3-spots. For all the possible spectral locations of the
three interference signals, the number of the IM3-free spots can be counted and a histogram of it can be
plotted as shown in Figure 2.7. The expected value is 13 white spots here, with a spread from 8 to 17.
15
Figure 2.7: Probability density of IM3-free spots for all possible scenarios with 3 interferers (Figure 2.6 depict one scenario)
Note that the fact that a spot is IM3-free does not mean that there is no 3rd order nonlinearity
requirement anymore for the receiver. Actually, the interferers will cause cross-modulation on top of
the desired signal as shown in Figure 2.8.
Figure 2.8: Cross-modulation products (XM3) due to interferers which cross-modulate the CR signal
To analyze the effect of cross-modulation, the mathematical expression for sin including a CR at spot 17
is written as follows:


s in = Re  2A CR e jω17 t + 2 ∑ A TVn e jω n t 
n =1,6,19


Equation 2.4
where ACR is the corresponding RMS value of the CR signal. Substituting sin in Equation 2.1 and retaining
just the terms that exist in the desired band gives:
16






 3

3
2
2
2 
Sout Desired Band = Re
2k1ACR + 2k3  ACR + 3ACR ATV + ATV + ATV  e jω
424
3
14444
4244444
3 
2 23
Desired1
 1
amplifiedCR signal
Cross-modulation
Self
interferen
ce


14444444442444444444
3


3 order distortion products
(
1
6
19
)
17
t
Equation 2.5
rd
where the first term on the right hand is the amplified version of CR signal, the second term is the selfinterference product [21] and the last three terms are the mentioned 3rd order cross-modulation [22]
products (XM3) of the DTV signals. The self-interference product can usually be neglected in comparison
to the XM3 products because the DTV signals are usually much stronger than the desired signal. The
same analysis is done to show the spectral location of the distortion products of 4 interference signals
with high power and the results can be seen in Figure 2.9 and Figure 2.10 (Expected value = 7).
Figure 2.9: Spectral location of the intermodulation products for a scenario with 4 interferers
Figure 2.10: Probability density of the existence of given number of free spots without intermodulation distortions for
scenarios with 4 interferers
17
As mention previously, both the DTV signals and the CR signal are wideband signals with about the same
bandwidth. Therefore their shape can be approximated by square blocks in the frequency domain. The
mathematical representation in the time domain of those signals and their distortion products is written
in the same way as mentioned previously with the exception that ATV and ACR now represent the
complex envelope of the DTV signal(s) and CR signal, respectively. In the time domain, the
intermodulation and the cross-modulation products are nothing else than a multiplication of three
signals, where the involved signals are indicated with arrows. E.g. three arrows constitute the inputs to
the distortion circles in Figure 2.5, Figure 2.6 and Figure 2.8, in some cases from two frequencies, in
other cases from three frequencies. In the frequency domain, this corresponds to two times convolution
of the signals. The result of this convolution is shown in Figure 2.11 for the interference scenario of
Figure 2.5, but now with more realistic block-shaped spectra instead of single tones.
Figure 2.11: Bell shapes for the distortion products of three wideband interferers
18
A Bell shape results, occupying three spots as shown in Figure 2.11. Analysis shows that more than 66%
of the distortion power is concentrated in the frequency spot in the center of the Bell, while the power
spilled over to the adjacent frequency spots is less than 17% of the distortion power each. Referring to
Equation 2.3, the bell shapes around spots 1, 6 and 19, each contain a self-interference product and two
cross-modulation products (e.g. the Bell shape around spot 6 contains its self-interference product and
the cross-modulation products of the other two interferers of spots 1 and 19 on top of the interferer of
spot 6). As a consequence of spreading, the distortion to the adjacent spots of 1, 6 and 19, the spots 2,
5, 7, 18 and 20 will be called adjacent-distortion spots (i.e. the distortion that is adjacent to the
interferers). The linearity requirements in those adjacent-distortion spots of Figure 2.11 will be higher
than that of Figure 2.5. Additionally, the bell shape of the IM3 products (i.e. around spot 11 and 14) will
reduce the number of IM3-free spots by 4 spots in comparison to Figure 2.5, as spots 10, 12, 13 and 15
will contain 17% of the total IM3 power of the Bell shape, hence those spots will be also called IM3 spots.
Paper [23] presents an analysis about the interference effects in Multi-Band - Orthogonal Frequency
Division Multiplexing (MB-OFDM) systems. It proves that the power of the IM3 resulting from two OFDM
signals (Quadrature Phase Shift Keyed (QPSK) modulation) each with a large number of sub-carriers is
higher than that of two tones that have the same power as the OFDM signals:
+
PIM = PIM + ∆Wideband [dBm]
3
Equation 2.6
3
where PIM3+ and PIM3 represent the intermodulation distortion of two wideband signals and two tones,
respectively. PIM3+ is distributed over three spots (e.g. see spots 10, 11, 12 of Figure 2.11). ∆Wideband
represents a correction term. Instead of dealing with wideband complex signals, Equation 2.6 provides
the possibility to analyze the nonlinearity with the traditional two-tone test method and simply add a
correction term.
As an example that can be applied to our case, a MatLab simulation setup has been built to simulate the
nonlinearity with two real OFDM signals at RF. QPSK modulation was used and the frequency difference
between the sub-channels is 4.46 kHz satisfying the Digital Video Broadcasting Terrestrial (DVB-T)
standard. The resulting IM3 was simulated. The simulation result is shown in Figure 2.12, where the
number of the sub-carriers per OFDM signal has gradually been increased from 1 to 900. Thus we can
compare the IM3 generated for two (narrowband) single tones with that for a truly wideband OFDM
signal with many sub-carriers, where the total signal power in all cases is kept constant.
19
Figure 2.12: Comparison of the IM3 power caused by two (single) tones and two multi-tone OFDM signals with a number
of sub-carriers from 1 to 900
Figure 2.12 shows a ∆Wideband of about 3 dB. Changing the modulation of the two OFDM signals to
16QAM, we find 3.5 dB. By further increasing the number of the QAM modulation to higher orders, the
signal characteristic of the two OFDM signals resembles the case of having two band limited noise
signals. In this case the factor is increased to 4.5 dB. Table 2.1 summarizes the simulation results:
Table 2.1: CORRECTION FACTOR ∆WIDEBAND FOR TWO OFDM SIGNALS WITH DIFFERENT MODULATION CHARACTERISTICS
Two Signals
Two OFDM: QPSK
Two OFDM: QAM16
Two Band limited Noise
Max. ∆Wideband [dB]
3
3.5
4.5
As the purpose of this paper is to compute the linearity requirements for CR, the highest value of the
maximum ∆Wideband will be taken, which is 4.5 dB.
2.4 Estimation of Linearity Requirements
The linearity of a receiver must be adequate to preserve a minimum Signal to Distortion Ratio needed at
the demodulator. As previously depicted in Figure 2.5, we distinguish two types of white spots: IM3-free
spots and IM3-spots. We will show that the CR requires higher linearity if it operates at IM3-spots than
for IM3-free spots with only self-interference and the XM3. This section begins with introducing an
estimation of the linearity requirement for the case that the CR is operating at an IM3-free spot and the
interferers are modeled by tones. Afterwards the linearity equation will be extended to cover the case
of wideband interference signals. Finally, this linearity requirement will be compared to that of a CR
operating at an IM3-spot and at an adjacent-distortion spot.
20
Referring to Figure 2.8 and Equation 2.5, the Signal to Distortion Ratio can be written as follows:
S
(k 1 A CR )
=
2
D 


 3k 3 A CR  ∑ A TV 2  
n



 n =1,6,19


2
Equation 2.7
where the self-interference has been neglected. The linearity requirement can be presented as the ratio
between k1 and k3:
k1
k3
= 3×
spot 17
S

D 
∑A
n =1,6,19
2
TV n



Equation 2.8
where S/D is the required Signal to Distortion Ratio for the CR receiver output. It is proportional to the
sum of the power of the DTV signals, which will cross-modulate the receiver CR-signal.
Equation 2.6 presents the correction factor ∆Wideband for the intermodulation power. The same ∆Wideband is
valid for the XM3 power (the ∆Wideband is just the relative difference between the modulation products of
the two tones in relation to that of the two wideband signals). Therefore Equation 2.6 provides a useful
tool to extend the derivation of Equation 2.8 to deal with wideband interference signals as explained in
Appendix B:
k1
k3
= 3×
spot 17
∆
S
× 10
D
Wideband
× 0.66
10
∑A
n =1,6,19
2
Equation 2.9
TVn
Using the mentioned maximum ∆Wideband of 4.5dB gives:
k1
k3
= 3×
spot 17
S
2
× 1.4 ∑ A TV
D
n =1,6,19
Equation 2.10
n
It is instructive to relate Equation 2.10 to the traditional IIP3 specification. According to the definition of
IIP3, the RMS value of IIP3 is related to k1/k3 as follows [15]:
AIIP =
3
4 k1
3 k3
2
Equation 2.11
4 k1
=
6 k3
[V]
21
Substituting Equation 2.10 in Equation 2.11 gives the following equation:
PIIP
3
spot 17
= 2×
S
× 1.4 ∑ PTV
D
n =1,6,19
n
Equation 2.12
[W]
We recognize that dealing with wideband signals instead of tones, the IIP3 requirement will be increased
by around 1.5 dB (i.e. 10*log(1.4)). To get a sense of what Equation 2.12 presents, let’s assume that the
interference signals of Figure 2.11 have a power of -10 dBm and the required Signal to Distortion Ratio is
10 dB. Then the required linearity to keep XM3 low enough expressed in terms of IIP3 would be 4 dBm.
Increasing the power of the interferences by 10 dB increases the requirement also by 10 dB. Equation
2.12 can be generalized as follows:
PIIP = 2 ×
3
S
× 1.4
D
∑P
TV n
[W]
Equation 2.13
n
where n refers to all the interference signals that exist in the RF bandpass filter of the CR receiver.
As a final step, let us compare IIP3 of Equation 2.12 (i.e. IIP3XM3) to that IIP3 requirement when the CR
operates in IM3-spot (IIP3IM3) like for example spot 11 in Figure 2.11. In this case the distortion products
will contain the IM3 plus the XM3 products:













3

2
2 
s out Desired Band = Re
2k1ACR
+ 2k 3  3ACR ∑ A TVn + ATV1 ATV6  e jω11 t 
1424
3


2 4243 

n =1,6,19
 1442443 1Intermodul

Desired amplified CR signal

ation
- modulation
 4Cross
144
44
4
424444444
3



3rd order distortion products


Equation 2.14
Writing down the Signal to Distortion Ratio gives:
(k1A CR )
S
=
2
D 




 3k 3A CR  ∑ A TV 2  +  3 k 3A TV A TV 2   
n
1
6

 2


n =1,6,19
2
4443 1442443  
 1444
Intermodul
ation

Cross- modulation

2
Equation 2.15
The resulting linearity requirement in terms of k1/k3 becomes:
k1
k3
= 3×
spot 11
S 
D 
∑ A TVn +
2
n =1,6,19
1 A TV1 A TV6
2
A CR
2




22
Equation 2.16
Expressing the linearity requirements in IIP3 and extending the derivation to wideband signals gives the
following two equations for spot 11 and 10, 12 (see Figure 2.11):


∆ Wdeband
∆ Wdeband
PTV1
S 
1

10
10
PIIP3
= 2×
10 424×4
0.66
0.66
4
3 ∑ PTVn + 2 110
44
424×4
4
3 P PTV6  [W]
spot 11
D  144
n =1,6,19
CR
1.4
1.4


PIIP
3
spot 10,12
= 2×
S
D
∆
∆


PTV
1
10
10
×
0.66
P
+
10
×
0.17
P
 110
TV
442443 n =∑
443 P TV  [W]
2 1442
1,6,19
CR

1.4
0.7

Wideband
Wideband
Equation 2.17
Equation 2.18
1
n
6
By comparing between Equation 2.17 and Equation 2.18, the IIP3 requirement of spots 10 or 12 (i.e. the
spots adjacent to the Bell center spot that contains 17% of the total IM3 power) is 3 dB lower than the
IIP3 requirement of spot 11 (i.e. the center spot of the Bell shape that contains 66% of the total IM3
0.17 / 0.66 = 1 / 2 .
power), hence
As mentioned previously, the CR requires also linearity requirement in the adjacent-distortion spots (see
Figure 2.11). Repeating the same analysis as previously done and using Equation 2.3, Equation 2.5 and
Equation 2.11 results:
PIIP
3
spot 2
= 2×
S
D
∆
∆

PTV 
1
10
10
×
0.66
P
+
10
×4
0.17
110
∑
442443 n=1,6,19TV 14424
3 ( 2 PTV + PTV + PTV ) P  [W]
CR 
1.4
0.7


Wideband
Wideband
1
n
1
6
19
Equation 2.19
To get a sense about the linearity requirement of Equation 2.12, Equation 2.17, Equation 2.18 and
Equation 2.19 consider the same numbers as for the previous example, while assuming that the CR
signal is equal to -60 dBm. Then the different IIP3 requirements are summarized in Table 2.2.
Table 2.2: LINEARITY REQUIREMENT OF EQUATIONS (12), (17), (18) AND (19) FOR THE SCENARIO OF Figure 2.11
Equation
Equation 2.12
Spot no.
17
Spot Description
IM3-free spot
IIP3 [dBm]
4
Equation 2.17
11
IM3-spot
21
Equation 2.18
10, 12
IM3-spot
18
Equation 2.19
2
Adjacent-distortion spot
25
Clearly, there is a very significant benefit in looking for IM3-free spots. Such spots like spot 17 can well
make a key difference in the feasibility (e.g. [28] and [31]) of CR.
23
2.5 Measurement Results
The spectral location of the distortion products is the consequence of modeling a memory-less weakly
nonlinear device with a Taylor series as presented in Equation 2.1. The model has been verified by
building a setup where a Low Noise Amplifier (LNA) from the Mini-Circuits (i.e. ZFL-1000LN) is tested
with three interference signals. The power and the spectral location of the distortion products are
measured. The LNA has gain of 18 dB, IIP3 of -6 dBm, Output 1 dB compression point of 3 dBm and its
operating region is from 0.1 to 1000 MHz (i.e. VHF/UHF band). The interference signals are three carriers
with QPSK modulation, which are generated by vector signal generators from Agilent technology. In the
first experiment, the power of each interferer is around -26 dBm and their frequency location is around
the following values: 548, 563 and 600 MHz. This case looks like the case shown in Figure 2.5 and the
measurement result is shown in Figure 2.13.
Figure 2.13: Spectral location of the distortion products in case of three interferers at 548, 563 and 600 MHz
Comparing to Figure 2.5, Figure 2.13 indeed contains two distortion products:
1. The first interferer at 548 MHz and the second interferer at 563 MHz generates the traditional
IM3 at 578 MHz (i.e. 2*563-548)
2. The combination of all the three interferers generates the previously mentioned new IM3 at
frequency 585 MHz (i.e. 600-(563-548) and (600-563)+548). Shown in Figure 2.13, this type of
IM3 is higher than the first IM3 by 6 dB as our theory has predicted (see Equation 2.3).
24
The power of those IM3 products are measured and compared to the results of our equations as shown
in Table 2.3.
Table 2.3: Comparison between the measurement and the theory results of the IM3 distortion products
IM3 frequency [MHz]
Theory [dBm]
578
585
-39.4
-33.4
Measurement
[dBm]
-38.7
-32.5
Another experiment has been done with a scenario that looks like the case shown in Figure 2.6. The
measurement result is shown in Figure 2.14.
Figure 2.14: Spectral location of the distortion products in case of three interferers at 554, 560 and 590 MHz
Consequently, the measurement results verify our theory about the effect of interference locations on
the spectral location of the distortions. The IM3 spots and the IM3-free spots are clearly shown in Figure
2.13 and Figure 2.14.
25
2.6 Conclusion
This chapter studied the 3rd order nonlinearity requirements for a CR receiver with a wideband RF
bandpass filter that operates in the DTV spectrum, where several high powers DTV signals fall in the RF
filter-band. Due to the nonlinearity of the RF receiver and the existence of strong interference signals,
different types of 3rd order distortion products will be produced. The spectral location of these distortion
products depends on the spectral locations of the interference signals. Analysis and measurements
show that white spots can be classified into two types:
1. IM3-spots, for which the distortion is dominated by IM3 products. The IM3 products are much
stronger than XM3 and self-interference if the interferers are much stronger than the desired
signal.
2. IM3-free spots where the CR signal will only be distorted by XM3 and self-interference.
If an IM3-free spot is selected for CR operation, its third order linearity is relaxed. In a scenario with
three interferers at -10dBm and a CR-signal at -60dBm, the resulting IIP3XM3=4dBm, while the IIP3IM3 is
18-25 dBm. Since the linearity of any radio design is limited, typically to not much more than 0dBm in
CMOS, spectral sensing information about the level and the spectral location of the strong interferers
can be used to predict the location of the IM3-free spots. Such a prediction mechanism seems to be
almost indispensable for a CR that can work in a wide RF-band. Next chapter is pursuing in this direction.
26
Chapter 3. Distortion Prediction Algorithm and Frequency Selection
This chapter investigates the second research question: “How to predict the distortion level in Dynamic
Spectrum Access, given Spectrum Sensing data?”
The coexistence of a Cognitive Radio signal among the primary users of the DTV spectrum increases the
linearity requirements of a wideband CR receiver considerably to very high numbers (see Table 2.2) as
derived in the previous chapter. The analysis in chapter 2 has classified the DTV white spots into IM3spots and IM3-free spots, depending on the existence of IM3 distortion products at a particular white
spot. It was concluded that the IIP3 requirements on the CR receiver will be different across the white
spots, and since the linearity of any radio design is limited, this means that not all the white spots are
equally appropriate for CR operation. This chapter aims at developing a Distortion Prediction Algorithm
(DPrA) that calculates the distortion level across the white spots. Based on the DPrA results in
combination with the theory in chapter 2, the linearity requirement per white spot can be quantified.
Section 3.1 presents the DPrA concept. Then, section 3.2 presents the verification of the DPrA concept
by applying it to the DSA1 scenarios of the previous chapter and some other application examples test its
ability to recognize the construction pattern of complex distortion products. Finally, the chapter ends
with conclusions in section 3.3.
1
Where the signals level across the bands/channels varies with respect to time for a given location
27
3.1 Distortion Prediction Algorithm (DPrA)
The analysis starts in representing the wideband1 DTV signals and the wideband CR signal (i.e. the CR
occupies the whole channel) as single frequency with a power equal to the integral of the wideband
signal power. Afterwards, the analysis will be extended to deal with the spectral shape of the wideband
signals.
The target of the DPrA is to derive the instantaneous level and spectral location of the 3rd order
distortion products across the white spots by using spectral sensing information. An example of the
instantaneous level and the spectral location of some DTV signals is shown in Figure 3.1. The calculation
of the output spectrum in this figure is already explained in chapter 2 (see Figure 2.5 and Equation 2.3).
rd
Figure 3.1: Example of the operation of the 3 order Distortion Prediction Algorithm (DPrA)
As indicated in Figure 3.1, the DPrA boils down to evaluate the following mathematical expression:
k 3s3in (t ) = k 3[sin (t )× sin (t )× sin (t )]
Equation 3.1
where k3 characterizes the third order nonlinearity term of the circuit. Note that sin contains all the input
signals passing through the RF bandpass filter. Equation 3.1 involves two times multiplication in the time
domain of sin. In this case, it consists of three signals at channel numbers: 1, 6 and 19, i.e.:


sin (t ) = Re 2 ∑ ATVn e jωn t 
 n =1,6,19

Equation 3.2
where n represents the corresponding channel number and ATV is the corresponding RMS value of the
DTV (single tone) signals.
1
Each wideband DTV signal occupies a channel bandwidth of 6-8 MHz depending on the different standards
28
Repeating the procedure of the previous chapter in substituting sin in Equation 3.1 produces the
following distortion products:
[
]

k 3 sin (t )× sin (t )× sin (t ) = Re


+ Re


+ Re


3

2  k 3ATV1 3 + 3k3ATV1 ATV6 2 + 3k3ATV1 ATV19 2 e jω1 t 

2


3

2  k 3ATV6 3 + 3k3ATV6 ATV1 2 + 3k3ATV6 ATV19 2 e jω6 t 
2



3

2  k 3ATV1 ATV6 2 e jω11 t 
2


{ (
)
+ Re 2 3k3ATV1 ATV6 ATV19 e
jω14 t
}
Equation 3.3
 3


+ Re 2  k 3ATV19 3 + 3k3ATV19 ATV1 2 + 3k3ATV19 ATV6 2 e jω19 t 

 2

Expanding the total number of the channels, and the number and spectral locations of the interferers
will increase the complexity of Equation 3.3. Thus it is not practical to implement this time domain
derivation, especially in case of a dynamic spectrum scenario that changes continuously. However, it will
be shown below that the analysis can also be done in the frequency domain in a systematic way using
matrix operations. In the frequency domain, the multiplications in Equation 3.1 correspond to
convolutions:
DPrA ⇒ Third Order Distortion Products in the Frequency Domain
Equation 3.4
⇒ [Sin (f ) ⊗ Sin (f ) ⊗ Sin (f )]
where Sin(f) is the frequency domain representations of sin(t).
To build an intuitive understanding of the matrix implementation of the DPrA, let’s start with just two
tones, as shown in Figure 3.2 and expressed as follows:


sin (t ) = Re 2 ∑ ATVn e jωn t  = ∑ Ame- jωm t + ∑ Ame+ jωm t
 n=2,3
 m=−2,-3
m=+2,+3
Equation 3.5
where m represents the channel number for double side spectrum (i.e. positive and negative
frequencies) and Am is the MAX value of the DTV (single tone) signals at channel m. For simplicity, let’s
assume that all the occupied channels contain equal signal amplitude, equal to A TV / 2 as follows:
A -3 = A - 2 = A + 2 = A + 3 = A TV / 2
Equation 3.6
29
Figure 3.2: Frequency spectral content of the of two tones (two single frequencies)
The frequency spectrum span of Sin(f) is equal to 7 channels (N=7 channels: m=-3,-2,-1,0,+1,+2,+3;)
around the zero-frequency (see Figure 3.2), each channel has a bandwidth equal to the DTV channel
bandwidth, namely fCH (for the clarity of later figures, the label fCh on x-axis will be removed). Assuming
pure tones, Sin(fm) can be written in the discrete form:
Sin (f m ) = [A + 2 δ (f m - 2 f ch ) + A − 2 δ (f m + 2 f ch )] + [A + 3δ (f m - 3 f ch ) + A − 3δ (f m + 3 f ch )]
=
2
2
A TV [δ (f m - 2 f ch ) + δ (f m + 2 f ch )] +
A TV [δ (f m - 3 f ch ) + δ (f m + 3 f ch )]
2
2
Equation 3.7
SCONV(fm) in Equation 3.8 presents the mathematical expression of the first convolution process between
the square brackets in Equation 3.4:
SCONV (f m ) = Sin (f m ) ⊗ Sin (f m )
=
=
∑ (S (ρ ) S (f
in
ρ m = − ∞ .. + ∞
m
in
ρ m = − ∞ .. + ∞
m
in
m
− ρ m ))
∑ (S (ρ ) S (m f
in
ch
Equation 3.8
− ρ m ))
Where ρm is a dummy variable required to perform the convolution operation [25]. This convolution
operation involves a multiplication between the Sin(ρm) spectrum and its flipped1 version, shifted to the
right by |m|fCH, written as Sin(m fCH - ρm). Figure 3.3 visualizes this convolution operation that calculates
the frequency content of SCONV(fm) at m=0..+6. For the negative frequency side of SCONV(fm) (i.e. at m=-6..1), a similar plot can be created, by multiplying the spectrum of Sin(ρm) with its flipped version, now
shifted to the left by |m|fCH. Figure 3.3 shows that seven shifting operations are required to calculate
the frequency content of SCONV(fm) at m=0..+6: one shifting operation for the frequency content at the
origin (i.e. at m=0) and another six shifting operations for the positive frequency content (i.e. at
m=+1..+6). To complete the spectrum content of SCONV(fm) with the negative frequency content (i.e. at
m=-6..-1), still extra 6 shifting operations are needed.
1
flipped version of the Sin(ρm) around the y-axis gives Sin( - ρm)
30
Figure 3.3: Visualizing the convolution operation of Equation 3.8
Hence for the total spectrum view, 13 shifting operation (i.e. in general form 2N-1, where N=7) from the
left to the right are required. The DPrA implements this convolution operation via Matrix multiplication
as shown in Figure 3.4, where the Qm elements of the SCONV (m f ch ) is derived. First, we form a vector,
called S in . This vector corresponds to Sin(fm) (see the 7 channels in Figure 3.2) and contains 7 elements.
The index and the value of those 7 elements correspond to the spectral location and the amplitude of
the signals that composite the Sin(fm) spectrum.
31
Figure 3.4: Matrix implementation of the convolution operation of Figure 3.3
The Sin(m fCH - ρm) is constructed in a matrix form as follows:
•
Its middle column (i.e. 7th column) is nothing else than the flipped version of S in shifted by zero
(i.e. Sin(- ρm)). Then multiplying S in elements by the elements of this column calculates the
signal SCONV (m f ch ) quantity at m=0.
•
The columns at the right of the 7th column (i.e. 8th - 13th columns) represent the flipped version
of S in , right shifted by m=+1 to +6, respectively (i.e. Sin(m fCH - ρm)). Multiplying S in with those
columns provides the signal SCONV (m f ch ) quantity at m=+1 to m=+6.
•
The columns at the left of the 7th column (i.e. 1st - 6th columns) represent the flipped version of
S in , left shifted by m=-6 to -1, respectively (i.e. Sin(m fCH - ρm)). Multiplying S in with those
columns provides the signal SCONV (m f ch ) quantity at m=-6 to m=-1.
32
The content of the Qm values of the signal SCONV (m f ch ) , which are calculated according to Figure 3.4 are
listed in the following table:
Table 3.1: Frequency content of
SCONV (m f ch ) at m=-6 to m=+6
Qm elements
Qm construction
Qm amplitude level
Q-6
Q-5
Q-4
Q-3
Q-2
Q-1
Q0
Q+1
Q+2
Q+3
Q+4
Q+5
Q+6
( A-3 )2
2 ( A-3 A-2 )
( A-2 )2
0
0
2 ( A-3 A+2 )
2 ( A-3 A+3 ) + 2 ( A-2 A+2 )
2 ( A-2 A+3 )
0
0
( A+2 )2
2( A+2 A+3 )
( A+3 )2
½ ( ATV )2
( ATV )2
½ ( ATV )2
0
0
( ATV )2
2 ( ATV )2
( ATV )2
0
0
½ ( ATV )2
( ATV )2
½ ( ATV )2
The second column of Table 3.1 specifies how the delta functions of the two tones will contribute to the
Qm elements of the signal SCONV (m f ch ) . Substituting the amplitude of the two tones, defined in Equation
3.6, the value of the distortion component is evaluated (see third column of Table 3.1). In general, the
vector Sin has a dimension 1x7 (i.e. 1xN) and the dimension of the resulting matrix is 7x13 (i.e. N x (2
N-1).
The second convolution operation in Equation 3.4 can be explained and plotted in the same way as we
did for the first convolution. The mathematical expression of this convolution is written as follows:
SCONV (f m ) ⊗ Sin (f m ) =
=
∑ (S
CONV
ρ m = − ∞ .. + ∞
(ρm ) Sin (f m − ρ m ))
∑ (S
CONV (ρ m ) Sin (m f ch − ρ m ))
Equation 3.9
ρ m = − ∞ .. + ∞
This convolution process is implemented by multiplying the result of the previous convolution
SCONV(mfCH) with the flipped and shifted version of S in . Figure 3.5 visualizes this convolution operation,
which requires in total 19 (i.e. 3N-2) shifting operations to produce the total spectrum view.
33
Figure 3.5: Visualizing the convolution operation of Equation 3.9
In this figure, one can recognize the conventional IM3 products [15] at m=±1 and m=±4 around the
spectral location of the two tones.
34
In the same way as in Figure 3.4, the convolution operation of Figure 3.5 can be mapped to a simple
matrix multiplication as shown in in Figure 3.6, where the Dm elements of the [S CONV (m f ch ) ⊗ S in (m f ch )] is
derived. Firstly, the vector SCONV is formed, which corresponds to SCONV (m f ch ) . Then the flipped
version of Sin(fm), shifted by m fCH , namely Sin(m fCH - ρm) is formed as a matrix as follows:
•
Its middle column (see Figure 3.6), the 10th column, is the flipped version of S in shifted by
zero and extrapolated with six zero’s, which is necessary to match the dimension of this
column to the dimension of SCONV (see Figure 3.5). Then multiplying SCONV by this
column provides the signal [S CONV (m f ch ) ⊗ S in (m f ch )] quantity for m=0.
•
The columns at the right of the 10th column (i.e. 11th - 19th columns) represent the flipped
version of S in , right shifted by m=+1 to +9, respectively (i.e. Sin(m fCH - ρm)). Multiplying
SCONV with those columns provides the signal [S CONV (m f ch ) ⊗ S in (m f ch )] quantity at m=+1 to
m=+9.
•
The columns at the left of the 10th column (i.e. 1st - 9th columns) represent the flipped version
of S in , left shifted by m=-9 to -1, respectively (i.e. Sin(m fCH - ρm)). Multiplying SCONV with
those columns provides the signal [S CONV (m f ch ) ⊗ S in (m f ch )] quantity at m=-9 to m=-1.
35
Figure 3.6: Matrix implementation of the convolution operation of Figure 3.5
36
The content of the Dm element values of [S CONV (m f ch ) ⊗ S in (m f ch )], which are calculated according to
Figure 3.6 are listed in the following table:
[
]
Table 3.2: Frequency content of S CONV (m f ch ) ⊗ S in (m f ch ) at m=-9 to m=+9
Dm elements
D-9
Dm construction
( A-3 )3
Dm amplitude level
( ATV )3 / (2 2 )
D-8
3 A-2 ( A-3 )2
3 ( ATV )3 / (2
D-7
3 A-3 ( A-2 )2
3 ( ATV )3 / (2
D-6
( A-2 )3
D-5
D-4
0
3 A+2 ( A-3 )2
)
2)
3
( ATV ) / (2 2 )
0
3 ( ATV )3 / (2
2
D-2
6 A+3 A-3 A-2 + 3 A+2 ( A-2 )
9 ( ATV )3 / (2
D-1
3 A+3 ( A-2 )2
3 ( ATV )3 / (2
)
2)
2)
2)
D0
D+1
0
3 A-3 ( A+2 )2
0
3 ( ATV )3 / (2
2
D-3
2
9 ( ATV ) / (2
2
6 A+2 A-2 A-3 + 3 A+3 ( A-3 )
3
6 A-3 A+3 A+2 + 3 A-2 ( A+2 )
2
9 ( ATV ) / (2
D+3
6 A-2 A+2 A+3 + 3 A-3 ( A+3 )
2
9 ( ATV )3 / (2
D+4
3 A-2 ( A+3 )2
3 ( ATV )3 / (2
D+5
D+6
0
( A+2 )3
0
( ATV )3 / (2
D+2
3 A+3 ( A+2 )
D+8
3 A+2 ( A+3 )2
D+9
( A+3 )3
)
2)
2)
2)
)
3 ( ATV ) / (2 2 )
3 ( ATV )3 / (2 2 )
( ATV )3 / (2 2 )
2
D+7
3
2
2
3
An interesting property of the DPrA matrix calculation process is its capability to track how a spectral
tones contributes to various distortion products as shown in column 2 of Table 3.2. Substituting the
amplitude of the two tones, defined in Equation 3.6, the value of the distortion component is evaluated
(see third column of Table 3.2). In summary, Table 3.2 givs insight in both the absolute and relative
contribution to the total distortion level. The content of the Dm elements in Table 3.2 is in agreement
with the conventional 3rd order distortion products derived in literature [15]) with using two tones in
time domain (i.e. two cos(ωt) signals).
37
Figure 3.7 visualizes the total 2-step matrix multiplication view of the DPrA process, which is built in the
Maple program of Appendix C. Observing the matrix dimensions of [A: 1 x N], [B: N x (2*N-1)] and [C:
(2*N-1) x (3*N-2)]. One can argue about the dimensions of the shown matrixes, which grow with the
number of the channels (i.e. N). However the intrinsic strength of the DPrA is in its matrix multiplication
form. Firstly, when the inner multiplications are between two elements where one or both of those
elements are “zeroes1”, then there is simply no need to apply a multiplication. Therefore the total
required number of multiplications may considerably be decreased. Besides that, the DPrA has another
interesting property in case the distortion in only one specific channel is of interest, e.g. when a
particular channel is proposed for use by a cognitive radio. In that case, the DPrA just needs to multiply
the Qm elements of SCONV with that associated column of the [C] matrix, defined in Figure 3.7, to that
channel number . For example, to calculate the distortion product at channel m=+4, then the DPrA just
needs to multiply the Qm elements of SCONV with the 14th column of [C] matrix.
1
Practically, zeros do not exist, as noise floor will always be there. A threshold may be specified to define this zero
content. This subject is out of our scope.
38
rd
Figure 3.7: The DPrA process to derive/predict the 3 order distortion products
Now that the distortion products across the white channels are calculated via DPrA, the equations in
chapter 2 can be directly adjusted to calculate the required IIP3 for any white Spot. The analysis of
chapter 2 derives the signal at the output of a nonlinear system when the CR signal operates at any
39
white spot among the occupied channels by the DTV primary service signals. Equation 3.10 reveals a
general form for such construction. The form consists of the amplified CR signal and the 3rd order
distortion products: Self-interference, XM3 and potential IM3 products.
Sout WhiteSpot


 jω

= Re Amplified CR signal + Self - interference + XM3 + IM 3 WhiteSpot  e
1444444424444444
3

3 order distortion products

(
)
rd
WhiteSpot t




Equation 3.10
As already concluded in Chapter 2, whether the CR operates in the IM3 spots or the IM3-free spots, the
self-interference and XM3 products will always exist. The self-interference product can be usually
neglected in comparison to the XM3 and the IM3 products because the DTV signals are usually much
stronger1 than the CR signal. The IM3 products of the DTV signals are calculated by multiplying the 3rd
order nonlinear term k3 with the DPrA results (e.g. see Figure 3.7) as follows:
Sout
WhiteSpot


 jω

= Re Amplified CR signal + XM 3 + k 3 DPrA WhiteSpot  e
14444244443 

3 order distortion products

(
(
rd
))
WhiteSpot t




Equation 3.11
Repeating the same procedure in chapter 2 by firstly writing the Signal to Distortion ratio of Equation
3.11 as follows:
S
(Amplified CR signal )
=
D
XM 3 + k 3 DPrA WhiteSpot
2
(
(
))
Equation 3.12
2
Then the amplified CR signal and the XM3 products of Equation 2.5 are substituted in Equation 3.12:


k 2A 
CR
 114243

 Amplified CR signal 
2
S
=
D 

 3k 3 2 A CR  ∑ A TV 2  + k 3 DPrA
WhiteSpot
 1444
 4n 444
42
3

XM

n
3
(
Equation 3.13






2
)
where n refers to all the channel numbers occupied by the DTV signal. ATV and ACR are the RMS values of
the DTV signal and CR signal, respectively. Then the linearity requirement can be presented as the ratio
between k1 and k3:
1
In case that the CR receiver is close to the CR transmitter, then the received CR signal is now larger than the TV
signals. An Automatic Gain Control (AGC) can be used at the CR receiver to attenuate the total received signals.
Now the attenuated TV signals will not affect the CR receiver nonlinearity anymore.
40
k1
k3
S
D
=
WhiteSpot

DPrA WhiteSpot
 3 ∑ A TV 2 +
 n
2 A CR

n




Equation 3.14
where S/D is the required Signal to Distortion Ratio for the CR receiver output. Then according to the
definition of IIP3, referred to 50Ω, the power value of IIP3 is related to k1/k3 as follows [15]:
PIIP =
3
2 k1
/ 50 [W]
3 k3
Equation 3.15
Substituting Equation 3.14 in Equation 3.15 gives the linearity requirement for any white spot:




PDPrA
S 
1
WhiteSpot
 [W]
PIIP3
= 2×
+
∑ PTV n
WhiteSpot

D 
150
P
n 23
CR44
1
14 4
424
3 

 CR - power independen t term CR - power dependent term 
Equation 3.16
2
Where:
PDPrA
WhiteSpot
 DPrA WhiteSpot 
 /50 [W]
= 


2


Equation 3.17
According to Equation 3.16, the linearity requirement depends on the following two terms:
•
A CR-power independent term: it is proportional to the sum of the power of the DTV signals,
which will cross-modulate the receiver CR-signal. This term exist in all the white spots.
•
A CR-power dependent term: it is proportional to the ratio between the IM3 distortion products
predicted by DPrA and the power value of the CR signal. This term can be minimized by
increasing the CR power signal that operates at the IM3-spots, as long as the third order
nonlinear model is still valid and no contribution yet from the 5th order nonlinear term.
The derivation of the linearity requirement can be further simplified by including the CR signal in the
vector S in , in case that its amplitude and potential spectral location is specified in advance. Then the
whole 3rd distortion products of Equation 3.10 can be directly derived by applying the DPrA to the whole
spectrum inclusive the CR signal. Hence the linearity requirements will be given as follows:
PIIP
3
Specific- WhiteSpot
=
2
S
×
150
D
PDPrA
Equation 3.18
Specific -WhiteSpot
[W]
PCR
Where PCR is now added to the construction of PDPrA
.
Specific - WhiteSpot
41
Application example As an example, let’s apply our analysis in predicting the linearity requirement to
channel 4 and 8 in Figure 3.7. From Table 3.2, the DPrA has predicted the following distortion levels:
Table 3.3: DPrA frequency content at m=±4 and m=±6
Dm elements
D-4 + D+4
Dm construction
3 A+2 ( A-3 )2 + 3 A-2 ( A+3 )2
3
D-6 + D+6
( A-2 ) + ( A+2 )
3
Dm amplitude level
( )
/ (2 2 )
2 x 3 ( ATV )3 / 2 2
2 x ( ATV )3
Then the linearity requirements at channels 4 and 6 are given as follows:


P


DPrA
S
1
Channel 4
 ∑ PTV +
 [W]
PIIP3
= 2×
n
Channel 4
D n = 2 and 3
150
PCR

4
4244
3 144
42444
3
1Fixed
term
Variable term


Equation 3.19
2
where : PDPrA
PIIP
3
Channel 4
  3 (A )3  
TV  
2
2
  2 2 
 3 (A )3 


TV  /50 [W]
 /50 = 
=


2
2













S 
1 PDPrA
= 2×
 ∑ PTV +
 [W]
D n1
150
P
= 2 and 3
CR
424
3 1442443 
 Fixed term

Variable term
Equation 3.20
Channel 6
Channel 6
n
  (A TV )3  
 2

2
 (A TV )3 
  2 2  
 /50 [W]
=
 /50 = 

2
2








2
where : PDPrA
Channel 6
Finally the analysis can be expanded to deal with wideband DTV signals instead of tones in the same
manner as in chapter 2 by introducing ∆Wideband (see Equation 2.6, maximum ∆Wideband=4.5dB) as a
correction factor to Equation 3.16 as follows:
PIIP
3
WhiteSpot
= 2×
S
D
∆
∆

PDPrA
 10 10 × 0.66 ∑ PTV + 1 10 10 × (0.66 or 0.17 )
1
4
4
2
4
4
3
244443
150 14444
PCR

n
1.4
1.4 or 0.7

Wideband
Wideband
WhiteSpot
n

 [W]


Equation 3.21
The choice between 0.66 or 0.17 depends on whether the white spot lies in the center spot or on one of
the two side spots of the Bell shape (see Figure 2.11).
42
3.2 DPrA Application Examples
To verify the DPrA concept, this section begins with applying DPrA on the chapter 2’s examples of Figure
2.5 and Figure 2.6 (i.e. with 3 interferers), which are verified by measurements shown in Figure 2.13 and
Figure 2.14, respectively. Further as the number of interferers increases, then the construction pattern
of 3rd order distortion products becomes more complex. This complexity is explored at the last example.
First spectrum scenario: There are three interferers at the channels/spots at 1, 6 and 19 (see Figure 2.5).
The 3rd order distortion products view is presented in Table 3.4 and visualized in Figure 3.8. The third
column (i.e. Signal amplitude level) is in exact consistent with Equation 3.3.
Figure 3.8: DPrA prediction of the distortion products for Interferers at Channels 1, 6 and 19
rd
Table 3.4: 3 order distortion products view
Ch.
No.
1
Signal amplitude level
Signal construction
+1 :=6 k A A A + 6 k A A A + 3 k A 2 A
3 -19 1 19
3 -6 1 6
3 1 -1
3
2
2
3
:=3 k3 2 ATV19 ATV1 + 3 k3 2 ATV6 ATV1 + k3 ATV1 2
2
-1 :=3 k A 2 A + 6 k A A A + 6 k A A A
3 -1 1
3 -19 -1 19
3 -6 -1 6
6
+6 :=6 k A A A + 3 k A 2 A + 6 k A A A
3 -19 6 19
3 6 -6
3 -1 6 1
3
2
3
2
:=3 k3 2 ATV19 ATV6 + k3 ATV6 2 + 3 k3 2 ATV1 ATV6
2
-6 := 3 k A 2 A + 6 k A A A + 6 k A A A
3 -6 6
3 -19 -6 19
3 -6 -1 1
11
14
19
2
3
2
:= k3 2 ATV6 ATV1
2
+11
:=3 k3 A-1 A6
-11
:=3 k3 A-6 A1
+14
:=6 k3 A-6 A1 A19
-14
:=6 k3 A-19 A-1 A6
+19
:=3 k3 A19 A-19 + 6 k3 A-6 A19 A6 + 6 k3 A-1 A19 A1
2
:=3 k3 2 ATV6 ATV1 ATV19
2
-19 :=3 k A 2 A + 6 k A A A + 6 k A A A
3 -19 19
3 -19 -6 6
3 -19 -1 1
43
3
3
2
2
:= k3 ATV19 2 + 3 k3 2 ATV6 ATV19 + 3 k3 2 ATV1 ATV19
2
Second spectrum scenario: There are three interferers at the channels/spots at 3, 5 and 16 (see Figure
2.6). The 3rd order distortion products view is presented in Table 3.5 and visualized in Figure 3.9. The
third column is in exact consistent with Equation 3.3.
Figure 3.9: DPrA prediction of the distortion products for Interferers at Channels 3, 5 and 16
rd
Table 3.5: 3 order distortion products view
Ch.
No.
1
3
Signal amplitude level
Signal construction
3
2
:= k3 2 ATV5 ATV3
2
2
+1
:=3 k3 A-5 A3
-1
:=3 k3 A-3 A5
2
+3 :=6 k A A A + 6 k A A A + 3 k A 2 A
3 -16 3 16
3 -5 3 5
3 3 -3
3
2
2
3
:=3 k3 2 ATV16 ATV3 + 3 k3 2 ATV5 ATV3 + k3 ATV3 2
2
-3 :=6 k A A A + 6 k A A A + 3 k A 2 A
3 -16 -3 16
3 -5 -3 5
3 -3 3
5
+5 :=6 k A A A + 3 k A 2 A + 6 k A A A
3 -16 5 16
3 5 -5
3 -3 5 3
3
2
3
2
:=3 k3 2 ATV16 ATV5 + k3 ATV5 2 + 3 k3 2 ATV5 ATV3
2
-5 :=6 k A A A + 6 k A A A + 3 k A 2 A
3 -16 -5 16
3 -5 -3 3
3 -5 5
7
14
16
:=3 k3 A-3 A5
-7
:=3 k3 A-5 A3
+14
:=6 k3 A-5 A3 A16
-14
:=6 k3 A-16 A-3 A5
2
:=3 k3 2 ATV5 ATV3 ATV16
+16 :=3 k A 2 A + 6 k A A A + 6 k A A A :=3 k ATV 3 2 + 3 k 2 ATV2 ATV + 3 k 2 ATV2 ATV
16
3
5
16
3
3
16
2 3
3 16 -16
3 -5 16 5
3 -3 16 3
-16 :=6 k A A A + 6 k A A A + 3 k A 2 A
3
18
3
2
:= k3 2 ATV5 ATV3
2
2
+7
-16 -5 5
3 -16 -3 3
+18
:=6 k3 A-3 A5 A16
-18
:=6 k3 A-16 A-5 A3
3
-16
16
:=3 k3 2 ATV5 ATV3 ATV16
44
Referring to the theory of chapter 2 and observing Table 3.5, the construction pattern of the 3rd order
distortion products across the channels are as follows:
1.
Channel 3 contains the cross-modulation (i.e XM3) products of the interferers of channel 5 and
channel 16 on top of the interferer of channel 3 plus its self-interference product:
3
2
2
3
:=3 k3 2 ATV16 ATV3 + 3 k3 2 ATV5 ATV3 + k3 ATV3 2 .
2
2. Channel 5 contains the cross-modulation (i.e XM3) products of the interferers of channel 3 and
channel 16 on top of the interferer of channel 5 plus its self-interference product:
3
2
3
2
:=3 k3 2 ATV16 ATV5 + k3 ATV5 2 + 3 k3 2 ATV5 ATV3 .
2
3. Channel 16 contains the cross-modulation (i.e XM3) products of the interferers of channel 3 and
channel 5 on top of the interferer of channel 16 plus its self-interference product:
3
3
2
2
:= k3 ATV16 2 + 3 k3 2 ATV5 ATV16 + 3 k3 2 ATV3 ATV16 .
2
4. Both channel 1 and channel 7 contain the conventional intermodulation products between
3
2
3
2
2
2
channel 3 and channel 5::= k3 2 ATV5 ATV3 and:= k3 2 ATV5 ATV3 , respectively.
5. Channel 14 contains the intermodulation product among channel 3, channel 5 and channel 16.
The distortion level in channel 14 is 2 times higher (i.e 6dB) than the conventional
intermodulation products of point 4:
:=3 k3 2 ATV5 ATV3 ATV16 , as presented in chapter 2.
6. Channel 18 contains the intermodulation product among channel 3, channel 5 and channel 16.
The distortion level in channel 18 is 2 times higher (i.e 6dB) than the conventional
intermodulation products of point 4:
:=3 k3 2 ATV5 ATV3 ATV16 , as presented in chapter 2.
45
Third spectrum scenario: There are four interferers at the channels/spots at 2, 6, 13 and 20. The 3rd
order distortion products view is quantified in Table 3.6 and visualized in Figure 3.10.
Figure 3.10: DPrA prediction of the distortion products for Interferers at Channels 2, 6, 13 and 20
rd
Table 3.6: 3 order distortion products view
Ch.
Signal construction
No.
2 + :=6 k3 A-20 A2 A20 + 6 k3 A-13 A2 A13 + 6 k3 A-6 A2 A6 + 3 k3 A22 A-2
Signal amplitude level
3
2
2
2
3
:=3 k3 2 ATV20 ATV2 + 3 k3 2 ATV13 ATV2 + 3 k3 2 ATV6 ATV2 + k3 ATV2 2
2
- :=6 k3 A-20 A-2 A20 + 6 k3 A-13 A-2 A13 + 6 k3 A-6 A-2 A6 + 3 k3 A-22 A2
6
9
10
13
16
17
2
:=6 k3 A-20 A6 A20 + 6 k3 A-2 A6 A2 + 6 k3 A6 A13 A-13 + 3 k3 A-20 A13 + 3 k3 A6 A-6
-
:= 6 k3 A-13 A-6 A13 + 6 k3 A-6 A-2 A2 + 6 k3 A-6 A20 A-20 + 3 k3 A20 A-13 + 3k3 A-6 A6
+
:= 6 k3 A-13 A2 A20 + 6 k3 A-6 A2 A13
-
:=6 k3 A-13 A-2 A6 + 6 k3 A-20 A-2 A13
+
:=3 k3 A-2 A6
-
:=3 k3 A-6 A2
+
-
:=6 k3 A-2 A13 A2 + 6 k3 A6 A13 A-6 + 6k3 A13 A20 A-20 + 6 k3 A20 A-13 A6 + 3 k3 A-13 A13
:= 6 k3 A-20 A-13 A20 + 6 k3 A-13 A-6 A6 + 6 k3 A-13 A-2 A2 + 6 k3 A13 A-20 A-6 + 3 k3 A13 A-13
+
:= 6 k3 A-6 A2 A20
-
:= 6 k3 A-20 A-2 A6
+
:= 6 k3 A-2 A6 A13
2
+
-
2
2
2
2 3
2
3 k3 2 ATV20 ATV6 + 3 k3 2 ATV2 ATV6 + 3 k3 2 ATV6 ATV13 + k3 2 ATV20 ATV13
2
3
3
+ k3 ATV6 2
2
:=3 k3 2 ATV13 ATV2 ATV20 + 3 k3 2 ATV6 ATV2 ATV13
2
3
2
:= k3 2 ATV6 ATV2
2
2
2
2
2
2
2
3k3 2 ATV2 ATV13 +3k3 2 ATV6 ATV13 +3k3 2 ATV13 ATV20 +3k3 2 ATV20 ATV13 ATV6
3
3
+ k3 2 ATV13
2
:=3 k3 2 ATV6 ATV2 ATV20
:=3 k3 2 ATV6 ATV2 ATV13
:= 6 k3 A-13 A-6 A2
20
2
+
2
2
:=3 k3 A20 A-20 + 6 k3 A-6 A20 A6 + 6 k3 A-2 A20 A2 + 6 k3 A13 A20 A-13 + 3 k3 A-6 A13
2
2
:=6 k3 A-20 A-13 A13 + 6 k3 A-20 A-2 A2 + 6 k3 A-20 A6 A-6 + 3 k3 A-20 A20 + 3 k3 A6 A-13
46
3
3
2
2
2
k ATV20 2 + 3 k3 2 ATV6 ATV20 + 3 k3 2 ATV2 ATV20 + 3 k3 2 ATV20 ATV13
2 3
3
2
+ k3 2 ATV6 ATV13
2
Observing Table 3.6, the construction pattern of the 3rd order distortion products across the channels
are as follows:
1.
Channel 2 contains the cross-modulation (i.e XM3) products of the interferers of channel 6,
channel 13 and channel 20 on top of the interferer of channel 2 plus its self-interference
2
2
2
3
2
3
:=3 k3 2 ATV20 ATV2 + 3 k3 2 ATV13 ATV2 + 3 k3 2 ATV6 ATV2 + k3 ATV2
product:
2
2. Channel 6 contains the cross-modulation (i.e XM3) products of the interferers of channel 2,
channel 13 and channel 20 on top of the interferer of channel 6 plus its self-interference
product plus the conventional intermodulation products between channel 13 and channel 20:
3. Channel 13 contains the cross-modulation (i.e XM3) products of the interferers of channel 2,
channel 6 and channel 20 on top of the interferer of channel 13 plus its self-interference
product plus the intermodulation products among channel 6, channel 13 and channel 20:
4. Channel 20 contains the cross-modulation (i.e XM3) products of the interferers of channel 2,
channel 6 and channel 13 on top of the interferer of channel 20 plus its self-interference
product plus the conventional intermodulation products between channel 6 and channel 13:
5. Channel 9 contains the intermodulation products among channel 2, channel 13 and channel 20
plus the intermodulation products among channel 2, channel 6 and channel 13:
:=3 k3 2 ATV13 ATV2 ATV20 + 3 k3 2 ATV6 ATV2 ATV13
6. Channel 10 contains the conventional intermodulation products between channel 2 and channel
3
2
2
:= k3 2 ATV6 ATV2
6:
7. Channel 16 contains the intermodulation products among channel 2, channel 6 and channel 20:
:=3 k3 2 ATV6 ATV2 ATV20
8. Channel 17 contains the intermodulation products among channel 2, channel 6 and channel 16:
:=3 k3 2 ATV6 ATV2 ATV13
47
3.3 Conclusions
The amplitude and the spectral location of the DTV signals are place and time dependent, and hence the
required linearity for a wideband CR receiver varies across the unoccupied DTV channels and varies over
place and time. The DPrA aims at predicting the instantaneous distortion levels in a dynamic way by
matrix multiplications on spectrum analysis measurements for each frequency channel. It is possible to
implement this DPrA in the digital domain of the receiver. The dimension of the required DPrA matrices
grows with the number of the covered channels, whereas the number of the inner multiplications grows
with the number of the occupied DTV channels.
The matrix algorithm used in the DPrA can be further extended to calculate the distortion products
generated by higher nonlinear order terms (e.g. 4th, 5th ... etc.). For example for the 4th order nonlinear
term, an extra matrix multiplication (i.e. due to an extra convolution) to the outcome of Figure 3.7 is
required. The outcome of this multiplication must be multiplied with k4 that characterize the 4th third
order nonlinearity of the circuit. Note that the 2nd distortion products is already calculated in Figure 3.4
by just one multiplication (i.e. due to one convolution operation) between vector S in and the
convolution matrix [A]. All the results of this chapter are verified by MatLab simulation.
48
Chapter 4. Analysis of a Very Linear Front-End [41]
Chapter 4 and 5 investigate the third research question: “How to increase the linearity of CR receiver
hardware, given less up-front filtering?”
This chapter conceptually presents the linearization technique via the application of negative
conductance on the receiver topology that is introduced in section 4.1. Then, section 4.2 gives an
intuitive model to understand the basic distortion cancellation concept. The optimum negative
conductance value is derived by mathematical analysis in section 4.3 and related to negative feedback
theory in section 4.4. Section 4.5 analyses stability issues related to this negative conductance
technique. Finally, some conclusions are given in section 4.6.
This chapter is a literal copy of the first part of the paper published in [41].
49
4.1 Motivation and Introduction
As discussed in the previous chapters, wideband receivers with high IIP3 are very important for
opportunistic DSA via a cognitive
radio.
Strong
interferers
(e.g.
incumbent TV signals) may present in
near or directly adjacent channels
inside the down converted band of
the receiver, making on-chip RF
filtering ineffective, as exemplified in
Figure 4.1 for a Digital TV band.
Figure 4.1: CR operates in an adjacent channel of Digital TV spectrum [25]
High linearity is required also to prevent cross-modulation effects [18] from desensitizing the receiver. A
part from the RF receivers, the spectrum sensing front-end also requires high in-band IIP3 in order to
minimize the errors in detecting the empty channels in the spectrum [26].
In general, linearity requirements on radio receivers become increasingly challenging, as the radio
spectrum becomes more and more crowded. Moreover, there is a trend towards more wideband and
more flexible radio hardware with less dedicated RF filtering (“Software Defined Radio”). As an example,
Figure 4.2 plots IIP3 requirements calculated
for E-UTRA for a wideband base station
receiver in three scenarios: wide area, local
area and home [27]. Apart from the high
100MHz bandwidth, note the sudden step in
IIP3 requirements at the band-edge. Also note
that less coverage area (home versus wide
area), corresponds to higher in-band IIP3 but a
smaller step to out-of-band IIP3 (i.e. around
Figure 4.2: Example IIP3 requirement for E-UTRA [29]
16dB for home area versus 40dB for wide
area). As a consequence of the lack of a reasonable transition band, on-chip analog filtering is ineffective
to relax the IIP3 requirement, and off-chip filters are expensive. Depending on the blocker scenario,
compression point requirements may or may not be affected. This chapter proposes a circuit technique
that can increase IIP3 simultaneously for in- and out-of-band, at roughly constant compression point.
50
Strong RF interference can easily clip baseband amplifiers, while higher required bandwidths limit the
amount of available loop-gain for negative feedback. When pushing linearity, avoiding voltage gain at RF
(see Figure 4.3) is instrumental [[28]-[34]]. Exploiting RF V-I conversion followed by passive down-mixing
and then simultaneous I-V conversion and filtering at IF/baseband with OpAmps, an out-of-band IIP3 of
up to +18dBm has been shown [[28],[29]]. Passive mixer-first architectures can even achieve up to
+25dBm out-of-band IIP3 [32]. However in-band IIP3 is much worse, certainly at high gain. The best inband IIP3 results that we found for receivers were +3.5dBm for [28] at 34dB gain and +11dBm for [31] at
19dB gain. Analysis shows that finite OpAmp gain can be a bottleneck, as a non-zero virtual ground node
voltage can result in distortion currents. In [35], we proposed to exploit a negative conductance
technique to cancel distortion currents. In this way, the design of the OpAmp is relaxed and its
performance no longer needs to be a bottleneck. The use of a negative conductance has been proposed
in [36] to realize TIA flicker noise shaping. Paper [36] also briefly mentions linearity improvement, but
linearity benefits were not the focus there.
Figure 4.3: High blocker tolerant linear receiver
51
4.2 Linearization Concept Analysis
To understand the OpAmp linearity limitation and the distortion cancellation technique intuitively, it is
instructive to follow a 4-step approach to analyze what happens at the virtual ground node “VGND”, as
illustrated in Figure 4.4 to Figure 4.8:
Step 1: Assume the RF V-I conversion and mixing are perfectly ideal (i.e. linear and infinite current
source resistance for GM), we can use the equivalent baseband model in Figure 4.4 (omitting the
downconversion for simplicity). Assuming a 2-tone input signal VS(f), the injected current IS(f) to the
VGND node is linear, so without IM3 tones. Now, if the OpAmp handles large signals at a high but finite
gain, its output stage will produce IM3 products at the OUT node, i.e. VOUT(f). However, as IS(f) has no IM3
and the feedback resistor RF is linear, the voltage over RF does not contain IM3 (assuming negligible
OPAMP input current). Consequently, the IM3 products of VVGND(f) are in absolute sense equal to those
of VOUT(f) both in magnitude and phase. Let’s denote this “IM3 copy” effect in Figure 4.4 as “problem A”.
Note that the two main tones of VVGND(f) are much smaller than that of VOUT(f), as the ratio VOUT(f)/VVGND
(f) for linear terms is equal to the loop gain. As a consequence the ratio between the linear terms and
the IM3 products at VGND node is much worse than at the OUT node, causing a more serious problem
discussed next.
Figure 4.4: OpAmp nonlinearity problem A: IM3 is copied from the OUT node to the VGND node
52
Step 2: Assume we add a finite output resistance RO as shown in Figure 4.5. The nonlinear voltage
VVGND(f) over RO now generates a nonlinear current IO(f), and hence IF (f) becomes nonlinear. This current
is absorbed by the OpAmp output stage and increases IM3 at both VOUT(f) and VVGND(f). We will denote
this “RO loading” effect on the VGND node in Figure 4.5 as “problem B”.
Figure 4.5: OpAmp nonlinearity problem B: RO loads the VGND node
Step 3: Once one realizes the main cause for distortion current is VVGND(f)/RO, it is easy to verify that
adding a negative conductance with value GO=1/RO between VGND and ground can be a solution (see
Figure 4.6). The negative conductance senses VVGND and generates a copy of the distorted current IO(f),
which now flows in a “local circle” via the ground. Consequently, the current injected to the VGND node
becomes linear again and we are back at the circuit of problem A, having solved problem B.
Figure 4.6: Solving problem B via negative conductance with GO=1/RO
53
Step 4: Still, the OpAmp output voltage contains some IM3, equal to that on the VGND node. By slight
overcompensation this IM3 contribution can also be cancelled. To show this, it is useful to model the
floating resistor RF with an equivalent network consisting of four single-ended linear transconductor
blocks GF (GF=1/RF), all referred to ground as shown in Figure 4.7(a). The two shorted GF blocks,
indicated with a dashed ellipse, can be replaced by a simple RF resistor to the ground (see Figure 4.7(b)).
Thus Figure 4.7(c) results with RF-VGND and RF-OUT, (loading resistances at the VGND node and the OUT
node, respectively), GF-VGND (the transconductance sensing VOUT and injecting current to the VGND node),
and GF-OUT (the transconductance sensing VVGND and injecting current to the OUT node). We assigned
different names to GF and RF blocks in order to distinguish between their effects on nonlinearity at the
VGND node and the OUT node separately. Figure 4.7(c) clearly shows the loading effect of RF (i.e. RF-VGND)
at the VGND node.
Figure 4.7: Equivalent model of the effect that RF has on the OUT node and the VGND node
54
Now, when the negative conductance cancels this loaded effect (see Figure 4.8), the injecting current of
GF-VGND becomes equal to the linear current IS. As VOUT=IS/GF-VGND=-IS.RF, the OpAmp output voltage VOUT
becomes linear. In this way, problem A is solved as well.
Figure 4.8: Solving problem A via negative conductance with GF = 1/RF
Overall, combining the solutions for problem A and B, the optimal total negative conductance is:
GTOTAL=1/RO+1/RF.
55
4.3 Mathematical Analysis
To mathematically prove this optimum cancellation condition, the OpAmp is modeled as an OTA with
nonlinear transconductance and also a nonlinear output resistance because we aim for high output
swing:
I F = gm1VIN + gm 3 VIN + go1VO + go 3 VO
3
3
Equation 4.1
In the model, we assume that the third order nonlinearities are more pronounced than the second order
nonlinear terms, which is reasonable considering the OpAmp will be implemented in fully differential
form.
In Appendix D, the nonlinear relation between VOUT and signal current IS is derived using the model in
Figure 4.9.
Figure 4.9: Baseband model with RO and the extended RF for nonlinearity derivations
It can be expressed in terms of a linear (Ω1) and third-order nonlinear (Ω3) coefficient:
VOUT = Ω 1 I S + Ω 3 I S
3
Equation 4.2
The linear coefficient Ω1 is the I/V conversion gain:
Ω1 =
1
1  1
1
+
 
 a  R O R F − VGND
Equation 4.3


 + G F − VGND 


56
Where (a) is a function of the linear terms of the OpAmp model (i.e. gm1, go1) and the RF effects at the
OUT node (i.e. RF-OUT and GF-OUT). For very high gm1, (a) reaches -∞. Consequently, the I/V conversion
gain of Equation 4.3 becomes 1/GF-VGND= -RF.
The third-order distortion coefficient (Ω3) is:
Ω3 =
 1

1

NL3 
+
 R O R F− VGND 
1  1


1
 + G F − VGND 
+
 
 a  R O R F− VGND 

Equation 4.4
4
where (NL3: see Appendix D) is related to the nonlinear terms of the OpAmp model and is a function of
(i.e. gm1, gm3, go1 and go3) and the effect of RF on the OUT node (i.e. RF-OUT and GF-OUT). Now, if the
negative conductance technique cancels 1/RO+1/RF-VGND from Equation 4.3 and Equation 4.4, we see that
Ω1 reaches 1/GF-VGND=-RF and Ω3 becomes zero (i.e. distortion is cancelled). Note that since the voltage
swing at the VGND node is small, the effect of negative conductance nonlinearity can be very small.
To verify the OpAmp model, we fitted the model derived above to simulations done for the OpAmp that
will be introduced in the next chapter. Figure 4.10 shows a close agreement.
Figure 4.10: OpAmp model Equation 4.1 verification
57
4.4 Feedback Control Analysis
The linearity benefit can also be verified by applying feedback theory to Figure 4.9 as shown in Figure
4.11, excluding GM. The feedback topology of the circuit is Voltage-Current Feedback [33]. The output
voltage (i.e. VOUT) is sensed and converted to a proportional feedback current βVOUT, where β=GF-VGND (in
Siemens). This feedback current is subtracted from the input current IS resulting in an error current
Ierror=IS-βVOUT to be amplified by the block A. Here, A=VOUT/Ierror, where A has the dimension of a
transimpedance [Ω]. It consists of all the blocks of Figure 4.9, excluding GM and GF-VGND. Actually for
finite A, there will be a non-zero Ierror due to the loading effect of RO and RF-VGND on the VGND node. Now
the negative conductance increases the input impedance of the A block to infinity by cancelling RO and
RF-VGND, so that Ierror becomes zero and A=VOUT/Ierror=infinity. Consequently, loopgain Aβ goes to infinity
and VOUT/IS achieves its ideal value 1/β=RF, i.e. perfect linearity. We conclude that the negative
conductance technique increases the loop gain by increasing the value of A. Also note that only a finite
value for GO is needed to make the loopgain theoretically approach infinity, which is not possible by
increasing gm1 in the gain block. Although the feedback theory puts the application of a negative
conductance technique in the right context, however the problem with control theory is that it assumes
blocks
with
operation,
unilateral
which
are
sometimes not easy to
identify (e.g. see Figure
4.11: feedback resistor RF
which
is
supposed
to
realize the β block also
becomes part of the A
block). In compare to the
feedback
analysis
analysis,
explains
our
in
a
simple way how IM3 is
Figure 4.11: Applying feedback theory to Figure 4.9, excluding GM
affected by RO and RF.
Now, before we finalize this chapter with conclusions, we will first deal with a potential caveat of
negative conductance: the risk of instability
58
4.5 Stability Analysis
We will consider two stability aspects: 1) the risk of oscillation, based on a small signal model, and 2) the
risk of latch-up. Let us first look at the small signal behavior, referring to Figure 4.12. As the low-pass
filtering is desired, CF is added as feedback capacitor. Capacitor CT models the total input capacitance to
ground of the OpAmp CIN-OpAmp and other capacitance CO at the VGND node (see Figure 4.3).
Figure 4.12: Circuit diagram for small signal stability analysis
For simplicity, the OTA is modeled as a frequency dependent transconductance with a dominant pole at
ωO and infinite output impedance:
gm(s) =
gm O
s
1+
ωO
Equation 4.5
Assuming no further loading at the OUT node, looking into the VGND node (see Figure 4.12), the
impedance (ZIN) consists of the reactance of CT in parallel to 1/gm(s):




1
1
1
1
s

Z IN =
//
=
// 
+
s C T gm(s) s C T  gm O
gm O ω O 
1
424
3 
 {
Resistance Inductance (L)
Equation 4.6
Therefore, a parallel RLC tank is seen looking into the VGND node. If the negative conductance would
both cancel 1/RO and gmO, then oscillation would happen at a resonance frequency that depends on the
value of CT and L (i.e. fres=1/ 2
)). However, note that the typical virtual ground impedance 1/gmO
will normally be much lower than RO and RF. Thus, as the negative conductance GTotal is designed to
cancel 1/RO and 1/RF, the point of small signal instability can be designed to be safely far away.
59
Let’s now look at the potential of latch-up of the OpAmp for a case that the negative conductance is too
strong, i.e. it produces more current than needed after compensating the current in RO. As shown in
Figure 4.13, the negative conductance injects current via RF (i.e. VVGNDGLatch-up-Risk) that needs to be
handled by the OpAmp output stage in addition to the main current coming from GM (i.e. IS):
IOpAmp-Latch-up-Risk = IS + VVGNDG Latch−up−Risk
= IS +
1



1
− G Latch −up−Risk  + a G F−VGND 


 R F−VGND

ISG Latch −up−Risk
Equation 4.7
Where the relation between VVGND and IS is derived in Appendix E.
Figure 4.13: Latch-up problem at the OUT node
Referring to Figure 4.13 and substituting GLatch-up-Risk=(1/RF)+∆G in Equation 4.7 gives the following
relation:
I OpAmp-Latch -up - Risk

1
∆G +

RF
= I S 1 −
 
a 

  ∆G +
R F 
 






Equation 4.8
The OpAmp output stage current flows throw RF and make a voltage drop. The peak of this voltage drop
is around VDD/2-VOpAmpOutputStage-OV, where VOpAmpOutputStage-OV is the over drive voltages of the OpAmp
output stage transistors. Hence, if very strong negative conductance has been used (i.e. high ∆G), then
the current of Equation 4.8 becomes higher than the OpAmp output stage current capability and the
latch-up occur.
60
4.6 Conclusion
Front-End topologies with RF V-I conversion followed by passive down-mixing and then simultaneous I-V
conversion and filtering at IF/baseband with OpAmps are quite linear. However their linearity is mainly
limited by the OpAmp where a very high swing appears at its output stage transistors. Analysis shows
that finite OpAmp gain can be a bottleneck, as a non-zero virtual ground node voltage can result in
distortion currents. Our concept proposes to exploit a negative conductance technique to cancel those
distortion currents. A mathematical derivation proves is the optimum/required negative conductance
value, which exactly agrees with our intuitive analysis. Additionally, from the feedback theory context,
this optimum negative conductance simply cancels the error current so that the loop gain theoretically
approaches infinity, which is not possible by increasing the OpAmp gain to infinity. In this way, the
design of the OpAmp is relaxed and its performance no longer needs to be a bottleneck. The risk of
stability is also analyzed, where the system is well stable around the optimum negative conductance.
61
62
Chapter 5. Chip Design: UTCRxVLi [41]
Chapter 5 follows chapter 4 in investigating the third research question: “How to increase the linearity
of CR receiver hardware, given less up-front filtering?”
A receiver prototype in 65nm STMicroelectronics is implemented to verify the concept of using a
negative conductance technique to linearize the whole receiver. This chapter discusses the receiver
design details in section 5.1. The receiver noise figure analysis including the negative conductance
contribution is discussed in section 5.2. This noise analysis and the analysis of the previous chapter are
all verified by measurements in section 5.3. The results are benchmarked to other highly linear receivers
in section 5.4. Finally, the chapter is ended with conclusions in section 5.5.
This chapter is a literal copy of the second part of the paper published in [41].
63
5.1 Receiver Implementation
We will now apply the negative conductance idea to highly linear zero-IF radio receiver architecture of
Figure 4.3 in the previous chapter, which is redrawn in Figure 5.1.
Figure 5.1: High blocker tolerant linear receiver
To demonstrate the linearity potential of this technique, we will replace the active V-I conversion by a
more linear fully passive mixer with resistors in series [29], as shown in Figure 5.2.
Figure 5.2: Replacing the LNTA (GM) of by a linear resistance RRF
64
Figure 5.3: Complete Receiver with distortion compensation by –GO
Figure 5.3 shows the complete front-end IC schematic including the negative conductance. Using the
equivalent model in Figure 4.5, we can model the RF part of each branch in I and Q as a grounded
resistor RO and a transconductor GM referred to ground as denoted in Figure 5.3. However, as resistor
RRF is in series with the mixer on-resistance RON-MIXER and the virtual ground impedance RVGND of the
OpAmp, the equivalent GM now equals 1/(RRF+RON-MIXER+RVGND). This is chosen to be 20mS to realize RF
input impedance matching of 50Ω, assuming perfect non overlapping 25% duty-cycle clocks, so the RFinput continuously sees a conduction path to ground. The equivalent output impedance of the mixer at
baseband now is RO=2(RBalUn+RRF+ RON-MIXER), where the factor 2 is due to the quadrature mixer with 25%
duty cycle, connecting each I and Q baseband part to RF two times per LO cycle. To understand this
point, let’s derive RO from the power that is delivered by a test voltage source (i.e. Vtest=Va cos(ωLOt))
“looking back” in RO as shown in Figure 5.4.
65
Figure 5.4: Derivation of RO
This source is connected to the first branch of the I-path. The current Itest will flow through
RON-MIXER+RRF+RBalUn two times LO-cycle, hence we get:
3TLO
 TLO

2
4
4
1 
Va

P=
V
I
dt
+
V
I
dt
=
test test
∫ test test  4 (R ON-MIXER + R BalUn + R RF
TLO  ∫0
TLO


2
)
Equation 5.1
This power must be equal to the power dissipation in RO:
P=
1
TLO
TLO
TLO
∫
0
2
2
Vtest
V
dt = a
RO
2R O
Equation 5.2
By equating Equation 5.1 and Equation 5.2, the following RO is derived:
RO = 2 (RON-MIXER+ RBalUn + RRF )
Equation 5.3
In the RO derivation, the power is only balanced with the fundamental, while the effect of the 3rd and
higher harmonics are neglected due to the existence of CO (see Figure 5.3). Now, the 50Ω input
impedance matching is implemented as a combination of series resistances RRF≈12Ω, the up-converted
impedances of the passive mixer switches RON-MIXER≈28Ω plus the VGND impedance RVGND≈7Ω. The
passive mixer consists of simple NMOS switches. CO=8pF effectively shorts the LO leakage and high IF
frequency components to ground. The TIA consists of a class-A input stage and a class-AB output stage,
to maximize output swing (see Figure 5.5, [37]). Common mode feedback ensures biasing at VDD/2.
66
Figure 5.5: Circuit Diagram of the fully differential OpAmp design [37]
The feedback impedance is RF=1.5kΩ and CF=8pF, to obtain 26dB voltage gain and a -3dB-bandwidth of
12MHz. The differential topology allows for a simple differential implementation of the negative
conductance (right part of Figure 5.3) and high IIP2. To be able to measure what is the effect of different
negative conductance values, -GO is implemented as a parallel array of identical “unit-transconductors”,
digitally controllable via multiplier M, with transconductance steps of 0.2mS. Thus M=28 renders
GO=5.6mS to compensate the nominal value of RO=180Ω (RO=2(RBalUn+RRF+RON-MIXER)=2(50+12+28)=180Ω).
We will now consider the noise degradation resulting from the introduction of the negative
conductance. Actually this noise can be cancelled by a noise cancellation path [[38],[29]], however this is
expected to result in a linearity bottleneck in the auxiliary noise cancellation path. Hence we will analyze
the noise figure degradation and aim for minimizing the noise penalty.
67
5.2 Noise Figure
Receiver topologies with a passive mixer and transimpedance amplifier (TIA), can suffer from
amplification of OpAmp noise [39]. The output referred OpAmp noise contribution can be written as:
2
2
n−OUT
V
 R 
= 1+ F  Vn2−OpAmp
 RO 
Equation 5.4
Where Vn-OpAmp refers to the (equivalent) input noise voltages of the OpAmp, RO and RF are as used in
Figure 4.5. For our design RF=1.5kΩ and RO=180Ω, then the amplification factor is equal to
(1+RF/RO)2=87.
Often a high RF V/I conversion (GM-value) is used to achieve an overall noise figure around or below
3dB. Here we will use 20mS, the value desired for input impedance matching.
Figure 5.6 shows a baseband model of Figure 5.2 with noise sources added. The noise of GM is
represented by the current noise source (In-Ro) of RO. The noise of RF is modeled via voltage noise source
Vn-RF, while In-GTotal represents the current noise source of the negative conductance.
Figure 5.6: Equivalent baseband model for Noise Figure analysis
For simplicity, the OpAmp is modeled as a simple Transconductance (gm). To analyze the noise
contributions of In-Ro and In-GTotal to the output voltage, Ω1 (i.e. the I/V conversion of the TIA (3)) is useful.
The straightforward NF analysis shows:


2
2

 Vn2− OpAmp 1 
 1

 1
 
2 1
 Ω1 

1 + R F 

NF = 1 +
2
 R + γ G Total  + 2
 R − G Total   + R F  gm R − 1 

R O 
I
1

O
O

F






n
−
R
O
3 144
 Ω1GM  R S  14442444
444424444443 1442443 
First Term
Third Term
2

Second Term


1
68
Equation 5.5
The first term between the square brackets in Equation 5.5 shows that the negative conductance GTotal
has a direct noise contribution to the output. Its noise contribution is scaled by /
. The
“noise excess factor” γ can be minimized to around 2/3 (i.e. theoretically) by choosing a non-minimum
channel length for the negative conductance transistors. Long-channel transistors are preferred for 1/f
noise. We used 1µm channel length in this design. The second term is the mentioned amplification
factor Equation 5.4 of OpAmp noise including the negative conductance effect (GTotal). It is interesting to
observe that this term reaches zero when the negative conductance reaches GTotal. However, the direct
noise contribution of the negative conductance is much higher than the canceled OpAmp noise
contribution. Hence the total noise figure of the circuit increases. We verified Equation 5.5 by noise
simulations using the OpAmp circuit of Figure 5.5.
The NF is increased from 6 to 7.5 dB given that GM is equal to 20mS. Note that it is also possible to apply
the negative conductance in combination with an LNTA with higher GM and hence lower NF. In that
case, the negative conductance can be lower, as RO>1/GM. However, then IIP3 of the LNTA becomes a
bottleneck.
5.3 Measurement Results
Figure 5.7 shows a photo of the implemented 65nm IC. The active area is < 0.2 mm2 including the clock
circuit. Thick metal was used for RRF for high linearity and low spread.
Figure 5.7: Die Photograph (65nm CMOS, 1.45mm x 1.45mm)
69
The front-end achieves 26 dB gain (BalUn losses are de-embedded) at 1 GHz LO, over 24MHz bandwidth
(BW), 12MHz on either side of LO. To demonstrate distortion cancelling, Figure 5.8(a) shows the
measured in-band IIP3 at 150kHz tone spacing (f1=1004.1MHz and f2=1004.25MHz) vs. M. IIP3 clearly
improves from around +9 dBm to +21 dBm!
Figure 5.8: Measurements: (a) In-band IIP3 vs. the number of parallel Negative Conductance Unit-Cells M (b) IM3 versus
input power for three M settings, with LO=1GHz
The optimum IIP3 of
+21 dBm is located at
M = 32,
which fits
to our
theory
GTotal=1/RO+1/RF=1/1500+1/180=6.22mS so M=6.22mS/0.2mS=31 very well. Figure 5.8(b) shows the IM3
curves versus power for three cases: M=0 (off), M=28 (cancelling of IO, Figure 4.6) and M=32 (overall
optimum IIP3). Up to -22dBm input power (note: this power is high for an in-band signal), IM3 improves.
The rise of distortion for high input powers > -23 dBm is due to the clipping of the OpAmp output stage
to its 1.2V supply. The negative conductance was pushed to instability (i.e. latch-up of OpAmp output
stage). This occurs at M=45 (see Equation 4.8 ∆G=∆Mx0.2mS), safely away from the optimum point by
∆M=45-32=13. This shows a close agreement with our explanation in section 4.4 and with the
simulations in Figure 5.9, which is done for the circuit of Figure 4.13.
Figure 5.9: Latch-up simulation of VOUT, input power of -16 dBm
70
One tone input signal with power of -16 dBm is used. Around this input power, the OpAmp output stage
begins to clip. According to our simulation, the latch-up occurs for ∆M ≥ 14. The same mechanism,
discussed in section 4.1, of this technique also improves IIP2 by more than 10 dB as shown in Figure 5.10.
Figure 5.10: Measurements: IM2 versus input power for three M settings, with LO=1GHz
Table 5.1 compares/summarizes the IIP2 and IIP3 improvement for three M settings 0, 28 en 32. Note
that the optimum linearity point will vary somewhat with Process, Voltage and Temperature (i.e. PVT).
The analysis in this paper gives the relation between the required negative conductance and the
resistance values RO and RF, which can be a basis for designing an automatic PVT correction circuit.
Table 5.1: IIP2 and IIP3 improvement
71
Figure 5.11 provides IIP3 curves versus the frequency offset Δf, with fixed 3.95MHz in-band IM3 position.
The negative conductance clearly increases the IIP3 both in- and out-of-band (all-Band) with a worst case
IIP3 >+10 dBm.
Figure 5.11: 2-tone IIP3 measured at IM3=3.95MHz versus tone spacing Δf, with LO=1GHz
The reason behind less linearity improvement in the transition band can be understood considering the
equivalent circuit earlier derived for stability analysis in Figure 4.12. The negative conductance cancels
only the loading of RO and RF. However, gm(s), CF and CT introduce frequency dependences.
Consequently, the “loading effect” on the VGND node (see Figure 4.5) becomes frequency dependent
and will introduce a phase shift compared with the (frequency independent) current generated by the
negative conductance. This results in imperfect cancellation, i.e. less linearity improvement at high
frequencies. This may be improved in the future by designing the negative conductance to be frequency
dependent as well. Up to 10MHz, in-band IIP3 is >+20dBm, i.e. >10dB improvement thanks to the
negative conductance. Then the IIP3 declines from 12MHz to 135MHz, on the one hand because the OTA
gain and hence its linearity degrades, but on the other hand also because the benefit from cancellation
drops (the top line in Figure 5.11 drops faster, versus Δf, than the bottom line). Note that the out-ofband IIP3 at Δf > 450 MHz is again high, +18 dBm. This is because at high Δf (i.e. spacing between the
carriers) the carriers are filtered due to the low pass filtering by CF, RF and CO, hence less IM3 products. In
this region the negative conductance doesn’t result in any benefit anymore.
72
Figure 5.12: Compression point
The compression point (CP) is around -13 dBm (hardly affected by M as shown in Figure 5.12). Due to
the virtual ground, S11 is hardly affected by the negative conductance and Figure 5.13(a) shows that
S11 < -25dB. Noise is more worrisome, but depending on the application some degradation may be
acceptable, provided that the overall SFDR still improves (i.e. IIP3 in dBm should improve more than NF
in dB degrades). Figure 5.13(b) shows that NF increases from 6.2 dB at M=0 to 7.5 dB at M=32. This
result is close to the NF prediction in the previous section. The 1/f corner was around 2MHz.
Figure 5.13: Measurements (a) S11 (b) Noise Figure, with LO=1GHz
The current consumption without the negative conductance at 1 GHz LO is 18 mA (including 8mA of
clock circuitry (i.e. on-chip drivers and divider)), and 1.6 mA more for M=32. The clock divider frequency
range (i.e. also the receiving RF frequency) is 0.2-2.6 GHz, where it consumes 2.8-19 mA. The maximum
Gate-Source voltage of the mixer switches is equal to the 1.2V supply. The LO leakage to the RF port is
less than -75 dBm. The optimum IIP3 has been measured for 5 samples. The optimum in-band IIP3 varies
±1 dB around +21 dBm and the corresponding M varies ±2 around M=32.
73
5.4 Benchmarking
Equation 5.2 benchmarks this work to other state-of-the-art receivers with high linearity and/or SFDR.
Our front-end is more linear than [[28],[30]] where active RF blocks are present. Even compared to the
mixer-first designs [[31],[32]] we achieve better in-band IIP3 while our SFDR in 1MHz of 85dB is the
highest.
Table 5.2: Summary of measurement results and comparison to other state-of-the-art receivers
Linearization
Technique
Matching
Mixer type
Baseband-stage
CMOS Techn.
Active Area
RF Frequency
Gain
In-band BW[1]
NF
In-band IIP3
SFDR @ 1MHz
bandwidth
Wide-Band IIP3
@2-tone Δf
Supply Voltage
Power
Consumption
This work
Negative GO
Ru [28]
Partial cancel
Noise/Distortion
Murphy [29]
Cancel Noise
Youssef [30]
Freq. Translated
Active feedback
Soer [31]
Feedback +
N-path filter
Andrews [32]
Feedback +
N-path filter
Switch-R
Switch-R
TIA + RC
65nm
< 0.2
0.2-2.6
26.5
24
7.5
Common-gate
Switch-I
TIA+RC
65nm
<1
0.4-0.9
34
24
4
+3.5
Switch-R
Switch-R&I
TIA + RC
40nm
1.2
0.08-2.7
70
4
2
-22
R
Gm + Switched-I
Inverter-RC
65nm
< 0.06
1.0-2.5
30
5
7.25-8.9
-20
Switch-RC
Voltage Amp
65nm
< 0.13
0.2-2.0
19
50
6.5
+11
via TIA
Switch-RC
TIA+RC
65nm
0.75
0.1-2.4
40-70
1.6
4
-67
mm2
GHz
dB
MHz
dB
dBm
75
60
57
79
29
dB
+18 @ Δf>800
+13.5
@Δf>40
1.3
15.6
> +12
@ Δf>60
1.2
62
Not
measured
1.2
60
+25
@ Δf>50
1.2 / 2.5
< 70[2]
dBm
@ MHz
V
mW
> +20
85
≥+18 @ >450
>+10 @ All Δf
1.2
13.9
1.2
39.6
[1] In-band BW is twice the zero-IF bandwidth around the LO frequency
units
[2] Includes the clock circuitry
5.5 Conclusions
Due to the strong relationship between linearity and voltage swing, it is challenging to improve linearity
in advanced CMOS technologies with low supply voltages. Architectures with RF V-I conversion followed
by a passive mixers and an OTA-RC Transimpedance Amplifier perform relatively well. In such
architectures, the OpAmp can become the bottleneck, especially for wide channel bandwidth, where the
amount of loop gain available for negative feedback is limited. Still high linearity is wanted, not only outof-band but also in-band, as RF-filtering often is ineffective for close-in interferers. The measurements in
this chapter of our receiver proves that virtual ground imperfections due to OTA nonlinearity lead to
distortion currents, which can be cancelled exploiting a negative conductance in parallel to the virtual
ground node. Although the technique results in slightly degraded noise figure from 6 to 7.5dB the inband IIP3 (and IIP2) is improved by much more (>10dB), resulting in-band SFDR=85dB in 1MHz
bandwidth.
74
Chapter 6. Conclusions and Recommendations
In this chapter, the contents and key conclusions of this thesis are summarized. Also, original
contributions are mentioned and recommendations for future work are given.
75
6.1 Summary and Conclusions
Wireless communication experiences an enormous growth, while there are hardly any new frequency
bands left for dedicated use by new services and applications. Moreover, the traditional static frequency
planning which assigns each application to a specific frequency band is inefficient, as it leaves Radio
Frequency (RF) spectrum locally and temporarily unused, as discussed in chapter 1. Although cellular
network bands are used intensively in most parts of the world, this is not true for many other frequency
bands including the Digital TV broadcasting (DTV) spectrum band. To improve the spectral utilization
efficiency, Dynamic Spectrum Access (DSA) can help. New unlicensed services can be defined as
secondary users that dynamically and opportunistically exploit bands/channels that are locally or
temporarily not used by the licensed primary users. A Software Defined Radio (SDR) that is aware of its
frequency spectrum, i.e. a Cognitive Radio (CR), can be a good candidate to implement the DSA concept.
The DTV spectrum is a well-organized channelized RF system and the chosen frequency spectrum
environment of the thesis.
The DTV interferers will often lead to distortion in practical receivers with a limited degree of up-front
filtering and linearity. Those distortion levels may limit the frequency spectral utilization efficiency that
the CR concept requires to solve via DSA. In traditional narrowband RF receivers, out of band strong
interferers are greatly suppressed by a Surface Acoustic Wave (SAW) RF band filter after the antenna.
This will greatly reduce the distortion level of the received signal at the Analog-to-Digital Convertor
(ADC) input, i.e. the required out of band linearity of the receiver will be also reduced. The SAW RF filter
with a fixed center frequency does not help to immune the CR receiver from the DTV interferers that
actually exists in-band, hence the in-band linearity becomes a critical specification for CR receiver. The
ultimate solution is a flexible and reconfigurable RF channel filter with high rejection (i.e. very high order
filtering) of the interferers from the adjacent channel(s), note that guard channel(s) are not desired to
maximize the spectral utilization efficiency. Such RF channel filters are very difficult to be implemented,
especially on-chip. As a consequence of the poor rejection performance of RF channel filter, the in-band
linearity of the receiver must be increased to tolerate the received DTV interferers that actually exists inband.
This thesis aims to mitigate the linearity requirements challenges by benefiting from the presence of
spectrum sensing device, which is assumed to be on board. The main research questions are as follows:
76
How to analyze the distortion level (i.e. linearity requirement) across the unoccupied DTV
band for Dynamic Spectrum Access CR radio receiver?
Due to the existence of in-band strong DTV interferers and the nonlinearity of a radio receiver that is
modelled by a Taylor memory-less model (focusing on the dominant 3rd order nonlinear term), different
types of distortion products are generated, namely intermodulation (IM3), cross-modulation (XM3) and
self-interference products. This point is analyzed in chapter 2, where the following two important
conclusions are drawn:
•
The level and spectral location of the intermodulation distortion products depends on the level
and spectral location of DTV signals. Given statistics of DTV signal scenarios lead to a certain
probability of existence of a “white spot” with no intermodulation products.
•
Although intermodulation products are (by definition) absent in a white spot, the CR receiver is
any always still limited by cross-modulation and self-interference products.
Therefore, true distortion free white spots do not exist if DTV signals are present after the wideband
Band Pass RF Filter. However, cross-modulation and self-interference distortion products are typically
much weaker than intermodulation products. Thus it makes sense to monitor the level and spectral
location of interferers by spectrum sensing and classifying the white spots into two types, assuming only
third order distortion products:
•
“IM3-spots”
•
“IM3-free spots”
The derived equations of chapter 2 quantify how much the linearity requirements of 3rd order Input
Intercept Point (IIP3) are relaxed when the CR operates at an IM3-free spot. For example: a scenario with
three interferers, each has -10 dBm power and a CR signal that has -60dBm power. The calculated
linearity requirements at IM3-free spots is IIP3XM3=4dBm. While the linearity requirements at IM3-spots
where the receiver is majorly distorted by intermodulation, is IIP3IM3 =18-25 dBm. The analysis is verified
by measurements. The analysis of chapter 2 is extended to not only take narrowband interferers (e.g.
narrowband complex modulated carriers) into account but also wideband Orthogonal Frequency
Division Multiplexing (OFDM) interferers by adding a correction factor (i.e. ∆Wideband).
Reviewing the literature before this research started, the in-band linearity of a Complementary Metal
Oxide Semiconductor (CMOS) radio receiver is typically not much more than 0dBm. As this is by far not
77
enough to keep distortion below the noise floor, relaxing IIP3 requirements is very useful, e.g. by
predicting where IM3-free spots will occur.
How to predict the distortion level in Dynamic Spectrum Access, given Spectrum Sensing
data?
Chapter 3 targets a Distortion Prediction Algorithm (DPrA) based on Spectrum Sensing data. The DPrA
calculates the level and the spectral locations of the distortion products per each white spot.
Theoretically this requires a long and complex multiplications in time domain with complex envelop
functions of the wideband OFDM signals. The following three steps are done to develop the DPrA in its
simple processing form:
1. By using the correction factor (i.e ∆Wideband), the wideband signals are replaced by tones, where
only two information are required about the signals, namely frequency and amplitude.
2. By replacing the multiplications in time domain with convolutions in the frequency domain.
3. By implementing the convolutions via matrix multiplication, where the position and the value of
each element of the matrix corresponds to the spectral location and the amplitude of the tones
that composite the input signals to the nonlinear receiver.
The DPrA output is substituted in the equations of chapter 2 to derive the linearity requirements across
the white spots. Accordingly, the CR receiver is now able to select the appropriate white spot(s), which
have linearity requirements below the linearity capability of the CR receiver hardware. Due to the simple
matrix multiplication of the DPrA, it can be implemented in the Digital Signal Processing (DSP) section of
a receiver, so that instantaneous prediction can probably easily be provided after every spectrum
sensing operation. Moreover, DPrA can be easily extended to be applied for higher order nonlinear
terms (4th, 5th ... etc.).
Although, the DPrA allows the CR receiver to select the spots where a not so high linearity capability is
required, however this will limit the total utilization of the white spots. By increasing the CR receiver
linearity performance, more white spots can be accessed and still achieving sufficient signal to distortion
level, even the near/adjacent channels (e.g. IIP3IM3=25dBm is required in the adjacent channels of three
interferers, each has -10dBm) to the interferers.
78
How to increase the linearity of CR receiver hardware, given less up-front filtering?
In chapter 4, it is shown that a CMOS radio receivers that exploit linear Voltage-to-Current (V-I)
conversion at RF, followed by passive down-mixing and an OpAmp-based Transimpedance Amplifier at
baseband shows high linearity potential. Due to nonlinearity and finite gain in the OpAmp, virtual
ground is imperfect, inducing distortion currents. The concept of a negative conductance is proposed in
chapter 4 to cancel such distortion currents. Through a simple intuitive analysis, the basic operation of
the technique is explained. By mathematical analysis the optimum negative conductance value is
derived and related to feedback theory. It is shown that a finite value of negative conductance is needed
to make the feedback loopgain theoretically approach infinity, which is not possible by increasing the
gain of the OpAmp block. Furthermore, stability concerns are analyzed, considering both potential
(small-signal) oscillatory behavior and (large-signal) latch-up. From the small signal behavior perspective,
the system is stable. However, for strong negative conductance value the OpAmp output stage cannot
sink the negative conductance current and latch-up can occur. Fortunately, this latch-up condition
happens at far away from the optimum negative conductance value, as quantified in chapter 4. This
point can be improved by carefully designing the output stage of the OpAmp. Another relevant
important point is that the optimum linearity point will vary somewhat with Process, Voltage and
Temperature. Hence an automatic PVT correction circuit is important to circumvent the latch-up state.
Our analysis about the optimum negative conductance value gives a basis for such circuit design.
The technique is applied to linearize an RF receiver. Chapter 5 presents a prototype that is implemented
in 65 nm technology. Measurement results show an increase of in-band IIP3 from 9dBm to >20dBm, and
IIP2 from 51 to 61dBm, at the cost of increasing the Noise Figure (NF) from 6 to 7.5dB and <10% power
penalty. The NF analysis and simulation shows that the OpAmp thermal noise contribution is totally
cancelled at the optimum negative conductance value, where the OpAmp linearity limitation is also
totally cancelled. However the negative conductance contributes now noise. In- and out-of-band (i.e. allband) linearity is also analyzed and measured, showed best performance in-band and far out-of-band. In
the transition band, linearity improvement is degraded due to phase shifts between the frequency
dependent distortion currents generated by the OpAmp versus the frequency-independent current
generated by the negative conductance. In summary, the chip achieves for in-band, in 1MHz bandwidth,
a Spurious-Free Dynamic Range (SFDR) of 85dB at <27mA up to 2GHz for 1.2V supply voltage.
79
6.2 Original Contributions
The research in this thesis contains the following original contributions:
•
The analysis of the in-band linearity requirements (IIP3) for a wideband Cognitive Radio,
considering the wideband multi-carrier property of OFDM signal; (Chapter 2 and [18])
•
The classification of white spots in a channelized DTV band into IM3 spots and IM3-free spots, for
which the linearity requirement for a CR receiver is relaxed; (Chapter 2 and [18])
•
The development and the application of Distortion Prediction Algorithm (DPrA) to calculate the
3rd order nonlinear distortion products across the white spots, by using the spectrum sensing
information. Then the translation of DPrA outcome to evaluate the linearity requirement across
the white spots; (Chapter 3)
•
The design of a highly linear receiver, with record in-band IIP3 > +21dBm, by canceling the
induced distortion current via a negative conductance at the virtual ground node of the
baseband Trans Impedance Amplifier; (Chapter 4 and 5 and [35]and [41]).
6.3 Recommendations for Future Work
The ultimate hardware implementation for CR is Software Radio, where a fast high-resolution ADC with
low power consumption and high dynamic range is required. Tunable RF channel filters with a high
quality factor help to lower the dynamic range of the ADC. Till the time that those solutions will be
reached, it is reasonable to focus on SDR receiver imperfections and mitigating them via the data given
by the spectrum sensing device that is assumed to be on board. While the spectrum sensing provides
information about the level and frequency spectral location, the cross modulation effect of a blocker
inside the SDR receiver that is presented in chapter 2 may
also provide information about the instantaneous amplitude
level of the blockers. The analyses of Chapter 2 (i.e. paper
published in [18]) did not address the cross modulation
effect of a blocker on the CR signal constellations.
Traditionally it is known that as the noise power increases,
the “clouds” around the constellations points extend
Figure 6.1: Blocker cross modulation inside the
beyond the decision boundary, thereby causing errors in the
CR received channel
received CR signal. Referring to the equations in chapter 2, the intermodulation products have the same
80
effect on the constellation points as the noise. However, that is not true for the cross modulation
products. Let’s assume a blocker with less up-front filtering is received in addition to the desired CR
signal as shown in Figure 6.1. Then as mathematically is represented in chapter 2, the RF band filter
output signal, which is the input signal to the nonlinear receiver front-end, can be written as follows:
{
sin = Re 2 sCR (t) e jωCR t + 2 sBlocker(t) e jωBlockert
}
Equation 6.1
Whereas the output signal of the nonlinear receiver at the CR channel frequency is written as follows:
sout (t) CR Channel = k1 sin (t) + k 3 sin (t)3
Equation 6.2




2  jωCR t 
= Re
2k1 sCR (t) + 3 2 k 3 sCR (t) s Blocker(t)  e

142
4 43
4
444
4244444
3
Desired amplified CR 14

Cross modulation product 


While Equation 6.2 is already represented in chapter 2 (i.e. Equation 2.5), Equation 6.2 is now written in
its instantaneous form for one blocker. Observing the cross modulation product, it can be seen that the
instantaneous CR signal (i.e. sCR(t)) is distorted with the
instantaneous square amplitude of the blocker signal. This reveals
that for example if the CR signal uses QPSK digital modulation, then
the cross modulation will just affect the amplitude of CR signal by
moving the CR constellation points towards the original on a
diagonal path as shown in Figure 6.2. Measurement verifies this
observation. Therefore the CR receiver may instantaneously track
and subtract the blocker(s) amplitude variations from its received
Figure 6.2: QPSK CR signal
signal by the mean of time domain filtering.
Other recommendation points are related to the linearization technique presented in chapter 4 and 5:
•
The receiver linearity is degraded in the transition band (i.e. the band between in and out of
band in Figure 5.11) due to phase shifts between the frequency dependent distortion currents
generated by the OpAmp versus the frequency-independent current generated by the negative
conductance. This may be improved by designing the negative conductance to be frequency
dependent as well by applying negative capacitance.
•
The receiver noise is high. This may be improved by investigating the possibility of applying the
noise cancelling technique.
81
82
Appendix A.
Digital TV broadcasting (DTV)
Most of the countries in the world has already shifted towards the Digital TV and leaving the classical
analog TV due to the efficient use of spectrum and provision of more capacity than analog, better
quality images, and lower operating costs for broadcast and transmission.
Frequency bandwidth: The DTV is implemented in different ways around the world, as example:
•
The EU has 7-8MHz bandwidth with Orthogonal Frequency Division Multiplexing (OFDM)
modulation. It follows the Digital Video Broadcasting Terrestrial (DVB_T) open standard. DVB_T
offers three different modulation schemes (QPSK, 16QAM, and 64QAM), using 1705 or 6817
carriers (2k or 8k mode, respectively). The difference between the sub-channels is 4.46kHz.
•
The USA has 6MHz bandwidth with 8-level vestigial sideband (8-VSB) modulation. This
modulation is outside the scope of this book.
Frequency range: In USA, the FCC allocated three bands of frequencies in the radio spectrum, chopped
into 6-MHz slices, to accommodate the following TV channels:
•
VHF band: 54-88 MHz for Channels 2-6 and 174-216 MHz for Channels 7-13
•
UHF band: 470-890 MHz for Channels 14-83
Each country in Europe has its own allocation of the channels across the DTV frequency spectrum. A
typical frequency range for west European countries is listed as follows:
•
VHF band: 41-82 MHz and 175-224 MHz
•
UHF band: 471-951 MHz
Statistical Distribution reported in [42]: Major cities have many occupied TV-channels, including some
with high power levels. As example: in New York, simulated TV-signals reach levels around -15 dBm,
while the strongest measured signal is at -20 dBm. In Chicago, measurements and simulations show
several channels above -10 dBm, with the highest measured power at 0 dBm. At the other extreme,
measurements in the West-Virginia National Radio Astronomy Observatory (NRAO) radio-quiet zone
show a very low occupancy, and only a single channel with a power above -80 dBm.
83
Appendix B.
Wideband Correction Factor (∆Wideband) Effect on CR
Linearity Requirement
Equation 2.6 can be rewritten for the cross-modulation products as follows:
+
PXM3 = PXM3 + ∆ Wideband [dBm]
Converting it to the linear range gives:
+
p XM3 = p XM3 ×10
∆ Wideband
10
[W]
The shape of this distortion looks like a Bell and covers three spots. However we are just interested in
the amount of it in the middle spot, which contains 66% of the total distortion power:
+
XM3
p
∆Wideband
10
spot17
= pXM3×10
×0.66 [W]
Rewriting it according to its RMS values and substituting the cross-modulation product of Equation 2.5 in
it gives:
AXM3+ = AXM3×
∆Wideband
10 10
∆
Wideband


2
× 0.66 = 3k3ACR ∑ ATVn  × 10 10 × 0.66
 n=1,6,19



where AXM3+ and AXM3 represent the RMS value of the cross-modulation products of two wideband
signals and two tones, respectively. The next step is substituting it in the Signal to Distortion Ratio,
which gives:
S
=
D 
 3k A 
 3 CR 

(k 1A CR )2
∆ Wideband

2

10

A
×
10
×
0.66
∑
TV n 

n =1,6,19


2
Then the linearity requirements can be rewritten as follows:
k1
k3
= 3 × 10
spot 17
∆ Wideband
10
× 0.66 ×
S
2
A TV n
∑
D n =1,6,19
84
Appendix C.
DPrA Maple Program Implementation
Begin
> restart;
> with(linalg): with(ListTools): with(VectorCalculus): with(plots):
> alpha:=[k[1],k[2],k[3]]:
Forming Vector Sin
SCRn:=sqrt(2)*ACR/2: SCRp:=sqrt(2)*ACR/2: # Cognitive Radio signal
> STV1n:=S1*A[-1]: STV1p:=S1*A[+1]: S1:=0: # Channel 1
> STV2n:=S2*A[-2]: STV2p:=S2*A[+2]: S2:=1: # Channel 2
> STV3n:=S3*A[-3]: STV3p:=S3*A[+3]: S3:=0: # Channel 3
> STV4n:=S4*A[-4]: STV4p:=S4*A[+4]: S4:=0: # Channel 4
> STV5n:=S5*A[-5]: STV5p:=S5*A[+5]: S5:=0: # Channel 5
> STV6n:=S6*A[-6]: STV6p:=S6*A[+6]: S6:=1: # Channel 6
> STV7n:=S7*A[-7]: STV7p:=S7*A[+7]: S7:=0: # Channel 7
> STV8n:=S8*A[-8]: STV8p:=S8*A[+8]: S8:=0: # Channel 8
> STV9n:=S9*A[-9]: STV9p:=S9*A[+9]: S9:=0: # Channel 9
> STV10n:=S10*A[-10]: STV10p:=S10*A[+10]: S10:=0: # Channel 10
> STV11n:=S11*A[-11]: STV11p:=S11*A[+11]: S11:=0: # Channel 11
> STV12n:=S12*A[-12]: STV12p:=S12*A[+12]: S12:=0: # Channel 12
> STV13n:=S13*A[-13]: STV13p:=S13*A[+13]: S13:=1: # Channel 13
> STV14n:=S14*A[-14]: STV14p:=S14*A[+14]: S14:=0: # Channel 14
> STV15n:=S15*A[-15]: STV15p:=S15*A[+15]: S15:=0: # Channel 15
> STV16n:=S16*A[-16]: STV16p:=S16*A[+16]: S16:=0: # Channel 16
> STV17n:=S17*A[-17]: STV17p:=S17*A[+17]: S17:=0: # Channel 17
> STV18n:=S18*A[-18]: STV18p:=S18*A[+18]: S18:=0: # Channel 18
> STV20n:=S20*A[-20]: STV20p:=S20*A[+20]: S20:=1: # Channel 20
> NZbegin:=15: # Number of zeroes added at the begin of the vector
> NZend:=1: # Number of zeroes added at the end of the vector
> input:=vector(1,[0]):
> for i from 1 to NZend-1 do input:=augment(input,[0]): od:
> input:=augment(input,[STV20n]): # Adding Channel -20
> input:=augment(input,[STV19n]): # Adding Channel -19
> input:=augment(input,[STV18n]): # Adding Channel -18
> input:=augment(input,[STV17n]): # Adding Channel -17
> input:=augment(input,[STV16n]): # Adding Channel -16
> input:=augment(input,[STV15n]): # Adding Channel -15
> input:=augment(input,[STV14n]): # Adding Channel -14
> input:=augment(input,[STV13n]): # Adding Channel -13
> input:=augment(input,[STV12n]): # Adding Channel -12
> input:=augment(input,[STV11n]): # Adding Channel -11
> input:=augment(input,[STV10n]): # Adding Channel -10
> input:=augment(input,[STV9n]): # Adding Channel -9
> input:=augment(input,[STV8n]): # Adding Channel -8
> input:=augment(input,[STV7n]): # Adding Channel -7
> input:=augment(input,[STV6n]): # Adding Channel -6
> input:=augment(input,[STV5n]): # Adding Channel -5
> input:=augment(input,[STV4n]): # Adding Channel -4
> input:=augment(input,[STV3n]): # Adding Channel -3
> input:=augment(input,[STV2n]): # Adding Channel -2
> input:=augment(input,[STV1n]): # Adding Channel -1
> for i from 1 to NZbegin do input:=augment(input,[0]): od:
> input:=augment(input,[0]):
> for i from 1 to NZbegin do input:=augment(input,[0]): od:
> input:=augment(input,[STV1p]): # Adding Channel +1
> input:=augment(input,[STV2p]): # Adding Channel +2
85
> input:=augment(input,[STV3p]): # Adding Channel +3
> input:=augment(input,[STV4p]): # Adding Channel +4
> input:=augment(input,[STV5p]): # Adding Channel +5
> input:=augment(input,[STV6p]): # Adding Channel +6
> input:=augment(input,[STV7p]): # Adding Channel +7
> input:=augment(input,[STV8p]): # Adding Channel +8
> input:=augment(input,[STV9p]): # Adding Channel +9
> input:=augment(input,[STV10p]): # Adding Channel +10
> input:=augment(input,[STV11p]): # Adding Channel +11
> input:=augment(input,[STV12p]): # Adding Channel +12
> input:=augment(input,[STV13p]): # Adding Channel +13
> input:=augment(input,[STV14p]): # Adding Channel +14
> input:=augment(input,[STV15p]): # Adding Channel +15
> input:=augment(input,[STV16p]): # Adding Channel +16
> input:=augment(input,[STV17p]): # Adding Channel +17
> input:=augment(input,[STV18p]): # Adding Channel +18
> input:=augment(input,[STV19p]): # Adding Channel +19
> input:=augment(input,[STV20p]): # Adding Channel +20
> for i from 1 to NZend do input:=augment(input,[0]): od:
> input:=convert(input,vector): # Sin vector or A vector
> N:=vectdim(input): # Sin vedctor dimension is [ 1 x N ]
First convolution implementation
> C1:=2*N-1: # SQ vector dimension is [ 1 x (2*N-1) ]
> M1:=transpose(convert(input,matrix)):
> for i from 1 to (2*N-1)-N do M1:=augment(M1,[0]); od:
> MM:=transpose(M1):
> v:=convert(M1,vector):
> v:=convert(v,list):
> for m from 1 to N-1 do vv:=Rotate(v,-1*m); vv:=convert(vv,vector): MM:=augment(MM,vv): od:
> M1:=transpose(MM):
> firstconv:=evalm(input&*M1): # SQ vector
Second convolution implementation
> C2:=3*N-2: # The distortion vector dimension is [ 1 x (3*N-2) ]
> M2:=transpose(convert(input,matrix)):
> for j from 1 to (3*N-2)-N do M2:=augment(M2,[0]); od:
> MM:=transpose(M2):
> v:=convert(M2,vector):
> v:=convert(v,list):
> for m from 1 to (2*N-1)-1 do vv:=Rotate(v,-1*m): vv:=convert(vv,vector): MM:=augment(MM,vv): od:
> M2:=transpose(MM):
> Outputt:=evalm(alpha[3]*firstconv&*M2):
> Output:=vector(vectdim(Outputt)):
> for n from 1 to vectdim(Outputt)do Output[n]:=expand(Outputt[n]): od:
> for m from 1 to vectdim(Output) do C[m-round(vectdim(Output)/2)]:=evalm(Output[m]): od:
Viewing the Distortion Products
> for Cn from -(NZbegin+20) to -(NZbegin+1) do Channel[Cn+NZbegin]:=evalm(C[Cn]): od:
> for Cp from NZbegin+1 to NZbegin+20 do Channel[Cp-NZbegin]:=evalm(C[Cp]): od:
> for CH from NZbegin+1 to NZbegin+20 do Channel[CH-NZbegin]:=evalm(C[CH]): Channel[-H+NZbegin]:=evalm(C[-CH]): od:
> A[-2]:=sqrt(2)/2*ATV[2]: A[+2]:=sqrt(2)/2*ATV[2]:
> A[-6]:=sqrt(2)/2*ATV[6]: A[+6]:=sqrt(2)/2*ATV[6]:
> A[-13]:=sqrt(2)/2*ATV[13]: A[+13]:=sqrt(2)/2*ATV[13]:
> A[-20]:=sqrt(2)/2*ATV[20]: A[+20]:=sqrt(2)/2*ATV[20]:
> for CH from NZbegin+1 to NZbegin+20 do Channel[CH-NZbegin]:=evalm(2*C[CH]): od:
> ATV[2]:=ATV: ATV[6]:=ATV: ATV[13]:=ATV: ATV[20]:=ATV:
> for CH from NZbegin+1 to NZbegin+20 do Channel[CH-NZbegin]:=evalm(2*C[CH]): od:
86
Appendix D.
Nonlinear Transimpedance Amplifier Analysis
In this section, a 3rd order Taylor approximation of VOUT versus IS (i.e. VOUT=VOUT(IS,IS3)) of the
transimpedance amplifier shown in Figure C will be derived. The following procedure will be applied:
j.
VOUT is derived as a function of VVGND, VVGND3 and VOUT3
VOUT=VOUT (VVGND,VVGND3,VOUT3).
k. The resulting relationship is rewritten as a function of VVGND and VVGND3, by using the definition of
the 3rd order Taylor coefficients
VOUT=VOUT(VVGND,VVGND3).
l. The inverse function, VVGND as a function of VOUT and VOUT3, is written as a 3rd order Taylor
function by using the procedure explained in [40]
m. IS is rewritten as a function of VVGND and VOUT
VVGND=VVGND (VOUT,VOUT3).
IS=IS(VVGND,VOUT).
n. Substituting VVGND of step 3 in IS of step 4 makes IS to be a function of VOUT and VOUT3
IS=IS(VOUT,VOUT3).
o. Finally, by repeating the procedure explained in [40], the function of step 5 is inversed to obtain
VO as a function of IS and IS3
VOUT=VOUT (IS, IS3).
Figure D: Baseband model with RO and the extended RF for nonlinearity derivations (redrawn Figure 4.9)
87
VOUT=VOUT (VVGND,VVGND3,VOUT3): We begin the derivation by expressing the feedback current IF at
Step 1
the VGND node and the OUT node (see Figure D) as follows:
At VGND node: IF =
VVGND
+ GF−VGNDVOUT
RF−VGND
At OUT node: IF = −
VOUT
− GF−OUTVVGND
RF−OUT
Equation D.0
Equation D.1
Referring to the OpAmp nonlinear model, we equate IF in Equation 4.1 to IF in Equation D.1 as follows:
gm1VVGND + gm3 VVGND + go1VOUT + go3 VOUT = −
3
3
VOUT
− G F−OUTVOUT
R F−OUT
VOUT = −
(gm1 + G F−OUT ) V
Step 2
VOUT=VOUT(VVGND,VVGND3): VOUT is defined as: VOUT=β1VVGND +β2VVGND2 +β3VVGND3, which is a 3rd

1 
 go1 +

R
F
−
OUT
144
42444
3
a
VGND
gm3
−
VVGND −
3

1 
 go1 +

R
F
−
OUT
144
42444
3
b
go3
3

1 
 go1 +

R
F
OUT
−
144
42444
3
c
VOUT
order Taylor approximation around VVGND=0, where β1, β2 and β3 are the Taylor coefficients:
β n =1,2,3 =
1  ∂ n VOUT 
n!  ∂ VVGND n 
VVGND =0
To derive β1, we differentiate Equation D.2 with respect to VVGND as follows:
∂ VOUT
∂ VVGND
= a + 3bVVGND + 3cVOUT
 ∂ VOUT
∴β1 = 
 ∂ VVGND
2
2
∂ VOUT
∂ VVGND
∂ VOUT
a + 3bVVGND
⇒
=
2
∂ VVGND
1 − 3cVOUT
2

(gm1 + G F−OUT )

=
a
=
−


1 
 VVGND=0
 go1 +

R F−OUT 

The same procedure is used to derive β2 and β3:
1  ∂ 2 VOUT 

β 2 = 
=0
2  ∂ VVGND 2 
VVGND = 0
and
1  ∂ 3 VOUT 

β 3 = 
= b + a3 c
6  ∂ VVGND 3 
VVGND = 0
88
Equation D.2
(
)
VOUT = {
a VVGND + b + a 3 c VVGND
1
424
3
β
1
3
Equation D.3
β3
VVGND=VVGND (VOUT, VOUT3): We write the inverse of Equation D.3 in the Taylor series form: VVGND
Step 3
=α1VOUT+ α2VOUT2 +α3VOUT3. Deriving α1, α2 and α3 can be done by the procedure below.
First, let’s substitute Equation D.3 into its abovementioned inversed form as follows:
(
VVGND = α 1 β 1 VVGND + β 3 VVGND
3
)+ α (β V
2
1
VGND
+ β 3 VVGND
)
3 2
(
+ α 3 β 1 VVGND + β 3 VVGND
)
3 3
By equating the right to the left side of the equation above [40], the coefficients α1, α2 and α3 are
derived:
(
)
1
b + a3 c
3
VVGND = VOUT −
VOUT
4
a
a 43
{
142
α1
Step 4
Equation D.4
α3
IS=IS(VVGND,VOUT): Referring to IS in Figure D, we substitute the IF Equation D.0 at the VGND node
in the following equation:
 1

1
VVGND + G F− VGND VOUT
IS = I O + I F = 
+
 R O R F − VGND 
Equation D.5
IS=IS(VOUT,VOUT3): By substituting Equation D.4 into Equation D.5, the following equation is
Step 5
obtained:
(
1  1


1
b + a3 c


IS =  
+
 + G F− VGND  VOUT − a 4
 a  R O R F−VGND 

Step 6
VOUT =
) 

1
1
 VOUT 3
R + R
F − VGND 
 O
Equation D.6
VOUT=VOUT (IS, IS3): Finally, by inversing Equation D.6, the following expression is obtained:
 1
1
NL 3 
+
R
R
F − VGND
 O



1
3
IS +
IS
4
1  1


1
1  1


1
 + G F− VGND 
+

 + G F− VGND 
 
+

R O R F− VGND 
R O R F− VGND 
 a4 4

 a4 4
3
1
444
42 4 4 4 4 44
3
1
4 4 4424 4 4 4 4 4
Ω1
Where: NL3 =
Ω3
(b + a c) is related to the nonlinear terms of the OpAmp model.
3
a4
89
Equation D.7
Appendix E.
Latch-Up Analysis
As shown in Figure D, the negative conductance injects current via RF (i.e. VVGNDGLatch-up-Risk) that needs to
be handled by the OpAmp output stage in addition to the main current coming from GM (i.e. IS):
I OpAmp- Latch- up - Risk = IS + VVGNDG Latch − up − Risk
Equation E.0
Hence, the relation between VVGND and IS is the next step in this derivation.
Figure E: Latch-up problem at the OUT node (redrawn Figure 4.13)
In order to simplify the latch-up analysis, let’s a linear OpAmp (i.e gm3=go3=0). Consequently, Equation
D.2 and Equation D.7 can be simplified as follows:
1
VVGND = VOUT
a
Equation E.1
VOUT = Ω1IS
Equation E.2
By combining the previous two equations, the following relation between VVGND and IS is derived:
VVGND =
Ω1
1
IS =
IS
a
 1


1
 + a G F−VGND 
+

R
R
O
F
−
VGND




Equation E.3
After that the negative conductance cancels the loading effect of RO on the VGND node, it injects current
via RF that needs to be handled by the OpAmp output stage (see Figure D and Figure 5.5). Now if the
negative conductance becomes too strong then the potential latch-up becomes a real risk. For the case
of latch-up, Equation E.3 can be further elaborated to obtain the following equation:
90
VVGND =
1



1
− G Latch −up  + a G F−VGND 


 R F−VGND

Equation E.4
IS
Finally, Equation E.4 is substituted in Equation E.0 to obtain the latch-up equation:
IOpAmp-Latch-up-Risk = IS + VVGND G Latch−up−Risk
= IS +
1



1
− G Latch −up−Risk  + a G F−VGND 


 R F−VGND

91
ISG Latch −up−Risk
Equation E.5
92
References
[j]
National Telecommunications and Information Administration (NTIA), available online:
http://www.ntia.doc.gov
[k] M. Reardon. (2010, Jan.) FCC chairman pushes policy agenda, available online: http:
//www.cnet.com/8301-31045_1-10430991-269.html.
[l] M. Lazarus, “The great spectrum famine,” IEEE Spectrum, vol. 47, no. 10, pp. 26–31, October
2010.
[m] Shared Spectrum Company (SSC), “Spectrum occupancy measurements on different locations at
different time lines”, available online: http://www. sharedspectrum.com.
[n] M. McHenry, P. Tenhula, D. McCloskey, D. Roberson, and C. Hood, “Chicago Spectrum
Occupancy Measurements & Analysis and a Long-term Studies Proposal”, Proceedings of the
First International Workshop on Technology and Policy for Accessing Spectrum, USA, 2006.
[o] R. B. Bacchus et al., “Long-term, wide-band spectral monitoring in support of dynamic spectrum
access networks at the IIT spectrum observatory,” in Proc. 3rd IEEE Symposium on New
Frontiers in Dynamic Spectrum Access Networks (DySPAN), Chicago, IL, USA, 14–17 Oct. 2008.
[q] Federal Communications Commission (FCC), “In the Matter of Unlicensed Operation in the TV
Broadcast Bands, Additional Spectrum for Unlicensed Devices below 900 MHz and in the 3 GHz
Band”, Technical Report, November 2008.
[r]
R. Tandra, S.M. Mishra, A. Sahai, “What is a spectrum hole and what does it take to recognize
one?”, Proc. IEEE, vol. 97, no. 5, pp. 824–848, May 2009.
[s] IEEE, “IEEE standard definitions and concepts for dynamic spectrum access: Terminology relating
to emerging wireless networks, system functionality, and spectrum management (IEEE Std
1900.1-2008),” Sep. 2008.
[jt] J. Mitola, “The software radio architecture”, IEEE Communications Magazine, vol. 33, no. 5, pp.
26–38, May 1995.
93
[jj] S. Subhan, “Cognitive Radio transmitter with a broadband clean frequency spectrum”, PhD
thesis,
The
Netherlands,
Twente
University,
2014,
available
online:
http://icd.ewi.utwente.nl/publications.
[jk] B. Murmann, “A/D converter trends: Power dissipation, scaling and digitally assisted
architectures”, IEEE Custom Integrated Circuits Conference, September 2008.
[jl] J. Mitola, Jr.G.Q. Maguire, “Cognitive radio: making software radios more personal”, IEEE
Personal Communications, vol. 6, no. 4, pp. 13–18, Augustus 1999.
[jm] J. Mitola, “Cognitive radio”, Licentiate thesis, KTH, Royal Institute of Technology, Sweden,
September 1999.
[jn] B. Razavi, “RF Microelectronics”, Prentice Hall, New Jersey, 1998.
[jo] A. Ghaffari, E.A.M. Klumperink, M.C.M.Soer, B. Nauta, “Tunable High-Q N-Path Band-Pass
Filters: Modeling and Verification" IEEE Journal of Solid-State Circuits, May 2011.
[jq] M. Darvishi, R.A.R. van der Zee, E.A.M. Klumperink, B. Nauta, “Widely Tunable 4th Order
Switched Gm -C Band-Pass Filter Based on N-Path Filters " IEEE Journal of Solid-State Circuits,
December 2012.
[jr] D. H. Mahrof, E. A. M. Klumperink, J. Haartsen and B. Nauta, "On the effect of spectral location
of interferers on linearity requirements for wideband cognitive radio receivers", IEEE
Symposium New Frontiers in Dynamic Spectrum Access Networks (DySPAN), pp. 1-9, April 2010.
[js] P.F. Marshall, “Dynamic Spectrum Management of Front End Linearity and Dynamic Range”,
IEEE Symposium New Frontiers in Dynamic Spectrum Access Networks (DySPAN), pp. 1-12, April
2008.
[kt] W. Sansen, “Distortion in elementary transistor circuits”, IEEE Transactions on Circuits and
Systems II: Analog and Digital Signal Processing, vol. 46, no. 3, pp. 315-325, March 1999.
[kj] C. van den Bos, M. Ksuwenhoven, and W. Serdijn, “Effect of smooth nonlinear distortion on
OFDM symbol error rate”, IEEE Transactions on Communications, vol. 49, no. 9, pp. 1510 – 1514,
September 2001.
94
[kk] Q. Zou, M. Mikhemar, and A.H. Sayed, “Digital Compensation of Cross-Modulation Distortion in
Software-Defined Radios”, IEEE Journal of Selected Topics in Signal Processing, vol. 3, no. 3, pp.
348-361, June 2009.
[kl] Y. Li, J.M. Rabaey, and A. Sangiovanni-Vincentelli, “Analysis of Interference Effects in MB-OFDM
UWB”, Wireless Communications and Networking Conference 2008 (WCNC), Las Vegas,
November 2008.
[km] M.S. Oude Alink, E. A. M. Klumperink, M.C.M. Soer, A.B.J. Kokkeler,B. Nauta, "A 50MHz-to1.5GHz Cross-Correlation CMOS Spectrum Analyzer for Cognitive Radio with 89dB SFDR in 1MHz
RBW", IEEE Symposium New Frontiers in Dynamic Spectrum Access Networks (DySPAN), pp. 1-9,
April 2010.
[kn] Samir S. Soliman and Mandyam D. Srinath, “Continuous and Discrete Signals and Systems”,
Prentice Hall, New Jersey, 1988.
[ko] M.S. OudeAlink, E.A.M. Klumperink, A.B.J. Kokkeler, Z. Ru, W. Cheng, B. Nauta, "Using
Crosscorrelation to Mitigate Analog/RF Impairments for Integrated Spectrum Analyzers", IEEE
Transactions on Microwave Theory and Techniques (MTT), March 2013.
[kq] 3GPP TS 36.104: "Evolved Universal Terrestrial Radio Access (E-UTRA); Base Station (BS) radio
transmission and reception", available online: http://www.3gpp.org.
[kr] Z. Ru, E.A.M. Klumperink, B. Nauta, “A Software-Defined Radio Receiver Architecture Robust to
Out-of-Band Interference”, ISSCC Digest of Technical Papers, pp. 230-231, February 2009.
[ks] D. Murphy, A. Hafez, A. Mirzaei, M. Mikhemar, H. Darabi, M.F. Chang, A. Abidi, "A BlockerTolerant Wideband Noise-Cancelling Receiver with a 2dB Noise Figure", ISSCC Digest of
Technical Papers, pp. 74-76, February 2012.
[lt] S.S.T. Youssef, R.A.R. van der Zee, B. Nauta, "Active Feedback Receiver with Integrated Tunable
RF Channel Selectivity, Distortion Cancelling, 48dB Stop-Band Rejection and > +12dBm
Wideband IIP3, Occupying < 0.06mm2 in 65nm CMOS", ISSCC Digest of Technical Papers, pp.
166-168, February 2012.
95
[lj] M.C.M. Soer, E.A.M. Klumperink, Z. Ru, F.E. van Vliet, B. Nauta, "A 0.2-to-2.0GHz 65nm CMOS
Receiver without LNA Achieving >11dBm IIP3 and <6.5dB NF", ISSCC Digest of Technical Papers,
pp. 222-223, February 2009.
[lk] C. Andrews, A.C. Molnar, “A Passive Mixer-First Receiver With Digitally Controlled and Widely
Tunable RF Interface”, IEEE Journal of Solid-State Circuits, vol. 45, no. 12, pp. 2696-2708,
December 2010.
[ll] B. Razavi, “Design of Analog CMOS Integrated Circuits”, McGraw-Hill, Singapore, 2001.
[lm] S.C. Blaakmeer, E.A.M. Klumperink, D.M.W. Leenaerts, B. Nauta, "The Blixer, a Wideband BalunLNA-I/Q-Mixer Topology", IEEE Journal of Solid-State Circuits, vol.43, pp.2706,2715, December
2008.
[ln] D. H. Mahrof, E. A. M. Klumperink, M. S. Oude Alink, and B. Nauta, “A receiver with in-band
IIP3>20dBm, exploiting cancelling of OpAmp finite-gain-induced distortion via negative
conductance”, IEEE Radio Frequency Integrated Circuits Symposium (RFIC), pp. 601–604, June
2013.
[lo] J. Deguchi et al, “A Fully Integrated 2x1 Dual-Band Direct-Conversion Mobile WiMAX Transceiver
With Dual-Mode Fractional Divider and Noise-Shaping Transimpedance Amplifier in 65 nm
CMOS”, IEEE Journal of Solid-State Circuits, vol. 45, pp. 2774 - 2784, December 2010.
[lq] Z. Ru, “Frequency Translation Techniques for Interference-Robust Software-Defined Radio
Receivers”, PhD thesis, The Netherlands, Twente University, 2009, available online:
http://icd.ewi.utwente.nl/publications.
[lr] F. Bruccoleri, E.A.M. Klumperink, and B. Nauta, ”Noise Cancelling in Wideband CMOS LNAs",
ISSCC Digest of Technical Papers, pp. 406-407, February 2002.
[ls] W. Redman-White, D.M.W. Leenaerts, “1/f Noise in Passive CMOS Mixers for Low and Zero IF
Integrated Receivers”, ESSCIRC, pp. 18-20, September 2001.
[mt] E.A.M. Klumperink, “Transconductance Based CMOS Circuits”, PhD thesis, The Netherlands,
Twente University, 1997, available online: http://jive.el.utwente.nl/home/erick/Klumperink-PhDThesis.pdf.
96
[mj] Dlovan Mahrof, Eric Klumperink, Zhiyu Ru, Mark Oude Alink, Bram Nauta,"Cancellation of
OpAmp Virtual Ground Imperfections by a Negative Conductance applied to improve RF
Receiver Linearity", IEEE Journal of Solid-State Circuits, May 2014. (Invited paper)
[mk] M.S. Oude Alink, “RF Spectrum Sensing in CMOS Exploiting Crosscorrelation”, PhD thesis, The
Netherlands, Twente University, 2013, available online: http://icd.ewi.utwente.nl/publications.
97
98
List of Publications
Dlovan Mahrof, Eric Klumperink, Zhiyu Ru, Mark Oude Alink, Bram Nauta, "Cancellation of OpAmp
Virtual Ground Imperfections by a Negative Conductance applied to improve RF Receiver Linearity"
IEEE Journal of Solid-State Circuits, May 2014. (Invited paper)
Dlovan Mahrof, Eric Klumperink, Mark Oude Alink, Bram Nauta, "A Receiver with in-band IIP3>20dBm,
exploiting Cancelling of OpAmp Finite-Gain-induced Distortion via Negative Conductance" IEEE Radio
Frequency Integrated Circuits Symposium (RFIC), Seattle, June 2013. (Best Student Paper Award – 1st)
Dlovan Mahrof, Eric Klumperink, Jaap Haartsen, Bram Nauta, "On the Effect of Spectral Location of
Interferers on Linearity Requirements for Wideband Cognitive Radio Receivers" 4th IEEE Symposium of
New Frontiers in Dynamic Spectrum Access Networks (DySPAN), Singapore, April 2010.
99
100
Acknowledgments
In my opinion, the Doctor of Philosophy (Ph D) is not just a work that needs to be done. It
affects/changes your life, your way of thinking to extend that you can never imagine till the point of
writing the thesis. Finding a gap in applied techniques by other researchers is trivial as nothing is perfect.
Being creative and viewing that gap in different ways than the standard way, which is written in books
and papers, is difficult and simple. It is difficult because we are actually biased around the known
knowledge. However it is simple if you just dare to take the first step towards how you think about the
problem in your own way and compare it to the rest. I could never go through this process without the
devoted guidance, patient and kindness from my daily supervisor Eric Klumperink and my promotor
Bram Nauta. Their continue encouragements and support to me and other colleges is really amazing.
Thanks for the technical support from Gerard Wienk and Henk de Vries . They helped me in solving a lot
of issues while I was designing my first chip. They treated me with a lot of patience knowing the
stress/mood around Tap Out and measurements. Thanks to our Management assistant Madam Gerdien
Lammers for all her kind support and notes. Thanks go to my Project partners, Mark Oude Alink, Amir
Ghaffari and Saqib Subhan.
I will never forget the nice talks and discussions with Mustafa, Milad, Shadi, Kasra, Ramen, Michiel,
Harish, Jasper, Erik, Harish, Haifeng, Remko and Johan.
Thanks Zu Ru, you really gave me the chance to learn from your experience in IC design. I really miss our
nice discussions about IC design at the late nights in EL/TN building.
Thanks to the Dutch Technology Foundation STW (i.e. the applied science division of the NWO, and the
Ministry of Economic Affairs Technology Program) for funding this research.
I thank STMicroelectronics for silicon donation and CMP, and of course a special thanks goes to Andreia
Cathelin for her kind advices through the receiver chip design in 65nm ST technology.
Dr. Anna-Johan Annema, Rien van Leeuwen and Dr. Ronan van der Zee brought my first steps into this
wonderful world of electronics at my bachelor study phase. Thanks a lot!!!
Prof. Dr. Frank van Vliet and Mr. Joe Tauritz opened my eyes to see the world in Electric and magnetic
waves and how to control/feel those waves with very basic and principle way. Thanks a lot!!!
101
I would like to take this chance to thank Prof. Dr. Miko Elwenspoek for teaching me to go beyond the
equations and inter the real understanding of electrodynamics. His way of teaching inspired all of his
students to love physics and material science.
My special thank goes to Mr. Reint Brink the study advisor of our faculty EWI and Mr. Said ben Ayad the
studentendecaan. They believed in me and guided me throw the whole study/life inside the University
of Twente. Their doors were always open. I can never pay them back for all their kindness and support.
Thanks, Thanks a lot to the University of Twente for giving me the chance in getting the knowledge, the
nice memories inside the classrooms and at the campus with my great ex-roommates in OchtendRust
building (Calslaan 38 – UT Campus), Gerard (& his girlfriend Aank), Pieter, Peter, Paul, Mark, Elco,
Wijnand (& his wife Nicolet), Jeroen, Marikka, Emile, Evo, Bas, Arnne, Johan, Marco, Pietra, Soublex, and
others. Thanks a lot!!!
From Bruco IC Company, I would like to thank Mr. Harry Diepenmaat, Mr. Edwin Veldman, Mr. David
Bradford, Madam Sharon Grunder and Mr. Charles Klaasen for their kind encouragements and support.
My special thank goes to Johan Koopman and Paul Bauis for their kind support in finishing my Ph D in
combination with my current work as RF IC design engineer in a flexible way. Of course, I would like to
thank all my colleges at Bruco IC, Frank, Marco H., Harry, Robert, Alexander, Eisse, Paul K., Paul V.,
Jurjen, Rinus, Hans B., Bram, Ellen, Erwin, Hans H., Jeroen, Jelle, Gerard, Rob, Marco W., Marc, Sander,
Tim, Tanya, Mirthe, Lizan, Wendy, Frans, Mariska, Michel, Marga, Caroline, Faust, Arie and Mathijs.
I would like to express my appreciation to Prof. Dr. Osama Shanaa from MediaTek USA. I learned a lot
from his comments on my RFIC 2013 paper, which have guided me through other aspects to be explored
in my JSSC2014 paper. Thanks a lot!!!
Finally, but foremost, I would like to thank my dear wife Almas for her support throughout the long
working hours and the associated ups and downs. She kept holding my hand in each step through all the
obstacles and made me to feel what is HOME and as you know West or East, HOME is the best. My real
apologies go to my dear son Varto to not be there in the first three years of his life. I can never bring the
days, months, and years back as I cannot bring one minute of the clock back. I hoop that I can teach you
from what I learned but foremost to leave you to enjoy the life as well. Last year, we completed our
family with our daughter flower Zhareen.
102
Biography
From left, Eric Klumperink, Bram Nauta and me at the IEEE Radio Frequency Integrated Circuits
Symposium (RFIC), Seattle, June 2013
Dlovan Hoshiar Mahrof was born on December 26, 1973, in Al-Sulaimaniya, Iraq. He received the B.Sc.
degree in power engineering from Baghdad University, Iraq, in 1995. He received the B.Sc. degree and
the M.Sc. degree in electrical engineering from Twente University, The Netherlands, in 2005 and 2008.
From 2008 to 2012, he did his Ph.D. research on “Distortion Mitigation in Cognitive Radio Receivers”
with the IC-Design group of Twente University, headed by Bram Nauta. During the RFIC 2013
Symposium, Seattle, WA, USA, he won the Best Student Paper Award—1st place.
From 2012 until the present, he has been working at Bruco IC Company, The Netherlands, as RFIC
system and circuit design engineer. His research interests include circuit design, and advanced
techniques applied to RF Front-Ends.
103
104