WO2017173913A1 - Method and circuit for forming pulse stream during sparse sampling of ultrasonic signal - Google Patents

Method and circuit for forming pulse stream during sparse sampling of ultrasonic signal Download PDF

Info

Publication number
WO2017173913A1
WO2017173913A1 PCT/CN2017/076691 CN2017076691W WO2017173913A1 WO 2017173913 A1 WO2017173913 A1 WO 2017173913A1 CN 2017076691 W CN2017076691 W CN 2017076691W WO 2017173913 A1 WO2017173913 A1 WO 2017173913A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
ultrasonic
frequency
module
signals
Prior art date
Application number
PCT/CN2017/076691
Other languages
French (fr)
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
Application filed by 江苏大学 filed Critical 江苏大学
Publication of WO2017173913A1 publication Critical patent/WO2017173913A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483Details of pulse systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation

Definitions

  • the invention belongs to the technical field of ultrasonic signal sparse sampling, and particularly relates to a method for forming a pulse stream in a sparse sampling of an ultrasonic signal based on a limited innovation rate and a hardware front end physical realization circuit.
  • Ultrasonic testing is an important non-destructive testing method.
  • multi-sensor array structure high detection frequency and long-time detection are often used.
  • Obtain rich ultrasound echo information of the target body The unfavorable consequence is that massive detection data is present, which brings great difficulties to the acquisition, transmission, storage and real-time processing of signals.
  • CS theory has been widely used in signal acquisition in the fields of communication and wireless sensor networks.
  • the hardware implementation models mainly include Analog-to-Information Convertor (AIC) and Modulated Wideband Conversion (MWC).
  • AIC Analog-to-Information Convertor
  • MWC Modulated Wideband Conversion
  • the concept of the new interest rate comes from the FRI (Finite Rate of Innovation) theory, which can be sampled for a specific pulse or pulse stream signal. It was first proposed by Vetterli et al.
  • the FRI theory is a sampling signal scheme combining traditional Shannon sampling theory and sub-Nyquist sampling, and is collected at a new interest rate. Signals with FRI properties can be represented by the known delay and amplitude of the pulse, thus reducing the signal sampling rate and reducing the amount of data.
  • the ultrasonic detection echo signal can be regarded as a superposition of a series of Gaussian pulse signals, so it has the FRI signal property, and the FRI theory can be used for signal acquisition and reconstruction.
  • the FRI sampling method for ultrasonic signals was first published in 2011 by Israel's YC. Eldar et al. (Innovation rate sampling of pulse streams with Application to ultrasound imaging, IEEE Transactions on Signal Process, 2011, 59(4), 1827–1842). Since the original ultrasound signal does not belong to the FRI signal, FRI sampling cannot be performed directly. A pulse stream must be formed from the original detection signal to perform limited Xinxin rate sparse sampling.
  • the existing ultrasonic pulse stream extraction method is mainly realized by a software algorithm, that is, the original signal is first collected at an acquisition rate higher than the Nyquist frequency, and then the acquired digital signal is processed by using a demodulation algorithm to obtain Pulse stream signal.
  • the data Li Zhenkun. Parameterized ultrasonic echo model and its parameter estimation [D]. South China University of Technology. 2013; Lin Weiyi. Ultrasonic low-frequency imaging based on FRI [D]. South China University of Technology.
  • the current software extraction methods for ultrasonic pulse flow mainly include Network detection method, Hilbert transform method and quadrature demodulation method.
  • the software method extracts the signal envelope with high precision, flexibility and versatility.
  • such methods are based on the traditional Shannon sampling theorem, which requires conventional Nyquist sampling of the original signal, which is not a physical implementation of sparse sampling.
  • the present invention is directed to an ultrasonic signal, and a physical forming method and circuit for the ultrasonic pulse stream are proposed.
  • the present invention proposes an ultrasonic pulse flow physical forming method and circuit based on the orthogonal demodulation principle, and uses the analog circuit to realize the direct extraction of the pulse stream from the ultrasonic signal. signal.
  • the method is easy to implement, has high demodulation accuracy and does not require an accurate carrier frequency.
  • the technical solution for implementing the present invention is as follows:
  • a method for pulse formation in ultrasonic signal sparse sampling comprising the following steps:
  • Step 2 the oscillating signal F LO (t) in step 1 is divided by two to form two square wave carrier signals with equal amplitude, same frequency and phase difference of 90 degrees, respectively F 1 (t), F 2 (t);
  • Step 3 multiplying and receiving the ultrasonic echo signals S(t) having the center frequency f 0 and the two square wave signals F 1 (t) and F 2 (t) in step 2, respectively,
  • the ultrasonic echo signal is spectrally shifted to form two orthogonal signals I(t), Q(t) including low frequency components and higher harmonic components;
  • Step 4 performing low-pass filtering on the two orthogonal signals I(t) and Q(t) in step 3, filtering out higher harmonic components, and forming two low-frequency component signals I'(t), Q' (t);
  • Step 5 performing modulo operation on the two orthogonal low-frequency component signals I'(t), Q'(t) in step 4, and acquiring an ultrasonic pulse stream signal A(t) for FRI sparse sampling,
  • the method comprises: squarely adding the two orthogonal low-frequency component signals I′(t) and Q′(t), and performing square root operation on the operation result to obtain an ultrasonic pulse stream signal A(t).
  • the step 2 further comprises: setting a frequency difference between the actual carrier signal frequency and the actual ultrasonic echo signal center frequency to be ⁇ f, and establishing the F 1 (t), F 2 (t) Fourier with the ultrasonic echo signal as a reference
  • the series expansion is as follows:
  • K the carrier amplitude coefficient
  • the step 3 further comprises: considering the ultrasonic echo signal S(t) as a Gaussian pulse stream signal amplitude-modulated by the probe center frequency signal, and establishing the following mathematical expression:
  • L is the number of echoes; Is the amplitude coefficient of the echo; Is the echo bandwidth factor; ⁇ is the arrival time of the echo; f 0 is the center frequency of the ultrasonic echo signal; Is the initial phase; A(t) represents a Gaussian pulse stream signal with a pulse number of L,
  • the parameters of the low pass filter in step 4 are determined according to the following formula:
  • f 0 represents the center frequency of the ultrasonic echo signal
  • f p represents the cut-off frequency of the filter pass band
  • f s represents the cut-off frequency of the filter stop band
  • BW A(t) represents the bandwidth of the ultrasonic pulse flow signal A(t).
  • the present invention also proposes a pulse forming circuit for ultrasonic signal sparse sampling, comprising: a local oscillation module, a quadrature mixing module, a low-pass filtering module, and a modulo operation module;
  • the output of the local oscillation module is connected to the orthogonal mixing module after being divided by two.
  • the orthogonal mixing module, the low-pass filtering module, and the modulo computing module are sequentially connected;
  • the quadrature mixing module is configured to multiply and mix the ultrasonic echo signals with the two carrier signals F 1 (t) and F 2 (t) to obtain two orthogonal signals I(t) and Q. (t);
  • the low pass filtering module is configured to filter out higher harmonic parts of the two orthogonal signals output by the quadrature mixing module to obtain low frequency signals I'(t), Q'(t);
  • the modulo operation module is configured to extract a DC component in the two signals output by the low-pass filter module, and synthesize the pulse stream signal to obtain an ultrasonic pulse stream signal A(t) for FRI sparse sampling.
  • the local oscillation module uses a varactor diode SVC321 and a high-speed CMOS Schmitt inverter 74HC14 to form a voltage-controlled square wave oscillating circuit, and the voltage-dividing circuit provides a control voltage of 0 to 10V, and is added by a DC bias resistor R1. Connect to the varactor diode and isolate the control voltage from the CMOS Schmitt inverter input with a large capacitor C1. The varactor diode capacitance varies from 20pF to 400pF, and the resistor R2 and the varactor diode VC1 are timed. The constant determines the frequency f c of the oscillating signal.
  • the quadrature mixing module is implemented by using an integrated quadrature demodulation chip RF2713, and the signal F LO (t) generated by the local oscillation module is input to the LO input pin of the RF2713 through the capacitor C2 in an AC coupling manner.
  • the RF2713 internally includes a digital frequency divider that divides F LO (t) into two square wave carrier signals with a phase difference of 90 degrees; the ultrasonic echo signal S(t) is AC coupled through capacitor C3.
  • the method is input into RF2713, and is mixed by two internal Gilbert mixing units and two square wave carrier signals respectively, and the two output signals I(t) and Q(t) after mixing are respectively passed through capacitors C7 and C8. Output in AC coupling.
  • the low-pass filter module comprises a third-order linear phase low-pass filter of a SallenKey structure by an integrated operational amplifier chip AD847 and a plurality of resistor capacitors.
  • the modulo operation module includes a two-way square operation circuit, an add-on operation circuit, and an open-square root operation circuit, and the two-way square operation circuits are all connected to the addition operation circuit, and the addition operation circuit and The square root computing circuit is connected;
  • the square operation circuit and the square root operation circuit are both realized by an integrated analog multiplier chip AD734 and a plurality of resistor capacitors, and the adder circuit is composed of an integrated op amp chip AD847 to form an in-phase adder.
  • the present invention is different from the existing ultrasonic signal FRI sampling method for extracting ultrasonic pulse flow by software method, and the ultrasonic pulse flow signal is directly extracted by the hardware circuit, and is suitable for the ultrasonic signal FRI sampling system.
  • the proposed ultrasonic pulse flow physical formation method based on the orthogonal demodulation principle can accurately extract the pulse flow signal from the original ultrasonic echo signal.
  • the method allows the carrier frequency to have a certain deviation from the center frequency of the actual echo signal, which reduces the design complexity.
  • the designed ultrasonic pulse stream physically forms a circuit, and the quadrature demodulation process is realized by the analog circuit built by the integrated chipset, and the circuit structure is simple and easy to implement.
  • the hardware implementation of the pulse stream is the key to the physical realization of the ultrasonic signal sparse sampling. While ensuring the information integrity in the signal, the data collection amount can be greatly reduced, and the problem that the conventional sampling method has too large data amount is solved.
  • FIG. 1 is a schematic diagram of a method for physically forming an ultrasonic pulse stream in the present invention
  • FIG. 5 is a schematic circuit diagram of a local oscillation module according to an embodiment of the present invention.
  • FIG. 6 is a schematic circuit diagram of a quadrature mixing module according to an embodiment of the present invention.
  • FIG. 7 is a schematic circuit diagram of a low pass filter module according to an embodiment of the present invention.
  • FIG. 8 is a schematic circuit diagram of a modulo operation module according to an embodiment of the present invention.
  • S1-local oscillation module S2-orthogonal mixing module, S3-low-pass filter module, S4-modulo operation module.
  • the actual ultrasonic echo signal used in this embodiment is as shown in FIG.
  • the center frequency of the ultrasonic probe is 5 MHz, and the actual measured echo center frequency f 0 is about 3.5 MHz.
  • the frequency f c of the oscillation output signal F LO (t) of the local oscillator is determined according to the center frequency of the actual echo signal:
  • the oscillation signal F LO (t) Since the actual ultrasonic echo signal center frequency f 0 is affected by external factors, a certain frequency fluctuation will occur, and the oscillation signal F LO (t) has a certain limit between the theoretical calculation value and the theoretical calculation value due to the performance limitation of the analog circuit. Frequency deviation. Therefore, there is a certain frequency difference ⁇ f between the actual carrier signal frequency generated by the frequency division and the actual echo signal center frequency. Taking the echo signal as a reference, the carrier signal frequency is (f 0 + ⁇ f), and the Fourier series expansion of the two carrier signals F 1 (t) and F 2 (t) are:
  • K the carrier amplitude coefficient
  • the ultrasonic echo signal S(t) can be regarded as a Gaussian pulse stream signal amplitude-modulated by the probe center frequency signal, which can be approximated by the following mathematical model:
  • L is the number of echoes; Is the amplitude coefficient of the echo; Is the echo bandwidth factor; ⁇ is the arrival time of the echo; f 0 is the center frequency of the echo signal; Is the initial phase; A(t) represents a Gaussian pulse stream signal with a pulse number of L,
  • the two signals I(t) and Q(t) formed after the mixing respectively contain low frequency components. And higher harmonic components with a frequency of 2f 0 + ⁇ f or more.
  • Low-pass filtering filters out higher-order harmonic components with frequencies above 2f 0 + ⁇ f in I(t) and Q(t), leaving only low-frequency components, resulting in low-frequency quadrature signals I'(t), Q' (t).
  • the setting of the filter parameters can be determined according to the following formula:
  • f 0 represents the center frequency of the ultrasonic echo signal
  • f p represents the cut-off frequency of the filter pass band
  • f s represents the cut-off frequency of the filter stop band
  • BW A(t) represents the bandwidth of the ultrasonic pulse flow signal A(t).
  • the ultrasonic pulse stream signal A(t) is thus obtained.
  • the method allows a certain frequency difference between the carrier signal frequency and the echo signal center frequency, which relaxes the limitation of the design of the hardware circuit.
  • FIG. 1 it is a schematic diagram of a physical formation method of ultrasonic pulse flow in the present invention. Including: local oscillation a module S1, a quadrature mixing module S2, a low-pass filtering module S3, and a modulo operation module S4; the output of the local oscillation module is connected to the orthogonal mixing module after being divided by two, the orthogonal mixing The frequency module, the low-pass filter module, and the modulo operation module are sequentially connected;
  • the quadrature mixing module is configured to multiply and mix the ultrasonic echo signals with the two carrier signals F 1 (t) and F 2 (t) to obtain two orthogonal signals I(t) and Q. (t);
  • the low pass filtering module is configured to filter out higher harmonic parts of the two orthogonal signals output by the quadrature mixing module to obtain low frequency signals I'(t), Q'(t);
  • the modulo operation module is configured to extract a DC component in the two signals output by the low-pass filter module, and synthesize the pulse stream signal to obtain an ultrasonic pulse stream signal A(t) for FRI sparse sampling.
  • the local oscillation module S1 in the present invention is as shown in FIG.
  • a square wave oscillating signal having a frequency twice the center frequency of the ultrasonic echo signal is generated.
  • the module uses a varactor diode SVC321 and a high speed CMOS Schmitt inverter 74HC14 to form a voltage controlled square wave oscillating circuit.
  • a voltage of 0 to 10 V is supplied from the voltage dividing circuit, and is applied to the varactor through a DC bias resistor R1, and the control voltage is isolated from the input end of the CMOS chip by the large capacitor C1.
  • the capacitance of the varactor diode ranges from about 20pF to 400pF.
  • the time constant of the resistor R2 and the capacitance VC1 of the varactor determines the frequency of the oscillating signal.
  • the oscillating frequency f c can be approximated by the following equation:
  • the orthogonal mixing module S2 in the present invention is as shown in FIG.
  • the local oscillation signal F LO (t) is divided by two to form two square wave carrier signals with a phase difference of 90 degrees, and the original ultrasonic signal is mixed with the two carrier signals to output two orthogonal signals I(t). Q(t).
  • the module uses the integrated quadrature demodulation chip RF2713 to achieve its function.
  • the local oscillating signal F LO (t) is ac-coupled to the LO input pin of the RF2713 via capacitor C2. It internally includes a digital divider that divides F LO (t) by two to a phase difference of 90. Square wave carrier signal.
  • the ultrasonic echo signal S(t) is input to the RF2713 by means of a capacitor C3 in an AC-coupled manner, and is separately mixed with the two carrier signals by an internal two-way Gilbert mixing unit.
  • the two output signals I(t) and Q(t) are outputted by AC coupling through capacitors C7 and C8, respectively.
  • the low pass filter module S3 in the present invention is as shown in FIG. Used to filter out the input of the quadrature mixing module S2
  • the two orthogonal signals I(t), Q(t) are the higher harmonic components.
  • the module includes two filtering circuits with the same parameters, as shown in (a) and (b), respectively filtering the I(t), Q(t).
  • Each channel uses an integrated operational amplifier chip AD847 to form a third-order linear phase low-pass filter of the SallenKey structure.
  • the modulo operation module S4 in the present invention is as shown in FIG. It is used to extract the DC component of the two signals outputted by the S3 low-pass filter module, and synthesize the pulse stream signal.
  • the module includes two square operation circuits, one addition circuit and an square root operation circuit.
  • the square operation and the square root operation circuit are composed of an integrated analog multiplier chip AD734, and the addition circuit uses an integrated operational amplifier chip AD847 to form an in-phase adder.
  • the ultrasonic pulse flow signal extracted from the ultrasonic echo signal shown in FIG. 2 by the ultrasonic pulse flow physical formation method and circuit proposed by the present invention is as shown in FIG. 4.
  • 3 is an ultrasonic pulse stream signal extracted from the ultrasonic echo signal shown in FIG. 2 using a conventional quadrature demodulation software method.
  • the orthogonal demodulation algorithm has the highest precision, which can meet the waveform requirements of the pulse stream in the sparse sampling of ultrasonic signals. It can be seen from the comparison between FIG. 3 and FIG. 4 that the pulse stream extracted by the hardware circuit of the present invention is very close to the pulse stream waveform extracted by the orthogonal demodulation algorithm, and can satisfy the waveform of the pulse stream in the hardware sparse sampling. Claim.
  • the existing algorithm needs to first obtain the signal according to the conventional sampling, and then form the pulse flow through the algorithm, thereby realizing the direct hardware formation of the pulse flow, and has real-time performance.
  • the invention is the key to physically achieving sparse sampling of ultrasonic signals.

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

