A Practical Guide To EMC Engineering (PDFDrive)
A Practical Guide To EMC Engineering (PDFDrive)
A Practical Guide To EMC Engineering (PDFDrive)
Levent Sevgi
Library of Congress Cataloging-in-Publication Data
A catalog record for this book is available from the U.S. Library of Congress.
All rights reserved. Printed and bound in the United States of America. No part of this book
may be reproduced or utilized in any form or by any means, electronic or mechanical, including
photocopying, recording, or by any information storage and retrieval system, without permission
in writing from the publisher.
All terms mentioned in this book that are known to be trademarks or service marks have been
appropriately capitalized. Artech House cannot attest to the accuracy of this information. Use of
a term in this book should not be regarded as affecting the validity of any trademark or service
mark.
10 9 8 7 6 5 4 3 2 1
Contents
Preface xiii
Introduction 1
v
vi A Practical Guide to EMC Engineering
2 Accreditation 33
2.1 Introduction 33
2.2 Accreditation 35
2.4 TÜRKAK 38
3 Electromagnetic Model 51
4 Circuit Model 65
4.5 Grounding 83
4.10 Crosstalk 99
References 100
Bibliography 101
Bibliography 289
Index 295
Preface
Every book is an accumulation of knowledge and/or experience.
The material for this book has been accumulated since 1996. At that
time, concern about cell phones, base stations, and public health was continu-
ously increasing at the same time as legislation of the necessity of Conformity
Europe (CE) marking (i.e., the Electromagnetic Compatibility Directive) on
every product entering the European Union market was taking effect. Health
concerns fall under bioelectromagnetics (BEM), while legislation is the subject
of electromagnetic compatibility (EMC) engineering. Both topics are equally
important. They cover a large number of subjects from public understanding of
science to problems of technological developments and from market control to
engineering education. In other words, they cover life itself.
During these years, while we were quietly studying in our labs or sitting in
front of our PCs, our phones started ringing. We were continuously bombarded
with questions from government offices, municipalities, chambers, and other
institutes and organizations as well as citizens, on EM pollution and possible
adverse health effects of base stations and cell phones. We did our best to an-
swer these (technical and nontechnical) questions. As time went on, our single
initial responses had turned into public seminars, panels, and radio and televi-
sion talks. We received a wide range of speaking invitations, from ministries
of health, environment, and industry to GSM companies, from military acad-
emies to medical faculties, and from schools to local associations. We wrote in
newspapers, magazines, and scientific journals. The aim was to provide people
with the right information.
In the same time slice, in September 1997, an older colleague of mine
knocked on my door and told me that in the next month, the Turkish Elec-
tronics Industry Association (TEIA) would organize two EMC training courses
at Istanbul Technical University. He said that nearly 60 attendees had already
xiii
xiv A Practical Guide to EMC Engineering
registered in each course; unfortunately, TEIA could not find a trainer for this.
I was going to ask “what is this EMC?” when he said that he gave my name as
the trainer because he thought I could teach these courses. You can imagine that
the next few weeks before the courses were a nightmare for me! Nevertheless, I
did my best and prepared myself and created course notes on time. Both courses
were highly successful and influenced many companies in Turkish industry.
Since then, I have repeated these two- to three-day EMC training courses at
many prestigious Turkish companies including Aselsan, Arçelik, Beko, Alcatel,
and Vestel, covering defense to information technology equipment (ITE) and
white goods producers. These courses enabled me to meet with engineers and
groups of people in different sectors from design to tests and measurements,
from planning to marketing. I have also organized these courses in universities
and research institutes open to personnel of small and midsize companies. There
were design/system/test engineers, quality control managers, directors, owners,
and even lawyers in these courses. There were also technical and nontechnical
representatives of ministries and municipalities, accreditation, and certification
authorities, who are responsible for market control. I have also learned much
more than I taught about the problems and challenges of different groups and
sectors. I have continuously updated and enriched my course notes. I published
three EMC-related books in Turkish.
This is what I tell the attendees at the beginning of these courses: “I know
you’re attending this course because you all face serious EMC problems. First,
I must say that that I’m not going to solve your problems; you’re going to do
that! Second, I’m not going to teach you something totally new. What I’ll do in
the following two to three days is mostly to fill in the gaps in your knowledge
and experience. I’m going to revisit everything you learned in your university
education and teach you how to look at them with an eye of an EMC engineer.
I’m quite sure that at the end of the course you will learn how to diagnose and
solve your problem; very often, you’ll suddenly realize where the problem is
during the course. That would be the success of the course.”
The title of this book is an interesting phrase that characterizes the worth
of EMC engineering. This book is unique in terms of its content. It is not a
specific EMC design book; it neither aims to be nor claims to be this! There are
several books on EMC (see the list of suggested readings) that aim to discuss
specific topics such as signal integrity, grounding, shielding, anechoic cham-
bers, and so on. The book covers essentials in EMC engineering that necessi-
tates a different look into many EM concepts and issues. As the title implies, it
contains all the procedures from design to market, including technical and non-
technical issues, which mostly have not been completely covered by any other
EMC books, such as market control, accreditation, calibration, EMC tests and
measurements, and EMC protection. There are many sections and subsections
as well as practical tips related to design. In addition, there are nearly two dozen
Preface xv
MATLAB scripts useful for EMC engineers. The purpose of giving these MAT-
LAB scripts is not to attack real-life EMC problems, but to establish the first
(most valuable) step for the beginners.
There are 10 chapters in the book. The logic of the book is as follows:
First, fundamental definitions, terms, and concepts, including processes in
achieving EMC compliance through testing and certification, are given in
Chapter 1. Accreditation is summarized in Chapter 2. Both EM and circuit
theories are revisited through the eyes of an EMC engineer in Chapters 3 and 4.
Chapter 5 is reserved for antenna and antenna calibration. Noise and frequency
analysis are given in Chapter 6. Chapter 7 focuses on EMC tests and measure-
ment environments and their calibration. Chapter 8 belongs to basic EMC
tests and measurement equipment. EMC tests and measurements methods and
implementation are presented in Chapter 9. Finally, major protection methods,
such as filtering, shielding/screening, grounding and cabling, and impedance
matching, are included in Chapter 10.
Note that modeling and numerical simulation in EMC engineering are
important topics that should be addressed in a separate book. There are many
sources on the Internet on these topics. Many commercial EM packages (e.g.,
CST, FEKO, COMSOL, EMA, REMCOM, ANSYS, SOLIDWORKS) have
continuously been modified to handle more and more complex EMC/EMI
problems. For example, several EMC/electromagnetic interference (EMI) mod-
eling and numerical simulation examples, from EM and lightning protection
to shielding, aircraft chambers to cabling for data and power lines, compos-
ite materials to EMC in radar systems, ESD to system-in-chip simulations,
which have been prepared using the commercial computer simulation technol-
ogy (CST) packages, may be found at www.cst.com/Applications/Category/
EMC. Another example is the site jointly sponsored by the IEEE EMC-S TC9
and ACES, www.ewh.ieee.org/cmte/tc9, where excellent practical examples of
EMC/EMI modeling, simulation, and validation with a focus on canonical,
benchmark, and standard problems are presented.
I have spent nearly three decades on EM and studied on topics such as
propagation through complex environments, radars and integrated surveil-
lance systems, EM scattering and diffraction, radar cross section (RCS) pre-
diction and reduction, EMC-BEM, antennas, EM modeling and simulation,
microstrip circuits, and so on. These topics certainly necessitate establishing
an intelligent balance between strong mathematical background (theory), en-
gineering experience (practice), and modeling and numerical computations
(simulation). Generally speaking, I know that EMC Society (EMC-S) people
find books written by Antennas and Propagation Society (AP-S) authors too
academic (theoretical), and that books written by EMC-S authors are found by
AP-S members to be too industrial (practical). Being a member of both IEEE
AP-S and EMC-S, I aimed, in this book, to balance theory and practice. As can
xvi A Practical Guide to EMC Engineering
be seen from the topics listed in the table of contents, they are the essentials of
EMC engineering that are mostly not covered by other published EMC books.
I believe this book will be a reference book for broad range of EMC-re-
lated technical and nontechnical people and can be used as a textbook in many
introduction and advanced-level EMC lectures. You want to have this in your
library! I hope you enjoy reading it and find it useful.
1
Introduction
The role of electromagnetic (EM) fields in our lives has been increasing [1–6].
Communication, remote sensing, integrated command/control/surveillance
systems, medicine, environment, education, marketing, and defense are only
a few areas where EM fields have critical importance. In 2015, the number of
first-time cell phone subscribers (7.5 billion) exceeded the world population
(7.4 billion) and the number of Internet users exceeded 3 billion. We are sur-
rounded by EM fields.
Figure 1.1 shows a typical recorded spectrum from DC to 1 GHz, show-
ing EM signals in our environment. These signals and levels can be recorded
anywhere on the Earth with quite similar spectra. Today, the spectrum is ex-
panded to 6 GHz.
All of these necessitate the design and production of less interfering, less
susceptible electronic devices on one hand; on the other hand, they create en-
gineering problems in terms of EM fields—human interactions. Two new en-
gineering disciplines were born: electromagnetic compatibility (EMC) and bio-
electromagnetics (BEM). EMC and BEM engineering, in general, deal with
identification of, respectively, adverse device-device and device-human inter-
actions (usually unintentional) and understanding the mechanisms occurring
there as well as their mitigation. Conventional EM engineers are concerned
with a range of product specific issues, whereas EMC engineers have to deal
with all possible external EM influences on the environment. They are both
multidisciplinary and cover electrical, electronics, communication, system en-
gineering and chemistry, physics, medicine, and biology and both have design,
test, production, quality, marketing, and legal implications.
At the same time, governments have tightened controls for the products
bought and sold in their markets. Military and aerospace systems must also
1
2 A Practical Guide to EMC Engineering
• 3 Hz–3 kHz range: Extremely low frequencies. (AC mains is the major
EMC problem in this region; high voltage power lines are the major
concern in BEM engineering.)
• 3 kHz–30 kHz range: Very low frequency (VLF) region.
• From a few hundred kHz up to 1012 Hz: Radiowaves.
• Around 1012 Hz: Infrared region. Everything radiates in this region.
For example, nearly 60% of our body heat is radiated out as infrared
energy. Infrared sensors (i.e., thermal cameras) are used in night vision
(detection in the dark).
• Visible light around 1014 Hz.
Electric and magnetic fields are defined as force lines surrounding elec-
trical devices. The unit of electric fields caused by voltage difference is [V/m]
and [A/m] for the current generated magnetic field. Electric fields are caused
by electric charges; magnetic fields are caused by moving charges (i.e., electric
currents).
The electric field is the force act on a charge. The magnetic field is the
force act on moving charges.
Electric and magnetic fields are vector quantities; they have both ampli-
tude and direction. One way of representing magnetic fields is magnetic flux
density (B); its unit is tesla [T] or gauss [G] (1T = 10000G). Magnetic fields
and magnetic flux density are related with the magnetic permeability (µ) of the
medium (B = µH). In free-space, in air and in living tissues it is µ = µ0 = 4π10–7
[Henry/m].
In general, DC magnetic fields do not induce body currents on humans.
But weak currents that can affect hearing devices and/or pacemakers may be in-
duced in industrial areas under very high DC magnetic fields. AC mains supply
electric power everywhere at 50/60 Hz frequency all around the world. Time
-varying EM fields cause EM energy absorption and in-body currents. The rela-
tion between the flowing current and electric field in a tissue is given by J=σE
[A/m2] where σ is the tissue conductivity [S/m].
We are surrounded in our environment by static (DC) and dynamic (AC)
electric and magnetic fields. Earth itself creates fields. Some typical values are:
• In air, under normal conditions, DC electric fields are around 120–150
V/m.
• On rainy days with severe thunderstorms and lightning, DC electric
field values may reach up to 10 kV/m.
• DC geomagnetic flux of Earth is around 0.5G; its AC value is around
10-5G.
• In our bodies, the DC magnetic flux is around heart 10-7G, in our ab-
domen 10-6G, and in our lungs 10-5G.
Introduction 7
ary principle lower their reference levels. For example, Italy uses 20 V/m
(at 900 MHz) value as a general level and further reduces it to 6 V/m
around schools, hospitals, and so on. Russia has 6 V/m and Sweden has
4 V/m at 900 MHz as reference levels.
Figure 1.4 shows ICNIRP reference levels versus frequency. These are ob-
tained as follows [8]:
The top and bottom curves in the figure correspond to occupational and
public exposures, respectively. The general public exposure limit is 300 V/m
from DC to 1 MHz, but it is 42 V/m at 900 MHz cell phone frequency band.
Note that the ICNIRP guides do not specify EM emission levels of the
devices we use in our daily lives. They are derived under normal conditions for
10 A Practical Guide to EMC Engineering
average humans [8]. They are given according to available scientific findings,
are dynamic, and are revised as more and more scientific evidence appears.
1.2.4 Epilogue
EM fields and public concern is a long-term, sensitive issue [9]. EM pollution
discussion has been raised in the public all around the world parallel to expo-
nential increase in cellular phone users. Because of the lack of information and
long-term experimental and epidemiological studies, contradictory opinions
exist on EM exposure safety levels in societies. Different governmental offices,
universities, media, and nonprofit organizations may have different approaches
to this problem. Since available scientific knowledge is far from being able to
end these discussions, this is normal. Beliefs may be totally different form one
group to another. It is the obligation of the scientists not to confuse beliefs with
scientific studies.
The debate on possible adverse effects has gained increasing attention.
Therefore, one needs to be careful even when choosing the right words. For
example, effect does not necessarily mean harm. Saying “no adverse effect has
been observed” does not mean “there is no adverse effect.” In order to say “there
is no . . . ” scientifically, all possible situations must be checked, which is almost
impossible in EM wave–human interaction.
The discussion about SAR (amount of EM power in watts/kg absorbed
by human tissues) is almost over. We now know that regulations based on SAR
are used in EM field measurements around base stations, and measured values
are well below specified threshold levels. The discussion has focused recently on
the question “What if very weak EM fields, well below the thresholds specified
in the regulations, trigger adverse effects in the long term (e.g., 10 years, 20
years)?” The answer will come only after long-term, accurate epidemiological
studies.
In May 2010, the results of the Interphone study [8] coordinated by
IARC were published. It was a 10-year multinational, epidemiologic study that
included nearly 13,000 human subjects from 13 countries around the world.
The study concerned the attempt to determine whether mobile phones increase
the risk of developing cancer involving the tissues around the ears and mouth in
humans. Final results showed no association between cell phone use (less than
10 years) and cancer.
In May 2011, the World Health Organization (WHO) [3] and Interna-
tional Agency for Research on Cancer (IARC) classified mobile phone use and
other RF EM fields as possible carcinogens (group 2B). Note that, IARC clas-
sifies substances to show whether or not they are suspected to cause cancer. It
places them into one of five categories depending on the strength of the evidence
for their carcinogenicity. These categories are Group 1 (carcinogenic to humans
Introduction 11
400 animals per study (50 animals per sex per group, with three dose groups
and a control) and at least two years of chronic study with complete necropsy
and histopathology [7]. In order to study nonthermal effects of cell phone radi-
ation, NTP requires (a) an ability to expose large numbers of rodents to a high
dose; (b) animals to be unconstrained and housed in standard laboratory cages;
(c) food and water to be available on demand; (d) excellent field and specific
absorption rate (SAR) homogeneity and detailed numerical and experimental
dosimetry; (e) the ability to discern a possible dose response; and (f ) third party
verification of the correct operation of the system. A study must satisfy these
criteria first (NTP spent multimillion US dollars for the establishment of a
RF radiation reverberation chamber exposure system for rodents). Note that, a
waste number of in vitro studies were conducted with at most 20–40 animals.
This is why WHO approaches their results with caution.
In view of the seriousness of the risk and the urgency of the situation, and
having regard to the objective of the decision, the Commission did not
act in a manifestly inappropriate manner by adopting the decision, on a
temporary basis and pending the production of more detailed scientific
information. Where there is uncertainty as to the existence or extent of
risks to human health, the Commission may take protective measures
without having to wait until the reality or seriousness of those risks
becomes apparent.
Pt
S= W m2 (1.1)
4 πR 2
The transmit antenna is usually directive with antenna gain Gt; therefore,
power density along the maximum radiation direction is:
Pt Gt
S= W m2 (1.2)
4 πR 2
Pt Gt λ2
Pr = Ae [W ] Ae = Gr m 2 (1.3)
4 πR 2 4π
2
Pt Gt Gr λ2 λ
Pr = = Pt Gt Gr [W ] (1.4)
2
4 πR 4 π 4 πR
Table 1.1
Basic Power Calculations in Communication Systems
Quantity Explanation Unit
Pt Transmit power Watt
Pr Receive power Watt
S Power density W/m2
Gt, Gr Antenna gains
R Transmit-receiver distance Meter
λ Wavelength Meter
f EM signal frequency Hz
Introduction 15
The power density at the receive point may also be found using
E2
S= W m2 (1.5)
120 π
if the electric field at the receive point is known. Combining (1.1)–(1.5) yields
30Pt Gt
E=
R
[ V m] (1.6)
Table 1.2
Typical Magnetic Fields Produced by Household Devices [mG]
Equipment d = 10 cm d = 30 cm d > 1m
Vacuum cleener 300–400 30–50 3–5
Iron 5–10 1 0.5
Washing machine 20–30 3–5 1
Hair dryer 400 10 1
Electrical shaver 200 5 0.5
Mixer 70–200 5–10 0.5
Dishing machine 20–40 5–10 0.5-1
Electric oven 5–10 0.5 0.5
Freezer 5–10 1 0.5
Music Set 3–5 1 0.5
TV 3–5 1 0.5
Digital watch 2–3 1 0.5
Air Conditioner 100 3–5 0.5
Photocopy machine 80–150 10–30 1–3
1992, after a three-year transition period. It delivers all necessary and sufficient
conditions that must be satisfied by the products. After years of feedback from
the market, its revised version, 2004/108/EC, was introduced. This new direc-
tive was based on the assumption of a healthy market with the most applicable
legal infrastructure. The new directive was published in the EU Official Journal
on December 31, 2004. On July 20, 2007, directive 2004/108/EC repealed the
long-standing EMC directive 89/336/EEC. A major difference is the change in
definition of apparatus to equipment. The term apparatus in the new directive
is now specific to devices or product only, and excludes fixed installations. The
term equipment covers both apparatus and fixed installations. It also simplified
the CE process. Finally, both component body and notified body referred to
the same institution. The last EU EMC directive was introduced in 2014 (EU
EMC directive 2014/30/EU). It covers products that are new to the EU market
and applies to all forms of supply, including distance selling. After April 2016,
the new directive was required for all EMC compliance files, and declarations
referencing 2004/108/EC were longer valid. In addition to those for produc-
ers and importers, it also brings obligations to law makers, market controllers,
and surveillance (i.e., to ministries and national accreditation institutions). The
EMC directive differs from most other directives in that its primary require-
ment is protection of the EM spectrum rather than safety of the equipment.
The EMC directive does not address safety!
The CE mark is neither a quality mark nor a guarantee that the product
meets all of the requirements of relevant EU product safety law.
• Basic standards: They define and describe the EMC problem, measure-
ment and test methods, principle measuring equipment, and test setup.
They contain no limits and state no interference criteria. Basic immu-
nity and emission EMC standards are specified by the IEC and CISPR,
respectively, (e.g., CISPR 16 and IEC 61000).
• Generic standards: These are standards that have to do with a particular
EM environment. They specify appropriate series of requirements and
tests that are used for all products and systems in this environment.
There are two major environments: (1) residential, commercial and light
industry; and (2) industry.
• Product family standards: These standards contain special limits for
emission and immunity for a specific category of devices (e.g., CISPR
22, information technologies). They contain many basic standards.
• Product-specific standards: They are for specific products (e.g.,
prEN50220 is the immunity standard for hearing aid products).
22 A Practical Guide to EMC Engineering
Table 1.3
EU EMC Standards
EN50081-1:1992— Electromagnetic compatibility—Generic emission standard—Part 1:
Residential, commercial, and light industry
EN50081-2:1994— Electromagnetic compatibility—Generic emission standard—Part 2:
Industrial environment.
EN50082-1:1998— Electromagnetic compatibility—Generic immunity standard—Part 1:
Residential, commercial, and light industry
EN50082-2:1995— Electromagnetic compatibility—Generic immunity standard—Part 2:
Industrial environment
EN55013:2001— Limits and methods of measurement of radio disturbance characteristics of
broadcast receivers and associated equipment
EN55014-1:2001— Electromagnetic compatibility—Requirements for household appliances,
electric tools, and similar apparatus—Part 1: Emission—Product family
standard
EN55014-2:1997— Electromagnetic compatibility—Requirements for household appliances,
electric tools, and similar apparatus—Part 2: Immunity—Product family
standard
Introduction 23
Table 1.4
IEEE/ANSI ASC-C63 EMC Standards Used in the US
C63.16-2016- ANS Guide for Electrostatic Discharge Test Methodologies and Acceptance
Criteria for Electronic Equipment
C63.7-2015- ANS Guide for Construction of Test Sites for Performing Radiated Emission
Measurements
C63.23-2012- ANS Guide for Electromagnetic Compatibility—Computations and Treatment
of Measurement Uncertainty
C63.7-2005- ANS Guide for Construction of Open-Area Test Sites for Performing Radiated
Emission Measurements
C63.5-1998- ANS for Calibration of Antennas Used for Radiated Emission Measurements
in Electromagnetic Interference (EMI) Control Calibration of Antennas (9 kHz
to 40 GHz)
C63.15-2010- ANS Recommended Practice for the Immunity Measurement of Electrical
and Electronic Equipment
C63.5-2006- ANS for Electromagnetic Compatibility Radiated Emission Measurements in
Electromagnetic Interference (EMI) Control Calibration of Antennas (9 kHz to
40 GHz)
In Figure 1.6, CE limiting voltage values, in dBµV, versus frequency for the
range of 150 kHz–30 MHz for both EN50081-1 and EN50081-2 standards
are shown. As observed, limiting CE voltage values vary between 40–60 dBµV
in residential regions, but it is between 65–80 dBµV in industrial areas.
Figure 1.7 shows FCC and CISPR RE limits as electric field (in dBµV/m)
versus frequency in the frequency range of 30 MHz–1 GHz. Here, Class A and
B, respectively, correspond to industrial and residential cases in EN standards.
The values in the figure belong to 10m distance between EUT and the EMI
receive antenna. As observed, FCC and CISPR values are similar but not ex-
actly the same.
Note that standardization is a dynamic process where EMC limiting values
and frequency ranges are continuously updated. For example, IEC/EN61000-
4-3 (1995) RI values were first given in the 80 MHz–1 GHz frequency range.
26 A Practical Guide to EMC Engineering
After 2001, tests were mandatory up to 2.5 GHz. Finally, in the 2006 issue of
IEC/EN61000-4-3 the maximum test frequency became 6 GHz.
These components are critical for EMC tests and measurements. For ex-
ample, effects and/or contributions of tests and measurement environments
shall clearly be known; otherwise, the difference between tests and measure-
ments conducted in a screened room and anechoic chamber shall not be un-
derstood and evaluated. Resonance effects of screened rooms, low and high
frequency contributions, and couplings in an anechoic chamber must be well
understood. Similarly, the isolation of the EUT and EMI receiver, the isolation
of mains, and the EUT must clearly be supplied and demonstrated before the
tests and measurements.
AC mains (power supply) is the worst interference source.
Introduction 27
EMC tests and measurements can be divided into two groups: emission
and immunity (or susceptibility).
Emission is measured; immunity/susceptibility is tested.
Both emission measurements and immunity/susceptibility tests can be
performed either radiated or conducted:
improper use of equipment, not from screened room, not from the mains, and
so on.
• Two EMC engineers do not agree with each other unless they are drunk.
• If an EMC engineer wants to be sure, he or she never measures twice.
• Nobody grows up wanting to be an EMC engineer; it usually just hap-
pens.
• Very often, what you predict falsifies you.
• The number of your mistakes is proportional to your experience.
• In general, SE of the tower case is better than AT case (this may change
with brand, with design, and so on).
• EMC problems of a computer with a 200-MHz processor may be im-
portant up to 10th harmonics (i.e., up to 2–3 GHz).
• Extra screening of both cases seems to modify SE in certain frequencies
but also degrade at some other frequency regions.
• SE versus frequency behaviors are too complicated.
Introduction 29
Two important conclusions that can be derived from this example are:
ered a black magic by those who do not understand EM. In reality, one can
solve some of the most complex aspects of EMC by understanding the funda-
mental or basic aspects of Maxwell’s equations.
Critical parameters/quantities of an EMC problem are frequency, signal
amplitude, time behavior, geometrical discontinuities, and boundary/edge con-
ditions. Specifically these interrelations are vital:
• Time-frequency;
• Wavelength (frequency)-geometry (size);
• Frequency-impedance.
Because:
• Crisis approach: Disregard EMC until the design is completed. Use add-
ons to fix the problems as they arise during testing and field experience.
This is also called the band-aid approach and tends to be expensive.
• Systems approach: Consider EMC throughout the design. Anticipate
all possible problems from the start. Do thorough EMC test of final
prototypes. This is more cost effective because EMC is designed into the
product and not added on when problems later arise.
References
[1] Sevgi, L., Elektomanyetik Uyumluluk—Elektromanyetik Kirlilik, Elektrik Mühendisleri
Odası (EMO) Yayınları, Istanbul, December 2000.
[2] Sevgi, L., Tasarımdan Piyasaya Elektromanyetik Uyumluluk, OKAN Üniversitesi Yayınları,
Istanbul, March 2015.
[3] Sevgi, L., Tasarımdan Üretime Elektromanyetik Uyumluluk, Eksen Yayıncılık, Istanbul,
June 2006.
[4] Sevgi, L., Ruediger Vahldieck (Guest Editor), Special issue on “From Design to Market of
Electronic Systems: EMC Engineering,” TJEECS, Turkish Journal of Electrical Engineering
and Computer Sciences, Vol. 17, No. 3, Nov. 2009, http://mistug.tubitak.gov.tr.
[5] Sevgi, L., “Modeling and Simulation Strategies in EMC Engineering: Problems, Chal-
lenges and Tools,” EMC Zurich 2009, 20th International Zurich Symposium on Electro-
magnetic Compatibility, Jan 12–16, 2009, Switzerland.
32 A Practical Guide to EMC Engineering
[6] Sevgi, L., “EMC Testing Centers and Labs: Problems Related to CE Marking,” EMC Zu-
rich 2009, 20th International Zurich Symposium on Electromagnetic Compatibility, Jan
12–16, 2009, Switzerland.
[7] WHO, World Health Organ�ization, http://www.who.int.
[8] ICNIRP, http:// http://www.icnirp.org.
[9] Sevgi, L., “Biostatistics and Epidemiology: Hypothetical Tests on Cell Phone Users,” IEEE
Antennas and Propagation Magazine, Vol. 52, No. 1, Feb. 2010, pp. 267–273.
[10] American Cancer Society, www.cancer.org.
[11] Ioannidis, John P. A., “Why Most Published Research Findings Are False,” PLoS Medicine,
Vol. 2, 8, Aug 2005, pp. 696–701.
[12] Sevgi, L.. “Biostatistics and Epidemiology: Hypothetical Tests on Cell Phone Users,” IEEE
Antennas and Propagation Magazine, Vol. 52, 1, Feb. 2010, pp. 267–273.
[13] CENELEC, http://www.cenelec.eu.
[14] Siemens, www.siemens.com
[15] IEC, http://www.iec.ch.
[16] CISPR, http://www.iec.ch/emc/iec_emc/iec_emc_players_cispr.htm.
[17] Military Standards, www.dsp.dla.mil.
Bibliography
EU Directive, http://ec.europa.eu/enterprise/sectors/electrical/emc.
2
Accreditation
Accreditation is discussed in this chapter. Worldwide accreditation institutes
are reviewed. Procedures and standards for accreditation, certification, and cali-
bration are summarized through EMC tests and measurements. Error, uncer-
tainty and statistical evaluations, and proficiency tests with typical results are
presented.
2.1 Introduction
Performance enhancing drugs (PEDs) have long been used in professional
sports and the Olympic Games. The first documented misuse was in the 1904
marathon by the winner Thomas Hicks. The International Olympic Commit-
tee (IOC) banned PEDs in 1967 and introduced the first drug use controls at
the 1968 Winter Olympics. These controls include systematic testing of ath-
letes’ urine and blood. The World Anti-Doping Agency (WADA) became the
authoritative body on the use of PEDs in 1999. The PEDs have become more
and more sophisticated, but WADA has also innovated new ways to detect these
drugs. The first positive test was at the 1968 Summer Olympics, and the Swed-
ish pentathlete Hans-Gunnar Liljenwall lost his bronze medal for alcohol use.
A well-known case was the Canadian sprinter Ben Johnson, who won the men’s
100 meters at the 1988 Seoul Olympics but lost his gold medal when he tested
positive in doping tests.
Diana Taurasi, a star of the Women’s National Basketball Association,
tested positive and was suspended for a banned substance in December 2010
while playing for the Turkish basketball team Fenerbahce, which could have
jeopardized not only her competition at the 2012 Olympic Games in Lon-
don but also her entire career. Tests done by the Doping Control Center at
33
34 A Practical Guide to EMC Engineering
How can we be sure on the results of the tests? How can we be sure on the preci-
sion of the measurement and accuracy of the result? How can we be sure that the
test/measurement devices/instruments are capable (calibrated)? How can we be
sure that the test/measurement operator is competent? How can we be sure that
everything is certified? The answer is accreditation.
2.2 Accreditation
Annual world trade (global export) approaches dozens of trillion US dollars.
Consumers buy products imported/exported all around from China to Africa,
South America to Australia. The safety and quality of the trade is controlled by
international standards and compliance programs that facilitate the movement
of goods and services across national borders. The ISO is the world’s largest
standard-setting body. Many of ISO’s standards include guidance on how to
demonstrate that a product, person, service, or system meets the requirements
contained within a standard. ISO/IEC 17011 (Conformity Assessment-Gener-
al Requirements for Accreditation Bodies Accrediting Conformity Assessment
Bodies) defines accreditation as “a third-party attestation related to a confor-
mity assessment body conveying formal demonstration of its competence to
carry out specific conformity assessment tasks.”
The philosophy of the accreditation is given in Figure 2.2. The aim is
to remove all the barriers of international trade without degrading the safety
and quality. It is an effective marketing tool that brings public confidence. It
gives confidence through ensuring consistently high standards in the quality of
products/services purchased, and minimizes failure or recalls. This is achieved
through repeatable tests/measurements carried out by a reliable independent
The activities and standards covered by most of the MLAs are listed in
Table 2.1. As listed there, members of the regional accreditation bodies provide
accreditation of labs (testing according to ISO/IEC 17025, medical examina-
tions according to ISO 15189, calibration according to ISO/IEC17025, in-
spection bodies according to ISO/IEC 17020), certification bodies (product
certification according to EN45011-ISO/IEC 17065, certification of persons
according to ISO/IEC 17024), management systems certification according to
ISO/IEC 17021.
The IAF oversees, evaluates, and formally recognizes bodies that accredit
organizations that certify management systems, products, or persons. The ILAC
oversees, assesses, and formally recognizes accreditation bodies that participate
in the accreditation of calibration and testing laboratories. It has regional co-
operations with EA, APLAC, and IAAC. The merging of IAF and ILAC in the
future is a common expectation.
Like all the other regional accreditation bodies, the EA is a nonprofit
association set up in November 1997 and registered in The Netherlands in
June 2000. EA has the overall strategic objective to safeguard the credibility of
accreditation services delivered by its member countries within the European
market.
The regional accreditation bodies’ missions include defining, harmoniz-
ing, and building consistency in accreditation in their region. They ensure com-
mon interpretation and application of the standards used by their members.
They also ensure transparency of the processes and procedures. They maintain
a multilateral agreement based on mutual recognition between accreditation
activities and reciprocal acceptance of accreditation services. They manage a
peer evaluation system consistent with international practices.
A typical structure of the EA is pictured in Figure 2.3 (structures of the
IAAC and APLAC are also the same). Some of the national institutes within
EA are the Danish (DANAK), Finnish (FINAS), French (COFRAC), German
(DAR-DKD), English (UKAS), Czech (CAI), Slovak (SNAS), Irish (INAB),
Table 2.1
MLA Accreditation Coverage
Accreditation of Scope/Activity: Standard
Laboratories: Testing and Medical examinations ISO/IEC 17025 ISO 15189
Laboratories: calibration ISO/IEC 17025
Certification bodies: product certification EN 45011 and ISO/IEC 17065
Certification bodies: certification of persons ISO/IEC 17024
Certification bodies: management systems certification ISO/IEC 17021
Inspection ISO/IEC 17020
Verification EN/ISO 14065
38 A Practical Guide to EMC Engineering
Spanish (ENAC), and Portuguese (IPAC) accreditation bodies. The Turkish ac-
creditation body (TÜRKAK) is also a full member of EA, IAF, and ILAC.
All accreditation bodies aim to provide a high level of service for regional
economies, industries, governments, and citizens. An accredited organization’s
competence, impartiality, and sustainable performance have to be demonstrat-
ed. An accreditation organization can ensure that everyone from specifiers, pur-
chasers, and suppliers to consumers have confidence in the safety and quality of
certified goods and services throughout the supply chain.
There are numerous mutual recognitions of conformity between EU and
other countries such as the United States, Japan, Canada, Australia, New Zea-
land, and Israel. Because of this, products from these countries may not have af-
fixed CE mark. Japan has its own marking known as the Technical Conformity
Mark. Switzerland and Turkey (which are not members of the EU) also require
products to bear CE marking as an affirmation of conformity.
2.4 TÜRKAK
Being related with Ministry of Industry and Trade and subject to private law
provisions, TÜRKAK was established as a legal entity having administrative
and financial autonomy with its headquarters in Ankara. (TÜRKAK Law No.
4457 was adopted on 27/10/1999 and published in the 4/11/1999/No.23866
issue of the Official Gazette.) The objective of TÜRKAK is to facilitate the free
movement of goods in the Turkish market by assuring the reliability of reports
and certifications granted by Conformity Assessment Boards (CABs) in the
Accreditation 39
Figure 2.5 A typical emission measurement setup: (1) eliminate external interference; (2)
minimize turntable contributions; (3) mutually isolate EUT and mains; (4) mutually isolate re-
ceive antenna and mains.
that there is no mutual disturbance between the mains and EUT and
that they are well isolated.
• Emissions captured by the receive antenna are recorded by an EMI re-
ceiver. The antenna, cables, and the EMI receiver must all be calibrated
[12, 13].
Only accredited labs can maintain all these requirements, and tests and
measurements can be traceable. National accreditation bodies do the accredi-
tation. Lab accreditation covers personnel management system accreditation,
the calibration of the test/measurement environments and instruments, and
accreditation of the test/measurement methods. Note that accreditation is a
dynamic process and the continuity is maintained through proficiency tests
summarized in Section 2.6.
Emissions are typically measured from a few hundred kHz up to several
GHz frequencies. Roughly speaking, CE along cables and wiring are performed
in the frequency range of 150 kHz–30 MHz and results are presented in dBµV.
Above 30 MHz CE is negligible. On the other hand, RE are performed in the
frequency range of 30 MHz–1 GHz (up to 6 GHz) and results are presented in
terms of EM field strength dBµV/m. Below 30 MHz RE is negligible.
Like all the others, EMC tests and measurements are challenging and
time consuming. For example, EN55016-2-3:2010 (CISPR 16-2-3:2010) cov-
ers RE measurements for household appliances in OATS and anechoic chambers
in the range of 30 MHz to 1 GHz. Measurements are performed with a resolu-
tion bandwidth (RBW) of 120 kHz with a scan step much less than this value.
Usually, a scan step of 50 MHz is chosen. This results in 970 MHz / 50 kHz =
42 A Practical Guide to EMC Engineering
19,400 steps. Note that peak (PK), quasi-peak (QP), and average (AV) detec-
tors are used in EMC measurements. PK measurements are a few order magni-
tudes faster than QP and AV measurements, but the limits are given in terms of
QP and AV values. EMC measurements are performed in two steps in order to
speed up the work: prescan and final. Prescan is done with the PK detector; QP
detector is used in final measurements. The typical measurement durations of
PK and QP detectors are 20s/MHz and 1 ms/MHz, which mean PK measure-
ments (one scan in 30 MHz–1 GHz band lasts 324 min in QP mode, but less
than a second with PK detector).
Note that CISPR 22 (2009) for ITE (also EN55022) necessitates condi-
tional RE measurements in the range of 30 MHz–6 GHz. If the internal EUT
sources are less than 108 MHz, then RE measurements are mandatory up to
1 GHz. If it is between 108 MHz and 500 MHz, then measurements shall be
performed up to 2 GHz. For EUTs with internal sources between 500 MHz
and 1 GHz, RE measurements are requested up to 5 GHz. Finally, if the EUT
has sources above 1 GHz, RE measurements shall be conducted up to 6 GHz or
5fmax, whichever is less (e.g., if EUT has a source at 1.1 GHz, then RE measure-
ments will be done in the range of 30 MHz–5.5 GHz). For RE measurements
above 1 GHz, the RWB is 1 MHz; measurements are requested with both PK
and AV detectors.
Both prescan and final measurements can be done in a 10m anechoic
chamber. First, the maximum emission direction of the EUT and its polariza-
tion are found with at least 4 prescans. These fast prescans yield critical frequen-
cy regions where recorded emissions are close to or above the limits given in
the standards. Then final measurements in these critical regions are done with
QP detector. the 10m anechoic chambers are expensive; therefore, a significant
number of EMC labs have 3m anechoic chambers. In this case, prescan is done
in the chamber, but final measurements are performed in an OATS. Note that
the reliability of an OATS is higher than an anechoic chamber.
2.5.2 Calibration
The CE marking of a product is possible if tests and measurements reveal that
the product meets all the conditions given in related standards. Tests and/or
measurements done in an accredited lab assuming the following:
the test with no data loss (criterion B), and the EUT either stops functioning
or its performance degrades and does recover after the test without intervention
(criterion C).
Once all the tests and measurements are over, an EMC compliance report
is issued by the accredited lab. This report contains the product info, lab ad-
dress, a copy of the lab’s accreditation certificate, info about lab facilities, mea-
surement instruments and their calibration certificates, all the standards used
in the tests and measurements, test/measurement configurations, tables and/or
figures of the test/measurement data and pass/fail info, and signatures of the
technical personnel.
Table 2.2
QP Values Recorded in Different Labs (CE Measurements of a
Hair Dryer Operating at Hot-Air Blowing Mode)
Lab|Freq 160 kHz 500 kHz 1 MHz 3 MHz 5 MHz 10 MHz 20 MHz 25 MHz
101 69.6 61.5 55.1 45.3 54.2 50.4 38.5 30.3
102 70.8 64.5 55.6 45.6 54.8 52.1 41 31.6
103 70.56 58.95 48.77 40.89 48.06 50.04 40.21 30.6
104 70.8 61.5 56.2 45.9 56.1 52.9 42.7 34
106 70.5 61 55 46 55 53 41 34
107 71.2 60.1 52 46 54 50 39 30
108 70.6 60.6 55 45.3 54.1 51 38 28
109 70.1 56 55.1 46.5 51 50 36 29
110 67 58.8 51.4 38.9 42.2 53.8 37.8 30.1
111 72.2 63.7 57.7 47.9 57.3 54.4 41.4 33.2
113 72.5 61.5 55.3 45.5 54.5 51.7 38.5 32.5
46 A Practical Guide to EMC Engineering
Figure 2.7 Typical results of a proficiency test (CE of a hair dryer operating at hot-air blowing
mode, QP values at 1 MHz).
Table 2.3
Suspicious Values of the CE Measurements Listed in Table 2.2
Conducted Em. –QP
Lab No 0.16 MHz 0.5 MHz 1 MHz 3 MHz 5 MHz
101 — — — — —
102 — — — — —
103 — — X — —
104 — — — — —
106 — — — — —
107 — — — — —
108 — — — — —
109 — X — — —
110 X — — X X
111 — — — — —
112 — — — — —
113 — — — — —
The proficiency tests map current status of accredited labs and measure-
ment centers. It increases the confidence level of the personnel of the participat-
ing labs in their abilities and knowledge of their capabilities. They can detect
Accreditation 47
any difficulties they may have with analyses, identify training needs, and have a
mechanism where these needs can be met.
The ISO/IEC 17025 standard used for laboratory accreditation necessi-
tates a laboratory to produce an estimation of the uncertainty of the conducted
measurements using accepted methods of analysis, through the production and
application of suitable uncertainty of measurement procedures [15]. Uncertain-
ty is a range that is likely to contain the true value of a quantity being measured
or calculated.
In general, no measurement or test is perfect, and the imperfections give
rise to error of measurement in the result. Consequently, the result of a mea-
surement is only an approximation to the value of the measurand and is only
complete when accompanied by a statement of the uncertainty of that approxi-
mation. It is related with the accuracy as well as the precision [16].
Accuracy is the closeness of the measured or calculated quantity to its ex-
act value and is expressed in terms of error. Error is the difference between a
measured or calculated value of a quantity and its exact value. The error may
be systematic or random. Systematic error is error that plagues experiments or
calculations caused by negative factors (e.g., a dc voltage component that is
unintentionally present, a drift because of the failure of a component, and so
on). Systematic errors can be complex but can be removed once understood or
discovered via careful controls and calibration. Random error is an error that is
always present but varies unpredictably in size and direction. Random errors are
related to the scatter in the data obtained under fixed conditions that determine
the repeatability (precision) of the measurement. Fortunately, random errors
usually follow well-behaved statistical rules. Their effects can be reduced by
repeating the measurement as often as possible.
Note that in addition to systematic and random errors, measurement er-
ror also includes drift error due to system performance changing after the cali-
bration, mostly because of climatic (e.g., temperature) variations.
Precision is a measure of the closeness of the value obtained via multiple
measurements to the true value. It is the total amount of random error present.
A very precise measurement means a small random error. Precision is given as
the percentage of the ratio of the value region to the true value of the quantity
being measured. The value region is the difference between the maximum and
minimum values in multiple measurements. As illustrated in Figure 2.8, high
accuracy means high precision, but high precision does not necessarily mean
high accuracy.
Confidence level [17] is the probability that the true value of the measure-
ment or calculation falls within a given range of uncertainty caused by the
inherent random nature of the measurement or calculation. Confidence levels
can be defined through a good understanding of the nature (probabilistic dis-
tributions) of the errors.
48 A Practical Guide to EMC Engineering
Note that P(–∞ < z < ∞). The probability of getting z between z1,2 =
±1.96 is 0.95. In other words, we’re 95% confident that z would be in the range
of –1.96 ≤ z ≤ 1.96 (i.e., the measured value is in the range of –1.96∆x ≤ x ≤
1.96∆x or approximately –2∆x ≤ x ≤ 2∆x).
The success of the world trade is directly related to the elimination of
technical/bureaucratic barriers. This can be achieved only with worldwide ac-
creditation and multilateral agreements. Accreditation is essential for the cor-
rect operation of a transparent and quality-oriented market.
The goal is to find the most effective way of achieving a single system that
will allow companies with an accredited conformity assessment certificate in
Accreditation 49
one part of the world to have that certificate recognized elsewhere. International
institutes and standards aim to ensure that accreditation bodies only accredit
bodies that are competent to do the work they undertake and reduce risks to
business and customers by ensuring that an accredited certificate may be relied
upon anywhere in the world.
References
[1] EU rapid alert system RAPEX, http://ec.europa.eu/consumers-/safety/rapex/alerts.
[2] IAF, International Accreditation Forum, http://www.iaf.nu.
[3] ILAC, International Laboratory Accreditation Cooperation, http://www.ilac.org.
[4] EU, The European co-operation for Accreditation, http://www.european-accreditation.
org.
[5] IAAC, Inter-American Accreditation Cooperation, http://www.iaac.org.mx.
[6] APLAC, Asia Pacific Laboratory Accreditation Cooperation, https://aplac.org.
[7] http://ec.europa.eu/enterprise/sectors/electrical/emc.
[8] Sevgi, L., and R. Vahldieck, special issue on “From Design to Market of Electronic Sys-
tems: EMC Engineering,” Turkish Journal of Electrical Engineering and Computer Sciences,
Vol. 17, No. 3, Nov. 2009, http://mistug.tubitak.gov.tr.
[9] Sevgi, L., “EMC Testing Centers and Labs: Problems Related to CE Marking,” EMC Zur-
ich 2009, 20th International Zurich Symposium on EMC, Jan. 12–16, 2009, Switzerland.
[10] Sevgi, L., “From Design to Market EMC,” IBCAST 2011, 8th International Bhurban
Conference on Applied Sciences and Technology, Jan. 13–16, 2011, Islamabad, Pakistan.
[11] Eser, S., and L. Sevgi, “Open Area Test Site (OATS) Calibration,” IEEE Antennas and
Propagation Magazine, Vol. 52, No. 3, Jun. 2010, pp. 204–212.
[12] Sevgi, L., S. Çakır, and G. Çakır, “Antenna Calibration for EMC Tests and Measurements,”
IEEE Antennas and Propagation Magazine, Vol. 50, No. 3, Jun. 2008, pp. 215–224.
[13] Çakır, S., R. Hamid, and L. Sevgi, “Loop Antenna Calibration,” IEEE Antennas and
Propagation Magazine, Vol. 53, No. 4, Oct. 2011, pp. 243–254.
[14] ISO/IEC Guide 43-1:1997 Proficiency Testing by Interlab Comparisons—Part 1:
Development and Operation of Proficiency Testing Schemes (revised by 17043:2010).
[15] The Expression of Uncertainty and Confidence in Measurement, UKAS Publication M3003,
First Edition, Dec. 1997.
[16] Sevgi, L., “Innumeracy: The Meaning of the Numbers We Use,” IEEE Antennas and
Propagation Magazine, Vol. 49, No. 2, Apr. 2007, pp. 195–190.
[17] Sevgi, L., “Hypothesis Testing and Decision Making: Constant-False-Alarm Rate,” IEEE
Antennas and Propagation Magazine, Vol. 51, No. 3, Jun. 2009, pp. 218–224.
[18] TSE Proficiency Test, http://global.tse.org.tr.
50 A Practical Guide to EMC Engineering
Bibliography
http://ec.europa.eu/enterprise/policies/single-market-goods/cemarking.
http://ec.europa.eu/enterprise/sectors/electrical/lvd.
http://ec.europa.eu/enterprise/sectors/toys.
3
Electromagnetic Model
EMC engineering necessitates a different look at two classical lectures, EM fields
and waves, taught to undergraduates. EM theory is briefly revisited through
they eyes of an EMC engineer. For further reading, a few suggested classical
books on EM are [1–5].
51
52 A Practical Guide to EMC Engineering
• EM waves interact with objects. Surface currents are generated during
this interaction, and objects radiate like secondary EM sources.
• Materials with nonzero total magnetic moments of their atoms (like a
magnet) generate magnetic fields around.
Table 3.1
Maxwell Equations and EM Model
Differential Form Integral Form
∇ ⋅D = ρ
∫ Dds = ∫ pdv
S V
∇⋅B = 0
∂B
∫ Bds = 0
S
∇×E = − +J
∂t ∂
D = εE , B = µH , J = σE
∫ Edl ∂t ∫S
=− Bds
C
∫ Hdl = ∫ Dds + ∫ Jds
C S S
Electromagnetic Model 53
Table 3.2
EM Quantities and Their Units
Symbol Quantity Unit
E Electric field Volt/m
H Magnetic field Amper/m
J Surface current density Amper/m2
ρ Volume charge density Coulomb/m3
ε Medium permittivity Farad/m
σ Medium conductivity Siemens/m
µ Medium permeability Henry/m
Cylindrical waves far from a line source may be assumed planar and are
called plane waves. The plane wave is the simplest wave object that can easily
be modeled mathematically. Both reflection and refraction occur when an EM
wave hits a penetrable surface. Reflection and refraction depend on the electri-
cal parameters of the surface as well as the polarization of the source.
In free space, EM waves propagate with the speed of light (c=3108 m/s).
The wavelength (λ) of the wave in free space is obtained from λ=c/f, where f is
the frequency of the wave. For example, at 10 GHz the wavelength is λ=3 cm.
At 1 GHz it is λ=30 cm. At 100 MHz it is λ=3m. A practical equation for the
wavelength is λ=300/fMHz [m].
∂2 ∂2 ∂2 2
∂x 2 ∂y 2 ∂z 2 + k ψ = 0
+ + (3.1)
The source term in the wave equation is given on the right-hand side.
Source-free means either the sources are very far, or that they are first applied,
the fields are excited, and then are removed. Source-free wave equation (3.1) is
used to investigate the environment. The solutions are called eigenfunctions, and
characteristic values that satisfy eigenfunctions are called eigenvalues. Therefore,
together with boundary/radiation conditions, (3.1) is called an eigenvalue
problem [6].
If a point source is located at the origin of Cartesian coordinate system,
(3.1) becomes
∂2 ∂2 ∂2 2
∂x 2 ∂y 2 ∂z 2 + k g = − δ ( x ) δ ( y ) δ ( z )
+ + (3.2)
be known. The Green’s function is the solution of (3.2) with the point source.
The Green’s function in EM theory corresponds to the impulse response in the
circuit theory.
Due to angular symmetry of the point source at the origin, it would be
better to deal with (3.2) in spherical coordinates. The full symmetry in both
horizontal and vertical angles the solution (i.e., the Green’s function would be
only radial distance dependent):
e − jkr e jkr
g= or g = (3.3)
4 πr 4 πr
These two, respectively, represent outgoing waves from the origin and
incoming waves toward the origin. The Green’s function (i.e., the solution of
the wave equation under the point source) decays with distance as 1/r and the
phase changes as (–jkr). Here, k = ω/c is the wave number).
If there is symmetry along one direction (e.g., along z), the wave equation
can be reduced to 2D (by assuming, ∂/∂z ≡ 0). In this case, the point source
reduces to a line source and can be given as δ(x)δ(y). The Green’s functions in
2D will then be:
e − jk ρ e jk ρ
g= or g = (3.4)
2π ρ 2π ρ
Here, ρ is the distance from the origin on the xy plane. As observed, waves
decay with one over square root of the distance (1 ρ) and the phase change
will be in the same form (exp(–jkρ)).
The Green’s function problem and the eigenvalue problem are not in-
dependent. One can be obtained from the other. This means EM behavior of
the structure at hand can be used to obtain the source effects, or vice versa.
The characteristic problem which connects source-driven and source-free solu-
tions (i.e., the Green’s function and the eigenfunction) is called Sturm-Liouville
problem [7, 8].
• Linear polarization: If the tip of the electric field vector moves along a
straight line, then the polarization is said to be linear. The two cases of
Electromagnetic Model 57
are three characteristic problem groups in EM, and every area in Figure 3.4 falls
in one of these three groups.
Table 3.3
EM Materials and Their Electrical Parameters (σCu=5.8 107
S/m, µCu= 1.0 H/m )
Relative
permittivity Relative Relative
(wrt free- Permeability Conductivity
Material space) (wrt copper) (wrt copper)
Air 1.0005 0.0000 0.0000
Glass 5.0000 0.0000 0.0000
Nylon 3.5000 0.0000 0.0000
Teflon 2.1000 0.0000 0.0000
Bakelite 4.9000 0.0000 0.0000
Porcelain 6.0000 0.0000 0.0000
Silicon 12.0000 0.0000 0.0000
Table 3.4
EM Materials and Their Electrical Parameters
(σCu=5.8 107 S/m, µCu= 1.0 H/m)
Relative Relative
Conductivity Permeability
Material (wrt copper) (wrt copper)
Copper 1.0000 1.0000
Brass 0.260 1.0000
Gold 0.700 1.0000
Silver 1.0800 1.0000
Nickel 0.2200 1.0000
Aluminum 0.6100 1.0000
Zinc 0.2900 1.0000
Steel 0.1000 500.000
Stainless Steel 0.0200 1000.0000
Iron 0.1700 1000.0000
Permaloy 0.0300 80,000.0000
Mu-metal 0.0300 30,000.0000
2
δ= [m] (3.5)
ωµσ
Electromagnetic Model 61
Figure 3.6 Basic radiators and their radiation patterns: (a) line source; (b) Poynting vector;
(c) vertical dipole; (d) horizontal dipole.
References
[1] Collin, R. E., Field Theory of Guided Waves, New York: McGraw-Hill, 1960.
[2] Jordan, E. C., and K. G. Balmain, Electromagnetic Waves and Radiating Systems, Engle-
wood Cliffs, NJ: Prentice-Hall, 1968.
[3] Dudley, D. G., Mathematical Foundations for Electromagnetic Theory, New York: IEEE
Press (Series on Electromagnetic Waves), 1994.
[4] Felsen, L. B., and N. Marcuvitz, Radiation and Scattering of Waves, Englewood Cliffs, NJ:
Prentice-Hall, 1973, Reissued IEEE Press, 1994.
[5] Sevgi, L., Electromagnetic Modeling and Simulation, New York: IEEE Press–John Wiley &
Sons, 2014.
[6] Sevgi, L., “Guided Waves and Transverse Fields: Transverse to What?” IEEE Antennas and
Propagation Magazine, Vol. 50, No. 6, Dec. 2008, pp. 221–225.
[7] Sevgi, L., “Guided Waves and Eigenvalue Extraction from Propagation Characteristics,”
IEEE Antennas and Propagation Magazine, Vol. 50, No. 2, Apr. 2008, pp. 222–234.
[8] Sevgi, L., “Sturm-Liouville Equation: The Bridge Between Eigenvalue and Green’s Func-
tion Problems,” ELEKTRIK, Turkish J. of Electrical Engineering and Computer Sciences,
Vol. 14, No. 2, 2006, pp. 293–311.
4
Circuit Model
One of the first technical lectures in electric, electronics, control, computer, and
communication engineering is circuit theory. Three important circuit elements
are resistor, inductor, and capacitor. Students learn definitions, their serial and
parallel connections, and functions in simple few-loop circuits, and they calcu-
late element currents (voltages) for given voltage (current) sources. For example,
resistance is defined as the ratio of voltage between its ends and current flowing
through and is the same for both DC and AC circuits. The impedance concept
is used to define inductor and capacitor in an AC circuit, which act as short-
circuited (SC) and open-circuited (OC) in DC circuits, respectively. Everything
works perfectly in circuit theory when definitions are used in ideal conditions.
However, EMC engineers have to go further and need to understand and model
these three elements and their circuits under realistic conditions. How does a
resistor/inductor/capacitor act at high frequencies? Is a resistor only a resistor
in practical circuits? How do EMC engineers define these elements in practice?
What is the difference between “lumped” and “scattered” element definitions?
65
66 A Practical Guide to EMC Engineering
Table 4.1
AWG Standards (1 mil = 10-3 in = 2.54
10-3 cm)
Standard Diameter Diameter
Wire [mil] [mm]
AWG 50 1.0 0.0254
AWG 44 2.0 0.0508
AWG 38 4.0 0.1016
AWG 32 8.0 0.2032
AWG 26 16.0 0.4064
AWG 20 32.0 0.8128
AWG 14 64.0 1.6256
Circuit Model 67
in gray areas. The unit resistances of the wire for DC and AC circuits are also
given in the figure. As shown there, unit resistance is inversely proportional
with conductivity and the cross section. Since the area A1 is much less than A,
RAC will be much greater than RDC.
Figure 4.3 shows RAC/RDC versus frequency of copper and iron wires for
different wire cross sections. As observed, this ratio is nearly 20 for a circular
cross-section copper wire with 20-mm radius at 1 kHz, but it increases to 700
at 1 MHz.
EMC engineers must pay attention to skin depth even in very simple
circuits in order to prevent circuit failures.
4l
L = 0.002 l 2.3 Log 10 − 0.75 (4.1)
d
bl
L=µ (4.2)
d
dl (4.3)
C =ε
b
Equations (4.2) and (4.3) yield L = µ and C = ε when the length, diam-
eter, and distance are assumed 1m. This helps to bring different interpretations
for EM parameters permittivity and permeability. The inductive and capacitive
effects are related with magnetic and electric energy storage capacities, respec-
tively. It should be noted that models given in Figure 4.1 belong to ideal case;
in practice element models are quite different.
Here are few circuit tips in terms of EMC engineering:
and load, and VL2 is the load voltage when the source is directly connected to
the load. IL in dB is then given as ILdB=20log10|VL2/VL1|. IL is negative with
the maximum value of 0 dB, which corresponds to VL2 = VL1. It should be
noted that the definition of IL may also be given as ILdB=20log10|VL1/VL2|. In
this case, IL would be positive with the minimum value of 0 dB. RL is defined
at the input port in terms of incident (V +) and reflected (V –) voltage waves as
RLdB=20log10| V –/V +|. The two-port circuit’s performance is then given as IL
versus frequency or RL versus frequency curves. Note also that voltage transfer
function (VTFdB=20log10|V2/V1|) is also used to show the performance of a
two-port circuit. Actually, the VTF curve is simply a 6-dB shifted (below) IL
curve.
Calculation of both IL and RL necessitates input impedance (Zin) and
input voltage (Vin). As shown in Figure 4.10, the input impedance (Zin = ZC in
this case) may be obtained, from load to source, by serial/parallel reduction of
the circuit impedances. Then, input voltage can directly be calculated from a
single loop circuit with the equivalent input impedance. Once the source cur-
rent is obtained, the currents through all branches, from source to load, can eas-
ily be obtained. This will yield all required voltages for IL and RL calculations.
In Figure 4.10, a π-circuit with three impedances is used as a two-port
example. The impedance Za is the parallel combination of RL and Z3; Zb is the
serial combination of Za and Z2; the input impedance Zc is the parallel com-
bination of Zb and Z1. Once they are calculated, currents flowing through all
branches will easily be obtained as shown in the figure. Table 4.2 lists a short
MATLAB code for the calculation of IL versus frequency of the π-circuit shown
in Figure 4.10. A typical result obtained with this code is shown in Figure 4.11.
In Figure 4.11, both IL and VTF versus frequency variations are plotted.
Z1 is a parallel resonance circuit with L1 = 0.949 nH and C1=11.867 pF. Z3 is
also a parallel resonance circuit with L3 = 0.964 nH and C3 =11.676 pF. Z2 is a
Table 4.2
Short MATLAB Code for the IL Calculations of a π-Circuit
%--------------------------------------------------------------------------
% Program: IL versus Frequency calculations for a PI circuit
%--------------------------------------------------------------------------
clc; clear all; rs=50; rl=50; vs=100; nn=500;
fmin = 1; fmax = 3e3; df=(fmax-fmin)/nn; % Frequency range
% Impedance values for L/C elements)
l1 = 0.949e-9; c1 = 11.867e-12; % Z1 parallel resonance circuit
l2 = 20.372e-9; c2 = 0.552e-12; % Z2 serial resonance circuit
l3 = 0.964e-9; c3 = 11.676e-12; % Z3 parallel resonance circuit
for k=1:nn
f(k) = (fmin+(k-1)*df); omega = 2*pi*f(k)*1e6;
zl1 = 1i*omega*l1; zc1 =-1i/(omega*c1); Z1=(zl1*zc1)/(zl1+zc1);
zl2 = 1i*omega*l2; zc2 =-1i/(omega*c2); Z2=zl2+zc2;
zl3 = 1i*omega*l3; zc3 =-1i/(omega*c3); Z3=(zl3*zc3)/(zl3+zc3);
za=(rl*Z3)/(rl+Z3); zb=Z2+za; zc=(zb*Z1)/(zb+Z1);
I1=vs/(rs+zc); I1B=I1/(1+zb/Z1); IY=I1B/(1+rl/Z3);
vl=IY*rl; vtf(k)=vl/vs; IL(k)=2.*vl/vs;
end
plot(f,20*log10(abs(vtf)),’b’,f,20*log10(abs(IL)),’r--’); ylim([-50 5]);
legend(‘GTF’, ‘AGK’); title(‘2-port behavior’)
ylabel(‘IL [dB]’); xlabel(‘Frequency [MHz]’) % Program End
serial resonance circuit with L2 = 20.372 nH and C2 = 0.552 pF. These values
belong to a three-element band pass filter. As observed, IL is 0 dB in the pass
band (i.e., between 1.3 and 1.8 GHz).
74 A Practical Guide to EMC Engineering
1
f0 = (4.4)
2 π LC
• Lowpass and highpass filters can be realized using only an inductor and/
or a capacitor. Bandpass and bandstop filters necessitate serial and/or
parallel resonance circuits.
• The inductor is SC for DC signals but acts OC at high frequencies.
On the contrary, the capacitor is SC for DC signals and acts SC at high
frequencies.
• An inductor (capacitor), connected in serial between input and output,
acts as a low (high) pass filter.
• An inductor (capacitor), connected in parallel between input and out-
put, acts as a high (low) pass filter.
fm
Q= ; B = fU − f L (4.5)
B
is the quality factor and show the selectivity of a resonance circuit. High-
Q means better frequency selectivity.
• Shape factor: The ratio of 60-dB and 3-dB bandwidths is the shape fac-
tor of a resonance circuit. It shows the degree of the amplitude decay be-
tween pass and stop bands. In ideal filters, the shape factor is equal to 1.
R + j ωL
Z0 = [Ω] (4.6)
G + j ωC
Circuit Model 79
−
V Z − Z0
Γ = ρ= = (4.8)
V + Z + Z0
• When ZL = Z0, reflection coefficient is zero (Γ = 0). The power along the
line totally transferred to the load. This is called matched TL.
• When ZL = 0 (SC), reflection coefficient is minus one (Γ = –1). In this
case, the power along the line totally reflects back (with a phase shift of
180°).
• When ZL →∞ (OC), reflection coefficient is one (Γ = 1). In this case,
the power along the line totally reflects back (without phase shift).
V max 1+ | ρ |
s= =
V min 1− | ρ | (4.9)
82 A Practical Guide to EMC Engineering
ized with a SC TL with length l can also be realized with OC line with
length l + λ/4.
• Electrical characteristics on a TL are periodic with half wavelength (λ/2).
4.5 Grounding
The voltage (potential difference) is the work capability between any two
points. In electrical engineering, voltage is the potential difference between any
two points. For example, VAB is the difference between potentials of the points
A and B. But, the potential of a single point (e.g., VA) is meaningful if a refer-
ence (zero potential) point is specified. This is called ground.
The classical definition of a ground is “an equipotential point or plane
which serves as a reference for a circuit or system” [12]. Unfortunately, this
definition is meaningless in the presence of ground current flow. Even where
signal currents are negligible, induced ground currents due to environmental
magnetic or electric fields will cause shifts in ground potential. An alternative
definition for a ground is “a low impedance path by which current can return
to its source” [12]. This emphasizes current flow and the consequent need for
low impedance and is more appropriate when high frequencies are involved. It
is important to remember that two physically separate “ground” points are not
at the same potential unless no current is flowing between them.
According to IEEE Std 142-2007, the definition of “ground” is “a con-
ducting connection, whether intentional or accidental, between an electrical
circuit or equipment and Earth, or to some other body that serves in place of
Earth.”
Grounding is one of the primary ways of minimizing unwanted noise and
of producing a safe system. That said, a noise-free system is not necessarily a safe
system, and, conversely, a safe system is not necessarily a noise-free system. It is
the responsibility of the designer to provide both a safe and a noise-free system.
One advantage of a well-designed ground system is that it can often provide
protection against unwanted interference and emission, without any additional
per-unit cost to the product.
Grounds fall into two categories: (1) safety grounds and (2) signal grounds.
The second category probably should not be called grounds at all, but rather
returns, and it could be further subdivided into either signal or power returns.
If a ground is connected to the enclosure or chassis of the equipment, then it
84 A Practical Guide to EMC Engineering
radius r, area A, distance to ground h, these impedances are Zin = l /Aσ at low
frequencies, Zin = jωL at medium frequencies, and Zin = jZ0tan(2πl /λ) at high
frequencies (the wire and ground form a kind of TL and its characteristic im-
pedance is Z 0 = 60 µr εr ln (2h r )). This input impedance becomes infinite
(zero) at frequencies when l = λ/4 (l = λ/2).
At low frequencies, a DC ground resistance gives a clue about ground im-
pedance, and single point grounding may be applied. The length of the ground-
ing conductor should be much less than a wavelength (as a rule of thumb, l<
λ/10). On the other hand, when the frequency is high, scattered inductance and
leakage capacitance effects create grounding problems.
Don’t forget: A nonzero (even as small as 0.3 Ω–0.5 Ω), imperfect ground
impedance may be negligible under normal operations but can create several
tens of kV potential differences under a 200 kA lightning. This in turn causes
high discharge and damage equipment. Grounding is also important in protect-
ing personnel, animals, structures, and buildings from harm because of electri-
cal shock or fire.
Grounding is a major problem in PCB circuits. Here are a few tips for a
good grounding:
luxury. EMC deals with real-life problems, and real-life problems are not linear.
Moreover, most of the systems have nonlinear circuit elements like a mixer,
amplifier, attenuator, and so on. Even, for a single-frequency input, they create
undesired signals at several frequencies. The signals with integer multiples of
the input signal are called harmonics. If there are multiple frequency inputs,
nonlinear circuits also create internal products (IPs) at their outputs. Genera-
tions of harmonics and internal products are pictured in Figure 4.22. Nonlinear
circuits have transfer functions with high-order terms; therefore, for an input
signal f1, they create signals f2 = 2f1, f3 = 3f1, f4 = 4f1, ..., and so on. These har-
monics cause serious EMC problems.
Similarly, a mixer, for example, with f1 and f2 signals at the input, creates
both harmonic signals (i.e., 2f1, 3f1, 4f1, …, 2f2, 3f2, 4f2, …) and internal prod-
ucts (f1 ± f2, 2f1 ± f2, f1 ± 2f2, 3f1 ± f2, ..., mf1 ± nf2) where m and n are integers.
The order of the internal product is important and would be (m + n + k) for
the internal product (mf1 ± nf2 ± kf3). Third order internal products (TOIP) are
usually critical because they fall in the signal band.
Figure 4.23 shows harmonic and internal products of two signals A and B
with close frequencies. As observed, the two third order internal products 2A-B
and 2B-A are very close to the signals A and B with significant amplitudes. The
other harmonics and internal products can be filtered out easily, but these two
TOIPs necessitate sharp filters.
Circuit Model 89
V1 = Z 11I 1 + Z 12 I 2
(4.10)
V 2 = Z 21I 1 + Z 22 I 2
Here, for example, Z11 is the input impedance extracted under a specific
condition (i.e., when the output port is OC). It is obtained from
V1
Z 11 = (I 2 = 0) (4.11)
I1
Vi +
ai = , i = 1,2,3,, N (normalized incident voltage wave ) (4.12)
2 Z 0i
Vi −
bi = i = 1,2,3,, N (normalized reflected voltage wave ) (4.13)
2 Z 0i
symmetry; therefore, S11 = S22. Using current continuity (i.e., I1 = –I2), the
other two parameters are found as S21 = S12 = 1 – ρ.
The circuit in Figure 4.26(b) is impedance connected in parallel between
input and output. The S-parameters of this circuit are S11 = ρ, S12 = 1 + ρ, S21
= 1 +, S22 = ρ, where the voltage reflection coefficient is ρ = ((Z1//Z0) – Z0)/
((Z1//Z0) + Z0).
Question: For the π–circuit given in the figure, SdB = 20log10S, show that
S11 = S22 = –53dB and S12 = S21 = –6dB. Here, Z1 = 150Ω, Z2 = 150Ω, Z3 =
37Ω, and Z0 = 50Ω.
Table 4.3 lists a short MATLAB code for the S-parameter calculations of
this two-port PI-circuit. The output of this code is given in Figure 4.27. The
dB values of the calculation results match with the values given in the question.
Another MATLAB code for the TEE circuit given in Figure 4.26(d) is
listed in Table 4.4. A sample output of this code is given in Figure 4.28.
Some important properties of S-parameters are as follows:
Table 4.3
MATLAB Code for Calculating S-Parameters of the PI-Circuit
%-----------------------------------
clear all; close all; clc;
z0 = 50; % Characteristic impedance of the TL;
z1 = 150; z2 = 150; z3 = 37; %z1 = 100-50j;
% Calculate Y-parameters
y1=1/z1; y2=1/z2; y3=1/z3; y11=y1+y3; y12=-y3; y21=-y3; y22=y2+y3;
% Calculate S-parameters
s11=((1/z0-(y1+y3))*(1/z0+(y2+y3))+y3^2)/((1/z0+(y1+y3))*(1/z0+(y2+y3))-y3^2);
s12=(2*y3/z0)/((1/z0+(y1+y3))*(1/z0+(y1+y3))-y3^2);
s21=(2*y3/z0)/((1/z0+(y1+y3))*(1/z0+(y1+y3))-y3^2);
s22=((1/z0+(y1+y3))*(1/z0-(y2+y3))+y3^2)/((1/z0+(y1+y3))*(1/z0+(y2+y3))-y3^2);
% Display Results
fprintf(‘ \n’); fprintf(‘ \n’);
fprintf(‘-------------Z3----------- \n’); fprintf(‘ | | \n’)
fprintf(‘ Z0 Z1 Z2 Z0 \n’); fprintf(‘ | | \n’)
fprintf(‘---------------------------- \n’); fprintf(‘ \n’);
fprintf(‘ Z1 = %5.1f Ohm; Z2 = %5.1f Ohm; Z3 = %5.1f Ohm \n’,z1,z2,z3);
fprintf(‘ \n’); fprintf(‘---S-Parameters of the above PI Circuit--- \n’);
fprintf(‘---------------------------------------------------- \n’); fprintf(‘ \n’);
fprintf(‘ S11 = ‘); disp(s11); fprintf(‘ S22 = ‘); disp(s22);
fprintf(‘ S12 = ‘); disp(s12); fprintf(‘ S21 = ‘); disp(s21); % Program END
Table 4.4
Short MATLAB Code for the TEE-Circuit in Figure 4.26(d)
%-----------------------------------
clear all; close all; clc;
z0 = 50; z1 = 150; z2 = 150; z3 = 37; %z1 = 100-50j;
% Calculate S-parameters
s11=((z1+z2-z0)*(z2+z3+z0)-z2^2)/((z1+z2+z0)*(z2+z3+z0)-z2^2);
s22=((z1+z2+z0)*(z2+z3-z0)-z2^2)/((z1+z2+z0)*(z2+z3+z0)-z2^2);
s12=(2*z2*z0)/((z1+z2+z0)*(z2+z3+z0)-z2^2);
s21=(2*z2*z0)/((z1+z2+z0)*(z2+z3+z0)-z2^2);
fprintf(‘ \n’); fprintf(‘ \n’); % Display Results
fprintf(‘--------Z1---------Z3--------- \n’);
fprintf(‘ | \n’);
fprintf(‘ Z0 Z2 Z0 \n’);
fprintf(‘ | \n’);
fprintf(‘----------------------------------- \n’); fprintf(‘ \n’);
fprintf(‘ Z1 = %5.1f Ohm; Z2 = %5.1f Ohm; Z3 = %5.1f Ohm \n’,z1,z2,z3);
fprintf(‘ \n’);
fprintf(‘---S-Parameters of the above TEE Circuit--- \n’);
fprintf(‘--------------------------------------------- \n’); fprintf(‘ \n’);
fprintf(‘ S11 = ‘); disp(s11); fprintf(‘ S22 = ‘); disp(s22);
fprintf(‘ S12 = ‘); disp(s12); fprintf(‘ S21 = ‘); disp(s21); % Program END
• S-matrix is unity in lossless circuits; that is, total power entering from
the ports is equal to the total reflected power (|a1|2 + |a2|2 = |b1|2 + |b2|2).
Why S-Parameters?
Because they are relatively easy to obtain at high frequencies; voltage traveling
waves can be measured with a vector network analyzer; they don’t need SC/OC
terminations, which can cause active devices to oscillate or self-destruct; they
are related to familiar measurements (gain, loss, reflection coefficient, and so
on); S-parameters of multiple devices can be cascaded to predict system perfor-
mance; H, Y, or Z parameters can be computed from S-parameters if desired;
and S-parameter files can easily be imported and used in numerical simulation
tools.
As shown in the figure, electric field lines go from the ribbon to the
ground plane, most of which are concentrated underneath the ribbon. Some
lines partially extended into the free space above the substrate. The magnetic
field lines encircle the microstrip ribbon and extend above the substrate. The
presence of electric field lines in the air reduces the effective dielectric constant
seen by the propagating wave. If all the fields existed between the ribbon and
ground plane, the dielectric constant would be that of the substrate. Instead, it
is somewhat less than that constant, and it depends on ribbon width, substrate
dielectric constant, and height.
EM energy is trapped under the top ribbon inside the dielectric substrate.
The trapped energy is proportional with the permittivity of the substrate. In
terms of EMC engineering, this reduces undesired emission. Characteristic im-
pedance of the microstripline is determined by the ratio W/h (0.1 ≤ W/h ≤
3.0).
In Figure 4.30, characteristic impedance versus frequency for different
W/h ratios is plotted. When the thickness of the strip is negligible (compared
to its width), characteristic impedances between 10–20Ω up to150–250Ω can
be obtained by changing the width, the substrate height, or substrate relative
permittivity.
Circuit Model 97
Table 4.5
Short MATLAB Code for the Calculating Z0 of a Microstripline
% Microstrip characteristic impedance plot
clc; clear all; close all;
er=9.6 % relative permittivity of the substrate
r=linspace(.01,1,100); % W/h <= 1
eff=(er+1)/2+((er-1)/2)*((1+12./r).^(-1/2)+0.04*(1-r).^2);
Zo=(60./sqrt(eff)).*log(8./r+r/4);
plot(r,Zo); xlabel(‘W/h ratio’); ylabel(‘Zo (Ohms)’); hold on;
grid on; title(‘Microstrip Characteristic Impedance versus W/h’)
r=linspace(1,5,400); % W/h > 1
eff=(er+1)/2+((er-1)/2)*(1+12./r).^(-1/2);
Zo=(120*pi)./(sqrt(eff).*(r+1.393+0.667*log(r+1.444)));
plot(r,Zo); % Program END
98 A Practical Guide to EMC Engineering
4.10 Crosstalk
Crosstalk is any phenomenon by which a signal transmitted on one cable cre-
ates an undesired effect in another cable. It is an internal EMC problem. It
is usually caused by undesired capacitive (electric field), inductive (magnetic
field), or conductive (common impedance) coupling from one cable to other
or one cable wire to another. Near end crosstalk (NEXT) is a measure of the
ability of cabling to reject crosstalk. Interference between two pairs in a cable
is measured at the same end of the cable as the interfering transmitter. Far end
crosstalk (FEXT) is the interference between two pairs of a cable measured at
the other end of the cable with respect to the interfering transmitter.
Generally speaking, a ground plane is a useful tool to combat crosstalk.
A good grounding may cause up to a 40–60 dB reduction in the common
100 A Practical Guide to EMC Engineering
impedance. The ground plane may also reduce mutual inductance coupling.
The capacitive coupling between cables or wires is not directly affected by the
ground plane, but the lowered common impedance indirectly reduces the ca-
pacitive crosstalk amplitude.
Crosstalk may be a major EMC problem in digital circuits.
References
[1] Electromagnetic Compatability Industry Association, journal articles of Armstrong, K.,
www.compliance-club.com. Last accessed, Oct. 21, 2016.
[2] Williams, T., EMC for Product Designers, 3rd Ed., Oxford: Newnes, 2001.
[3] Hong, J. G., and M. J. Lancaster, Microstrip Filters for RF/Microwave Applications, New
York: John Wiley and Sons, 2001.
[4] Sevgi, L., Electromagnetic Modeling and Simulation, New York: IEEE–John Wiley and
Sons, 2014.
[5] Sevgi, L., Complex Electromagnetic Problems and Numerical Simulation Approaches, New
York: IEEE–John Wiley Press, June 2003.
[6] Sevgi, L,. “Electromagnetic Modeling and Simulation: Challenges in Validation, Verifica-
tion and Calibration,” IEEE Trans on EMC, Vol. 56, No. 4, Aug. 2014, pp. 750–758.
[7] Çakır, G., S. Gündüz, and L. Sevgi, Geni Bantlı Filtre Tasarımı, Eksen Yayıncılık, Istanbul,
2006.
Circuit Model 101
[8] Gündüz, S., G. Çakır, and L. Sevgi, “A Generic Microstrip Structure for the Realization of
All-Type Broadband Filters,” Microwave and Optical Technology Letters, Vol. 48, No. 12,
Dec. 2006, pp. 2390–2393.
[9] Çakır, G., and L. Sevgi, “A Design of Microstrip Electromagnetic Double-Band Band-
gap (DB-EBG),” Microwave and Optical technology Letters, Vol. 48, No. 7, Jul. 2006, pp.
1405–1408.
Bibliography
Bowick, C., Circuit Design, Boston: Newness, 1982.
Besser, L., and R. Gilmore, Practical RF Circuit Design for Modern Wireless Systems Volume I,
Norwood, MA: Artech House, 2003.
Ludwig, R., and P. Bretchko, RF Circuit Design Theory and Applications, Englewood Cliffs, NJ:
Prentice Hall, 2000.
5
Antennas and Antenna Calibration
An antenna is a transducer [1]. A transducer is a device that transforms one type
of energy to another (e.g., mechanical to heat or kinetic to electrical, and so
on). Some may say that an antenna is not a transducer since it converts electri-
cal energy to electrical energy again. But it is appropriate to call an antenna a
transducer because it transforms electrical energy to EM energy and vice versa.
An antenna is a reciprocal device and behaves the same on transmit and
on receive [2]. The transmit antenna converts the voltage between its input ter-
minal, or the current flowing through to EM fields, while the receive antenna
does the opposite. It is a device that provides a transition from a guided wave on
a TL to a free-space EM wave (or vice versa). The IEEE definition of an antenna
is given as [3] “that part of a transmitting or receiving system that is designed to
radiate or receive EM waves.”
Antennas can focus and shape the radiated power in space by enhancing
the power in some desired directions and suppressing it in others. Antenna
theory is a well-established topic, so there are wonderful classical books (see
e.g., [4–15]).
Loosely speaking, antennas may be grouped into two classes depending
on their usage: communication antennas and EMC antennas. The fundamental
reason for this grouping is that communication and EMC antennas are used
in the far and near field regions, and mostly for narrowband and broadband
applications, respectively. Nobody deals with a TV and/or a radio broadcast an-
tenna once it is installed, except during periodic maintenance and failures. But
an EMC antenna is a device that is used many times a day in screened/anechoic
rooms on OFTS during emission measurements and immunity tests as either
receive or transmit elements.
A typical RE measurement setup is pictured in Figure 5.1. These mea-
surements are performed with broadband log-periodic, spiral, and biconic
103
104 A Practical Guide to EMC Engineering
antennas. EMC engineers perform tests and measurements almost every day
and write their reports. They record unwanted radiations from the EUT and
compare results with the threshold values in a specified standard. The objective
of immunity tests is to find out whether or not the EUT yields its designated
performance under the influence of EM interference (fast transients, ESD,
surges, voltage dips, and so on) at a specified level. No matter what the tests are,
EMC engineers usually do not pay much attention to fundamental antenna pa-
rameters such as standing wave ratio, gain, radiation pattern, input impedance,
and so on. The antenna factor (AF) and the transmit antenna factor (TAF) are
the two most important antenna parameters for the EMC engineers.
The primary design parameter for an EMC antenna is bandwidth;
VSWR, gain, efficiency, input impedance, and so on, are of secondary impor-
tance. This is mainly because EMC tests and measurements are performed over
broad frequency ranges specified by the standards. For example, RE measure-
ments of most commercial electrical and electronic devices are required in the
frequency range of 30 MHz–6 GHz. On the other hand, antenna engineers
who design communication antennas are accustomed to work with those sec-
ondary parameters.
The antenna is an interesting device. It combines well-established EM
and circuit theories. Therefore, an antenna engineer should be familiar with
both its EM and circuit equivalents.
Figure 5.2 depicts one of the simple models of a transmit/receive antenna
pair. The basic circuit representation of a transmit antenna consists of a com-
plex impedance ZA connected across a voltage source with an internal imped-
ance Zs . The radiated power Pradiated at the transmit side is determined by the
antenna input impedance and the voltage across its input terminal supplied
by the transmit circuit. The mismatch at the transmit antenna input deter-
mines the antenna input power PA extracted from the supplied power Ps (under
Antennas and Antenna Calibration 105
Here Ps is the source power in dB, Gt is the transmit antenna gain in dBi,
and Ltot is the total loss that accounts for impedance mismatch as well as cable
losses, insertion losses, and so on, in dB. Note that effective radiated power
(ERP) instead is used when the antenna gain is given in dBd.
The free-space path loss (or power dissipation, LFree) of an antenna range
determines the difference in power levels between the output of the transmit
antenna and the output of an isotropic (0 dBi) antenna located at the receive
site. This free-space loss is due to the dispersive nature of a transmitting antenna
(not atmospheric losses). A transmitting antenna radiates EM waves with a
spherical wavefront; only a portion of this spherical wavefront is captured by
the receiving antenna. For a free space, at a far-field range d and a given fre-
quency f, the path loss in dB is easily determined as follows:
(5.3)
PA = EIRP − LFree + Gr
All of these are derived from the well-known Friis formula:
2 2
λ 4π d
PA = Pt Gt Gr ; L =
4 π d
Free
λ (5.4)
G r λ2
Ae = (5.5)
4π
is used as the effective aperture of the receive antenna. Equation (5.4) is widely
used in signal power calculations in communication links assuming antennas
108 A Practical Guide to EMC Engineering
λ
PA = Pt
4 π d
(
Gt 1 − Γt
2
)G (1 − Γ )
r r
2
(5.6)
Note that the gains in (5.4) are as defined in the standard IEEE defini-
tions, while (modified) gains (used in the EMC society) derived from [14]
(
Gt' = Gt 1 − Γt
2
), (
Gr' = Gr 1 − Γr
2
) (5.7)
1 ( VA )rms
2 2
1 VA
(PDelivered )max =
8 RA
or (PDelivered )max =
4 RA
(5.9)
when RA = RRec and XA = –XRec (here VA and (VA)rms are the peak and root-
mean-square values of the source and are related via (V A )rms = V A 2 for the
sinusoidal case). In many cases, the antenna impedance is not known; therefore,
VSWR is measured and used.
V Re fl Z A − Z0 1+ | Γ |
Γ= = , VSWR = s = (5.10)
Vinc Z A + Z0 1− | Γ |
Antennas and Antenna Calibration 111
Note that Γ in (5.10) is for the direction toward the antenna; it becomes
Γ = (Z0 – ZA)/(Z0 + ZA) if the direction is toward the cable (connecting to the
source). But, the percentage of the reflected power does not change with direc-
tion, since it is proportional with |Γ|2 (reflection loss in dB is defined as RL =
–20log10|Γ|). If the power delivered to the antenna by the source is Ps, then
antenna power (the sum of radiated power and the power dissipated by the
antenna itself because of the losses) is PA = Ps(1 – |Γ|2)
Ae =
( )
R Rec I A2
rms
(5.11)
Pd
where (IA)rms is the rms value of the current flowing through the circuit shown
in Figure 5.4:
(V )
(I )
A rms
=
A rms
(R A + RRe c ) + j ( X A + X Re c ) (5.12)
Combining (5.9), (5.11), and (5.12), and using the fact that the maxi-
mum delivered power is the transferred power of the Friis transmission equa-
tion one obtains Ae as:
G A λ2 4 R A R Re c
Ae =
4 π (R A + R Re c )2 + ( X A + X Re c )2 (5.13)
Ae =
G A λ2
4π
(
1 − Γ2 ) (5.14)
where the reflection coefficient is as defined in (5.10). Note that (5.10) is dif-
ferent than (5.4) by a factor of impedance mismatch; in other words, (5.5) does
not include the mismatch factor.
112 A Practical Guide to EMC Engineering
It presents the amount of voltage captured by the antenna (note that VRec
and VA are at the same junction if the antenna is connected directly to the
receiver without any intervening cable). As a logarithmic value, AF is given as
G r λ2
PA = Pd × Ae [ W ], Ae = m 2 (5.17)
4π
The power density can also be calculated if the electric field in the vicinity
of the antenna is known:
E2 E2
Pd = = W m2 (5.18)
η0 120 π
where η0 = 120π = 377Ω is the free-space wave impedance. This is valid un-
der far-field conditions (i.e., transmit antenna that causes this EM wave is far
from the receiver; this could be taken as the distance which satisfies d ≥ λ/2π
for simple and low-gain antennas. Use d ≥ 2D2/λ for the high-gain complex
antennas; D is the longest dimension of the antenna).
The antenna is usually connected to a receiver (e.g., an EMI device) over
a finite-length TL. Assuming (1) the frequency-independent receiver input im-
pedance is 50 Ω; (2) the characteristic impedance of loss-free TL is real and
Antennas and Antenna Calibration 113
equal to Z0 = 50Ω, allow the use of the simple circuit model depicted in Figure
5.2. In this case, the output power of the receive antenna (equivalently, the
input power of the receiver) PA can be given as
V A2 V A2 V A2 V A2
PA = = = = [W ] (5.19)
R Re c R A Z 0 50
V A2 E2 G λ2
== × r [W ] (5.20)
R Re c 120 π 4π
9.73
AF = m −1 (5.22)
λ Gr
or, in dB:
η0
AF = m −1 (5.25)
R Re c Ae
114 A Practical Guide to EMC Engineering
Pt Gt
Pd = W m2 (5.26)
2
4πd
The same power density can be derived from the electric field at that
point and the wave impedance as given in (5.18); therefore, the combination of
(5.18) and (5.26) yields
30Pt Gt
E=
d
[ V m] (5.27)
The input power may be derived from Ohm’s law using voltage across its
terminals and its input resistance as
V A2
Pt = [ W ], or V A = Pt R A (5.28)
RA
E 1 30Gt
TAF = = m −1 (5.29)
VA d RA
Antennas and Antenna Calibration 115
1
TAF = 0.6Gt m −1 (5.30)
d
Note that (5.32) is valid only when one of AF-TAF is known and both
are obtained under the same conditions. For example, if AF is measured over a
ground plane, then TAF is requested over ground plane, too.
A typical AF graphic is given in Figure 5.6. Here, AFs of two different
models against frequency are plotted. As observed, AFs of different antennas
may be totally different. The important point is that AF (i.e., electric field–re-
ceive voltage transformation ratio) is known at every single frequency. If the
recorded voltage value is X dB and AF at this frequency is Y dB, then the mea-
sured electric field in the environment will be X + Y dB.
CALTS is an open field test site with metallic ground plane and tightly
specified site attenuation performance in horizontal and vertical polarization.
A CALTS is used for determining free-space antenna factor of an antenna.
Site attenuation measurements of a CALTS are used for comparison to
corresponding site attenuation measurement of a compliance test site, in
order to evaluate the performance of the compliance test site.
between the two antennas is recorded. Typically, three antennas are needed to
perform such a calibration, and they are measured in three pairings. Calcula-
tions are then performed to derive the AFs. Although a ground plane is used,
the aim of the ANSI/IEEE C63.5 SSM is to obtain the free-space AF by theo-
retically removing the ground plane effect. The conducting ground plane is
there to establish a repeatable calibration environment.
SSM assumes that (1) the antennas under test have radiation patterns of
a short dipole (i.e., horizontally uniform and figure “∞” vertically), (2) anten-
nas and the metal ground plane are all uncoupled and, (3) antennas are in the
far field so that the physical size of an antenna has negligible effect. However,
these are not always acceptable. Depending on types of antennas, error intro-
duced from the SSM can be as large as 2.5 dB [16]. This means the total error
could be 5 dB (because an antenna pair is used for the site attenuation). A new
ANSI/IEEE C63.5-2006 addresses this limitation by providing correction fac-
tors (CFs) [28], which are based on numerical simulations.
Reference antenna method (RAM) is basically a substitution method speci-
fied in ANSI/IEEE C63.5 [16]. Here, the responses between two known an-
tennas are measured, and then one is replaced by the antenna under test. The
AF is derived from the difference. Again, the mutual coupling between the two
standard antennas and the antenna under test, as well as ground plane effects,
can be significantly different, leading to significant errors.
The equivalent capacitance substitution method (ECSM) is used to cali-
brate monopole antennas from 9 kHz to 30 MHz. In this method, a dummy
antenna consisting of a capacitor equal to the self capacitance of the monopole
is used in place of the actual element (see 2.4 of IEEE Std 291-1991).
The other methods are standard field method (SFM) (IEEE Std 291-
1991), standard antenna method (SAM) (IEEE Std 149-1979), and standard
transmitting loop method (STLM), a specialized version of the SAM (IEEE Std
291-1991) [29]. One important note is that even if a perfect free-space AF has
been obtained, one should still apply the CFs provided by ANSI/IEEE C63.5
for the NSA test [13], because NSA tests are not in free space. The CFs are
used to correct the influences from the test setup (i.e., the differences between
free space and the specific geometry of an NSA setup). Finally, antennas with
high gain/narrow beamwidth (such as horns) do not see the ground plane when
placed in a close distance; therefore, calibration is carried out under free-space
conditions.
Antenna calibration is simply measuring AF.
High gain antennas have narrow radiation beams and therefore may be
calibrated in free space when elevated slightly above the ground because an el-
evated narrowbeam antenna does not suffer from ground reflections.
118 A Practical Guide to EMC Engineering
NSA has become the standard method for determining the adequacy of
OFTS to perform EM emission and immunity measurements. NSA is defined
as the ratio of the power input of a matched, balanced, lossless, tuned dipole ra-
diator to that at the output of a similarly matched, balanced, lossless, tuned di-
pole receiving antenna for specified polarization, separation, and heights above
a flat electromagnetically reflecting surface [16]. It is a measure of the transmis-
sion path loss between two antennas.
The theoretical method proposed by ANSI/IEEE for the evaluation of
NSA is obtained with a geometrical-optics (GO) approximation based on
the Friis equation since Norton surface waves are negligible for these antenna
heights and frequency bands [21, 29]. Standards define the NSA as the available
input power to a short dipole for the field strength of 100 µV/m at a distance d
= 3m over a leveled conducting ground screen [29].
279.1
NSA = (5.33)
f (MHz )E D max ( µV /m )
where
are horizontal and vertical electric field components [24]. The other parameters
in these equations are:
h + hr
g = tan −1 t
d
(reflection angle ) (5.38)
( εr − j 60σλ)
( εr − j 60σλ) sin g −
− cos 2 g
ρv = ( ver reflection coefficient ) (5.40)
( εr − j 60σλ)
( εr − j 60σλ) sin g +
− cos 2 g
2π Im { ρh ,v }
β= ; φh ,v = tan −1
Re { ρh ,v }
λ (5.41)
In (5.33), EDmax is the maximum electric field strength (in µV/m) for the
receive antenna height scanned between 1m and 4m from a theoretical half-
wave dipole with 1 pW radiated power and 1.64 dB antenna gain calculated
from (5.34) and (5.35) for horizontal and vertical polarizations, respectively.
The calibration steps are as follows:
Table 5.1
MATLAB Code for Theoretical NSA Calculations
%---------------------- Program : LS_NSA.m--------------------------
clear all; clc;
fprintf(‘============================== \n’);
fprintf(‘(1) Ver Pol (TM) \n’); fprintf(‘(2) Hor Pol (TE) \n’);
fprintf(‘============================== \n’);
nx = input (‘Your Choice? : ‘); ht = 1; d = 10; er = 15; sigma = 0.0; k=0;
for freq=3e7:1e7:1e9 % Frequency loop
k=k+1; k0=2*pi*freq/3e8; lambda=3e8/freq; Edmax(k)=1e-6;
Hrmax(k)=1; ep0=8.854e-12; n2=er-j*60*lambda*sigma;
for hr=1:.02:4 % Receive antenna height
d1=sqrt(d^2+(ht-hr)^2); d2=sqrt(d^2+(ht+hr)^2); theta=(atan(ht+hr)/d);
if nx==1 % VER Polarizasyon
rho=(n2*sin(theta)-sqrt(n2-cos(theta)^2))/(n2*sin(theta)+sqrt(n2-cos(theta)^2));
rhom=abs(rho); rhoa=angle(rho); rh(k)=rho;
Ed=sqrt(49.2)*d^2*(d2^6+d1^6*rhom^2+2*d1^3*d2^3*...
rhom*cos(rhoa-k0*(d2-d1)))^(0.5)/(d1^3*d2^3);
elseif nx==2 % HOR Polarizasyon
rho=(sin(theta)-sqrt(n2-cos(theta)^2))/(sin(theta)+sqrt(n2-cos(theta)^2));
rhom=abs(rho); rhoa=angle(rho); rh(k)=rho;
Ed=sqrt(49.2)*(d2^2+d1^2*rhom^2+2*d1*d2*...
rhom*cos(rhoa-k0*(d2-d1)))^(0.5)/(d1*d2);
end
if Edmax(k)<Ed; Edmax(k)=Ed; Hrmax(k)=hr; end
end
f(k)=freq/1e6; nsa(k)=279.1/(f(k)*Edmax(k));
end
figure(1); plot(f,20*log10(nsa),’LineWidth’,2);
xlim([30,1000]); ylim([-20,30]); grid on;
xlabel(‘Frequency [MHz]’); ylabel(‘NSA [dB]’); % Program END
Table 5.2
NSA Versus Frequency over Metal
Ground (Ver Polarization, d = 10m,
ht = 1m)
Frequency Edmax NSA
[MHz] [dBµ/m] [dB]
30 2.671 16.702
50 2.641 12.295
80 2.567 8.286
100 2.499 6.416
200 1.917 0.978
400 1.009 -4.135
500 1.668 -6.732
600 2.021 -8.669
800 2.368 -11.515
900 2.457 -12.627
1000 2.523 -13.608
122 A Practical Guide to EMC Engineering
Figure 5.8 A typical NSA versus frequency variation over a metal ground.
5Z 0 1 d
NSA =
2 π f (MHz ) −4 (5.42)
1 − ( βd ) + ( βd )
−2
for short antennas. In the far field, 1/d2 and 1/d4 terms in (5.42) can be ne-
glected and the free-space theoretical NSA equation becomes:
5Z d
NSA = 0
2 π f (MHz ) (5.43)
The logarithmic form of (5.43) for impedance matched (Z0 = 50Ω) sys-
tems will then be:
Table 5.3 lists NSA values computed from (5.9), (5.10), and (5.11) in the
free space. The free-space NSA from (5.9) may be obtained directly by setting
relative permittivity and conductivity values to one and zero, respectively.
As observed, using (5.11) instead of (5.10) introduces significant errors
that are not acceptable [20] for short range measurements. For d = 3m, this
error is 1 dB and maximum at 30 MHz. This error is less than 0.1 dB above
60 MHz at d = 5m, and above 110 MHz at d = 10m. Therefore, the free-space
NSA equation with only term can be used for EMC and AF measurements in
the frequency range of 30 MHz–1 GHz at d = 10m and d = 30m, but a CF
Antennas and Antenna Calibration 123
Table 5.3
NSA Versus Frequency in Free Space NSA1
(5.33), NSA2 (5.42), NSA3 (5.43), Ver-pol,
d=3m, hv=2m)
Frequency NSA1 NSA2 NSA3
[MHz] [dB] [dB] [dB]
30 12.00 12.98 12.00
50 7.56 7.97 7.56
70 4.64 4.94 4.64
90 2.45 2.59 2.46
110 0.71 0.80 0.71
130 –0.74 –0.68 –0.74
150 –1.98 –1.94 –1.98
170 –3.07 –3.03 –3.07
190 –4.04 –4.01 –4.03
200 –4.48 –4.46 –4.48
should be added at shorter ranges for lower end of this frequency band. The
CF includes not only the effects of 1/d 2 and 1/d 4 terms, but also nonuniform
radiation of receive antenna, mutual antenna couplings, and mutual coupling
between antennas and the ground [20]. Further investigation of theoretical
NSA models including Norton surface wave contributions is presented in [19].
where AFT and AFR are antenna factors of the transmit and receive antennas (in
dBm–1), and AFTOT is the mutual impedance CF (in dB). Note that the dif-
ference VDirect – VSite is equal to the classic site attenuation.
124 A Practical Guide to EMC Engineering
The single-point NSA measurement stated earlier is valid only for OFTS.
For the NSA measurements in a semi-anechoic chamber, a single-point NSA
measurement is insufficient to pick up possible reflections from the construc-
tion or RF-absorbing materials comprising the walls and ceiling of the test site.
For anechoic chambers, a test volume is defined [27] as that volume traced out
by the largest equipment or system to be tested as it is rotated about its center
location through 360°, such as by a turntable. In evaluating the site, the trans-
mit antenna should be placed at various points within the test volume with
both horizontal and vertical polarizations, such as illustrated in Figure 5.9. This
may require a maximum of 20 separate site attenuation measurements—that is,
five positions in the horizontal plane (center, left, right, front, and rear, mea-
sured with respect to the center and a line drawn from the center to the position
of the measuring antenna), for two polarizations (horizontal and vertical), and
for two heights (1m and 2m, horizontal) (1m and 1.5m, vertical) [24, 27].
These measurements are carried out with broadband antennas. The dis-
tance d is measured from the center of transmit antenna to the center of receive
antenna. This distance shall be maintained for all measurements, which requires
that the receive antenna be moved along the line in the directions shown in the
figure. Also, transmit and receive antennas shall be aligned with the antenna
elements parallel to each other and orthogonal to the measurement axis. Fur-
thermore, the lower tip of the antenna shall be at a distance greater than 25 cm
from the floor, which may require the center of the antenna to be slightly higher
than 1m for the lowest height measurement.
Both CISPR and ANSI/IEEE necessitates NSA of an OFTS must be
within ±4 dB of the theoretical values in order to accredit a CALTS. It is rec-
ommended that this should be within ±1 dB to reduce the measurement un-
certainty significantly. There are currently no CALTS validation requirements
above 1 GHz. However, facilities suitable for measurements in the frequency
range 30 MHz to 1000 MHz are considered suitable for the frequency range 1
GHz to 40 GHz, including the presence of the reference (metal) ground plane.
• There are three or more antennas with the same frequency band; how-
ever, they may or may not be identical.
• There are two antennas with the same frequency band; however, they
may or may not be identical.
• There are two nonidentical antennas with the same frequency band. In
addition, one of them already has a recent calibration certificate from an
accredited lab.
• Antenna-cable mutual coupling and scattering from the cable are neg-
ligible.
• Horizontally polarized ground reflection is less sensitive to differences
in the ground plane, as well as to differences in the ground conductivity
and permittivity.
• Ground screen edge diffractions are smaller for horizontal polarization.
126 A Practical Guide to EMC Engineering
The TAC method requires three site attenuation (A) measurements under
identical geometries (ht, hr, d), using three antennas at hand in pairs. The three
equations associated with these measurements are:
Here, AF1, AF2, and AF3 are the antenna factors of antennas 1, 2, and 3
(in dBm-1) to be calculated; A1, A2, and A3 are the measured site attenuations
(in dB). Rearranging these equations yields the AF (in dBm-1) of the three an-
tennas as:
E D max + A1 + A2 − A3
AF1 = 10 log ( f MHz ) − 24.46 + (5.49)
2
E D max + A1 + A3 − A2
AF2 = 10 log ( f MHz ) − 24.46 + (5.50)
2
E D max + A2 + A3 − A1
AF3 = 10 log ( f MHz ) − 24.46 + (5.51)
2
E D max + A
AF = 10 log ( f MHz ) − 24.46 + (5.52)
2
In practice, two antennas are never identical, and the AF calculated from
(5.20) is the geometric mean (in linear units) of the individual factor for each of
the two antennas. Certain antennas can be constructed to be so nearly identical
that their AFs are different by much less than the measurement uncertainty.
If the two antennas are not identical, an antenna that has the same frequency
range can be borrowed from another EMC lab to do the calibration by using
the TAC method. If there are two antennas with the same frequency band and
one of them already has a recent calibration certificate from an accredited lab,
then site attenuation A can be measured and then the AF of unknown antenna
can be calculated from
where AF2 is the AF of the certificated antenna. Table 5.4 lists a typical antenna
calibration result of a Schaffner BHA9118 horn antenna.
Table 5.4
Typical Antenna Calibration Results
ANSI/IEEE
BHA 9118 ANSI/IEEE C63.5 149|G| [dB]
Horn Schaffner Specs TAC Method Method Difference
Freq. [GHz] AF [dB] |G| [dB] AF [dB] |G| [dB] |G| [dB] ∆G [dB]
1 24.37 5.83 24.50 5.71 5.73 0.022
2 27.28 8.94 27.09 9.14 9.16 0.022
3 29.94 9.80 29.87 9.87 9.91 0.024
4 32.62 9.62 32.45 9.80 9.83 0.023
5 33.01 11.16 33.76 10.43 10.45 0.022
6 34.27 11.49 34.48 11.29 11.31 0.022
utilized. However, the advantage of the transmit mode is to find the possibility
to measure the absolute radiated power in addition to the pattern, which in
turn also yields directivity and gain.
In addition to the instruments listed earlier in the TAC method, a cali-
brated receive antenna in the same frequency band, a positioning system for θ
and φ variations, a recording system, and data processing software are also re-
quired for the APM method. Radiation patterns of most practical antennas are
so complex that closed form mathematical expressions are not available.
The radiated power of an antenna can be calculated from the double inte-
gral over a sphere enclosing the antenna if the 3D directional radiation function
F(θ, φ) is known. In discrete form, this integral reduces to a double summation:
M N
Prad = Bo ∆θ∆φ ∑ ∑ F ( θi , ϕ j )sin θi (5.54)
j =1 i =1
Antennas and Antenna Calibration 129
(B0 is constant) and can be used to calculate radiated power once the pat-
tern F(θi, φi) is measured at a number of points θi and φi:
π π π
N = , θi = i ∆θ = + (i − 1) , i = 1,2,, N (5.55)
∆θ N N
2π 2π 2π
M= , φ j = j ∆φ = + (i − 1) , j = 1,2,, M (5.56)
∆φ M M
It has been tested with a number of different antenna patterns that radia-
tion pattern discretization error (i.e., error introduced using (5.54) instead of
using the mathematical exact double-integral) is less than 0.5 dB if N and M are
between 10–15 and 20–30, respectively. Once the radiated power is found the
directivity at any direction can be calculated from [12]
4 πU ( θ, φ)
D ( θ, φ) = (5.57)
Prad
where PFwd and PRvs are the powers measured at ports 1 and 2 as shown in the
figure, respectively. Absolute gain of the antenna (in dB) is then equal to
130 A Practical Guide to EMC Engineering
Table 5.5
Uncertainty Budget for Antenna Calibration
Ci Ui
Source of Probability (Sensitivity (Uncertainty
Sym Uncertainty Unit Value Distribution Divisor Coefficient) Component) U i2
dM Model Error dB 0.2 Rectangular 1,732 1 0.115 0.0132
dY Site reflections dB 0.2 Rectangular 1,732 1 0.115 0.0132
dIG Ap dB 0.1 Normal 2 2 1 0.05 0.0025
dIF dB 0.1 Rectangular 1,732 1 0.06 0.0036
dAG dB 0.2 Normal 2 2 1 0.100 0.01
dFH dB 0.2 Rectangular √3 1 0.115 0.0132
dP dB 0.1 Rectangular √3 1 0.060 0.004
dCP dB 0.1 Rectangular √3 1 0.060 0.004
dMS dB 0.3 Rectangular √3 1 0.173 0.030
dA dB 0.2 Rectangular √3 1 0.115 0.0132
dB dB 0.2 Rectangular 1 0.115 0.0132
Antennas and Antenna Calibration
√3
dAF dB 1 Normal 2 2 1 0.500 0.250
dKK dB 0.2 Normal 2 2 1 0.100 0.010
dU dB 0.5 U-Shaped √2 1 0.353 0.125
dYK dB 0.1 Rectangular √3 1 0.060 0.004
dUYK dB 0.1 U-Shaped √2 1 0.070 0.005
dT dB 0.5 Normal 1 1 1 0.500 0.250
Combined standard uncertainty: 0.87 dB. Expanded uncertainty: 1.74 dB
131
132 A Practical Guide to EMC Engineering
experimentally shown in the following sections. The factors of the RE101 an-
tenna are often presented in unit of dBpT/µV since the limits in the standard
are in the unit of dBpT. The RE101 test setup is shown in Figure 5.13.
Another test in MIL-STD 461E/F that requires the use of loop antennas
is the RS101 test. This test defines two loop antennas. Unlike RE101 test, this
is an immunity test and aims to expose the EUT to magnetic radiated emis-
sions in the frequency range 30 Hz–100 kHz. One of required antennas by
the RS101 test is a transmitting loop antenna with 20 turns and a diameter
of 12 cm for exposing EUTs to magnetic fields. This loop is also chosen as a
transmitter by SAE-ARP 958 standard for the calibrations of RE101 receiving
loop antennas. The other antenna required by the RS101 test is a receiver loop
with 51 turns and a diameter of 4 cm for verifying the test setup before the
actual test. Unlike other immunity tests like RS103, during the RS101 test, no
feedback process is applied to measure the produced magnetic field, as shown
in Figure 5.14.
The current that flows on the transmitter loop is measured by means of
a current probe, and produced magnetic field on the EUT is then calculated
using this current. The precise measurement of the current therefore becomes
essential for the reliability of the test. Before the actual RS101 test, a verification
process must be applied to verify all the test setup by using the receiving loop
as shown in Figure 5.15.
In addition to military applications, loop antennas are commonly utilized
in SE measurements below 30 MHz [37] (see Figure 5.16). In SE measure-
ments, the calibration factors are not needed. The comparison of two measure-
ment cases with and without the target shield, as in Figure 5.16, yields the
desired SE value.
where A, N, f, and B are the area of the loop in square meters, the number of
turns on the coil, the frequency in Hz, and the magnetic flux in Tesla, respec-
tively. The voltage VL across the load RL , when the loop antenna (with the coil
impedance ZC = RC + jXC) is loaded, will be
136 A Practical Guide to EMC Engineering
VL VOC
= (5.63)
RL RL + RC + jX C
B (1 + RC / RL )2 + (2 πfL / RL )2
AFH = = (5.64)
VL 2 πfAN
B [T ] = µo H [ A / m ] = 4 π 10 −7 H [ A / m ]
B [ pT ] = 4 π 10 −7 H [ A / m ]1012 = 4 π H [ A / m ]105
B [ pT ] H [A / m ] S
AFH [dBpT / µV ] = = 1.256 = 1.256( )
V [ µV ] V [V ] m
Hence, one subtracts 1.98 dB from dBpT/µV and obtains dBS/m. Some-
times, AFH is expressed in terms of equivalent electric field intensity (AFE =
E[µV/m]/V[µV]). In this case, the unit of the AF is dBm–1 (i.e.,AFE[dBm–1]
= AFE[dBS/m] + 51.5). A typical AF graph is given in Figure 5.18. In order
to obtain magnetic flux density at a given frequency, AFH value read from the
Antennas and Antenna Calibration 137
graphics and the dBµV value read from the EMI receiver are added. For ex-
ample, at 10 kHz, AFH is read from the graphics as 35 dBpT/µV. If the EMI
receiver shows 23 dBµV, then the measured magnetic flux density will be 58
dBpT.
Details of the loop antenna calibration according to IEEE 291 [33] and
SAE-ARP 958 [34] are given in [38].
horizontal and vertical. For example, in order to plot horizontal radiation pat-
tern, θ is kept constant (θ = 90°) and φ is changed from 0° to 360° with a given
angular resolution ∆φ. Similarly, φ angle is kept constant (φ = 0° for xz-plane
or φ = 90° for yz-plane) and θ angle is rotated 360° with an angular resolution
of ∆θ.
N
E ( θ, ϕ) = ∑ I i e j ψi ( θ, ϕ) (5.65)
i =1
Here, Ii is the complex current for the ith element and ψi is phase con-
tribution at the far field point from the ith radiator with respect to the origin,
which can be expressed as:
Antennas and Antenna Calibration 139
r
ψi = k ri ⋅ (5.66)
r
For planar arrays located on the xy-plane, the radiators’ coordinates are xi = ri
cosφι, yi = ri sinθi, and zi = 0, where ri = x i2 + yi2 and φi = a tan–1(yi /xi). For
the planar arrays (i.e., θi = 90°) the total electric field, which is known as the
array factor, can be expressed as:
N
E ( θ, ϕ) = ∑ I i exp ( jkri sin θ cos( ϕi − ϕ)) (5.68)
i =1
N
E ( θ, ϕ) = ∑ I i exp jkri [ sin θ cos( ϕi − ϕ) − sin θ0 cos( ϕi − ϕ0 )] (5.69)
i =1
140 A Practical Guide to EMC Engineering
if maximum radiation is desired to point (θ0, φ0). Figure 5.21 shows broadside
and endfire radiation patterns obtained with the same array of five isotropic
radiators; only the phases of the radiators are changed accordingly.
and Theta/Phi value box. If the θ-plane is selected (vertical pattern), the angle
θ varies between [0°, 360°] at a fixed φ value (specified by the Theta/Phi box).
If the φ-plane (horizontal pattern) is selected, the angle φ varies between [0°,
360°] at a fixed θ-value specified by the user. All the 2D radiation patterns are
plotted with an angular resolution of ∆θ = ∆φ = 1°. Figure 5.22 presents a
circular array with maximum radiation along (θ0 = 45°, φ0 = 45°) direction.
Figure 5.23 shows the vertical radiation pattern of 10×10 element planar
array of isotropic radiators. Note that radiation patterns may also be plotted in
3D [39].
Another exciting feature of ARRAY is that the user can use sliding bars to
change N, f, θ0, φ0. For instance, by pressing the sliding bar of the frequency
continuously, the user can increase or decrease the frequency, see the corre-
sponding radiation pattern immediately, and may easily understand the effect
of frequency change on the radiation pattern. The package also allows the user
to save the radiation pattern data (field intensity versus angle) to a file named
Arraypattern.dat. The first column of this file corresponds to 361 observation
angle values (in radians), and the corresponding array factors are in the second
column.
Figure 5.23 A 10×10 planar array and its radiation pattern on the yz-plane.
Extra line phasing capability has been included in ARRAY to change FBR
characteristics in the planar array option. The user may change the phase angle
∆β of the lined-up radiators and observe how FBR improves or gets worse. In
this case complex currents Ii in (5.65) is replaced with Iin as
M
I i → I in = ∑ I i exp { j (n − 1)∆β} (5.70)
n =1
large FBR is required. This may be achieved by phasing the second row radia-
tors (and the others if there are more than two rows) with respect to the first
row. This is illustrated at right in Figure 5.24, where a 90° phase-difference is
given for the second row radiators. As observed in the figure, a very good FBR
is achieved.
A short MATLAB code is also given in Table 5.6. This code plots radia-
tion pattern of a circular array of isotropic radiators. The user supplies num-
ber of radiators, radius of the circle, the frequency, and the beam direction. A
sample pattern plotted with this code is shown in Figure 5.25.
Note that side lobes of an array may also be reduced by changing the
radiators’ feeding amplitudes [8–12]. Two examples are given in Figures 5.26
and 5.27. Here, 1D and 2D arrays of isotropic radiators are fed with amplitude
values like in binomial series [12]. Binomial feeding is best in eliminating side
lobes; the payoff is the broader main lobe.
Table 5.6
MATLAB Code for a Circular Array Pattern Plot
%---------------------------------------------------------------------
% Program: Circular antenna array
%---------------------------------------------------------------------
clc; close all; clear all;
N = 4; % input(‘Number of radiators = ?’);
r = 0.25; %input(‘Radius [m]= ?’);
fr = 150; % input(‘Frequency [MHz] = ?’);
phi0 = 90; %input(‘Beam direction [degree]= ?’);
c=3e8; fr=fr*1e6; k=2*pi*fr/c; dtr = pi/180;
rtd = 180/pi; phi0r = phi0*dtr;
for i=1:N
phi=2*pi*(i-1)/N; phin(i)=phi; xn(i)=r*cos(phi); yn(i)=r*sin(phi);
end
i = 0;
for phi = 0:1:359
i = i+1; phir(i) = phi*dtr; Af(i) = 0;
for n = 1:N
Af(i)=Af(i)+exp(1i*k*r*(cos(phir(i)-phin(n))-cos(phi0r-phin(n))));
end
Af(i) = abs(Af(i));
end
Array = Af/max(Af);
polar(phir,Array); hold on;
plot(xn*0.25/max(xn),yn*0.25/max(xn),’ro’)
%--------------------PROGRAM END-------------------------
The other two basic antennas that are widely used in antenna engineering
and important in EMC engineering are the electric and magnetic dipoles.
Antennas and Antenna Calibration 145
Figure 5.26 Side lobe reduction using binomial feeding: (left) 1-1-1-1-1-1-1-1; (right) 1-7-21-
35-35-21-7-1.
Figure 5.27 Side and back lobe reduction using binomial feeding.
Characteristics of electric and magnetic dipoles are just the opposite in the
near field. In the far field (once the contributions of r –2 and r –3 terms become
negligible), they are both assumed as a plane wave. The impedances of electric
and magnetic dipoles versus normalized distance are plotted in Figure 5.32. As
observed, the wave impedance of an electric dipole is very high (5 kΩ–10 kΩ),
but it is very low for a magnetic dipole (10Ω–30Ω) in the near field. Around
and beyond the normalized distance of λ/2π [m], wave impedances of both
types become free-space wave impedance (η0 = Z0 = |Eθ|/|Hφ| = 377 Ω).
In EMC engineering, the theoretical boundary of the near field is ac-
cepted as λ/2π.
Table 5.7
Short MATLAB Code of an LPDA Design
%------------------------------------------------
% Program: LPDA.m
%------------------------------------------------
clc; clear all;
sig = [.139.149.158.160.163.168.172.174.176.177.180];
ta = [.782.822.843.868.882.900.920.952.948.956.949];
% Supply design parameters
fmin = input(‘Lower Freq [MHz]: ‘); fmax = input(‘Upper Freq [MHz]: ‘);
G = input(‘Enter Gain [6-11 dB]: ‘);
while G<6 || G>11; G = input(‘Enter GAIN [6-11 dB]: ‘); end
c = 300; sigma = sig(2*G-11); tau = ta(2*G-11);
alpha = 2 * atan((1-tau) / (4 * sigma)) * 180/pi;
lambda_lower = c / fmin; lambda_upper = c / fmax;
N = 1; L(N) = 0.5 * lambda_lower; Lu = 0.5 * lambda_upper;
fprintf(1, ‘\n Alpha degree is: %2.1f\n’, alpha);
fprintf(1, ‘ Tau for the optimum design: %1.3f\n’, tau);
fprintf(1, ‘ Sigma for the optimum design: %1.3f\n\n’, sigma);
while L(N) > Lu; L(N + 1) = tau * L(N); N = N + 1; end; fprintf(1, ‘ \n’)
for n = 1 : N-1
d(n) = 2 * sigma * L(n);
fprintf(1, ‘Dipole # %2.0f-- L = %3.3f-- d = %3.3f \n’, n, L(n), d(n)); end % Program END
Figure 5.36 A sample LPDA designed with the code in Table 5.7.
Figure 5.37 (a) Vertical, and (b) horizontal radiation patterns of the LPDA at 800 MHz.
5.38). Horns are widely used in EMC tests and measurements above 250–300
MHz frequencies. They can operate over a wide range of frequencies (e.g., from
1 GHz to 20 GHz) because they have no resonant elements. Their gains range
up to 25 dBi with 10–20 dBi being typical. Horns have very little loss; their
directivity are roughly equal to their gain and therefore are also used in antenna
calibration.
References
[1] Sevgi, L., “The Antenna as a Transducer: Simple Circuit and Electromagnetic Models,”
IEEE Antennas and Propagation Magazine, Vol. 49, 6, Dec. 2007, pp. 211–218.
[2] Sevgi, L., “Reciprocity: Some Remarks from Field Point of View,” IEEE Antennas and
Propagation Magazine, Vol. 52, 2, Apr. 2010, pp. 205–210.
[3] IEEE Standard Definitions of Terms for Antennas, IEEE Std 145–1983.
[4] King, R. P. W., The Theory of Linear Antennas, Cambridge, MA: Harvard University Press,
1956.
[5] Jordan, E., and K. Balmain, Electromagnetic Waves and Radiating Systems, 2nd Ed., Engle-
wood Cliffs, NJ: Prentice Hall, 1968.
[6] Elliot, R. S., Antenna Theory and Design, Englewood Cliffs, NJ: Prentice Hall, 1981.
[7] Colin, R. E., Antennas and Radiowave Propagation, New York: McGraw Hill, 1985.
[8] Milligan, Thomas A., Modern Antenna Design, New York: McGraw Hill, 1985.
[9] Kraus, J. D., Antennas, Second Ed., New York: McGraw Hill, 1988.
[10] Brookner, E., Practical Phased Array Systems, Norwood, MA: Artech House, 1991.
[11] Manser, G., Antenna Design: A Practical Guide, New York: McGraw Hill, 1996.
[12] Balanis, C. A., Antenna Theory, 2nd Ed., New York: John Wiley and Sons, 1997.
Antennas and Antenna Calibration 153
[13] Hansen, R. C., Phased Array Antennas, New York: John Wiley and Sons, 1997.
[14] Shutzman, W. L., and G. A. Thiele, Antenna Theory and Design, 2nd Edition, New York:
John Wiley and Sons, 1998.
[15] Macnamara, Thereza M., Handbook of EMC Antennas, Norwood, MA: Artech House,
1995.
[16] ANSI C63.5-1998, American National Standard for Electromagnetic Compatibility-
Radiated Emission Measurements in Electromagnetic Interference (EMI) Control
Calibration of Antennas (9kHz–40 GHz).
[17] Sevgi, L., S. Çakır, and G. Çakır, “Antenna Calibration for EMC Tests and Measurements,”
IEEE Antennas and Propagation Magazine, Vol. 50, No. 3, Jun. 2008, pp. 215–224.
[18] Smith, A. A., R. F. German, and J. B. Pate, “Calculation of Site Attenuation from Antenna
Factors,” IEEE Trans on EMC, Vol. 24, No. 3, Aug. 1982, pp. 301–316.
[19] Trakadas, P. T., and C. N. Capsalis, “A Mixed Model for the Determination of Normalized
Site Attenuation in OATS,” IEEE Trans on EMC, Vol. 43, No. 1, Feb. 2001, pp. 29–36.
[20] Gavenda, J. D., “Near Field Corrections to Site Attenuation,” IEEE Trans on EMC, Vol.
36, No. 3, Aug. 1994, pp. 213–220.
[21] Askri, A., C. Vollaire, L. Nicolas, and D. Prebet, “Normalized Site Attenuation Standard
Correction from Numerical Computing,” IEEE Trans on EMC, Vol. 38, No. 2, Mar.
2002, pp. 693–696.
[22] Agilent Technologies, http://www.agilent.com.
[23] TDK RF Solutions, www.tdkrfsolutions.com.
[24] CISPR16-1-5:2003 Antenna calibration test sites (CALTS) for 30 MHz to 1000 MHz.
[25] CISPR 16-1-4, “Specification for Radio Interference Measuring Apparatus and
Measurement Methods,” IEC 2004.
[26] CISPR16-4-2:2003 Measuring Apparatus and Methods—Uncertainties, Statistics, and
Limit Modeling—Uncertainty in EMC Measurements.
[27] ANSI/IEEE Standard C63.4-1992, American National Standard for Methods of
Measurement of Radio-Noise Emissions from Low-Voltage Electrical and Electronic
Equipment in the Range of 9 kHz to 40 GHz.
[28] ANSI/IEEE C63.5-2006 (Revision of ANSI C63.5-2003), American National Standard
for Electromagnetic Compatibility–Radiated Emission Measurements in Electromagnetic
Interference (EMI) Control–Calibration of Antennas (9 kHz to 40 GHz).
[29] ANSI/IEEE Std 149-1979 (Revision of IEEE Std 149-1965) IEEE Standard Test
Procedures for Antennas.
[30] IEEE Standard 145-1983, IEEE Standard Definitions of Terms for Antennas.
[31] ISO/IEC/OIML/BIPM, Guide to the Expression of Uncertainty in Measurement, Prepared
by ISO/TAG 4/WG 3, January 1993.
[32] Eser, S., and L. Sevgi, “Open Area Test Site (OATS) Calibration,” IEEE Antennas and
Propagation Magazine, Vol. 52, No. 3, Jun. 2010, pp. 204–212.
154 A Practical Guide to EMC Engineering
[33] IEEE Std 291-1991, IEEE Standard Methods for Measuring Electromagnetic Field
Strength of Sinusoidal Continuous Waves, 30 Hz to 30 GHz.
[34] SAE-ARP 958 Rev. D, Electromagnetic Interference Measurement Antennas; Standard
Calibration Method.
[35] MIL-STD-461E-1999, Interference Standards Requirements for the Control of
Electromagnetic Interference Characteristics of Subsystems and Equipment, US
Department of Defense.
[36] MIL-STD-461F-2007, Interference Standards Requirements for the Control of
Electromagnetic Interference Characteristics of Subsystems and Equipment, US
Department of Defense.
[37] CENELEC EN50147-1: 1996, Anechoic Chambers Part 1: Shield Attenuation
Measurement.
[38] Çakır, S., R. Hamid, and L. Sevgi, “Loop Antenna Calibration,” IEEE Antennas and
Propagation Magazine, Vol. 53, No. 4, Oct. 2011, pp. 243–254.
[39] Sevgi, L., and Ç. Uluışık, “A MATLAB-Based Visualization Package for Planar Arrays of
Isotropic Radiators,” IEEE Antennas and Propagation Magazine, Vol. 47, No. 1, Feb. 2005,
pp. 156–163.
Bibliography
Paul, C. R., Introduction to Electromagnetic Compatibility, New York: Wiley Interscience, 1992.
McLeon, J., R. Sutton, and R. Hoffman, “Interpreting Antenna Performance Parameters for
EMC Applications,” Application Notes, TDK RF Solutions Inc., www.tdkrfsolutions.com.
IEEE Std 1309-2005, IEEE Standard for Calibration of Electromagnetic Field Sensors and
Probes, Excluding Antenna, from 9 kHz to 40 GHz.
6
Noise and Frequency Analysis
Suppose you are a 100m sprinter. You run on your own track; therefore, you
do not to worry about being obscured by the next runners. If you are a 3000m
runner, then the situation is different: after a short start on your track, like the
other runners, you prefer to continue on the innermost track. In that case, a
runner can obscure others. EMC problems are similar to this example. Two
or three decades ago, there was no EMI problem between a long-wave radio
transmitter and defense radar because their locations and frequencies were far
enough apart. Today, the electrical and electronics systems we use in our daily
lives (e.g., a laptop, a cell phone, hundreds of systems in automobiles, and fast-
trains) are all EMI-vulnerable because they are digital, broadband, and have to
operate close to each other.
Two parameters are important for all equipment operating in close prox-
imity: signal strength and frequency. In order not to interfere with each other
here are two simple EMC engineering tips:
• If frequencies of two pieces of equipment are very close, then the dis-
tance in between them should be far enough apart or some sort of isola-
tion should be applied.
• If two pieces of equipment need to operate very close to each other, they
must use different frequencies.
155
156 A Practical Guide to EMC Engineering
6.2 Noise
Noise [5] is an uncorrelated signal that limits signal detection. Figure 6.1 pres-
ents typical noise spectra with manmade or natural components (sun burst,
magnetic storms, cosmic noise, atmospheric effects, etc.).
In terms of EMC engineering, noise can be grouped into two types: envi-
ronmental noise and thermal noise (also called internal noise, electronic noise,
Gaussian noise, and white noise). The thick line in the figure belongs to ther-
mal noise. As observed, above a few hundred MHz frequencies thermal noise
is dominant. Roughly speaking, EMC engineers do not need to worry about
environmental noise if frequencies are above 300 MHz–500 MHz; the only
noise they have to pay attention to is thermal noise. On the contrary, for the
frequencies below 200 MHz–300 MHz, thermal noise is negligible; EMC en-
gineers have to take care of the environmental noise. At high frequencies (above
500 MHz) thermal noise is critical in EMC engineering.
Thermal noise exists in all circuits and devices caused by random varia-
tions in current or voltage by random movement of charges. Thermal noise can
be reduced by reducing the temperature. This phenomenon limits the mini-
mum signal level that any receiver can usefully respond to, because there will
Noise and Frequency Analysis 157
always be a small but significant amount of thermal noise arising in its input
circuits. There are several other sources of noise in electronic circuits such as
shot noise, seen in very low-level signals where the finite number of energy-
carrying particles becomes significant, or flicker noise (1/f noise) in semicon-
ductor devices.
The sensitivity (minimum detectable signal) of an EMI receiver is deter-
mined by the noise level. Thermal noise is uncorrelated; this means it appears at
all frequencies. It is a random signal and therefore necessitates stochastic mod-
eling. Thermal noise is Gaussian with zero average value. It fluctuates around
useful signal. Noise power can roughly be predicted from
Pn = kTB [ W ] (6.1)
Vm2 V2
SNR = or SNRdB = 10Log 10 m (6.2)
2kTB 2kTB
Any of these three (SNR, signal voltage, and noise voltage) can be cal-
culated if the other two are known or given. The noise voltage, for example, is
obtained from
Vm
Vn = (6.3)
2SNR
Figure 6.2 SNR versus time of a sinusoidal signal with Gaussian noise (SNR=10 dB, Vm = 3V,
f = 100 Hz, T = 0.25s).
Noise and Frequency Analysis 159
1 ∞ − j ωt
S ( ω) = ∫ s (t )e dt
2 π −∞
(6.4)
Noise and Frequency Analysis 161
1 ∞
s (t ) = ∫
2 π −∞
S ( ω)e j ωt dt (6.5)
Here, s(t) is the time signal, S(ω) is its FT, and ω=2πf is the angular fre-
quency. FT is defined for continuous time (analog) signals. To obtain FT of a
time function at every frequency with infinitesimal frequency steps, the time
signal must be observed forever.
A few functions and their FT are given in Figure 6.4. FT of a triangle
function is square of a Sinc() function. FT of a Gauss pulse is also a Gauss pulse.
As the pulse gets narrower in the time domain, its FT gets broader. FT of an
impulse train is also an impulse train. As the time interval between impulses
gets narrower, the frequency intervals gets broader. For a single impulse, the FT
is constant. In communication society, it is called impulse response.
Knowing impulse response of a circuit means knowing all frequency
responses.
Fourier transform is used for energy signals (i.e., for signals which are non-
∞
2
zero in a certain period of time and goes to zero when t → ±∞ ( ∫ s (t ) dt is
−∞
finite)). Fourier series (FS) representation is used for periodic (power) signals. A
v(t) function in a finite T1 ≤ t ≤ T2can be represented with FS in terms of sine
and cosine functions:
A0 ∞ 2 πnt 2 πnt
v (t ) ≈ + ∑ An cos + Bn sin
T2 − T1
(6.6)
2 n =1 T2 − T1
T2 − T1 T∫
A0 = v (t )dt (6.7a)
1
T2
2 2 πnt
An =
T2 − T1 T∫ v (t )cos dt
T2 − T1 (6.7b)
1
T2
2 2 πnt
Bn =
T2 − T1 T∫ v (t )sin dt
T2 − T1 (6.7c)
1
Table 6.1
MATLAB Code for Fourier Series Summation of a Trapezoidal Pulse
�%-----------------------------------------------------------------------------------
% PROGRAM: Fourier Series representation of a trapezoidal function
%-----------------------------------------------------------------------------------
clear all; clc;
N=400; xf=0; xl=10; P=abs(xl-xf); delx=P/N; x0=2.5; xT=.4;
tr=0.5; tau=2; % Pulse rise time and duration
for j=1:N
xx(j)=xf+(j-1)*delx;
% Trapezoidal pulse
if xx(j)<=2*tau ; fnr(j)=0; elseif xx(j)<=2*tau+tr; fnr(j)=(xx(j)-2*tau)/tr;
elseif xx(j)<=3*tau+tr ; fnr(j)=1; elseif xx(j)<=3*tau+2*tr
fnr(j)=3-(xx(j)-3*tau+tr)/tr; else; fnr(j)=0; end; end
A0=trapz(xx,fnr)*2./P; % Calculating A0
for j=1:N; fn(j)=A0/2; end
OK=1; k=0; Nt=0;
fprintf(1, ‘--CoefficientNo---------An----------Bn---\n’);
fprintf(1, ‘------------------------------------------------------\n ‘);
while OK>0 % Fourier series of the function
Nt=Nt+1;
for k=1:N % Calculate A(k) and B(k)
g1(k)=fnr(k)*cos(2*pi*Nt*xx(k)/P); g2(k)=fnr(k)*sin(2*pi*Nt*xx(k)/P);
end
An=trapz(xx,g1)*2./P; Bn=trapz(xx,g2)*2./P; % Use trapz integral function
fprintf(1, ‘------------------------------------------------------\n ‘);
fprintf(1, ‘ %3d %12.8f %12.8f \n’, Nt, An, Bn);
fprintf(1, ‘------------------------------------------------------\n ‘);
for j=1:N % Calculate the function
fn(j)=fn(j)+An*cos(2*Nt*pi*xx(j)/P)+Bn*sin(2*Nt*pi*xx(j)/P);
end
plot(xx,fnr,’k’,xx,fn, ‘r--’,’LineWidth’,2); grid; xlabel(‘Time [s])’); ylabel(‘f(t)’);
title(ttle); legend(‘Function’,’FSeries’); pause; end % Program END
even-symmetric trapezoidal pulse with sharper rise and fall times. As observed,
the FS representation is still insufficient even after the first 80 terms.
Figure 6.5 Trapezoidal pulse and its FS representation: (a) 1 term; (b) 10 terms; (c) 18 terms.
Noise and Frequency Analysis 165
Figure 6.6 Symmetric trapezoidal pulse and its Fourier series representation.
frequency of the sinusoid differs in value from the frequency samples, which is
termed scalloping loss.
The DFT of a continuous time signal sampled over the period of T, with
a sampling rate of t can be given as
N −1
T
S (m ∆f ) =
N
∑ s (n ∆t )e − j 2 πm ∆f n ∆t (6.8)
n =0
1 1
f max = , ∆f = (6.9)
2 ∆t T
Since sampling interval and signal record lengths are finite in numeri-
cal computations maximum frequency and the resolution are also finite. This
means:
168 A Practical Guide to EMC Engineering
Tables 6.2 and 6.3 list short MATLAB codes for the computations of
DFT and FFT of given functions. Note that the code in Table 6.2 uses (6.4) in
discrete form, so any number of frequency samples () can be used to plot a fre-
quency spectrum. Unfortunately, (6.8) still holds for the DFT. In other words,
for example, if one wants to discriminate two sinusoids with 50 Hz and 52 Hz
in the frequency domain, the frequency resolution must be less than 2 Hz. This
necessitates minimum T = 0.5s. On the other hand, the DFT of a signal can be
taken using any ∆f, but the capability of discriminating two nearby sinusoids in
the frequency domain is limited by the length of the time record. The FFT code
in Table 6.3 uses built-in MATLAB commands fft, ifft, and fftshift.
A typical spectrum obtained with these codes (either DFT or FFT) is
presented in Figure 6.9. Here, two sine functions with 35 Hz and 50 Hz (s(t) =
sin(100πt) + sin(70πt)) are used as the time function. The sampling parameters
are ∆t = 10 ms and T = 1 s. As observed, only 35-Hz components appear in
the spectrum. The other (±50 Hz) components cannot be observed since they
Table 6.2
Short MATLAB DFT Code
%----------------------------------------------------------
% Program: DFT.m
%----------------------------------------------------------
clear all; close all; clc;
% Supply f1, f2, T, dt, df, fmax
a1=1; a2=0.5; N=T/dt; w1=2*pi*fr1; w2=2*pi*fr2; M=fmax/df;
% build the input time series
for k=1:N; t(k)=k*dt; st(k)=a1*sin(w1*dt*k)+a2*sin(w2*dt*k); end
for k=1:M % Apply the DFT with M points
Sf(k)=complex(0,0);
for n=1:N; Sf(k)=Sf(k)+st(n)*exp(-i*2*pi*n*dt*k*df); end
Sf(k)=Sf(k)*dt;
End
% Prepare the frequency samples
for k=1:M; f(k)=k*df; end
plot(t,st)); title(‘The Sum of two Sinusoids’); xlabel(‘Time [s]’); ylabel(‘Amplitude’)
figure(2); plot(f,abs(Sf)); title(‘The DFT of the Sum of two Sinusoids’)
xlabel(‘Frequency [Hz]’); ylabel(‘Amplitude’); % Program END
Noise and Frequency Analysis 169
Table 6.3
Short MATLAB FFT Code
%-----------------------------------------------
% Program: FFT.m
%-----------------------------------------------
clear all; close all; clc;
% Get the input parameters f1, f2, dt, T
N=floor(T/dt); df=1/T; fmax=1/(2*dt); a1=1.0; a2=0.5;
for n=1:N
t(n)=(n-1)*dt; f(n)=-fmax+(n-1)*df;
st(n)=a1*sin(2*pi*f1*t(n))+a2*sin(2*pi*f2*t(n));
end
Sf=fftshift(fft(st)*dt); % Apply FFT Swap the left and right halves
plot(t,st); % plot the output in the time domain
title(‘Sum of two Sinusoids’); xlabel(‘Time [s]’); ylabel(‘Amplitude’)
figure(2); plot(f,abs(Sf)); title(‘The DFT of the function’)
xlabel(‘Frequency [Hz]’); ylabel(‘Amplitude’); % Program END
exactly appear at ±fmax. One can run either of these codes and see that the cor-
rect spectrum will be obtained by just reducing the sampling rate, for example,
to ∆t = 5 ms (which means fmax = ±100 Hz).
Any function may be approximated by series summation of sine and co-
sine functions. The number of terms required in the FS representation depends
on the smoothness of the function and the specified accuracy. The degree of
smoothness of the function determines the number of terms in its FS represen-
tation. In addition, only sine or cosine terms contribute the function if it is odd
or even symmetric.
The frequency content of the rectangular pulse should be well analyzed in
order to understand digital communication concepts. A symmetric rectangular
pulse is defined as
T T
1 if − ≤t ≤ +
t 2 2
Re ct = (6.10)
T T
0 if t>
2
occur at both edges of the pulse. These are called overshoot or ringing (i.e., Gibbs
effect in communication, after mathematical physicist Josiah Gibbs, who ex-
plained the phenomenon in 1899). They occur even with a very high number
of terms in the FS representation.
Full FS representation of a rectangular pulse requires an infinite number
of terms. In other words, a rectangular pulse has infinite harmonics.
Digital signals have an infinite number of harmonics and create severe
EMC problems. Because of this, do not upgrade your processor unless it is
essential.
As noted earlier, the FS representation represents the function only in
the given finite period; beyond that, the signal and its FS representation may
totally be different. Figure 6.11 illustrates this. Here, fmax = (x2sin(3x))/ex/2 is
chosen as an example. Its FS representation is given between 0 ≤ x ≤ 10. Both
the function and its 17-term FS representation are plotted for 0 ≤ x ≤ 2 0. As
observed, they match very well in the first period, but are totally different in
the second period.
A few signals widely used in electronic systems and number of significant
FS terms are shown in Table 6.4. The table shows why a triangular pulse is
preferred against square or rectangular pulse in order to mitigate EMC prob-
lems. The last column presents the number of FS terms required for 10% error
between the function and its FS representation.
Question: What is the maximum frequency of a triangular pulse with 50-
µs pulse repetition period assuming 10% and more harmonic contributions?
Noise and Frequency Analysis 171
Table 6.4
Waveforms and Significant Fourier Components
Harmonics
Waveform First Second Third Fourth Fifth Sixth Total (10%)
Sine 1 0 0 0 0 0 1
Square 1.273 0.0 0.424 0.0 0.255 0.0 9
Triangle 0.811 0.0 0.090 0.0 0.032 0.0 3
Pulse (50%) 0.637 0.0 0.212 0.0 0.127 0.0 9
Pulse (25%) 0.450 0.318 0.150 0.0 0.090 0.105 14
Pulse (10%) 0.197 0.187 0.172 0.151 0.127 0.101 26
harmonics of a sinusoidal signal and total waveforms when they add up are
shown in Figure 6.12.
• Power disturbance: Any deviation from the nominal value (or from
some selected thresholds based on load tolerance) of the input AC power
characteristics.
• Power factor, total: The ratio of the total power input in [W], to the
total voltampere input.
• Transient: A subcycle disturbance in the AC waveform that is evidenced
by a sharp, brief discontinuity of the waveform. May be of either polar-
ity and may be additive to, or subtractive from, the nominal waveform.
• Interrupt: The complete loss of voltage for a time period.
• Outage: The complete loss of voltage for at least half a cycle. For exam-
ple, in 50-Hz AC mains, outage means loss of voltage for at least 10 ms.
• Voltage distortion: Any deviation from the nominal sine waveform of
the AC line voltage.
• Voltage impulse: From 50V to 5 kV, voltage pulses with durations be-
tween 0.5 ms and 20 ms.
• Notches: A switching (or other) disturbance of the normal power volt-
age waveform, lasting less than a half cycle, which is initially of oppo-
site polarity than the waveform and is thus subtractive from the normal
waveform in terms of the peak value of the disturbance voltage. This
includes complete loss of voltage for up to a half cycle.
• Harmonic distortion: The mathematical representation of the distortion
of the pure sine waveform.
• Distortion factor: The ratio of the root square value of the harmonic
content to the root square value of the fundamental quantity, expressed
as a percent of the fundamental.
• Sag: An rms reduction in the AC voltage, at the power frequency, for du-
rations from a half cycle to a few seconds. In other words, voltage drops
up to 80% of nominal value for a period longer than a full cycle (20 ms).
• Undervoltage: When used to describe a specific type of long duration
variation, refers to an rms decrease in the AC voltage, at the power fre-
quency, for a period of time greater than 1 min.
• Voltage regulation: The degree of control or stability of the rms voltage
at the load. Often specified in relation to other parameters, such as input
voltage changes, load changes, or temperature changes.
• Swell: An increase in rms voltage or current at the power frequency for
durations from a half cycle to 1.0 min. In other words, more than 110%
Noise and Frequency Analysis 175
increase of the voltage or current for a period longer than a full cycle
(i.e., longer than 20 ms).
• Frequency change: The change in nominal frequency (i.e., frequency
fluctuations around 50 Hz or 60 Hz AC mains).
• Flicker: A variation of input voltage, either magnitude or frequency,
sufficient in duration to allow visual observation of a change in electric
light source intensity. In other words, it is defined as a periodic voltage
decrease in nearly a 6–7 full cycle (i.e., 8–9 Hz). The human eye is sensi-
tive at 8.2 Hz and suffers from signals around this frequency.
• DC voltage component: The inequality of positive and negative halves
of the AC voltage.
• Common-mode noise (longitudinal): The noise voltage that appears
equally, and in phase, from each current carrying conductor to ground.
• Differential-mode noise (transverse-mode noise): Noise signals mea-
surable between or among active circuit conductors feeding the subject
load, but not between the equipment grounding conductor or associated
signal reference structure and the active circuit conductors.
• Crest factor (of a periodic function): The ratio of the peak value of a
periodic function to the rms value.
• Ground: (a) A conducting connection, whether intentional or acciden-
tal, by which an electric circuit or equipment is connected to Earth or
to some conducting body of relatively large extent that serves in place of
the Earth; (b) high-frequency reference or approximately that potential
on conductors connected to it and for conducting ground currents to
and from Earth (or the conducting body).
frequency shift, glitch, power surge, shared ground, spike, and subcycle outages,
have different histories of usage, and some may have specific definitions for
other applications. Therefore, they are avoided in [9].
V (t ) = V 0 + V1 sin( ωt + Φ1 ) + V 2 sin( ωt + Φ 2 )
(6.11)
+ V 3 sin( ωt + Φ 3 ) + ...... + Vn sin( ωt + Φn )
where
The EU standard for the protection of the AC mains from the connected
equipment is EN61000-3-2 (1998) [11]. This standard gives the limits for har-
monic current emissions for equipment with input current less than 16A per
phase. According to this standard, voltage harmonic distortion [HD]V and total
harmonic distortion [THD]V are defined as:
Vk
∑ Vk2
(6.12)
[HD ]V = , [THD ]V = k =2
V1 V1
I
∑ I k2
(6.13)
[HD ]I = k , [THD ]I = k =2
I1 I1
Noise and Frequency Analysis 177
• Class A, B, C, and D limits are given in the tables for steady state opera-
tions (see Table 6.5).
• Harmonics observed/recorded for a period less than 10s, no matter what
their amplitudes are, are not taken into account.
Table 6.5
Harmonic Current Limits in EN61000-3-2
Class A Class B Class C Class D
Harmonic No [A rms] [A rms] (%) [mA/W]
2 1.08 1.62 2 —
3 2.30 3.45 30 3.4
4 0.43 0.645 — —
5 1.14 1.71 10 1.9
6 0.30 0.45 — --
7 0.77 1.115 7 1.0
8 ≤ n ≤ 40 1.84 / n 2.76 / n — —
9 0.40 0.60 5 0.5
11 0.33 0.495 3 0.35
13 0.21 3.15 3 0.296
15 ≤ n ≤ 39 2.25 / n 3.375 / n 3 3.85 / n
178 A Practical Guide to EMC Engineering
• The installation does not disturb AC mains if [THD]V < 2.5% and
[THD]I < 10%.
• The installation may disturb AC mains if [THD]V ≅ 2.5%–3 % and
[THD]I > 10%. It might be necessary to use harmonic filters.
• The installation does disturb AC mains if [THD]V > 3%. The use of
harmonic filters would be essential technically and financially.
Noise and Frequency Analysis 179
• Current harmonics:
• Increase in power losses, decrease in usable power, low power factor;
• Excess currents in triphase systems;
• Excess heating in transformers and generators;
• Increase in acoustics noise;
• Interference with telephone lines.
• Voltage harmonics:
• Excess heating in transformers and generators;
• Excess heating in capacitors and motors;
• Stress decrease in isolation systems;
• Resonance occurrence or high-voltage breakdown;
• Mechanical oscillations or induction problems in motors.
The best way in harmonic protection is to use active and/or passive filters.
Filtration is a very effective EMC protection approach, not only for harmonics
but also for several other EMC problems.
Signal behaviors determine EMC/EMI problems and protection ap-
proaches. The time-domain EM transients are modeled in terms of Laplace
transforms. Steady-state EM signals are well understood by means of Fourier
transforms. A Fourier transform is a mathematical definition. Electronic sys-
tems are digital; therefore, DFT is used for spectrum analyzes. FT and DFT are
totally different operations; one needs to be very careful in dealing with discrete
signals. DFT introduces periodicity in both the time and frequency domains.
Aliasing and spectrum leakage are two artificial effects of the DFT.
References
[1] Baron Jean Baptiste Fourier, http://bartleby.com/65/fo/Fouriers.html.
[2] Hewlett Packard, The Fundamentals of Signal Analysis, Application Note 243, 1994.
[3] Sevgi, L., “Numerical Fourier Transforms DFT and FFT,” IEEE Antennas and Propagation
Magazine, Vol. 49, No. 3, Jun. 2007, pp. 238–243.
[4] Sevgi, L. “Synthetic Radar Signal Environment: Computer Generation of Signal, Noise,
and Clutter,” IEEE Antennas and Propagation Magazine, Vol. 49, No. 5, Oct. 2007, pp.
192–198.
180 A Practical Guide to EMC Engineering
[5] Sevgi, L., Electromagnetic Modeling and Simulation, New York: IEEE Press–John Wiley
and Sons, April 2014.
[6] Sevgi, L., F. Akleman, L. B. Felsen, “Ground Wave Propagation Modeling: Problem-
Matched Analytical Formulations and Direct Numerical Techniques,” IEEE Antennas and
Propagation Magazine, Vol. 44, No. 1, Feb. 2002, pp. 55–75.
[7] Sevgi, L., Complex Electromagnetic Problems and Numerical Simulation Approaches, New
York: IEEE Press–John Wiley and Sons, June 2003.
[8] M. Levy, Parabolic equation methods for electromagnetic wave propagation, IEE, Institution
of Electrical Engineers, 2000
[9] IEEE Std 1100-1999 -Recommended Practice for Powering and Grounding Electronic
Equipment.
[10] IEC 61000-4-30 Ed. 3.0 b(2015) Electromagnetic Compatibility (EMC)-Part 4-30:
Testing and Measurement Techniques -Power Quality Measurement Methods.
[11] EN 61000-3-2 (1998): Electromagnetic Compatibility (EMC)-Part 3: Limits-Section 2:
Limits for Harmonic Current Emissions (Equipment Input Current <=16 A Per Phase).
Bibliography
F. J. Harris, “On the Use of Windows for Harmonic Analysis with the Discrete Fourier Trans-
form,” Proc. IEEE, Vol. 66, No. 1, Jan 1978, pp. 51–83.
Sevgi, L., and Ç. Uluışık, “A Labview-Based Virtual Instrument for Engineering Education: A
Numerical Fourier Transform Tool,” ELEKTRIK, Turkish J. of Electrical Engineering and Comput-
er Sciences (Special issue on Electrical and Computer Engineering Education in the 21st Century:
Issues, Perspectives and Challenges), Vol. 14, No. 1, 2006, pp. 129–152.
7
EMC Test and Measurement
Environments
EMC tests and measurements are performed in two different environments:
inside and outside. The place where outside tests and measurements are per-
formed is called open area test site (OATS). Inside tests and measurements are
performed in screened or shielded rooms and anechoic chambers. Inside or out-
side, the main problem in tests and measurements is the control/elimination of
environmental interference: up to a certain extent, free-space conditions must
be satisfied. Screened rooms isolate inside and outside but do not eliminate
inside reflections. Moreover, tests and measurements may yield values 30–50
dB higher than EUT’s emission because of resonance effects. Anechoic cham-
bers both isolate inside/outside and eliminate internal reflections (emulate free-
space conditions). OATS must be located in an electromagnetically quiet zone.
Whether or not such tests and measurements areas satisfy these requirements is
controlled and traced via calibration.
181
182 A Practical Guide to EMC Engineering
hand, is also an OATS used for determining the free-space AF and can be cali-
brated/validated according to [2].
OATS calibration is performed through site measurements and compari-
sons against a reference site. The ideal or reference site is simply a perfectly con-
ducting, large plane over leveled ground in an electrically quiet area. In [1], the
reference or standard site is defined as “A site comprised of a flat, open-area,
devoid of nearby scatterers such as trees, power lines, and fences that has a large
metallic ground plane.” There are many papers on OATS calibration in the
literature; a few of them may be found in [1–3].
A typical picture of an OATS is sketched in Figure 7.1. The site is an
elliptical open, leveled area with the major and minor axes of 2R and R√3, re-
spectively. The distance between the antenna and EUT is R. The measurement
distance R may be 3m, 10m, or 30m.
The standard site calibration is based on the measurement of the IL be-
tween transmit and receive antennas, on a large, flat, and unobstructed con-
ducting ground plane. As specified in Clause 5.8.2 of CISPR 16-1-4 [1], site
performance may be validated by the NSA method. One antenna is set to be
at a fixed height (e.g., 1m, 2m), while the other one is scanned from 1m to 4m
in height. The maximum response between the two antennas is recorded. The
conducting ground plane is there to ensure the repeatability of the calibration
measurements.
The critical parameter in OATS calibration is NSA.
volume and the receive antenna is outside this test volume at a prescribed
orientation and position. The transmit antenna shall have an approximately
omnidirectional H-plane pattern. Typical receive antennas are hybrid anten-
nas (biconicaI/LPD combination) for 30 to 1000 MHz, or separate biconical
antennas (for 30 MHz to 200 MHz) and LPDA (for 200 MHz to 1000 MHz).
The equipment used during the NSA measurements is listed in Table 7.1. Due
to their large dimensions and phase center problems hybrid antennas are not
recommended for NSA measurements.
As described in [1] the NSA method is used for the calibration of OATS
having interantenna distance greater than 5m. Free-space AFs of the antennas
Table 7.1
Equipment Used During the Calibration
Equipment Firm Model Characteristics
Biconic antenna Schwarzbeck VHA 9103 30–300 MHz
Biconic antenna EMCO 3109 30–300 MHz
Log-periodic antenna Schwarzbeck VUSPL 9111 200–2000 MHz
Bi-log antenna Schaffner CBL 6141A 30–2000 MHz
Spectrum analyzer HP 8560E 30 Hz–2.9 GHz
Coaxial cable – RG213U 18m, 4m, 10m
Signal generator HP 8657B 0.1–2000 MHz
Signal generator Anritsu MG3633A 0.01–2700 MHz
EMC Test and Measurement Environments 185
used in the calibration are needed in this method. For each frequency site, cali-
bration is performed with the following steps:
• Adjust the output level of the signal generator to give a received voltage
display well above ambient and measured receiver or spectrum analyzer
noise.
• Raise the receiving antenna on the mast and scan 1–4m.
• Record the maximum signal level. This value is VSite in (5.45) in Chapter
5.
• Disconnect transmit and receive cables from their antennas. Directly
connect these cables with a straight-through adapter.
• Record the signal level when transmit and receive cables connected. This
value is VDirect in (5.45) of Chapter 5.
• At each frequency and for each polarization record VSite and VDirect .
• Insert transmit and receive AFs at the measurement frequency as shown
in (5.45) of Chapter 5.
• Insert the mutual impedance correction factor ∆AFtot , which applies
only for the specific geometry of horizontal polarization using tunable
dipoles separated by 3m. ∆AFtot = 0 for all other geometries.
• Solve (5.45) of Chapter 5 for AN, which is the NSA for the measure-
ment frequency and polarization used. Then, subtract this AN from the
NSA value.
• If the result is less then ±4 dB, the site is considered validated at that
frequency and polarization.
• Repeat these steps for the next frequency and polarization combination.
Table 7.2
NSA Measurement Table (d=10m, Hor-Pol)
EMC Test and Measurement Environments
187
188 A Practical Guide to EMC Engineering
Figure 7.4 NSA versus frequency and uncertainty limits (Hor-Pol, d=10m).
Figure 7.5 NSA versus frequency and uncertainty limits (Ver-Pol, d=10m).
EMC Test and Measurement Environments 189
Figure 7.6 NSA versus frequency and uncertainty limits (Hor-Pol, d=3m).
Figure 7.7 NSA versus frequency and uncertainty limits (Ver-Pol, d=3m).
190 A Practical Guide to EMC Engineering
choosing proper locations for your source and receiver, you may excite a reso-
nance or you may get rid of that resonance. This is critical in EMC engineering.
The resonance frequencies of a screened room with dimensions a, b, and
c can be calculated from
2 2 2
c m n p
f mnp = + + (7.1)
2 a b c
EMC Test and Measurement Environments 193
Table 7.3 lists a short MATLAB code that calculates the first several reso-
nances of a screened room. A sample output of this code, for a screened room
with a = 4.5m, b = 2.5m, and c = 2.5m, is given in Figure 7.10.
It is possible to locate absorbers (e.g., ferrite tiles) in a screened room and
damp some of the resonances that are critical [14]. The damped screened room
may then be used as a precompliance test facility for both radiated emission and
immunity testing.
Table 7.3
MATLAB Code Calculating Resonances in a Screened Room
%--------------------------------------
% Program: RezDamping.m
%--------------------------------------
clc; clear all; close all;
fmin = 30e6; fmax = 230e6; c=3e8;k=1; a=4.5; b=2.5; d=2.5;
fprintf(‘ \n’);
fprintf(‘---Lowest Order Resonances of the screened room--- \n’);
fprintf(‘ (a = %5.2f [m] b = %5.2f [m] d = %5.2f [m]) \n’,a,b,d);
fprintf(‘---------------------------------------------------- \n’);
fprintf(‘ \n’);
for m=0:4; for n=0:2; for p=1:1
f(k)=(c/2)*sqrt((m/a)^2+(n/b)^2+(p/d)^2); f(k)=f(k)/1e6;
fprintf(‘ m = %2d n = %2d p = %2d f [MHz] = %8.2f \n’,m,n,p,f(k));
k=k+1;
end; end; end % Program END
Figure 7.11 Recorded field strengths before and after resonance damping.
EMC Test and Measurement Environments 195
All four side walls, the ceiling, and the floor are covered with absorbing
materials in a fully anechoic chamber (FAC). As shown in Figures 7.12 and
7.13, tests and measurements are performed for heavy vehicles where the floor
is covered with metal plates. These types of chambers are semi-anechoic (SAC).
CISPR16-1-4 [1] explains all the requirements required to be satisfied by an-
echoic chambers.
An anechoic chamber that operates down to 30 MHz is either very large
with up to 2.5m depth of carbon-loaded absorber over all surfaces or is smaller
but requires greater strength to be lined with ferrite tiles. The ferrite tiles are
expensive, and the room strength required to support the weight of the tiles
adds to the cost. Ferrite tiles on their own are unlikely to give the performance
196 A Practical Guide to EMC Engineering
of an open area test site, especially in small rooms, although for precompliance
work this is not necessarily a problem. Composite absorbers (carbon-loaded
foam with a ferrite backing) fall midway between the extremes of weight and
size and can also be used.
EMC Test and Measurement Environments 197
Figure 7.16 Side view of a TEM cell. The height of the cell (h) is approximately three times
of the height of the EUT. The width (cross section) of the cell is d. The maximum frequency of
the cell depends on d and h.
198 A Practical Guide to EMC Engineering
Figure 7.18 Photos of commercial (a) TEM, and (b) GTEM cells (www.teseq.com).
References
[1] CISPR 16-1-4: 2003 Radio Disturbance and Immunity Measuring Apparatus—Ancillary
Equipment—Radiated Disturbances.
[2] CISPR16-1-5: 2003 Antenna Calibration Test Sites (CALTS) for 30 MHz to 1000 MHz.
[3] Federal Communications Commission, “Calibration of a Radiation Measurement Site-
Site Attenuation,” in Bull. OCE 44. Washington, DC: US Government. Printing Office,
Sept. 1977, Docket 21371.
[4] ANSI/IEEE C63.4-1992, American National Standard Guide for Methods of Measure-
ment of Radio-Noise Emissions from Low-Voltage Electrical and Electronic Equipment
in the Range of 9 kHz to 40 GHz.
[5] ANSI/IEEE C63.5-2006 (Revision of C63.5-2003) American National Standard Guide
for Electromagnetic Compatibility—Radiated Emission Measurements in Electromag-
netic Interference (EMI) Control—Calibration of Antennas (9 kHz to 40 GHz).
[6] Smith, A. A., R. F. German, and J. B. Pate, “Calculation of Site Attenuation from Antenna
Factors,” IEEE Trans on EMC, Vol. 24, No. 3, Aug. 1982, pp. 301–316.
[7] Trakadas, P. T., and C. N. Capsalis, “A Mixed Model for the Determination of Normalized
Site Attenuation in OATS,” IEEE Trans. on EMC, Vol. 43, No. 1, Feb. 2001, pp. 29–36.
[8] Askri, A., C. Vollaire, L. Nicolas, and D. Prebet, “Normalized Site Attenuation Standard
Correction from Numerical Computing,” IEEE Trans. on EMC, Vol. 38, No. 2, Mar.
2002, pp. 693–696.
[9] Sevgi, L., “The Antenna as a Transducer: Simple Circuit and Electromagnetic Models,”
IEEE Antennas and Propagation Magazine, Vol. 49, No. 6, Dec. 2007, pp. 211–218.
200 A Practical Guide to EMC Engineering
[10] Sevgi, L., S. Çakır, and G. Çakır, “Antenna Calibration for EMC Tests and Measurements,”
IEEE Antennas and Propagation Magazine, Vol. 50, No. 3, Jun. 2008, pp. 215–224.
[11] Eser, S., and L. Sevgi, “Open Area Test Site (OATS) Calibration,” IEEE Antennas and
Propagation Magazine, Vol. 52, No. 3, Jun. 2010, pp. 204–212.
[12] ANSI C63.6-1996 (Revision of C63.6-1988), American National Standard Guide for the
Computation of Errors in Open Area Test Site Measurements.
[13] ANSI C63.7-1992 (Revision of C63.7-1988), American National Standard Guide for
Construction of Open Area Test Sites for Performing Radiated Emission Measurements.
[14] Dawson, L., J. F. Dawson, A. C. Marwin, and D. Welsh, “Damping Resonances Within
a Screened Enclosure,” IEEE Trans. on Electromagnetic Compatibility, Vol. 43, No. 1, Mar.
2001, pp. 45–55.
Bibliography
EN 61000-4-20:2010: Electromagnetic Compatibility (EMC)—Testing and Measurement
Techniques. Emission and Immunity Testing in Transverse Electromagnetic (TEM) Waveguides.
8
EMC Test and Measurement Devices
Fundamental measuring equipment used in electrical and electronics engineer-
ing are a multimeter (MM), oscilloscope, spectrum analyzer (SA), and network
analyzer (NA). A MM measures resistance, voltage, and current. It is also used
to check a fault in a circuit, and to find out whether the source supplies energy.
An oscilloscope displays time variations of DC and AC signals. SA and EMI
receivers display signal versus frequency (an EMI receiver is a slightly modified
SA). An NA is used to investigate input/output characteristics (e.g., insertion
loss, reflection loss, or voltage transfer functions of n-port circuits), impedance
versus frequency variation on Smith chart, and so on.
Generally speaking, MMs and EMI receivers are mostly used in EMC
tests and measurements.
201
202 A Practical Guide to EMC Engineering
f
N e = N 0 − 12.6 ln MHz +10Log 10 (B ) [dBW ] (8.1)
3
where, No is equal to the manmade noise level of the measurement site (typi-
cally, –136, –148, or –164 in residential, rural, or remote sites, respectively) and
the term inside the first parenthesis is the noise density given in dBW/Hz. The
noise power of an HF receiver operating at 5 MHz with 1-MHz bandwidth
at residential and remote sites will be −82.4 dBW and −110.4 dBW, respec-
tively (which are nearly 55–70 dB higher than thermal noise levels). They are
dominant at different frequencies. At microwaves, environmental noise can be
neglected, but below, thermal noise is negligible.
Note that both thermal and environmental noise powers depend on sig-
nal bandwidth (i.e., the bandwidth of the measurement). Therefore, one may
expect DMM’s DC voltage measurement noise to be zero (because the fre-
quency of the signal being measured is zero). This is not true. The measurement
204 A Practical Guide to EMC Engineering
Traceability error is the absolute error relative to the standards used. For DC
voltage measurements, the factory traceability error of this product is given to
be 2 ppm [7].
Question 1: What is the relative accuracy if 24-hr operating temperature
is Tcal ± 1°C?
Answer: The accuracy specification in Table 8.1 for this case is ±(0.5 ppm
of reading + 0.05 ppm of range). Using these specifications, relative measure-
ment error is obtained as ±(0.5/1,000,000 × 8) + (0.05/1,000,000 × 10) = ±4.5
µV or 0.45 ppm of 10V.
Question 2: What is the relative accuracy if the measurement is performed
after 24 hours and the operating temperature is 27°C?
Answer: As shown in the table, the optimum technical specifications of
this DMM are based on auto-calibration (ACAL) of the instrument within the
previous 24 hours and for the operating temperatures of Tcal ± 1°C. After 24
hours, the accuracy specification is ±(4.1 ppm of reading + 0.05 ppm of range).
This yields a total relative error of ±(4.1 ppm × 8V) + (0.05 ppm × 10V) =
±33.3 µV. There is no temperature correction for this example.
Table 8.1
Technical Specs of Agilent 3458 DMM
DC Voltage
Temperature Coefficient (ppm
Maximum Input or Reading + ppm of Range)/°C
Range Full Scale Resolution Impedance Without ACAL1 With ACAL2
100 mV 120.00000 10 nV >10 GΩ 1.2 + 1 0.15 + 1
1V 1.20000000 10 nV >10 GΩ 1.2 +0.1 0.15 + 0.1
10 V 12.0000000 100 nV >10 GΩ 0.5 + 0.01 0.15 +0.01
100 V 120.0000000 1 µV 10 MΩ ± 1% 2 + 0.4 0.15 + 0.1
1000 V 1050.00000 10 µV 10 MΩ ± 1% 2 + 0.04 0.15 + 0.01
Accuracy3 [ppm of Reading (ppm of Reading for Option 002) + ppm of Range]
Range 24 Hour4 90 Day5 1 Year5 2 Year5
100 mv 2.5 + 3 5.0 (3.5) + 3 9 (5) + 3 14 (10) + 3
1V 1.5 + 0.3 4.6 (3.1) + 0.3 8 (4) + 0.3 14 (10) + 0.3
10 V 0.5 + 0.05 4.1 (2.6) + 0.05 8 (4) + 0.05 14 (10) + 0.05
100 V 2.5 + 0.3 6.0 (4.5) + 0.3 10 (6) + 0.3 14 (10) + 0.3
1000 V6 2.5 + 0.1 6.0 (4.5) + 0.1 10 (6) + 0.1 14 (10) + 0.1
1Additional error for Tcal or last ACAL±5°C.
2Additional error for Tcal±5°C.
(Tcal: calibration temperature; Acal: Auto calibration in the last 24 hours when the temperature change is less
than ±1°C).
EMC Test and Measurement Devices 207
Question 3: What is the relative accuracy with and without ACAL if the
measurement is performed after 24 hours and the operating temperature is
35°C?
Answer: The total relative accuracy for this case is the same with the pre-
vious example: ±(4.1 ppm × 8V) + (0.05 ppm × 10V) = ±33.3 µV. The tem-
perature coefficient, without ACAL, is equal to ±61.2 µV (±((0.5 ppm × 8V)
+ (0.01 ppm × 10V)) ×12°C). The total relative error is the addition of these
two: (33.3+61.2) = ±94.5 µV. Assuming the same conditions, but using ACAL,
significantly reduce the error due to temperature difference from calibration
temperature. Since the measurement temperature is 7°C beyond the standard
range of Tcal ± 5°C the temperature coefficient is then equal to ±((0.15 ppm
× 8V) + (0.01 ppm × 10V)) ×7°C = ±9.1 µV. Finally, the total relative error is
(33.3+9.1) = ±42.4 µV. If the absolute error is of interest for the last two ex-
amples, the traceability error of (2 ppm × 8) = ±16 µV should be added.
Absolute errors of the past two examples can be calculated by adding the
traceability error (2 ppm × 8) = ±16 µV to the relative errors.
DMM resolution is inversely proportional with the measurement speed.
Speeding up the measurements makes the resolution worse. For example, the
DMM given in [7] supplies 100,000 readings per second with only 4.5-digit
resolution; this is reduced to 6 readings per second if the resolution is raised to
8.5-digit.
Aperture (or integration) time is a period during which the ADC reads the
input signal. Increasing the aperture time potentially improves the accuracy of
the measurement. This is because more samples are internally averaged during a
longer aperture time (noise is sample-to-sample uncorrelated signal; therefore,
simple averaging reduces the noise level and increases SNR). Configuring the
DMM’s aperture time to be an integer multiple of power line cycles also reduces
power line noise. This is called normal mode rejection (NMR). Aperture time
is often expressed in number of power line cycles (NPLCs), where 1 PLC for
60 Hz (50 Hz) is 16.67 ms (20 ms). The noise can be greatly reduced by set-
ting the aperture time to 1 PLC or greater. The larger the value of NPLC, the
greater the line noise reduction, but the longer the measurement period. If frej
and are NS the signal frequency to be rejected and the number of samples to be
averaged, respectively, the aperture time can be set to Ta = NPLC/(NS × frej) in
order to eliminate frej and DC voltage noise. For example, averaging 40 samples
for a NPLC of 3 acquired at a sampling interval of 1.25 ms (or 800 samples/s)
rejects 60 Hz power line noise.
Question 4: What do you know about the accuracy, resolution, and sensi-
tivity of the DMM in Table 8.2?
Answer: The DMM’s resolution is the smallest number called count that
can be displayed. Suppose a DC voltage is measured with the 3V scale of the
DMM described in Table 8.2. Since a DC voltage with 4-digit resolution can
208 A Practical Guide to EMC Engineering
Table 8.2
Technical Specs of a 4-Digit DMM
DC Voltage
Scale: 300 mV, 3V, 30V, 300V, 1000V
Accuracy: ±(0.75% of reading + 3 counts)
Input resistance: Min 10 MΩ
AC Voltage
Scale: 3V, 30V, 300V, 750V
Accuracy: ±(1.2% of reading + 4 count)
Max frequency 500 Hz
DC current
Scale: 300 µA, 3 mA, 30 mA, 300 mA, 10A
Accuracy: ±(1.5% of reading + 5 count)
Resistance
Scale: 300 Ω, 3 kΩ, 30 kΩ, 300 kΩ, 3 MΩ, 30 MΩ
Accuracy: ±(0.75% of reading + 3 count)
be displayed as, for example, 2.543V, or 1.902V, and so on, one count that can
be displayed (i.e., the resolution) is 0.001V. If, on the other hand, a DC voltage
is being measured by the 300V scale, then one count will be 0.1V (the DMM
displays a DC voltage of, for example, 267.3V, or 112.8V). This shows that the
resolution is a relative quantity.
Specifically, for the DMM described here, sensitivity of the device can be
assumed as the resolution under the lowest scale. For example, in this 4-digit
DMM, minimum DC voltage measurement scale is 300mV and one count is
equal to 0.1 mV or 100 µV, therefore sensitivity is said to be 100 µV.
Question 5: What do you expect to see on the display if you want to mea-
sure a 20V DC voltage source having an internal resistance of 5 kΩ, and 150
mV having an internal resistance of 1 kΩ? What are the resolutions and the
absolute and relative errors of these measurements?
Answer: The accuracy of this DMM for different quantity measurements
are specified as in the table. For example, for a DC voltage measurement it is
given to be “±(0.75% of reading + 3 count).”
Error introduced because of loading effect and total measurement error
because of scaling and quantization errors should be calculated for each case. In
general, error introduced because of loading effect is much less than the limited
accuracy of the device and is negligible.
To measure a DC voltage of 20V, the DMM should be best used in 30V
scale. Since the internal and input resistances are given to be 5 kΩ and 10 MΩ,
respectively, the loading error of this measurement will be 20×5/(10005)=0.01V.
The scaling error is calculated to be 20 × 0.0075 = 0.150V. The 4-digit DMM
displays this measurement result as, for example, 19.55. Since, one count is
0.01V, the quantization errors will be 0.03V. As a result, the total error will be
the addition of the loading, scaling, and quantization errors, ±∆Vm = 0.01 +
EMC Test and Measurement Devices 209
0.15 + 0.03 = 0.19V. The value of the measured DC voltage is then in the range
of 19.81V ≤ Vmeas ≤ 20.19V.
To measure a DC voltage of 150 mV, the DMM should be best used in
300 mV scale. In this case the three error values are:
Table 8.3
Basic Technical Specs of the Spectrum Analyzer
Frequency accuracy ±0.2% or 5 MHz + 20% scan resolution
Long-term drift 3 kHz/10 min after the first 30 min
Noise side-bands –75 dBc for 30 × resolution drift from center frequency
Residual FM 1 kHz peak to peak for 2ms
Input sensitivity –115 dBm for 1 kHz bandwidth
3rd order intermodulation (TOIP) 70 dB below 2 full screen signal
Resolution bandwidth –6 dB points from 1 MHz to 1 kHz
212 A Practical Guide to EMC Engineering
• Deflection factor: The ratio of the input signal amplitude to the resul-
tant output indication. The ratio may be in terms of rms volts per divi-
sion, dB per division, watts per division, or any other specified factor.
• Resolution bandwidth: The width in Hz of the SA’s response to a CW
signal. This width is usually defined as the frequency difference at speci-
fied points on the response curve, such as the 3 dB or 6 dB down points.
The manufacturer will specify the dB down points to be used.
• Residual response: A spurious response in the absence of an input, not
including noise and zero pip.
• Residual FM: Short-term displayed frequency instability (jitter) of the
SA caused by instability of the local oscillators (given in terms of PP
frequency deviation).
• Noise sidebands: Undesired response caused by noise internal to the SA
appearing on the display around a desired response.
• Shape factor (skirt selectivity): A measure of the asymptotic shape of the
resolution bandwidth response curve of a SA. Shape factor is defined
as the ratio between bandwidths at two widely spaced points on the
response curve, such as the 3-dB and 60-dB down points.
• Spurious response (spurii, spur): A characteristic of a SA wherein the
displayed frequency does not conform to the input frequency.
Span in this formula is the frequency range between the first and last
frequencies. Span/(sweep points-1) is defined as the horizontal resolution. Al-
though internal reference accuracy is given as “≤ 5% ppm/year (within two
years of adjustment),” it may need to be calculated, for example, from:
214 A Practical Guide to EMC Engineering
(1 × 10 ) × (±1.8 × 10
9 −7
year ) = 180
400 kHz × %0.25 = 1000
3kHz × %5 = 150
2Hz + 0.5 × 400 EkHz (1000 − 1) = 202
Total uncertainty = ±1532 Hz
improved dynamic range. Figure 8.8 shows a photo of Rohde and Schwarz EMI
receiver.
Typical parameters of SAs are first and last frequencies of a sweep, resolu-
tion bandwidth of 3 dB or 6 dB, detectors, sweep time, and video bandwidth.
Although manual control is allowed, resolution and video bandwidths, and
sweep time of SAs are kept constant during a measurement. Typical parameters
of EMI receivers are first and last frequencies of a sweep, RWB of 3 dB or 6 dB,
detectors, dwell time, and frequency step.
General SAs measure either power or voltage; EMI receivers measure and
display only field strengths. EMC standards specify limit values in terms of field
strength at some distance (i.e., 1m, 3m, 10m, 30m) from the EUT. To arrive
at a measurement that is comparable to a specified limit value, one needs to be
sure that the measured value belongs to only the EUT. This is why the AF of the
EMI receive antenna also includes cables and connectors [12].
A typical EMI receiver display is shown in Figure 8.9. Here, emissions
measured in the frequency range of 30 MHz – 1 GHz with RBW = 120 kHz
are displayed as field amplitude, in dBµV/m, versus frequency. A PK detector
is used with 23.67 µs measurement at each frequency. The table at the bottom
of the display lists the PK values marked on the graphics. At the right of the
display, recorded PK, QP, and AV field values of a selected frequency are dis-
played. On the top of the display, PASS/FAIL message according to the related
CISPR is also given. As observed, this EUT fails to comply with the related
CISPR standard.
Note that it is important to be an experienced user in EMC test and
measurements. Measurement devices from simple DMM to EMI receivers have
many parameters that are critical in any test and measurement [3, 12, 13].
The CISPR 16-1-1 standard [14] gives the requirements such as detector
type, resolution bandwidth, selectivity, harmonic suppression, input imped-
ance, dynamic range, noise threshold, and so on, of voltage, current, and EM
field measuring devices in the frequency range of 9 kHz–18 GHz. Commercial
EMI receivers fully comply with CISPR 16-1-X standard family.
Figure 8.10 (a) Ampere’s law, (b) fluke current clamp, (c) Teseq current-injection clamp.
220 A Practical Guide to EMC Engineering
References
[1] Witte, R. A., Electronic Test Instruments: Theory and Applications, Englewood Cliffs, NJ:
Prentice Hall, 1993.
[2] Dunn, P. F., Measurement and Data Analysis for Engineering and Science, New York: Mc-
Graw-Hill, 2005.
[3] Sevgi, L., “Innumeracy: The Meaning of the Numbers We Use,” IEEE Antennas and Prop-
agation Magazine, Vol. 49, No. 2, Apr. 2007, pp.195–190.
[4] National Instruments Tutorial, “Digital Multimeter Measurement Fundamentals,” Feb.
2006, http://zone.ni.com/devzone.
[5] Rec.P–372 “Radio Noise,” ITU (formerly CCIR)-Recs., 2003, or Rec. 322-3, “Character-
istics and Applications of� Atmospheric Radio Noise Data,” CCIR Recs., 1988.
[6] Agilent Technologies Application Note, “System Cabling Errors and DC Voltage Mea-
surement Errors in Digital Multimeters,” AN-1389-1, Jan. 2005, www.agilent.com.
[7] Agilent Technologies 3458 Multimeter Data Sheet, 5965-4971E, 2005.
[8] Fluke Application Note “Understanding Specifications for Precision Multimeters,”
2547797 A-EN-N, Rev A, 2006, www.fluke.com/library.
[9] National Instruments Tutorial, “Digital Multimeter Measurement Cycle,” Feb. 2006,
http://zone.ni.com/devzone.
[10] Cawley, Kevin, “Choosing DMMs and More for High-Performance Applications,”
Keithley Instruments, Inc., No 25260604, Jun. 2004, www.keithley.com.
[11] IEEE Std 748-1979: IEEE Standard for Spectrum Analyzers.
[12] L. Sevgi, S. Çakır, and G. Çakır, “Antenna Calibration for EMC Tests and Measurements,”
IEEE Antennas and Propagation Magazine, Vol. 50, No. 3, Jun. 2008, pp. 215–224.
[13] Sevgi, L., “Digital Multimeters and Basic Measurements,” IEEE Antennas and Propagation
Magazine, Vol. 49, No. 4, Aug. 2007, pp. 232–237.
[14] CISPR 16-1-4: 2003, Radio Disturbance and Immunity Measuring Apparatus—
Measuring Apparatus.
9
EMC Tests and Measurements
EMC requires that systems/equipment must be able to tolerate a specified de-
gree of interference and not generate more than a specified amount of interfer-
ence. These are monitored via EMC tests and measurements. EMC tests and
measurements are performed to certify that a product complies with all the
requirements specified in related standards. A product has to be designed and
fabricated accordingly; it must not cause unintentional, undesired emissions
that may affect or be susceptible to other products in the environment.
The three major components of EMC tests and measurements of EUT
are (1) finding out related standards; (2) knowing operational status during tests
and measurements; and (3) understanding the methods, measurement devices,
and procedures of tests and measurements. EMC tests and measurements must
be reliable, realizable, and repeatable. These are the issues that are controlled via
traceability. Traceability is the possibility of a test or measurement result avail-
able for a comparison within a chain of national/international references within
a specified uncertainty.
Figure 9.1 shows main EMC tests and measurements. As shown in the
figure, emission is measured but susceptibility (immunity) is testes.
Emission measurements are divided into two parts: conducted emissions
(CE) and radiated emissions (RE). CE measurements are performed over EUT
cables up to 30 MHz frequencies, and the voltage values are recorded in dBµV.
A line impedance stabilization network is used during all measurements to sta-
bilize impedances. CE measurements are performed inside a screened room. RE
measurements, performed either in anechoic chamber or on OATS, are usually
requested between the 30-MHz and (up to) 6-GHz frequency range. RE val-
ues are given in terms of electric field in dBµV/m. In addition to CE and RE
measurement, harmonics, flicker, and clicks are also measured. Immunity tests
are also performed conducted and radiated: conducted immunity (CI) below
223
224 A Practical Guide to EMC Engineering
30 MHz and radiated immunity (RE) above 30 MHz. In CI tests, high cur-
rents are injected through cables and EUT’s performance is observed. RI tests
are performed under high homogeneous electric and magnetic fields described
in related standards. Immunity tests also include against ESD, EFT, and surge.
EMC emission and immunity standards are prepared by internat‑ional
expert institutions. These standards specify conditions to be satisfied by equip-
ment that operates in an environment including other equipment. EMC tests
and measurements are complex and time consuming, but not complicated. Ev-
ery detail of the test and measurement procedures are given in related standards.
Table 9.1 lists mandatory EU EMC standards. Some of the fundamental
CISPR standards referred in EU standards are given in Table 9.2.
The CISPR standards given in Table 9.2 are about test and measurement
methods, environments and devices [1–15]. The EU standards in Table 9.1,
on the other hand, present EMC tests and measurements for various products
in our lives [16–29]. For example, the EU standards for conducted emissions
EN55022, EN55011, EN550141, EN61000−6−3, and EN61000−6−4 refer
to CISPR 14 and CISPR 16 as reference standards (see, for example, EMC
standards for IT products in Figure 9.2.)
The starting point of EMC tests and measurements is to identify related
standards.
Table 9.1
EU Standards for Fundamental EMC Tests and Measurements
Conducted EN55022, EN55011, EN55014−1, EN61000−6−3, EN61000−6−4
Emission
EN55022, EN55011, EN61000−6−3, EN61000−6−4, 95/54/EC,
Radiated Emission 72/245/EEC, 75/322/EEC
Power Parasitics EN55014−1
Click EN 55014−1
EN61000−4−3, EN55024, EN61000−6−1, EN61000−6−2,
Radiated Immunity EN55014−2, 95/54/EC, 72/245/EEC, 75/322/EEC
Conducted EN61000−4−6, EN55024, EN55014−2, EN61000−6−1,
Immunity EN61000−6−2
EN61000−4−4, EN55024, EN55014−2, EN61000−6−1,
EFT / Burst EN61000−6−2
EN61000−4−5 , EN55024, EN55014−2, EN61000−6−1,
Surge EN61000−6−2
EN61000−4−2 , EN55024, EN55014−2, EN61000−6−1,
ESD EN61000−6−2
Magnetic Field EN61000−4−8 , EN55024, EN55014−2, EN61000−6−1,
Immunity EN61000−6−2
Voltage outage EN61000−4−11 , EN55024, EN55014−2, EN61000−6−1,
and dip EN61000−6−2
Harmonic Currents EN61000−3−2
Harmonic Immunity EN61000−4−13
Flicker EN61000−3−3
Table 9.2
Basic CISPR EMC Test and Measurement Standards
CISPR 11 EMC thresholds and measurement methods—ISM RF devices
CISPR 14-1 EMC requirements—emission (Household appl, electrical equipment)
CISPR 14-2 EMC requirements—immunity (Household appl, electrical equipment)
CISPR 16-1 EMC requirements 1—emission and immunity measurement setups
CISPR 16-2 EMC requirements 2—emission and immunity measurement methods
CISPR 16-3 EMC requirements 3—reporting and CISPR recommendations
CISPR 16-4 EMC requirements 4—uncertainty and statistical evaluations
CISPR 24 Information technology devices—immunity, thresholds, and test methods
positive events are called opportunities. Risks may come from various sources
including uncertainty, threats at any phase in design, development, production,
legal liabilities, accidents, natural causes and disasters, deliberate attacks, and
226 A Practical Guide to EMC Engineering
so on. Living with those risks is called risk perception and risk management
(see Section 1.2.5). Risk management is the identification, assessment, and pri-
oritization of risks (defined in [30] as the effect of uncertainty on objectives)
followed by coordinated and economical application of resources to minimize,
monitor, and control the probability and/or impact of unfortunate events or to
maximize the realization of opportunities.
Risk is never zero. It is expressed with a real number between zero and one
if calculated using a kind of risk model. If not, risk is expressed as “high,” “low,”
“negligibly low,” “acceptable,” and so on. According to the definition of risk, a
risk is the possibility that an event will occur and adversely affect the achieve-
ment of an objective. Therefore, risk itself has uncertainty. Risk management
aims to control the risk.
An important issue in EMC tests and measurements is the risk analysis.
The degrees of risks in EMC tests and measurements are listed in Table 9.3.
Unless EMC tests and measurements are repeated for every product, there
will always be a risk. The important issue is to minimize this. Depending on
the size of the company and type of the product, it may not be economic to es-
tablish an in-house EMC lab. But it may also not be economic to continuously
send a sample product to a commercial lab. As in all engineering problems,
here, an optimization is necessary.
Table 9.3
Risks in EMC Tests and Measurements
Least risk Establish a complete EMC test center and test every product
Low risk Establish a complete EMC test center and test by sampling
Medium risk (1) Establish a complete EMC test center and test one product
Medium risk (2) Establish a pre-compliance EMC test center and test a product
High risk Do not test—pray!
Table 9.4
Bandwidths and Measurement Times for Emission Measurements
Frequency Band Step # of Meas Meas
[MHz] [kHz] [kHz] steps time (PK) time (QP)
0.15–30 9 5 5970 2 dak 100 dak
30–1000 120 50 19400 6.5 dak 323 dak
EMC Tests and Measurements 229
measurements (i.e., in the frequency range of 150 kHz–30 MHz). On the other
hand, signals having bandwidths more than 120 kHz are assumed broadband
in RE measurements up to 1 MHz. Above this, signals with more than 1 MHz
bandwidth are broadband.
CE measurements are performed between 150 kHz–30 MHz. Mandatory
standards are EN55016-2-1, EN55022, EN55011, and EN55014-1. Results
are given as voltage in dBµV versus frequency. Both the EUT and EMI receiv-
ers are connected via a LISN. AC power supply is also connected to LISN. The
three mission of the LISN are (1) to filter out AC mains and supply the EUT
clean enery; (2) to direct EUT-generated emissions to reach only to the EMI
receiver; (3) to prevent other reflections by impedance matching between AC
mains and the EUT in the frequency range of 150 kHz–30 MHz. Table 9.5 lists
equipment necessary for CE measurements.
To satisfy traceability, CE measurements are performed under the follow-
ing conditions specified in the standards:
• A vertical/horizontal ground plane shall extend at least 0.5m beyond the
projection of the test arrangement, but shall nevertheless have a mini-
mum size of 2m × 2m.
• Bottom and rear of the EUT shall be at a controlled distance of 40 cm
from the ground plane. This ground plane is normally the wall or floor
of a shielded room.
• The EUT is placed on an 80-cm-high nonconducting table. Floor-
standing EUTs are subject to the same provisions as tabletop equipment
with the exception that they should be placed on a floor, with the points
of contact being consistent with normal use.
• Interconnecting cables that hang closer than 40 cm to the ground plane
shall be folded back and forth forming a bundle 40 cm long or less,
Table 9.5
Equipment Necessary for CE Measurements
Equipment Specification
EMI receiver CISPR 16-1-1 compliant
LISN 9 KHz-30 MHz, CISPR 16-1-2
Coaxial cable —
Reference ground plane At least 2m by 2 m.
Shielded enclosure (optional) Depends on the EUT dimension
RF mains filters (optional) Broadband RF filter
Voltage probe (CISPR 16-1-2) 9 KHz–30 MHz, high impedance
230 A Practical Guide to EMC Engineering
hanging approximately in the middle between the ground plane and the
table. The minimum bend radius of the cable should not be exceeded.
• I/O cables that are connected to a peripheral should be bundled in
the center. The end of the cable may be terminated if required using
correct terminating impedance. The total length should not exceed 1m,
if possible.
• If the EUT has only one port, that port is connected to the CDN used
for injection.
Radiated immunity tests are performed with the typical RI test setup given
in EN 61000−4−3 standard (see Figure 9.6). For most of the commercial prod-
ucts, RI tests are mandatory from 80 MHz up to 1 GHz frequency. Table 9.6
lists necessary test equipment. Test levels for commercial products are listed in
Table 9.7 (levels of special and/or military products may be higher sometimes
mandatory by the standards, very often is requested by the producer).
The vertical area where homogeneous high-level electric fields will be gen-
erated is shown in Figure 9.6. This is called calibration plane or uniform field
area (UFA). The calibration plane should be larger than the EUT. Immunity
tests can also be performed by shifting and relocating the EUT if it is larger than
UFA. The critical issue in RI tests is the generation of this field. This should
be verified before every RE test. The generated fields is measured at 16 points
(on a 4 × 4 grid). The four anomalous values are left out, and the field values
on the other 12 points are averaged. According to EN61000-4-3, this average
must be within ± 6 dB of the requested level. If this is satisfied, then RI tests
will be performed.
The following are described in the EN61000-4-3 standard:
• During the test, it should be verified that the EUT has sufficient im-
munity against radiated EM fields. All testing of equipment shall be
Table 9.6
RI Test Equipment
Equipment Specification
Signal generator 80–1000/2700 MHz or higher
Power amplifier Broadband, 80–1000 MHz
Power amplifier Broadband, 1000–2700 MHz
Power meter 80–2700 MHz or higher
Power sensor(s) –60 dBm +20 dBm
Directional coupler(s) 40/60 dB
Function generator (for PM) —
Broadband antenna(s), bilog 80–1000 MHz or higher
Double-ridged/horn antenna 1000–2700 MHz or higher
Anechoic chamber Dimension depends on the EUT.
Isotropic EM probe 80–2700 MHz or higher
RF immunity software —
Coaxial cable(s) —
Computer —
Table 9.7
RI Test Levels Given in EN
6100-4-3 Standard
Electric Field
Test Level Level
1 1 V/m
2 3 V/m
3 10 V/m
4 30 V/m
• The EUT is initially placed with one face coincident with the calibra-
tion plane.
• The EUT face being illuminated shall be contained within the UFA.
• The dwell time of the amplitude modulated carrier at each frequency
shall not be less than the time necessary for the EUT to be exercised and
to respond, but shall in no case be less than 0.5s.
• The test shall normally be performed with the generating antenna facing
each side of the EUT.
• When equipment can be used in different orientations, all sides shall be
exposed to the field during the test.
• The test shall be performed for both horizontal and vertical antenna
polarization for each side of the EUT.
and valid for all AC mains (50−60 Hz 220/380V, 230/400V, and 240/415V).
EN61000-3-2 standard makes harmonic current measurements mandatory up
to 40th harmonics. If the amplitudes of the harmonics monotonically decrease,
then measurements may be stopped at the 19th harmonics.
fluctuations and flicker, limits, and test conditions. Flicker meters are used in
flicker tests.
[8]. ESD immunity tests are performed with ESD simulators/guns manufac-
tured with this model.
ESD pulses (between 4 kV and 8 kV) during ESD immunity tests, as
described in details in EN61000-4-2 standards, are applied directly (through
contacts) and indirectly (over air). The standards describe everything including
how an ESD gun is held, at what distance, how many consecutive pulses will
be applied, and so on. For the sake of repeatability (traceability), ESD immu-
nity tests are performed in a screened room on a conducting/metal plate. The
conducting plate must be, at all sides, 50 cm wider than the size of the EUT.
ESD pulses are applied only to those points and surfaces of the EUT that are
expected to be touched during usual operation, including user access, as speci-
fied in the user manual (e.g., cleaning or adding consumables when the EUT is
powered). The EUT is exposed to at least 200 pulses, 100 each at negative and
positive polarity, at a minimum of four test points. For tabletop equipment,
one of the test points should be the center-front edge of the horizontal coupling
plane, which should be subjected to at least 50 indirect discharges (25 of each
polarity). All other test points should each receive at least 50 direct contact dis-
charges (25 of each polarity). All areas normally touched by the user should be
tested. A photo of a typical ESD immunity test is given in Figure 9.8.
ESD occurs because of potential differences; therefore, the main ESD
protection approach is to ground or bring all elements to the same potential.
ANSI ESD S20.20 and EN100-015/1-1991 standards cover the requirements
for designing, establishing, implementing, and maintaining ESD protection.
Figure 9.9 Parasitic signals (1: useful signal, 2: out-of-band signals, 3: parasitics).
always present, but varies unpredictably in size and direction. They are related
to scatter in the data obtained under fixed conditions, which determines the
repeatability (precision) of the measurement. Random errors (fortunately) fol-
low well-behaved statistical rules. Random error can be reduced by repeating
the measurement as often as possible. Error can be given in one of two ways.
Absolute error is an error that is expressed in physical units. It is the absolute
value of the difference between the measured value and the true value (or the
average value if the true value is not known) of a quantity. Relative error is an
error expressed as a fraction of the absolute error to the true (or average) value
of a quantity. It is always given as a percentage. Uncertainty is a range that is
likely to contain the true value of a quantity being measured or calculated. Un-
certainty can be expressed in absolute or relative terms.
Note that the difference between error and uncertainty [30] should always
be kept in mind. For example, the result of a measurement after correction can
unknowably be very close to the unknown value of the measurand, and thus
have negligible error, even though it may have a large uncertainty.
Error and uncertainty analysis of EMC tests and measurements are done
statistically. Two important parameters in statistics are the average value and
standard deviation. The average value of a measurement performed N times is
given as:
1 x1 + x 2 + x 3 + x 4 + ... + x N
∑ k =1
N
x av = xk = (9.1)
N N
Very often, this is assumed to be the true value of the measurand. The
standard deviation (σ) (i.e., absolute error –∆x) of this measurement repeated
N times will be:
∑ k =1 (x k − x av )
N 2
(9.2)
σ = ∆x =
N −1
C = ( A ± B ) ± ∆c ⇒ ∆c = ∆a + ∆b (9.3)
∆c ∆a ∆b
C = A ×B or C = A / B ⇒ = + (9.4)
C A B
References
[1] IEC/EN61000-1-1, Electromagnetic Compatibility (EMC)—Part 1: General-Section 1:
Application and Interpretation of Fundamental Definitions and Terms.
[2] IEC/EN61000-2-1, Electromagnetic Compatibility (EMC)—Part 2: Environment-Sec-
tion 1: Description of the Environment-Electromagnetic Environment for Low-Frequen-
cy Conducted Disturbances and Signaling in Public Power Supply Systems.
[3] IEC/EN61000-2-3, Electromagnetic compatibility (EMC)—Part 2: Environment-Sec-
tion 3: Description of the Environment-Radiated and Non-Network-Frequency-Related
Conducted Ohenomena.
[4] IEC EN61000-3-2, Electromagnetic Compatibility (EMC)—Part 3-2-Limits—Limita-
tion of Harmonic Current Emissions (equipment input current ≤ 16A per phase).
[5] IEC EN61000-3-3, Electromagnetic Compatibility (EMC)—Part 3-3-Limits—Limita-
tion of Voltage Changes, Voltage Fluctuations and Flicker in Public Low-Voltage Supply
Systems (equipment with rated current ≤ 16A per phase).
[6] IEC/EN61000-3-4, Electromagnetic Compatibility (EMC)—Part 3-4: Limits—Limita-
tion of Emission of Harmonic Currents in Low-Voltage Power Supply Systems (equip-
ment with rated current greater than 16A).
[7] IEC/EN61000-3-5, Electromagnetic Compatibility (EMC)—Part 3: Limits—Section 5:
Limitation of Voltage Fluctuations and Flicker in Low-Voltage Power Supply Systems
(equipment with rated current greater than 16A).
[8] IEC/EN61000-4-2, Electromagnetic Compatibility (EMC)—Part 4-2: Testing and Mea-
surement Techniques-ESD Immunity Test.
[9] IEC/EN61000-4-3, Electromagnetic Compatibility (EMC)—Part 4-3: Testing and Mea-
surement Techniques-Radiated, Radio-Frequency, EM Field Immunity Test.
[10] IEC/EN61000-4-4, Electromagnetic Compatibility (EMC)—Part 4-4: Testing and
Measurement Techniques-Electrical Fast Transient / Burst Immunity Test.
[11] IEC/EN61000-4-5, Electromagnetic Compatibility (EMC)—Part 4-5: Testing and
Measurement Techniques-Surge Immunity Test.
[12] IEC/EN61000-4-6, Electromagnetic Compatibility (EMC)—Part 4-6: Testing and
Measurement Techniques-Immunity to Conducted Disturbances, Induced by Radio-
Frequency Fields.
[13] IEC/EN61000-4-7, Electromagnetic Compatibility (EMC)—Part 4-7: Testing
and Measurement Techniques-General Guide on Harmonics and Inter-Harmonics
Measurements and Instrumentation, for Power Supply Systems and Equipment Connected
Thereto.
EMC Tests and Measurements 245
Bibliography
NIST Technical Note 1297-1997: Guidelines for Evaluating and Expressing the Uncertainty of
NIST Measurement Results.
ITU-Rec SM.328-10 Spectra and Bandwidth of Emissions.
10
EMC and Protection
In order to avoid EMC problems, (1) EMI sources must be identified and sup-
pressed, (2) EMI victims must be strengthened, (3) finally, the coupling means
and mechanisms must be removed. EMC protection should take into account
all of these approaches. Well-known EMC protection approaches are (1) filter-
ing, (2) screening, (3) cabling and grounding, and (4) impedance matching.
247
248 A Practical Guide to EMC Engineering
designed according to the given specifications can then be given via frequency
responses. This may be presented as either IL versus frequency or RL versus
frequency. They are defined in Section 4.2. Serial and parallel LC resonance
circuits are also discussed in Section 4.3. As shown in Figure 4.13, four types
of filters can be realized with only inductors (L) and capacitors (C). While
parallel combinations of L and C elements are enough for LPF and HPF, LC
resonance circuits are required for the realization of BPF and BSF. Any L (C)
element, connected serially between input and output of a two-port circuit, acts
as LPF (HPF). They act oppositely when connected in parallel between input
and output. Similarly, any serial (parallel) LC resonance circuit connected seri-
ally between input and output of a two-port circuit acts as BPF (BSF). They act
oppositely when connected in parallel between input and output.
f
AdB = 10 log 1 + χ 2n ; χ = (10.1)
fc
The user needs only to specify the cut-off frequency, and either the at-
tenuation or the order of the filter. The order of filter n corresponds to the
number of poles of the transfer function and changes with the bandwidth as
well as the attenuation slope. A (four-element) Butterworth prototype LPF is
given in Figure 10.2, which is valid when the source and load impedances are
normalized to 1 Ohm and 3 dB corner frequency to 1/2π Hz. The prototype
lumped element values Ak can be calculated from (A1 → L1 , A2 → C2, A3 →
L3, and so on):
Ak = 2 sin
(2k − 1) π ; k = 1,2,3,...,n (10.2)
2n
Table 10.1 lists prototype element values for the first several filter orders.
Equations (10.1) and (10.2) are enough to design a LPF. Generally, the user
doesn’t have to deal with the order of the filter, unless there is strict space and
cost requirements. All he or she wants is a passband characteristics and the de-
gree of out-of-band attenuation. For example, a statement such as a cut-off fre-
quency of 1000 MHz and 25 dB attenuation at 1500 MHz is enough to design
a LPF. Once the prototype values are known, the real inductor and capacitor
values of the desired filter are then obtained from
RL An An
Ln = H, C n = F (10.3)
2 πf c 2 πf c RL
Table 10.1
Butterworth Prototype LPF Coefficients
n A1 A2 A3 A4 A5 A6 A7 A8 A9 A10
1 2.0000 1.0 — — — — — — — —
2 1.4142 1.4142 1.0 — — — — — — —
3 1.0000 2.0000 1.0000 1.0 — — — — — —
4 0.7654 1.8478 1.8478 0.7654 1.0 — — — — —
5 0.6180 1.6180 2.0000 1.6180 0.6180 1.0 — — — —
6 0.5176 1.4142 1.9318 1.9318 1.4142 0.5176 1.0 — — —
7 0.4450 1.2470 1.8019 2.0000 1.8019 1.2470 0.4450 1.0 — —
8 0.3902 1.1111 1.6629 1.9616 1.9616 1.6629 1.1111 0.3902 1.0 —
9 0.3473 1.0000 1.5324 1.8794 2.0000 1.8794 1.5321 1.0000 0.3473 1.0
Table 10.2
Short MATLAB Code for Butterworth LPF Design
%-------------------------------------------------------
% Butterworth LPF Design
%-------------------------------------------------------
fc = input (‘Cut-off Frequency? ‘); AdB = input (‘Filter Attenuation, AdB? ‘);
f1 = input (‘Attenuation frequency? ‘);
n=ceil(0.5*log10(10^(AdB/10)-1)/log10(f1/fc)) ;
for k=1:n; g(k)=2*sin(((2*k)-1)*pi/(2*n)); end;
if mod(n,2) ~= 0; mm = n-1; else; mm = n; end; nn=200; df=(fmax-fmin)/nn;
for k=1:nn % Frequency loop
f(k) = (fmin+(k-1)*df); omega = 2*pi*f(k);
for kk=1:n
if mod(kk,2) ~= 0
ll=g(kk).*rload./(2*pi*fc); z(kk)=i*omega*ll;
else
cc=g(kk)./(2*pi*fc*rload); z(kk)=-i/(omega*cc);
end; end; end % Program END
EMC and Protection 251
RL 1
Ln = , H, C n = F (10.4)
An 2 πf c An 2 πf c RL
The BPF and BSF design procedures are also straightforward. Figure 10.4
shows the transformation from LPF to BPF and/or BSF. Here, inductors and
capacitors of the LPF are replaced with resonance circuits according to the re-
quirements. For the BPF, serial elements must show low-impedance and parallel
elements high-impedance inside the passband B. This may be achieved by using
serial resonance circuits at serial arms and parallel resonance circuits at parallel
arms of the circuit. For the BSF, serial elements must show high-impedance and
parallel elements low-impedance inside the stopband B. Serial resonance circuit
elements of the real BSF are obtained via
An RL B
Cn = F, Ln = H (10.5)
2 π RL B 2 π f c 2 An
and parallel resonance circuit elements of the real BSF are obtained
B R A
Cn = 2
F, Ln = L n H (10.6)
2 π f c RL An 2 πB
B R A
Cn = 2
F, Ln = L n F (10.7)
2 π f c RL An 2 πB
An RL B
Cn = F, Ln = H (10.8)
2 π RL B 2 π f c 2 An
The components in a filter are classified for a range of climatic and me-
chanical requirements, including lower temperature limit, over temperature
limit, and relative humidity. Characteristics of CM and DM currents are also
important in classification.
A filter should be positioned as close as possible to the connection to the
supply—ideally before the AC mains switch in the front panel and fuses. Oth-
erwise, the connecting cables could provide a coupling path via stray induction
to the unfiltered cables. A good ground connection (with large contact surfaces,
without paint) is also essential.
Z SC = jZ 0 tan θ (10.9)
where θ = βl is the electrical length, and β = 2π/λ is the phase constant (λ is the
wavelength). This impedance becomes
if the line is OC. These two form the basis of LC to TL filter transformation.
Transforming lumped elements to distributed elements is desired espe-
cially when the physical length of lumped element becomes comparable with
the signal wavelength [2, 6]. Richards’ transformation may be used to convert
lumped prototype elements into ideal TLs that have equal commensurate
length. Commensurate lengths can be λ/4 or λ/8, but generally λ/8 is being
used to obtain a more compact filter in size. Using Richards’ transformation,
capacitors can be replaced with OC TLs and inductors with SC TLs. The char-
acteristic impedances of the TLs with SC and OC terminations are given as
RL An R tan θ
Z SC = , ZOC = L (10.11)
tan θ An
Table 10.3
Filter Order Versus Number of Required Kuroda Transformations
Filter Order 2 3 4 5 6 7 8 9 10 11 12
Kuroda 1 2 5 8 13 18 25 32 41 50 61
Transformations
increases. When converting lumped element prototype values into TLs using
Kuroda identities, one needs to separate TL elements via unit elements. Add-
ing a unit element having the same impedance of the scaled filter and the same
commensurate length at both sides of the filter does not change the frequency
response. Kuroda identities are heuristic, and the number of transforms re-
quired in converting lumped element prototype into a distributed form is a
disadvantage [7].
Example: Design a LPF with a 500-MHz cut-off frequency and a 30-dB
attenuation at 700 MHz. The order of the corresponding LC filter is found
from (10.1) at n = 5. The TL equivalent has five stubs. Both LC and TL filters
are given in Figure 10.7.
Tables 10.4 and 10.5 list short MATLAB scripts which plot IL versus
frequency characteristics of LC and TL LPFs, respectively, for another LPF with
100-MHz cut-off frequency and 30-dB attenuation at 200 MHz (see Figure
10.8).
As observed, LC LPF has the exact desired performance. On the oth-
er hand, TL LPF is equivalent to LC LPF in the passband and at the cut-off
Table 10.4
Short MATLAB Butterworth LC LPF Code
%-------------------------------------------------------------------------------------
% Program: LC_LPF.m (fc = 100 MHz; 30dB attenuation at 200 MHz)
%-------------------------------------------------------------------------------------
clc; clear all; tpi=2*pi; Vs = 100; Zs = 50; Zl = 50; % source/ load impedance
f = 0:1e6:250e6; IL = zeros(1,length(f)); % define vectors for freq/IL
% Define low-pass filter component values
L1 = 49.2e-9; C2 = 51.5e-12; L3 = 159.2e-9; C4 = C2; L5 = L1;
for n=1:length(f)
Z1=1i*tpi*f(n)*L1; Z2=-1i/( tpi*f(n)*C2); Z3=1i*tpi*f(n)*L3; Z4=Z2; Z5=Z1;
Ze=Z5+Zl; Zd=Ze*Z4/(Ze+Z4); Zc=Z3+Zd; Zb=Z2*Zc/(Z2+Zc);
Za=Z1+Zb; Zeq=Zs+Za; % Calculate equivalent impedances
I = Vs/Zeq; Ib = I/(1+Zc/Z2); Ie = Ib/(1+Ze/Z4); % Calculate currents
% Calculate voltage at the load with or without the filter present
V2 = Zl*Ie; V2bar = Vs*Zl/(Zl+Zs); IL(n)=20*log10(abs(V2)/abs(V2bar));
end; plot(f/1e6,IL); xlabel(‘Freq [MHz]’); ylabel(‘IL [dB]’); ylim([-50 10]); % END
Table 10.5
Short MATLAB Butterworth TL LPF Code
%------------------------------------------------------------------------------------
% Program: TL_LPF.m (fc = 100 MHz; 30dB attenuation at 200 MHz)
%------------------------------------------------------------------------------------
clc; clear all; Vs = 100; Zs = 50; Zl = 50; % source/load impedance
f = 0:1e6:250e6; IL = zeros(1,length(f)); % define vectors for freq/IL
% Define low-pass filter component values
fc=1e8; lambdac=3e8/fc; A1 = 0.618; A2 = 1.618; A3 = 2.000;
A4 = 1.618; A5 = 0.618; teta=2*pi*fc/3e8*lambdac/8;
Z01 = A1*Zl/tan(teta); Z03 = A3*Zl/tan(teta); Z05 = A5*Zl/tan(teta);
Z02 = Zl*tan(teta)/A2; Z04 = Zl*tan(teta)/A4;
for n=1:length(f)
teta=2*pi*f(n)/3e8*lambdac/8;
Z1 = 1i*Z01*tan(teta); Z3 = 1i*Z03*tan(teta); Z5 = 1i*Z05*tan(teta);
Z2 =-1i*Z02/tan(teta); Z4 =-1i*Z04/tan(teta);
Ze = Z5+Zl; Zd = Ze*Z4/(Ze+Z4); Zc = Z3+Zd;
Zb = Z2*Zc/(Z2+Zc); Za = Z1+Zb; Zeq = Zs+Za;
I = Vs/Zeq; Ib = I/(1+Zc/Z2); Ie = Ib/(1+Ze/Z4); % currents
V2 = Zl*Ie; V2bar = Vs*Zl/(Zl+Zs); % Load voltages
IL(n)=20*log10(abs(V2)/abs(V2bar));
end; plot(f/1e6,IL); xlabel(‘Freq [MHz]’); ylabel(‘IL [dB]’); ylim([-50 10]); % END
Figure 10.8 IL versus frequency of the LPF obtained with the codes in Tables 10.4 and 10.5.
Table 10.6
Short MATLAB Code for IL Calculation of a LC BPF
%---------------------------------------------------------------
% Program : BW_BPF.m
%---------------------------------------------------------------
clc; clear all; rs = 50; rload=rs; vl2=rload/(rload+rs);
fc = input (‘BPF--Cut-off Frequency [MHz] : ‘); fc = fc*1e6;
Band = input (‘Pass Band [MHz] : ‘); Band = Band*1e6;
SBand= input (‘Stop Band (should be more larger Pass Band) [MHz] : ‘);
SBand = SBand*1e6; ff=SBand/Band; AdB = input (‘Out of band atten [dB] : ‘);
display(‘Frequency range for the graphics:’);
fmin = input(‘Min Frequency [kHz] : ‘); fmin=fmin*1e3;
fmax = input(‘Max Frequency [MHz] : ‘); fmax=fmax*1e6;
n=ceil(0.5*log10(10^(AdB/10)-1)/log10(ff)) ; % n=Filter order
display(‘The order of the Filter: ‘); n
%.........Prototype Filter elements..............
for k=1:n; g(k)=2*sin(((2*k)-1)*pi/(2*n)); s(k)=1/g(k); end
display(‘Element values of the LPF prototype filter: ‘); g
if mod(n,2) ~= 0; mm = n-1; else; mm = n; end; nn=200; df=(fmax-fmin)/nn;
for k=1:nn % Frequency loop
f(k) = (fmin+(k-1)*df); omega = 2*pi*f(k);
for kk=1:n
if mod(kk,2) ~= 0
cc=Band/(2*pi*fc*fc*rload*g(kk)); ll=rload*g(kk)/(2*pi*Band);
z(kk)=i*omega*ll-i/(omega*cc);
else
cc=g(kk)/(2*pi*rload*Band); ll=rload*Band/(2*pi*fc*fc*g(kk));
z(kk)=(i*omega*ll)*(-i/(omega*cc))/(i*omega*ll-i/(omega*cc));
end; end
% Checkt the last element and load
if mod(n,2) ~= 0; ZA = rload + z(n); else; ZA = rload; end;
ZLast = ZA; % Calculate total input impedance
for kk=mm:-1:1
if mod(kk,2) == 0; ZT(kk)=z(kk)*ZLast/(z(kk)+ZLast); else;
ZT(kk)=z(kk)+ZLast; end;
ZLast = ZT(kk); end
iload = 1/(rs+ZLast); % Calculate source current
% Distribute and find out load current
for kk=2:2:mm-1; iload=iload*z(kk)/(z(kk)+ZT(kk+1)); end;
iload=iload*z(mm)/(z(mm)+ZA); vl1=iload*rload;
iloss(k)=-20*log10(abs(vl2/vl1));
end % End of frequency loop
plot(f/1e6,iloss); ylim([-40 10]); ylabel(‘IL[dB]’); xlabel(‘Freq [MHz]’) % END
The user only selects the filter type and inputs the two design parameters.
The ML material properties (substrate height of the dielectric material and its
relative permittivity) are also supplied by the user. The rest is handled in the
program automatically, and the LC filter, the TL filter, the ML filter, and all
element values and dimensions are given by the program. The layout of the ML
filter and the frequency responses are also plotted.
The ML equivalent of the filter given in Figure 10.7, obtained with this
package, is shown in Figure 10.10 Both the top view and dimensions are given
EMC and Protection 259
Table 10.7
Short MATLAB Code for IL Calculation of a TL with Stubs
%------------------------------------------------------------------------------------------
% Program : LS_IL_TRLine.m
% Parameters:
% ls(ns) : lengths of ns stubs
% d(ns) : lengths of ns nodes that ls(ns) are connected
% type(ns) : ‘oc’, ‘sc’ for open and short circuit stubs
% Z0 : Characteristic impedance; ZL : Load impedance
% Zin(ns) : Equivalent impedance of the ns node looking towards the load
% Zs(ns) : parallel effects of the ns stubs at the nodes
% | d3 | d2 | d1 |
%---------------\----------\-----------\---------------|
% line Zo \ \ \ ZL
%------------\----\-----\----\------\----\-------------|
% l3 \ \ l2 \ \ l1 \ \
% \____\ \____\ \____\
%-----------------------------------------------------------------------------------------------
clear all; clc; close all; % ZL=75+j*100;
Z0=50; ZL=50; Zsource=50; Vsource=1; VL1=Vsource*ZL/(Zsource+ZL);
fmin = input(‘\n Min Freq [MHz] : ‘); fmax = input(‘\n Max Freq [MHz] : ‘);
fmax=fmax*1e6; fmin=fmin*1e6; nf=200; df=(fmax-fmin)/nf;
ns = input(‘\n Number of stubs [] : ‘);
for k=1:ns
fprintf(‘Give d(%1d) in [cm] \n’,k); d(k) = input(‘\n ‘)*1e-2;
fprintf(‘Give stub length ls(%1d) in [cm] \n’,k);
ls(k) = input(‘\n ‘)*1e-2; tip = input(‘Stub termination [oc] or [sc] \n’, ‘s’);
if tip ==’oc’; type(k)=1; else; type(k)=0; end
end
for n=1:200 % Frequency Loop
freq(n)=fmin+(n-1)*df; beta=2*pi*freq(n)/3e8;
for k=1:ns % Calculate equivalent stub impedances
if type(k) == 1; Zs(k) =-j*Z0/tan(beta*ls(k)); else; Zs(k) = j*Z0*tan(beta*ls(k));
end; end
Ztot=ZL; % Calculate equivalent input impedances at each node
for k=1:ns
Zin(k) = Z0*(Ztot+j*Z0*tan(beta*d(k)))/(Z0+j*Ztot*tan(beta*d(k)));
Ztot = Zin(k)*Zs(k)/(Zin(k)+Zs(k));
end
Isource = Vsource/(Zsource+Ztot); Vleft = Isource*Ztot; Iright = Isource;
for k=ns:-1:1
Ileft = Iright*Zs(k)/(Zs(k)+Zin(k));
Vright = cos(beta*d(k))*Vleft-j*Z0*sin(beta*d(k))*Ileft;
Iright =-j*sin(beta*d(k))*Vleft/Z0+cos(beta*d(k))*Ileft; Vleft = Vright; end
VL2 = Vright; IL(n) = 20*log10(abs(VL2/VL1));
end; plot(freq/1e6,IL); xlabel(‘Freq [MHz]’); ylabel(‘IL [dB]’); % Program END
in this figure. IL versus frequency of all these (LC, TL, and ML) filters is plotted
in Figure 10.11.
Note that frequency responses of LC and TL filters are automatically cal-
culated with this package. On the other hand, the ML structure needs a full
wave simulation package. The GDSII-format file may be used together with the
fullwave FDTD-based package MSTRIP [10] and filter characteristics can be
260 A Practical Guide to EMC Engineering
obtained. The MS curve in Figure 10.10 is obtained in this way. Both MWFil-
terDesigner and MSTRIP tools can be downloaded from www.leventsevgi.net.
10.2 Ferrites
A ferrite bead is a passive electric component used to suppress high-frequency
noise in electronic circuits. It is a specific type of electronic choke. Ferrite beads
employ the mechanism of high dissipation of high-frequency currents in a fer-
rite to build high-frequency noise suppression devices. Ferrite beads may also
EMC and Protection 261
be called ferrite rings, ferrite EMI filters, or ferrite chokes. Ferrite beads are
used as passive bandpass filters. The geometry and EM properties of coiled
wire over the ferrite bead result in a high resistive impedance (resistance) for
high-frequency signals, attenuating high-frequency EMI/RFI electronic noise.
Typical ferrite filters are shown in Figure 10.12.
Ferrite, with high permeability (µ) and very low conductivity (σ), is a
ceramic material, made like pottery or bricks. It contains iron and other fer-
romagnetic elements. If a cable is surrounded by ferrite, then the magnetic
field encircling the cable due to CM current in the cable magnetizes the fer-
rite. Because of its high permeability, the amount of energy stored magnetically
in the ferrite is very high. Thus, the inductance per unit length of cable sur-
rounded by ferrite is very high. A CM choke is an inductor having a high value
of inductance.
The energy is either reflected back up the cable or absorbed resistively
within the ferrite core and dissipated as heat. Only in extreme cases will the
heat be noticeable. A pure inductor does not dissipate energy; it only absorbs
energy from the circuit and returns it at a later time. A ferrite bead, by design,
filters out the high-frequency noise in the circuit by dissipating it as heat. The
ferrite bead is effectively an inductor with a very small Q factor. When electrical
noise within the target frequency range travels in the signal cable, an opposite
electromotor force is induced in the ferrite bead because of its high inductance.
The material used to construct the ferrite bead becomes highly resistive
at the design frequency range, and the induced current inside the bead is dissi-
pated as heat instead of inducing an opposing current back in the signal cable. It
is for this reason that the specific circuit characteristics as well as the frequency
band of the noise need to be taken into account when the ferrite bead is in-
stalled as a noise filter.
The right ferrite filter is selected according to the frequency range where
maximum IL is needed, the amount of attenuation, the cable type and struc-
ture, temperature requirements, and DC currents. It may be mono-block or in
pieces.
In two-port circuits, ferrites are modeled as serial impedances.
10.3 Screening/Shielding
An effective EM protection method is screening or shielding. Shielding is the
degree of isolating two pieces of equipment, devices, or environments [11–14].
This could be satisfied in two ways: (1) EM fields exponentially decay with the
distance, so increasing the distance between an EMI source and victim may sup-
ply required isolation; (2) when this is impossible and/or undesired, then the
isolation is supplied by putting obstacles between the EMI source and victim.
Enclosing devices and/or circuitry in a shielded enclosure is a good way of
controlling EM emissions.
Figure 10.13 shows a scenario for the shielding effectiveness (SE) defini-
tion. SE is a degree of isolation of an enclosure between EM interfering source
and a victim:
E H
SE = 20Log 10 0 = 20Log 10 0 (10.12)
E1 H1
EMC and Protection 263
where E0(H0) and E1(H1) are electric (magnetic) field strength values measured
or calculated at the victim position without and with the enclosure, respectively.
If the enclosure is absent or fully transparent in terms of EM waves, then E1 =
E0 and the SE is 0 dB; otherwise, positive and negative SE corresponds to the
attenuation and amplification, respectively, of the EM interfering signal.
The value of SE depends upon a number of factors: EM interfering source
characteristics, such as wave impedance; frequency; amplitude; polarization;
enclosure characteristics such as thickness, material permittivity, permeability,
conductivity, and so on. Among all, the factor related to the type of EM inter-
fering source is usually dominant.
Choosing the right shielding approach depends on identification of the
EMI source (i.e., the distance and type of disturbing signal). In the far field,
EM waves act as plane waves; therefore, electrical shielding may always be ap-
plied. In the near field, shielding approaches are totally different for electrical
and magnetic disturbing sources. As discussed in Section 3.5, all disturbing
signals act either as an electric dipole or a magnetic dipole. Electric dipole has
strong electric fields in the near field; magnetic fields are dominant in the near
field of magnetic dipole. This is explained in Figure 10.14, where a nonmag-
netic (µ = µ0), mono-block metal sphere with thickness d is used as a shield in
the vicinity of both strong electric and magnetic disturbing signals.
Simple SE expressions are also given in the figure. As observed, SE of a
nonmagnetic metal conductor enclosure is quite different for electric and mag-
netic fields in the near field. For high electric fields, electrical shielding is high at
low frequencies; as the frequency increases, SE decreases, reaching a minimum
at medium frequencies, and then increases again toward high frequencies. For
high magnetic fields, SE is almost zero (the enclosure is almost transparent) at
low and medium frequencies and starts increasing at high frequencies.
264 A Practical Guide to EMC Engineering
Reflection Losses
The SE from reflections can be calculated in terms of wave impedance Zw and
conductor impedance Zs:
E ( Z + Z )2
SE R = 20Log 10 i = 20Log 10 s w
(10.13)
Et 4 Z Z
s w
where Ei and Et are field strengths incident upon and transmitted through the
screen, respectively. The screen impedance is Z s = ωµ σ Ω. For the far field,
the (plane-) wave impedance in free space is Zw = Z0 = 120π Ω. Therefore,
(10.13) reduces to
EMC and Protection 265
1 σr σr
SE R = 20Log 10 = 168 + 10Log 10 (10.14)
4 µr f µr f
But in the near field, the type of interfering source is important. The wave
impedance of electrical (i.e., dipole-like emitting) source is Zwe = Z0(l/2πd) Ω,
while the wave impedance of magnetic (i.e., loop-like emitting) source is Zwm =
Z0(l/2πd) Ω ( (see Section 5.9.1 and Figure 5.31). Therefore, the SER in [dB]
in the near field is
σ r
SE Re = 332 + 10Log 10 3 2 (10.15)
µr f d
for electrical sources, and for the magnetic sources it is given as:
fd σ r
2
SE Rm = 14.6 + 10Log 10 (10.16)
µr
Reflection loss in the near field is different for electric and magnetic
dipoles.
266 A Practical Guide to EMC Engineering
Absorption Losses
EM waves barely penetrate conductors. The degree of this penetration is speci-
fied in terms of skin (penetration) depth δ = 2 ωσµ[m] (see Section 3.4), and
absorption losses are calculated in terms of skin depth as:
t
SE A = 20Log 10 exp = 131.4 t µr σ r f (10.17)
δ
The total SE is then obtained from the addition of SER and SEA. A sim-
ple, MATLAB-based SE calculator code based on (10.12)–(10.17) is given in
Table 10.8. A simple virtual SE calculator tool is also designed for this pur-
pose [15]. The user only needs to choose shielding material type and supply
its thickness and the interfering source distance. The virtual SE tool yields SE
versus frequency plots of both electrical and magnetic interfering sources. After
specifying input parameters and the type of annoying source, one calculates SE
caused by the absorption using (10.17), then checks whether or not the inter-
fering source is in the near or far field region and calculates SE caused by the
reflections accordingly either from (10.14) or (10.15)–(10.16), and finally finds
out the total SE by adding these two contributions.
Two typical examples are given in Figures 10.16 and 10.17 for 0.01-mm-
thick screen from a 30-cm distant interfering source. Figures 10.16 and 10.17
belong to copper (µr = 1, σr = 1.0) and mumetal (µr = 30,000, σr = 0.03) shields,
Table 10.8
Short MATLAB Code for SE Calculations
%--------------------------------------------------------------------
% Program : Shielding Effectiveness (SE) Calculator
%--------------------------------------------------------------------
clear all; clc; format long; sigmar=0.1; mur=500; N=250;
dt=input(‘Material thickness [mm] = ? ‘); d=input(‘\n Screen Distance [cm] = ? ‘);
dt=dt/1e3; d=d/1e2; eps0=10^(-9)/(36*pi); mu0=4*pi*10^(-7); z0=sqrt(mu0/eps0);
fmin=1e3; fmax=1e9; df=(fmax-fmin)/N; df=(fmax-fmin)/N;
for k=1:2*N % Frekans çevrimi
fr(k)=fmin+k*df; wavel=3e8/fr(k);
if d < wavel/(2*pi) % Near Field
SRe=322+10*log10(sigmar/(mur*fr(k)^3*d^2));
SRm=14.6+10*log10(fr(k)*d^2*sigmar/mur);
SA=131.4*dt*sqrt(fr(k)*sigmar*mur);
else % Far Field
SRf=168-10*log10(mur*fr(k)/sigmar); SA=131.4*dt*sqrt(fr(k)*sigmar*mur);
SEf(k)=SRf+SA;
end
SEe(k)=SRe+SA; SEm(k)=SRm+SA;
end
semilogx(fr,SEe,’k’,fr,SEm, ‘r--’); grid; xlabel(‘Freq [MHz])’); ylabel(‘SE [dB]’);
legend(‘Electrical Shielding’,’Magnetic Shielding’); % Program END
EMC and Protection 267
Magnetic Screening/Filtering
Magnetic screening is done via creating magnetic fields opposite to magnetic
field lines of the disturbing signals. Enclosures like Faraday cage do not work
for magnetic fields, since field lines form closed loops. Roughly speaking, mag-
netic screening/filtering is important below 30 MHz. Holes and apertures are
barely important.
Electrical Screening
Electrical screening is applied by blocking field lines that start at positive charg-
es and end at negative charges. Electrical screening is important for frequencies
above 30 MHz. Note that isolated pieces of conductors act as an antenna. Cable
connectors, holes, and apertures reduce electrical screening performance. Con-
tact resistance between different conductors is also important. Several empirical
SE expressions can be used to get a feeling of a degree of screening. One of them
given for an aperture is SEdB = 20log10(λ/2d ), where λ and d are the wavelength
and the longest aperture size. As a rule of thumb, d ≤ λ/20 is necessary to reach
at least 20 dB SE. Experience tells us that in electrical shielding, material type
and connections/bonding are important below 10 MHz; above this frequency,
holes and apertures are critical. Shielding materials can be grouped into three
categories:
right side) enclosure. The solid lines represent two-aperture cases. Some obser-
vations from this figure may be listed as follows:
a voltage on the inside dipole. This is a typical RI problem. Are they the same
for the SE characteristics of the enclosure? Does it matter which dipole radi-
ates? The answer of this question is given in Figure 10.19, where SE versus
frequency of both cases is simulated [15]. Here, a 10-cm cubical enclosure with
2-cm by 6-cm rectangular aperture is used. Two, infinitely thin, 6-cm-long,
center-fed dipoles are taken into account. The inside antenna is located right
across the aperture vertically at the center of the enclosure. The outside antenna
is located vertically 15 cm away from the aperture. Simulations are performed
up energy from a disturbing signal, and they radiate undesired EM energy and
disturb others.
The capacitive or electric coupling results from the interaction of electric
fields between circuits. The inductive, or magnetic, coupling results from the
interaction between the magnetic fields of two circuits. The combination of
electric and magnetic fields is appropriately called EM coupling or radiation.
To prevent radiation of a magnetic field from a conductor grounded at both
ends, the conductor should be shielded, and the shield should be grounded at
both ends.
Common-Impedance Coupling
When a disturbing signal shares a ground connection with a victim, then its
current flowing through the common impedance may develop an unwanted
victim voltage. The common impedance increases with high-frequency signals
or signals with high di/dt. A solution is to separate the connections (no com-
mon current path), and eliminate common impedance. Grounds are the most
usual source of common impedance. The inductance of a straight length of wire
is L = 0.0051/2.3ln [mH], where l and d are length and diameter in inches (a
useful rule of thumb is 20 nH/in).
Magnetic Induction
Even with no direct connection, the fields resulting from circuit operation will
allow coupling between two adjacent circuits. This is known as reactive or near-
field coupling. An AC current flowing in one conductor creates a magnetic
field that may couple with a nearby conductor and induce a voltage in it. The
voltage induced in the victim conductor is V = –M diL/dt, where M is the mu-
tual inductance in [H]. M depends on the areas of the source and victim cur-
rent loops, their orientation and separation, and the presence of any magnetic
screening (note that this is not the same as electric field screening). Typical
values for short lengths of cable loomed together are 0.1 mH to 3mH.
Electrical Induction
Changing voltage on one conductor creates an electric field that may couple
with a nearby conductor and induce a voltage on it. The voltage induced on
the victim conductor is V = –CdVL /dtZin , where C is the coupling capacitance
and Zin is the impedance (to ground) of the victim circuit. The value of C is a
function of the distance between the conductors, their effective areas, and the
presence and point of connection of any electric screening material. Typically,
two parallel wires, 1–3 mm apart, show a coupling capacitance of about 50
pF/m. Note that stray capacitance will complete the coupling path even if the
two circuits are not directly referenced to each other.
274 A Practical Guide to EMC Engineering
Transfer Impedance
The quality of a cable at low frequencies is the transfer impedance and at high
frequencies is screening attenuation. Transfer impedance is defined as the ratio
of the voltage induced on the inner conductor to the current flowing along the
outer conductor because of the external effects (see Figure 10.20).
10.4.1 Cables
Cables in a circuit guide signals from one end to the other. Performances of
cables depend on their screen material and type of the grounding used for their
screen. Coaxial cables are naturally screened. The critical issue for the coaxial
cables is their connection and grounding. Pair or twisted-pair cables must be
screened. Two typical cables are pictured in Figure 10.21. On top, a multicon-
ductor cable with copper screen and ferrite isolator layers is shown. Electrical
screening of this cable is good. Even though ferrite cover is used, magnetic
screening of this cable is poor because of multiconductor. On the other hand,
the cable shown at the bottom is good in terms of both electrical and magnetic
screenings. As much as 100–110 dB electric and magnetic isolation could be
obtained with this cable if used properly in a circuit.
Undesired currents are induced on the cable screen because of environ-
mental EMI. These currents must not couple internal conductors. The perfor-
mance of the cable screen depends on preventing this coupling. A voltage differ-
ence also occurs along the screen of the cable if there is a magnetic field around.
The amount of this voltage difference depends on the unit-length resistance
and inductance of the cable. If the cable is screenless, then this voltage increases
as the frequency increases and may cause significant noise voltages. This occurs
only up to the cable cut-off frequency if the cable is screened. Electrical screen-
ing is effective above cable cut-off frequency. Magnetic screening becomes effec-
tive for frequencies several times higher than the cut-off frequency. Table 10.9
lists a few cables and their cut-off frequencies.
Table 10.9
Cut-Off Frequencies of Some Cables
Cable Type Imp [Ω] fc [kHz] 5fc [kHz] Explanation
RG-6A coaxial 75 0.6 3.0 Double shielded
RG-213 coaxial 50 0.7 3.5
RG-214 coaxial 50 0.7 3.5 Double shielded
assemblies should be used for analog/digital signals, and they should be located
as far apart as possible.
Lightning and pulsed sources may cause high-energy transients into
cables.
10.4.2 Connectors
Connectors are also important in EMC engineering. A tight connection be-
tween two cables and surfaces mitigate EMI. EM signals penetrate materials
like air, water, light, and sound. Leakage points where air, sound, water, or light
leak should be closed. This may not be enough, because EM waves also pen-
etrate solids where air, sound, water, or light cannot leak.
Discontinuities in system/device, such as caps, input panels, ventilation
holes and apertures, and so on, are inevitable. In such cases, connectors, and
bonding materials (i.e., gaskets, seals) are used to prevent EMI. Tight contact is
critical when connecting cables with each other as well as to a panel or surface.
In addition, wide strips are best for grounding conductors. Better/tight contact
for the doors or windows may be achieved by metal doped plastics, gaskets, and
fabrics.
Corrosion is also important in surface mounts. Galvanic corrosion is an
electrochemical process in which one metal corrodes preferentially to another
when both metals are in electrical contact and immersed in an electrolyte. Dis-
similar metals and alloys have different electrode potentials, and when two or
EMC and Protection 277
more come into contact in an electrolyte a galvanic couple is set up. Some met-
als such as gold, silver, copper, and nickel are cathodic; others (e.g., aluminum,
zinc, magnesium, and iron) are anodic. One way in reducing or preventing the
corrosion is to electrically insulate the two metals from each other. This can be
done using plastic or another insulator. Another way is to keep the metals dry
and/or shielded from salts, acids, and so on. Painting, coating, or choosing met-
als with similar potentials or the same metal for all construction also reduces
the corrosion.
Note that BNC connectors developed in Bell Labs back in the 1990s have
50Ω–75Ω characteristic impedances and can be used up to 500–600 MHz
frequencies. SMA connectors extend this range up to 20–25 GHz frequencies.
Power planes can also reduce the loop areas of signal and power traces,
causing a decrease of EMI emissions/susceptibility.
• A ground plane can lower the overall ground impedance and can re-
duce high-frequency ground bounce. Also, the impedance between the
ground and voltage planes is lowered at the high frequencies, and this
reduces power bus ringing.
• The capacitors must be properly selected in their frequency response to
deliver the energy needed. However, as the number of decoupling paths
increase, so do the number of voltage drops across them, and this can
result in power bus transients along with associated CM emissions. This
problem can be minimized with proper power plane design.
• The power plane acts as an effective high-frequency capacitor and, con-
sequently, as an additional energy source needed for cleaner IC outputs.
• The delay for each subsection should be no larger than 1/10th the signal
rise time.
• Multilayer PCBs are better than single-layer PCBS in terms of EMI
mitigation. If a single-layer board must be used, a ground plane should
be utilized to help reduce radiation.
• Top and bottom ground planes can help reduce radiation from multi-
layer boards by at least 10 dB.
• Segmented PC board ground planes are useful for reducing cable radia-
tion due to common mode currents.
280 A Practical Guide to EMC Engineering
1
Z 0 = jX + (10.18)
jB + 1 (RL + jX L )
X y ± RL / Z 0 RL2 + X L2 − Z 0RL
B= (10.19)
RL2 + X L2
1 X LZ0 Z0
X = + − (10.20)
B RL BRL
Equation (10.19) yields two solutions for B and X, which means both B
and X may be inductive (positive) or capacitive (negative).
The circuit in Figure 10.25(b) (i.e., for RL < Z0) may be solved in a similar
approach and the following equations may be obtained:
1 1
= jB +
Z0 RL + j ( X + X L ) (10.21)
X = ± RL ( Z 0 − RL ) − X L , (10.22)
( Z 0 − RL ) / RL
B=± (10.23)
Z0
Table 10.10
Short MATLAB Code for a Single-Stub Impedance Matching
%------------------------------------------------------------------------------------
% Program: ImpedanceMatch.m; Usage: ImpedanceMatch(50,15+10j)
% Result = [distance, length] (in wavelengths)
%-------------------------------------------------------------------------------------
function stub = ImpedanceMatch(Z0,ZL)
if ZL == Z0
disp(‘Load is already matched.’);
else
% Voltage reflection coefficient at the LOAD
GammaL = (ZL-Z0)/(ZL+Z0); thL = angle(GammaL);
bl = thL/2 + [1;-1]*acos(-abs(GammaL))/2;
bd = acot(tan(2*bl-thL)/2);
end
position = bl/2/pi; stublength = bd/2/pi; disp(‘ ‘);
disp(‘Distance and Stub length in wavelengths:’);
stub = mod([position,stublength], 0.5);
end %------------------PROGRAM END------------------------
The Smith chart is a graphic showing the relation between voltage reflec-
tion coefficient (Γ) and normalized TL impedance (z = Zy /Z0 = r ± x). It is
plotted on the complex reflection coefficient plane in 2D and scaled with the
normalized impedance.
Polar coordinates on the chart are the amplitude and phase angle of the re-
flection coefficient. The Cartesian coordinates are the real and imaginary parts
of the reflection coefficient. The entire chart lies within the unit circle |Γ| = 1.
The chart also shows curves of normalized constant resistances and constant
reactance. Normalized impedance is defined as the ratio of impedance over
characteristic impedance of the line used. The horizontal axis, which divides
the chart into top and bottom semi-circles, is the normalized resistance axis.
It starts from zero at left (SC point) and goes to infinity at right (OC point).
The circles that are tangent to the right boundary are normalized resis-
tance = constant circles. The center of the chart is 1, where ZL = Z0. On the left,
ZL < Z0; on right, ZL > Z0. The circle that is tangent l to both center point and
right boundary is the loci, where real part of the normalized impedance is 1.
The arcs above the resistance axis (initiated from right) are normalized induc-
tance values (i.e., XL /Z0 = ωL/Z0). The arcs below the resistance axis correspond
to normalized capacitances (i.e., YL /Z0 = 1/ωCZ0). Therefore, normalized com-
plex impedance with (positive) negative imaginary parts are located at the top
(bottom) semi-circle.
The center of |Γ| = constant circle is the origin. The phase angle of the
reflection coefficient may be directly read on the chart measured from the right
end of the resistance axis up to the line connecting the origin to the normalized
impedance point. For a lossless TL, moving from load toward the source along
the line correspond to turning clockwise on the |Γ| = constant circle, or vice
versa. A complete tour on |Γ| = constant circle corresponds to half wavelength.
That is, impedances at λ/2, λ, 3λ/2, 2λ, and so on, away from the load are also
ZL . The point on |Γ| = constant circle, which is symmetric with respect to the
origin, is the normalized admittance at the same point or normalized imped-
ance at λ/4 away. At the origin, the reflection coefficient is zero (VSWR is equal
to 1). The philosophy of impedance matching is to bring any point on the chart
to the center.
the distance between the normalized impedance point and the center of
the chart. Moving from the load toward the source along the TL cor-
responds to turning on this circle in the clockwise direction.
• Step 3: Mark the normalized admittance y L = a ′ + jb ′, which is just the
symmetric of normalized impedance with respect to the origin.
• Step 4: Turn on the |Γ| = constant circle and go to the intersection point
between |Γ| = constant circle and the unit circle. The real part of nor-
malized impedance is equal to 1 (z′ = 1 + jc). Read the distance in terms
of wavelength between these two points on the outer scale of the chart.
This is the distance d, where series stub will be connected.
• Step 5: Add an open-ended parallel stub at this point, which has a ca-
pacitive effect of jc. Since parallel admittances are directly added, the
total normalized admittance at this point will then be yt = 1.
• Step 6: The stub is open ended (OC) at one end and jc at the other. The
problem is reduced to find out the length of the stub with infinite nor-
malized admittance at one and jc at the other. Read this length on the
chart as the distance from SC (OC impedance means SC admittance) to
jc in the clockwise direction.
References
[1] Bowick, C., Circuit Design, Newness, Boston, 1982.
EMC and Protection 287
Bibliography
Nguyen, C., Analysis Methods for RF, Microwave, and Millimeter-Wave Planar Transmission Line
Structures, New York: John Wiley and Sons, 2000.
Sobol, H., “Applications of Integrated Circuit Technology to Microwave Frequencies,” Proceed-
ings of the IEEE, 1971.
288 A Practical Guide to EMC Engineering
Suggested Readings
The author has been the associate editor and a writer of the “Testing Ourselves”
column of the IEEE Antennas and Propagation Magazine since 2007. The fol-
lowing tutorials have been published, and several free EM virtual tools have
been shared with international EM community. Some of these are very basic,
but some are highly comprehensive. All of them have been used all around the
world in undergraduate- and graduate-level courses as well as in research.
Sevgi, L., “Electromagnetic Modeling and Simulation: Challenges in Validation, Verification and
Calibration,” IEEE Trans on EMC, Vol. 56, Aug. 2014, pp. 750–758.
Uslu, A., and L. Sevgi, “MATLAB-Based Filter Design Program: From Lumped Elements to Mi-
crostriplines,” IEEE Antennas and Propagation Magazine, Vol. 53, No. 1, Feb. 2011, pp. 213–224.
Eser, S., and L. Sevgi, “Open Area Test Site (OATS) Calibration,” IEEE Antennas and Propagation
Magazine, Vol. 52, No. 3, Jun. 2010, pp. 204–212.
Sevgi, L., “Biostatistics and Epidemiology: Hypothetical Tests on Cell Phone Users,” IEEE An-
tennas and Propagation Magazine, Vol. 52, No. 1, Feb. 2010, pp. 267–273.
Sevgi, L., “Electromagnetic Screening and Shielding Effectiveness (SE) Modeling,” IEEE Anten-
nas and Propagation Magazine, Vol. 51, No. 1, Feb. 2009, pp. 211–216.
Sevgi, L., “Electromagnetic Compatibility Engineering Education: Problems, Challenges and
Perspectives,” (special issue) TJEECS, Turkish Journal of Electrical Engineering and Computer Sci-
ences, Vol. 17, No. 3, Dec. 2009, pp. 273–278.
Sevgi, L., “The Antenna as a Transducer: Simple Circuit and Electromagnetic Models,” IEEE
Antennas and Propagation Magazine, Vol. 49, No. 6, Dec. 2007, pp. 211–218.
289
290 A Practical Guide to EMC Engineering
Sevgi, L., “Digital Multimeters and Basic Measurements,” IEEE Antennas and Propagation Maga-
zine, Vol. 49, No. 4, Aug. 2007, pp.232–237.
Sevgi, L., “Numerical Fourier Transforms DFT and FFT,” IEEE Antennas and Propagation Maga-
zine, Vol. 49, No. 3, Jun. 2007, pp.238–243.
Sevgi, L., “Innumeracy: The Meaning of the Numbers We Use,” IEEE Antennas and Propagation
Magazine, Vol. 49, No. 2, Apr. 2007, pp.195–190.
Sevgi, L., and Ç. Uluıık, “A MATLAB-based Transmission Line Virtual Tool: Finite-Difference
Time-Domain Reflectometer,” IEEE Antennas and Propagation Magazine, Vol. 48, No. 1, Feb.
2006, pp. 141–145.
Sevgi, L., “Transmission Line Fault Analysis Using a MATLAB-Based Virtual Time Domain
Reflectometer,” IEEE EMC Newsletter, Fall 2006, pp. 67–72.
Sevgi, L., “EMC and BEM Engineering Education: Physics based Modeling, Hands-on Train-
ing and Challenges,” IEEE Antennas and Propagation Magazine, Vol. 46, No. 2, Apr. 2004, pp.
140–145.
Suggested Books
The followings are several useful EMC books for the EMC engineers who
would like to go beyond the scope of this book.
Suggested Sites
ANSI, American National Standardization Institute, www.ansi.org
CENELEC, European Committee for Electrotechnical Standardization, www.cenelec.eu
CISPR, International Special Committee on Radio Interference, www.iec.ch/cispr
FCC, Federal Communications Commission, www.fcc.gov
EMC Academy, www.emcacademy.org
EMC Information Center, www.compliance-club.com
EMC Banana Skins, www.compliance-club.com/archive/old_archive/Bananaskins.htm
EN Standards, http://ec.europa.eu/growth/single-market/european-standards/
harmonised-standards/
EU EMC Directive, http://ec.europa.eu/growth/single-market/european-standards/
harmonised-standards/electromagnetic-compatibility/
IAF, International Accreditation Forum, www.iaf.nu
ICNIRP, International Commission on Non-Ionizing Radiation Protection, www.icnirp.org
IEC, International Electrotechnical Commission, www.iec.ch
IEC ACEC, Advisory Committee on EMC, www.iec.ch/acec
IEEE, Institute for Electrical and Electronics Engineers, www.ieee.org
ILAC, International Laboratory Accreditation Cooperation, www.ilac.org
292 A Practical Guide to EMC Engineering
293
294 A Practical Guide to EMC Engineering
and teaching electromagnetics via virtual tools. He also teaches popular science
lectures in science, technology, and society.
He is a Fellow of the IEEE, a member of Turkish Chamber of Electrical
Engineers (EMO), an associate editor of the IEEE Antennas and Propagation
magazine’s “Testing Ourselves” column, a member of the IEEE Antennas and
Propagation Society Education Committee, a column writer in the IEEE Re-
gion 8 Newsletter’s “Scientific Literacy” column, on the editorial boards of IEEE
Access, ELEKTRIK, the Turkish Journal of Electrical Engineering and Computer
Sciences, Wiley’s RF and Microwave Computer Aided Engineering, the Journal of
Applied Computational Electromagnetics (ACES), the American Journal of Food
Technology, IJAP, the International Journal of Antennas and Propagation, and
Asian Journal of Management.
He is the author or coauthor of more than 175 journal and 100 inter-
national conference papers. He has published several books in English and
Turkish.
His two books, Complex Electromagnetic Problems and Numerical Simula-
tion Approaches and Electromagnetic Modeling and Simulation, were published
by the IEEE Press (Wiley) in 2003 and 2014, respectively.
Index
Absolute errors, 242–43 RE measurement setup, 196
Absorption losses, 266–68 RI measurement setup, 196
AC circuits, 65 use of, 194
Accreditation, 33–49 Antenna arrays
benefits of, 36 ARRAY package, 140–43
certification versus, 36 broadside, 137, 140, 143
as dynamic process, 36 circular, 144
EMC tests and measurements, 40–44 defined, 137
institutions, 36–38 endfire, 137, 140, 143
interlab comparisons, 44–49 illustrated, 137
introduction to, 33–35 with isotropic radiators, 138–40
lab, 41 phased, 137
MLA coverage, 37 planar, 139
philosophy of, 35 radiation patterns, 137–38
proficiency testing, 44–49 Antenna calibration, 116–17
regional bodies, structure of, 38 AF measurement, 117, 126
TÜRKAK, 38–39 antenna pattern measurement (APM),
Accuracy, proficiency testing and, 47, 48 127–31
AC mains, power quality issues of, 173–76 CALTS, 116, 124, 125
Admittance equivalent capacitance substitution
high frequency equivalent, 71 method (ECSM), 117
parameters, 90 horizontal polarization and, 125
Advisory Committee of Electromagnetic performing, 125–27
Compatibility (ACEC), 20 referencce antenna method (RAM), 116
Aliasing, 166 standard site method (SSM), 116
American wire gauge (AWG) standards, 66 three-antenna calibration (TAC)
Ampere’s law, 52 method, 125, 126
Analog-to-digital converter (ADC), 201–2 uncertainty budget for, 131
Anechoic chamber Antenna effective aperture, 111
composite absorbers, 196 Antenna factor (AF), 112–14
illustrated, 196 defined, 112
operation, 195 effective aperture, 113–14
295
296 A Practical Guide to EMC Engineering
Error effects, 85
absolute, 242–43 function of, 83–86
analysis, 240–44 good, tips for, 85
defined, 241 ground categories, 83–84
propagation, 242–43 lightening effects and, 84
random, 241–42 myths, 86
uncertainty and, 242–43 problem, 84, 85
EU EMC directives, 16–17, 21 wire, 84–85, 86
European Accreditation (EA), 36, 37
Expanded uncertainty, 130, 243 Half-power beamwidth, 106
Extremely low frequency (ELF) radiation, 5 Harmonic distortions, 174, 176
Harmonic generation, 88, 89, 171–72
Faraday’s law, 52 Harmonics
Far end crosstalk (FEXT), 99 analysis, 178
Ferrites, 260–62 current and voltages effects, 179
Filters defined, 169, 171
Butterworth LC, 248–54 effects of, 173
design parameters, 248 measurement of, 236–37
ferrite, 260–62 measurements, 177–78
transmission line, 254–56 Horn antennas, 151–52
types of, 247
See also EMC protection IEEE/ANSI ASC-C63 EMC standards, 23
Flicker Immunity, 3, 31
defined, 175, 237 conducted (CI), 223–24
tests, 237–38 measurements, 27, 43
Fourier coefficients, 161 radiated (RI), 224
Fourier series (FS), 163–71 standards, 224
of a function, 162 Immunity tests, 27, 43
rectangular pulse representation, 165, conducted, 232, 233–34
171 difficulty in performing, 232
summation, MATLAB code for, 163 performing, 236
trapezoidal pulse representation, 164 radiated, 232, 234–36
Fourier transform (FT) types of, 232
defined, 160–61 Impedance
discrete (DFT), 163–71 dipole, 147
for energy signals, 161 high frequency equivalent, 72
numerical computation difficulties, 165 in lumped element model, 71
performing, 166 mismatch, 81
Frequency analysis, 160 parameters, 90
Impedance matching
Gain, 106 bandwidth and, 286
Gain pattern, 106 connectivity and, 286
Galvanic corrosion, 276 importance, 280
Gauss’s law, 52 lumped element circuits, 281–82
Gigahertz transverse electromagnetic parallel stub, 285–86
(GTEM) cells, 197–98 practical realization, 286
Green’s function problem, 55–56, 143 scattered parameter circuits, 282–83
Grounding, 83–86 simplicity and, 286
cable, 275–76 single-stub, MATLAB code, 283
Index 301