Biomedical Engineering Reference
In-Depth Information
EXAMPLE PROBLEM 4.15
Solve Eq. (4.75) and plot ventricular pressure
) for one human heartbeat. Use initial ventri-
cular volume of 150 ml and the parameter values in Table 4.3.
p
v
(
t
Solution
The following MATLAB m-file will perform the required computation and plot the results,
shown in Figure 4.42.
Isovolumic Ventricular Pressure
350
300
250
200
150
100
50
0
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Time [s]
FIGURE 4.42
Isovolumic ventricular pressure computed for a human heartbeat.
% ventricle.m
%
% MATLAB m-file to compute isovolumic pressure using ventricle model
%
% Initial conditions:
%
delt
¼
0.001; % The iteration time step delta t
a
¼
7e-4;
¼
b
20.;
¼
c
2.5;
¼
d
80.;
¼
tc
0.264;
¼
tp
0.371;
¼
tr
0.299;
¼
tb
0.258;
alpha
¼
2.88;
Continued
