25.01.2015 Views

Download Full Issue in PDF - Academy Publisher

Download Full Issue in PDF - Academy Publisher

Download Full Issue in PDF - Academy Publisher

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!