Description
Biological Signals and Systems
Analyzing a biological system modeled by an LCCDE: As described in the lecture, the relationship
between cardiac output and aortic blood pressure can be represented with a second-order LCCDE
model of the arterial system. However, blood pressure is more easily measured in a peripheral
artery than the central aorta, and the peripheral blood pressure signal (ππ
(π‘)) differs
substantially from the central blood pressure signal (ππ
(π‘)). A third-order LCCDE model of the
arterial system may be more useful in that it can represent the relationship between the left
ventricular blood flow rate signal (πΜ
π
(π‘) whose average gives the cardiac output) and both ππ
(π‘)
and ππ
(π‘). Figure 1 shows the model in electrical analog form, where voltage corresponds to
blood pressure; charge, to blood volume; and current, to blood flow rate.
Figure 1: Third-order LCCDE model of arterial system.
The parameters of the model are πΆπ and πΆπ, which represent the central and peripheral
βcomplianceβ or blood volume storage ability; πΏπ
, which represents the central βinertanceβ and
indicates the pressure required to accelerate blood; and π
π, which represents the peripheral
resistance to blood flow. This circuit model is governed by the following two third-order LCCDEs.
πΏππ
ππΆππΆπ
π
3ππ
(π‘)
ππ‘
3 + πΏππΆπ
π
2ππ
(π‘)
ππ‘
2 + π
π(πΆπ + πΆπ)
πππ
(π‘)
ππ‘ + ππ
(π‘) = π
ππΜ
π
(π‘)
πΏππ
ππΆππΆπ
π
3ππ
(π‘)
ππ‘
3 + πΏππΆπ
π
2ππ
(π‘)
ππ‘
2 + π
π(πΆπ + πΆπ)
πππ
(π‘)
ππ‘ + ππ
(π‘)
= πΏππ
ππΆπ
π
2πΜ
π
(π‘)
ππ‘
2 + πΏπ
ππΜ
π
(π‘)
ππ‘ + π
ππΜ
π
(π‘)
2
1. Apply Laplace Transform techniques with βpaper and pencilβ to find the system functions,
π»π
(π ) = ππ
(π ) πΜ
π β (π ) and π»π
(π ) = ππ
(π ) πΜ
π β (π ).
2. Let πΆπ = 1.6 ml/mmHg, πΆπ = 0.2 ml/mmHg, πΏπ = 0.015 mmHg/(ml/s
2
), and π
π = 1.1
mmHg/(ml/s). Use the built-in functions, freqs, tf, bode, and impulse, to plot the frequency
response, Bode plot, and impulse response for π»π
(π ). Does this system show overdamped,
underdamped, or critically damped characteristics? What MATLAB function can you use to
verify your answer to this question?
3. Figure 2 shows a simple sine model of one beat of the input signal πΜ
π
(π‘). Create a vector to
define the sine beat signal at a sampling interval of 0.001 sec. Set ππ (stroke volume, which
is the amount of blood ejected by the left ventricle per beat) to 80 ml, and π (the beat length)
to 1 sec. Now create another vector to define a train of at least 20 unit-impulses spaced apart
by T with zeros inserted in between. Use the same sampling interval for this βimpulse trainβ.
Finally, form a periodic or βpulsatileβ input signal by convolving the sine beat signal with the
impulse train. Plot the input signal. Why does this method work?
Figure 2: Sine model of one beat of the left ventricular blood flow rate.
4. Use the built-in function lsim to determine the output ππ
(π‘) of the system with system
function π»π
(π ) in response to the periodic sine flow rate waveform. Plot the output signal.
How do you think this built-in function works? What is the cause of the transient part of the
output?
5. Determine the output ππ
(π‘) of the system with system function π»π
(π ) to the same input. Plot
the output signal. Compare the two steady-state outputs. Are the results what you expected?
6. Vary the model parameters (ππ, π, πΆπ
, and π
π) individually by Β±50% and compute the output
for each case. How do the cardiac output, mean blood pressure, and peripheral pulse pressure
(systolic minus diastolic blood pressure) change in response to each parametric variation?
Show exemplary plots to support your answers. What are biological models good for?
Deliverables: Submit a single pdf file containing answers to the questions, properly labeled
plots, and the source code.