0% found this document useful (0 votes)
23 views24 pages

SS03 - 2023S1 New

This document discusses representations of continuous-time linear time-invariant systems in the time domain. It covers modeling CT LTI systems using linear constant-coefficient differential equations and solving these equations using calculus approaches. It also discusses the zero-input response, zero-state response, and impulse response of CT LTI systems.

Uploaded by

hm.quy.9223
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
23 views24 pages

SS03 - 2023S1 New

This document discusses representations of continuous-time linear time-invariant systems in the time domain. It covers modeling CT LTI systems using linear constant-coefficient differential equations and solving these equations using calculus approaches. It also discusses the zero-input response, zero-state response, and impulse response of CT LTI systems.

Uploaded by

hm.quy.9223
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 24

ELT2035 Signals & Systems

Lesson 3: Representations of CT LTI


systems in time domain

Hoang Gia Hung


Faculty of Electronics and Telecommunications
University of Engineering and Technology, VNU Hanoi
Last lecture review
❑ Fundamental concepts of system
➢ System model: mathematical description of the rules of operation
➢ System analysis vs system synthesis (design)

❑ Classification of system
➢ By configurations (SISO, MIMO, …)
➢ By properties of the system (CT/DT, Instantaneous/dynamic,
causal/noncausal, time varying/time invariant, …)

❑ Properties of systems
➢ Stability: BIBO criterion
➢ Memory/Memoryless
➢ Causality: cannot “see” future inputs
➢ Invertibility: ability to restore the input from the output
➢ Time invariance: system behaviours do not change with time
➢ Linearity: superposition + homogeneity

❑ Today’s lecture: Modelling of CT LTI systems in time domain


Introduction

Input signal Output signal


(excitation) System (response)
x(t) H[.] y(t)

❑ Objectives: design/synthesis a system which will produce the


desired output for the given input. How?
➢ The first step is the construction of a system model (i.e. determine an
operator H[.] that satisfactorily approximates the behaviour of the system).
➢ Next: use the constructed system to analyze the system rules of operation

❑ System modelling techniques


➢ Internal description approach: H[.] is derived from a knowledge of the
internal structure of the system
➢ External description (“black box”) approach: H[.] is obtained from the
measurements of the system’s input and the corresponding output
A circuit modelling example
❑ Input: 𝑥(𝑡), output: 𝑦 𝑡 , what’s the system’s model?
➢ Model the resistor (Ohm’s law): 𝑣𝑅 𝑡 = 𝑅𝑦(𝑡)
𝑑𝑦(𝑡)
➢ Model the inductor (the 𝑖 − 𝑣 eq.): 𝑣𝐿 𝑡 = 𝐿
𝑑𝑡
➢ Model the circuit (Kirchoff’s valtage law):
𝑑𝑦(𝑡)
𝑅𝑦 𝑡 + 𝐿 = 𝑥(𝑡) (1)
𝑑𝑡
➢ Eq. (1): linear constant-coefficient ordinary differential equation (LCCODE)

❑ Generalization: a CT LTI system is described by the LCCODE


𝑛 𝑑 𝑖 𝑦(𝑡) 𝑚 𝑑 𝑗 𝑥(𝑡)
σ𝑖=0 𝑎𝑖 = σ𝑗=0 𝑏𝑗 , (2)
𝑑𝑡 𝑖 𝑑𝑡 𝑗
where 𝑎𝑖 and 𝑏𝑗 are constants.
𝑑
➢ Denote 𝐷 = , eq. (2) is rewritten as: 𝐷𝑛 + 𝑎𝑛−1 𝐷𝑛−1 + ⋯ + 𝑎0 𝑦 𝑡 =
𝑑𝑡
𝑏𝑚 𝐷 + 𝑏𝑚−1 𝐷𝑚−1 + ⋯ + 𝑏0 𝑥 𝑡 , or
𝑚