Disclosed are a method and a circuit for forming a pulse stream during sparse sampling of an ultrasonic signal. The method comprises: generating, by a local oscillator, an oscillating signal having a frequency double a center frequency of an ultrasonic signal S(t); dividing the oscillating signal by two and performing modulation and mixing on the respective signals and the S(t); passing two modulated and mixed signals through a low-pass filter to generate two output signals; and squaring the two output signals respectively, adding the squares, and calculating a square root of the sum, so as to obtain a final output signal A(t), wherein A(t) is a pulse stream obtained by detecting the ultrasonic signal S(t). The circuit comprises: a local oscillating module (S1), an orthogonal mixing module (S2), a low-pass filtering module (S3), and a modulo operation module (S4). The method and the circuit of the present invention are applicable to a sparse sampling system for an ultrasonic signal based on finite rate of innovation, so as to achieve sparse sampling at a rate much lower than the conventional Nyquist sampling rate, and thus resolving the issue in conventional sampling methods of overly large amounts of data, while retaining original signals and information, thereby enabling real-time formation of an ultrasonic pulse stream.

Description

一种超声信号稀疏采样中脉冲流形成的方法与电路Method and circuit for forming pulse flow in ultrasonic signal sparse sampling 技术领域Technical field
本发明属于超声信号稀疏采样技术领域,特别涉及一种基于有限新息率的超声信号稀疏采样中脉冲流的形成方法及硬件前端物理实现电路。The invention belongs to the technical field of ultrasonic signal sparse sampling, and particularly relates to a method for forming a pulse stream in a sparse sampling of an ultrasonic signal based on a limited innovation rate and a hardware front end physical realization circuit.
背景技术Background technique
超声波检测是一种重要无损检测手段,为了提高检测效率、提高其对材料内部缺陷的检出率和成像分辨率,往往采用多传感器的阵列结构、高的检测频率、以及长时间的检测,以获得目标体的丰富的超声回波信息。但是,其带来的不利后果是出现海量检测数据,给信号的采集、传输、存储以及实时处理带来极大的困难。为了解决这一问题,科研人员探索采用常规香农-奈奎斯特采样方法以外的新型数据采样方法,在不丢失或尽量少丢失目标体信息的前提下,减少采集数据量。Ultrasonic testing is an important non-destructive testing method. In order to improve the detection efficiency and improve the detection rate and imaging resolution of internal defects of materials, multi-sensor array structure, high detection frequency and long-time detection are often used. Obtain rich ultrasound echo information of the target body. However, the unfavorable consequence is that massive detection data is present, which brings great difficulties to the acquisition, transmission, storage and real-time processing of signals. In order to solve this problem, researchers explored a new data sampling method other than the conventional Shannon-Nyquist sampling method to reduce the amount of data collected without losing or losing the target information as little as possible.
近几年来,围绕降低超声信号采样速率,减少采集数据量的新方法不断涌现。其代表性的采样方法主要有压缩传感技术和有限新息率采样技术。这些方法均属于信号的稀疏采样方法。尽管这些方法均可有效地降低对信号的采样速率,但目前还主要停留在理论分析阶段,在理论上实现了对信号的稀疏采样,并验证了经这些采样方法对信号进行采样后,信号中的有用信息可以重建。压缩传感(CompressiveSensing,CS)理论的出现打破了常规奈奎斯特采样对最低采样频率的限制,为实现高速信号的低速采集奠定了基础。目前,CS理论已广泛应用于通信、无线传感网络等领域的信号采集中。其硬件实现模型主要有模拟/信息转换器(Analog-to-Information Convertor,AIC)及调制宽带转换器(Modulated Wideband Conversion,MWC)等。新息率的概念来源于FRI(Finite Rate of Innovation)理论,可针对特定的脉冲或脉冲流信号进行采样,它最早由Vetterli等人提出来。FRI理论是传统的香农采样理论和次奈奎斯特采样相结合的采样信号方案,以新息率速率进行采集。具有FRI性质的信号可以由已知脉冲的时延和幅值完备的表示,因此可以降低信号采样速率,减少数据量。超声波检测回波信号可以看作是由一系列高斯脉冲信号叠加组成,因此具有FRI信号性质,可以采用FRI理论进行信号采集和重建。超声信号的FRI采样方法最早由以色列的YC.Eldar等人在2011年发表的论文(Innovation rate sampling of pulse streams with  application to ultrasound imaging,IEEE Transactions on Signal Process,2011,59(4),1827–1842)当中提出。由于原始超声信号不属于FRI信号,不能直接进行FRI采样。必须由原始检测信号形成脉冲流才能进行有限新息率稀疏采样。In recent years, new methods have been emerging around reducing the sampling rate of ultrasonic signals and reducing the amount of data collected. Its representative sampling methods mainly include compression sensing technology and limited innovation rate sampling technology. These methods are all sparse sampling methods of signals. Although these methods can effectively reduce the sampling rate of the signal, at present, it mainly stays in the theoretical analysis stage, theoretically realizes the sparse sampling of the signal, and verifies that after the signal is sampled by these sampling methods, the signal is Useful information can be rebuilt. The emergence of Compressive Sensing (CS) theory breaks the limitation of conventional Nyquist sampling on the lowest sampling frequency, which lays a foundation for low-speed acquisition of high-speed signals. At present, CS theory has been widely used in signal acquisition in the fields of communication and wireless sensor networks. The hardware implementation models mainly include Analog-to-Information Convertor (AIC) and Modulated Wideband Conversion (MWC). The concept of the new interest rate comes from the FRI (Finite Rate of Innovation) theory, which can be sampled for a specific pulse or pulse stream signal. It was first proposed by Vetterli et al. The FRI theory is a sampling signal scheme combining traditional Shannon sampling theory and sub-Nyquist sampling, and is collected at a new interest rate. Signals with FRI properties can be represented by the known delay and amplitude of the pulse, thus reducing the signal sampling rate and reducing the amount of data. The ultrasonic detection echo signal can be regarded as a superposition of a series of Gaussian pulse signals, so it has the FRI signal property, and the FRI theory can be used for signal acquisition and reconstruction. The FRI sampling method for ultrasonic signals was first published in 2011 by Israel's YC. Eldar et al. (Innovation rate sampling of pulse streams with Application to ultrasound imaging, IEEE Transactions on Signal Process, 2011, 59(4), 1827–1842). Since the original ultrasound signal does not belong to the FRI signal, FRI sampling cannot be performed directly. A pulse stream must be formed from the original detection signal to perform limited Xinxin rate sparse sampling.
现有的超声脉冲流的提取方法主要是通过软件算法来实现,即先将原始信号以高于奈奎斯特频率的采集速率进行采集,然后利用解调算法对采集的数字信号进行处理,得到脉冲流信号。根据资料(卢振坤.参数化的超声回波模型及其参数估计[D].华南理工大学.2013;林伟毅.基于FRI的超声波低频成像[D].华南理工大学.2013;曹文,刘春梅,胡莉.一种超声回波信号的数字正交检波方法及FPGA实现[J].西南科技大学学报.2006,21.第3期,56-60),目前超声脉冲流的软件提取方法主要有包络检波法,希尔伯特变换法和正交解调法等。软件方法提取信号包络具有精度高,灵活性、通用性强等特点。但是此类方法是建立在传统的香农采样定理的基础上,需要对原始信号进行常规奈奎斯特采样,不属于稀疏采样的物理实现方法。The existing ultrasonic pulse stream extraction method is mainly realized by a software algorithm, that is, the original signal is first collected at an acquisition rate higher than the Nyquist frequency, and then the acquired digital signal is processed by using a demodulation algorithm to obtain Pulse stream signal. According to the data (Lu Zhenkun. Parameterized ultrasonic echo model and its parameter estimation [D]. South China University of Technology. 2013; Lin Weiyi. Ultrasonic low-frequency imaging based on FRI [D]. South China University of Technology. 2013; Cao Wen, Liu Chunmei, Hu Li.A digital quadrature detection method for ultrasonic echo signals and FPGA implementation[J].Journal of Southwest University of Science and Technology.2006,21.3, 56-60), the current software extraction methods for ultrasonic pulse flow mainly include Network detection method, Hilbert transform method and quadrature demodulation method. The software method extracts the signal envelope with high precision, flexibility and versatility. However, such methods are based on the traditional Shannon sampling theorem, which requires conventional Nyquist sampling of the original signal, which is not a physical implementation of sparse sampling.
要想使稀疏采样方法真正用于实际,必须从根本上解决脉冲流的物理实现问题,本发明专门针对超声信号,提出了一种超声脉冲流的物理形成方法和电路。In order to make the sparse sampling method be used in practice, the physical realization of the pulse stream must be fundamentally solved. The present invention is directed to an ultrasonic signal, and a physical forming method and circuit for the ultrasonic pulse stream are proposed.
发明内容Summary of the invention
为了解决现有的软件算法不能实现超声信号FRI直接物理采样的问题,本发明提出一种基于正交解调原理的超声脉冲流物理形成方法和电路,利用模拟电路来实现超声信号直接提取脉冲流信号。该方法具有易于实现,解调精度高且无需精确的载波频率的特点。实现本发明的技术方案如下:In order to solve the problem that the existing software algorithm cannot realize the direct physical sampling of the ultrasonic signal FRI, the present invention proposes an ultrasonic pulse flow physical forming method and circuit based on the orthogonal demodulation principle, and uses the analog circuit to realize the direct extraction of the pulse stream from the ultrasonic signal. signal. The method is easy to implement, has high demodulation accuracy and does not require an accurate carrier frequency. The technical solution for implementing the present invention is as follows:
一种超声信号稀疏采样中脉冲形成的方法,包括如下步骤:A method for pulse formation in ultrasonic signal sparse sampling, comprising the following steps:
步骤1,由本地振荡器产生固定频率的方波振荡信号FLO(t),其振荡频率fc满足:fc=2f0;其中,f0为超声回波信号中心频率;Step 1: The local oscillator generates a square wave oscillating signal F LO (t) of a fixed frequency, and the oscillating frequency f c satisfies: f c = 2f 0 ; wherein f 0 is the center frequency of the ultrasonic echo signal;
步骤2,将步骤1中所述振荡信号FLO(t)进行二分频,形成两路幅值相等、频率相同、相位差为90度的方波载波信号,分别为F1(t)、F2(t); Step 2, the oscillating signal F LO (t) in step 1 is divided by two to form two square wave carrier signals with equal amplitude, same frequency and phase difference of 90 degrees, respectively F 1 (t), F 2 (t);
步骤3,将接收到的中心频率为f0的超声回波信号S(t)分别与步骤2中所述两路方波信号F1(t)、F2(t)相乘混频,对超声回波信号进行频谱搬移,形成两路包含低频分量和高次谐波分量的正交信号I(t)、Q(t); Step 3, multiplying and receiving the ultrasonic echo signals S(t) having the center frequency f 0 and the two square wave signals F 1 (t) and F 2 (t) in step 2, respectively, The ultrasonic echo signal is spectrally shifted to form two orthogonal signals I(t), Q(t) including low frequency components and higher harmonic components;
步骤4,对步骤3中所述两路正交信号I(t)、Q(t)进行低通滤波,滤除高次谐波分量,形成两路低频分量信号I′(t)、Q′(t); Step 4, performing low-pass filtering on the two orthogonal signals I(t) and Q(t) in step 3, filtering out higher harmonic components, and forming two low-frequency component signals I'(t), Q' (t);
步骤5,对步骤4中所述两路正交的低频分量信号I′(t),Q′(t)进行取模运算,获取用于FRI稀疏采样的超声脉冲流信号A(t),具体包括:对所述两路正交的低频分量信号I′(t)、Q′(t)分别平方相加,再对运算结果进行开平方根运算,得到超声脉冲流信号A(t)。 Step 5, performing modulo operation on the two orthogonal low-frequency component signals I'(t), Q'(t) in step 4, and acquiring an ultrasonic pulse stream signal A(t) for FRI sparse sampling, The method comprises: squarely adding the two orthogonal low-frequency component signals I′(t) and Q′(t), and performing square root operation on the operation result to obtain an ultrasonic pulse stream signal A(t).
进一步优选方案,步骤2还包括:设实际载波信号频率与实际超声回波信号中心频率存在频差为Δf,以超声回波信号为参考,建立F1(t)、F2(t)的傅立叶级数展开式分别如下:Further preferably, the step 2 further comprises: setting a frequency difference between the actual carrier signal frequency and the actual ultrasonic echo signal center frequency to be Δf, and establishing the F 1 (t), F 2 (t) Fourier with the ultrasonic echo signal as a reference The series expansion is as follows:
Figure PCTCN2017076691-appb-000001
Figure PCTCN2017076691-appb-000001
Figure PCTCN2017076691-appb-000002
Figure PCTCN2017076691-appb-000002
式中,K为载波幅值系数,n=1,3,5,…。Where K is the carrier amplitude coefficient, n = 1, 3, 5, ....
进一步优选方案,步骤3还包括:将超声回波信号S(t)看成是由探头中心频率信号进行幅值调制的高斯脉冲流信号,并建立如下数学表达式:Further preferably, the step 3 further comprises: considering the ultrasonic echo signal S(t) as a Gaussian pulse stream signal amplitude-modulated by the probe center frequency signal, and establishing the following mathematical expression:
Figure PCTCN2017076691-appb-000003
Figure PCTCN2017076691-appb-000003
式中,L为回波个数;
Figure PCTCN2017076691-appb-000004
为回波幅值系数;
Figure PCTCN2017076691-appb-000005
为回波带宽因子;τ为回波到达时刻;f0为超声回波信号中心频率;
Figure PCTCN2017076691-appb-000006
为初相位;A(t)表示脉冲数为L的高斯脉冲流信号,
Figure PCTCN2017076691-appb-000007
Where L is the number of echoes;
Figure PCTCN2017076691-appb-000004
Is the amplitude coefficient of the echo;
Figure PCTCN2017076691-appb-000005
Is the echo bandwidth factor; τ is the arrival time of the echo; f 0 is the center frequency of the ultrasonic echo signal;
Figure PCTCN2017076691-appb-000006
Is the initial phase; A(t) represents a Gaussian pulse stream signal with a pulse number of L,
Figure PCTCN2017076691-appb-000007
进一步优选方案,步骤4中所述低通滤波器的参数根据下式确定:Further preferably, the parameters of the low pass filter in step 4 are determined according to the following formula:
Figure PCTCN2017076691-appb-000008
Figure PCTCN2017076691-appb-000008
其中,f0表示超声回波信号中心频率;fp表示滤波器通带截止频率;fs表示滤波器阻带截止频率;BWA(t)表示超声脉冲流信号A(t)的带宽。Where f 0 represents the center frequency of the ultrasonic echo signal; f p represents the cut-off frequency of the filter pass band; f s represents the cut-off frequency of the filter stop band; BW A(t) represents the bandwidth of the ultrasonic pulse flow signal A(t).
为了将上述方法在硬件底层上实现,本发明还提出了一种超声信号稀疏采样中脉冲形成的电路,包括:本地振荡模块、正交混频模块、低通滤波模块、取模运算模块;所述本地振荡模块的输出经二分频后与所述正交混频模块相连接,所 述正交混频模块、所述低通滤波模块、所述取模运算模块依次相连接;In order to implement the above method on the hardware underlayer, the present invention also proposes a pulse forming circuit for ultrasonic signal sparse sampling, comprising: a local oscillation module, a quadrature mixing module, a low-pass filtering module, and a modulo operation module; The output of the local oscillation module is connected to the orthogonal mixing module after being divided by two. The orthogonal mixing module, the low-pass filtering module, and the modulo computing module are sequentially connected;
所述本地振荡模块用于产生方波信号,所述方波信号频率fc=2f0,f0为超声回波信号中心频率;所述方波信号经二分频后产生频率为f0、相位差为90度的两路方波载波信号F1(t)、F2(t);The local oscillation module is configured to generate a square wave signal, the square wave signal frequency f c = 2f 0 , f 0 is the ultrasonic echo signal center frequency; the square wave signal is divided by two to generate a frequency f 0 , Two square wave carrier signals F 1 (t), F 2 (t) with a phase difference of 90 degrees;
所述正交混频模块用于将超声回波信号分别与所述两路载波信号F1(t)、F2(t)相乘混频,得到两路正交信号I(t)、Q(t);The quadrature mixing module is configured to multiply and mix the ultrasonic echo signals with the two carrier signals F 1 (t) and F 2 (t) to obtain two orthogonal signals I(t) and Q. (t);
所述低通滤波模块用于滤除正交混频模块输出的两路正交信号中的高次谐波部分,得到低频信号I′(t)、Q′(t);The low pass filtering module is configured to filter out higher harmonic parts of the two orthogonal signals output by the quadrature mixing module to obtain low frequency signals I'(t), Q'(t);
所述取模运算模块用于提取所述低通滤波模块输出两路信号中的直流分量,并合成脉冲流信号,得到用于FRI稀疏采样的超声脉冲流信号A(t)。The modulo operation module is configured to extract a DC component in the two signals output by the low-pass filter module, and synthesize the pulse stream signal to obtain an ultrasonic pulse stream signal A(t) for FRI sparse sampling.
进一步优选方案,所述本地振荡模块采用变容二极管SVC321以及高速CMOS施密特反相器74HC14构成压控方波振荡电路,由分压电路提供0~10V控制电压,通过直流偏置电阻R1加到变容二极管上,并利用大电容C1将控制电压与CMOS施密特反相器输入端隔离,变容二极管电容量的变化范围为20pF~400pF,电阻R2和变容二极管的电容VC1的时间常数决定振荡信号的频率fcFurther preferably, the local oscillation module uses a varactor diode SVC321 and a high-speed CMOS Schmitt inverter 74HC14 to form a voltage-controlled square wave oscillating circuit, and the voltage-dividing circuit provides a control voltage of 0 to 10V, and is added by a DC bias resistor R1. Connect to the varactor diode and isolate the control voltage from the CMOS Schmitt inverter input with a large capacitor C1. The varactor diode capacitance varies from 20pF to 400pF, and the resistor R2 and the varactor diode VC1 are timed. The constant determines the frequency f c of the oscillating signal.
进一步优选方案,所述正交混频模块采用集成正交解调芯片RF2713实现,所述本地振荡模块产生的信号FLO(t)通过电容C2以交流耦合的方式输入到RF2713的LO输入引脚,所述RF2713内部包括一个数字分频器,将FLO(t)二分频为两路相位差为90度的方波载波信号;超声回波信号S(t)通过电容C3以交流耦合的方式输入RF2713,并通过内部两路吉尔伯特混频单元与两路方波载波信号分别进行混频,混频后的两路输出信号I(t)、Q(t)分别通过电容C7、C8以交流耦合的方式输出。Further preferably, the quadrature mixing module is implemented by using an integrated quadrature demodulation chip RF2713, and the signal F LO (t) generated by the local oscillation module is input to the LO input pin of the RF2713 through the capacitor C2 in an AC coupling manner. The RF2713 internally includes a digital frequency divider that divides F LO (t) into two square wave carrier signals with a phase difference of 90 degrees; the ultrasonic echo signal S(t) is AC coupled through capacitor C3. The method is input into RF2713, and is mixed by two internal Gilbert mixing units and two square wave carrier signals respectively, and the two output signals I(t) and Q(t) after mixing are respectively passed through capacitors C7 and C8. Output in AC coupling.
进一步优选方案,所述低通滤波模块由集成运放芯片AD847和若干电阻电容构成SallenKey结构的三阶线性相位低通滤波器。In a further preferred solution, the low-pass filter module comprises a third-order linear phase low-pass filter of a SallenKey structure by an integrated operational amplifier chip AD847 and a plurality of resistor capacitors.
进一步优选方案,所述取模运算模块包括两路平方运算电路、一路加法运算电路以及一路开平方根运算电路,所述两路平方运算电路均与所述加法运算电路相连,所述加法运算电路与所述开平方根运算电路相连;In a further preferred embodiment, the modulo operation module includes a two-way square operation circuit, an add-on operation circuit, and an open-square root operation circuit, and the two-way square operation circuits are all connected to the addition operation circuit, and the addition operation circuit and The square root computing circuit is connected;
所述平方运算电路和开平方根运算电路均由集成模拟乘法器芯片AD734和若干电阻电容实现,所述加法运算电路由集成运放芯片AD847构成同相加法器。 The square operation circuit and the square root operation circuit are both realized by an integrated analog multiplier chip AD734 and a plurality of resistor capacitors, and the adder circuit is composed of an integrated op amp chip AD847 to form an in-phase adder.
本发明的有益效果:The beneficial effects of the invention:
1、本发明异于现有的超声信号FRI采样中利用软件方法提取超声脉冲流的方法,利用硬件电路直接提取超声脉冲流信号,适合于超声信号FRI采样系统当中。1. The present invention is different from the existing ultrasonic signal FRI sampling method for extracting ultrasonic pulse flow by software method, and the ultrasonic pulse flow signal is directly extracted by the hardware circuit, and is suitable for the ultrasonic signal FRI sampling system.
2、提出的基于正交解调原理的超声脉冲流物理形成方法,能够从原始超声回波信号中准确提取出脉冲流信号。2. The proposed ultrasonic pulse flow physical formation method based on the orthogonal demodulation principle can accurately extract the pulse flow signal from the original ultrasonic echo signal.
3、该方法允许载波频率与实际回波信号中心频率有一定的偏差,降低了设计复杂度。3. The method allows the carrier frequency to have a certain deviation from the center frequency of the actual echo signal, which reduces the design complexity.
4、所设计的超声脉冲流物理形成电路,利用集成芯片组建成的模拟电路实现正交解调过程,电路结构简单易于实现。4. The designed ultrasonic pulse stream physically forms a circuit, and the quadrature demodulation process is realized by the analog circuit built by the integrated chipset, and the circuit structure is simple and easy to implement.
5、该脉冲流的硬件实现是超声信号稀疏采样物理实现的关键,在保证信号中信息完整性的同时,可极大地减少数据采集量,解决常规采样方法数据量太大的问题。5. The hardware implementation of the pulse stream is the key to the physical realization of the ultrasonic signal sparse sampling. While ensuring the information integrity in the signal, the data collection amount can be greatly reduced, and the problem that the conventional sampling method has too large data amount is solved.
附图说明DRAWINGS
图1为本发明中超声脉冲流物理形成方法原理图;1 is a schematic diagram of a method for physically forming an ultrasonic pulse stream in the present invention;
图2为本发明实施例中实际超声回波信号;2 is an actual ultrasonic echo signal in an embodiment of the present invention;
图3为现有的软件方法提取的超声脉冲流信号;3 is an ultrasonic pulse stream signal extracted by an existing software method;
图4为本发明实施例中实际电路提取的超声脉冲流信号;4 is an ultrasonic pulse stream signal extracted by an actual circuit according to an embodiment of the present invention;
图5为本发明实施例中本地振荡模块电路原理图;5 is a schematic circuit diagram of a local oscillation module according to an embodiment of the present invention;
图6为本发明实施例中正交混频模块电路原理图;6 is a schematic circuit diagram of a quadrature mixing module according to an embodiment of the present invention;
图7为本发明实施例中低通滤波模块电路原理图;7 is a schematic circuit diagram of a low pass filter module according to an embodiment of the present invention;
图8为本发明实施例中取模运算模块电路原理图。FIG. 8 is a schematic circuit diagram of a modulo operation module according to an embodiment of the present invention.
图中标记:S1-本地振荡模块,S2-正交混频模块,S3-低通滤波模块,S4-取模运算模块。Marked in the figure: S1-local oscillation module, S2-orthogonal mixing module, S3-low-pass filter module, S4-modulo operation module.
具体实施方式detailed description
下面结合附图和实施例对本发明的技术方案作进一步描述。The technical solutions of the present invention are further described below in conjunction with the accompanying drawings and embodiments.
本实施例中所用实际超声回波信号如图2所示。超声探头中心频率为5MHz,实际测定回波信号中心频率f0约为3.5MHz。The actual ultrasonic echo signal used in this embodiment is as shown in FIG. The center frequency of the ultrasonic probe is 5 MHz, and the actual measured echo center frequency f 0 is about 3.5 MHz.
根据实际回波信号的中心频率确定本地振荡器的振荡输出信号FLO(t)的频率fcThe frequency f c of the oscillation output signal F LO (t) of the local oscillator is determined according to the center frequency of the actual echo signal:
fc=2f0=7MHzf c =2f 0 =7MHz
由于实际超声回波信号中心频率f0受外界因素影响,会产生一定的频率波动,而所述振荡信号FLO(t)由于模拟电路的性能限制,与理论计算值之间也会存在一定的频率偏差。因此,由其分频后产生的实际载波信号频率与实际回波信号中心频率存在一定的频差Δf。以回波信号为参考,则载波信号频率为(f0+Δf),两路载波信号F1(t)、F2(t)的傅立叶级数展开式分别为:Since the actual ultrasonic echo signal center frequency f 0 is affected by external factors, a certain frequency fluctuation will occur, and the oscillation signal F LO (t) has a certain limit between the theoretical calculation value and the theoretical calculation value due to the performance limitation of the analog circuit. Frequency deviation. Therefore, there is a certain frequency difference Δf between the actual carrier signal frequency generated by the frequency division and the actual echo signal center frequency. Taking the echo signal as a reference, the carrier signal frequency is (f 0 +Δf), and the Fourier series expansion of the two carrier signals F 1 (t) and F 2 (t) are:
Figure PCTCN2017076691-appb-000009
Figure PCTCN2017076691-appb-000009
Figure PCTCN2017076691-appb-000010
Figure PCTCN2017076691-appb-000010
式中,K为载波幅值系数,n=1,3,5,…。Where K is the carrier amplitude coefficient, n = 1, 3, 5, ....
所述超声回波信号S(t)可以看成是由探头中心频率信号进行幅值调制的高斯脉冲流信号,可由如下数学模型近似表示:The ultrasonic echo signal S(t) can be regarded as a Gaussian pulse stream signal amplitude-modulated by the probe center frequency signal, which can be approximated by the following mathematical model:
Figure PCTCN2017076691-appb-000011
Figure PCTCN2017076691-appb-000011
式中,L为回波个数;
Figure PCTCN2017076691-appb-000012
为回波幅值系数;
Figure PCTCN2017076691-appb-000013
为回波带宽因子;τ为回波到达时刻;f0为回波信号中心频率;
Figure PCTCN2017076691-appb-000014
为初相位;A(t)表示脉冲数为L的高斯脉冲流信号,
Figure PCTCN2017076691-appb-000015
Where L is the number of echoes;
Figure PCTCN2017076691-appb-000012
Is the amplitude coefficient of the echo;
Figure PCTCN2017076691-appb-000013
Is the echo bandwidth factor; τ is the arrival time of the echo; f 0 is the center frequency of the echo signal;
Figure PCTCN2017076691-appb-000014
Is the initial phase; A(t) represents a Gaussian pulse stream signal with a pulse number of L,
Figure PCTCN2017076691-appb-000015
将所述超声回波信号S(t)与所述两路载波信号F1(t),F2(t)相乘,得到两路正交信号I(t),Q(t):Multiplying the ultrasonic echo signal S(t) by the two carrier signals F 1 (t), F 2 (t) to obtain two orthogonal signals I(t), Q(t):
Figure PCTCN2017076691-appb-000016
Figure PCTCN2017076691-appb-000016
Figure PCTCN2017076691-appb-000017
Figure PCTCN2017076691-appb-000017
因此,混频过后形成的两路信号I(t),Q(t)中分别包含低频分量
Figure PCTCN2017076691-appb-000018
以及频率为2f0+Δf以上的高次谐波分量。
Therefore, the two signals I(t) and Q(t) formed after the mixing respectively contain low frequency components.
Figure PCTCN2017076691-appb-000018
And higher harmonic components with a frequency of 2f 0 + Δf or more.
通过低通滤波将I(t),Q(t)中频率为2f0+Δf以上的高次谐波分量滤除,只保留低频分量,从而得到低频正交信号I′(t),Q′(t)。Low-pass filtering filters out higher-order harmonic components with frequencies above 2f 0 +Δf in I(t) and Q(t), leaving only low-frequency components, resulting in low-frequency quadrature signals I'(t), Q' (t).
滤波器参数的设定,可根据下式确定:The setting of the filter parameters can be determined according to the following formula:
Figure PCTCN2017076691-appb-000019
Figure PCTCN2017076691-appb-000019
其中,f0表示超声回波信号中心频率;fp表示滤波器通带截止频率;fs表示滤波器阻带截止频率;BWA(t)表示超声脉冲流信号A(t)的带宽。Where f 0 represents the center frequency of the ultrasonic echo signal; f p represents the cut-off frequency of the filter pass band; f s represents the cut-off frequency of the filter stop band; BW A(t) represents the bandwidth of the ultrasonic pulse flow signal A(t).
滤除高频分量后,形成两路所述低频正交信号I′(t),Q′(t),其表达式为:After filtering out the high frequency component, two low frequency quadrature signals I'(t), Q'(t) are formed, and the expression is:
Figure PCTCN2017076691-appb-000020
Figure PCTCN2017076691-appb-000020
Figure PCTCN2017076691-appb-000021
Figure PCTCN2017076691-appb-000021
对I′(t),Q′(t)分别平方后相加,再进行开平方根运算:For I'(t), Q'(t) are squared and added separately, and then the square root operation is performed:
Figure PCTCN2017076691-appb-000022
Figure PCTCN2017076691-appb-000022
由此得到超声脉冲流信号A(t)。The ultrasonic pulse stream signal A(t) is thus obtained.
可见,通过合理设计滤波器的参数,可以使得由频差Δf引入的误差在取模运算过程中被消除掉。因此该方法允许载波信号频率与回波信号中心频率存在一定的频差,这为硬件电路的设计放宽了限制。It can be seen that by reasonably designing the parameters of the filter, the error introduced by the frequency difference Δf can be eliminated in the modulo operation. Therefore, the method allows a certain frequency difference between the carrier signal frequency and the echo signal center frequency, which relaxes the limitation of the design of the hardware circuit.
如图1所示,为本发明中超声脉冲流物理形成方法原理图。包括:本地振荡 模块S1、正交混频模块S2、低通滤波模块S3、取模运算模块S4;所述本地振荡模块的输出经二分频后与所述正交混频模块相连接,所述正交混频模块、所述低通滤波模块、所述取模运算模块依次相连接;As shown in FIG. 1, it is a schematic diagram of a physical formation method of ultrasonic pulse flow in the present invention. Including: local oscillation a module S1, a quadrature mixing module S2, a low-pass filtering module S3, and a modulo operation module S4; the output of the local oscillation module is connected to the orthogonal mixing module after being divided by two, the orthogonal mixing The frequency module, the low-pass filter module, and the modulo operation module are sequentially connected;
所述本地振荡模块用于产生方波信号,所述方波信号频率fc=2f0,f0为超声回波信号中心频率;所述方波信号经二分频后产生频率为f0、相位差为90度的两路方波载波信号F1(t)、F2(t);The local oscillation module is configured to generate a square wave signal, the square wave signal frequency f c = 2f 0 , f 0 is the ultrasonic echo signal center frequency; the square wave signal is divided by two to generate a frequency f 0 , Two square wave carrier signals F 1 (t), F 2 (t) with a phase difference of 90 degrees;
所述正交混频模块用于将超声回波信号分别与所述两路载波信号F1(t)、F2(t)相乘混频,得到两路正交信号I(t)、Q(t);The quadrature mixing module is configured to multiply and mix the ultrasonic echo signals with the two carrier signals F 1 (t) and F 2 (t) to obtain two orthogonal signals I(t) and Q. (t);
所述低通滤波模块用于滤除正交混频模块输出的两路正交信号中的高次谐波部分,得到低频信号I′(t)、Q′(t);The low pass filtering module is configured to filter out higher harmonic parts of the two orthogonal signals output by the quadrature mixing module to obtain low frequency signals I'(t), Q'(t);
所述取模运算模块用于提取所述低通滤波模块输出两路信号中的直流分量,并合成脉冲流信号,得到用于FRI稀疏采样的超声脉冲流信号A(t)。The modulo operation module is configured to extract a DC component in the two signals output by the low-pass filter module, and synthesize the pulse stream signal to obtain an ultrasonic pulse stream signal A(t) for FRI sparse sampling.
本发明中的本地振荡模块S1,如图5所示。产生频率为超声回波信号中心频率两倍的方波振荡信号。该模块采用变容二极管SVC321以及高速CMOS施密特反相器74HC14构成压控方波振荡电路。由分压电路提供0~10V控制电压,通过直流偏置电阻R1加到变容二极管上,并利用大电容C1将控制电压与CMOS芯片输入端隔离。变容二极管电容量的变化范围约为20pF~400pF,电阻R2和变容二极管的电容VC1的时间常数决定振荡信号的频率,振荡频率fc可由下式近似表示:The local oscillation module S1 in the present invention is as shown in FIG. A square wave oscillating signal having a frequency twice the center frequency of the ultrasonic echo signal is generated. The module uses a varactor diode SVC321 and a high speed CMOS Schmitt inverter 74HC14 to form a voltage controlled square wave oscillating circuit. A voltage of 0 to 10 V is supplied from the voltage dividing circuit, and is applied to the varactor through a DC bias resistor R1, and the control voltage is isolated from the input end of the CMOS chip by the large capacitor C1. The capacitance of the varactor diode ranges from about 20pF to 400pF. The time constant of the resistor R2 and the capacitance VC1 of the varactor determines the frequency of the oscillating signal. The oscillating frequency f c can be approximated by the following equation:
Figure PCTCN2017076691-appb-000023
Figure PCTCN2017076691-appb-000023
本发明中的正交混频模块S2,如图6所示。对本地振荡信号FLO(t)二分频形成两路相位差90度的方波载波信号,并将原始超声信号与两路载波信号进行混频,输出两路正交信号I(t),Q(t)。该模块采用集成正交解调芯片RF2713来实现其功能。本地振荡信号FLO(t)通过电容C2以交流耦合的方式输入到RF2713的LO输入引脚,其内部包括一个数字分频器,将FLO(t)二分频为两路相位差为90度的方波载波信号。超声回波信号S(t)通过电容C3以交流耦合的方式输入RF2713,并通过内部两路吉尔伯特混频单元与两路载波信号分别进行混频。两路输出信号I(t),Q(t)分别通过电容C7,C8以交流耦合的方式输出。The orthogonal mixing module S2 in the present invention is as shown in FIG. The local oscillation signal F LO (t) is divided by two to form two square wave carrier signals with a phase difference of 90 degrees, and the original ultrasonic signal is mixed with the two carrier signals to output two orthogonal signals I(t). Q(t). The module uses the integrated quadrature demodulation chip RF2713 to achieve its function. The local oscillating signal F LO (t) is ac-coupled to the LO input pin of the RF2713 via capacitor C2. It internally includes a digital divider that divides F LO (t) by two to a phase difference of 90. Square wave carrier signal. The ultrasonic echo signal S(t) is input to the RF2713 by means of a capacitor C3 in an AC-coupled manner, and is separately mixed with the two carrier signals by an internal two-way Gilbert mixing unit. The two output signals I(t) and Q(t) are outputted by AC coupling through capacitors C7 and C8, respectively.
本发明中的低通滤波模块S3,如图7所示。用于滤除正交混频模块S2中输 出的两路正交信号I(t),Q(t)中高次谐波分量。该模块包括两路参数一致的滤波电路,如(a)(b)所示,分别实现对I(t),Q(t)的滤波。每一路采用集成运放芯片AD847构成SallenKey结构的三阶线性相位低通滤波器。The low pass filter module S3 in the present invention is as shown in FIG. Used to filter out the input of the quadrature mixing module S2 The two orthogonal signals I(t), Q(t) are the higher harmonic components. The module includes two filtering circuits with the same parameters, as shown in (a) and (b), respectively filtering the I(t), Q(t). Each channel uses an integrated operational amplifier chip AD847 to form a third-order linear phase low-pass filter of the SallenKey structure.
本发明中的取模运算模块S4,如图8所示。用于提取S3低通滤波模块输出两路信号中的直流分量,并合成脉冲流信号。该模块包括两路平方运算电路,一路加法运算电路以及开平方根运算电路。其中平方运算和开平方根运算电路采用集成模拟乘法器芯片AD734构成,加法运算电路采用集成运放芯片AD847构成同相加法器。The modulo operation module S4 in the present invention is as shown in FIG. It is used to extract the DC component of the two signals outputted by the S3 low-pass filter module, and synthesize the pulse stream signal. The module includes two square operation circuits, one addition circuit and an square root operation circuit. The square operation and the square root operation circuit are composed of an integrated analog multiplier chip AD734, and the addition circuit uses an integrated operational amplifier chip AD847 to form an in-phase adder.
利用本发明提出的超声脉冲流物理形成方法与电路从图2所示超声回波信号中提取的超声脉冲流信号如图4所示。图3为利用现有的正交解调软件方法从图2所示超声回波信号中提取的超声脉冲流信号。在现有的脉冲流形成算法中,正交解调算法的精度最高,可以满足超声信号稀疏采样中对脉冲流的波形要求。从图3和图4的对比中可以看出,通过本发明中的硬件电路提取的脉冲流与正交解调算法中提取的脉冲流波形非常接近,能满足硬件稀疏采样中对脉冲流的波形要求。并且避免了现有算法需要先按常规采样得到信号,再通过算法形成脉冲流的环节,实现了脉冲流的直接硬件形成,具有实时性。本发明是物理实现超声信号稀疏采样的关键。The ultrasonic pulse flow signal extracted from the ultrasonic echo signal shown in FIG. 2 by the ultrasonic pulse flow physical formation method and circuit proposed by the present invention is as shown in FIG. 4. 3 is an ultrasonic pulse stream signal extracted from the ultrasonic echo signal shown in FIG. 2 using a conventional quadrature demodulation software method. In the existing pulse stream formation algorithm, the orthogonal demodulation algorithm has the highest precision, which can meet the waveform requirements of the pulse stream in the sparse sampling of ultrasonic signals. It can be seen from the comparison between FIG. 3 and FIG. 4 that the pulse stream extracted by the hardware circuit of the present invention is very close to the pulse stream waveform extracted by the orthogonal demodulation algorithm, and can satisfy the waveform of the pulse stream in the hardware sparse sampling. Claim. Moreover, the existing algorithm needs to first obtain the signal according to the conventional sampling, and then form the pulse flow through the algorithm, thereby realizing the direct hardware formation of the pulse flow, and has real-time performance. The invention is the key to physically achieving sparse sampling of ultrasonic signals.
应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,本领域技术人员应当将说明书作为一个整体,各实施方式中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。It should be understood that, although the present specification is described in terms of embodiments, not every embodiment includes only one independent technical solution, and those skilled in the art should have the specification as a whole, and the technical solutions in the respective embodiments may also be combined as appropriate. Other embodiments that can be understood by those skilled in the art are formed.
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。 The series of detailed descriptions set forth above are merely illustrative of the possible embodiments of the present invention, and are not intended to limit the scope of the present invention. Changes are intended to be included within the scope of the invention.

