DynamicSystems II SignalsSystems
DynamicSystems II SignalsSystems
Control 1
1
Signals and Systems
Outline
1. Introduction
2. Discrete systems, the z-transform
3. Continuous systems, the Laplace transform
4. Laplace and z-transform, in practice
5. Frequency domain
6. Continuous and discrete systems
● Examples
• Aircraft flying at 𝑈 = 10,000𝑓𝑡 : we will consider altitude small variations
𝑢 =+/−100𝑓𝑡
• Autocruise at 𝑈 = 100𝑘𝑚/ℎ : we will consider
𝑢 =+/−1𝑘𝑚/ℎ
𝑢 𝑦 0, ∀𝑛 ≠ 0
𝑢 𝑛 = 𝛿0 𝑛 = ቊ
ℎ 1, 𝑛=0 ℎ0 𝑛
⇒ 𝑦 𝑛 = ℎ0 (𝑛)
𝛿0 𝑛 𝛿1 𝑛 𝛿2 𝑛
1. « Time Invariant »
property:
ℎ0 𝑛 ℎ1 𝑛 ℎ2 𝑛
= + + + +⋯
𝑢 𝑛 = 𝑢0 𝛿0 𝑛 + 𝑢1 𝛿1 𝑛 + 𝑢2 𝛿2 𝑛 + 𝑢3 𝛿3 𝑛 +…
2. « Linear » property: = + + + +…
𝑦 𝑛 = 𝑢0 ℎ0 𝑛 + 𝑢1 ℎ1 𝑛 + 𝑢2 ℎ2 𝑛 + 𝑢3 ℎ3 𝑛 +⋯
𝑦 𝑛 = 𝑢 𝜏 ℎ 𝑛 −𝜏
𝜏=0
𝑧 ∈ ℂ, 𝑈 𝑧 = 𝑢 𝑛 𝑧 −𝑛
𝑛=−∞
𝑦 𝑛 = 𝑢 𝜏 ℎ 𝑛−𝜏
𝜏=0
(𝑑𝑒𝑚𝑜)
𝑌(𝑧) = 𝐻 𝑧 𝑈 𝑧
● Definition
𝑛=∞
𝑧 ∈ ℂ, 𝑈 𝑧 = 𝑢 𝑛 𝑧 −𝑛
𝑛=−∞
● Convergence
– We should always give the Region of Convergence (ROC): the set of points in the complex
plane where the z-transform summation converge…
𝑈 𝑧 = 𝑢 𝑛 𝑧 −𝑛 𝑈 𝑧 = 𝑢 𝑛 𝑧 −𝑛
𝑛=−∞ 𝑛=0
– When we use the z-transform for the impulse response of a causal system ℎ(𝑛) is always
equal to 0 for every negative value of 𝑛
● Inverse z-transform
1
𝑢 𝑛 = ර 𝑈 𝑧 𝑧 𝑛−1 𝑑𝑧
2𝜋𝑗 𝐶
– Where C is a counterclockwise path encircling the origin and entirely the ROC
● Linearity
𝛼𝑢 𝑛 + 𝛽𝑣 𝑛 ⇒ 𝛼𝑈 𝑧 + 𝛽𝑉 𝑧
● Delay
𝑢 𝑛 − 𝑑 ⇒ 𝑢 −𝑑 + 𝑢 −𝑑 + 1 𝑧 −1 + ⋯ 𝑢 −1 𝑧 −𝑑−1 + 𝑧 −𝑑 𝑈 𝑧
● Convolution
𝑢 𝑛 ∗𝑣 𝑛 ⇒𝑈 𝑧 𝑉 𝑧
● Difference
𝑢 𝑛 − 𝑢 𝑛 − 1 ⇒ 1 − 𝑧 −1 𝑈 𝑧
● Integration
𝑛 1
𝑢(𝑘) ⇒ 𝑈 𝑧
𝑘=0 1 − 𝑧 −1
𝛿(𝑛) 1
𝛿(𝑛 − 𝑑) 𝑧 −𝑑
1
𝑢(𝑛)
1 − 𝑧 −1
1
𝛼 𝑛 𝑢(𝑛)
1 − 𝛼𝑧 −1
𝑧 −1
𝑛𝑢(𝑛)
1 − 𝑧 −1 2
1 − cos 𝜔 𝑧 −1
cos 𝜔𝑛 𝑢(𝑛)
1 − 2 cos 𝜔 𝑧 −1 + 𝑧 −2
sin 𝜔 𝑧 −1
sin 𝜔𝑛 𝑢(𝑛)
1 − 2 cos 𝜔 𝑧 −1 + 𝑧 −2
1 − 𝑎 cos 𝜔 𝑧 −1
an cos 𝜔𝑛 𝑢(𝑛)
1 − 2 𝑎 cos 𝜔 𝑧 −1 + 𝑎2 𝑧 −2
𝑎 sin 𝜔 𝑧 −1
an sin 𝜔𝑛 𝑢(𝑛)
1 − 2 𝑎 cos 𝜔 𝑧 −1 + 𝑎2 𝑧 −2
𝐼𝑚𝑎𝑔
Stability region
𝑅𝑒𝑎𝑙
1
0, ∀𝑡 ≠ 0
𝑢 𝑡 =𝛿 𝑡 =ቊ 𝑤𝑖𝑡ℎ න 𝛿 𝜏 𝑑𝜏 = 1
∞, 𝑡 = 0
𝑦 𝑡 = ℎ(𝑡)
𝑦 𝑡 = න 𝑢 𝜏 ℎ 𝑡 − 𝜏 𝑑𝜏
𝜏=0
𝑠 ∈ ℂ, 𝑈 𝑠 = න 𝑢 𝑡 𝑒 −𝑠𝑡 𝑑𝑡
𝑡=0−
𝑦 𝑡 = න 𝑢 𝜏 ℎ 𝑡 − 𝜏 𝑑𝜏
𝜏=0
(𝑑𝑒𝑚𝑜)
𝑌(𝑠) = 𝐻 𝑠 𝑈 𝑠
● Definition
𝑡=∞
𝑠 ∈ ℂ, 𝑈 𝑠 = න 𝑢 𝑡 𝑒 −𝑠𝑡 𝑑𝑡
𝑡=0−
● Convergence
– The set of values for which U(𝑠) converges absolutely is for 𝑅𝑒 𝑠 > 𝑎 where 𝑎 is a
constant…
● Linearity
𝛼𝑢 𝑡 + 𝛽𝑣 𝑡 ⇒ 𝛼𝑈 𝑠 + 𝛽𝑉 𝑠
● Delay
𝑢 𝑡 − 𝑡0 ⇒ 𝑒 −𝑡0 𝑠 𝑈 𝑠
● Convolution
𝑢 𝑡 ∗𝑣 𝑡 ⇒𝑈 𝑠 𝑉 𝑠
Derivation
● Derivation
𝑑𝑢 𝑡
⇒ 𝑠𝑈 𝑠 − 𝑢(0− )
𝑑𝑡
● Integration Pre-initial value
𝑡
1
න 𝑢 𝜏 𝑑𝜏 ⇒ 𝐹 𝑠
0 𝑠
𝛿 𝑡 (unit impulse) 1
𝛿(𝑡 − 𝑡0 ) 𝑒 −𝑠𝑡0
1
𝑢(𝑛) (unit step)
𝑠
1
𝑡𝑢(𝑡)
𝑠2
1
𝑒 −𝑎𝑡 𝑢 𝑡
𝑠+𝑎
𝜔
sin 𝜔𝑡 𝑢 𝑡
𝑠2 + 𝜔2
𝑠
cos 𝜔𝑡 𝑢 𝑡
𝑠2 + 𝜔2
𝜔
𝑒 −𝛼𝑡 sin 𝜔𝑡
𝑠 + 𝛼 2 + 𝜔2
𝑠+𝑎
𝑒 −𝛼𝑡 cos 𝜔𝑡
𝑠 + 𝛼 2 + 𝜔2
This example catches most practical examples. To be more precise read the Partial Fraction
Decomposition page of Wikipedia
https://en.wikipedia.org/wiki/Partial_fraction_decomposition
Stability region
𝑅𝑒𝑎𝑙
● Stability
– Poles says if the system is stable or not (real part for 𝑠, radius for 𝑧)
𝑏1 𝑠 + 𝑏0 𝑠𝑦 0− + 𝑦ሶ 0− + 𝑎1 𝑦 0− − 𝑏0 𝑢(0− )
𝑌 𝑠 = 2 𝑈 𝑠 +
𝑠 + 𝑎1 𝑠 + 𝑎0 𝑠 2 + 𝑎1 𝑠 + 𝑎0
𝑢 𝑡 𝑦 0 , 𝑦ሶ 0 , 𝑢 0
𝑦 𝑡
𝑎1 , 𝑎0 , 𝑏1
Laplace Inverse
Transform Laplace Transform
𝑈 𝑠 𝑌 𝑠
𝑦 𝑡 = 𝐴 × 𝑀 𝜔 cos 𝜔𝑡 + 𝜙 𝜔
Output is also sinusoidal but:
• Amplitude amplified by magnitude of 𝐇 𝐢𝛚
• Phase shifted by phase of 𝑯 𝒊𝝎
• Input was sinusoidal forever: not transient response
● Bode diagram
– Vertical axis:
magnitude of 𝐻(𝑖𝜔) in dB
Phase of 𝐻(𝑖𝜔) in 𝑑𝑒𝑔
– Horizontal axis: frequency 𝑓 (𝐻𝑧) or pulsation 𝜔 (𝑟𝑎𝑑. 𝑠 −1 ) in log scale
20 log10 𝐻 𝑖𝜔
arg 𝐻 𝑖𝜔
– Combinations
𝑇
𝑦(𝑡) 𝑦𝑠 (𝑡)
× =
𝑘=+∞
𝑦𝑠 𝑡 = 𝑦 𝑡 × 𝛿 𝑡 − 𝑘𝑇
𝑘=−∞
𝑘=+∞
= 𝑦 𝑘𝑡 𝛿 𝑡 − 𝑘𝑇
𝑘=−∞
1/𝑇𝑠 ok
• Shannon theorem
• If the signal spectrum is limited in a band [0 𝑓𝑚𝑎𝑥 ] the sampling frequency must be chosen
such as 𝑓𝑠 > 2 × 𝑓𝑚𝑎𝑥
• If the sampling rate is large enough the continuous signal it can be fully
reconstructed by correctly filtering the sampled signal
• If the sampling rate is too small there is loss of information by aliasing
𝑦𝑠 (𝑡) 𝑦𝑠 (𝑡)
ZOH