𝑄 𝐷 𝑦 𝑡 =𝑃 𝐷 𝑥 𝑡 . (3)
➢ Eq. (3) is obtained by applying natural laws (e.g. Ohm/Kirchoff laws) to the
internal structure of the system (e.g. components and topology in a circuit).
Solving for LCCODE – the calculus
approach
❑ System analysis: use the model to determine the response of
the system to any given input. Two methods:
➢ Solution to linear constant-coefficients differential equations (today)
➢ Laplace transform (later in the course)
𝑑𝑦(𝑡)
❑ The R-L circuit example: 𝑅𝑦 𝑡 + 𝐿 = 𝑥(𝑡). Given: 𝐿 =
𝑑𝑡
1 𝐻 , 𝑅 = 1 [Ω], 𝑥 𝑡 = cos(𝑡), 𝑦 0 = 2 [𝐴], find 𝑦 𝑡 .
𝑑𝑦𝑛 (𝑡)
➢ Homogeneous solution (natural response): 𝑦𝑛 𝑡 + = 0 → 𝑦𝑛 𝑡 =
𝑑𝑡
𝐴𝑛 𝑒 𝜆𝑡 →𝐴𝑛 𝑒 𝜆𝑡 1 + 𝜆 = 0 → 𝜆 = −1 → 𝑦𝑛 (𝑡) = 𝐴𝑛 𝑒 . −𝑡

𝑑𝑦𝑓 (𝑡)
➢ Particular solution (forced response): 𝑦𝑓 𝑡 + = cos(𝑡) → 𝑦𝑓 𝑡 =
𝑑𝑡
1
𝐵 cos(𝑡 + 𝐶) → 𝐵 cos 𝑡 + 𝐶 − 𝐵 sin 𝑡 + 𝐶 = cos(𝑡) → 𝑦𝑓 𝑡 = cos(𝑡) +
2
1
sin(𝑡).
2
1
➢ General solution (total response): 𝑦 𝑡 = 𝑦𝑛 𝑡 + 𝑦𝑓 𝑡 = 𝐴𝑛 𝑒 −𝑡 + cos(𝑡) +
2
1 1 3 3 1 1
sin(𝑡) → 𝑦 0 = 𝐴𝑛 + → 𝐴𝑛 = , 𝑦 𝑡 = 𝑒 −𝑡 + cos(𝑡) + sin(𝑡).
2 2 2 2 2 2
❖ Note: the initial condition is imposed to the total response, which begins right
after the input is applied at 𝑡 = 0 → the i.c. is required at 𝑡 = 0+ .
Solving for LCCODE – the engineering
approach
❑ In practice, it is difficult to measure initial conditions at 𝑡 = 0+ .
➢ It’s much easier to determine 𝑦 0− .
➢ 𝑦 0− = internal condition (e.g. 𝑡=0
𝑦 0−
energy storage) of the system prior
to the application of the input.

❑ With respect to initial conditions at 𝑡 = 0− , the total response is


the sum of two independent components:
Total response = zero-input response + zero-state response

➢ Z-I response is due to 𝑦 0− alone & independent of external input: 𝑦𝑍𝐼 𝑡 +


𝑑𝑦𝑍𝐼 (𝑡)
= 0 → 𝑦𝑍𝐼 (𝑡) = 𝐴𝑍𝐼 𝑒 −𝑡 → 𝑦𝑍𝐼 0 = 2 → 𝐴𝑍𝐼 = 2, 𝑦𝑍𝐼 (𝑡) = 2𝑒 −𝑡 .
𝑑𝑡
➢ Z-S response is caused by external input only, i.e. the internal condition is 0:
𝑑𝑦 (𝑡) 1 1
𝑦𝑍𝑆 𝑡 + 𝑍𝑆 = cos(𝑡) → 𝑦𝑍𝑆 𝑡 = 𝐴𝑍𝑆 𝑒 −𝑡 + cos(𝑡) + sin(𝑡) → 𝑦 0 =
𝑑𝑡 2 2
1 1 1 1
𝐴𝑍𝑆 + = 0 → 𝑦𝑍𝑆 𝑡 = − 𝑒 −𝑡 + cos(𝑡)
+ sin(𝑡).
2 2 2 2
3 1 1
➢ Total response: 𝑦 𝑡 = 𝑦𝑍𝐼 𝑡 + 𝑦𝑍𝑆 𝑡 = 𝑒 −𝑡 + cos(𝑡) + sin(𝑡).
2 2 2
Relations between the two approaches
❑ Natural response vs zero-input response:
➢ Both are solutions to 𝑄 𝐷 𝑦 𝑡 = 0 (with different initial conditions) → same
form, differs in constant coefficient of 𝑒 −𝑡 .
➢ natural response = zero-input response + all characteristic modes (i.e. 𝑒 −𝑡
term in this case) from zero-state response.

❑ Forced response vs. zero-state response:


➢ Both are solutions to 𝑄 𝐷 𝑦 𝑡 = 𝑃 𝐷 𝑥 𝑡 (different initial conditions) →
same form, differs in constant coefficient of 𝑒 −𝑡 .
➢ forced response = zero-state response - characteristic mode terms

𝟑 −𝑡
𝑦𝑍𝐼 𝑡 = 𝟐𝑒 −𝑡 𝑦𝑛 𝑡 = 𝑒
𝟐
𝟏 1 1 1 1
𝑦𝑍𝑆 𝑡 = − 𝑒 −𝑡 + cos(𝑡) + sin(𝑡) 𝑦𝑓 𝑡 = 𝟎. 𝑒 −𝑡 + cos(𝑡) + sin(𝑡)
𝟐 2 2 2 2
Zero-input response of CT-LTI systems
❑ Solution to 𝐷 𝑛 + 𝑎𝑛−1 𝐷 𝑛−1 + ⋯ + 𝑎0 𝑦0 𝑡 = 0 (4)
➢ Linear combination of 𝑦0 𝑡 and its 𝑛 successive derivatives are zero for all
t → 𝑦0 𝑡 and its derivatives must have the same form → 𝑦0 𝑡 = 𝑐𝑒 𝜆𝑡 .
➢ Substitute 𝑐𝑒 𝜆𝑡 for 𝑦0 𝑡 yields: 𝑐 𝜆𝑛 + 𝑎𝑛−1 𝜆𝑛−1 + ⋯ + 𝑎0 𝑒 𝜆𝑡 = 0 ∀𝑡 →
nontrivial solution satisfies 𝜆𝑛 + 𝑎𝑛−1 𝜆𝑛−1 + ⋯ + 𝑎0 = 0 or 𝑄 𝜆 = 0.
➢ Eq. (4) has 𝑛 possible solutions 𝑐1𝑒 𝜆1𝑡 , … , 𝑐𝑛𝑒 𝜆𝑛𝑡 → the general solution is
𝑦0 𝑡 = 𝑐1 𝑒 𝜆1 𝑡 + ⋯ + 𝑐𝑛 𝑒 𝜆𝑛 𝑡 , where 𝑐1:𝑛 are defined by 𝑛 initial conditions.

❑ 𝑄 𝜆 is referred to as the characteristic polynomial of the


system, while 𝑄 𝜆 = 0 is called the characteristic equation.
➢ 𝜆1:𝑛 are called characteristic roots, or eigenvalues, or natural frequencies.
➢ 𝑒 𝜆𝑖𝑡 are called characteristic/natural modes of the system.

❑ The zero-input response is a linear combination of the


system’s characteristic modes.
➢ Non-characteristic modes are lumped together in the forced response.
Zero-state response of CT LTI systems
to the unit impulse input
h(t) = H[δ(t)]

System
H[.]

❑ When 𝛿(𝑡) enters the input of the system:


➢ All the initial conditions are 0 at 𝑡 = 0− , then at 𝑡 = 0 the impulse creates
nonzero initial conditions at 𝑡 = 0+ → ℎ 𝑡 must consists of the system’s
characteristic modes for 𝑡 ≥ 0+ .
➢ 𝛿(𝑡) vanishes for 𝑡 > 0, but what happens on the system’s output at the
single moment 𝑡 = 0? There can at most be an impulse.