Claims (9)

  1. 一种超声信号稀疏采样中脉冲流形成的方法,其特征在于,包括如下步骤:A method for forming a pulse stream in ultrasonic signal sparse sampling, comprising the steps of:
    步骤1,由本地振荡器产生固定频率的方波振荡信号FLO(t),其振荡频率fc满足:fc=2f0;其中,f0为超声回波信号中心频率;Step 1: The local oscillator generates a square wave oscillating signal F LO (t) of a fixed frequency, and the oscillating frequency f c satisfies: f c = 2f 0 ; wherein f 0 is the center frequency of the ultrasonic echo signal;
    步骤2,将步骤1中所述振荡信号FLO(t)进行二分频,形成两路幅值相等、频率相同、相位差为90度的方波载波信号,分别为F1(t)、F2(t);Step 2, the oscillating signal F LO (t) in step 1 is divided by two to form two square wave carrier signals with equal amplitude, same frequency and phase difference of 90 degrees, respectively F 1 (t), F 2 (t);
    步骤3,将接收到的中心频率为f0的超声回波信号S(t)分别与步骤2中所述两路方波信号F1(t)、F2(t)相乘混频,对超声回波信号进行频谱搬移,形成两路包含低频分量和高次谐波分量的正交信号I(t)、Q(t);Step 3, multiplying and receiving the ultrasonic echo signals S(t) having the center frequency f 0 and the two square wave signals F 1 (t) and F 2 (t) in step 2, respectively, The ultrasonic echo signal is spectrally shifted to form two orthogonal signals I(t), Q(t) including low frequency components and higher harmonic components;
    步骤4,对步骤3中所述两路正交信号I(t)、Q(t)进行低通滤波,滤除高次谐波分量,形成两路低频分量信号I′(t)、Q′(t);Step 4, performing low-pass filtering on the two orthogonal signals I(t) and Q(t) in step 3, filtering out higher harmonic components, and forming two low-frequency component signals I'(t), Q' (t);
    步骤5,对步骤4中所述两路正交的低频分量信号I′(t),Q′(t)进行取模运算,获取用于FRI稀疏采样的超声脉冲流信号A(t),所述取模运算具体包括:对所述两路正交的低频分量信号I′(t)、Q′(t)分别平方相加,再对运算结果进行开平方根运算,得到超声脉冲流信号A(t)。Step 5, performing modulo operation on the two orthogonal low-frequency component signals I'(t), Q'(t) in step 4, and acquiring an ultrasonic pulse stream signal A(t) for FRI sparse sampling, Specifically, the modulo operation includes: square-adding the two orthogonal low-frequency component signals I′(t) and Q′(t), and performing square root operation on the operation result to obtain an ultrasonic pulse stream signal A ( t).
  2. 根据权利要求1所述的一种超声信号稀疏采样中脉冲流形成的方法,其特征在于,步骤2还包括:设实际载波信号频率与实际超声回波信号中心频率存在频差为Δf,以超声回波信号为参考,建立F1(t)、F2(t)的傅立叶级数展开式分别如下:The method for forming a pulse stream in an ultrasonic signal sparse sampling according to claim 1, wherein the step 2 further comprises: setting the frequency difference between the actual carrier signal frequency and the actual ultrasonic echo signal center frequency to Δf to the ultrasound The echo signal is used as a reference, and the Fourier series expansion of F 1 (t) and F 2 (t) is established as follows:
    Figure PCTCN2017076691-appb-100001
    Figure PCTCN2017076691-appb-100001
    Figure PCTCN2017076691-appb-100002
    Figure PCTCN2017076691-appb-100002
    式中,K为载波幅值系数,n=1,3,5,…。Where K is the carrier amplitude coefficient, n = 1, 3, 5, ....
  3. 根据权利要求1所述的一种超声信号稀疏采样中脉冲流形成的方法,其特征在于,步骤3还包括:将超声回波信号S(t)看成是由探头中心频率信号进行幅值调制的高斯脉冲流信号,并建立如下数学表达式: The method for forming a pulse stream in ultrasonic signal sparse sampling according to claim 1, wherein the step 3 further comprises: treating the ultrasonic echo signal S(t) as amplitude modulation by the probe center frequency signal. The Gaussian pulse stream signal and establishes the following mathematical expression:
    Figure PCTCN2017076691-appb-100003
    Figure PCTCN2017076691-appb-100003
    式中,L为回波个数;
    Figure PCTCN2017076691-appb-100004
    为回波幅值系数;
    Figure PCTCN2017076691-appb-100005
    为回波带宽因子;τ为回波到达时刻;f0为超声回波信号中心频率;
    Figure PCTCN2017076691-appb-100006
    为初相位;A(t)表示脉冲数为L的高斯脉冲流信号,
    Figure PCTCN2017076691-appb-100007
    Where L is the number of echoes;
    Figure PCTCN2017076691-appb-100004
    Is the amplitude coefficient of the echo;
    Figure PCTCN2017076691-appb-100005
    Is the echo bandwidth factor; τ is the arrival time of the echo; f 0 is the center frequency of the ultrasonic echo signal;
    Figure PCTCN2017076691-appb-100006
    Is the initial phase; A(t) represents a Gaussian pulse stream signal with a pulse number of L,
    Figure PCTCN2017076691-appb-100007
  4. 根据权利要求1所述的一种超声信号稀疏采样中脉冲流形成的方法,其特征在于,步骤4中所述低通滤波器的参数根据下式确定:The method for forming a pulse stream in ultrasonic signal sparse sampling according to claim 1, wherein the parameter of the low-pass filter in step 4 is determined according to the following formula:
    Figure PCTCN2017076691-appb-100008
    Figure PCTCN2017076691-appb-100008
    其中,f0表示超声回波信号中心频率;fp表示滤波器通带截止频率;fs表示滤波器阻带截止频率;BWA(t)表示超声脉冲流信号A(t)的带宽。Where f 0 represents the center frequency of the ultrasonic echo signal; f p represents the cut-off frequency of the filter pass band; f s represents the cut-off frequency of the filter stop band; BW A(t) represents the bandwidth of the ultrasonic pulse flow signal A(t).
  5. 一种超声信号稀疏采样中脉冲流形成的电路,其特征在于,包括:本地振荡模块S1、正交混频模块S2、低通滤波模块S3、取模运算模块S4;所述本地振荡模块的输出经二分频后与所述正交混频模块相连接,所述正交混频模块、所述低通滤波模块、所述取模运算模块依次相连接;A circuit for forming a pulse stream in an ultrasonic signal sparse sampling, comprising: a local oscillation module S1, a quadrature mixing module S2, a low-pass filtering module S3, a modulo operation module S4, and an output of the local oscillation module After being divided by two, the orthogonal mixing module is connected, and the orthogonal mixing module, the low-pass filtering module, and the modulo computing module are sequentially connected;
    所述本地振荡模块用于产生方波信号,所述方波信号频率fc=2f0,f0为超声回波信号中心频率;所述方波信号经二分频后产生频率为f0、相位差为90度的两路方波载波信号F1(t)、F2(t);The local oscillation module is configured to generate a square wave signal, the square wave signal frequency f c = 2f 0 , f 0 is the ultrasonic echo signal center frequency; the square wave signal is divided by two to generate a frequency f 0 , Two square wave carrier signals F 1 (t), F 2 (t) with a phase difference of 90 degrees;
    所述正交混频模块用于将超声回波信号分别与所述两路载波信号F1(t)、F2(t)相乘混频,得到两路正交信号I(t)、Q(t);The quadrature mixing module is configured to multiply and mix the ultrasonic echo signals with the two carrier signals F 1 (t) and F 2 (t) to obtain two orthogonal signals I(t) and Q. (t);
    所述低通滤波模块用于滤除正交混频模块输出的两路正交信号中的高次谐波部分,得到低频信号I′(t)、Q′(t);The low pass filtering module is configured to filter out higher harmonic parts of the two orthogonal signals output by the quadrature mixing module to obtain low frequency signals I'(t), Q'(t);
    所述取模运算模块用于提取所述低通滤波模块输出两路信号中的直流分量,并合成脉冲流信号,得到用于FRI稀疏采样的超声脉冲流信号A(t)。The modulo operation module is configured to extract a DC component in the two signals output by the low-pass filter module, and synthesize the pulse stream signal to obtain an ultrasonic pulse stream signal A(t) for FRI sparse sampling.
  6. 根据权利要求5所述的一种超声信号稀疏采样中脉冲流形成的电路,其特征在于,所述本地振荡模块采用变容二极管SVC321以及高速CMOS施密特反相器74HC14构成压控方波振荡电路,由分压电路提供0~10V控制电压,通过直流偏置电阻R1加到变容二极管上,并利用大电容C1将控制电压与CMOS施密特反相器输入端隔离,变容二极管电容量的变化范围为20pF~400pF,电阻 R2和变容二极管的电容VC1的时间常数决定振荡信号的频率fcThe circuit for forming a pulse stream in ultrasonic signal sparse sampling according to claim 5, wherein the local oscillation module comprises a varactor diode SVC321 and a high speed CMOS Schmitt inverter 74HC14 to form a voltage controlled square wave oscillation. The circuit is provided with a 0~10V control voltage from the voltage dividing circuit, and is applied to the varactor through the DC bias resistor R1, and the control voltage is isolated from the input end of the CMOS Schmitt inverter by the large capacitor C1, and the varactor diode is electrically connected. The capacity varies from 20 pF to 400 pF, and the time constant of the resistor R2 and the capacitance VC1 of the varactor determines the frequency f c of the oscillating signal.
  7. 根据权利要求5所述的一种超声信号稀疏采样中脉冲流形成的电路,其特征在于,所述正交混频模块采用集成正交解调芯片RF2713实现,所述本地振荡模块产生的信号FLO(t)通过电容C2以交流耦合的方式输入到RF2713的LO输入引脚,所述RF2713内部包括一个数字分频器,将FLO(t)二分频为两路相位差为90度的方波载波信号;超声回波信号S(t)通过电容C3以交流耦合的方式输入RF2713,并通过内部两路吉尔伯特混频单元与两路方波载波信号分别进行混频,混频后的两路输出信号I(t)、Q(t)分别通过电容C7、C8以交流耦合的方式输出。The circuit for forming a pulse stream in the ultrasonic signal sparse sampling according to claim 5, wherein the orthogonal mixing module is implemented by using an integrated quadrature demodulation chip RF2713, and the signal generated by the local oscillation module is F LO (t) is ac-coupled to the LO input pin of RF2713 via capacitor C2. The RF2713 internally includes a digital divider that divides F LO (t) by two to a phase difference of 90 degrees. The square wave carrier signal; the ultrasonic echo signal S(t) is input to the RF2713 by means of the capacitor C3 in an AC coupling manner, and is separately mixed by the internal two-way Gilbert mixing unit and the two square wave carrier signals, after mixing The two output signals I(t) and Q(t) are outputted by AC coupling through capacitors C7 and C8, respectively.
  8. 根据权利要求5所述的一种超声信号稀疏采样中脉冲流形成的电路,其特征在于,所述低通滤波模块由集成运放芯片AD847和若干电阻电容构成SallenKey结构的三阶线性相位低通滤波器。The circuit for forming a pulse stream in ultrasonic signal sparse sampling according to claim 5, wherein the low-pass filter module comprises a third-order linear phase low-pass of the SallenKey structure by an integrated operational amplifier chip AD847 and a plurality of resistor capacitors. filter.
  9. 根据权利要求5所述的一种超声信号稀疏采样中脉冲流形成的电路,其特征在于,所述取模运算模块包括两路平方运算电路、一路加法运算电路以及一路开平方根运算电路,所述两路平方运算电路均与所述加法运算电路相连,所述加法运算电路与所述开平方根运算电路相连;The circuit for forming a pulse stream in ultrasonic signal sparse sampling according to claim 5, wherein the modulo operation module comprises two square operation circuits, one add operation circuit and one open square root operation circuit, Two square operation circuits are connected to the addition circuit, and the addition circuit is connected to the square root operation circuit;
    所述平方运算电路和开平方根运算电路均由集成模拟乘法器芯片AD734和若干电阻电容实现,所述加法运算电路由集成运放芯片AD847构成同相加法器。 The square operation circuit and the square root operation circuit are both realized by an integrated analog multiplier chip AD734 and a plurality of resistor capacitors, and the adder circuit is composed of an integrated op amp chip AD847 to form an in-phase adder.
PCT/CN2017/076691 2016-04-08 2017-03-15 Method and circuit for forming pulse stream during sparse sampling of ultrasonic signal WO2017173913A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201610219295.0A CN105738885B (en) 2016-04-08 2016-04-08 The method and circuit that stream of pulses is formed in a kind of ultrasonic signal sparse sampling
CN201610219295.0 2016-04-08

Publications (1)

Publication Number Publication Date
WO2017173913A1 true WO2017173913A1 (en) 2017-10-12

Family

ID=56253959

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2017/076691 WO2017173913A1 (en) 2016-04-08 2017-03-15 Method and circuit for forming pulse stream during sparse sampling of ultrasonic signal

Country Status (2)

Country Link
CN (1) CN105738885B (en)
WO (1) WO2017173913A1 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105738885B (en) * 2016-04-08 2018-04-24 江苏大学 The method and circuit that stream of pulses is formed in a kind of ultrasonic signal sparse sampling
CN108896980B (en) * 2018-04-09 2022-06-21 飞依诺科技股份有限公司 Method and device for generating RF data based on IQ data
CN110261492B (en) * 2019-05-14 2021-11-23 江苏大学 Method and device for driving ultrasonic sensor by digital sequence
CN110501429B (en) * 2019-07-24 2022-05-20 江苏大学 Sparse sampling method for array ultrasonic signals
CN110808753B (en) * 2019-11-13 2021-02-02 中国矿业大学 Three-path signal aliasing processing method applied to band-pass sampling
CN111130513B (en) * 2019-12-13 2024-01-05 江苏大学 Encoding ultrasonic signal pulse stream forming method and circuit based on high-order moment
CN113037303B (en) * 2021-03-02 2022-04-12 清华大学 Quasi-coherent pulse ultra-wideband receiver and signal demodulation method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102195617A (en) * 2010-03-14 2011-09-21 技术研究及发展基金有限公司 Low-rate sampling of pulse streams
US20130038479A1 (en) * 2010-04-27 2013-02-14 Technion Research & Development Foundation Ltd. Multi-Channel Sampling of Pulse Streams at the Rate of Innovation
CN105738885A (en) * 2016-04-08 2016-07-06 江苏大学 Method and circuit for forming pulse flow in ultrasonic signal sparse sampling

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104931584B (en) * 2015-05-08 2017-08-15 哈尔滨工业大学 Based on the ultrasound computed tomography detection method that compression sampling is theoretical
CN104914440B (en) * 2015-06-23 2017-06-16 重庆大学 A kind of compressed sensing ultrasonic beam synthetic method for merging dynamic aperture

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102195617A (en) * 2010-03-14 2011-09-21 技术研究及发展基金有限公司 Low-rate sampling of pulse streams
US20130038479A1 (en) * 2010-04-27 2013-02-14 Technion Research & Development Foundation Ltd. Multi-Channel Sampling of Pulse Streams at the Rate of Innovation
CN105738885A (en) * 2016-04-08 2016-07-06 江苏大学 Method and circuit for forming pulse flow in ultrasonic signal sparse sampling

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CAO, WEN ET AL.: "Design of Digital I/Q Demodulation System in Echoed Ultrasonic Waves Based on FPGA", JOURNAL OF SOUTHWEST UNIVERSITY OF SCIENCE AND TECHNOLOGY, vol. 21, no. 3, 30 September 2006 (2006-09-30), pages 57 - 59, ISSN: 1671-8755 *
CHANG, J.H. ET AL.: "ANovel Envelope Detector for High-Frame Rate, High-Frequency Ultrasound Imaging", IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL, vol. 54, no. 9, 30 September 2007 (2007-09-30), XP055428994, ISSN: 0885-3010 *
DEMIRLI, R. ET AL.: "Model-Based Estimation of Ultrasonic Echoes Part I: Analysis and Algorithms", IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL, vol. 48, no. 3, 31 May 2001 (2001-05-31), pages 789, XP011438299, ISSN: 0885-3010 *
HUANG, ZHIWEI, MODEM CIRCUIT DESIGN, 30 April 2009 (2009-04-30), pages 310 - 313 *
PENG, CHENGQING: "Research on Pipeline Ultrasonic Testing Signal Sampling Method and Application Based on Finite Rate of Innovation", WANFANG DATA, 12 October 2015 (2015-10-12), pages 49 - 52 *

Also Published As

Publication number Publication date
CN105738885A (en) 2016-07-06
CN105738885B (en) 2018-04-24

Similar Documents

Publication Publication Date Title
WO2017173913A1 (en) Method and circuit for forming pulse stream during sparse sampling of ultrasonic signal
CN104122444B (en) All-digital IF spectrum analyzer and frequency spectrum analysis method
CN106645949A (en) Heterodyne sweep-frequency type spectrum analyzer based on low frequency detection
CN111766597B (en) Frequency measuring method and device for frequency-modulated continuous wave laser ranging beat signal
CN105259410B (en) The frequency estimating methods and its device of a kind of lack sampling waveform under very noisy interference
CN205427032U (en) Weak signal detection device based on phase -locked loop
CN103983849B (en) A kind of Electric Power Harmonic Analysis method of real-time high-precision
CN112067130A (en) Amplitude extraction method and system for multi-frequency components
CN106353594A (en) Spectrum analyzer system with fast multi-resolution and method
CN104677486B (en) The aero-engine vibration signal Method for Phase Difference Measurement reconstructed based on tacho-pulse
CN107807276B (en) A kind of Real-time Spectrum Analyzer tracking source design method
CN105471525B (en) A kind of four-way compressed sensing digital receiver signal processing method of vector network analyzer
CN205670191U (en) A kind of echo wave signal acquisition device of LDV technique
CN107765085A (en) A kind of signal supervisory instrument and method
CN109633243B (en) Beam signal peak amplitude accurate extraction method based on multiphase sampling
CN110161310B (en) Weak signal detection method based on difference frequency modulation phase locking
CN113627398B (en) Signal characteristic detection method based on adaptive reconstruction filtering
CN107091956A (en) A kind of method and circuit for measuring network phase shift parameters
Shoupeng et al. Quadrature demodulation based circuit implementation of pulse stream for ultrasonic signal FRI sparse sampling
CN105515694A (en) Method for implementing real-time FFT frame averaging of digital signal processing
CN112666392A (en) Carrier frequency measuring circuit and method of high-speed pulse modulation signal
Karelin et al. An FPGA-based fourier FFTS-160 spectrometer for atmospheric molecular radiation research
CN104198813A (en) Device and method for measuring impedance angle of ultrasonic energy transducer through orthogonal correlation method
CN110794208A (en) Broadband rapid modulation domain frequency measurement method based on digital phase calculation
Ma et al. A signal processing practice method based on LFMCW radar

Legal Events

Date Code Title Description
NENP Non-entry into the national phase

Ref country code: DE

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17778577

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 17778577

Country of ref document: EP

Kind code of ref document: A1