004 Analog Signal Processing en
004 Analog Signal Processing en
Where does analog signal and processing fit into a mechatronic system?
A general transfer system with the input variable x(t) and the output variable y(t) is considered
Problem definition
How can the relationship between input x(t) and output y(t) be described?
In the following, some key functions for the analysis of transmission systems will be considered
=1 x≥0
∞
∫ δ (x) dx = 1
−∞
In the following, a linear time-invariant system (LTI) will be considered as a special case of the general transfer system
Linear System
f (a ⋅ x) = a ⋅ y f (x1 + x2 ) = y1 + y2
Time-invariant system
An event at time t1 is results in the same behavior as an event at t2 if the system is in the same state.
In general, the input function can be represented with the help of the Dirac function
∞
x(t) = ∫ x(τ ) ⋅ δ(t − τ ) dτ
−∞
Transfers functions from the time domain into the image domain
Enables the transformation of time-dependent differential equations into easily solvable algebraic functions
Convolution in the time domain can be represented as a simple multiplication in the image domain
∞
F (s) = L(f (t)) = ∫ f (t) ⋅ e −st
dt s = σ + jω
0
e ⋅ F (s) ds = {
γ+j∞
1 f (t) t ≥ 0
−1
f (t) = L (F (s)) = ⋅∫ st
2πj γ−j∞ 0 t<0
Example - Displacement y(t) for damped spring-mass pendulum with excitation u(t)
1 2D
2
⋅ ÿ (t) + ⋅ ẏ (t) + y(t) = u(t)
ω0 ω0
1 2 2D
2 ⋅ s ⋅ Y (s) + ⋅ s ⋅ Y (s) + Y (s) = U (s)
ω0 ω0
1 s e−sT
Ähnlichkeitssatz f (at) ⋅F ( )
δ(t − T )
a a
1 −sT
Shifting the Image Function e −tT
⋅ f (t) F (s + T ) σ(t − T ) ⋅e
s
2
tn−1 1
df
2nd derivation of the original function = ¨(t) s2 ⋅ F (s) − s ⋅ f (0) − f˙(0) (n − 1)!
sn
f
dt2
n −at 1
d f e
nth Derivation of the original function = f (n)
(t) sn ⋅ F (s) − sn−1 ⋅ f (0) − sn−2 ⋅ f˙(0) − ⋯ − s ⋅ f (n−2) (0) − f (n−1) (0) s+a
dt n
a
d F (s) n sin(at)
nth Derivation of the image function (−1) ⋅ t ⋅ f (t) n n s2 + a2
dsn
a
∂f (t,a) ∂F (s,a) sinh(at)
s2 − a2
s
cos(at)
t
1 s2 + a2
0 s
s
cosh(at)
s2 − a2
∞
f (t)
Integration of the image function
∫ F (ω) dω
t s
−at
a
1−e
s ⋅ (s + a)
a2 a2
∫ ∫
−at 1
t⋅e
(s + a)2
t
Convolution ∫ f1 (τ ) ⋅ f2 (t − τ ) dτ
F1 (s) ⋅ F2 (s)
0 e−at − e−bt 1
=b
a
(s + a) ⋅ (s + b)
b−a
γ+j∞
1
⋅∫
Multiplication of Original Functions f1 (t) ⋅ f2 (t) f1 (w) ⋅ F2 (s − w)dw
2πj γ+j∞
1 1
⋅ e−Dω0 t ⋅ sin(ω0
1− D2 t)
1 − D2 s2 + 2Dω0 s + ω0
ω0
9
Converts functions from the time domain into the frequency domain
Enables the transformation of time-dependent differential equations into easily solvable algebraic functions
∞
1
F (jω) = F(f (t)) = ∫ f (t) ⋅ e −jωt
dt
2π −∞
Unlike the Laplace transformation, the transformed system does not describe a transient system behaviour
The system is viewed in steady state
The transfer function specifies the relationship between the output value Ua and the input value Ue of a component.
The transfer behavior of a component can be described by the complex transfer function
The general transfer function with transient behavior is defined as F (s) with s = σ + jω
Ua
F (s) =
Ue
Ua
F (jω) =
Ue
If the Fourier transform is used instead of the Laplace transform to determine the transfer function, the description of the transient
behavior (e.g. transient response) is lost.
Poles and zeros of the transfer function form an important basis for the analysis
and synthesis of transfer systems.
The representation is done in the pole-zero diagram in the complex number
plane
Example Calculation
2
s
F (s) = 2
s + 2s + 1
Zeros
2
s =0 → sn1,2 = 0 (Double zero)
Poles
2
s + 2s + 1 = 0 → sp1 = −1 ± j
(Conjugate complex pole pair)
The goal, in most cases, is for the output of the assembly to be greater than the input to the assembly. This process is called
amplification.
The input and output quantities of the module are usually voltage signals that are continuous in time.
In most cases, the module should have a linear amplification behaviour. A multiplication of the input variable by the factor n should
therefore also result in a multiplication of the output variable by the factor n.
In order to achieve a signal amplification, active components are used in the amplifier circuits.
Commonly used components are
Bipolar transistors
Field-effect transistors
Operational amplifiers
Output: The output current has no influence on the output voltage, i.e.
you can drive infinitely high loads
Low output resistance
Ra = 0
Ua = A ⋅ (Up − Um ) = A ⋅ Ud
1000 ≤ Ra ≤ 1000Ω
For an ideal operational amplifier the input current is I =0
Very high gain factor
−Ue + Ud = 0
Ua = a ⋅ Ud
This results in
Ua
=a
Ue
For the output value y of a general operational amplifier with amplification factor A the following applies
y = A ⋅ (Up − Um )
If we now define the positive input Up as input value x and a general feedback transfer function f we get
y = A ⋅ (x − f y)
y
Convert to transfer function F = x
y A
F = =
x 1 + Af
y A
F = =
x 1 + Af
A f
This corresponds to the transfer function for an operational amplifier with negative feedback loop.
For very high amplification factors, the transfer function is thus independent of the amplification factor itself.
As long as the differential voltage is not zero, the voltage level Um of the negative input is increased. Asymptotically, the input voltage
difference tends to zero. For a stable system in equilibrium with a negative feedback loop, Ud = 0 can therefore be assumed.
I1 = 0
I2 = 0
−Ue + Ud + Ua = 0
Ud = 0
This results in
Ua
=1
Ue
I1 − I2 = 0
−Ue + I1 ⋅ R1 + Ud = 0
−Ud + I2 ⋅ R2 + Ua = 0
Ud = 0
This results in
Ua R2
=−
Ue R1
−Ue + Ud + R2 ⋅ I = 0
−(R1 + R2 ) ⋅ I + Ua = 0
Ud = 0
This results in
Ua R1 + R2 R1
= =1+
Ue R2 R2
−Ue + I ⋅ R1 + Ud = 0
→ I = (Ue − Ud )/R1
Ud
=0
The operational amplifier switches when Ud =0
The system knows only two stable states depending on the positive
Ue − Ud Ud − Ua Ue −Ua R1
supply voltage UCC and negative supply voltage −UCC of the = → = → Ue = −Ua ⋅
R1 R2 R1 R2 R2
operational amplifier
To determine the threshold values for the two stable states, a case
distinction is necessary
R1
Uth1 = −UCC ⋅
R2
R1
Uth2 = +UCC ⋅
R2
Reset
RUN / Stop
Simulation Speed
Current Speed
Power Brightness
Current Circuit:
−Ue − Ud + I ⋅ R2 = 0
→ I = (Ue + Ud )/R2
Ud
=0 The operational amplifier switches when Ud =0
The system knows only two stable states depending on the positive (Ue + Ud )/R2 = Ua /(R1 + R2 ) → Ue /R2 = Ua /(R1 + R2 )
To determine the threshold values for the two stable states, a case
distinction is necessary
R2
Uth1 = −UCC ⋅
R1 + R2
R2
Uth2 = +UCC ⋅
R1 + R2
Reset
RUN / Stop
Simulation Speed
Current Speed
Power Brightness
Current Circuit:
I1 − I2 = 0
−Ue + Uc + Ud = 0
−Ud + I2 ⋅ R + Ua = 0
dUc dUe
I1 = C ⋅ =C⋅
dt dt
This results in
dUe
Ua = −RC ⋅
dt
I1 − I2 = 0
−Ue + I1 ⋅ R + Ud = 0
−Ud + Uc + Ua = 0
dUc dUa
I2 = C ⋅ = −C ⋅
dt dt
dUa
Ue = −RC ⋅
dt
1
Ua = − ⋅ ∫ Ue dt + Ua (t = 0)
RC
I1 + I2 − I3 = 0
−Ue1 + I1 ⋅ R1 + Ud = 0
−Ue2 + I2 ⋅ R2 + Ud = 0
−Ud + I3 ⋅ R3 + Ua = 0
Ua = − ( ) ⋅ R3
Ue1 Ue2
+
R1 R2
I1 − I3 = 0
−Ue1 + I1 ⋅ R1 + Ud + I2 ⋅ R4 = 0
−Ue2 + I2 ⋅ (R3 + R4 ) = 0
−I2 ⋅ R4 − Ud + I3 ⋅ R2 + Ua = 0
R2
Ua = ⋅ (Ue2 − Ue1 ) mit R1 = R3 und R2 = R4
R1
9 12
Very high input impedance (10 − 10 Ω)
Robust against noise/interference
High common mode rejection ratio
Low common mode gain
High amplification achievable
Gain factor is predefined by integrated resistors
Mostly used as integrated circuit
2R1
Ua = ( 1 + )⋅
R3
⋅ (Ue2 − Ue1 )
RG R2
Task
Elimination of unwanted frequency components in a signal
Frequency-dependent weighting of a signal
Application
Low-pass filtering before analog/digital conversion (aliasing)
Elimination of interfering signals and noise
Filters are circuits that can be used to "isolate" frequency ranges
Implementation options
The magnitude ∣F (jω)∣ of the transfer function F (jω) is The phase response φjω of the transfer function F (jω) is often
often also called amplitude gain or amplitude response also called phase shift or phase difference.
For a general complex transfer function of 1st order the For a general complex transfer function of 1st order the following
following applies applies
a + jb r1 ⋅ e jφ1
r1 j(φ1 −φ2 )
F (jω) = = = ⋅ e = r ⋅ e jφ
F (ω) = = = ⋅ e = r ⋅ e jφ
c + jd r2 ⋅ e 2
jφ r2 c + jd r2 ⋅ e 2
jφ r2
φ1 = arctan ( ) φ2 = arctan ( )
r1 = a2 + b2 r2 =
c 2 + d 2
b d
a c
a + jb ∣a + jb∣ a2 + b2
∣F (jω)∣ = = =
c + jd ∣c + jd∣ 2 2 φ(ω) = φ1 − φ2
+
c d
= arctan ( ) − arctan ( )
b d
a c
ωg = 2π ⋅ fg
and
1
∣F (jωg )∣ = ⋅ ∣F (jωmax )∣
2
Approximation of the amplitude response by applying the tangent to the filter flank and the maximum value of the amplitude
response.
Intersection of the two straight lines at the position of the corner angular frequency
Approximation of the phase response by determining the maximum and minimum values one decade above ω = 10ωg and below
Variant 2
Approximation of the phase response by approximating the tangent at the location of the corner angular frequency ωg
= 20 ⋅ log10 ( )
Ua
∣A(ω)∣dB
Ue
Ua
Ue
A(ω)dB
10 20dB
1 0dB
1/ 2 −3dB
1/2 −6dB
1/10 −20dB
1
∣F (jωg )∣ = ⋅ ∣F (jωmax )∣
2
and the equation for converting the amplitude response into decibels
results
1 1
∣F (jωg )∣dB = 20 ⋅ log ( ⋅ ∣F (jωmax )∣) = 20 ⋅ log ( ) + 20 ⋅ log(∣F (jωmax ∣)
2 2
1
= 20 ⋅ log ( ) + ∣F (jωmax )∣dB = −3dB + ∣F (jωmax )∣dB
2
At the point of the corner angular frequency ωg , the magnitude of the amplitude of the transfer function ∣F (jωg )∣ is approximately
3dB less than the maximum magnitude value of the transfer function ∣F (jωmax )∣
Considering the magnitude of a general transfer function ∣F (jω)∣ of 1st order at the filter slope, the following applies for the
proportionality of the magnitude depending on the angular frequency
1
∣F (jω)∣ ∝ ω or ∣F (jω)∣ ∝
ω
The following applies for the conversion of the absolute value into decibels
From this follows for the magnitude change per decade (ω2 /ω1 = 10)
The filter slope for a 1st order filter is +20dB/decade for the rising filter slope and −20dB for the falling filter slope.
According to the derivation above, slopes of ±40dB/decade result for 2nd order filters or +60dB/decade for 3rd order filters and
so on.
IR = IC
1
−Ue + I ⋅ (R + )=0
jωC
1
−I ⋅ + Ua = 0
jωC
Ua 1
F (jω) = =
Ue 1 + jωRC
1
∣F (jω)∣ =
1+ (ωRC)2
1
∣F (jωg )∣ =
2
1
ωg =
RC
IR = IC
1
−Ue + I ⋅ (R + )=0
jωC
−I ⋅ R + Ua = 0
Ua 1
F (jω) = =
1
Ue 1 + jωRC
1
∣F (jω)∣ =
1
1+ (ωRC)2
1
∣F (jωg )∣ =
2
1
ωg =
RC
1
−Ue + I2 ⋅ (R + )=0 (1)
jωC
1
−Ue + I3 ⋅ (R + )=0 (2)
jωC
1
−Ua + I2 ⋅ − I3 ⋅ R = 0 (3)
jωC
1
Ua jωC −R 1 − jωRC
F (jω) = = =
1 1 + jωRC
Ue +R
jωC
(1)2 + (ωCR)2
∣F (jω)∣ = =1
(1)2 + (ωRC)2
−ωRC
φ(jω) = arctan ( ) − arctan ( ) = −2 ⋅ arctan ( )
ωRC ωRC
1 1 1
1
−Ue + I ⋅ (R + ) + Ua = 0
jωC
−R ⋅ IR + Ua = 0
1
− ⋅ I C + Ua = 0
jωC
I − IR − IC = 0
This results in
Ua jωRC
=
Ue 1 + 3jωRC − (ωRC)2
1
ω0 = ωl ⋅ ωh =
RC
I = IR + IC
Ua
IR =
R
IC = Ua ⋅ jωC
1
−Ue + ( + Ua ⋅ jωC ) ⋅ (R + ) + Ua = 0
Ua
R jωC
Extracting Ua gives
1 1
−Ue + Ua ⋅ [( + jωC ) ⋅ (R + ) + 1] = 0
R jωC
Ua
Converting to U results in
e
Ua 1
=
Ue ( R + jωC ) ⋅ (R + )
1 1
jωC +1
1
= 1
1+ + jωRC + 1 + 1
jωRC
1
= 1
jωRC + 3 +
jωRC
Ua jωRC
=
Ue 1 + 3jωRC − (ωRC)2
−Ue + I ⋅ (jωL + R) + Ua = 0
1
−I ⋅ + Ua = 0
jωC
This results in
Ua 1
=
Ue 1 + jωRC − ω LC
2
1
F1 (jω) = 2
1 + jωRC − ω LC
jωRC
F2 (jω) =
1 + jωRC
The group delay τgr is an important parameter of an LTI-System and describes the delay of the envelope of a narrowband signal.
The group delay τgr can be calculated from the phase response φ(ω) of the transfer function
dφ(ω)
τgr (ω) = −
dω
Butterworth-Filter:
Flat amplitude-frequency response in the passband
Amplitude response drops off sharply
Overshoots with step response (increases with increasing filter order)
Tschebyscheff-Filter:
Steeper falling frequency response
Ripple in the passband of the filter
Bessel (Thomson-Filter):
Optimal behavior of the frequency response (maximum smooth)
Low slope
Constant group delay in the passband
Cauer-Filter
Also known as elliptical filters
Very steep amplitude response between passband and stopband
Oscillating behaviour in the course of the transfer function
Strong phase distortion
To understand the difference between analog signal processing and digital signal processing, a clear distinction between signal and
information is necessary
Signal
In reality, deviations of the discrete values also occur with digital signals. If these deviations are too high, the transmitted information is
lost!
Information
Analog/digital converters are electrical components for converting an analog signal into a digital signal
Since digital signals can only take previously defined discrete levels, values between these discrete levels must be rounded
The values determined from this are also referred to as code words
The number of valid discrete values can be defined for the representation in a binary number system via the number of necessary
bits
The number of bits required to represent all allowed states is also referred to as the resolution
Sample-and-hold circuit enables a defined voltage value to be "held" for a short duration at a specific point in time.
The analog/digital converter needs a certain amount of time to quantize the signal
Sample-and-hold circuitry is therefore necessary as part of the analog/digital conversion process
can be represented
The voltage ULSB (Least Significant Bit) describes the input voltage
distance between two output code words and is related to the FS (Full
Scale) value.
The term ‘Full Scale’ is used to describe the highest possible digital
output code of an ADC with n-bit resolution.
UFSR = URef +
(URef − = 0V)
UFSR = URef + − URef −
UFSR
ULSB = n
2
The Dirac comb Ш(t) represents the sum of a periodic sequence of shifted Dirac-functions
Ш(t) = ∑ δ (t − nT )
T >0
n∈N
With the help of the Dirac comb, the sampling of an arbitrary function y(t) with a sampling rate T can be described. The sampling
process is represented as a multiplication with the Dirac-Comb
Resolution
Smallest quantization of the input voltage range that can be digitally represented by the converter.
Speed (sample rate)
Number of measured values per time unit
Input voltage range, cost
The selection of the required converter depends on the task at hand
See exercises
−q q
e(t) = st <t<
2s 2s
−q/2s 2
s q
eMSE = (e(t)) = ⋅ ∫
2 2
(st) dt =
q q/2s 12
2
q
eRMSE = (e(t)) =
12
The graph shows an example of the quantization error for a sinusoidal input signal
Resolution
The resolution r is defined as the number n of bits used to represent the output code words.
Dynamic Range
The dynamic range DR of an AD-converter is defined as the ratio between the largest and smallest signals that can be precisely
measured and is closely linked to the resolution r and is usually specified in decibels.
2 ⋅ LSB
n
DR = 20 ⋅ log10 ( ) = 6.02 ⋅ n dB
LSB
In addition to the described characteristics such as resolution and dynamic range, there are various characteristics for describing the
dynamic behaviour of AD-converters.
The usual metrics such as SNR, THD, SINAD, ENOB, SFDR for comparing the dynamic behaviour are based on a Fast Fourier
Transformation (FFT) following the AD-converter for spectral analysis of the output.
The following equations show the theoretical relationships for determining the dynamic characteristics. These can be calculated for
the ideal AD-converter and then compared with the actual values of a real AD-converter to obtain a measure of the performance of
the analysed AD-converter.
q⋅2 n
2
q ⋅ 2 n
v(t) = ⋅ sin(ωt) (v(t)) =
2 2 2
Signal/Noise-Ratio
The signal-to-noise ratio (SNR) is a measure of the noise behaviour of an AD-converter and is defined as the ratio of the signal power
PSignal to the noise power PN oise .
2
RMSSignal
=( ) = 20 ⋅ log ( )
PSignal USignal,RMS
SNR =
RMSN oise q/ 12 2
2
∑ Uharmonics,RMS
∑ PHarmonics
THD = =
PSignal USignal,RMS
PSignal
SINAD = SNDR = = 20 ⋅ log ( 10 −SNR/10 + 10 THD/10 )
PN oise + PDistortion
RMSSignal
SFDR = 20 ⋅ log ( )
RMSSpur,Largest
For an optimal AD-converter, the effective number of bits should also correspond to the theoretical resolution of the AD-converter
Noise and distortion reduce the effective number of bits in reality
ENOB specifies the remaining effective number of bits
SINADMeasured − 1.76dB + 20 ⋅
log ( Input Amplitude )
Full Scale Amplitude
SINAD − 1.76dB
ENOB = ENOBFull Scale Input =
6.02 6.02
The digital representation of a signal only allows sampling at discrete time values
The possible incorrect determination of a signal frequency in this case is called aliasing
For a sampling frequency that is too low fs lead to a too low assumed signal frequency
spectrum
n ⋅ fs
n∈N
In order to avoid undersampling and thus the development of the alias effect, the Nyquist-Shannon sampling theorem describes a
minimum sampling frequency that should not be undercut
fs ≥ 2 ⋅ b
The sampling frequency fs should correspond to at least twice the bandwidth b of the signal (Nyquist criterion)
The half of the sampling frequency fs /2 is called Nyquist frequency, which should not be exceeded by the signal
By adhering to the Nyquist-Shannon sampling theorem, an overlap of the frequency spectra is avoided
Function
Parameters
−URef + I ⋅ Re = 0
I ⋅ R + Ua = 0
1 1 1 1 1 16R
= z3 ⋅ + z2 ⋅ + z1 + z0 → Re =
Re 2R 4R 8R 16R z0 + 2z1 + 4z2 + 8z3
n−1
1 1 2 R
n
= ∑ zi ⋅ n−i → Re =
Re 2 R n−1
2 ⋅ zi
i=0
∑i=0
i
URef
Ua = −R ⋅ I = −R ⋅
Re
1 1 1 1
= −URef ⋅ (z3 ⋅ + z2 ⋅ + z1 + z0 )
2 4 8 16
Z
Ua = −URef ⋅
ZM ax + 1
ZM ax = 2 − 1
n
Reconstruction filters are often used to smooth out the abrupt behaviour of the output signal
−τ
t
UA,TP = UA ⋅ (1 − e
Initial situation
Rampe abgetastet mit Sample-and-Hold-Schaltung
The "optimal" setting
UA = UDA ⋅ (1 − e −τ
)
t
τ = RC
1 1
fg,opt = ⋅ fa (Shannon) fg =
2 2πRC
1 1 Ta
= → τopt =
2πτopt 2Ta π
Characteristic Parameters
tH tH
D= = [D ] = %
tL + tH T
1
f=
T
tH
UM = ⋅ UP W M
T
A digital/analogue-converter using the counting method is a simple and cost-effective alternative to the binary weighted method
Generation of an analogue output signal from a PWM signal using a low-pass filter
Decoupling of input and output voltage using an operational amplifier circuit with the transfer function F (jω) =1
When the switch is closed, the capacitor C is charged via the resistor R
The time constant τ and therefore the charging curve depends on the dimensioning of the components
τ =R⋅C
Ideal result
Determined by the on-resistance of the switches and the accuracy of the negative feedback resistor
Non-linearity
Indicates by how much a step is greater or less than 1 LSB in the most unfavourable case
Non-linearity is usually below 0.5 LSB
Phone
Fax
E-Mail [email protected]
Internet www.ika.rwth-aachen.de