❑ For LTI systems: ℎ 𝑡 = 𝑏𝑛 𝛿 𝑡 + 𝑃 𝐷 𝑦𝑛 𝑡 𝑢(𝑡).


➢ 𝑦𝑛 𝑡 is a linear combination of the characteristic modes that satisfies:
𝑛 = 1: 𝑦𝑛 0 = 1
𝑛 = 2: 𝑦𝑛 0 = 0 and 𝑦ሶ𝑛 0 = 1
𝑛 = 3: 𝑦𝑛 0 = 𝑦ሶ𝑛 0 = 0 and 𝑦ሷ𝑛 0 = 1
𝑛 = 4: 𝑦𝑛 0 = 𝑦ሶ𝑛 0 = 𝑦ሷ𝑛 0 = 0 and 𝑦ഺ𝑛 (0) = 1
➢ ℎ 𝑡 is a descriptor of the system and is called the impulse response.
Impulse response example
Determine the unit impulse response ℎ(𝑡) for a system specified
by the equation 𝐷 2 + 3𝐷 + 2 𝑦 𝑡 = 𝐷𝑓(𝑡).
=?
f(t) = δ(t)

𝐷2 + 3𝐷 + 2 𝑦 𝑡 = 𝐷𝑓(𝑡)

Solution
➢ 𝑄 𝐷 = 𝐷2 + 3𝐷 + 2, 𝑃 𝐷 = 𝐷 ⟹ 𝑏2 = 0 ⟹ ℎ 𝑡 = 𝑏2 𝛿 𝑡 +
𝑃 𝐷 𝑦𝑛 𝑡 𝑢 𝑡 = 𝑦ሶ𝑛 𝑡 𝑢(𝑡).
➢ Characteristic polynomial: 𝜆2 + 3𝜆 + 2 = 0 ⟹ eigenvalues 𝜆1 =
− 1, 𝜆2 = −2 ⟹ 𝑦𝑛 𝑡 = 𝑐1 𝑒 −𝑡 + 𝑐2 𝑒 −2𝑡 .
➢ Since 𝑛 = 2, the initial conditions are: 𝑦ሶ𝑛 0 = 1, 𝑦𝑛 0 = 0.
➢ Differentiate 𝑦𝑛 𝑡 and substitute to the initial conditions yields
0 = 𝑐1 + 𝑐2
ቊ ⟹ 𝑐1 = 1, 𝑐2 = −1 ⟹ 𝑦𝑛 𝑡 = 𝑒 −𝑡 − 𝑒 −2𝑡 .
1 = −𝑐1 − 2𝑐2
➢ Hence, ℎ 𝑡 = 𝐷𝑦𝑛 𝑡 𝑢 𝑡 = −𝑒 −𝑡 + 2𝑒 −2𝑡 𝑢(𝑡).
Zero-state response of CT LTI systems
to an arbitrary input
❑ Given ℎ(𝑡), what’s the Z-S response to an arbitrary input 𝑓(𝑡)?
System
𝛿(𝑡) → ℎ(𝑡)

❑ Approximate 𝑓(𝑡) as a sum of narrow rectangular pulses

1
𝛿መ 𝑡 = ቐ∆𝜏 0 ≤ 𝑡 < ∆𝜏
0 otherwise


ℎ(𝑡)
𝛿መ 𝑡

Impulse approximation

