CN105433931A - Processing device of light volume change description waveform and method thereof - Google Patents

Processing device of light volume change description waveform and method thereof Download PDF

Info

Publication number
CN105433931A
CN105433931A CN201510565554.0A CN201510565554A CN105433931A CN 105433931 A CN105433931 A CN 105433931A CN 201510565554 A CN201510565554 A CN 201510565554A CN 105433931 A CN105433931 A CN 105433931A
Authority
CN
China
Prior art keywords
volume
signal
digital filtering
describe
smooth change
Prior art date
Application number
CN201510565554.0A
Other languages
Chinese (zh)
Inventor
庄政达
罗展鹏
潘德城
Original Assignee
义明科技股份有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Priority to US201462051957P priority Critical
Priority to US62/051,957 priority
Priority to TW104125755 priority
Priority to TW104125755A priority patent/TWI558375B/en
Application filed by 义明科技股份有限公司 filed Critical 义明科技股份有限公司
Publication of CN105433931A publication Critical patent/CN105433931A/en

Links

Abstract

The present invention discloses a processing device of light volume change description waveform and a method thereof. The processing device comprises a light volume change description waveform (PPG) signal capture unit captures a PPG signal, and during capturing the PPG signal, an operation information capture unit captures motion information, and the PPG signal and the motion information are output to a first adaptive digital filtering unit to remove a motion artifact interference signal from the PPG signal. The PPG signal after removing is output to a second adaptive digital filtering unit to capture the period signals in the second PPG signal and improve the signal to noise ratio of the period signals to be taken as the basis of reading out accurate heartbeat signals. Therefore, the PPG signal processing device of the invention is applicable to a wearable heartbeat detection device for provide correct heartbeat information.

Description

Light change in volume describes blood processor and the method thereof of waveform

Technical field

The present invention describes blood processor and the method thereof of waveform about a kind of smooth change in volume, and espespecially a kind of light change in volume being applied to Wearable heartbeat detection device describes blood processor and the method thereof of waveform.

Background technology

Wearable electronic installation provides monitoring user motor function (as step function) usually, or provides the function such as detection and record human body physiological characteristics (as heart beating, blood pressure or blood oxygen), therefore quite likes by consumer.But, still more general fixed sphygomanometer detection heart beating result is poor for the heartbeat detection result of Wearable electronic installation, main cause is that this Wearable electronic installation is in time capturing the physiological feature signal of human body, and the limb activity of user can produce interfering signal to this physiological feature signal.

Generally the physiological feature signal captured can be converted to frequency-domain waveform (as shown in Figure 9 B) from time domain waveform (as shown in Figure 9 A), again using the periodic signal in frequency domain as the foundation measuring palmic rate, as measured correct palmic rate namely as shown in fig. 9d.But, but in frequency-region signal except heartbeat signal S heartoutward, as shown in Figure 9 B, also include because hand swings the low frequency signal caused with body-sway motion, or comprise and rock caused high-frequency signal because hand swings fast fast with health, when wherein the low frequency of part or the energy of high-frequency signal can higher than the energy of heartbeat signal.Owing to including motion disturbance signals S in frequency signal motion(or claim motion artifacts; MotionArtifact), can be surveyed as heartbeat signal, the heartbeat signal therefore detected turns S by mistake heartchange to time-domain signal, namely there will be non-heartbeat signal S motion, as shown in Figure 9 C.

Therefore, the accuracy of current Wearable electronic installation detection heart beating is necessary to improve it further.

Summary of the invention

Because the interference that Wearable electronic installation detection physiological signal is vulnerable to motion artifacts reduces its accuracy in detection, the blood processor that main purpose system of the present invention provides a kind of light change in volume to describe waveform (PPG) and method thereof, contribute to the accuracy of measurement improving physiological signal.

The blood processor that the technical way system used for reaching above-mentioned purpose makes this light change in volume describe waveform (PPG) includes:

One smooth change in volume describes waveform (PPG) signal acquisition unit, is capture a PPG signal;

One movable information acquisition unit is capture a movable information;

One first adaptive digital filtering unit, system is connected to this PPG signal acquisition unit and this movable information acquisition unit, to receive a PPG signal and this movable information, and according to this movable information, the motion artifacts interfering signal that one PPG signal comprises is eliminated, and exports the second smooth change in volume and describe waveform (PPG) signal;

One second adaptive digital filtering unit, system is connected to this first adaptive digital filtering unit, to receive the 2nd PPG signal, and capture the cyclical signal in the 2nd PPG signal, and after improving the signal to noise ratio of this cyclical signal, export a physiology characteristic signal.

The processing method that the technical way system used for reaching above-mentioned purpose makes this light change in volume describe waveform (PPG) includes:

A () obtains one first smooth change in volume and describes waveform (PPG) signal and movable information;

B one PPG signal, according to this movable information, is comprised motion artifacts interfering signal and is eliminated by (), and export one the 2nd PPG signal; And

C () captures the cyclical signal in the 2nd PPG signal, and export this physiological feature signal after improving the signal to noise ratio of this cyclical signal.

As shown in the above description, the present invention is mainly during acquisition PPG signal, with an operation information acquisition unit acquisition movable information, because the PPG signal captured from human body can include the interfering signal of motion artifacts, therefore this first adaptive digital filtering unit is according to movable information, by the filtering in this PPG signal of the motion artifacts interfering signal of major part, export the PPG signal after filtering to one second adaptive digital filtering unit again, the signal to noise ratio of cyclical signal is improved with the cyclical signal captured in the 2nd PPG signal, as the foundation reading accurate heartbeat signal.Therefore, PPG signal processing apparatus of the present invention can be applicable to, on Wearable heartbeat detection device, provide correct heartbeat message.

Accompanying drawing explanation

Fig. 1 is the functional block diagram of the first preferred embodiment of PPG signal processing apparatus of the present invention.

Fig. 2 is the use schematic diagram that the present invention is applied to a Wearable electronic installation.

Fig. 3 A to Fig. 3 G is the signal output waveform figure of Fig. 1 part function block.

Fig. 4 A is a part of functional block diagram of the second preferred embodiment of PPG signal processing apparatus of the present invention.

Fig. 4 B is another part functional block diagram of the second preferred embodiment of PPG signal processing apparatus of the present invention.

Fig. 5 A is a functional block diagram of the first adaptive digital filtering unit of Fig. 1 of the present invention.

Fig. 5 B is a functional block diagram of the first adaptive digital filtering unit of Fig. 4 A of the present invention.

Fig. 6 is a functional block diagram of the second adaptive digital filtering unit of Fig. 1 of the present invention.

Fig. 7 A is that Fig. 3 A and Fig. 3 D are converted to a frequency-domain waveform figure.

Fig. 7 B is that Fig. 3 A and Fig. 3 G are converted to a frequency-domain waveform figure.

Fig. 7 C is that Fig. 3 A, Fig. 3 D and Fig. 3 G are converted to a frequency-domain waveform figure.

Fig. 8 is the frequency-domain waveform figure according to changing Fig. 3 G in Fig. 7 C, obtains the time domain beamformer of palmic rate.

Fig. 9 A is the time domain beamformer measuring one of human body physiological feature signal.

Fig. 9 B is that Fig. 9 A is converted to a frequency-domain waveform figure.

Fig. 9 C is the time domain beamformer obtaining palmic rate according to Fig. 9 B.

Fig. 9 D is a time domain ripple figure of normal heartbeat frequency.

Wherein, Reference numeral:

10PPG signal acquisition unit 11 first band filter

12 first normalizer 20 movable information acquisition units

20a tri-axle electrodynamic induction device 21 second direct current level adjustment unit

21a second normalizer 22 delay circuit

30 first adaptive digital filtering unit

301 first FIR digital filter unit

302 subtractor 303 weight adjustment unit

The automatic adjustment unit 32 of 31 first error convergence coefficient removes extremum unit

33 second band filter 34 the 3rd normalizer

40 second adaptive digital filtering unit

401 second FIR digital filter unit

402 subtractor 403 weight adjustment unit

The automatic adjustment unit of 41 second error convergence coefficient

50 wrist 51 blood vessels

60 Wearable electronic installations

Detailed description of the invention

System of the present invention describes waveform (PPG) signal mainly for a kind of light change in volume of reacting physiology signal, and to its blood processor carried out and method thereof, with embodiment, the blood processor of PPG signal of the present invention and the technology contents of method are described in detail below.

First refer to shown in Fig. 1, for the functional block diagram of a preferred embodiment of PPG signal processing apparatus of the present invention, it includes PPG signal acquisition unit 10, movable information acquisition unit 20,1 first adaptive digital filtering unit 30 and an one second adaptive digital filtering unit 40; Wherein this first adaptive digital filtering unit 30 is electrically connected to this PPG signal acquisition unit 10 and this movable information acquisition unit 20, and this second adaptive digital filtering unit 40 is electrically connected to this first adaptive digital filtering unit 30.

Please refer to Fig. 2, if above-mentioned PPG signal acquisition unit 10 is applied to a kind of Wearable electronic installation 60 as accutron, the light-emitting components such as infrared light element, green glow element, HONGGUANG element or laser diode can being adopted, obtaining should one the one PPG signal PPG (as shown in Figure 3A) that shrinks of blood vessel 51 by irradiating the blood vessel 51 of wrist 50.But, in detection wrist, physical activity is non-is simultaneously static, therefore a PPG signal PPG includes the interfering signal of motion artifacts, namely as shown in Figure 7 A, as the PPG signal PPG of Fig. 3 A through fast fourier transform (FFT) to a frequency-domain waveform figure, comprise two obvious crests in frequency-domain waveform, be namely respectively heartbeat signal S heartand the interfering signal S of motion artifacts motion.

As shown in Figure 1, when a PPG signal PPG is applied to specific physiological signal (as the heart beating) of human body, heartbeat signal due to the mankind has specific frequency range, therefore first the PPG signal PPG that this PPG signal acquisition unit 10 exports can be exported to one first band filter 11, noise beyond the frequency range of this first band filter 11 meeting filtering heartbeat signal in a PPG signal PPG, export the PPG signal PPG_F (as shown in Figure 3 B) crossed through bandpass filtering to one first normalizer 12 more again, adjust the PPG signal PPG_C (as shown in Figure 3 C) after the range scale of the PPG signal PPG_F crossed through bandpass filtering, export this first adaptive digital filtering unit 30 again to.In addition, this first normalizer 12 also can use the signal adjustment elements such as one first direct current level adjustment unit, to adjust the direct current level of the PPG signal PPG_F crossed through bandpass filtering.

As Fig. 1 shows, above-mentioned movable information acquisition unit 20 is in order to capture the movable information of this wrist, capture this movable information M in particularly during the above-mentioned PPG signal PPG of acquisition simultaneously, this movable information acquisition unit 20 exports one second direct current level adjustment unit 21 and a delay circuit 22 to, to adjust direct current level and after postponing a period of time, then export this first adaptive digital filtering unit 30 to.In the present embodiment, then coordinate shown in Fig. 4 A, this movable information acquisition unit 20 adopts three axle electrodynamic induction device 20a, therefore this movable information M includes the first axle, the second axle and three-axis moving information X, Y, Z.In addition, the first axle, the second axle and the 3rd axle operation information X, Y, Z that this three axles electrodynamic induction device 20a exports preferably sequentially can export one second direct current level adjustment unit 21 and a delay circuit 22 respectively to, to adjust its range scale direct current level and after postponing a period of time, then export this first adaptive digital filtering unit 30 to.In addition, this the second direct current level adjustment unit 21 and a delay circuit 22 can replace it by one second normalizer 21a, namely as shown in Figure 4 B, the first axle, the second axle and the 3rd axle operation information X, Y, Z that this three axles electrodynamic induction device 20a exports preferably sequentially can export this second normalizer 21a to respectively, to adjust its range scale, then export this first adaptive digital filtering unit 30 to.

The first adaptive digital filtering unit 30 as shown in Figure 1 obtains to be crossed and direct current level or the adjusted PPG signal PPG_C of range scale through bandpass filtering, and the movable information M_C after direct current level or range scale adjustment, this the first adaptive digital filtering unit can according to this movable information M_C, most motion artifacts interfering signal that one PPG signal PPG_C comprises is eliminated, and is exported one the 2nd PPG signal Error_C.At the present embodiment, this first adaptive digital filtering unit 30 can preferably use least square algorithm, and its cost function is: Error_C (k)=PPG_C (k)-m (k); Wherein Error_C (k) error amount that is this; PPG_C (k) is a current PPG signal; M (k) is the numerical value of current movable information after digital filtering.The error amount of this first adaptive digital filtering unit can minimize by this least square algorithm, to eliminate the most motion artifacts interfering signal in a PPG signal.

This first adaptive digital filtering unit 30 realizes the framework of above-mentioned least square algorithm again, and as illustrated in figs. 1 and 5a, it includes one first FIR digital filter unit 301, subtractor 302 and a weight adjustment unit 303.Wherein this first FIR digital filter unit 301 is connected to this movable information acquisition unit 20, this subtractor 302 is exported to after sequentially received this movable information m is given digital filtering, by this subtractor 302, a PPG signal PPG_C and this movable information M_C after digital filtering are subtracted each other, to obtain this error value E rror_C.Because this weight adjustment unit 303 is connected to this first FIR digital filter unit 301 and this subtractor 302, therefore the error value E rror_C after can being subtracted each other according to a PPG signal PPG_C and this movable information M_C after digital filtering, adjust this first FIR digital filter unit 301 in the weighted value W of each digital filtering, until error amount minimizes rear output the 2nd PPG signal Error_C, namely as shown in Figure 3 D.Again please refer to shown in Fig. 7 A, for the 2nd PPG signal Error_C of Fig. 3 D is converted to frequency-domain waveform, known after compared with being converted to frequency-domain waveform equally with Fig. 3 A, the interfering signal S of motion artifacts motion(1) effectively suppressed, the summit in waveform is heartbeat signal S heart(1).

The renewal function of above-mentioned first finite impulse response digital filter 301 is:

M (k+1)=m (k)+w (k) × M_C (k); Wherein m (k+1) is the numerical value of this movable information after digital filtering next time captured, and w (k) is this weighted value of this first FIR digital filter unit, and M_C (k) is current movable information.This weighted value is by following formulae discovery:

W (k+1)=w (k)+2 × μ × Error_C (k) × X_C (k); Wherein w (k+1) is the weighted value next time of this first FIR digital filter unit, and μ is error convergence coefficient.

Refer to shown in Fig. 4 A and Fig. 5 B again, the present embodiment is the movable information acquisition unit with the use of three axle electrodynamic induction device 20a, this three axles electrodynamic induction device 20a is connected with three the first FIR digital filter unit 301, and this error amount of this first adaptive digital filtering unit 30a is:

Error_C (k)=PPG_C (k)-y_0 (k)-y_1 (k)-y_2 (k), and respectively the renewal function of this first FIR digital filter unit 301 is as follows respectively:

y_0(k+1)=y_0(k)+w_0(k)×X_C(k);

y_1(k+1)=y_1(k)+w_1(k)×Y_C(k);

Y_2 (k+1)=y_2 (k)+w_2 (k) × Z_C (k); Wherein:

X_C (k) is current first axle movable information;

Y_0 (k+1) is the numerical value of this first axle movable information after digital filtering next time;

Y_C (k) is current second axle movable information;

Y_1 (k+1) is the numerical value of this second axle movable information after digital filtering next time;

Z_C (k) is current three-axis moving information;

Y_2 (k+1) is the numerical value of this three-axis moving information of next record after digital filtering.

In addition, adjust for making the error convergence coefficient dynamic of this first adaptive digital filtering unit 30 tool, accelerate the speed of error minimize, as shown in Figure 1, can include the automatic adjustment unit 31 of one first error convergence coefficient further, it comprises the step of following this error convergence coefficient of dynamic generation:

A () is by its vector value of PPG signal acquisition;

B one PPG signal is converted to matrix by (), and obtain the inner product of this matrix and transposed matrix, to obtain the intensity of a PPG signal;

C () sets a first-in first-out type relief area;

D () obtains the product of step (b) continuously in this first-in first-out type relief area; And

E () selects to be greater than 0 and the numerical value of the inverse of maximum in relief area is this error convergence coefficient; Wherein this maximum is the max product in this first-in first-out type relief area.

Again as shown in Figure 1, although above-mentioned first adaptive digital filtering unit 20 export one the 2nd PPG signal Error_C can the interfering signal of the most motion artifacts of filtering, but still have other interfering signal, therefore can sequentially export the 2nd PPG signal Error_C to one further and remove extremum unit 32, one second band filter 33 and one the 3rd normalizer 34, extremum unit 32 is gone to set a marginal value by this, exceptional value energy in the 2nd PPG signal Error_C under time domain being greater than marginal value is eliminated, again by the noise (as shown in FIGURE 3 E) beyond the frequency range of this second band filter 33 filtering heartbeat signal, finally by the 3rd normalizer 34, the energy of the 2nd PPG signal is given normalization, can be distributed between 1 and-1, as shown in FIGURE 3 E, to increase the operating dynamic range of system, as illustrated in Figure 3 F.

Then, export the 2nd PPG signal PPG_N after normalization to this second adaptive digital filtering unit 40, after the signal to noise ratio (SNR) of cyclical signal wherein being improved, export a physiology characteristic signal.Shown in Fig. 1 and Fig. 6, this second adaptive digital filtering unit 40 includes one second FIR digital filter unit 401, subtractor 402 and a weight adjustment unit 403 in the present embodiment; Wherein this second FIR digital filter unit 401 receives the 2nd PPG signal PPG_N after normalization through a time delay unit 404, to receive the 2nd PPG signal after delay, and export this subtractor 402 to after the 2nd PPG signal is given digital filtering, this is received the 2nd PPG signal PPG_N after normalization by this subtractor 402, subtract each other with normalization the 2nd PPG signal PPG_N that Δ before this is secondary, to export an error value E rror_P.This error value E rror_P inputs to this weight adjustment unit 403, to adjust the weighted value W_P of the next digital filtering of this second FIR digital filter unit 401 according to this error amount; When Error_P level off to 0 time, output signal after this second FIR digital filter unit 401 digital filtering, be this cyclical signal; Namely as shown in Figure 7 B, the cyclical signal of Fig. 3 G is converted to frequency-domain waveform, and waveform summit is still heartbeat signal S heart(2), then consult shown in Fig. 7 C, find after the frequency-domain waveform two that Fig. 3 D of Fig. 7 A changes is compared, Fig. 3 G cyclical signal S heart(2) peak value in frequency-domain waveform exceeds S compared with Fig. 3 D the 2nd PPG signal Error_C in the peak value of frequency-domain waveform d; Therefore the 2nd PPG signal Error_C after normalization is after the second adaptive digital filtering unit 40 process, really can by cyclical signal S heart(2) signal to noise ratio (SNR) improves.

In the present embodiment, this cyclical signal that this second FIR digital filter unit 401 above-mentioned exports is: PPG_P (k+1)=PPG_P (k)+w_P (k) × PPG_N (k-Δ); Wherein: PPG_P (k+1) is the numerical value of this cyclical signal after digital filtering next time; PPG_P (k) is the numerical value of current cyclical signal after digital filtering; W_P (k) is the weighted value of this this FIR digital filter unit; PPG_N (k-Δ) is this 2nd PPG signal that Δ is secondary before.And the weighted value of this second adaptive digital filtering unit is: w_P (k+1)=w_P (k)+2 × μ × Error_P (k) × PPG_N (k-Δ); Wherein w_P (k+1) is the weighted value of this FIR digital filter unit next time; μ is this error convergence coefficient; And Error_P (k) is this error amount.The error amount of this second adaptive digital filtering unit is: Error_P (k)=PPG_N (k)-PPG_P (k); Wherein PPG_N (k) is current 2nd PPG signal.Therefore when error is close to 0, then this cyclical signal that this second FIR digital filter unit exports is close to actual physiological feature signal (beats), therefore can calculate beats according to cyclical signal.

As shown in Figure 1, for making the error convergence coefficient dynamic of this second adaptive digital filtering unit 30 tool adjust, to accelerate the speed of error minimize, can include the automatic adjustment unit 41 of one second error convergence coefficient equally further, it comprises following steps:

A () is by the 2nd its vector value of PPG signal acquisition;

B 2nd PPG signal is converted to matrix by (), and the inner product to this matrix and a transposed matrix, to obtain the intensity of the 2nd PPG signal;

C () sets a first-in first-out type relief area;

D () obtains the product of step (b) continuously in this first-in first-out type relief area;

E () selects to be greater than 0 and the numerical value reciprocal of maximum is this error convergence coefficient in relief area; Wherein this maximum is the max product in this first-in first-out type relief area.

Beats is calculated for using the cyclical signal of Fig. 3 G, then after passing fast vertical leaf conversion (FFT) namely frequency-domain waveform as shown in Figure 7 B (X-axis represents it with index (index), Y-axis is then through the energy frequency spectrum after FFT), calculate beats F per minute according to frequency-domain waveform heart=index × (fs/n) × 60, wherein fs is sampling frequency, and n is that FFT counts.Please coordinate shown in Fig. 8, be calculated this palmic rate oscillogram under time domain, compare Fig. 9 C and be subject to the heartbeat waveform figure of motion artifacts interference known, palmic rate of the present invention measures the interference being namely no longer subject to motion artifacts.

As shown in the above description, the present invention is mainly during acquisition PPG signal, with an operation information acquisition unit acquisition movable information, because the PPG signal captured from human body can include the interfering signal of motion artifacts, therefore this first adaptive digital filtering unit is according to movable information, by the filtering in this PPG signal of the motion artifacts interfering signal of major part, export the PPG signal after filtering to one second adaptive digital filtering unit again, the signal to noise ratio of cyclical signal is improved with the cyclical signal captured in the 2nd PPG signal, as the foundation reading accurate heartbeat signal.Therefore, PPG signal processing apparatus of the present invention can be applicable to, on Wearable heartbeat detection device, provide correct heartbeat message.

The above is only preferred embodiment of the present invention, not any pro forma restriction is done to the present invention, although the present invention with preferred embodiment openly as above, but and be not used to limit the present invention, any person of ordinary skill in the field, not departing from the scope of technical solution of the present invention, make a little change when above-mentioned disclosed technology contents can be utilized or be revised as the Equivalent embodiments of equivalent variations, in every case be the content not departing from technical solution of the present invention, according to any simple modification that technical spirit of the present invention is done above embodiment, equivalent variations and amendment, all still belong in the scope of technical solution of the present invention.

Claims (21)

1. light change in volume describes a blood processor for waveform, it is characterized in that, comprising:
One smooth change in volume describes waveshape signal acquisition unit, describes waveshape signal in order to capture the first smooth change in volume;
One movable information acquisition unit, in order to capture a movable information;
One first adaptive digital filtering unit, be electrically connected to this light change in volume and describe waveshape signal acquisition unit and this movable information acquisition unit, to receive this first smooth change in volume, waveshape signal and this movable information are described, and according to this movable information, this first smooth change in volume is described the motion artifacts interfering signal that waveshape signal comprises to be eliminated, and export the second smooth change in volume and describe waveshape signal;
One second adaptive digital filtering unit, be electrically connected to this first adaptive digital filtering unit, to receive this second smooth change in volume, waveshape signal is described, and capture this second smooth change in volume and describe cyclical signal in waveshape signal, and after improving the signal to noise ratio of this cyclical signal, export a physiology characteristic signal.
2. light change in volume as claimed in claim 1 describes the blood processor of waveform, and it is characterized in that, this first adaptive digital filtering unit includes:
One first FIR digital filter unit, is electrically connected to this movable information acquisition unit, exports after this movable information of received in sequence is given digital filtering;
One subtractor, connect this light change in volume and describe waveshape signal acquisition unit and this first FIR digital filter unit, this first smooth change in volume is described waveshape signal signal to be subtracted each other with this movable information after digital filtering, describe waveshape signal to export this second smooth change in volume; And
One weight adjustment unit, be connected to this first FIR digital filter unit and this subtractor, according to this first smooth change in volume describe waveshape signal subtracted each other with this movable information after digital filtering after error amount, adjust this first FIR digital filter unit in the weighted value of each digital filtering.
3. light change in volume as claimed in claim 2 describes the blood processor of waveform, it is characterized in that, this first adaptive digital filtering unit uses least square algorithm, and its cost function is:
E [Error_C (k) 2], wherein Error_C (k)=PPG_C (k)-m (k); Wherein:
The error amount that Error_C (k) is this;
PPG_C (k) describes waveshape signal for current first smooth change in volume; And
M (k) is the numerical value of current movable information after digital filtering.
4. light change in volume as claimed in claim 3 describes the blood processor of waveform, and the renewal function of this first finite impulse response digital filter is: m (k+1)=m (k)+w (k) × M_C (k); Wherein:
M (k+1) is the numerical value of this movable information after digital filtering next time captured;
W (k) is this weighted value of this first FIR digital filter unit; And
M_C (k) is current movable information.
5. light change in volume as claimed in claim 4 describes the blood processor of waveform, it is characterized in that, this movable information acquisition unit is three axle electrodynamic induction devices, this movable information includes the first axle, the second axle and three-axis moving information, and this three axles electrodynamic induction device is connected with three the first FIR digital filter unit; Wherein:
This error amount of this first adaptive digital filtering unit is:
Error_C(k)=PPG_C(k)-y_0(k)-y_1(k)-y_2(k);
Respectively the renewal function of this first finite impulse response digital filter is respectively:
y_0(k+1)=y_0(k)+w_0(k)×X_C(k);
y_1(k+1)=y_1(k)+w_1(k)×Y_C(k);
Y_2 (k+1)=y_2 (k)+w_2 (k) × Z_C (k); Wherein:
X_C (k) is current first axle movable information;
Y_0 (k+1) is the numerical value of this first axle movable information after digital filtering next time;
Y_C (k) is current second axle movable information;
Y_1 (k+1) is the numerical value of this second axle movable information after digital filtering next time;
Z_C (k) is current three-axis moving information;
Y_2 (k+1) is the numerical value of this three-axis moving information of next record after digital filtering.
6. light change in volume as claimed in claim 4 describes the blood processor of waveform, and it is characterized in that, this weighted value is by following formulae discovery:
W (k+1)=w (k)+2 × μ × Error_C (k) × X_C (k); Wherein:
W (k+1) is the weighted value next time of this first FIR digital filter unit; And
μ is error convergence coefficient.
7. light change in volume as claimed in claim 6 describes the blood processor of waveform, it is characterized in that, this the first adaptive digital filtering unit includes the automatic adjustment unit of one first error convergence coefficient further, and it comprises the step of following this error convergence coefficient of dynamic generation:
A, describe waveshape signal by this first smooth change in volume and obtain its vector value;
B, this first smooth change in volume is described waveshape signal be converted to matrix, and obtain the inner product of this matrix and transposed matrix, to obtain the intensity that this first smooth change in volume describes waveshape signal;
C, set a first-in first-out type relief area;
D, in this first-in first-out type relief area, obtain continuously the product of step b;
E, selection is greater than 0 and the numerical value reciprocal of maximum is this error convergence coefficient in relief area; Wherein this maximum is the max product in this first-in first-out type relief area.
8. light change in volume as claimed in claim 1 describes the blood processor of waveform, and it is characterized in that, this second adaptive digital filtering unit includes:
One second FIR digital filter unit, this the first adaptive digital filtering unit is connected to through a time delay unit, describe waveshape signal to receive this second smooth change in volume, and this second smooth change in volume is described after waveshape signal gives digital filtering export this physiological feature signal;
One subtractor, connect this first adaptive digital filtering unit and this second FIR digital filter unit, sequentially this first smooth change in volume is described waveshape signal to describe waveshape signal with corresponding the second smooth change in volume through digital filtering and subtracted each other, to export an error amount; And
One weight adjustment unit, be connected to this second FIR digital filter unit and this subtractor, according to this second smooth change in volume describe waveshape signal and received wherein second smooth change in volume before this describe waveshape signal subtracted each other after error amount, adjust the weighted value of each digital filtering of this second FIR digital filter unit.
9. light change in volume as claimed in claim 8 describes the blood processor of waveform, it is characterized in that, this cyclical signal that this second FIR digital filter unit exports is:
PPG_P (k+1)=PPG_P (k)+w_P (k) × PPG_N (k-Δ); Wherein:
PPG_P (k+1) is the numerical value of this cyclical signal after digital filtering next time;
PPG_P (k) is the numerical value of current cyclical signal after digital filtering;
W_P (k) is the weighted value of this this FIR digital filter unit; And
Second smooth change in volume that PPG_N (k-Δ) is Δ before this time describes waveshape signal.
10. light change in volume as claimed in claim 9 describes the blood processor of waveform, and it is characterized in that, the weighted value of this second adaptive digital filtering unit is:
W_P (k+1)=w_P (k)+2 × μ × Error_P (k) × PPG_N (k-Δ); Wherein
W_P (k+1) is the weighted value of this FIR digital filter unit next time;
μ is this error convergence coefficient; And
Error_P (k) is this error amount.
11. light change in volume as claimed in claim 10 describe the blood processor of waveform, and the error amount of this second adaptive digital filtering unit is: Error_P (k)=PPG_N (k)-PPG_P (k); Wherein:
PPG_N (k) describes waveshape signal for current second smooth change in volume.
12. light change in volume as claimed in claim 11 describe the blood processor of waveform, it is characterized in that, this the first adaptive digital filtering unit includes the automatic adjustment unit of one first error convergence coefficient further, and it comprises following steps and produces dynamic error convergence coefficient:
A, describe waveshape signal by this first smooth change in volume and obtain its vector value;
B, this first smooth change in volume is described waveshape signal be converted to matrix, and the inner product to this matrix and a transposed matrix, to obtain the intensity that this first smooth change in volume describes waveshape signal;
C, set a first-in first-out type relief area;
D, in this first-in first-out type relief area, obtain continuously the product of step b;
E, select to be greater than 0 and the numerical value of the inverse of maximum in relief area is this error convergence coefficient; Wherein this maximum is the max product in this first-in first-out type relief area.
13. light change in volume as claimed in claim 11 describe the blood processor of waveform, it is characterized in that, this the second adaptive digital filtering unit includes the automatic adjustment unit of one second error convergence coefficient further, and it comprises following steps and produces dynamic error convergence coefficient:
A, describe waveshape signal by this second smooth change in volume and obtain its vector value;
B, this second smooth change in volume is described waveshape signal be converted to matrix, and the inner product to this matrix and a transposed matrix, to obtain the intensity that this second smooth change in volume describes waveshape signal;
C, set a first-in first-out type relief area
D, in this first-in first-out type relief area, obtain continuously the product of step b;
E, selection is greater than 0 and the numerical value of the inverse of maximum is this error convergence coefficient in relief area; Wherein this maximum is the max product in this first-in first-out type relief area.
14. describe the blood processor of waveform as the light change in volume as described in arbitrary in claim 1 to 13, it is characterized in that, this light change in volume describes waveshape signal acquisition unit and is sequentially connected to this first adaptive digital filtering unit through one first band filter and one first direct current level adjustment unit further, or is sequentially connected to this first adaptive digital filtering unit through this first band filter and one first normalization further.
15. describe the blood processor of waveform as the light change in volume as described in arbitrary in claim 1 to 13, it is characterized in that, this movable information acquisition unit is connected to this first adaptive digital filtering unit through one second direct current level adjustment unit and a delay circuit further, or is connected to this first adaptive digital filtering unit through one second normalizer further.
16. describe the blood processor of waveform as the light change in volume as described in arbitrary in claim 1 to 13, it is characterized in that, this first adaptive digital filtering unit sequentially goes extremum unit, one second band filter and one the 3rd normalizer to be connected to this second adaptive digital filtering unit through one further.
17. 1 kinds of light change in volume describe the processing method of waveform, it is characterized in that, comprising:
A, acquisition one first smooth change in volume describe waveshape signal and movable information;
B, according to this movable information, this first smooth change in volume is described waveshape signal and comprise motion artifacts interfering signal and eliminated, and export one second smooth change in volume and describe waveshape signal;
C, capture this second smooth change in volume and describe a cyclical signal in waveshape signal, and after improving the signal to noise ratio of this cyclical signal, export a physiology characteristic signal.
18. light change in volume as claimed in claim 17 describe the processing method of waveform, it is characterized in that, in this step b, use one first adaptive digital filter to eliminate the motion artifacts interfering signal that this first smooth change in volume describes waveshape signal.
19. light change in volume as claimed in claim 18 describe the processing method of waveform, it is characterized in that, use this first adaptive digital filter to eliminate in this step b motion artifacts interfering signal that this first smooth change in volume describes waveshape signal, comprises following steps:
B1, this movable information is given digital filtering after export;
B2, this first smooth change in volume is described waveshape signal subtracted each other with the corresponding movable information through digital filtering, describe waveshape signal to export this second smooth change in volume; And
B3, to describe according to this first smooth change in volume waveshape signal subtracted each other with the corresponding movable information through digital filtering after error amount, adjust this first adaptive digital filter in the weighted value of each digital filtering.
20. describe the processing method of waveform as the light change in volume as described in arbitrary in claim 17 to 19, it is characterized in that, the noise using this second smooth change in volume of this second adaptive digital filter numeral filtering to describe waveshape signal in this step c to comprise, describe this cyclical signal in waveshape signal to capture this second smooth change in volume, comprise following steps:
C1, this second smooth change in volume described after waveshape signal gives digital filtering export this physiological feature signal;
C2, this first smooth change in volume described waveshape signal and describe waveshape signal with corresponding the second smooth change in volume through digital filtering and subtracted each other, to export an error amount; And
C3, to describe according to this second smooth change in volume waveshape signal and wherein one second smooth change in volume received before this describe waveshape signal subtracted each other after error amount, adjust this second adaptive digital filter in the weighted value of each digital filtering.
21. light change in volume as claimed in claim 20 describe the processing method of waveform, it is characterized in that, respectively an error convergence coefficient of this first and second adaptive digital filter is dynamic error convergence coefficient.
CN201510565554.0A 2014-09-18 2015-09-08 Processing device of light volume change description waveform and method thereof CN105433931A (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US201462051957P true 2014-09-18 2014-09-18
US62/051,957 2014-09-18
TW104125755 2015-08-07
TW104125755A TWI558375B (en) 2014-09-18 2015-08-07 Ppg signal processing device and method thereof

Publications (1)

Publication Number Publication Date
CN105433931A true CN105433931A (en) 2016-03-30

Family

ID=55544898

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510565554.0A CN105433931A (en) 2014-09-18 2015-09-08 Processing device of light volume change description waveform and method thereof

Country Status (1)

Country Link
CN (1) CN105433931A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106691425A (en) * 2016-12-30 2017-05-24 北京工业大学 Wrist heart rate monitoring method of sports bracelet
CN107223036A (en) * 2017-04-10 2017-09-29 深圳市汇顶科技股份有限公司 Object wearing device, the adaptive method and device for filtering out motion artifacts
US10466783B2 (en) 2018-03-15 2019-11-05 Sanmina Corporation System and method for motion detection using a PPG sensor

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1989895A (en) * 2003-05-09 2007-07-04 三星电子株式会社 Ear type apparatus for measuring a bio signal and measuring method therefor
WO2012134395A1 (en) * 2011-03-25 2012-10-04 Nitto Denko Corporation A method of measuring an artefact removed photoplethysmographic (ppg) signal and a measurement system
CN103099615A (en) * 2013-01-23 2013-05-15 深圳市理邦精密仪器股份有限公司 Method and device for eliminating exercise electrocardiosignal interference
CN103596491A (en) * 2011-01-21 2014-02-19 伍斯特理工学院 Physiological parameter monitoring with a mobile communication device
US20140058217A1 (en) * 2012-08-25 2014-02-27 The Board Of Trustees Of The Leland Stanford Junior University Motion artifact mitigation methods and devices for pulse photoplethysmography
CN103781414A (en) * 2011-09-16 2014-05-07 皇家飞利浦有限公司 Device and method for estimating the heart rate during motion
US20140288436A1 (en) * 2012-06-22 2014-09-25 Fitbit, Inc. Wearable heart rate monitor

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1989895A (en) * 2003-05-09 2007-07-04 三星电子株式会社 Ear type apparatus for measuring a bio signal and measuring method therefor
CN103596491A (en) * 2011-01-21 2014-02-19 伍斯特理工学院 Physiological parameter monitoring with a mobile communication device
WO2012134395A1 (en) * 2011-03-25 2012-10-04 Nitto Denko Corporation A method of measuring an artefact removed photoplethysmographic (ppg) signal and a measurement system
CN103781414A (en) * 2011-09-16 2014-05-07 皇家飞利浦有限公司 Device and method for estimating the heart rate during motion
US20140288436A1 (en) * 2012-06-22 2014-09-25 Fitbit, Inc. Wearable heart rate monitor
US20140058217A1 (en) * 2012-08-25 2014-02-27 The Board Of Trustees Of The Leland Stanford Junior University Motion artifact mitigation methods and devices for pulse photoplethysmography
CN103099615A (en) * 2013-01-23 2013-05-15 深圳市理邦精密仪器股份有限公司 Method and device for eliminating exercise electrocardiosignal interference

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106691425A (en) * 2016-12-30 2017-05-24 北京工业大学 Wrist heart rate monitoring method of sports bracelet
CN106691425B (en) * 2016-12-30 2019-06-21 北京工业大学 A kind of wrist rhythm of the heart method of motion bracelet
CN107223036A (en) * 2017-04-10 2017-09-29 深圳市汇顶科技股份有限公司 Object wearing device, the adaptive method and device for filtering out motion artifacts
WO2018187895A1 (en) * 2017-04-10 2018-10-18 深圳市汇顶科技股份有限公司 Wearable device, method and device for self-adaptive filtering of motion interference
CN107223036B (en) * 2017-04-10 2019-03-29 深圳市汇顶科技股份有限公司 Object wearing device, the method and device for adaptively filtering out motion artifacts
US10466783B2 (en) 2018-03-15 2019-11-05 Sanmina Corporation System and method for motion detection using a PPG sensor

Similar Documents

Publication Publication Date Title
Luo et al. A review of electrocardiogram filtering
US5521851A (en) Noise reduction method and apparatus
Ruha et al. A real-time microprocessor QRS detector system with a 1-ms timing accuracy for the measurement of ambulatory HRV
KR20120065399A (en) Method and apparatus for artifacts mitigation with multiple wireless sensors
JP2009131628A (en) System, method and program for vital sign estimation
Peng et al. Extraction of heart rate variability from smartphone photoplethysmograms
WO1997038626A1 (en) Arrhythmia detector
JP2008284352A (en) Coherent signal rejection in ecg
CN103501694B (en) Apparatus and method for measuring physiological signal quality
US8880576B2 (en) Systems and methods for determining respiration information from a photoplethysmograph
EP2757944B1 (en) Systems and methods for determining respiration information from a photoplethysmograph
US9119597B2 (en) Systems and methods for determining respiration information from a photoplethysmograph
KR20130060317A (en) Method and apparatus for low complexity compression of signals
CN103549950B (en) Improved difference threshold detection algorithm for mobile ECG (electrocardiogram) monitoring
WO2006034024A2 (en) Method for adaptive complex wavelet based filtering of eeg signals
Chandrakar et al. Denoising ECG signals using adaptive filter algorithm
Soria-Olivas et al. Application of adaptive signal processing for determining the limits of P and T waves in an ECG
Van Helleputte et al. 18.3 A multi-parameter signal-acquisition SoC for connected personal health applications
ES2559263T3 (en) Handheld device with ability to monitor health
Zou et al. An energy-efficient design for ECG recording and R-peak detection based on wavelet transform
Zhang et al. VLSI friendly ECG QRS complex detector for body sensor networks
Sornmo Time-varying filtering for removal of baseline wander in exercise ECGs
CN1977767B (en) Method for improving respiratory wave identification rate
Tang et al. ECG de-noising based on empirical mode decomposition
Gradolewski et al. Wavelet-based denoising method for real phonocardiography signal recorded by mobile devices in noisy environment

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160330

WD01 Invention patent application deemed withdrawn after publication