% The function read_param.m reads a file which contains the
% parameter values characterizing the human cardiovascular model
% and its execution.
%
% Note that if the parameters characterizing the model are modified,
% then this file must be accordingly modified.
%
% Function argument
% inputfileid - id of the file containing the parameter values
%
% Function output
% th - 113x1 MATLAB vector containing the read-in parameter values
% waveform - string of numerical values corresponding to the waveforms
% to be viewed on-line
%
function [th,waveform] = read_param(inputfileid)
% Pre-allocating all of the parameters.
Cls = [];
Cld = [];
Ca = [];
Cv = [];
Crs = [];
Crd = [];
Cpa = [];
Cpv = [];
Qlo = [];
Qao = [];
Qvo = [];
Qro = [];
Qpao = [];
Qpvo = [];
Rl = [];
Ra = [];
Rv = [];
Rr = [];
Rpa = [];
Rpv = [];
Qtot = [];
F = [];
Qlmax = [];
Qrmax = [];
Ppa = [];
Pa = [];
Pv = [];
Plmax = [];
Prmax = [];
Pms = [];
Qfrs = [];
againr = [];
bgain = [];
pgain = [];
Pasp = [];
Pratrsp = [];
Tr = [];
fcp = [];
Sa = [];
stdwr = [];
stdwf = [];
dgainp = [];
dgains = [];
againq = [];
cgainr = [];
cgainq = [];
cgainfp = [];
Sc = [];
Qt = [];
Pvs = [];
Crds = [];
Pas = [];
Pvc = [];
again = [];
againc = [];
cgainc = [];
againfs = [];
cgainfs = [];
againfp = [];
Qfr = [];
Qds = [];
Rair = [];
Clu = [];
Qluo = [];
window = [];
annotations = [];
numerics = [];
alpha = [];
step = [];
% Reading in each of the values of each of the above parameters
% from the input file according to the prescribed format of the
% file.
line = fgetl(inputfileid);
while (strcmp(line,'%%%END OF FILE%%%') == 0)
if (isempty(findstr(line,'%')) == 1)
if (isempty(line) == 1)
;
elseif (isempty(findstr(line,':')) == 1)
;
else
index = findstr(line,':');
tstr = char(line(1:index-1));
count = 1;
x = isspace(tstr(count));
mbintscalar(x)
while (x)
count=count+1;
x = isspace(tstr(count));
mbintscalar(x);
end;
start = count;
count = index-1;
x = isspace(tstr(count));
mbintscalar(x)
while (x)
count=count-1;
x = isspace(tstr(count));
mbintscalar(x);
end
finish = count;
tstr = char(tstr(start:finish));
if (strcmp(tstr,'Cls') == 1)
tmp = sscanf(line,'%s %f');
Cls = tmp(length(tmp));
elseif (strcmp(tstr,'Cld') == 1)
tmp = sscanf(line,'%s %f');
Cld = tmp(length(tmp));
elseif (strcmp(tstr,'Ca') == 1)
tmp = sscanf(line,'%s %f');
Ca = tmp(length(tmp));
elseif (strcmp(tstr,'Cv') == 1)
tmp = sscanf(line,'%s %f');
Cv = tmp(length(tmp));
elseif (strcmp(tstr,'Crs') == 1)
tmp = sscanf(line,'%s %f');
Crs = tmp(length(tmp));
elseif (strcmp(tstr,'Crd') == 1)
tmp = sscanf(line,'%s %f');
Crd = tmp(length(tmp));
elseif (strcmp(tstr,'Cpa') == 1)
tmp = sscanf(line,'%s %f');
Cpa = tmp(length(tmp));
elseif (strcmp(tstr,'Cpv') == 1)
tmp = sscanf(line,'%s %f');
Cpv = tmp(length(tmp));
elseif (strcmp(tstr,'Qlo') == 1)
tmp = sscanf(line,'%s %f');
Qlo = tmp(length(tmp));
elseif (strcmp(tstr,'Qao') == 1)
tmp = sscanf(line,'%s %f');
Qao = tmp(length(tmp));
elseif (strcmp(tstr,'Qvo') == 1)
tmp = sscanf(line,'%s %f');
Qvo = tmp(length(tmp));
elseif (strcmp(tstr,'Qro') == 1)
tmp = sscanf(line,'%s %f');
Qro = tmp(length(tmp));
elseif (strcmp(tstr,'Qpao') == 1)
tmp = sscanf(line,'%s %f');
Qpao = tmp(length(tmp));
elseif (strcmp(tstr,'Qpvo') == 1)
tmp = sscanf(line,'%s %f');
Qpvo = tmp(length(tmp));
elseif (strcmp(tstr,'Rl') == 1)
tmp = sscanf(line,'%s %f');
Rl = tmp(length(tmp));
elseif (strcmp(tstr,'Ra') == 1)
tmp = sscanf(line,'%s %f');
Ra = tmp(length(tmp));
elseif (strcmp(tstr,'Rv') == 1)
tmp = sscanf(line,'%s %f');
Rv = tmp(length(tmp));
elseif (strcmp(tstr,'Rr') == 1)
tmp = sscanf(line,'%s %f');
Rr = tmp(length(tmp));
elseif (strcmp(tstr,'Rpa') == 1)
tmp = sscanf(line,'%s %f');
Rpa = tmp(length(tmp));
elseif (strcmp(tstr,'Rpv') == 1)
tmp = sscanf(line,'%s %f');
Rpv = tmp(length(tmp));
elseif (strcmp(tstr,'Qtot') == 1)
tmp = sscanf(line,'%s %f');
Qtot = tmp(length(tmp));
elseif (strcmp(tstr,'F') == 1)
tmp = sscanf(line,'%s %f');
F = tmp(length(tmp));
% elseif (strcmp(tstr,'Pth') == 1)
% tmp = sscanf(line,'%s %f');
% Pth = tmp(length(tmp));
% elseif (strcmp(tstr,'Ppac') == 1)
% tmp = sscanf(line,'%s %f');
% Ppac = tmp(length(tmp));
elseif (strcmp(tstr,'Qlmax') == 1)
tmp = sscanf(line,'%s %f');
Qlmax = tmp(length(tmp));
elseif (strcmp(tstr,'Qrmax') == 1)
tmp = sscanf(line,'%s %f');
Qrmax = tmp(length(tmp));
elseif (strcmp(tstr,'Ppa') == 1)
tmp = sscanf(line,'%s %f');
Ppa = tmp(length(tmp));
elseif (strcmp(tstr,'Pa') == 1)
tmp = sscanf(line,'%s %f');
Pa = tmp(length(tmp));
elseif (strcmp(tstr,'Pv') == 1)
tmp = sscanf(line,'%s %f');
Pv = tmp(length(tmp));
elseif (strcmp(tstr,'Plmax') == 1)
tmp = sscanf(line,'%s %f');
Plmax = tmp(length(tmp));
elseif (strcmp(tstr,'Prmax') == 1)
tmp = sscanf(line,'%s %f');
Prmax = tmp(length(tmp));
elseif (strcmp(tstr,'Pms') == 1)
tmp = sscanf(line,'%s %f');
Pms = tmp(length(tmp));
elseif (strcmp(tstr,'Qfrs') == 1)
tmp = sscanf(line,'%s %f');
Qfrs = tmp(length(tmp));
elseif (strcmp(tstr,'againr') == 1)
tmp = sscanf(line,'%s %f');
againr = tmp(length(tmp));
elseif (strcmp(tstr,'bgain') == 1)
tmp = sscanf(line,'%s %f');
bgain = tmp(length(tmp));
elseif (strcmp(tstr,'pgain') == 1)
tmp = sscanf(line,'%s %f');
pgain = tmp(length(tmp));
elseif (strcmp(tstr,'Pasp') == 1)
tmp = sscanf(line,'%s %f');
Pasp = tmp(length(tmp));
elseif (strcmp(tstr,'Pratrsp') == 1)
tmp = sscanf(line,'%s %f');
Pratrsp = tmp(length(tmp));
elseif (strcmp(tstr,'Tr') == 1)
tmp = sscanf(line,'%s %f');
Tr = tmp(length(tmp));
elseif (strcmp(tstr,'fco') == 1)
tmp = sscanf(line,'%s %f');
fco = tmp(length(tmp));
elseif (strcmp(tstr,'Sa') == 1)
tmp = sscanf(line,'%s %f');
Sa = tmp(length(tmp));
elseif (strcmp(tstr,'stdwr') == 1)
tmp = sscanf(line,'%s %f');
stdwr = tmp(length(tmp));
elseif (strcmp(tstr,'stdwf') == 1)
tmp = sscanf(line,'%s %f');
stdwf = tmp(length(tmp));
elseif (strcmp(tstr,'dgainp') == 1)
tmp = sscanf(line,'%s %f');
dgainp = tmp(length(tmp));
elseif (strcmp(tstr,'dgains') == 1)
tmp = sscanf(line,'%s %f');
dgains = tmp(length(tmp));
elseif (strcmp(tstr,'againq') == 1)
tmp = sscanf(line,'%s %f');
againq = tmp(length(tmp));
elseif (strcmp(tstr,'cgainr') == 1)
tmp = sscanf(line,'%s %f');
cgainr = tmp(length(tmp));
elseif (strcmp(tstr,'cgainq') == 1)
tmp = sscanf(line,'%s %f');
cgainq = tmp(length(tmp));
elseif (strcmp(tstr,'cgainfp') == 1)
tmp = sscanf(line,'%s %f');
cgainfp = tmp(length(tmp));
elseif (strcmp(tstr,'Sc') == 1)
tmp = sscanf(line,'%s %f');
Sc = tmp(length(tmp));
elseif (strcmp(tstr,'Qt') == 1)
tmp = sscanf(line,'%s %f');
Qt = tmp(length(tmp));
elseif (strcmp(tstr,'Pvs') == 1)
tmp = sscanf(line,'%s %f');
Pvs = tmp(length(tmp));
elseif (strcmp(tstr,'Crds') == 1)
tmp = sscanf(line,'%s %f');
Crds = tmp(length(tmp));
elseif (strcmp(tstr,'Pas') == 1)
tmp = sscanf(line,'%s %f');
Pas = tmp(length(tmp));
elseif (strcmp(tstr,'Pvc') == 1)
tmp = sscanf(line,'%s %f');
Pvc = tmp(length(tmp));
elseif (strcmp(tstr,'again') == 1)
tmp = sscanf(line,'%s %f');
again = tmp(length(tmp));
elseif (strcmp(tstr,'againc') == 1)
tmp = sscanf(line,'%s %f');
againc = tmp(length(tmp));
elseif (strcmp(tstr,'cgainc') == 1)
tmp = sscanf(line,'%s %f');
cgainc = tmp(length(tmp));
elseif (strcmp(tstr,'againfs') == 1)
tmp = sscanf(line,'%s %f');
againfs = tmp(length(tmp));
elseif (strcmp(tstr,'cgainfs') == 1)
tmp = sscanf(line,'%s %f');
cgainfs = tmp(length(tmp));
elseif (strcmp(tstr,'againfp') == 1)
tmp = sscanf(line,'%s %f');
againfp = tmp(length(tmp));
elseif (strcmp(tstr,'time') == 1)
tmp = sscanf(line,'%s %f');
time = tmp(length(tmp));
elseif (strcmp(tstr,'preparation') == 1)
tmp = sscanf(line,'%s %f');
preparation = tmp(length(tmp));
elseif (strcmp(tstr,'breathing') == 1)
tmp = sscanf(line,'%s %f');
breathing = tmp(length(tmp));
elseif (strcmp(tstr,'dncm') == 1)
tmp = sscanf(line,'%s %f');
dncm = tmp(length(tmp));
elseif (strcmp(tstr,'baro') == 1)
tmp = sscanf(line,'%s %f');
baro = tmp(length(tmp));
elseif (strcmp(tstr,'dra') == 1)
tmp = sscanf(line,'%s %f');
dra = tmp(length(tmp));
elseif (strcmp(tstr,'df') == 1)
tmp = sscanf(line,'%s %f');
df = tmp(length(tmp));
elseif (strcmp(tstr,'Fs') == 1)
tmp = sscanf(line,'%s %f');
Fs = tmp(length(tmp));
elseif (strcmp(tstr,'Qfr') == 1)
tmp = sscanf(line,'%s %f');
Qfr = tmp(length(tmp));
elseif (strcmp(tstr,'Qds') == 1)
tmp = sscanf(line,'%s %f');
Qds = tmp(length(tmp));
elseif (strcmp(tstr,'Rair') == 1)
tmp = sscanf(line,'%s %f');
Rair = tmp(length(tmp));
elseif (strcmp(tstr,'Clu') == 1)
tmp = sscanf(line,'%s %f');
Clu = tmp(length(tmp));
elseif (strcmp(tstr,'Qluo') == 1)
tmp = sscanf(line,'%s %f');
Qluo = tmp(length(tmp));
% elseif (strcmp(tstr,'timeql') == 1)
% tmp = sscanf(line,'%s %f');
% timeql = tmp(length(tmp));
elseif (strcmp(tstr,'window') == 1)
tmp = sscanf(line,'%s %f');
window = tmp(length(tmp));
elseif (strcmp(tstr,'annotations') == 1)
tmp = sscanf(line,'%s %f');
annotations = tmp(length(tmp));
elseif (strcmp(tstr,'numerics') == 1)
tmp = sscanf(line,'%s %f');
numerics = tmp(length(tmp));
elseif (strcmp(tstr,'waveform') == 1)
lnum = findstr(line,':');
waveform = line(lnum+2:length(line));
elseif (strcmp(tstr,'alpha') == 1)
tmp = sscanf(line,'%s %f');
alpha = tmp(length(tmp));
elseif (strcmp(tstr,'step') == 1)
tmp = sscanf(line,'%s %f');
step = tmp(length(tmp));
end
end
end
tstr = [];
line = [];
line = fgetl(inputfileid);
end
% Assigning these values to the th vector.
theta = [Cls*((Plmax)/(Qlmax-Qlo)) Cld*((Plmax)/(Qlmax-Qlo)) Ca Cv Crs*((Prmax)/(Qrmax-Qro)) Crd*((Prmax)/(Qrmax-Qro)) Cpa Cpv Qlo Qao Qvo Qro Qpao Qpvo Rl Ra Rv Rr Rpa Rpv Qtot F Pvc-((Qfr-Qluo)/Clu) 1/F Pvc Qlmax Qrmax Ppa Pa Pv Plmax Prmax Pms 3 5 Qfrs againr bgain pgain Pasp Pratrsp Tr 50 fco 2 Sa stdwr stdwf 2 2 2 dgainp dgains againq cgainr cgainq 0 cgainfp alpha Sc 0.05 0.05 6 0.0575 0.1 0.00016 6.75 numerics 0.025 -0.01 865.4 1.45 0.15 647.9687 67.0312 784.2687 Qt step 0 1.25 2.5 5 0.003 2 4 5 0.003 Pvs]';
theta = [theta; Crds; Pas; Pvc; again; againc; cgainc; againfs; cgainfs; againfp; Qfr; Qds; Rair; Clu; Qluo; 300; window; annotations];
th = [theta; time*Fs; preparation; breathing; dncm; baro; dra; df; Fs];
% If 113 parameters have not been assigned numerical values,
% then return a scalar -1 which will cause rcvsim to abort.
if (length(th) < 113)
th = [];
th = -1;
end