෠ − 𝑛∆𝜏)
𝑓(𝑛∆𝜏)∆𝜏 ℎ(𝑡
መ − 𝑛∆𝜏)
𝑓(𝑛∆𝜏)∆𝜏 𝛿(𝑡

LTI properties
Zero-state response of CT LTI systems

𝑓 0 ∆𝜏 𝛿(𝑡) ෠
𝑓 0 ∆𝜏 ℎ(𝑡)
System
H[.]

መ − ∆𝜏) ෠ − ∆𝜏)
𝑓 ∆𝜏 ∆𝜏 ℎ(𝑡
𝑓 ∆𝜏 ∆𝜏 𝛿(𝑡

መ − 𝑛∆𝜏)
𝑓(𝑛∆𝜏)∆𝜏 𝛿(𝑡 ෠ − 𝑛∆𝜏)
𝑓(𝑛∆𝜏)∆𝜏 ℎ(𝑡

Superposition


𝑓(𝑡)
≈ 𝑓(𝑡)
𝑦(𝑡)

≈ 𝑦(𝑡)
Convolution

❑ In the limit, approximations become exact representations


➢ The operator that performs superposition of the scaled, shifted impulse
responses to produce the system’s output is called convolution.
Derivation of the convolution integral

System
H[.]

❑ Applying linearity and time-invariant properties:



𝛿(𝑡) → ෠
ℎ(𝑡)
መ − 𝑛∆𝜏)
𝑓(𝑛∆𝜏)∆𝜏 𝛿(𝑡 → ෠ − 𝑛∆𝜏)
𝑓(𝑛∆𝜏)∆𝜏 ℎ(𝑡
∞ ∞
መ − 𝑛∆𝜏) →
෍ 𝑓(𝑛∆𝜏)∆𝜏 𝛿(𝑡 ෠ − 𝑛∆𝜏)
෍ 𝑓(𝑛∆𝜏)∆𝜏 ℎ(𝑡
𝑛=−∞ 𝑛=−∞

𝑓(𝑡)≈ 𝑓(𝑡) ො
𝑦(𝑡)≈ 𝑦(𝑡)

❑ Taking the limit:


∞ ∞
መ − 𝑛∆𝜏)∆𝜏
lim ෍ 𝑓(𝑛∆𝜏)𝛿(𝑡 → ෠ − 𝑛∆𝜏)∆𝜏
lim ෍ 𝑓 𝑛∆𝜏 ℎ(𝑡
∆𝜏→0 ∆𝜏→0
𝑛=−∞ 𝑛=−∞
∞ ∞
න 𝑓 𝜏 𝛿 𝑡 − 𝜏 𝑑𝜏 → න 𝑓 𝜏 ℎ 𝑡 − 𝜏 𝑑𝜏
−∞ −∞
= 𝑓(𝑡) = 𝑦(𝑡)

❑ Convolution: 𝑦 𝑡 = ‫׬‬−∞ 𝑓 𝜏 ℎ 𝑡 − 𝜏 𝑑𝜏 .
Convolution integral evaluation
❑ Find 𝑦 𝑡 = 𝑥1 𝑡 ∗ 𝑥2 𝑡 ➢ Step 1

= ‫׬‬−∞ 𝑥1 𝜆 𝑥2 𝑡 − 𝜆 𝑑𝜆.

➢ Steps 2 → 4 for 2 ≤ 𝑡 < 3

❑ Steps in convolution operation


1. Fold/Reflect 𝑥2 𝑡 → 𝑥2 (−𝜆).
2. Shift/Slide 𝑥2 (−𝜆) by 𝒕 (seen as a
➢ Steps 2 → 4 for 3 ≤ 𝑡 < 4
constant inside the integration).
3. Find 𝑔𝑡 𝜆 = 𝑥1 𝜆 𝑥2 𝑡 − 𝜆 .
4. Integrate 𝑔𝑡 𝜆 over 𝝀.
5. Repeat 2 → 4 for all values of 𝑡.
❑ Result:
0, 𝑡<2 ➢ Steps 2 → 4 for 4 ≤ 𝑡 ≤ 5
𝑡−2 2≤𝑡 <3
𝑦 𝑡 = 1 3≤𝑡<4
5−𝑡 4≤𝑡 <5
0 𝑡≥5
Impulse response - a descriptor of LTI
systems
❑ Given ℎ(𝑡), we can determine the system’s output 𝑦(𝑡) for any
input 𝑥(𝑡) by performing the convolution 𝑦 𝑡 = ℎ(𝑡) ∗ 𝑥(𝑡).
❑ Both impulse response and differential equation are used to
describe the input-output relationship of CT LTI systems, but:
➢ The impulse response description applies only to systems that are initially
at rest → can determine the system’s zero−state response only.
➢ Differential equation representation is more flexible: it applies to systems
which are either at rest or with non-zero initial conditions.

❑ Relationship between the responses from each description:


Total response = σ𝑛𝑗=1 𝑐𝑗 𝑒 𝜆𝑗𝑡 + ℎ 𝑡 ∗ 𝑥(𝑡)
zero−input response zero−state response

❑ Both descriptions do NOT provide insight into the system’s


internal mechanism to produce the output from the input.
➢ They are declarative (i.e. tell rules that the system obeys), while the block
diagram is imperative (i.e. tells what the system does with signals).
Properties of the convolution integral
(proofs: homework)
❑ Commutative:
𝑓1 𝑡 ∗ 𝑓2 𝑡 = 𝑓2 𝑡 ∗ 𝑓1 𝑡

❑ Distributive:
𝑓1 𝑡 ∗ 𝑓2 𝑡 + 𝑓3 𝑡 = 𝑓1 𝑡 ∗ 𝑓2 𝑡 + 𝑓1 𝑡 ∗ 𝑓3 𝑡

❑ Associative:
𝑓1 𝑡 ∗ 𝑓2 𝑡 ∗ 𝑓3 𝑡 = 𝑓1 𝑡 ∗ 𝑓2 𝑡 ∗ 𝑓3 𝑡

❑ The shift property: if 𝑓1 𝑡 ∗ 𝑓2 𝑡 = 𝑐(𝑡) then


𝑓1 𝑡 ∗ 𝑓2 𝑡 − 𝑇 = 𝑐 𝑡 − 𝑇 , and
𝑓1 𝑡 − 𝑇 ∗ 𝑓2 𝑡 = 𝑐 𝑡 − 𝑇 , and
𝑓1 𝑡 − 𝑇1 ∗ 𝑓2 𝑡 − 𝑇2 = 𝑐 𝑡 − 𝑇1 − 𝑇2 .

❑ Convolution with an impulse


𝑓 𝑡 ∗ 𝛿 𝑡 = 𝑓(𝑡)
Impulse response as an indicator of
system’s properties

❑ The output of a CT LTI system: 𝑦 𝑡 = ‫׬‬−∞ ℎ 𝜏 𝑥 𝑡 − 𝜏 𝑑𝜏
➢ Memoryless:
ℎ 𝑡 = 𝑐𝛿 𝑡 , where c is an arbitrary constant
➢ Causal:
ℎ 𝑡 = 0 for all 𝑡 < 0.
➢ Stable: ℎ 𝑡 is absolutely integrable, i.e.

‫׬‬−∞ ℎ 𝜏 𝑑𝜏 < ∞.

❑ Invertibility and deconvolution


➢ A system is invertible if the input can be recovered from the output except
for a constant scale factor.
➢ An inverve system (of an invertible system) recovers 𝑥 𝑡 from ℎ(𝑡) ∗ 𝑥(𝑡)
→ an inverse system performs deconvolution: ℎ 𝑡 ∗ ℎinv 𝑡 = 𝛿(𝑡).
➢ Real-life application: equalizers are inverse systems that remove signal
distortion in communication modems or noise-cancelling headphones.
Recap: block diagram representation
❑ Use basic network elements to represent basic operations
performed along the signal flow paths.
𝑡
න 𝑦 𝑡 = න 𝑥 𝑡 𝑑𝑡
−∞

❑ Allows explicit representation of the system.


➢ Impulse responses or differential equations are mathematically compact
while block diagrams are more “natural” as they illustrate signal flow paths.

❑ Note: the block diagram description of a system is not unique.


❑ Example: find a block diagram representation for the system
described by the following differential equation
𝑑
𝑦 𝑡 + 10𝑦 𝑡 = 2𝑥 𝑡
𝑑𝑡
Example solution
𝑑 1 𝑑 1
❑ Form I: 𝑑𝑡 𝑦 𝑡 + 10𝑦 𝑡 = 2𝑥 𝑡 ⟺ 2 𝑑𝑡 𝑦 𝑡 = 𝑥 𝑡 − 10 𝑦 𝑡
2

𝑑
❑ Form II: 𝑑𝑡 𝑦 𝑡 + 10𝑦 𝑡 = 2𝑥 𝑡 ⟺ 𝑦 𝑡 = 2 ‫ 𝑡𝑑 𝑡 𝑥 ׬‬− 10 ‫𝑡𝑑 𝑡 𝑦 ׬‬
State variable description of LTI
systems

න න

❑ The state of a system is a minimal set of signals that


represents the system’s entire memory of the past.
➢ Each signal in the selected set is called a state variables & the selection of
the state of a system is not unique.
➢ A convenient way to select state variables are picking the outputs of the
integrators in the block diagram.

❑ An LTI system can be represented by a set of 1st-order


differential equations describing operations in its block diagram
➢ The state equations characterize the dynamics of the system’s state.
➢ The output equation relates the output to the system’s state and the input.
Example

න න

4. State variable description is often


1. Select state variables: 𝑞1 𝑡 , 𝑞2 𝑡 written in matrix form as follows
➢ The inputs to the integrators are the 𝒒ሶ 𝑡 = 𝑨𝒒 𝑡 + 𝒃𝑥 𝑡
derivatives of the state variables 𝑦 𝑡 = 𝒄𝒒 𝑡 + 𝐷𝑥(𝑡)
2. Write down the state equations 5. The state equation in matrix form
𝑑
𝑞 𝑡 = 2𝑞1 𝑡 − 𝑞2 𝑡 + 𝑥(𝑡) 𝑞ሶ 1 (𝑡) 2 −1 𝑞1 𝑡 1
𝑑𝑡 1 = + 𝑥(𝑡)
𝑑 𝑞ሶ 2 (𝑡) 1 0 𝑞2 𝑡 0
𝑞 𝑡 = 𝑞1 𝑡
𝑑𝑡 2
6. The output equation in matrix form
3. Write down the output equation
𝑞 𝑡
𝑦 𝑡 = 3𝑞1 𝑡 + 𝑞2 𝑡 𝑦 𝑡 = 3 1 1 + 0 𝑥(𝑡)
𝑞2 𝑡
Exploring concepts with Octave
❑ Simulating CT convolution
Ts=0.001;
x=[ones(1/Ts,1);zeros(3/Ts,1)];
y=[zeros(4/Ts,1);ones(2/Ts,1)];
z=conv(x,y)*Ts;


Exploring concepts with Octave
❑ Impulse response
➢ E.g. 𝐷2 + 3𝐷 + 2 𝑦 𝑡 = 𝐷𝑓(𝑡)
pkg load control;
sys=tf([1 0],[1 3 2]);
impulse(sys);
Summary
❑ Modelling of CT LTI systems
➢ Internal description approach: applying natural laws to the internal
structure of the system to obtain LCCODE.
➢ External description (“black box”) approach: measure the response of the
system to the unit impulse response.
➢ Other models: block diagram & state-variable representations.

❑ System analysis: given the system & input, what’s the output?
➢ Solving LCCODE: total response = natural response + forced response =
zero-input response + zero-state response.

➢ Convolution of the input and impulse response: 𝑦 𝑡 = ‫׬‬−∞ 𝑓 𝜏 ℎ 𝑡 − 𝜏 𝑑𝜏.

❑ System’s properties can be inferred from its impulse response


➢ Memoryless if ℎ(𝑡)=𝑐𝛿(𝑡), where 𝑐 is an arbitrary constant.
➢ Causal if ℎ 𝑡 = 0 for all 𝑡 < 0.

➢ Stable if ‫׬‬−∞ ℎ 𝜏 𝑑𝜏 < ∞.

❑ Next lecture: Modelling of DT LTI systems in time domain

You might also like