Download Full Issue in PDF - Academy Publisher
Download Full Issue in PDF - Academy Publisher
Download Full Issue in PDF - Academy Publisher
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
JOURNAL OF COMPUTERS, VOL. 8, NO. 6, JUNE 2013 1473<br />
B. Cardiovascular Function Parameters K Comput<strong>in</strong>g<br />
Based on Pulse Wave<br />
Needless to say, the characteristic <strong>in</strong>formation of pulse<br />
wave is closely related with the physiological factors. To<br />
study the relationship each other, many researchers get<br />
<strong>in</strong>formation from the time doma<strong>in</strong> or frequency doma<strong>in</strong><br />
characteristics based on pulse wave <strong>in</strong> cl<strong>in</strong>ical trials or<br />
model. In the time doma<strong>in</strong>, usually the pulse extracts<br />
some po<strong>in</strong>t with a clear physical mean<strong>in</strong>g (such as the<br />
ma<strong>in</strong> wave peak, heavy pump wave height, etc.). The<br />
comb<strong>in</strong>ation with the characteristic po<strong>in</strong>ts and the<br />
correspond<strong>in</strong>g physiological factors may get much<br />
cl<strong>in</strong>ical value. Some researchers have used simulation<br />
models to measure pulse wave different model<br />
parameters, to determ<strong>in</strong>e the person's physical condition<br />
accord<strong>in</strong>g to different parameters. Facts have proved that<br />
this method is more effective, but the simulation models<br />
and the extracted characteristic parameters must be<br />
proper, can effectively dist<strong>in</strong>guish the pathological state.<br />
In many studies, because the extracted parameters are<br />
too complicated to make the dist<strong>in</strong>ction between pulse<br />
waves, it often occurs the misjudg<strong>in</strong>g phenomenon.<br />
Therefore, the extraction of the pulse wave parameter is<br />
the critical research. Professor Luo Zhichang used the<br />
exist<strong>in</strong>g two-chamber model of elastic wave pulse to<br />
extract the characteristics of K (called form factor) which<br />
represents changes of the pulse wave’s area [11].<br />
Through the model theoretical analysis, thousands of<br />
animal experiments and cl<strong>in</strong>ical test<strong>in</strong>g with different<br />
age’s healthy people and patients with cardiovascular<br />
disease, confirmed that caused the pulse wave map<br />
features and the correspond<strong>in</strong>g changes <strong>in</strong> the area by<br />
physiological and pathological cardiovascular changes,<br />
and then reflect on the changes <strong>in</strong> K value. Determ<strong>in</strong>e the<br />
body's physical condition with the K value, although it<br />
can not achieve accurate quantitative analysis, but a<br />
simple calculation, differentiation, and the advantages of<br />
high sensitivity, which is important <strong>in</strong> the cl<strong>in</strong>ical<br />
reference value, is an important physiological <strong>in</strong>dicators<br />
of the cardiovascular cl<strong>in</strong>ical exam<strong>in</strong>ation.<br />
The K value reflects the characteristic quantities which<br />
changes <strong>in</strong> the amount of area of the pulse wave [11],<br />
which is def<strong>in</strong>ed as the average of the relative position of<br />
the pulse wave, which is def<strong>in</strong>ed by type (1) and figure 1.<br />
Pm<br />
− Pd<br />
K = . (1)<br />
P − P<br />
where K is the form factor; P<br />
m<br />
is mean arterial pressure,<br />
P<br />
d<br />
is the diastolic blood pressure; P s<br />
is the systolic<br />
blood pressure.<br />
Thus, the form factor K value have noth<strong>in</strong>g to do with<br />
the absolute value of systolic and diastolic blood pressure,<br />
it only depends on wave map area of the pulse wave, is a<br />
dimensionless parameter. Pulse waveform and area will<br />
have a great change <strong>in</strong> different physiological and<br />
pathological conditions, these changes can be expressed<br />
as K value.<br />
Because the pulse wave is difficult to accurately<br />
measure and solve, this paper propose a solution based on<br />
s<br />
d<br />
the averag<strong>in</strong>g method that K value of the pulse wave can<br />
be computed, and then through the specific network<br />
simulation of the cerebral circulation, the results is<br />
co<strong>in</strong>cided with cl<strong>in</strong>ical measurement.<br />
Ⅲ. CARDIOVASCULAR NETWORK HEMODYNAMIC<br />
ANALYSIS AND AVERAGING COMPUTATION<br />
Accord<strong>in</strong>g to the aforementioned study, <strong>in</strong> order to<br />
solve the pulse wave of blood circulation network<br />
diagram, first, analyze its network model [12-18]. In<br />
order to build blood circulation network's model, at first<br />
establishes the dynamic equation of one blood vessel<br />
branch. For simplicity, we make the follow<strong>in</strong>g<br />
assumptions: Al. the blood is <strong>in</strong>compressible; A2.the<br />
temperatures <strong>in</strong> all branches are identical. Under<br />
assumptions Al and A2, one branch of the blood network<br />
is described with the follow<strong>in</strong>g equations [12-18]:<br />
dQ<br />
j<br />
T<br />
j<br />
= −R<br />
j<br />
Q<br />
j<br />
Q<br />
dt<br />
TQ<br />
2<br />
= −Q<br />
R + H<br />
D<br />
j<br />
+ H<br />
j<br />
. (2)<br />
where Q<br />
j<br />
is flow through a branch j , R<br />
j<br />
are<br />
hemodynamic resistances, H<br />
j<br />
are pressure drops of the<br />
branches, Tj = ρl j<br />
/ S are <strong>in</strong>ertia coefficients, j = 1,<br />
,<br />
n<br />
j<br />
and n is the number of network branches (exclud<strong>in</strong>g the<br />
generator branch). T = diag T } , R = col R } and<br />
{ j<br />
{ j<br />
2<br />
Q = diag{<br />
Q Q } . (3)<br />
D<br />
Let nc<br />
denote the number of nodes. Then l = n − nc<br />
+ 1<br />
is the number of l<strong>in</strong>ks (exclud<strong>in</strong>g the generator branch)<br />
and n − l is the number of tree branches.<br />
Like an electrical network, a fluid network must satisfy<br />
Kirchhoff's current law, i.e., the flow out of any node is<br />
equal to the flow <strong>in</strong>to that node. Mathematically,<br />
Kirchhoff's current law for fluid flow networks can be<br />
expressed as:<br />
n<br />
∑<br />
j=<br />
1<br />
E<br />
Qij<br />
E<br />
⎡Q<strong>in</strong><br />
⎤<br />
⎢ ⎥ = 0<br />
⎣ Q ⎦<br />
Q<strong>in</strong><br />
or<br />
Q + e<br />
j<br />
Q<strong>in</strong>i<br />
Q<br />
<strong>in</strong><br />
j<br />
j<br />
= 0,<br />
i = 1, , n − l . (4)<br />
where n − l + 1 is the number of nodes (of which one is a<br />
“reference” node), Q is a vector of flows,<br />
E = [ e E ] Q<br />
, and E = E ] is a full rank matrix of<br />
Q<strong>in</strong><br />
Q<strong>in</strong><br />
Q<br />
[<br />
Qij<br />
order ( n − l)<br />
× n where E = 1 if branch j is connected<br />
Qij<br />
to node i and the flow goes away from node i ,<br />
E = −1 if it goes <strong>in</strong>to node i , E = 0 if branch j is<br />
Qij<br />
not connected to node i ; e<br />
Q<strong>in</strong><br />
is an (n-l)×1 vector such<br />
that, if the generator is connected to node i and the flow<br />
goes away from node i then e = 1 , if the flow goes<br />
Q<strong>in</strong>i<br />
Qij<br />
© 2013 ACADEMY PUBLISHER