BIOENG1320 Biological Signals and Systems MATLAB Project 2 -Solved

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/s2 ), 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.