CN106388863A - A method and a device for extracting ultrasonic Doppler frequency responses - Google Patents
A method and a device for extracting ultrasonic Doppler frequency responses Download PDFInfo
- Publication number
- CN106388863A CN106388863A CN201610853841.6A CN201610853841A CN106388863A CN 106388863 A CN106388863 A CN 106388863A CN 201610853841 A CN201610853841 A CN 201610853841A CN 106388863 A CN106388863 A CN 106388863A
- Authority
- CN
- China
- Prior art keywords
- frequency
- response
- setting
- signal
- state model
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
Abstract
The invention provides a method and a device for extracting ultrasonic Doppler frequency responses. The method comprises the steps of configuring a state model according to the LCR circuit response principles; determining a frequency response center through setting of an LC in the state model; determining an amplitude bandwidth and response through setting of an R in the state model. The ultrasonic Doppler signal extracting method is simple, can realize analysis and extraction of any frequency components, and can obtain Doppler envelopes through integral operation directly. The method has the advantages of high frequency estimation accuracy, great noise suppressing capability, low algorithm complexity and low calculating resource consumption, and can meet the requirement for instantaneity; the method is free of double frequency multiplication interference of Doppler flow, which is the biggest problem in ultrasonic continuous Doppler envelope extracting and heart rate detection.
Description
Technical field
The present embodiments relate to the technical field of medical science algorithm, more particularly, to a kind of extraction ultrasonic doppler frequency response
Method and device.
Background technology
Time frequency analysis, as the powerful of analysis time-varying non-stationary signal, are focuses of modern signal processing research
One of, describe the time dependent relation of signal frequency.In ultrasonic doppler medical diagnosiss, be by time frequency analysis and
Envelope extraction, physiological parameter needed for realizing extracting from DOPPLER ULTRASOUND SIGNAL, detects blood flow state with no damage, and then is blood
The diagnosis of fluid circulation and angiopathy provides foundation, is widely used on clinical medicine.
For DOPPLER ULTRASOUND SIGNAL, in actual blood flow detection, except intrinsic non-stationary in addition to, low noise
The characteristic such as stronger than, diversified background noise, multiple frequence interference and randomness has a lot of demands to the Time-Frequency Analysis Method of signal.In early days
The frequency domain transform method that uses, digital filtering, the space overlapping mode such as averagely are eliminated one after another.At present, the time frequency analysis of main flow are calculated
Method is with Hilbert-Huang transform method, short time discrete Fourier transform and study hotspot wavelet transformation in recent years.But, this
The complexity of a little main stream approach is high, and required computing resource is big, time-consuming it is difficult to meet supersonic blood and fetal rhythm Doppler is defeated in real time
The demand going out.
Content of the invention
The purpose of the embodiment of the present invention be to propose a kind of method and device extracting ultrasonic doppler frequency response it is intended to
Solve the problems, such as how to realize optional frequency composition is extracted.
For reaching this purpose, the embodiment of the present invention employs the following technical solutions:
In a first aspect, a kind of method extracting ultrasonic doppler frequency response, methods described includes:
According to LCR circuit response theory configuration status model;
By the setting to the LC in described state model, determine frequency response center;
By the setting to the R in described state model, determine bandwidth and the response of amplitude.
Preferably, described according to LCR circuit response theory configuration status model, including:
Ultrasonic echo complex signal after described LCR circuit response theory obtains demodulation
Wherein, described A is the amplitude of signal, described ωdFor the Doppler frequency deviation of signal, describedFor the phase shift of signal, institute
State n and be noise.
Preferably, described frequency response center is determined by the setting to the LC in described state model, including:
According to the frequency in described LCR circuitObtain frequency parameter, and determined according to described frequency parameter
Frequency response center.
Preferably, described the bandwidth of amplitude is determined by the setting to the R in described state model, including:
Determine the truncate length of data according to the points of single computing input data, and width is determined according to described truncate length
The bandwidth of degree.
Preferably, described the response of amplitude is determined by the setting to the R in described state model, including:
If described frequency is f it is determined that described response H is:
Second aspect, a kind of device of extraction ultrasonic doppler frequency response, described device includes:
Configuration module, for according to LCR circuit response theory design point model;
First determining module, for by the setting to the LC in described state model, determining frequency response center;
Second determining module, for by the setting to the R in described state model, determining bandwidth and the response of amplitude.
Preferably, described configuration module, specifically for:
Ultrasonic echo complex signal after described LCR circuit response theory obtains demodulation
Wherein, described A is the amplitude of signal, described ωdFor the Doppler frequency deviation of signal, describedFor the phase shift of signal, institute
State n and be noise.
Preferably, described first determining module, specifically for:
According to the frequency in described LCR circuitObtain frequency parameter, and determined according to described frequency parameter
Frequency response center.
Preferably, described second determining module, specifically for:
Determine the truncate length of data according to the points of single computing input data, and width is determined according to described truncate length
The bandwidth of degree.
Preferably, described second determining module, specifically for:
If described frequency is f it is determined that described response H is:
The embodiment of the present invention provides a kind of method and device extracting ultrasonic doppler frequency response, is rung according to LCR circuit
Answer principle configuration state model;By the setting to the LC in described state model, determine frequency response center;By to described
The setting of the R in state model, determines bandwidth and the response of amplitude.The present invention proposes a kind of simple DOPPLER ULTRASOUND SIGNAL and carries
Take method, it is possible to achieve the analysis to optional frequency composition is extracted, and in digital signal processing, data is discrete, directly tires out
Plus can get Doppler's envelope;Because the setting of LC can accurately extract required for frequency, in theory for Nai Kuisi
Optional frequency under distinct frequence, including the negative frequency of complex signal, because Frequency Estimation accuracy is high;Noise inhibiting ability is strong;
Because the LC concussion that the theory adopting is very common in circuit is theoretical, algorithm complex is low;Computing resource consumption is few, compared to it
He filters the algorithm extracting frequency etc, and LC can once be multiplied accumulating realizes computing, and it is all intended to repeatedly with method
Shift to multiply accumulating and just enable, computing faster, can guarantee that the requirement of real-time;In the setting of LC, do not extract double frequency composition,
Do not disturbed by double frequency, very freedom and flexibility are variable, the method is not disturbed by the double frequency of flow Doppler, and double frequency is done
Disturbing is ultrasonic continuous wave Doppler envelope extraction, the maximum difficult point in heart rate detection.
Brief description
Fig. 1 is a kind of schematic flow sheet of method extracting ultrasonic doppler frequency response provided in an embodiment of the present invention;
Fig. 2 is a kind of structural representation based on LCR circuit model provided in an embodiment of the present invention;
Fig. 3 is that schematic diagram is extracted in a kind of optional frequency response based on LCR model provided in an embodiment of the present invention;
Fig. 4 is that a kind of functional module of device extracting ultrasonic doppler frequency response provided in an embodiment of the present invention is illustrated
Figure.
Specific embodiment
With reference to the accompanying drawings and examples the embodiment of the present invention is described in further detail.It is understood that this
The described specific embodiment in place is used only for explaining the embodiment of the present invention, rather than the restriction to the embodiment of the present invention.In addition also
It should be noted that for the ease of description, illustrate only the part related to the embodiment of the present invention rather than entire infrastructure in accompanying drawing.
With reference to Fig. 1, Fig. 1 is a kind of flow process of method extracting ultrasonic doppler frequency response provided in an embodiment of the present invention
Schematic diagram.
As shown in figure 1, the described method extracting ultrasonic doppler frequency response includes:
Step 101, according to LCR circuit response theory configuration status model;
Wherein, the setting according to LC can be with location response frequency, you can only to obtain the ultrasound echo signal of this frequency, filter
Except remaining interference signal.
Alternatively, described according to LCR circuit response theory configuration status model, including:
Ultrasonic echo complex signal after described LCR circuit response theory obtains demodulation
Wherein, described A is the amplitude of signal, described ωdFor the Doppler frequency deviation of signal, describedFor the phase shift of signal, institute
State n and be noise.
Specifically, the present invention refer to LCR circuit response theory as the design reference to state model.For LCR electricity
For road, by the setting to LC, can determine that frequency response center, the setting to R, can determine that bandwidth and the response of amplitude.Its
It similarly is a radio in fact, old-fashioned radio devices capture radio station by way of coil, to adjust the side of volume
Formula obtains signal bandwidth and amplitude.Mathematical formulae is extracted in the optional frequency response that the present invention is drawn, is based on LCR circuit mould
The thought of type, as shown in Figure 2.
Step 102, by the setting to the LC in described state model, determines frequency response center;
Specifically, this frequency response center is equal to by circuit Meter
Calculation obtains corresponding frequencies parameter.
Being specifically calculated as frequency parameter:The parameter that LC is arranged for model, corresponding with the principle that Fig. 2 is shown,
If for striked frequency f0=10kHz, corresponding Then in Digital Signal Processing
In, as long as it is 1.6 × 10 that corresponding LC meets multiplication income value-5, specific L-value and C value can arbitrarily be arranged in pairs or groups.
Step 103, by the setting to the R in described state model, determines bandwidth and the response of amplitude.
Preferably, described the bandwidth of amplitude is determined by the setting to the R in described state model, including:
Determine the truncate length of data according to the points of single computing input data, and width is determined according to described truncate length
The bandwidth of degree.
Preferably, described the response of amplitude is determined by the setting to the R in described state model, including:
If described frequency is f it is determined that described response H is:
Specifically, in conjunction with LCR model, during actual algorithm is realized, the ultrasonic echo complex signal x after demodulation is such as
Shown in following formula:
Wherein, A is the amplitude of signal, ωdFor the Doppler frequency deviation of signal,For the phase shift of signal, n is noise.For
Frequency f of required frequency response, its system responds H expression formula and is:
From there through simple multiplying, you can obtain signal in this f Frequency pointResponse.For many
For the Doppler signal of frequency content, by arranging multiple frequency response points, integral operation and weighted accumulation can be managed
The Doppler's envelope curve thought, i.e. multifrequency signal extraction is realized in multiple LC weightings, if having demand to large-scale signal,
Expand this frequency response range in the setting of R, because parameter all can adjust one by one, therefore can adjust acquisition preferable
Doppler's envelope, and realize to frequency multiplication, drift, distorted signal suppression.Therefore, it is corresponding LC frequency-selecting to the setting of f, logarithm
Correspond to the wide control of R-band according to the setting of length, as shown in Figure 3.The optional frequency that this algorithm finally achieves based on LCR model rings
Should extract.
Specifically, as shown in figure 3, arriving in the range of Nyquist fs/2 0, fo=fs/4 frequency signal is carried
Take it can be seen that when data length (R value) setting is not both, its frequency response bandwidth is also different, and R value is less, and bandwidth is wider, R
Value is bigger, and bandwidth is narrower.Extraction for remaining frequency is also that, as illustrated, R value is less, bandwidth is wider, R for its responsive bandwidth
Value is bigger, and bandwidth is narrower.
The embodiment of the present invention provides a kind of method extracting ultrasonic doppler frequency response, according to LCR circuit response theory
Configuration status model;By the setting to the LC in described state model, determine frequency response center;By to described state mould
The setting of the R in type, determines bandwidth and the response of amplitude.The present invention proposes a kind of simple DOPPLER ULTRASOUND SIGNAL extraction side
Method, it is possible to achieve the analysis to optional frequency composition is extracted, and in digital signal processing, data is discrete, directly adds up i.e.
Can get Doppler's envelope;Because the setting of LC can accurately extract required for frequency, in theory for Nyquist frequency
Optional frequency under rate, including the negative frequency of complex signal, because Frequency Estimation accuracy is high;Noise inhibiting ability is strong;Because
Using theory be that very common in circuit LC concussion is theoretical, algorithm complex is low;Computing resource consumption is few, compared to other filters
Ripple extracts the algorithm of frequency etc, and LC can once be multiplied accumulating realizes computing, and it is all intended to multi-shift with method
Multiply accumulating and just enable, computing faster, can guarantee that the requirement of real-time;In the setting of LC, do not extract double frequency composition, that is, not
Disturbed by double frequency, very freedom and flexibility are variable, the method is not disturbed by the double frequency of flow Doppler, and the interference of double frequency is
Ultrasonic continuous wave Doppler envelope extraction, the maximum difficult point in heart rate detection.
With reference to Fig. 4, Fig. 4 is a kind of function of device extracting ultrasonic doppler frequency response provided in an embodiment of the present invention
Module diagram.
As shown in figure 4, described device includes:
Configuration module 401, for according to LCR circuit response theory design point model;
First determining module 402, for by the setting to the LC in described state model, determining frequency response center;
Second determining module 403, for by the setting to the R in described state model, determining bandwidth and the sound of amplitude
Should.
Preferably, described configuration module 401, specifically for:
Ultrasonic echo complex signal after described LCR circuit response theory obtains demodulation
Wherein, described A is the amplitude of signal, described ωdFor the Doppler frequency deviation of signal, describedFor the phase shift of signal, institute
State n and be noise.
Preferably, described first determining module 402, specifically for:
According to the frequency in described LCR circuitObtain frequency parameter, and determined according to described frequency parameter
Frequency response center.
Preferably, described second determining module 402, specifically for:
Determine the truncate length of data according to the points of single computing input data, and width is determined according to described truncate length
The bandwidth of degree.
Preferably, described second determining module 403, specifically for:
If described frequency is f it is determined that described response H is:
The embodiment of the present invention provides a kind of device extracting ultrasonic doppler frequency response, according to LCR circuit response theory
Configuration status model;By the setting to the LC in described state model, determine frequency response center;By to described state mould
The setting of the R in type, determines bandwidth and the response of amplitude.The present invention proposes a kind of simple DOPPLER ULTRASOUND SIGNAL extraction side
Method, it is possible to achieve the analysis to optional frequency composition is extracted, and in digital signal processing, data is discrete, directly adds up i.e.
Can get Doppler's envelope;Because the setting of LC can accurately extract required for frequency, in theory for Nyquist frequency
Optional frequency under rate, including the negative frequency of complex signal, because Frequency Estimation accuracy is high;Noise inhibiting ability is strong;Because
Using theory be that very common in circuit LC concussion is theoretical, algorithm complex is low;Computing resource consumption is few, compared to other filters
Ripple extracts the algorithm of frequency etc, and LC can once be multiplied accumulating realizes computing, and it is all intended to multi-shift with method
Multiply accumulating and just enable, computing faster, can guarantee that the requirement of real-time;In the setting of LC, do not extract double frequency composition, that is, not
Disturbed by double frequency, very freedom and flexibility are variable, the method is not disturbed by the double frequency of flow Doppler, and the interference of double frequency is
Ultrasonic continuous wave Doppler envelope extraction, the maximum difficult point in heart rate detection.
Describe the know-why of the embodiment of the present invention above in association with specific embodiment.These descriptions are intended merely to explain this
The principle of inventive embodiments, and the restriction to embodiment of the present invention protection domain can not be construed to by any way.Based on herein
Explanation, those skilled in the art do not need to pay performing creative labour and can associate the other concrete of the embodiment of the present invention
Embodiment, these modes fall within the protection domain of the embodiment of the present invention.
Claims (10)
1. a kind of method extracting ultrasonic doppler frequency response is it is characterised in that methods described includes:
According to LCR circuit response theory configuration status model;
By the setting to the LC in described state model, determine frequency response center;
By the setting to the R in described state model, determine bandwidth and the response of amplitude.
2. method according to claim 1 it is characterised in that described according to LCR circuit response theory configuration status model,
Including:
Ultrasonic echo complex signal after described LCR circuit response theory obtains demodulation
Wherein, described A is the amplitude of signal, described ωdFor the Doppler frequency deviation of signal, describedFor the phase shift of signal, described n
It is noise.
3. method according to claim 1 it is characterised in that described by the setting to the LC in described state model,
Determine frequency response center, including:
According to the frequency in described LCR circuitObtain frequency parameter, and frequency is determined according to described frequency parameter
Responsing center.
4. method according to claim 1 it is characterised in that described by the setting to the R in described state model, really
The bandwidth of tentering degree, including:
Determine the truncate length of data according to the points of single computing input data, and amplitude is determined according to described truncate length
Bandwidth.
5. method according to claim 1 it is characterised in that described by the setting to the R in described state model, really
The response of tentering degree, including:
If described frequency is f it is determined that described response H is:
6. a kind of device extracting ultrasonic doppler frequency response is it is characterised in that described device includes:
Configuration module, for according to LCR circuit response theory design point model;
First determining module, for by the setting to the LC in described state model, determining frequency response center;
Second determining module, for by the setting to the R in described state model, determining bandwidth and the response of amplitude.
7. device according to claim 6 is it is characterised in that described configuration module, specifically for:
Ultrasonic echo complex signal after described LCR circuit response theory obtains demodulation
Wherein, described A is the amplitude of signal, described ωdFor the Doppler frequency deviation of signal, describedFor the phase shift of signal, described n
It is noise.
8. device according to claim 6 is it is characterised in that described first determining module, specifically for:
According to the frequency in described LCR circuitObtain frequency parameter, and frequency is determined according to described frequency parameter
Responsing center.
9. device according to claim 6 is it is characterised in that described second determining module, specifically for:
Determine the truncate length of data according to the points of single computing input data, and amplitude is determined according to described truncate length
Bandwidth.
10. device according to claim 6 is it is characterised in that described second determining module, specifically for:
If described frequency is f it is determined that described response H is:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610853841.6A CN106388863A (en) | 2016-09-27 | 2016-09-27 | A method and a device for extracting ultrasonic Doppler frequency responses |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610853841.6A CN106388863A (en) | 2016-09-27 | 2016-09-27 | A method and a device for extracting ultrasonic Doppler frequency responses |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106388863A true CN106388863A (en) | 2017-02-15 |
Family
ID=57996997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610853841.6A Pending CN106388863A (en) | 2016-09-27 | 2016-09-27 | A method and a device for extracting ultrasonic Doppler frequency responses |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106388863A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110720905A (en) * | 2019-09-20 | 2020-01-24 | 中南大学 | Fetal electrocardiosignal interference suppression system based on array processing |
-
2016
- 2016-09-27 CN CN201610853841.6A patent/CN106388863A/en active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110720905A (en) * | 2019-09-20 | 2020-01-24 | 中南大学 | Fetal electrocardiosignal interference suppression system based on array processing |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101138507B (en) | Doppler bloodstream aural signal processing method and device thereof | |
CN105007130B (en) | The detection method of LFM signals under a kind of low signal-to-noise ratio | |
CN102508031B (en) | Fourier series based measurement method of phase angle of partial discharge pulse | |
CN102355273B (en) | Digital channelizing method and circuit | |
CN102999473A (en) | Detection and parameter estimation method for linear frequency modulation signals | |
CN103944535B (en) | A kind of method of all phase DFT filter group utilizing Frequency Response to configure and device thereof | |
Zeng et al. | Normalized complex Teager energy operator demodulation method and its application to fault diagnosis in a rubbing rotor system | |
CN103519847A (en) | Doppler blood flow velocity estimation method and system based on ultrasonic echo radio frequency signals | |
CN106037743A (en) | Method and device for extracting respiratory rate | |
CN107595247A (en) | A kind of monitoring method and system of the depth of anesthesia based on EEG signals | |
Ma et al. | Analysis and design of modified window shapes for S-transform to improve time–frequency localization | |
CN105337587A (en) | Non maximum extraction system integrated filter bank construction method based on DFT | |
CN103344948A (en) | Method for computing external illuminator radar cross-ambiguity function utilizing sparse Fourier transform | |
CN107359903A (en) | Spread spectrum communication code capture method and device | |
CN106388863A (en) | A method and a device for extracting ultrasonic Doppler frequency responses | |
CN104049246A (en) | Delay inequality estimation method with unknown frequency | |
CN111224672A (en) | Multi-harmonic signal undersampling method based on multi-channel time delay | |
CN107495939A (en) | Live biometric monitoring method, device and system | |
CN104281775B (en) | A kind of method for detecting weak signals that multistable stochastic resonance system is compensated based on parameter | |
CN108919200A (en) | A kind of any radar waveform Photographing On-line method based on FPGA | |
CN106073822A (en) | A kind of superonic spectrum adaptive noise estimation envelope extraction method | |
CN104049247A (en) | Delay inequality estimating method based on frequency variance weight | |
CN103929143A (en) | Method for acquiring amplitude-frequency response characteristic of filter | |
CN114675253A (en) | Water surface micro-amplitude wave frequency estimation method and device, electronic equipment and storage medium | |
CN103580646B (en) | A kind of method for estimating analog filter frequency response characteristic |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination |