CN107766845A - A kind of breathing and BCG method for extracting signal based on light shock sensor - Google Patents
A kind of breathing and BCG method for extracting signal based on light shock sensor Download PDFInfo
- Publication number
- CN107766845A CN107766845A CN201711159859.7A CN201711159859A CN107766845A CN 107766845 A CN107766845 A CN 107766845A CN 201711159859 A CN201711159859 A CN 201711159859A CN 107766845 A CN107766845 A CN 107766845A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- msup
- bcg
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
Abstract
The invention belongs to sensor signal processing and signal extraction technical field, disclose a kind of breathing based on light shock sensor and BCG method for extracting signal and system, by sample, filter, extract frequency reducing, Wavelet Denoising Method, FFT spectrum analysis obtain breathing and BCG numerical value;IIR LPFs extraction respiratory curve, IIR high-pass filterings extraction BCG curves.The present invention using STM32L476RGT6 processors realizes dsp system, suits low cost, clear in structure, reliability and cost-effective, and scalability, the portable the characteristics of modern product such as getting well.Inventive algorithm simplifies circuit structure, and by filtering and high performance pre-filtering is realized in extraction, guarantee meets time domain sampling theorem;The method for first passing through Wavelet Denoising Method FFT again enhances the stationarity of signal, and breathing and BCG are determined using the two big value spectral lines that can pass through amplitude-frequency after FFT.
Description
Technical field
The invention belongs to sensor signal processing and signal extraction technical field, more particularly to a kind of shaken based on light to pass
The breathing of sensor and BCG method for extracting signal.
Background technology
Respiratory rate and BCG are important human body physiological parameters.Concept detection based on microbend fiber and energy loss is exhaled
Inhale and BCG theory it is verified that.In the motion of mechanical disturbance quasi-periodic (breathing chest and body kinematics), anamorphoser plate
(grenadine) extrudes optical fiber and induced a series of along fiber axis.Microbend causes optical coupling to enter radiation mode from kernel boot pattern
Formula, irreversible light loss is caused, and reduce the luminous intensity of radio transceiver.The sensor mat of embedded multimode fibre, one
The embeded processor of individual optoelectronics transceivers and a DSP algorithm, maximum micro-bend sensitiveness by appropriate optical configuration come
Realize.TOSA and ROSA is respectively transmitting and the receiver of optical signal, with this come optical fiber corresponding to matching and hardware circuit.
But because sensor is typically located on bed, or even under mattress, also have personal sleeping position, breathing
Situations such as state, heart condition, snoring, body weight, is all different, signal intensity that sensor detects, signal to noise ratio, signal it is flat
Situations such as stability, is different, belongs to blind source signal detection category.The prior art that the field can be used for be difficult adapt to it is various in the case of
Detection, even ordinary circumstance testing result deviation is also very big.Such as analysis corrugation pitch, FFT is based purely on to realize,
Correlation detection, all can not effectively solve the problem the methods of Wiener filter.
In summary, the problem of prior art is present be:Prior art can not all well adapt to be based on fibre optical sensor
The blind source specificity analysis of sophisticated signal, be particularly due to individual and environmental difference caused by measurement signal difference can not effectively solve
Certainly.Prior art is more using the correlation of signal or stationarity etc. as premise, and actual sensor signal is non-stable, and
Autocorrelation is poor, is not prove effective based on one-dimensional single method processing.At present, industry does not solve the reason for problem and is that
Blind signal analysis is also in period of expansion, and technology is immature, and the degree of accuracy of signal analysis is not high.
The content of the invention
The problem of existing for prior art, the invention provides a kind of breathing based on light shock sensor and BCG
Method for extracting signal.
The present invention is achieved in that a kind of breathing based on light shock sensor and BCG method for extracting signal pass through
Sampling, filtering, extract frequency reducing, Wavelet Denoising Method, FFT spectrum analysis acquisition breathing and BCG numerical value, the extraction breathing of IIR LPFs
Curve, IIR high-pass filterings extraction BCG curves.
Further, described to use high-speed sampling, the filtering uses LPF.Here linear phase fir is used
Lowpass digital filter is extracted to realize.Calculating process is as follows:
Further, the algorithm of the Wavelet Denoising Method is realized to signal denoising using wavelet package transforms multiresolution analysis, is carried
The stationarity of signal is risen, then carries out the extraction of breathing and BCG waveforms and numerical value again.It is specific as follows:
The Mallat algorithms of wavelet package transforms, for l=0,1,2 ..., introduce mark
Moreover,Then, f (t) is in wavelet packet space
On projection can be write as
Mallat algorithmic formulas based on WAVELET PACKET DECOMPOSITION
By primary signal f (t) in wavelet packet space
With
On rectangular projection be designated as respectivelyWithSo,In corresponding specific wavelet packet basis
{μ2l,j,n(t);N ∈ Z } under the coefficient that deploys beIn wavelet packet basis { μ2l+1,j,n(t);N ∈ Z } under open up
The coefficient opened is
Decompose and synthesized again after removing noise coefficient.Wavelet packet synthesis Mallat algorithms be
Synthesize the result of gainedActually primary signal f (t) is in wavelet packet space
On rectangular projectionIn base { μl,j+1,n(t);N ∈ Z } under coefficient.
Further, respiratory waveform and BCG waveforms are extracted respectively after IIR LPFs and IIR high-pass filterings is carried out.
Realized using direct II types structure
Further, the FFT spectrum analysis obtain two high-power spectrum spectral lines, including assert breathing spectral line and BCG spectrums
Line.
Another object of the present invention is to provide a kind of breathing based on light shock sensor and BCG signal extraction systems.
Advantages of the present invention and good effect are:
1. the present invention constructs complete breathing and BCG infomation detection schemes, algorithm structure is compact, on room and time all
It is optimized so that in general embeded processor can complete processing task, it is not necessary to accessory system or high-end processors
Calculated.The present invention using STM32L476RGT6 processors realizes DSP systems, suits that cost is low, clear in structure, reliability
With it is cost-effective, and the characteristics of the modern product such as scalability, portable good, whole electronic system can be integrated into small
In fibre optical sensor mat, external structure and installation are saved, especially enhances stability and anti-destructive, reduction takes after sale
Business at least 30%, reduce production cost at least 15%.
2. inventive algorithm simplifies circuit structure, by filtering and high performance pre-filtering is realized in extraction, ensure to meet
Time domain sampling theorem.The method for first passing through Wavelet Denoising Method FFT again enhances the stationarity of signal, remove the noises such as muscle vibrations
Influence, and snoring and sleep apnea vibration influence etc., it can be determined after using FFT by the two big value spectral line of amplitude-frequency
Breathing and BCG provide premise, and the adaptability of various crowds and measuring condition, including personal sleeping position, breathing is substantially improved
Adaptability of state, heart condition, snoring, body weight etc. etc..
3. after Wavelet Denoising Method, respiratory waveform and BCG waveforms are extracted respectively by IIR digital lowpasses and high-pass filter.
4.FFT algorithms using directly calculate and table look-up be combined by the way of lift calculating speed.
5.FIR filtering sampling linear phase structures, filtering data use circle queue, Lifting Convey speed.
Brief description of the drawings
Fig. 1 is the breathing provided in an embodiment of the present invention based on light shock sensor and BCG method for extracting signal
Fig. 2 is the breathing provided in an embodiment of the present invention based on light shock sensor and BCG signal extraction system diagrams.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to embodiments, to the present invention
It is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not used to
Limit the present invention.
A variety of method integrated uses are avoided jejune blind signal processing method, will differentiated more by the present invention completely
Rate is analyzed, signal denoising, and frequency-selective filtering and FFT power spectrumanalysises etc. carry out organic cooperation, and structure is suitably based on fibre optical sensor
Breathing and BCG detection methods.The application of result of this method, can be by fibre optical sensor, under non-direct contact, in not shadow
In the case of ringing people's normal life and rest, breathing and BCG parameters under accurate measurement sleep quality, it can be widely applied to cure
Institute's public ward patient care, the field such as community endowment monitoring and sleep monitor, is with a wide range of applications.
Below in conjunction with the accompanying drawings 1 and specific embodiment to the present invention application principle be further described.
Breathing and BCG method for extracting signal provided in an embodiment of the present invention based on light shock sensor, including it is following
Step:
Sampling, filtering, extract frequency reducing, Wavelet Denoising Method, FFT and spectrum analysis acquisition breathing and BCG numerical value, IIR low pass filtereds
Ripple extraction respiratory curve, IIR high-pass filterings extraction BCG curves.
As the preferred embodiment of the present invention, described to use high-speed sampling, the filtering uses LPF.
As the preferred embodiment of the present invention, the algorithm of the Wavelet Denoising Method uses wavelet package transforms multiresolution analysis
Realize the extraction for signal denoising, the stationarity of promotion signal, then carrying out breathing and BCG waveforms and numerical value again.
As the preferred embodiment of the present invention, extract and exhale respectively after IIR LPFs and IIR high-pass filterings is carried out
Inhale waveform and BCG waveforms.
Described to use high-speed sampling, the filtering uses LPF.Here extracted using linear phase fir low
Logical digital filter is realized.Calculating process is as follows:
The algorithm of the Wavelet Denoising Method uses wavelet package transforms multiresolution analysis to realize to signal denoising, promotion signal
Stationarity, the extraction of breathing and BCG waveforms and numerical value is then carried out again.It is specific as follows:
The Mallat algorithms of wavelet package transforms, for l=0,1,2 ..., introduce mark
Moreover,Then, f (t) is in wavelet packet space
On projection can be write as
Mallat algorithmic formulas based on WAVELET PACKET DECOMPOSITION
By primary signal f (t) in wavelet packet space
With
On rectangular projection be designated as respectivelyWithSo,In corresponding specific wavelet packet basis
{μ2l,j,n(t);N ∈ Z } under the coefficient that deploys beIn wavelet packet basis { μ2l+1,j,n(t);N ∈ Z } under open up
The coefficient opened is
Decompose and synthesized again after removing noise coefficient.Wavelet packet synthesis Mallat algorithms be
Synthesize the result of gainedActually primary signal f (t) is in wavelet packet space
On rectangular projectionIn base { μl,j+1,n(t);N ∈ Z } under coefficient.
Respiratory waveform and BCG waveforms are extracted respectively after IIR LPFs and IIR high-pass filterings is carried out.Using direct
II types structure is realized
The FFT spectrum analysis obtain two high-power spectrum spectral lines, including assert breathing spectral line and BCG spectral lines.
Fig. 2 is the breathing provided in an embodiment of the present invention based on light shock sensor and BCG signal extraction system diagrams.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention
All any modification, equivalent and improvement made within refreshing and principle etc., should be included in the scope of the protection.
Claims (6)
1. a kind of breathing and BCG method for extracting signal based on light shock sensor, it is characterised in that described to be shaken based on light
The breathing of dynamic sensor and BCG method for extracting signal, by sampling, filtering, extracting frequency reducing, Wavelet Denoising Method, FFT spectrum analysis obtain
Take breathing and BCG numerical value;IIR LPFs extraction respiratory curve, IIR high-pass filterings extraction BCG curves.
Respiratory waveform and BCG waveforms are extracted respectively after IIR LPFs and IIR high-pass filterings is carried out.
2. breathing and BCG method for extracting signal as claimed in claim 1 based on light shock sensor, it is characterised in that
Described to use high-speed sampling, described filter uses LPF, and calculating process is as follows:
<mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>h</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>z</mi>
<mrow>
<mo>-</mo>
<mi>n</mi>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>z</mi>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msup>
<mi>z</mi>
<mi>n</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>h</mi>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&rsqb;</mo>
<msup>
<mi>z</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
</msup>
</mrow>
Wherein, h [n] is the coefficient of wave filter, and N is the exponent number of wave filter, and H (z) is the system function of wave filter.
3. breathing and BCG method for extracting signal as claimed in claim 1 based on light shock sensor, it is characterised in that
The algorithm of Wavelet Denoising Method realized to signal denoising using wavelet package transforms multiresolution analysis, the stationarity of promotion signal, then
The extraction of breathing and BCG waveforms and numerical value is carried out again, it is specific as follows:
The Mallat algorithms of wavelet package transforms, for l=0,1,2 ..., introduce mark
<mrow>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mo>&Integral;</mo>
<mi>R</mi>
</msub>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>&mu;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>l</mi>
<mo>;</mo>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>;</mo>
</mrow>
Moreover,Then, f (t) is in wavelet packet space
<mrow>
<msubsup>
<mi>U</mi>
<mi>j</mi>
<mi>l</mi>
</msubsup>
<mo>=</mo>
<mi>C</mi>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>e</mi>
<mi>s</mi>
<mi>p</mi>
<mi>a</mi>
<mi>n</mi>
<mo>{</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mn>2</mn>
<mfrac>
<mi>j</mi>
<mn>2</mn>
</mfrac>
</msup>
<msub>
<mi>&mu;</mi>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<mi>j</mi>
</msup>
<mi>t</mi>
<mo>-</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>Z</mi>
<mo>}</mo>
<mo>;</mo>
</mrow>
On projection write as
<mrow>
<msubsup>
<mi>g</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>&Element;</mo>
<mi>Z</mi>
</mrow>
</munder>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Mallat algorithmic formulas based on WAVELET PACKET DECOMPOSITION
<mrow>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>&Element;</mo>
<mi>Z</mi>
</mrow>
</munder>
<msub>
<mover>
<mi>h</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>&Element;</mo>
<mi>Z</mi>
</mrow>
</munder>
<msub>
<mover>
<mi>g</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>;</mo>
</mrow>
By primary signal f (t) in wavelet packet space
<mrow>
<msubsup>
<mi>U</mi>
<mi>j</mi>
<mrow>
<mn>2</mn>
<mi>l</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mi>C</mi>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>e</mi>
<mi>s</mi>
<mi>p</mi>
<mi>a</mi>
<mi>n</mi>
<mo>{</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mn>2</mn>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mn>2</mn>
<mfrac>
<mi>j</mi>
<mn>2</mn>
</mfrac>
</msup>
<msub>
<mi>&mu;</mi>
<mrow>
<mn>2</mn>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<mi>j</mi>
</msup>
<mi>t</mi>
<mo>-</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>Z</mi>
<mo>}</mo>
</mrow>
With
<mrow>
<msubsup>
<mi>U</mi>
<mi>j</mi>
<mrow>
<mn>2</mn>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mi>C</mi>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>e</mi>
<mi>s</mi>
<mi>p</mi>
<mi>a</mi>
<mi>n</mi>
<mo>{</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mn>2</mn>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mn>2</mn>
<mfrac>
<mi>j</mi>
<mn>2</mn>
</mfrac>
</msup>
<msub>
<mi>&mu;</mi>
<mrow>
<mn>2</mn>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<mi>j</mi>
</msup>
<mi>t</mi>
<mo>-</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>Z</mi>
<mo>}</mo>
</mrow>
On rectangular projection be designated as respectivelyWith In corresponding specific wavelet packet basis { μ2l,j,n(t);n∈
Z } under the coefficient that deploys be In wavelet packet basis { μ2l+1,j,n(t);N ∈ Z } under the coefficient that deploys be
Decompose and synthesized again after removing noise coefficient, the Mallat algorithms of wavelet packet synthesis are
<mrow>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>&Element;</mo>
<mi>Z</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
<msubsup>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Synthesize the result of gainedFor:Primary signal f (t) is in wavelet packet space
<mrow>
<msubsup>
<mi>U</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</msubsup>
<mo>=</mo>
<mi>C</mi>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>e</mi>
<mi>s</mi>
<mi>p</mi>
<mi>a</mi>
<mi>n</mi>
<mo>{</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mn>2</mn>
<mfrac>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<msub>
<mi>&mu;</mi>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mn>2</mn>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>t</mi>
<mo>-</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>Z</mi>
<mo>}</mo>
</mrow>
On rectangular projectionIn base { μl,j+1,n(t);N ∈ Z } under coefficient.
4. breathing and BCG method for extracting signal as claimed in claim 1 based on light shock sensor, it is characterised in that
Respiratory waveform and BCG waveforms are extracted respectively after IIR LPFs and IIR high-pass filterings is carried out, using direct II types structure
Realize
<mrow>
<mi>y</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
<mi>x</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>-</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mrow>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<mi>y</mi>
</mrow>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>-</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
<mo>.</mo>
</mrow>
Wherein, akAnd bkIt is the coefficient of wave filter, N is the exponent number of wave filter, and M is the input delay quantity of wave filter.
5. breathing and BCG method for extracting signal as claimed in claim 1 based on light shock sensor, it is characterised in that
The FFT spectrum analysis obtain two high-power spectrum spectral lines, including assert breathing spectral line and BCG spectral lines.
6. a kind of breathing and BCG method for extracting signal as claimed in claim 1 based on light shock sensor based on light
The breathing of shock sensor and BCG signal extraction systems.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711159859.7A CN107766845A (en) | 2017-11-20 | 2017-11-20 | A kind of breathing and BCG method for extracting signal based on light shock sensor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711159859.7A CN107766845A (en) | 2017-11-20 | 2017-11-20 | A kind of breathing and BCG method for extracting signal based on light shock sensor |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107766845A true CN107766845A (en) | 2018-03-06 |
Family
ID=61279987
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711159859.7A Pending CN107766845A (en) | 2017-11-20 | 2017-11-20 | A kind of breathing and BCG method for extracting signal based on light shock sensor |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107766845A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108712808A (en) * | 2018-04-04 | 2018-10-26 | 湖南城市学院 | A kind of intelligent environment artistic decoration lamp light control system |
CN108836283A (en) * | 2018-04-17 | 2018-11-20 | 六盘水市人民医院 | It is a kind of based on the Respiratory Medicine of big data patient monitor control method |
CN111160090A (en) * | 2019-11-22 | 2020-05-15 | 新绎健康科技有限公司 | BCG signal noise reduction method and system |
CN112512412A (en) * | 2020-10-29 | 2021-03-16 | 香港应用科技研究院有限公司 | Microbend optical fiber sensor for vital sign monitoring and respiratory and heart rate co-extraction |
CN113080918A (en) * | 2021-03-10 | 2021-07-09 | 杭州澳芯科技有限公司 | BCG-based non-contact heart rate monitoring method and system |
CN113499059A (en) * | 2021-06-01 | 2021-10-15 | 武汉理工大学 | BCG signal processing system and method based on optical fiber sensing non-contact |
CN113509169A (en) * | 2021-08-05 | 2021-10-19 | 成都乐享智家科技有限责任公司 | Multi-parameter-based non-contact sleep apnea detection system and method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005065231A (en) * | 2003-07-31 | 2005-03-10 | Matsushita Electric Ind Co Ltd | Signal processor and signal processing method |
CN103070687A (en) * | 2013-02-06 | 2013-05-01 | 南京理工大学 | Signal processing algorithm of non-contact type vital sign monitoring system |
CN103684348A (en) * | 2013-12-31 | 2014-03-26 | 中国人民解放军国防科学技术大学 | Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter |
US20160051156A1 (en) * | 2013-03-24 | 2016-02-25 | Seoul National University R&Db Foundaton | Film-type biomedical signal measuring apparatus, blood pressure measuring apparatus using the same, cardiopulmonary fitness estimating apparatus, and personal authentication apparatus |
CN105919568A (en) * | 2016-05-24 | 2016-09-07 | 北京千安哲信息技术有限公司 | Gabor transformation based method and device for extracting and analyzing breathing and heartbeat signals |
CN106264502A (en) * | 2016-10-13 | 2017-01-04 | 杭州电子科技大学 | A kind of contactless bio-signal acquisition method |
CN106568988A (en) * | 2016-11-10 | 2017-04-19 | 北京清环智慧水务科技有限公司 | Drainage pipeline water body speed measuring system |
-
2017
- 2017-11-20 CN CN201711159859.7A patent/CN107766845A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005065231A (en) * | 2003-07-31 | 2005-03-10 | Matsushita Electric Ind Co Ltd | Signal processor and signal processing method |
CN103070687A (en) * | 2013-02-06 | 2013-05-01 | 南京理工大学 | Signal processing algorithm of non-contact type vital sign monitoring system |
US20160051156A1 (en) * | 2013-03-24 | 2016-02-25 | Seoul National University R&Db Foundaton | Film-type biomedical signal measuring apparatus, blood pressure measuring apparatus using the same, cardiopulmonary fitness estimating apparatus, and personal authentication apparatus |
CN103684348A (en) * | 2013-12-31 | 2014-03-26 | 中国人民解放军国防科学技术大学 | Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter |
CN105919568A (en) * | 2016-05-24 | 2016-09-07 | 北京千安哲信息技术有限公司 | Gabor transformation based method and device for extracting and analyzing breathing and heartbeat signals |
CN106264502A (en) * | 2016-10-13 | 2017-01-04 | 杭州电子科技大学 | A kind of contactless bio-signal acquisition method |
CN106568988A (en) * | 2016-11-10 | 2017-04-19 | 北京清环智慧水务科技有限公司 | Drainage pipeline water body speed measuring system |
Non-Patent Citations (3)
Title |
---|
SONIA GILABERTEJOAN GOMEZ-CLAPERSRAMON CASANELLARAMON PALLAS-ARENY: "Heart and respiratory rate detection on a bathroom scale based on the ballistocardiogram and the continuous wavelet transform" * |
王艳华: "基于CPLD的心冲击图信号实时处理方法研究" * |
王衍学: "基于小波和支持向量机的风机故障趋势预测研究" * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108712808A (en) * | 2018-04-04 | 2018-10-26 | 湖南城市学院 | A kind of intelligent environment artistic decoration lamp light control system |
CN108836283A (en) * | 2018-04-17 | 2018-11-20 | 六盘水市人民医院 | It is a kind of based on the Respiratory Medicine of big data patient monitor control method |
CN111160090A (en) * | 2019-11-22 | 2020-05-15 | 新绎健康科技有限公司 | BCG signal noise reduction method and system |
CN111160090B (en) * | 2019-11-22 | 2023-09-29 | 新绎健康科技有限公司 | BCG signal noise reduction method and system |
CN112512412A (en) * | 2020-10-29 | 2021-03-16 | 香港应用科技研究院有限公司 | Microbend optical fiber sensor for vital sign monitoring and respiratory and heart rate co-extraction |
CN113080918A (en) * | 2021-03-10 | 2021-07-09 | 杭州澳芯科技有限公司 | BCG-based non-contact heart rate monitoring method and system |
CN113499059A (en) * | 2021-06-01 | 2021-10-15 | 武汉理工大学 | BCG signal processing system and method based on optical fiber sensing non-contact |
CN113509169A (en) * | 2021-08-05 | 2021-10-19 | 成都乐享智家科技有限责任公司 | Multi-parameter-based non-contact sleep apnea detection system and method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107766845A (en) | A kind of breathing and BCG method for extracting signal based on light shock sensor | |
US9901306B2 (en) | Tracking mechanism for heart rate measurements | |
Wang et al. | Contactless respiration monitoring using ultrasound signal with off-the-shelf audio devices | |
Jia et al. | Monitoring a person's heart rate and respiratory rate on a shared bed using geophones | |
US9247346B2 (en) | Apparatus, system and method for noise cancellation and communication for incubators and related devices | |
CN106618497A (en) | Method for monitoring sleep in complicated environment based on channel state information | |
US20190038216A1 (en) | Methods for detecting a sleep disorder and sleep disorder detection devices | |
CN103584847B (en) | Non-contact magnetic induction heart rate and respiration rate synchronous detection method and system | |
CN109998549A (en) | A kind of human body respiration detection method based on WiFi channel state information | |
CN106725488A (en) | A kind of wireless field density breathing detection method, device and breathing detection instrument | |
CN101632588B (en) | Measurement method of blood oxygen saturation for eliminating motion interference | |
CN107530001A (en) | Integrated wearable noise dose meter | |
CN111568425B (en) | Non-contact type multi-person respiration detection method | |
CN111696575A (en) | Low ventilation and apnea detection and identification system based on hybrid neural network model | |
CN106796258A (en) | Slowly varying frequency and the application in health care is tracked in noisy environment | |
Meng et al. | A kind of integrated serial algorithms for noise reduction and characteristics expanding in respiratory sound | |
Chen et al. | Apneadetector: Detecting sleep apnea with smartwatches | |
Li et al. | Mo-sleep: Unobtrusive sleep and movement monitoring via Wi-Fi signal | |
CN113848544A (en) | Human body existence induction detection method and device based on Doppler radar and storage medium | |
CN112244794A (en) | Vital sign detection method and device based on periodic characteristics and storage medium | |
Liu et al. | Human biometric signals monitoring based on wifi channel state information using deep learning | |
Vasireddy et al. | K-band Doppler radar for contact-less overnight sleep marker assessment: a pilot validation study | |
CA2585824A1 (en) | Breathing sound analysis for detection of sleep apnea/hypopnea events | |
Fonseca et al. | Estimating actigraphy from motion artifacts in ECG and respiratory effort signals | |
WO2016176067A1 (en) | Time-domain interference removal and improved tracking mechanism for heart rate measurements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |