CN109412763B - Digital signal existence detection method based on signal energy-entropy ratio - Google Patents

Digital signal existence detection method based on signal energy-entropy ratio Download PDF

Info

Publication number
CN109412763B
CN109412763B CN201811356297.XA CN201811356297A CN109412763B CN 109412763 B CN109412763 B CN 109412763B CN 201811356297 A CN201811356297 A CN 201811356297A CN 109412763 B CN109412763 B CN 109412763B
Authority
CN
China
Prior art keywords
signal
frame
energy
entropy
data
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.)
Active
Application number
CN201811356297.XA
Other languages
Chinese (zh)
Other versions
CN109412763A (en
Inventor
曹蕾
董彬虹
张存林
陈延涛
李芊饶
赵宇轩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201811356297.XA priority Critical patent/CN109412763B/en
Publication of CN109412763A publication Critical patent/CN109412763A/en
Application granted granted Critical
Publication of CN109412763B publication Critical patent/CN109412763B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/20Arrangements for detecting or preventing errors in the information received using signal quality detector
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/18Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being spectral information of each sub-band
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/21Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being power information
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/78Detection of presence or absence of voice signals
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/78Detection of presence or absence of voice signals
    • G10L25/87Detection of discrete points within a voice signal
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/78Detection of presence or absence of voice signals
    • G10L2025/783Detection of presence or absence of voice signals based on threshold decision

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Computational Linguistics (AREA)
  • Health & Medical Sciences (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Quality & Reliability (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

The invention discloses a digital signal existence detection method based on signal energy-entropy ratio, belonging to the field of signal detection of wireless communication systems, and relating to a digital signal existence detection method which is low in calculation complexity and does not need any priori knowledge, wherein the method can be used as a frame synchronization method of non-cooperative communication, is used for solving the problem of coarse synchronization of a receiving end signal under the condition of unknown frame structure of a sender, replaces frame synchronization, and performs coarse extraction of information; the invention discloses a digital signal existence detection method based on signal energy-entropy ratio. Since the CME iterative algorithm proposed by the background art requires multiple iterations to converge, it will bring about extremely high computational complexity. Therefore, the invention provides the method for judging whether a signal exists in one frame by using the energy-entropy ratio as the measurement index of the existence of the signal and comparing the energy-entropy ratio with the threshold value through setting the threshold value, so that the iteration is not needed, the complexity of calculation is simplified, and the accuracy of the detection of the existence of the signal is higher under the condition of higher signal-to-noise ratio.

Description

Digital signal existence detection method based on signal energy-entropy ratio
Technical Field
The invention belongs to the field of signal detection of a wireless communication system, and relates to a digital signal existence detection method which is low in calculation complexity and does not need any priori knowledge.
Background
Due to the multipath effect and Doppler frequency shift influence in short wave channel transmission, signal fading is obvious, and the diversity combining technology has obvious advantages in resisting fading. Before combining the multi-channel signals, the multi-channel signals are synchronized, and only the useful signals are aligned to the maximum extent, so that a more remarkable combining effect can be achieved. It is therefore important to find the approximate location of the desired signal from the faded single-channel signal. The traditional digital signal detection method generally adopts energy detection, and provides a CME interference detection method in the literature PERTTI HENTTU, SAMI AROMAA. Consequential Mean evaluation Algorithm [ C ], and the method adopts an iterative idea to estimate the signal-to-noise ratio by multiple iterations so as to separate noise and interference signals. The algorithm has the defects of high complexity of multiple iterations, uncertain convergence, experience-based setting of threshold parameters and poor robustness. The endpoint detection is an extremely important method in speech signal processing, is used for determining the starting and ending positions of voice, and is developed more mature. The end point detection method of the voice signal utilizes the characteristics of voice, such as short-time energy of unvoiced consonants, short-time zero crossing rate and a specific design method of difference of signal spectrum entropy to distinguish silent sections without signals from voice sections, and in the sixth chapter of book MATLAB application in voice signal analysis and synthesis, several end point detection methods of the voice signal are taught, and the detection method based on energy-entropy ratio is particularly good in detection effect of signals under fading channels. The disadvantage of this scheme is that the sum of the squares of the amplitudes of all data in one frame is taken as the energy in the energy-entropy ratio, and the energy also contains the energy of noise, so that the error caused by low signal-to-noise ratio is large. The threshold setting uses double thresholds, which are helpful for voice signal detection, but also for digital signals, they are analyzed specifically. Due to the time variability of the voice signal, the short-time analysis can be performed only assuming that the short time (about 30ms) is unchanged, but for the digital signal, the requirement of the short time can be relaxed, that is, different frame lengths can be selected according to specific situations when the frames are divided.
Disclosure of Invention
The signal detection thought of the invention is derived from an end point detection method of a voice signal, the voice signal can also be regarded as a digital signal, the method of end point detection based on the energy-entropy ratio in the voice signal is also suitable for the existence detection of the digital signal, and only the spectral line energy near the modulation frequency point can be taken as the energy in the energy-entropy ratio under the condition that the modulation frequency point of the modulation signal is known, so that the accuracy is improved. Because the performance evaluation modes of the digital signal and the voice signal are different, little noise can be allowed at the head and the tail of a signal segment for the detection of the voice signal, but the length of a noise segment at the head and the tail of the digital signal is not allowed to exceed 1 symbol length, otherwise, error codes can be caused. The invention adopts twice signal existence detection when detecting digital signal end point, firstly, the frame length of a little bit larger can be selected for detection to determine the approximate position of the signal, thus ensuring no signal loss, secondly, the frame length is further reduced on the basis of the signal extracted for the first time, and the position of starting and stopping is further accurately combined with a merging algorithm. For the end point detection of digital signals, the threshold can be set by double thresholds or single thresholds (two thresholds of the double thresholds are set to be the same, namely the single threshold), the double thresholds are better for the first coarse signal detection, and the signal loss is avoided; the single threshold effect is superior to the double thresholds in the second fine detection, and the noise from the head to the tail can be removed as much as possible. Simulation shows that the method also has a good detection effect, is low in complexity and can accurately detect the position of the digital signal.
The technical scheme of the invention is as follows:
framing the received signal, and setting an overlapped section with a half frame length for smooth transition between frames; performing Fast Fourier Transform (FFT) on the framed data in each frame to obtain a magnitude spectrum of each frame, and obtaining the energy and the spectrum entropy value of each frame of data according to the magnitude spectrum, wherein the ratio of the energy to the spectrum entropy value is the energy entropy ratio; the energy entropy ratio is used as a basis for judging whether one frame of data is noise or digital signal; the invention adopts double-threshold judgment, so that the judgment of the existence of the signal is more accurate, and the double-threshold is more favorable for finding out the start-stop position of the signal, thereby reducing the error of the start-stop position judgment; generally speaking, the method for finding the starting and ending positions of a digital signal segment from a segment of data containing noise and having discontinuous signal segments is to traverse all frames, find out the starting frame and the ending frame of a segment of digital signal by comparing the energy-entropy ratio and the threshold value of each frame data, and then continue the search of the next signal segment; finally, recording the starting and ending positions of each digital signal segment in the segment of data;
the invention relates to a digital signal existence detection method based on signal energy-entropy ratio, which comprises the following steps:
step 1, sampling a received signal r (t) by a sampling frequency Fs to obtain r (N), setting the data length of r (N) as N, and framing r (N); the specific framing method and parameter setting are as follows: taking t seconds of data as a frame, namely the frame length wlen is t multiplied by Fs; the overlap ratio between the frame and the interframe is set to be n%, that is, the frame shift (the number of moving sampling points of the next frame relative to the current frame) is inc ═ wlen × n%; total number of minutes of N long data
Figure GDA0002806493260000021
So that the received sequence of length N, r (N), becomes [ r ] after framing1(n),r2(n),...,ri(n),...,rfn(n)]Wherein each element represents a frame of data;
step 2, performing Fast Fourier Transform (FFT) on each frame of data, and obtaining a modulus value to obtain an energy spectrum:
Yi(k),i=1,2,...,fn;k=1,2,...,wlen;
and 3, calculating the energy of each frame of data, wherein the time domain energy calculating method comprises the following steps: el (electro luminescence)i=log10(1+AMPiA), i ═ 1,2,. fn, where
Figure GDA0002806493260000031
1,2, fn, a is a set constant; due to the presence of a, when AMP is usediAmplitude hasWill be in EL with larger variationiThe fading is eased, so that signals and noise at the serious fading position can be more easily distinguished; or the energy can also be obtained from the frequency domain, and for the received data of the known modulation frequency point, the energy is calculated as follows:
Figure GDA0002806493260000032
1,2.., fn, where Φ is the set of several locations near where the modulation bins are located.
Step 4, calculating the spectrum entropy value of each frame data: the normalized spectral probability density function for each frequency component within a frame of data is:
Figure GDA0002806493260000033
1,2, fn; 1,2., wlen; the spectral entropy is:
Figure GDA0002806493260000034
i=1,2,...,fn;
for noise, the distribution of the normalized spectrum probability density function is uniform, so the spectrum entropy value is larger; for digital signals, due to the fact that energy is concentrated in modulation, the normalized spectral probability density function of the digital signals is not uniformly distributed, and therefore the spectral entropy of the digital signals is generally lower than that of noise;
step 5, calculating the energy-entropy ratio of each frame of data:
Figure GDA0002806493260000035
i=1,2,...,fn;
step 6, setting double thresholds: the setting of the threshold is related to the energy-entropy ratio of each frame, and the energy-entropy ratio is obtained through calculation; the method comprises the following specific steps: solving the maximum value Me and the initial average value eth of the energy-entropy ratios of all the frames, wherein eth is the average value of the energy-entropy ratios of the previous NIS frames, and NIS is a constant and represents the length of the noise frame before the signal segment starts; setting the threshold Th1 as a noise threshold: th1 ═ alpha1X (Me-eth) + eth, setting threshold Th2 as the signal threshold: th2 ═ alpha2×(Me-eth)+eth,α12Is constant and alpha1<α2
Step 7, setting the maximum allowable noise frame length MaxSilence and the minimum signal frame length MinLen, and starting the end point detection of the signal section;
and 8: detecting an end point of a signal segment;
step 8.1, initialization: the signal frame Count is 0, the status flag Case is 0, and the noise frame Count is 0;
step 8.2, reading the energy-entropy ratio of one frame of data, and if Case is 0 or 1, turning to step 8.3; if Case is 2, go to step 8.6;
step 8.3, comparing the energy-entropy ratio with the signal threshold Th 2; if the current frame number is greater than Th2, the signal is detected, the Count is increased by 1, the current frame number is recorded as the initial frame position of the signal, Case is set to be 2, and the step is switched to step 8.2; otherwise, turning to step 8.4;
8.4, comparing the energy entropy ratio with the noise threshold Th1, if the energy entropy ratio is larger than Th1, indicating that a signal possibly exists, increasing the Count by 1, setting the Case to be 1, and turning to the step 8.2; otherwise, turning to step 8.5;
step 8.5, switching to step 8.1 when no signal is detected;
8.6, comparing the magnitude of the energy entropy ratio and the noise threshold Th1, wherein the magnitude is larger than Th1, and the data of the frame is still in the signal section, the Count is increased by 1, the Case is unchanged, and the step 8.2 is switched; otherwise, turning to step 8.7;
step 8.7, noise exists, and the Silence is increased by 1; judging whether the Silence reaches the maximum allowable noise frame length, if not, continuing to serve the current frame as a part of the signal section, increasing the Count by 1, keeping the Case unchanged, and turning to the step 8.2; otherwise, judging whether the signal frame Count is greater than the minimum signal frame length, if the signal frame Count is less than the minimum signal frame length, determining that the number of the detected signal section frames is too short, and turning to the step 8.1 if no signal exists; otherwise, turning to step 8.8;
and 8.8, judging the signal section, recording the position of the ending frame, adding 1 to the signal section mark, preparing to detect the next signal section, and turning to the step 8.1.
The invention has the beneficial effects that: the invention discloses a digital signal existence detection method based on signal energy-entropy ratio. Since the CME iterative algorithm proposed by the background art requires multiple iterations to converge, it will bring about extremely high computational complexity. Therefore, the invention provides the method for judging whether a signal exists in one frame by using the energy-entropy ratio as the measurement index of the existence of the signal and comparing the energy-entropy ratio with the threshold value through setting the threshold value, so that the iteration is not needed, the complexity of calculation is simplified, and the accuracy of the detection of the existence of the signal is higher under the condition of higher signal-to-noise ratio.
Drawings
FIG. 1 is a general flow of data processing;
FIG. 2 is a flow chart of signal segment dual threshold endpoint detection in accordance with the present invention;
FIG. 3 is a flow chart of signal segment single threshold endpoint detection according to the present invention;
FIG. 4 is a diagram of the effect of simulation detection on a signal segment containing noisy data;
Detailed Description
The technical scheme of the invention is detailed below by combining the accompanying drawings and the embodiment. It should be understood that the scope of the present invention is not limited to the following examples, and any techniques implemented based on the present disclosure are within the scope of the present invention.
According to the program, the following parameters are initially set:
the sampling frequency Fs is 9600Hz, and the frequency point of the known signal is f1800Hz and f21400Hz, a frame data length is 30ms, i.e. the frame length wlen is 288 sampling points, the frame-to-frame overlap rate is set to 50%, the frame shift inc is 144 sampling points, and the threshold setting parameter α is set120.05 and 0.15 respectively, the maximum noise section length Maxsilence is set to 5 frames, and the minimum signal section length MinLen is set to 5 frames.
As shown in fig. 1, for a received digital sequence r (N), N is 1,2, N, the signal presence detecting step is as follows:
step 1, dividing r (n) into frames. The specific framing method and parameter setting are as follows: take 30ms of data as a frame, i.e. the frame is 288 samples long. The overlap ratio between the frame and the frame is set to 50%, i.e. the frame shift (the number of samples of the next frame relative to the current frame) is 144 samples. Total number of minutes of N long data
Figure GDA0002806493260000051
So that the received sequence of length N, r (N), becomes [ r ] after framing1(n),r2(n),...,ri(n),...,rfn(n)];
Step 2, FFT is carried out on each frame of data, and a modulus value is obtained to obtain an energy spectrum Yi(k),i=1,2,...,fn;k=1,2,...,288;
Step 3, calculating the energy of each frame of data: here, since the modulation frequency point is known, the energy is obtained from the frequency domain:
Figure GDA0002806493260000052
i=1,2,...,fn,
Figure GDA0002806493260000053
step 4, calculating the spectrum entropy value of each frame data: the normalized spectral probability density function for each frequency component is defined as:
Figure GDA0002806493260000054
1,2, fn; the 288 spectral entropy is defined as:
Figure GDA0002806493260000055
1,2., fn, the normalized spectral probability density function distribution is relatively uniform for noise, so the spectral entropy value is relatively large; for digital signals, because the frequency spectrum has formant frequency spectrum characteristics and the normalized spectral probability density function of the frequency spectrum is unevenly distributed, the spectral entropy of the digital signals is generally lower than that of noise;
step 5, calculating the energy-entropy ratio of each frame of data:
Figure GDA0002806493260000056
i=1,2,...,fn;
step 6, setting double thresholds: and calculating the maximum value Me of the energy-entropy ratios of all the frames and an initial mean value eth, wherein eth is the mean value of the energy-entropy ratios of the previous 15 frames, and considering that the signal segment has noise of at least 15 frames before the beginning. Setting the threshold Th1 as a noise threshold: th1 ═ 0.05 × (Me-eth) + eth, with a threshold Th2 set as the signal threshold: th2 ═ 0.15 × (Me-eth) + eth;
step 7, setting the maximum allowable noise frame length MaxSilence and the minimum signal frame length MinLen, and starting the end point detection of the signal section;
as shown in fig. 2, the end point detection algorithm for the signal segment has the following steps:
step 1, initialization: the signal frame Count is 0, the status flag Case is 0, and the noise frame Count is 0;
step 2, reading the energy-entropy ratio Eer of one frame of data, and turning to step 3 if Case is 0 or 1; if Case is 2, turning to step 6;
and 3, comparing the energy-entropy ratio with the signal threshold Th 2. If the current frame number is larger than Th2, the signal is detected, Count is increased by 1, the current frame number is recorded as the initial frame position of the signal, Case is set to be 2, and step 2 is switched. Otherwise, turning to the step 4;
and 4, comparing the energy-entropy ratio Eer with the noise threshold Th1, wherein the energy-entropy ratio Eer is larger than Th1, indicating that a signal possibly exists, increasing the Count by 1, setting the Case to be 1, and turning to the step 2. Otherwise, turning to the step 5;
step 5, turning to step 1 when no signal is detected;
and 6, comparing the energy-entropy ratio Eer with the noise threshold Th1, wherein the energy-entropy ratio is larger than Th1, indicating that the frame data is still in the signal section, the Count is increased by 1, the Case is unchanged, and turning to the step 2. Otherwise, turning to step 7;
step 7, noise exists, and the Silence is increased by 1 (possibly due to signal fading, the energy entropy of a signal segment is lower); and judging whether the Silence reaches the maximum allowable noise frame length, if not, continuing to serve as a part of the signal section by the current frame, increasing the Count by 1, keeping the Case unchanged, and turning to the step 2. Otherwise, judging whether the signal frame Count is greater than the minimum signal frame length, if so, judging that the number of the detected signal section frames is too short, considering that no signal exists, and turning to the step 1. Otherwise, turning to the step 8;
and 8, judging the signal segment, recording the position of the ending frame, adding 1 to the signal segment mark, preparing to detect the next signal segment, and turning to the step 1.
End point detection is carried out on a section of noisy data by Matlab simulation software, the detection effect is shown in figure 4 when the signal-to-noise ratio is 15dB, and a detected signal section is arranged between two sections of marked lines, so that the position of the signal section can be accurately found.

Claims (1)

1. A digital signal existence detection method based on signal energy-entropy ratio comprises the following steps:
step 1, sampling a received signal r (t) by a sampling frequency Fs to obtain r (N), setting the data length of r (N) as N, and framing r (N); the specific framing method and parameter setting are as follows: taking t seconds of data as a frame, namely the frame length wlen is t multiplied by Fs; the overlap ratio between the frame and the interframe is set as n%, that is, the frame shift is inc ═ wlen-wlen × n%; total number of minutes of N long data
Figure FDA0002806493250000011
So that the received sequence of length N, r (N), becomes [ r ] after framing1(n),r2(n),...,ri(n),...,rfn(n)]Wherein each element represents a frame of data;
step 2, performing fast Fourier transform on each frame data, and obtaining a modulus value to obtain an energy spectrum Yi(k),i=1,2,...,fn;k=1,2,...,wlen;
And 3, calculating the energy of each frame of data, wherein the time domain energy calculating method comprises the following steps: el (electro luminescence)i=log10(1+AMPiA), i ═ 1,2,. fn, where
Figure FDA0002806493250000012
a is a set constant;
step 4, calculating the spectrum entropy value of each frame data: the normalized spectral probability density function for each frequency component within a frame of data is:
Figure FDA0002806493250000013
the spectral entropy is:
Figure FDA0002806493250000014
step 5, calculatingEnergy-entropy ratio of each frame data:
Figure FDA0002806493250000015
step 6, setting double thresholds: the setting of the threshold is related to the energy-entropy ratio of each frame, and the energy-entropy ratio is obtained through calculation; the method comprises the following specific steps: solving the maximum value Me and the initial average value eth of the energy-entropy ratios of all the frames, wherein eth is the average value of the energy-entropy ratios of the previous NIS frames, and NIS is a constant and represents the length of the noise frame before the signal segment starts; setting the threshold Th1 as a noise threshold: th1 ═ alpha1X (Me-eth) + eth, setting threshold Th2 as the signal threshold: th2 ═ alpha2×(Me-eth)+eth,α12Is constant and alpha1<α2
Step 7, setting the maximum allowable noise frame length MaxSilence and the minimum signal frame length MinLen, and starting the end point detection of the signal section;
and 8: detecting an end point of a signal segment;
step 8.1, initialization: the signal frame Count is 0, the status flag Case is 0, and the noise frame Count is 0;
step 8.2, reading the energy-entropy ratio of one frame of data, and if Case is 0 or 1, turning to step 8.3; if Case is 2, go to step 8.6
Step 8.3, comparing the energy-entropy ratio with the signal threshold Th 2; if the current frame number is greater than Th2, the signal is detected, the Count is increased by 1, the current frame number is recorded as the initial frame position of the signal, Case is set to be 2, and the step is switched to step 8.2; otherwise, turning to step 8.4;
8.4, comparing the energy entropy ratio with the noise threshold Th1, if the energy entropy ratio is larger than Th1, indicating that a signal possibly exists, increasing the Count by 1, setting the Case to be 1, and turning to the step 8.2; otherwise, turning to step 8.5;
step 8.5, switching to step 8.1 when no signal is detected;
8.6, comparing the magnitude of the energy entropy ratio and the noise threshold Th1, wherein the magnitude is larger than Th1, and the data of the frame is still in the signal section, the Count is increased by 1, the Case is unchanged, and the step 8.2 is switched; otherwise, turning to step 8.7;
step 8.7, noise exists, and the Silence is increased by 1; judging whether the Silence reaches the maximum allowable noise frame length, if not, continuing to serve the current frame as a part of the signal section, increasing the Count by 1, keeping the Case unchanged, and turning to the step 8.2; otherwise, judging whether the signal frame Count is greater than the minimum signal frame length, if the signal frame Count is less than the minimum signal frame length, determining that the number of the detected signal section frames is too short, and turning to the step 8.1 if no signal exists; otherwise, turning to step 8.8;
and 8.8, judging the signal section, recording the position of the ending frame, adding 1 to the signal section mark, preparing to detect the next signal section, and turning to the step 8.1.
CN201811356297.XA 2018-11-15 2018-11-15 Digital signal existence detection method based on signal energy-entropy ratio Active CN109412763B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811356297.XA CN109412763B (en) 2018-11-15 2018-11-15 Digital signal existence detection method based on signal energy-entropy ratio

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811356297.XA CN109412763B (en) 2018-11-15 2018-11-15 Digital signal existence detection method based on signal energy-entropy ratio

Publications (2)

Publication Number Publication Date
CN109412763A CN109412763A (en) 2019-03-01
CN109412763B true CN109412763B (en) 2021-03-30

Family

ID=65473269

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811356297.XA Active CN109412763B (en) 2018-11-15 2018-11-15 Digital signal existence detection method based on signal energy-entropy ratio

Country Status (1)

Country Link
CN (1) CN109412763B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110233687B (en) * 2019-04-30 2021-04-30 电子科技大学 Multi-channel data signal joint decision detection method
CN110336764B (en) * 2019-07-04 2021-03-30 电子科技大学 Short wave channel blind symbol synchronization method based on diversity signal decoding feedback
CN110290084B (en) * 2019-07-04 2021-03-30 电子科技大学 Short wave channel blind symbol synchronization method based on data frequency energy peak value
CN110706693B (en) * 2019-10-18 2022-04-19 浙江大华技术股份有限公司 Method and device for determining voice endpoint, storage medium and electronic device
CN112201269B (en) * 2020-10-19 2021-09-07 深圳市车宝信息科技有限公司 MMSE-LSA speech enhancement method based on improved noise estimation
CN112102851B (en) * 2020-11-17 2021-04-13 深圳壹账通智能科技有限公司 Voice endpoint detection method, device, equipment and computer readable storage medium
CN113838476B (en) * 2021-09-24 2023-12-01 世邦通信股份有限公司 Noise estimation method and device for noisy speech
CN117250603B (en) * 2023-11-10 2024-01-09 西北工业大学深圳研究院 Multichannel entropy detection method for underwater weak target signal

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1426048A (en) * 2001-12-13 2003-06-25 中国科学院自动化研究所 End detection method based on entropy
CN104125579A (en) * 2014-08-07 2014-10-29 桂林电子科技大学 Spectrum sensing method and device based on time-domain energy and frequency-domain spectral entropy
CN106653062A (en) * 2017-02-17 2017-05-10 重庆邮电大学 Spectrum-entropy improvement based speech endpoint detection method in low signal-to-noise ratio environment
GB2554943A (en) * 2016-10-16 2018-04-18 Sentimoto Ltd Voice activity detection method and apparatus

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103117067B (en) * 2013-01-19 2015-07-15 渤海大学 Voice endpoint detection method under low signal-to-noise ratio
CN104934032B (en) * 2014-03-17 2019-04-05 华为技术有限公司 The method and apparatus that voice signal is handled according to frequency domain energy
CN107833582B (en) * 2017-11-20 2021-02-09 南京财经大学 Arc length-based voice signal endpoint detection method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1426048A (en) * 2001-12-13 2003-06-25 中国科学院自动化研究所 End detection method based on entropy
CN104125579A (en) * 2014-08-07 2014-10-29 桂林电子科技大学 Spectrum sensing method and device based on time-domain energy and frequency-domain spectral entropy
GB2554943A (en) * 2016-10-16 2018-04-18 Sentimoto Ltd Voice activity detection method and apparatus
CN106653062A (en) * 2017-02-17 2017-05-10 重庆邮电大学 Spectrum-entropy improvement based speech endpoint detection method in low signal-to-noise ratio environment

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于子带能熵比的语音端点检测算法;张毅等;《计算机科学》;20170531;第304-307页 *

Also Published As

Publication number Publication date
CN109412763A (en) 2019-03-01

Similar Documents

Publication Publication Date Title
CN109412763B (en) Digital signal existence detection method based on signal energy-entropy ratio
US8699635B2 (en) Frame boundary detection
JP6377862B2 (en) Encoder selection
US20180262832A1 (en) Sound Signal Processing Apparatus and Method for Enhancing a Sound Signal
EP3252771B1 (en) A method and an apparatus for performing a voice activity detection
CN113132287B (en) Synchronous detection method and device in OFDM system
US20140337039A1 (en) Frame Loss Compensation Method And Apparatus For Voice Frame Signal
CN106653062A (en) Spectrum-entropy improvement based speech endpoint detection method in low signal-to-noise ratio environment
JP4384814B2 (en) How to select a synchronization pattern from a test signal
CN107835036B (en) Non-cooperative frequency hopping signal cracking method
CN103632681B (en) A kind of spectral envelope silence detection method
CN114268528B (en) Frame capturing synchronization method and device
CN109061591B (en) Time-frequency line spectrum detection method based on sequential clustering
CN113447893B (en) Radar pulse signal frequency spectrum automatic detection method, system and medium
CN103618575A (en) Automatic real-time frame synchronization method of frequency-modulated hydro-acoustic communication system
CN109087657B (en) Voice enhancement method applied to ultra-short wave radio station
JP2005516247A (en) Voice activity detector and enabler for noisy environments
CN101611606A (en) The symbol timing synchronization method of the digital signal that receives and the digit signal receiver that utilizes this method
CN110233687B (en) Multi-channel data signal joint decision detection method
CN109507675B (en) Method for realizing underwater multi-target time delay estimation processing based on frequency division system
CN107612864A (en) Burst mode simplifies the decision threshold value setting method of symbol timing synchronization method
JP5944565B2 (en) Method and apparatus for initial gain control in a communication system
CN105703850B (en) Based on short-time Fourier transform according to chain signal edge detection method
JP2008085894A (en) Symbol speed detection device, and program
CN114944975B (en) Signal modulation type real-time identification method based on FPGA processing

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant