Programa 1. Introdução aos circuitos eléctricos 2. Grafos e circuitos resistivos lineares 3. Circuitos dinâmicos lineares 4. Regime forçado sinusoidal 5. Análise no domínio da frequência complexa – Funções de rede H(s): pólos e zeros – Diagramas de Bode de amplitude e de fase – Traçado assimptótico 6. Circuitos resistivos não-lineares
112
Embed
Funções de rede H(s): pólos e zeros Diagramas de Bode de ... · Programa 1. Introdução aos circuitos eléctricos 2. Grafos e circuitos resistivos lineares 3. Circuitos dinâmicos
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Programa
1. Introdução aos circuitos eléctricos2. Grafos e circuitos resistivos lineares3. Circuitos dinâmicos lineares4. Regime forçado sinusoidal5. Análise no domínio da frequência complexa
– Funções de rede H(s): pólos e zeros– Diagramas de Bode de amplitude e de fase– Traçado assimptótico
6. Circuitos resistivos não-lineares
1. Variable-Frequency Response AnalysisNetwork performance as function of frequency.Transfer function
2. Sinusoidal Frequency AnalysisBode plots to display frequency response data
3. Resonant CircuitsThe resonance phenomenon and its characterization
4. ScalingImpedance and frequency scaling
5. Filter NetworksNetworks with frequency selective characteristics:low-pass, high-pass, band-pass
VARIABLE-FREQUENCY NETWORKPERFORMANCE
LEARNING GOALS
°∠== 0RRZRResistor
VARIABLE FREQUENCY-RESPONSE ANALYSIS
In AC steady state analysis the frequency is assumed constant (e.g., 60Hz).Here we consider the frequency as a variable and examine how the performancevaries with the frequency.
Variation in impedance of basic components
°∠== 90LLjZL ωωInductor
Capacitor °−∠== 9011CCj
Zc ωω
Frequency dependent behavior of series RLC network
CjRCjLCj
CjLjRZeq ω
ωωω
ω 1)(1 2 ++=++=
CLCjRC
jj
ωωω )1( 2 −+
=−−
×
CLCRCZeq ω
ωω 222 )1()(|| −+= ⎟⎟
⎠
⎞⎜⎜⎝
⎛ −=∠ −
RCLCZeq ω
ω 1tan2
1
sCsRCLCssZ
sj
eq1)(
notation"in tion Simplifica"2 ++
=
≈ω
For all cases seen, and all cases to be studied, the impedance is of the form
011
1
011
1
......)(
bsbsbsbasasasasZ n
nn
n
mm
mm
++++++++
= −−
−−
sCZsLsZRsZ CLR
1,)(,)( ===
Simplified notation for basic components
Moreover, if the circuit elements (L,R,C, dependent sources) are real then theexpression for any voltage or current will also be a rational function in s
LEARNING EXAMPLE
sL
sC1
R
So VsCsLR
RsV/1
)(++
= SVsRCLCs
sRC12 ++
=
So VRCjLCj
RCjV
js
1)( 2 ++=
≈
ωωω
ω
°∠+××+××
××= −−
−
0101)1053.215()1053.21.0()(
)1053.215(332
3
ωωω
jjjVo
MATLAB can be effectively used to compute frequency response characteristics
USING MATLAB TO COMPUTE MAGNITUDE AND PHASE INFORMATION
011
1
011
1
......)(
bsbsbsbasasasasV n
nn
n
mm
mm
o ++++++++
= −−
−−
),(];,,...,,[
];,,...,,[
011
011
dennumfreqsbbbbden
aaaanum
nn
mm
>>=>>=>>
−
− MATLAB commands required to display magnitudeand phase as function of frequency
NOTE: Instead of comma (,) one can use space toseparate numbers in the array
This sequence will alsowork. Must be careful notto insert blanks elsewhere
GRAPHIC OUTPUT PRODUCED BY MATLAB
Log-logplot
Semi-logplot
LEARNING EXAMPLE A possible stereo amplifier
Desired frequency characteristic(flat between 50Hz and 15KHz)
Postulated amplifier
Log frequency scale
Frequency domain equivalent circuit:
Frequency Analysis of Amplifier
)()(
)()(
sVsV
sVsV
in
o
S
in= )(/1
)( sVsCR
RsV Sinin
inin +
= ]1000[/1
/1)( inoo
oo V
RsCsCsV+
=
⎥⎦
⎤⎢⎣
⎡+⎥
⎦
⎤⎢⎣
⎡+
=ooinin
inin
RsCRsCRsCsG
11]1000[
1)( ⎥⎦
⎤⎢⎣⎡+⎥⎦
⎤⎢⎣⎡+
=π
ππ 000,40
000,40]1000[100 sss
( ) ( )( ) ( ) π
π
000,401001058.79
100101018.3191
1691
≈××=
≈××=−−−
−−−
oo
inin
RC
RC
required
actual
ππππ
000,40000,40]1000[)(000,40||100
sssGs ≈⇒<<<<
Frequency dependent behavior iscaused by reactive elements
)()()(
sVsVsG
S
o=
Voltage Gain
)50( Hz
)20( kHz
NETWORK FUNCTIONS
INPUT OUTPUT TRANSFER FUNCTION SYMBOLVoltage Voltage Voltage Gain Gv(s)Current Voltage Transimpedance Z(s)Current Current Current Gain Gi(s)Voltage Current Transadmittance Y(s)
When voltages and currents are defined at different terminal pairs we define the ratios as Transfer Functions
If voltage and current are defined at the same terminals we defineDriving Point Impedance/Admittance
Some nomenclature
EXAMPLE
⎩⎨⎧
=admittanceTransfer
tanceTransadmit)()()(
1
2
sVsIsYT
gain Voltage)()()(
1
2
sVsVsGv =
To compute the transfer functions one must solve the circuit. Any valid technique is acceptable
LEARNING EXAMPLE
⎩⎨⎧
=admittanceTransfer
tanceTransadmit)()()(
1
2
sVsIsYT
gain Voltage)()()(
1
2
sVsVsGv =
The textbook uses mesh analysis. We willuse Thevenin’s theorem
sLRsC
sZTH ||1)( 1+=1
11RsL
sLRsC +
+=
)()(
1
112
RsLsCRsLLCRssZTH +
++=
)()( 11
sVRsL
sLsVOC +=
−+
)(sVOC
)(sZTH
−
+)(2 sV
2R)(2 sI
=+
=)(
)()(2
2 sZRsVsI
TH
OC
)(
)(
1
112
2
11
RsLsCRsLLCRsR
sVRsL
sL
++++
+
121212
2
)()()(
RCRRLsLCRRsLCssYT ++++
=
)()(
)()()()( 2
1
22
1sYR
sVsIR
sVsVsG T
sv ===
)()(
1
1
RsLsCRsLsC
++
×
POLES AND ZEROS (More nomenclature)
011
1
011
1
......)(
bsbsbsbasasasasH n
nn
n
mm
mm
++++++++
= −−
−− Arbitrary network function
Using the roots, every (monic) polynomial can be expressed as aproduct of first order terms
))...()(())...()(()(
21
21
n
m
pspspszszszsKsH
−−−−−−
=
function network the of polesfunctionnetworktheofzeros
==
n
m
pppzzz
,...,,,...,,
21
21
The network function is uniquely determined by its poles and zerosand its value at some other value of s (to compute the gain)
EXAMPLE
1)0(22,22
,1
21
1
=−−=+−=
−=
Hjpjp
z :poles:zeros
=++−+
+=
)22)(22()1()(
jsjssKsH
841
2 +++
sssK
⇒== 181)0( KH
8418)( 2 ++
+=
ssssH
LEARNING EXTENSIONFind the driving point impedance at )(sVS
)()()(
sIsVsZ S=
)(sI
)(1)()(: sIsC
sIRsVin
inS +=KVL
=+=in
in sCRsZ 1)(
Replace numerical values
Ω⎥⎦⎤
⎢⎣⎡ + M
sπ1001
LEARNING EXTENSION
π)104(
000,20,50 :poles0 :zero
721
1
×=
−=−==
K
HzpHzpz
⎥⎦
⎤⎢⎣
⎡+⎥
⎦
⎤⎢⎣
⎡+
=ooinin
inin
RsCRsCRsCsG
11]1000[
1)( ⎥⎦
⎤⎢⎣⎡+⎥⎦
⎤⎢⎣⎡+
=π
ππ 000,40
000,40]1000[100 sss
For this case the gain was shown to be
))...()(())...()(()(
21
21
n
m
pspspszszszsKsH
−−−−−−
= Zeros = roots of numeratorPoles = roots of denominator
Find the pole and zero locations and value K for the voltage gain G(s) = Vo(s) / Vs(s)
Formas da função de transferência:
1. Forma geral
2. Forma de ganho, pólos e zeros
3. Forma das constantes de tempo
Ganho estático (K0) é o ganho que se observa quando ω=0, i.e. a resposta a sinais DC
Obtém-se a partir de qualquer razão entre tensões e correntes (ex: Vo/Vi, Vo/Ii, Io/Ii, Io/Vi)
zi = zero
pi = pólo
τZi=1/zi constante de tempo associada ao zero zi
τPi=1/pi constante de tempo associada ao pólo pi
SINUSOIDAL FREQUENCY ANALYSIS
)(sH
Circuit represented bynetwork function
⎭⎬⎫
+
+
)cos(0
)(0
θω
θω
tBeA tj
( )⎩⎨⎧
∠++
+
)(cos|)(|)(
0
)(0
ωθωωω θω
jHtjHBejHA tj
)()()(
)()(|)(|)(
ωφωω
ωωφωω
jeMjH
jHjHM
=
∠==
Notation
stics.characteri phase and magnitudecalledgenerally are offunctionas of Plots ωωφω ),(),(M
)(log)(
))(log20PLOTS BODE 10
10 ωωφ
ωvs
(M
⎩⎨⎧
. of function a as function network theanalyzewefrequency theof functionaasnetworkaofbehavior thestudy To
ωω)( jH
HISTORY OF THE DECIBEL
Originated as a measure of relative (radio) power
1
2)2 log10(|PPP dB =1Pover
21
22
21
22)2
22 log10log10(|
II
VVP
RVRIP dB ==⇒== 1Pover
By extension
||log20|||log20|||log20|
10
10
10
GGIIVV
dB
dB
dB
===
Using log scales the frequency characteristics of network functionshave simple asymptotic behavior.
The asymptotes can be used as reasonable and efficient approximations
]...)()(21[)1(]...)()(21[)1()()( 2
233310
bbba
N
jjjjjjjKjH
ωτωτςωτωτωτςωτωω
++++++
=±
General form of a network function showing basic terms
Frequencyindependent
...|)()(21|log20|1|log20
...|)()(21|log20|1|log20
||log20log20
21010
233310110
10010
−++−+−
++++++
±=
bbba jjj
jjj
jNK
ωτωτςωτ
ωτωτςωτ
ω DNDN
BAAB
loglog)log(
loglog)log(
−=
+=|)(|log20|)(| 10 ωω jHjH dB=
212
1
2121
zzzz
zzzz
∠−∠=∠
∠+∠=∠
...)(1
2tantan
...)(1
2tantan
900)(
211
23
3311
1
−−
−−
+−
++
°±=∠
−−
−−
b
bba
NjH
ωτωτςωτ
ωτωτςωτ
ω
Idea: Display each basic termseparately and add the results to
obtain final answer.
Let’s examine each basic term >>
]...)()(21)[1(]...)()(21)[1()( 2
233310
bbba
N
sssssssKsH
ττςτττςτ
++++++
=±
Poles/zerosat the origin
First order terms Quadratic terms for complex conjugate poles/zeros
Constant Term
Poles/Zeros at the origin
⎩⎨⎧
°±=∠×±=
→±
±±
90)()(log20|)(|)( 10
NjNjj N
dBN
N
ωωωω
linestraight a is thislogisaxis-xthe 10ω
0
1a
1bSimple pole or zero ωττ
ω
jsjs
+=+=
11 ⎪⎩
⎪⎨⎧
=+∠+=+
− ωτωτωτωτ
1
210
tan)1()(1log20|1|
jj dB
asymptotefrequency low 0|1| ≈+ dBjωτ
(20dB/dec)asymptotefrequency high ωτωτ 10log20|1| ≈+ dBjfrequency)akcorner/bre1whenmeet asymptotes two The (=ωτ
Behavior in the neighborhood of the corner:
FrequencyAsymptoteCurvedistance to asymptote Argument
corner 0dB 3dB 3 45
octave above 6dB 7db 1 63.4
octave below 0dB 1dB 1 26.6
1=ωτ2=ωτ5.0=ωτ
°≈+∠ 0)1( ωτj
°≈+∠ 90)1( ωτj
⇒<<1ωτ
⇒>>1ωτ
Asymptote for phase
High freq. asymptoteLow freq. Asym.
Simple zero
Simple pole
1b
1b
Quadratic pole or zero ])()(21[ 22 ωτωτς jjt ++= ])()(21[ 2ωτωτς −+= j
Evaluation of frequency response using MATLAB User controlled
>> clear all; close all %clear workspace and close any open figure
>> figure(1) %open one figure window (not STRICTLY necessary)
>> w=logspace(-1,3,200);%define x-axis, [10^{-1} - 10^3], 200pts total
[ ]1004)()5.0(25)( 2 +++
=ωωω
ωωjjj
jjG
>> G=25*j*w./((j*w+0.5).*((j*w).^2+4*j*w+100)); %compute transfer function>> subplot(211) %divide figure in two. This is top part>> semilogx(w,20*log10(abs(G))); %put magnitude here
>> grid %put a grid and give proper title and labels>> ylabel('|G(j\omega)|(dB)'), title('Bode Plot: Magnitude response')
>> semilogx(w,unwrap(angle(G)*180/pi)) %unwrap avoids jumps from +180 to -180>> grid, ylabel('Angle H(j\omega)(\circ)'), xlabel('\omega (rad/s)')>> title('Bode Plot: Phase Response')
Evaluation of frequency response using MATLAB User controlled Continued
Repeat for phase
No xlabel here to avoid clutter
USE TO ZOOM IN A SPECIFIC REGION OF INTEREST
Compare with default!
LEARNING EXTENSION Sketch the magnitude characteristic
DETERMINING THE TRANSFER FUNCTION FROM THE BODE PLOT
This is the inverse problem of determining frequency characteristics. We will use only the composite asymptotes plot of the magnitude to postulatea transfer function. The slopes will provide information on the order
A
A. different from 0dB.There is a constant Ko
B
B. Simple pole at 0.11)11.0/( −+ωj
C
C. Simple zero at 0.5
)15.0/( +ωj
D
D. Simple pole at 3
1)13/( −+ωj
E
E. Simple pole at 20
1)120/( −+ωj
)120/)(13/)(11.0/()15.0/(10)(
++++
=ωωω
ωωjjj
jjG
20|
00
0
1020|dBK
dB KK =⇒=
If the slope is -40dB we assume double real pole. Unless we are given more data
LEARNING EXTENSIONDetermine a transfer function from the composite magnitude asymptotes plot
A
A. Pole at the origin. Crosses 0dB line at 5
ωj5
B
B. Zero at 5
C
C. Pole at 20
D
D. Zero at 50
E
E. Pole at 100
)1100/)(120/()150/)(15/(5)(++
++=
ωωωωωωjjj
jjjG
)1100/)(120/()150/)(15/(5)(
++++
=ssssssG
Resumo blocos diagrama de Bode
2. Pólo na origem:
H(s)=1/s
1. Pólo simples:
H(s)=a/(s+a)3. Pólo complexo
conjugado (assimpt.)
5. Zero na origem:
H(s)= s
4. Zero simples:
H(s)= (s+a)/a6. Zero complexo
conjugado (assimpt.)
Sistemas de 2ª ordem, sobre-amortecimento vs sub-amortecimento / ressonância
s → jω
RESONANT CIRCUITS - SERIES RESONANCE
Im { } 0Z⇒ =
⇒RESONANT FREQUENCY
PHASOR DIAGRAM
QUALITY FACTOR
RESONANT CIRCUITS
These are circuits with very special frequency characteristics…And resonance is a very important physical phenomenon
CjLjRjZ
ωωω 1)( ++=
circuit RLC Series
LjCjGjY
ωωω 1)( ++=
circuitRLCParallel
LCCL 11
0 =⇒= ωω
ω
whenzeroiscircuit eachof reactanceThe
The frequency at which the circuit becomes purely resistive is calledthe resonance frequency
Properties of resonant circuits
At resonance the impedance/admittance is minimal
Current through the serial circuit/voltage across the parallel circuit canbecome very large (if resistance is small)
ξωω
211 :FactorQuality
0
0 ===CRR
LQ
222 )1(||
1)(
CLRZ
CjLjRjZ
ωω
ωωω
−+=
++=
222 )1(||
1)(
LCGY
CjLj
GjY
ωω
ωω
ω
−+=
++=
Given the similarities between series and parallel resonant circuits, we will focus on serial circuits
40dB/dec- of asymptotefrequency High0dB/dec of asymptotefrequency low
rad/s1000 :frequencyer Break/corn
( )2
4/)/( 20
2 ωω
++−=
LRLRLO
( )2
4/)/( 20
2 ωω
++=
LRLRHI
100010 ==
LCω
srad /618=
srad /1618=
decdB /40−
ACTIVE FILTERS
Passive filters have several limitations
1. Cannot generate gains greater than one
2. Loading effect makes them difficult to interconnect
3. Use of inductance makes them difficult to handle
Using operational amplifiers one can design all basic filters, and more, with only resistors and capacitors
The linear models developed for operational amplifiers circuits are valid, in amore general framework, if one replaces the resistors by impedances
Ideal Op-Amp
These currents arezero
Basic Inverting Amplifier
0=+V
+− =⇒ VV gain Infinite
0=−V
0==⇒ +II-impedanceinput Infinite
02
2
1
1 =+ZV
ZV
11
22 V
ZZV −=
Linear circuit equivalent
0=−I
1
2
ZZG −=
1
11 Z
VI =
Basic Non-inverting amplifier
1V
1V
0=+I
1
1
2
10
ZV
ZVV
=−
11
120 V
ZZZV +
=
1
21ZZG +=
01 =I
Basic Non-inverting Amplifier
Due to the internal op-amp circuitry, it haslimitations, e.g., for high frequency and/orlow voltage situations. The OperationalTransductance Amplifier (OTA) performswell in those situations
Operational Transductance Amplifier (OTA)
∞== 0RRin :OTAIdeal
COMPARISON BETWEEN OP-AMPS AND OTAs – PHYSICAL CONSTRUCTION
Comparison of Op-Amp and OTA - Parameters
Amplifier Type Ideal Rin Ideal Ro Ideal Gain Input Current input VoltageOp-Amp 0 0 0
OTA gm 0 nonzero∞∞ ∞
∞
Basic Op-Amp Circuit Basic OTA Circuit
⎥⎦
⎤⎢⎣
⎡+⎥
⎦
⎤⎢⎣
⎡+
==
+=
+=
inS
inv
L
L
S
SinS
inin
invL
L
RRRA
RRR
VvA
VRR
Rv
vARR
Rv
0
0
00
⎥⎦
⎤⎢⎣
⎡+⎥
⎦
⎤⎢⎣
⎡+
==
+=
+=
inS
inm
Linm
SinS
inin
inmL
RRRg
RRR
viG
VRR
Rv
vgRR
Ri
0
00
0
00
∞== vAAAmp-Op Ideal mm gG =
OTAIdeal
Basic OTA Circuits
)0()(10
000
10
vdxxiC
v
vgit
m
+=
=
∫
)0()()( 00
10 vdxxvCgtv
tm += ∫
Integrator
In the frequency domain
10 VCj
gV m
ω=
00
0
=+−=ii
vgi
in
inm polarity)(notice
meq
in
in
gR
iv 1
==
ResistorSimulated
Basic OTA Adder
11vgm
22vgm
21 21vgvg mm +
ResistorSimulated
Equivalent representation
OTA APPLICATION
)(121
30 21
vgvgg
v mmm
+=
mg ofility Programmab
)10
10.,.(
10
7mSgge
gmSg
m
m
m
≥
≤decades 7 - 3 :range
valuesTypical
ABCm IASg ⎥⎦⎤
⎢⎣⎡= 201
Controlling transconductance
LEARNING EXAMPLE
ABCm
m
m
Ig
SmSg
mSg
20
1041044
74
=
×=≥
≤
−
resistoraProduce Ωk25
SSgg m
m
753 10410411025 −− ×>×=⇒=×
meq
in
in
gR
iv 1
==
Resistor Simulated
)(20104 5 AIASS ABC⎥⎦⎤
⎢⎣⎡=× −
AAI ABC μ2102 6 =×= −
LEARNING EXAMPLE Floating simulated resistor
1101 vgi m−= 1202 vgi m=
inmvgi −=0
One grounded terminal
011 ii −= 102 ii =
21 mm gg = operationproper For
ABCm
m
m
Ig
SmSg
mSg
20
1041044
74
=
×=≥
≤
−
resistor10MaProduce Ω
Sgm7
6 101010
1 −=×
= S7104 −×<
The resistor cannot be producedwith this OTA!
LEARNING EXAMPLE
210
210
321
210210
,,,
vvvvvv
ggg mmm
−=+=
b) a)
producetoSelect
ABCm
m
m
Ig
SmSg
mSg
20
1041044
74
=
×=≥
≤
−
)(121
30 21
vgvgg
v mmm
+=
Case a
2;103
2
3
1 ==m
m
m
m
gg
gg
Two equations in three unknowns.Select one transductance
)(102011.0 4
33 AImSg ABCm−×=⇒= Aμ5=
AImSg ABCm μ102.0 22 =⇒=
AImSg ABCm μ501 11 =⇒=
Case b
Reverse polarity of v2!
ANALOG MULTIPLIER
ASSUMES VG IS ZERO
Based on ‘modulating the control current
AUTOMATIC GAIN CONTROL
For simplicity of analysiswe drop the absolute value
IN O IN
IN O
v small v AvAv big vB
⇒ ≈
⇒ ≈
OTA-C CIRCUITS
Circuits created using capacitors, simulated resistors, adders and integrators
integrator
resistor
Frequency domain analysis assumingideal OTAs
1101 im VgI = 0202 VgI m−=
0201 IIIC +=CICj
Vω1
0 =
[ ]021101 VgVgCj
V mim −=ω
1
2
21
0 1 i
m
mm
Vg
Cjg
gV
ω+=
Magnitude Bode plot
1
0
iv V
VG =
2
1
m
mdc g
gA =
Cgf
Cg
mC
mC
2
2
2 =
=
π
ω
LEARNING EXAMPLE
)10(21
4
51
0
πωjV
VGi
v+
== :Desired
ABCm
m
m
Ig
SmSg
mSg
20
101011
63
=
=≥
≤
−
4=dcA kHzfCC 100)10(2 5 =⇒= πω
2
1
m
mdc g
gA =
Cgf
Cg
mC
mC
2
2
2 =
=
π
ω
Two equations in three unknowns.Select the capacitor value
=frequency notch To design, pick one, e.g., C and determine the other
LEARNING BY DESIGN Notch filter to eliminate 60Hz hum
Notch filter characteristic
DESIGN EXAMPLE ANTI ALIASING FILTER FOR MIXED MODE CIRCUITS
Visualization of aliasing
Signals of differentfrequency and the samesamples
Ideally one wants to eliminate frequency components higher than twice the sampling frequency and make sure that all useful frequencies as properly sampled Design specification
Simplifying assumption
Infinite input resistance (no load on RC circuit)
Design equation
15.9R k∴ = Ω
(non-inverting op-amp)DESIGN EXAMPLE “BASS-BOOST” AMPLIFIER