CN102525431A - Cardiovascular physiology signal detection device and method - Google Patents
Cardiovascular physiology signal detection device and method Download PDFInfo
- Publication number
- CN102525431A CN102525431A CN2010105828034A CN201010582803A CN102525431A CN 102525431 A CN102525431 A CN 102525431A CN 2010105828034 A CN2010105828034 A CN 2010105828034A CN 201010582803 A CN201010582803 A CN 201010582803A CN 102525431 A CN102525431 A CN 102525431A
- Authority
- CN
- China
- Prior art keywords
- signal
- cardiovascular physiology
- pulse
- self
- adaptive processing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
The invention provides a cardiovascular physiology signal detection device and a method, wherein the cardiovascular physiology signal detection device comprises an acquisition unit and a self-adaptive processing unit, the acquisition unit is used for collecting cardiovascular physiology signals, and the self-adaptive processing unit is used for treating the cardiovascular physiology signals in self adapting to obtain the cardiovascular physiology signals with maximum signal to interference and noise ratio. The cardiovascular physiology signal detection device in the embodiment obtains high-gain cardiovascular physiology signals by treating the cardiovascular physiology signals in self adapting, reduces influence of environment noise, and increases detection precision of the cardiovascular physiology signals. Users can complete detection of the cardiovascular physiology signals under unrestrained situation.
Description
Technical field
The present invention relates to technical field of medical, particularly, relate to a kind of checkout gear and method of cardiovascular physiology signal.
Background technology
Symptoms such as cardiopalmus, dizziness, chest pain and dyspnea can appear in some heart rate unusual (atrial fibrillation) patient regularly; But also symptom seldom or never appears in some patient; So only rely on subjective sensation, therefore the patient usually can not get diagnosis in time and treatment, thereby has increased the apoplexy risk.
In order to promote the awareness of patient for heart rate unusual (atrial fibrillation) and apoplexy risk thereof; For the patient can in time be diagnosed and treat and reduce the apoplexy risk a kind of objective guarantee means are provided; Therefore need in time to find the variation of cardiovascular physiology signals such as cardiac rhythm and pulse wave velocily, as rapid heart rate, slow excessively, uneven, pulse deficit (frequency of pulse is less than the heartbeat number of times) appear suddenly, the heartbeat power does not wait and disease such as blood pressure.Therefore need a kind of subjective sensation of not dependent patient, can in daily routines, in time detect the device that cardiovascular physiology signals such as changes in heart rate, pulse wave velocily and blood pressure change as big as you please.
Yet; The checkout gear of cardiovascular physiology signals such as the conventional sense heart rate that uses in the prior art, pulse wave, blood pressure; Usually need checkout gear be laid, clamping or be worn on detected person's the health, cause the uncomfortable or psychological anxiety on the user's body easily, and detected cardiovascular physiology signal receives the influence of environment and noise signal is too much easily; Make detected cardiovascular physiology signal inaccurate, cause user's the state of an illness to incur loss through delay easily.
Summary of the invention
For addressing the above problem, the present invention provides a kind of checkout gear and method of cardiovascular physiology signal, is used for solving prior art and detects the inaccurate problem of cardiovascular physiology signal.
For this reason, the present invention provides a kind of checkout gear of cardiovascular physiology signal, wherein, comprising: collecting unit and self-adaptive processing unit;
Said collecting unit is used to gather the cardiovascular physiology signal;
Said self-adaptive processing unit is connected with said collecting unit, is used for said cardiovascular physiology signal is carried out self-adaptive processing, obtains the maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio;
Wherein, also comprise analytic unit and output unit;
Said analytic unit is connected with said self-adaptive processing unit, is used for according to the maximum cardiovascular physiology signal of said Signal to Interference plus Noise Ratio the cardiovascular physiology parameter being made analytical judgment;
Said output unit is connected with said analytic unit, is used for the judged result output with said analytic unit.
Wherein, said collecting unit comprises: cardiac signal acquisition module and at least one pulse signal acquisition module;
Said cardiac signal acquisition module is used for gathering the cardiac signal of cardiovascular physiology signal, and said cardiac signal acquisition module comprises at least two cardiac signal harvesters;
Said pulse signal acquisition module is used for gathering the pulse signal of cardiovascular physiology signal, and each said pulse signal acquisition module comprises at least two pulse signal harvesters.
Wherein, comprise the analog digital conversion passage in the said cardiac signal harvester, said analog to digital conversion passage is used for converting said cardiac signal into digital signal;
Comprise the analog digital conversion passage in the said pulse signal harvester, said analog to digital conversion passage is used for converting said pulse signal into digital signal.
Wherein, said self-adaptive processing unit comprises: heartbeat self-adaptive processing module and pulse self-adaptive processing module;
Said heartbeat self-adaptive processing module is used for said cardiac signal is carried out self-adaptive processing;
Said pulse self-adaptive processing module is used for said pulse signal is carried out self-adaptive processing.
Wherein, said heartbeat self-adaptive processing module comprises: phase place adjustment submodule, equilibrium treatment submodule and signal synthon module;
Said equilibrium treatment submodule comprises the filtering channel group, is used for the cardiac signal of said cardiac signal acquisition module collection is carried out matched filtering with the filter factor of being set by said phase adjusting module output; The input of said phase place adjustment submodule is connected with outfan with the input of the cardiac signal of said bank of filters respectively; Be used to calculate the filter factor of said cardiac signal, said equalization filter group is carried out matched filtering according to said filter factor to said cardiac signal;
Said signal synthon module is used for obtaining the maximum cardiac signal of Signal to Interference plus Noise Ratio with carrying out anded through the said cardiac signal after the matched filtering.
Wherein, said phase place adjustment submodule comprises sef-adapting filter and channel selecting portion;
Wherein, said pulse self-adaptive processing module comprises: parameter estimation algorithm submodule and spatially adaptive filtering submodule;
Said parameter estimation algorithm submodule is used to calculate said pulse signal conduction delay time estimated value and filter weights;
Said spatially adaptive filtering submodule carries out weighted sum according to said filter weights to said pulse signal and calculates, and obtains the maximum pulse signal of Signal to Interference plus Noise Ratio.
Wherein, said parameter estimation algorithm submodule comprises: calculation channel selecting portion, passage gating portion and algorithm for estimating portion;
Said calculation channel selecting portion is used for selecting the calculation passage;
Said algorithm for estimating portion is used to calculate pulse signal conduction delay time estimated value and filter weights.
Wherein, the cardiovascular physiology parameter of said analytic unit analysis comprises at least a in cardiac signal intensity, pulse signal intensity, pulse wave conduction velocity and the blood pressure value.
Wherein, said collecting unit is the cardiac signal acquisition module, and said cardiac signal acquisition module is used for gathering the cardiac signal of cardiovascular physiology signal, and said cardiac signal acquisition module comprises at least two cardiac signal harvesters.
Wherein, said collecting unit is made up of at least one pulse signal acquisition module, and said pulse signal acquisition module is used for gathering the pulse signal of cardiovascular physiology signal, and each said pulse signal acquisition module comprises at least two pulse signal harvesters.
The present invention also provides a kind of detection method of cardiovascular physiology signal, wherein, comprising:
Collecting unit is gathered the cardiovascular physiology signal;
The self-adaptive processing unit carries out self-adaptive processing to said cardiovascular physiology signal, obtains the maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio.
Wherein, analytic unit is made analytical judgment according to the maximum cardiovascular physiology signal of said Signal to Interference plus Noise Ratio to the cardiovascular physiology parameter;
Output unit is with the judged result output of said analytic unit.
The present invention has following beneficial effect:
The checkout gear of embodiment of the invention cardiovascular physiology signal; For example cardiac signal, pulse signal etc. carry out self-adaptive processing to the cardiovascular physiology signal through the self-adaptive processing unit; Obtain the cardiovascular physiology signal of high-gain; Reduce influence of environmental noise, capacity of resisting disturbance is strong, has improved cardiovascular physiology signal detection precision; The user can be in bed or seat detect, guaranteeing when the user detects can be in the detection of loosening, accomplish as big as you please under the situation of no psychology or physiological load the cardiovascular physiology signal.
The detection method of embodiment of the invention cardiovascular physiology signal; Through for example cardiac signal, pulse signal etc. carry out self-adaptive processing to the cardiovascular physiology signal; Obtain the cardiovascular physiology signal of high-gain; Improved cardiovascular physiology signal detection precision, the user can be in bed or seat detect, guaranteeing when the user detects can be in the detection of loosening, accomplish under the situation of no psychology or physiological load the cardiovascular physiology signal.
Description of drawings
Fig. 1 is the structural representation of checkout gear first embodiment of cardiovascular physiology signal provided by the invention;
Fig. 2 is the structural representation of checkout gear second embodiment of cardiovascular physiology signal provided by the invention;
Fig. 3 is the structural representation of checkout gear the 3rd embodiment of cardiovascular physiology signal provided by the invention;
Fig. 4 is the structural representation of the heartbeat self-adaptive processing module among Fig. 3;
Fig. 5 is the equilibrium treatment sub modular structure sketch map among Fig. 4;
Fig. 6 is the structural representation of the phase place synthesis module among Fig. 4;
Fig. 7 is the workflow diagram of heartbeat self-adaptive processing module shown in Figure 4;
Fig. 8 is the structural representation of the pulse self-adaptive processing module among Fig. 3;
Fig. 9 is the workflow diagram of pulse self-adaptive processing module shown in Figure 8;
Figure 10 is the structural representation of checkout gear the 4th embodiment of cardiovascular physiology signal provided by the invention;
Figure 11 is the structural representation of checkout gear the 5th embodiment of cardiovascular physiology signal provided by the invention;
Figure 12 is the flow chart of detection method first embodiment of cardiovascular physiology signal provided by the invention;
Figure 13 is the flow chart of detection method second embodiment of cardiovascular physiology signal provided by the invention.
The specific embodiment
For making those skilled in the art understand technical scheme of the present invention better, be described in detail below in conjunction with the checkout gear and the method for accompanying drawing to cardiovascular physiology signal provided by the invention.
Fig. 1 is the structural representation of checkout gear first embodiment of cardiovascular physiology signal provided by the invention.As shown in Figure 1; The checkout gear of embodiment of the invention cardiovascular physiology signal comprises: collecting unit 10 and self-adaptive processing unit 20, and wherein, collecting unit 10 is used to gather the cardiovascular physiology signal; The cardiovascular physiology signal comprises pulse, heartbeat etc.; Self-adaptive processing unit 20 is connected with collecting unit 10, and the cardiovascular physiology signal that 20 pairs of collecting units in self-adaptive processing unit 10 are gathered carries out self-adaptive processing, and self-adaptive processing helps to reduce influence of environmental noise; It is strong to improve capacity of resisting disturbance, obtains the maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio.
Fig. 2 is the structural representation of checkout gear second embodiment of cardiovascular physiology signal provided by the invention.As shown in Figure 2; On the basis of first embodiment, the checkout gear of embodiment of the invention cardiovascular physiology signal also comprises analytic unit 30 and output unit 40, wherein; The maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio that analytic unit 30 obtains according to self-adaptive processing unit 20; Cardiovascular physiology parameter and health are made analytical judgment, and whether the cardiovascular physiology parameter of analyzing human body is unusual, and whether as arrhythmia, the heart quiver etc.; Can also analyze detected person's sleep state, for example: symptoms such as sleep disorder, apnea disease.Output unit 40 is for example exported to user side through digital terminals such as mobile phone, computer monitor or printers with the analysis result output of analytic unit 30.
Fig. 3 is the structural representation of checkout gear the 3rd embodiment of cardiovascular physiology signal provided by the invention.As shown in Figure 3; In embodiments of the present invention, collecting unit 10 comprises cardiac signal acquisition module 11 and at least one pulse signal acquisition module 12, wherein; Cardiac signal acquisition module 11 is used for gathering the cardiac signal of cardiovascular physiology signal; Comprise at least two cardiac signal harvesters 111 in the cardiac signal acquisition module 11, pulse signal acquisition module 12 is used for gathering the pulse signal of cardiovascular physiology signal, comprises at least two pulse signal harvesters 121 in the pulse signal acquisition module 12; In embodiments of the present invention; Comprise N cardiac signal harvester 111 in the cardiac signal acquisition module 11, the quantity of pulse signal acquisition module 12 is D (D＞0), comprises M pulse signal harvester 121 in each pulse signal acquisition module 12.All include the analog digital conversion passage in cardiac signal harvester 111 and the pulse signal harvester 121; The analog digital conversion passage is used for cardiac signal and pulse signal are carried out amplification filtering and analog digital conversion; Convert cardiac signal and pulse signal into can carry out digital processing cardiac signal and pulse signal respectively, cardiac signal harvester 111 can adopt optical fiber type harvester, piezoelectric type harvester, electret condenser type harvester, gasbag pressure formula harvester, MEMS (Micro-Electro-Mechanical SystemsMEMS) acceleration formula harvester or ultrasound wave harvester etc. with pulse signal harvester 121.
When actual acquisition cardiovascular physiology signal such as heartbeat signal or pulse signal; Can N cardiac signal harvester 111 and the flow direction of M*D pulse signal harvester 121 according to blood be arranged in order, with user's cardia contact with cardiac signal harvester 111, with user's heart in addition the tremulous pulse at position contact with pulse signal harvester 121 and get final product.In embodiments of the present invention; N cardiac signal harvester 111 and M*D pulse signal harvester 121 are arranged in order laying in bed or on the seat according to the flow direction of blood in the descending large artery trunks; The user can lie on a bed or be sitting on the seat; Only need user's cardia be contacted with cardiac signal harvester 111 with the following position of user's heart is contacted with pulse signal harvester 121 and get final product; Need the cardiac signal harvester 111 in the collecting unit or pulse signal harvester 121 not worn or be clamped on the user's body; Guarantee that the user can not have the psychology burden, detect as big as you please, helps obtaining user's the cardiovascular physiology of cardiac signal or pulse signal etc. more accurately signal under the state that loosens.In practical application; The signals collecting frequency f can be set to 45000 times/second; The spacing distance that cardiac signal harvester 111 or pulse signal harvester 121 are arranged can be set to 10mm; Cardiac signal harvester 111 or pulse signal harvester 121 also can not be equidistant arrangements, as long as spread length is respectively less than the wavelength of cardiac signal or pulse signal.
Self-adaptive processing unit 20 in the checkout gear of cardiovascular physiology signal shown in Figure 3 comprises heartbeat self-adaptive processing module 21 and pulse self-adaptive processing module 22; Heartbeat self-adaptive processing module 21 is connected with cardiac signal acquisition module 11; The cardiac signal that 21 pairs of cardiac signal acquisition modules of heartbeat self-adaptive processing module 11 are gathered carries out self-adaptive processing; Obtain the maximum cardiac signal of Signal to Interference plus Noise Ratio; And the analytic unit 30 that sends to of the cardiac signal that Signal to Interference plus Noise Ratio is maximum, pulse self-adaptive processing module 22 is connected with D pulse signal acquisition module 12, and pulse self-adaptive processing module 22 is used for the pulse signal of D pulse signal acquisition module collection is carried out self-adaptive processing; Obtain the maximum pulse signal of Signal to Interference plus Noise Ratio, and the analytic unit 30 that sends to of the pulse signal that Signal to Interference plus Noise Ratio is maximum.
Fig. 4 is the structural representation of the heartbeat self-adaptive processing module among Fig. 3.As shown in Figure 4 and combine Fig. 3; Heartbeat self-adaptive processing module 21 comprises equilibrium treatment submodule 211, phase place adjustment submodule 212 and phase place synthesis module 213; Wherein, Phase place adjustment submodule 212 is used to set the filter factor of cardiac signal; The input of phase place adjustment submodule 212 is connected to the input and the outfan of the cardiac signal of equilibrium treatment submodule 211, and the outfan of phase place adjustment submodule 212 is connected with the filter factor input of equilibrium treatment submodule 211, and equilibrium treatment submodule 211 matees cardiac signal according to the filter factor that phase place adjustment submodule 212 obtains; Signal synthon module 213 is connected with equilibrium treatment submodule 211, and the cardiac signal of 213 pairs of equilibrium treatment submodules of signal synthon module, 211 inputs synthesizes anded to obtain the maximum cardiac signal of Signal to Interference plus Noise Ratio.
Fig. 5 is the equilibrium treatment sub modular structure sketch map among Fig. 4.As shown in Figure 5 and combine Fig. 3; Equilibrium treatment submodule 211 comprises that impulse response length is the N road filtering channel of L; Filtering channel is made up of adder 2111 grades; The input of N road filtering channel is connected respectively to N cardiac signal harvester 111 in the cardiac signal acquisition module 11, is used to receive cardiac signal X
1～X
N, the outfan Y of N road filtering channel
1～Y
NThe cardiac signal that will pass through equilibrium treatment outputs to signal synthon module 213.
Fig. 6 is the structural representation of the phase place synthesis module among Fig. 4.As shown in Figure 6 and combine Fig. 5; Phase place adjustment submodule 212 comprises sef-adapting filter 2121 and channel selecting portion 2122; Wherein, sef-adapting filter 2121 can adopt adaptive filter methods such as Wei Na (Wiener) optimal filter method, recurrence least square (RLS) adaptive filter method.In the embodiment of the invention, sef-adapting filter 2121 adopts lowest mean square, and (Least MeanSquade, LMS) adaptive filter method constitute lowest mean square (Least Mean Squade, LMS) sef-adapting filter.The convergence coefficient μ of LMS sef-adapting filter satisfies the condition of convergence 0＜μ＜2/R=λ
Max, R=λ wherein
MaxAutocorrelation matrix eigenvalue of maximum for input cardiac signal Xj.
Fig. 7 is the workflow diagram of heartbeat self-adaptive processing module shown in Figure 4.As shown in Figure 7 and combine Fig. 5 and Fig. 6, the workflow of the heartbeat self-adaptive processing module in the embodiment of the invention specifically comprises following job step:
In an embodiment, the impulse response length of N road filtering channel is L in the equilibrium treatment submodule 211, and phase place adjustment submodule 212 is set the initial filter coefficients of N road filtering channel, and the initial filter coefficients of N road filtering channel is: W
11～W
1L... W
N1～W
NL
In embodiments of the present invention, the channel selecting submodule comes the 2122 i road filtering channel filtering channels as a reference selected in the equilibrium treatment submodules 211, and the filter factor of i road filtering channel is W
I1～W
IL, the filtering signal of the input and output of i road filtering channel is X
iAnd Y
i, get into step 603 then
In this step, if i＞N then selects the 1 road filtering channel filtering channel as a reference.
Channel selecting portion 2122 selects j road filtering channel as the adjustment filtering channel, and the filtering signal of the input and output of j road filtering channel is X
jAnd Y
j, wherein,,, then, get into step 604 if j＞N then selects the 1 road filtering channel as the adjustment filtering channel if j=i then selects j+1 bar filtering channel as the adjustment filtering channel.
Whether the deviate of the cardiac signal of step 604, the cardiac signal of judging the output of reference filtering passage and the output of adjustment filtering channel restrains.
The cardiac signal Y of reference filtering passage output
iCardiac signal Y with the output of adjustment filtering channel
jDeviate e, preset maximum deflection difference value is ε, if e＜ε explains convergence, to be ε disperse the astringent empirical value according in the past deviate e (n) to preset maximum deflection difference value obtains, and gets into step 605 then, otherwise, then get into step 602.
The filter factor of step 605, calculating adjustment filtering channel.
In embodiments of the present invention; Sef-adapting filter 2121 adopts the LMS sef-adapting filter; The cardiac signal Yj of the cardiac signal Xi of the adjustment filtering channel input that 2121 pairs of sef-adapting filters receive and the cardiac signal Yi of output and the output of reference filtering passage carries out time-domain adaptive and handles the filter factor W of the adjustment filtering channel that obtains
J1～W
JL, get into step 606 then.
The filter factor of step 606, renewal adjustment filtering channel.
The filter factor W of the j road filtering channel of the conduct adjustment filtering channel that calculates in the step 605
J1～W
JLBe updated in the filtering channel of j road, get into step 607 then.
Select j+1 bar filtering channel as the adjustment filtering channel, circulation is carried out the time-domain adaptive of step 603-606 and is handled.
The cardiac signal of 213 pairs of equilibrium treatment submodules of signal synthon module, 211 inputs synthesizes anded to obtain the maximum cardiac signal of Signal to Interference plus Noise Ratio.
The workflow that it is pointed out that heartbeat self-adaptive processing module also can realize through the mode of other adaptive-filtering, no longer the mode of other adaptive-filtering is discussed one by one at this.
Fig. 8 is the structural representation of the pulse self-adaptive processing module among Fig. 3.As shown in Figure 8 and combine Fig. 3; Pulse self-adaptive processing module 22 comprises: parameter estimation algorithm submodule 221 and spatially adaptive filtering submodule 222; Wherein, The filter weights of parameter spectrum algorithm for estimating submodule 221 calculating pulse signals and pulse conduction velocity (pulse wave averagevelocity, PWV), the filter weights that spatially adaptive filtering submodule 222 calculates according to parameter estimation algorithm submodule 221 is carried out weighted sum to pulse signal and is calculated; Obtain the maximum pulse signal of Signal to Interference plus Noise Ratio, then that Signal to Interference plus Noise Ratio is maximum pulse signal and pulse conduction velocity output to analytic unit 30.Parameter estimation algorithm submodule 221 comprises: calculation channel selecting portion 2211, passage gating portion 2212 and algorithm for estimating portion 2213; Wherein, Calculation channel selecting portion 2211 is used for Cyclic Selection calculation passage; Algorithm for estimating portion 2213 is used to calculate the filter weights and the pulse wave conduction velocity of pulse signal, and the pulse wave conduction velocity of d road pulse signal is PWVd, and passage gating portion 2212 is used for the gating transmission channel.Algorithm for estimating portion 2213 can adopt least mean-square error (Minimum Mean Squared Error; MMSE) method of estimation, constant modulus algorithm (Constant Modulus Algorithm; CMA) multiple signal classification (the Multi SignalClassification of method of estimation, super-resolution spectrum algorithm for estimating; MUSIC) method of estimation and invariable rotary technology (Estimating signalparameter via rotational invariance techniques, ESPRIT) method for parameter estimation of spatially adaptive filtering such as method of estimation.Spatially adaptive filtering submodule 222 comprises D road filtering channel, and D road filtering channel is connected respectively to D pulse signal acquisition module, and each road filtering channel has M road input, is connected respectively to M pulse signal harvester 121.
Fig. 9 is the workflow diagram of pulse self-adaptive processing module shown in Figure 8.As shown in Figure 9, the workflow of pulse self-adaptive processing module specifically comprises the steps:
Step 801, calculating pulse signal conduction delay time estimated value.
In embodiments of the present invention, utilize multiple signal classification (Multi Signal Classification, MUSIC) the computational methods calculating pulse signal conduction delay time estimated value of super-resolution spectrum algorithm for estimating through algorithm for estimating portion 2213.At first the mid frequency from pulse signal extracts the needed narrow band signal of MUSIC computational methods, and the wavelength of the mid frequency of this pulse signal is greater than spread length between two pulse signal harvesters 121 when detecting the user.The d road filtering channel in the filtering channel of D road is selected to connect by calculation channel selecting portion 2211, the pulse signal Y that d the pulse signal acquisition module 12 that 2213 receptions of algorithm for estimating portion are transmitted through d road filtering channel collects
D1～Y
DM, extract M pulse signal Y in d the pulse signal acquisition module 12 then
D1～Y
DMMid frequency under narrow band signal Q and I value, generate the signal phasor X (t) of M dimension, pulse signal vector X and covariance matrix R thereof are respectively shown in formula (1) and formula (2):
X(t)＝[X
1(t)，...X
m(t)，...X
M(t)]
T (1)
R＝E[X(t)X(t)
H] (2)
Wherein, T representes the vector transposition, and H representes complex-conjugate transpose, the element X among the signal phasor X (t)
m(t) imaginary part and the real part of (1≤m≤M, m are natural number) are respectively Q and I, X
m(t) be m pulse signal Y
DmComplex values, t is the sampling instant at interval of pulse signal certain hour.
In practical application, the maximum of covariance matrix R is intended right function shown in formula (3),
χ in the formula (3)=[X (1) ... X (T)], χ is the vector matrix of pulse signal, and T is the hop count of the pulse signal that collects, and T can be 1280.
Calculate the maximum of covariance matrix R and intend right function R
XxCharacteristic value λ
1～λ
M(wherein, λ
1>=λ
2>=...>=λ
M), according to the R that surpasses preset noise power
XxCharacteristic value pulse signal incidence wave quantity is set, calculate corresponding to characteristic value λ
1～λ
MEigenvector e
1～e
MObtain noise eigenvectors E then
N, E
NBe the corresponding eigenvector of (M-G) that be lower than noise power individual characteristic value, G is the incidence wave quantity of pulse signal, i.e. signal source quantity, wherein, E
NShown in formula (4),
E
N＝(e
G+1，...e
M) (4)
According to the mathematical model in far field, arrowband, can the expression of the signal phasor X (t) shown in the formula (1) be represented with formula (5),
X(t)＝AF(t)+N(t) (5)
Wherein, the noise signal vector of N (t) expression M dimension, the spacing wave vector of F (t) expression G dimension, A is the steering vector matrix of M*G dimension, A is shown in formula (6):
A＝[a(τ
1)，a(τ
2)，...，a(τ
G)] (6)
Wherein, steering vector is shown in formula (7):
A(τ
i)＝[exp(-jω
0τ
i2)，...，exp(-jω
0τ
iM)]T (7)
Wherein, ω
0=2 π f, f are the frequency of pulse signal, the delay time T of pulse signal harvester 121
iFor shown in formula (8):
Wherein, j is the serial number of pulse signal harvester, s
jBe j the distance between signal picker and first pulse signal harvester of fighting of feeling the pulse, i is by the detected pulse signal source sequence of pulse signal harvester number, V
iBe the signaling rate in pulse signal source, θ is the direction of pulse signal, and in the present embodiment, the orientation of pulse signal harvester is identical with the direction of propagation of pulse signal, and getting θ is 90 °, and formula (8) can use formula (9) and (10) to represent:
τ
i＝[(s
2+…+s
j)/((j-1)*V
i)] (10)
Can obtain evaluation function P according to formula (4) and (7)
MU(τ), evaluation function P
MU(τ) shown in formula (11),
Wherein, a (τ) expression pulse signal harvester is about the complex response of delay time T.
Find out evaluation function P according to formula (11)
MUBe exactly the pulse signal conduction delay time estimated value τ of adjacent pulse signal harvester the time delay that maximum point (τ) is corresponding
P, get into step 802 then.
Step 802, the filter weights of calculating pulse signal.
The computing formula of filter factor shown in formula (12),
W
H＝(A
HA)
-1A
H≡[W
1，W
2，...，W
M]
H (12)
Then according to formula (6), formula (7) and formula (12) and pulse signal conduction delay time estimated value τ
PCalculate the filter weights W of d road filtering channel
K1～W
KM, get into step 803 then.
Step 803, obtain the maximum pulse signal of Signal to Interference plus Noise Ratio.
The pulse signal Y that collects according to the pulse signal acquisition module
1～Y
MAnd filter weights W
K1～W
KMCarry out the complex conjugate weighted sum computing of multiplying each other, obtain the maximum pulse signal of Signal to Interference plus Noise Ratio, then that Signal to Interference plus Noise Ratio is maximum pulse signal is transported to analytic unit.
Checkout gear with cardiovascular physiology signal shown in Figure 3 is an example; Receive maximum cardiac signal and the pulse signal of Signal to Interference plus Noise Ratio that obtains like Fig. 7 with like Fig. 9 when analytic unit 30; The peak value number of times that extracts interior cardiac signal of unit interval and pulse signal appearance is respectively as cardiac signal number of times and pulse signal number of times; Calculate cardiac signal and pulse signal signal intensity in the unit interval and and extract the pulse wave conduction velocity PWV of d road pulse signal as cardiac signal intensity and pulse signal intensity
d, then according to PWV calculate cardiac signal and pulse signal average pulse wave conduction velocity (pulse wave average velocity, PWAV); Unit is m/s, calculates the difference of PWV and PWAV, calculates blood pressure value (Press Diff; PD); Unit is mmHg, wherein, and shown in average pulse wave conduction velocity PWAV and blood pressure value PD difference formula (13) and the formula (14):
PWAV(t+1)＝β*PWAV(t)+(1-β)*PWV(t)
PD＝γ*(PWV(t)-PWAV(t)
Wherein, β is the rolling average coefficient (0.01＜β＜0.98) of pulse wave conduction velocity; γ changes proportionality coefficient (10＜γ＜60) for the pulse wave conduction velocity changes with pressure value, and t detects constantly for the pulse wave conduction velocity, carries out the detection of pulse wave signal usually at a certain time interval.
Analytic unit is made analytical judgment according to the maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio to cardiovascular physiology parameter and health.The embodiment of the invention is analyzed the detected person according to one or more cardiovascular physiology parameters in cardiac signal intensity, pulse signal intensity, pulse wave conduction velocity and the blood pressure value etc. whether cardiovascular disease is taken place, and for example heart rate is unusual, dysarteriotony etc.Then analysis result is sent to display terminal, display terminal can be digital terminals such as mobile phone, printer or computer.
For example cardiac signal, pulse signal etc. carry out self-adaptive processing to the embodiment of the invention to the cardiovascular physiology signal through the self-adaptive processing unit; Obtain the cardiovascular physiology signal of high-gain; Reduce influence of environmental noise, capacity of resisting disturbance is strong, has improved cardiovascular physiology signal detection precision; The user can be in bed or seat detect, guaranteeing when the user detects can be in the detection of loosening, accomplish under the situation of no psychology or physiological load the cardiovascular physiology signal.
Figure 10 is the structural representation of checkout gear the 4th embodiment of cardiovascular physiology signal provided by the invention.Shown in figure 10, the checkout gear of embodiment of the invention cardiovascular physiology signal can include only cardiac signal acquisition module 11 and heartbeat self-adaptive processing module 21, and wherein, the workflow of heartbeat self-adaptive processing module 21 is as shown in Figure 7, repeats no more at this.
Figure 11 is the structural representation of checkout gear the 5th embodiment of cardiovascular physiology signal provided by the invention.Shown in figure 11; The checkout gear of embodiment of the invention cardiovascular physiology signal also can include only pulse signal acquisition module 12 and pulse self-adaptive processing module 22; Wherein, concrete workflow such as Fig. 9 of pulse self-adaptive processing module 22 are said, repeat no more at this.
Figure 12 is the flow chart of detection method first embodiment of cardiovascular physiology signal provided by the invention.Shown in figure 12, the concrete workflow diagram of the detection method of embodiment of the invention cardiovascular physiology signal comprises the steps:
In embodiments of the present invention, adopt the checkout gear of cardiovascular physiology signal as shown in Figure 1, collecting unit 10 is used to gather the cardiovascular physiology signal, and the cardiovascular physiology signal comprises pulse, heartbeat etc., gets into step 122 then.
The cardiovascular physiology signal that 20 pairs of collecting units in self-adaptive processing unit 10 are gathered carries out self-adaptive processing, and self-adaptive processing helps to reduce influence of environmental noise, and it is strong to improve capacity of resisting disturbance, obtains the maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio.
Figure 13 is the flow chart of detection method second embodiment of cardiovascular physiology signal provided by the invention.Shown in figure 13, the concrete workflow diagram of the detection method of embodiment of the invention cardiovascular physiology signal comprises the steps:
In embodiments of the present invention, adopt the checkout gear of cardiovascular physiology signal as shown in Figure 2, collecting unit 10 is used to gather the cardiovascular physiology signal, and the cardiovascular physiology signal comprises pulse, heartbeat etc., gets into step 132 then.
The cardiovascular physiology signal that 20 pairs of collecting units in self-adaptive processing unit 10 are gathered carries out self-adaptive processing, obtains the maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio.
The embodiment of the invention is through for example cardiac signal, pulse signal etc. carry out self-adaptive processing to the cardiovascular physiology signal; Obtain the cardiovascular physiology signal of high-gain; Improved cardiovascular physiology signal detection precision; The user can be in bed or seat detect, guaranteeing when the user detects can be in the detection of loosening, accomplish under the situation of no psychology or physiological load the cardiovascular physiology signal.
It is understandable that above embodiment only is the illustrative embodiments that adopts for principle of the present invention is described, yet the present invention is not limited thereto.For the one of ordinary skilled in the art, under the situation that does not break away from spirit of the present invention and essence, can make various modification and improvement, these modification also are regarded as protection scope of the present invention with improving.
Claims (14)
1. the checkout gear of a cardiovascular physiology signal is characterized in that comprising: collecting unit and self-adaptive processing unit;
Said collecting unit is used to gather the cardiovascular physiology signal;
Said self-adaptive processing unit is connected with said collecting unit, is used for said cardiovascular physiology signal is carried out self-adaptive processing, obtains the maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio.
2. the checkout gear of cardiovascular physiology signal according to claim 1 is characterized in that, also comprises: analytic unit and output unit;
Said analytic unit is connected with said self-adaptive processing unit, is used for according to the maximum cardiovascular physiology signal of said Signal to Interference plus Noise Ratio the cardiovascular physiology parameter being made analytical judgment;
Said output unit is connected with said analytic unit, is used for the judged result output with said analytic unit.
3. the checkout gear of cardiovascular physiology signal according to claim 1 and 2 is characterized in that, said collecting unit comprises: cardiac signal acquisition module and at least one pulse signal acquisition module;
Said cardiac signal acquisition module is used for gathering the cardiac signal of cardiovascular physiology signal, and said cardiac signal acquisition module comprises at least two cardiac signal harvesters;
Said pulse signal acquisition module is used for gathering the pulse signal of cardiovascular physiology signal, and each said pulse signal acquisition module comprises at least two pulse signal harvesters.
4. the checkout gear of cardiovascular physiology signal according to claim 3 is characterized in that:
Comprise the analog digital conversion passage in the said cardiac signal harvester, said analog to digital conversion passage is used for converting said cardiac signal into digital signal;
Comprise the analog digital conversion passage in the said pulse signal harvester, said analog to digital conversion passage is used for converting said pulse signal into digital signal.
5. according to the checkout gear of claim 3 or 4 described cardiovascular physiology signals, it is characterized in that said self-adaptive processing unit comprises: heartbeat self-adaptive processing module and pulse self-adaptive processing module;
Said heartbeat self-adaptive processing module is used for said cardiac signal is carried out self-adaptive processing;
Said pulse self-adaptive processing module is used for said pulse signal is carried out self-adaptive processing.
6. the checkout gear of cardiovascular physiology signal according to claim 5 is characterized in that, said heartbeat self-adaptive processing module comprises: phase place adjustment submodule, equilibrium treatment submodule and signal synthon module;
Said equilibrium treatment submodule comprises the filtering channel group, is used for the cardiac signal of said cardiac signal acquisition module collection is carried out matched filtering with the filter factor of being set by said phase adjusting module output; The input of said phase place adjustment submodule is connected with outfan with the input of the cardiac signal of said bank of filters respectively; Be used to calculate the filter factor of said cardiac signal, said equalization filter group is carried out matched filtering according to said filter factor to said cardiac signal;
Said signal synthon module is used for obtaining the maximum cardiac signal of Signal to Interference plus Noise Ratio with carrying out anded through the said cardiac signal after the matched filtering.
7. the checkout gear of cardiovascular physiology signal according to claim 6 is characterized in that, said phase place adjustment submodule comprises sef-adapting filter and channel selecting portion.
8. the checkout gear of cardiovascular physiology signal according to claim 5 is characterized in that, said pulse self-adaptive processing module comprises: parameter estimation algorithm submodule and spatially adaptive filtering submodule;
Said parameter estimation algorithm submodule is used to calculate said pulse signal conduction delay time estimated value and filter weights;
Said spatially adaptive filtering submodule carries out weighted sum according to said filter weights to said pulse signal and calculates, and obtains the maximum pulse signal of Signal to Interference plus Noise Ratio.
9. the checkout gear of cardiovascular physiology signal according to claim 8 is characterized in that, said parameter estimation algorithm submodule comprises: calculation channel selecting portion, passage gating portion and algorithm for estimating portion;
Said calculation channel selecting portion is used for selecting the calculation passage;
Said algorithm for estimating portion is used to calculate pulse signal conduction delay time estimated value and filter weights.
10. the checkout gear of cardiovascular physiology signal according to claim 2; It is characterized in that the cardiovascular physiology parameter of said analytic unit analysis comprises at least a in cardiac signal intensity, pulse signal intensity, pulse wave conduction velocity and the blood pressure value.
11. the checkout gear of 2 described cardiovascular physiology signals as requested; It is characterized in that; Said collecting unit is the cardiac signal acquisition module; Said cardiac signal acquisition module is used for gathering the cardiac signal of cardiovascular physiology signal, and said cardiac signal acquisition module comprises at least two cardiac signal harvesters.
12. the checkout gear of 2 described cardiovascular physiology signals as requested; It is characterized in that; Said collecting unit is made up of at least one pulse signal acquisition module; Said pulse signal acquisition module is used for gathering the pulse signal of cardiovascular physiology signal, and each said pulse signal acquisition module comprises at least two pulse signal harvesters.
13. the detection method of a cardiovascular physiology signal is characterized in that, comprising:
Collecting unit is gathered the cardiovascular physiology signal;
The self-adaptive processing unit carries out self-adaptive processing to said cardiovascular physiology signal, obtains the maximum cardiovascular physiology signal of Signal to Interference plus Noise Ratio.
14. the detection method of cardiovascular physiology signal according to claim 13 is characterized in that, also comprises:
Analytic unit is made analytical judgment according to the maximum cardiovascular physiology signal of said Signal to Interference plus Noise Ratio to the cardiovascular physiology parameter;
Output unit is with the judged result output of said analytic unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010582803.4A CN102525431B (en) | 2010-12-10 | 2010-12-10 | Cardiovascular physiology signal detection device and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010582803.4A CN102525431B (en) | 2010-12-10 | 2010-12-10 | Cardiovascular physiology signal detection device and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102525431A true CN102525431A (en) | 2012-07-04 |
| CN102525431B CN102525431B (en) | 2014-01-08 |
Family
ID=46334405
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201010582803.4A Active CN102525431B (en) | 2010-12-10 | 2010-12-10 | Cardiovascular physiology signal detection device and method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102525431B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106333655A (en) * | 2015-07-09 | 2017-01-18 | 三星电子株式会社 | Apparatus And Method For Analyzing Living Body Information |
| CN108567417A (en) * | 2018-04-24 | 2018-09-25 | 深圳还是威健康科技有限公司 | A kind of cardiovascular monitoring method and system based on intelligent folding fan |
| CN110960199A (en) * | 2019-12-24 | 2020-04-07 | 中国人民解放军陆军军医大学第一附属医院 | System for double-variable measurement of arteriosclerosis degree |
| CN113257418A (en) * | 2021-03-29 | 2021-08-13 | 广州科克里特生命科技有限公司 | Risk detection system and method for low back pain |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1111121A (en) * | 1994-08-30 | 1995-11-08 | 中国科学院上海技术物理研究所 | Self-adaptation analytical method and apparatus for electrocardiac and pulse signal |
| CN1121798A (en) * | 1994-08-16 | 1996-05-08 | 北京工业大学 | Cardiovascular function dynamic parameter testing analysis method and apparatus |
| US6647293B2 (en) * | 2000-01-19 | 2003-11-11 | Biotronik Mess-und Therapiegeräte GmbH & Co. Ingenieurbüro Berlin | Method for predictively calculating a cardiac signal course and for controlling the delivery of a stimulation pulse by an implantable cardiac device |
| CN1665443A (en) * | 2001-09-13 | 2005-09-07 | 康曼德公司 | A signal processing method and device for signal-to-noise improvement |
| JP4120100B2 (en) * | 1999-07-09 | 2008-07-16 | オムロンヘルスケア株式会社 | Non-invasive continuous blood pressure estimation device and non-invasive continuous blood pressure prediction device |
| CN101732040A (en) * | 2009-12-24 | 2010-06-16 | 中国科学院力学研究所 | Non-invasive multipath pulse wave detection device, system and analytical system |
-
2010
- 2010-12-10 CN CN201010582803.4A patent/CN102525431B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1121798A (en) * | 1994-08-16 | 1996-05-08 | 北京工业大学 | Cardiovascular function dynamic parameter testing analysis method and apparatus |
| CN1111121A (en) * | 1994-08-30 | 1995-11-08 | 中国科学院上海技术物理研究所 | Self-adaptation analytical method and apparatus for electrocardiac and pulse signal |
| JP4120100B2 (en) * | 1999-07-09 | 2008-07-16 | オムロンヘルスケア株式会社 | Non-invasive continuous blood pressure estimation device and non-invasive continuous blood pressure prediction device |
| US6647293B2 (en) * | 2000-01-19 | 2003-11-11 | Biotronik Mess-und Therapiegeräte GmbH & Co. Ingenieurbüro Berlin | Method for predictively calculating a cardiac signal course and for controlling the delivery of a stimulation pulse by an implantable cardiac device |
| CN1665443A (en) * | 2001-09-13 | 2005-09-07 | 康曼德公司 | A signal processing method and device for signal-to-noise improvement |
| CN101732040A (en) * | 2009-12-24 | 2010-06-16 | 中国科学院力学研究所 | Non-invasive multipath pulse wave detection device, system and analytical system |
Non-Patent Citations (1)
| Title |
|---|
| 张昕等: "《用自适应AR模型提高血氧饱和度测量中的脉搏波信号的检出率》", 《生物医学工程学杂志》, vol. 17, no. 3, 31 December 2000 (2000-12-31), pages 285 - 287 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106333655A (en) * | 2015-07-09 | 2017-01-18 | 三星电子株式会社 | Apparatus And Method For Analyzing Living Body Information |
| CN106333655B (en) * | 2015-07-09 | 2021-06-15 | 三星电子株式会社 | Apparatus and method for analyzing living body information |
| CN108567417A (en) * | 2018-04-24 | 2018-09-25 | 深圳还是威健康科技有限公司 | A kind of cardiovascular monitoring method and system based on intelligent folding fan |
| CN110960199A (en) * | 2019-12-24 | 2020-04-07 | 中国人民解放军陆军军医大学第一附属医院 | System for double-variable measurement of arteriosclerosis degree |
| CN110960199B (en) * | 2019-12-24 | 2022-05-27 | 中国人民解放军陆军军医大学第一附属医院 | System for double-variable measurement of arteriosclerosis degree |
| CN113257418A (en) * | 2021-03-29 | 2021-08-13 | 广州科克里特生命科技有限公司 | Risk detection system and method for low back pain |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102525431B (en) | 2014-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Venkatesan et al. | A novel LMS algorithm for ECG signal preprocessing and KNN classifier based abnormality detection | |
| CN101658425B (en) | Device and method for detecting attention focusing degree based on analysis of heart rate variability | |
| CN110236518B (en) | Electrocardio and heart-shock signal combined classification method and device based on neural network | |
| CN102488503B (en) | Continuous blood pressure measurer | |
| CN110179643A (en) | A kind of neck rehabilitation training system and training method based on annulus sensor | |
| CN114748080B (en) | Method and system for detecting and quantifying sensory-motor function | |
| CN103038772A (en) | Method of predicting the survivability of a patient | |
| CN103610447B (en) | A kind of Mental Workload online test method based on forehead EEG signal | |
| Barros et al. | Filtering noncorrelated noise in impedance cardiography | |
| CN102525431B (en) | Cardiovascular physiology signal detection device and method | |
| KR102075503B1 (en) | System of Predicting Dementia and Operating Method The Same | |
| CN106333652B (en) | A kind of sleep state analysis method | |
| CN104173043A (en) | Electrocardiogram (ECG) data analysis method suitable for mobile platform | |
| CN108742613A (en) | Orient coupling analytical method between the flesh of coherence partially based on transfer entropy and broad sense | |
| CN109871808A (en) | Atrial fibrillation model training and detecting method and device | |
| WO2023093770A1 (en) | Millimeter-wave radar-based noncontact electrocardiogram monitoring method | |
| CN106137167A (en) | A kind of motion artifacts detection method based on photoplethysmographic signal | |
| CN107582040B (en) | Method and device for monitoring heart rhythm | |
| CN109288515A (en) | Periodical monitoring method and device based on premature beat signal in wearable ECG signal | |
| WO2005000108A2 (en) | Radiation stress non-invasive blood pressure method | |
| Zhang et al. | Method of diagnosing heart disease based on deep learning ECG signal | |
| Kelwade et al. | Prediction of heart abnormalities using particle swarm optimization in radial basis function neural network | |
| CN113040738B (en) | Blood pressure detecting device | |
| Anuradha et al. | Classification of cardiac signals using time domain methods | |
| CN116226624A (en) | Channel selection method of motor imagery brain-computer interface based on tensor decomposition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |