CN112151065A - Method, device and equipment for detecting single tone signal frequency and computer storage medium - Google Patents
Method, device and equipment for detecting single tone signal frequency and computer storage medium Download PDFInfo
- Publication number
- CN112151065A CN112151065A CN201910578565.0A CN201910578565A CN112151065A CN 112151065 A CN112151065 A CN 112151065A CN 201910578565 A CN201910578565 A CN 201910578565A CN 112151065 A CN112151065 A CN 112151065A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signals
- sampling points
- signal
- tone
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000005070 sampling Methods 0.000 claims abstract description 87
- 238000001514 detection method Methods 0.000 claims abstract description 52
- 238000009825 accumulation Methods 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 238000006243 chemical reaction Methods 0.000 claims description 44
- 238000001914 filtration Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 11
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 2
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 238000012958 reprocessing Methods 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 19
- 238000004891 communication Methods 0.000 description 10
- 238000003775 Density Functional Theory Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000011664 signaling Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000012850 discrimination method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/48—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Human Computer Interaction (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Computational Mathematics (AREA)
- Signal Processing (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
The embodiment of the invention discloses a method, a device and equipment for detecting single tone signal frequency and a computer storage medium, wherein the method comprises the following steps: obtaining the single tone signals of m sampling points, carrying out frequency estimation on the single tone signals of the m sampling points to obtain m frequency estimation values, and carrying out at least one accumulation average calculation on the m frequency estimation values to determine the frequency value of the single tone signals. By adopting the detection method, the frequency blind detection can be realized without acquiring the frequency of the received single tone signal in advance, and the signal-to-noise ratio and the detection precision are improved by reprocessing the frequency estimation value obtained after the frequency estimation processing, so that the detection error reaches the detection precision smaller than 2 Hz.
Description
Technical Field
The invention belongs to the field of signal detection, and particularly relates to a method, a device and equipment for detecting the frequency of a single tone signal and a computer storage medium.
Background
In a wireless communication system, a receiving device needs to perform frequency detection on a single-frequency real signal and identify signaling, address numbers or characters according to a detection result, for example, analog tone transmission signaling, dial numbers or characters are often used in analog frequency modulation-based trunking communication, analog talkback and analog wired telephone transmission.
Currently, common methods for detecting tone (frequency) signals are the gratzel algorithm (Goertzel algorithm) and the discrete fourier transform DFT algorithm. However, in the conventional method, it is necessary to detect a signal when a frequency point is known and confirm whether a received signal is a known frequency signal, so that frequency blind detection cannot be realized, and the accuracy of the detected result is not high.
Disclosure of Invention
The embodiment of the invention provides a method, a device and equipment for detecting the frequency of a single tone signal and a computer storage medium, which can realize the blind detection of the single tone signal and improve the effect of detecting the frequency precision.
In one aspect, an embodiment of the present invention provides a method for detecting a single-tone signal frequency, where the method includes:
acquiring single tone signals of m sampling points;
carrying out frequency estimation on the single tone signals of the m sampling points to obtain m frequency estimation values;
and performing at least one accumulation average calculation on the m frequency estimation values to determine the frequency value of the single tone signal.
In another aspect, an embodiment of the present invention provides a single tone signal frequency detection apparatus, including:
the single tone signal acquisition module is used for acquiring single tone signals of m sampling points;
the frequency discrimination module is used for carrying out frequency estimation on the single-tone signals of the m sampling points to obtain m frequency estimation values;
and the tone signal frequency determining module is used for performing at least one accumulated average calculation on the m frequency estimated values to determine the frequency value of the tone signal.
In another aspect, an embodiment of the present invention provides a single tone signal frequency detection apparatus, where the apparatus includes:
a processor and a memory storing computer program instructions;
the processor, when executing the computer program instructions, implements the method as described in any of the above.
In yet another aspect, an embodiment of the present invention provides a computer storage medium having computer program instructions stored thereon, where the computer program instructions, when executed by a processor, implement the method as described in any one of the above.
The method, the device, the equipment and the computer storage medium for detecting the frequency of the single tone signal acquire the single tone signals of m sampling points, perform frequency estimation on the single tone signals of the m sampling points to obtain m frequency estimation values, perform at least one accumulation average calculation on the m frequency estimation values, and determine the frequency value of the single tone signal. By adopting the detection method, the frequency blind detection can be realized without acquiring the frequency of the received single tone signal in advance, and the frequency estimation values obtained after the frequency estimation processing are sequentially subjected to accumulation, elimination and accumulation processing, so that the signal-to-noise ratio and the detection precision are improved, and the detection error reaches the detection precision smaller than 2 Hz.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required to be used in the embodiments of the present invention will be briefly described below, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of a prior art Getz algorithm;
FIG. 2 is a flow chart of a tone signal frequency detection method according to an embodiment of the present invention;
FIG. 3 is a schematic structural diagram of a single tone signal frequency detection method according to an embodiment of the present invention;
FIG. 4 is a schematic structural diagram of a single tone signal frequency detection method according to another embodiment of the present invention;
FIG. 5 is a schematic structural diagram of a single tone frequency detection apparatus according to an embodiment of the present invention;
fig. 6 is a schematic structural diagram of a single tone signal frequency detection apparatus according to another embodiment of the present invention.
Detailed Description
Features and exemplary embodiments of various aspects of the present invention will be described in detail below, and in order to make objects, technical solutions and advantages of the present invention more apparent, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not to be construed as limiting the invention. It will be apparent to one skilled in the art that the present invention may be practiced without some of these specific details. The following description of the embodiments is merely intended to provide a better understanding of the present invention by illustrating examples of the present invention.
It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
At present, the commonly used methods for detecting the single frequency signal include a gratzer algorithm (Goertzel algorithm) and a discrete fourier transform DFT algorithm.
1. Discrete Fourier transform algorithm (DFT)
Setting the time domain samples of the received signal asTaking I +1 time domain samples x (I), I ═ 0,1, …, and I as samples, the frequency domain response of the target frequency point can be calculated according to the following formula:
in the above formula, f is the target frequency point, fsIs the sampling frequency, and I +1 is the number of time domain sampling points used for calculating the frequency domain response of the target frequency point. Discrete Fourier transform algorithm needs to calculate the rotation factor at each stepIs finally according to XI(f) The amplitude can be detected whether the received signal x (i) contains the target frequency point signal.
2. Grazier algorithm
The principle of the lattice algorithm is consistent with the discrete fourier transform, but the implementation is different from the discrete fourier transform. The lattice algorithm implements the discrete fourier transform with the structure of an IIR filter, as shown in fig. 1.
The algorithm represented in fig. 1 can be represented by the following equation:
according to the above formula wherein yI(f) The amplitude can be detected whether the received signal x (i) contains the target frequency point signal.
The common feature of the discrete fourier transform and the lattice algorithm is that the signal needs to be detected under the condition of a known frequency point to determine whether the received signal is a known frequency signal. When the received signal is one of several possible frequencies, several sets of detection modules are required to operate simultaneously, and the frequency of the received signal is determined according to the maximum value of the amplitude of the calculation result. Or both algorithms cannot blindly detect the signal frequency without receiving the signal frequency information.
In order to solve the problems in the prior art, embodiments of the present invention provide a method, an apparatus, a device, and a computer storage medium for single tone frequency detection. The method for detecting the frequency of the monophonic signal according to the embodiment of the present invention will be described first.
The scheme framework of the invention comprises three parts: the frequency discrimination module + M-point accumulated average + N-point comparison + elimination deviation value + N-point accumulated average.
The first frequency discrimination mode:
the frequency discrimination module 1: digital frequency converter, low-pass filter and frequency discriminator.
The second frequency discrimination mode:
and the frequency discrimination module 2: difference, low-pass filtering, received signal amplitude estimation, amplitude normalization, table lookup and frequency calculation.
The method can be used in scenes needing to carry out frequency detection on the analog real single-frequency signal, such as dial tone in a wired telephone or analog talkback equipment, analog sub-tone signaling detection and the like.
Fig. 2 is a flowchart illustrating a method for detecting the frequency of a single tone signal according to an embodiment of the present invention. As shown in fig. 2, the method includes:
step S1: and acquiring single tone signals of m sampling points.
Step S2: and carrying out frequency estimation on the single tone signals of the m sampling points to obtain m frequency estimation values.
Step S3: and performing at least one accumulation average calculation on the m frequency estimation values to determine the frequency value of the single tone signal.
Specifically, in steps S1-S3, the monophonic signal is also called a single-frequency signal, and the monophonic signals of m (m is a positive integer greater than 1) sampling points are collected, wherein the frequency of each sampling point is unknown. This application adopts two kinds of modes of discriminating right the monophonic signal of m sampling points carries out frequency estimation to obtain m frequency estimation value.
The method, the device, the equipment and the computer storage medium for detecting the frequency of the single tone signal acquire the single tone signals of m sampling points, perform frequency estimation on the single tone signals of the m sampling points to obtain m frequency estimation values, perform at least one accumulation calculation on the m frequency estimation values, and determine the frequency value of the single tone signal. According to the detection method, blind frequency detection can be achieved without acquiring the frequency of the received single tone signal in advance. The frequency estimation values obtained after the frequency estimation processing are sequentially subjected to accumulation, rejection and accumulation processing, so that the signal-to-noise ratio and the detection precision can be improved, and the detection error can reach the detection precision of less than 2 Hz.
Referring to fig. 3, a first frequency discrimination method is as follows:
in one embodiment, step S2 includes:
step S21: and carrying out digital frequency conversion on the single-tone signals of the m sampling points to obtain difference frequency signals and sum frequency signals of the m sampling points after the single-tone signals are subjected to frequency conversion.
Step S22: and filtering the sum frequency signals obtained after the frequency conversion of the single-tone signals of the m sampling points to obtain the difference frequency signals obtained after the frequency conversion of the single-tone signals of the m sampling points, and performing frequency discrimination processing on the difference frequency signals obtained after the frequency conversion of the single-tone signals of the m sampling points to obtain m frequency estimation values.
the preset frequency refers to a frequency value set inside the system. In order to obtain the desired difference frequency signal, signals of other frequency bands (i.e. sum frequency signals) are filtered out.
In one embodiment, step S21 includes:
step S211: and acquiring an I path signal corresponding to the m sampling points after the frequency conversion of the single tone signals and a Q path signal corresponding to the m sampling points after the frequency conversion of the single tone signals.
Step S212: and carrying out digital frequency conversion on the single-tone signals of the m sampling points to respectively obtain a difference frequency signal and a sum frequency signal of an I-path signal corresponding to the single-tone signals of the m sampling points after frequency conversion, and a difference frequency signal and a sum frequency signal of a Q-path signal corresponding to the single-tone signals of the m sampling points after frequency conversion.
In steps S211-S212, the monophonic signal of m sample points is represented as: the single-tone signal after frequency conversion comprises an I path signal and a Q path signal, and two algorithms are introduced here, wherein the first algorithm is specifically as follows:
the second algorithm:
in the formulas (3) to (6), f is the received signal frequency, foIs the central frequency, f, adopted during digital frequency conversionsIs the sampling frequency of the received signal,For a white noise value at the sampling point,is an in-phase white noise component,is an orthogonal white noise component.
In one embodiment, step S22 includes:
step S221: respectively filtering the sum frequency signal of the I-path signal corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion, and the sum frequency signal of the Q-path signal corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion, and determining the difference frequency signal of the I-path signal corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion and the difference frequency signal of the Q-path signal corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion.
Step S222: and carrying out frequency discrimination processing on the difference frequency signal of the I path signal corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion and the difference frequency signal of the Q path signal corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion to obtain m frequency values.
In steps S221 to S222, the sum frequency component, the in-phase white noise component and the quadrature white noise component of the I/Q signal are filtered out by the low pass filter from the signal after the digital frequency conversion, and the difference frequency signal of the I/Q signal is input into the frequency discriminator for frequency discrimination to determine m frequency values.
In one embodiment, step S222 includes:
step S2221: and acquiring a central frequency value adopted during digital frequency conversion.
Step S2222: and calculating each frequency value in the m frequency values and the central frequency value to obtain m frequency estimation values.
In steps S2221 to S2222, for equations (3) and (4), the output after frequency discrimination is: for equations (5) and (6), the output after frequency discrimination is:where f isdIs the frequency of the output of the frequency discriminator,Is the single tone signal frequency estimated by the present algorithm.
Referring to fig. 4, the second frequency discrimination method is as follows:
in one embodiment, step S2 further includes:
step S23: and carrying out differential processing on the single-tone signals of the m sampling points to obtain difference frequency signals and sum frequency signals of the single-tone signals of the m sampling points.
Step S24: and filtering the sum frequency signal in the single-tone signals of the m sampling points, and determining the difference frequency signal of the single-tone signals of the m sampling points.
Step S25: and calculating the square of the single tone signals of the m sampling points to obtain the square value of the single tone signals of the m sampling points.
Step S26: and carrying out low-pass filtering on the square values of the single tone signals of the m sampling points to obtain the amplitude estimation values of the single tone signals of the m sampling points.
Step S27: and sequentially carrying out normalization and inverse cosine operation on the difference frequency signals of the single tone signals of the m sampling points by using the amplitude estimation values of the single tone signals of the m sampling points to obtain m frequency estimation values.
In steps S23-S27, the monophonic signals of m sampling points are: the first path carries out differential processing on the single-tone signals of the m sampling points, and the output signals of the differentiator are as follows:
assuming a single tone signal as a constant envelope signal, Ai=Ai-1Then the above equation can be further expanded as:
further low pass filtering equation (8) yields the signal:
the second path performs square operation on the single tone signals of the m sampling points to obtain the following formula:
low-pass filtering the formula (10) to obtain the amplitude estimation value of the received signalUsing the amplitude estimate to pair the differentiated output signalsCarrying out amplitude normalization to obtain an output signalAnd further looking up a table to realize an inverse cosine operation, and obtaining frequency estimation values of the m single tone signals as follows:
in one embodiment, step S3 includes:
step S31: and sequentially selecting a first preset number of frequency estimation values from the m frequency estimation values, and performing accumulation averaging on the first preset number of frequency estimation values to obtain n first accumulation average results.
Step S32: and eliminating abnormal first accumulated average results in the n first accumulated average results to obtain k first accumulated average results.
Step S33: sequentially selecting a second preset number of first accumulated average results from the k first accumulated average results, performing accumulated averaging on the second preset number of first accumulated average results to obtain at least one second accumulated average result, and taking the at least one second accumulated average result as the frequency value of the single tone signal; wherein m, n and k are integers, m > n > k, and k > 1.
In steps S31-S33, m frequency estimates are obtained after the mono signal is frequency discriminated in the above two ways. Further, m-point accumulation is carried out on the m frequency estimation values, and the signal-to-noise ratio is improved; and for the signals after the m-point accumulation, n points are taken for observation, the values with obviously deviated values are removed from the n-point signals, and then accumulation is carried out, so that the signal-to-noise ratio is further improved. And finally, taking the residual frequency estimation value as the frequency value of the single tone signal. The first preset number and the second preset number refer to fixed numbers set inside the system, the first preset number is represented by M, the second preset number is represented by N, and M, N are integers greater than 1.
In one embodiment, the step S32 includes:
step S321: and if the n first accumulated average results have first accumulated average results exceeding a preset threshold range, rejecting the first accumulated average results exceeding the preset threshold range to obtain k first accumulated average results.
Specifically, the preset range threshold refers to a value range set in advance with respect to the frequency estimation value. The obtaining of the k first accumulation results is not limited to the above embodiment, and a preset number of first accumulation results may be selected first, then an average value of the preset number of first accumulation results is obtained, the obtained average value is compared with the preset average value, and if there is a first accumulation result significantly deviating from the preset average value, the first accumulation result significantly deviating from the average value is removed to determine the k first accumulation results.
The advantages of the invention include:
1. the frequency detection is carried out by adopting a digital frequency conversion and frequency discrimination method, so that the blind detection of a single-frequency real number signal can be realized, the frequency of a received signal needs to be known in the prior art, and the blind detection of the frequency cannot be realized.
2. The method of accumulating M points and accumulating N points after eliminating deviation values makes full use of the characteristic of the frequency discriminated signal, ensures the detection performance under the condition of low signal-to-noise ratio and is flexible to realize.
In the prior art, the signal-to-noise ratio of the detection result can be improved by increasing the DFT length or the filtering time length, and the frequency detection accuracy is improved. However, the cost of increasing the DFT length or increasing the filtering time is to increase the detection time, so that the detection time cannot meet the requirement of the communication system in the case where the detection time is sensitive, for example, in the case of using a single-tone signal as the communication signaling. Especially, when the known frequency information of the receiving end is inconsistent with the received actual frequency information, the DFT calculation in the first technology needs to be restarted, or the intermediate value of the filter in the second technology needs to be cleared, so that the purpose of detecting a new frequency point is achieved, and the response to the mutation situation is not flexible. The technology adopted by the invention can quickly obtain a preliminary detection result by observing the output value after M-point accumulation, and can further improve the signal-to-noise ratio and improve the detection precision by a method of removing the deviation value by N-point comparison and then accumulating. In addition, the change of the frequency of the received signal can be conveniently found through observation, and the emergency response is flexible.
Referring to fig. 5, the present application also provides a tone signal frequency detection apparatus, including:
and the signal acquisition module 10 is used for acquiring the single tone signals of the m sampling points.
And the frequency discrimination module 20 is configured to perform frequency estimation on the single-tone signals of the m sampling points to obtain m frequency estimation values.
And the signal frequency determining module 30 is configured to perform at least one cumulative average calculation on the m frequency estimation values to determine the frequency value of the monophonic signal.
Fig. 6 is a schematic diagram illustrating a hardware structure of tone information frequency detection according to an embodiment of the present invention.
The tone information frequency detection device may include a processor 301 and a memory 302 having stored computer program instructions.
In particular, the processor 301 may include a Central Processing Unit (CPU), or an Application Specific Integrated Circuit (ASIC), or may be configured as one or more Integrated circuits implementing embodiments of the present invention.
The processor 301 may implement any of the above-described tone signal frequency detection methods in the embodiments by reading and executing computer program instructions stored in the memory 302.
In one example, the tone information frequency detection device may also include a communication interface 303 and a bus 310. As shown in fig. 3, the processor 301, the memory 302, and the communication interface 303 are connected via a bus 310 to complete communication therebetween.
The communication interface 303 is mainly used for implementing communication between modules, apparatuses, units and/or devices in the embodiment of the present invention.
In addition, in combination with the method for detecting the frequency of the monophonic signal in the foregoing embodiments, the embodiments of the present invention may be implemented by providing a computer storage medium. The computer storage medium having computer program instructions stored thereon; the computer program instructions, when executed by a processor, implement any one of the tone signal frequency detection methods in the above embodiments.
It is to be understood that the invention is not limited to the specific arrangements and instrumentality described above and shown in the drawings. A detailed description of known methods is omitted herein for the sake of brevity. In the above embodiments, several specific steps are described and shown as examples. However, the method processes of the present invention are not limited to the specific steps described and illustrated, and those skilled in the art can make various changes, modifications and additions or change the order between the steps after comprehending the spirit of the present invention.
The functional blocks shown in the above-described structural block diagrams may be implemented as hardware, software, firmware, or a combination thereof. When implemented in hardware, it may be, for example, an electronic circuit, an Application Specific Integrated Circuit (ASIC), suitable firmware, plug-in, function card, or the like. When implemented in software, the elements of the invention are the programs or code segments used to perform the required tasks. The program or code segments may be stored in a machine-readable medium or transmitted by a data signal carried in a carrier wave over a transmission medium or a communication link. A "machine-readable medium" may include any medium that can store or transfer information. Examples of a machine-readable medium include electronic circuits, semiconductor memory devices, ROM, flash memory, Erasable ROM (EROM), floppy disks, CD-ROMs, optical disks, hard disks, fiber optic media, Radio Frequency (RF) links, and so forth. The code segments may be downloaded via computer networks such as the internet, intranet, etc.
It should also be noted that the exemplary embodiments mentioned in this patent describe some methods or systems based on a series of steps or devices. However, the present invention is not limited to the order of the above-described steps, that is, the steps may be performed in the order mentioned in the embodiments, may be performed in an order different from the order in the embodiments, or may be performed simultaneously.
As described above, only the specific embodiments of the present invention are provided, and it can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working processes of the system, the module and the unit described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again. It should be understood that the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive various equivalent modifications or substitutions within the technical scope of the present invention, and these modifications or substitutions should be covered within the scope of the present invention.
Claims (11)
1. A method for single tone signal frequency detection, the method comprising:
acquiring single tone signals of m sampling points;
carrying out frequency estimation on the single tone signals of the m sampling points to obtain m frequency estimation values;
and performing at least one accumulation average calculation on the m frequency estimation values to determine the frequency value of the single tone signal.
2. The method of claim 1, wherein the frequency estimating the monophonic signals of the m sampling points to obtain m frequency estimated values comprises:
carrying out digital frequency conversion on the single-tone signals of the m sampling points to obtain difference frequency signals and sum frequency signals after frequency conversion of the single-tone signals of the m sampling points;
and filtering the sum frequency signals obtained after the frequency conversion of the single-tone signals of the m sampling points to obtain the difference frequency signals obtained after the frequency conversion of the single-tone signals of the m sampling points, and performing frequency discrimination processing on the difference frequency signals obtained after the frequency conversion of the single-tone signals of the m sampling points to obtain m frequency estimation values.
3. The method of claim 2, wherein the step of digitally converting the mono-tone signals of the m sampling points to obtain the difference frequency signals and the sum frequency signals after the frequency conversion of the mono-tone signals of the m sampling points comprises:
acquiring an I path signal corresponding to the m sampling points after the frequency conversion of the single tone signals and a Q path signal corresponding to the m sampling points after the frequency conversion of the single tone signals;
and carrying out digital frequency conversion on the single-tone signals of the m sampling points to respectively obtain a difference frequency signal and a sum frequency signal of an I-path signal corresponding to the single-tone signals of the m sampling points after frequency conversion, and a difference frequency signal and a sum frequency signal of a Q-path signal corresponding to the single-tone signals of the m sampling points after frequency conversion.
4. The method of claim 3, wherein the filtering out the sum frequency signals obtained by frequency conversion of the mono-tone signals of the m sampling points, determining the difference frequency signals obtained by frequency conversion of the mono-tone signals of the m sampling points, and performing frequency discrimination on the difference frequency signals obtained by frequency conversion of the mono-tone signals of the m sampling points to obtain m frequency estimation values comprises:
respectively filtering out sum frequency signals of the I-path signals corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion, and sum frequency signals of the Q-path signals corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion, and determining difference frequency signals of the I-path signals corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion and difference frequency signals of the Q-path signals corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion;
and carrying out frequency discrimination processing on the difference frequency signal of the I path signal corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion and the difference frequency signal of the Q path signal corresponding to the m sampling points after the single-tone signals are subjected to frequency conversion to obtain m frequency values.
5. The method of claim 4, wherein the frequency discriminating between the difference frequency signal of the I-path signal corresponding to the m-sampled point mono-tone signal after frequency conversion and the difference frequency signal of the Q-path signal corresponding to the m-sampled point mono-tone signal after frequency conversion to obtain m frequency values comprises:
acquiring a central frequency value adopted during digital frequency conversion;
and calculating each frequency value in the m frequency values and the central frequency value to obtain m frequency estimation values.
6. The method of claim 1, wherein the frequency estimating the monophonic signals of the m sampling points to obtain m frequency estimated values further comprises:
carrying out differential processing on the single-tone signals of the m sampling points to obtain difference frequency signals and sum frequency signals after the single-tone signals of the m sampling points are differentiated;
filtering out the sum frequency signal after the single tone signals of the m sampling points are differentiated, and determining the difference frequency signal after the single tone signals of the m sampling points are differentiated;
calculating the square of the single tone signals of the m sampling points to obtain the square value of the single tone signals of the m sampling points;
carrying out low-pass filtering on the square values of the single tone signals of the m sampling points to obtain amplitude estimation values of the single tone signals of the m sampling points;
and sequentially carrying out normalization and inverse cosine operation on the difference frequency signals after the single tone signals of the m sampling points are differentiated by utilizing the amplitude estimation values of the single tone signals of the m sampling points to obtain m frequency estimation values.
7. The method of claim 5 or 6, wherein said performing at least one cumulative average calculation on said m frequency estimates and determining the frequency values of said monophonic signal comprises:
sequentially selecting a first preset number of frequency estimation values from the m frequency estimation values, and performing accumulation averaging on the first preset number of frequency estimation values to obtain n first accumulation average results;
rejecting abnormal first accumulated average results in the n first accumulated average results to obtain k first accumulated average results;
sequentially selecting a second preset number of first accumulated average results from the k first accumulated average results, performing accumulated averaging on the second preset number of first accumulated average results to obtain at least one second accumulated average result, and taking the at least one second accumulated average result as the frequency value of the single tone signal; wherein m, n and k are integers, m > n > k, and k > 1.
8. The method of claim 7, wherein said culling outlier first accumulated average results from the n first accumulated average results to obtain k first accumulated average results comprises:
and if the n first accumulated average results have first accumulated average results exceeding a preset threshold range, rejecting the first accumulated average results exceeding the preset threshold range to obtain k first accumulated average results.
9. An apparatus for single tone frequency detection, the apparatus comprising:
the single tone signal acquisition module is used for acquiring single tone signals of m sampling points;
the frequency discrimination module is used for carrying out frequency estimation on the single-tone signals of the m sampling points to obtain m frequency estimation values;
and the tone signal frequency determining module is used for performing at least one accumulated average calculation on the m frequency estimated values to determine the frequency value of the tone signal.
10. A tone signal frequency detection apparatus, the apparatus comprising: a processor and a memory storing computer program instructions;
the processor, when executing the computer program instructions, implements the method of any of claims 1-8.
11. A computer storage medium having computer program instructions stored thereon which, when executed by a processor, implement the method of any one of claims 1-8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910578565.0A CN112151065B (en) | 2019-06-28 | 2019-06-28 | Method, device, equipment and computer storage medium for detecting single-tone signal frequency |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910578565.0A CN112151065B (en) | 2019-06-28 | 2019-06-28 | Method, device, equipment and computer storage medium for detecting single-tone signal frequency |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112151065A true CN112151065A (en) | 2020-12-29 |
CN112151065B CN112151065B (en) | 2024-03-15 |
Family
ID=73891491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910578565.0A Active CN112151065B (en) | 2019-06-28 | 2019-06-28 | Method, device, equipment and computer storage medium for detecting single-tone signal frequency |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112151065B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023197933A1 (en) * | 2022-04-11 | 2023-10-19 | 维沃移动通信有限公司 | Signal processing method and apparatus, and electronic device and medium |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1130577A2 (en) * | 2000-03-02 | 2001-09-05 | Volkswagen Aktiengesellschaft | Method for the reconstruction of low speech frequencies from mid-range frequencies |
CN101031814A (en) * | 2004-09-28 | 2007-09-05 | 秦内蒂克有限公司 | Frequency modulated continuous wave (FMCW) radar having improved frequency sweep linearity |
US20080240301A1 (en) * | 2006-12-19 | 2008-10-02 | Soumyajit Mandal | Architectures for universal or software radio |
US20110103573A1 (en) * | 2008-06-30 | 2011-05-05 | Freescale Semiconductor Inc. | Multi-frequency tone detector |
CN202177705U (en) * | 2011-05-03 | 2012-03-28 | 上海通用卫星导航有限公司 | High-precision caesium optical pump signal processor |
CN102435842A (en) * | 2011-09-20 | 2012-05-02 | 武汉理工大学 | Device for comparing double-mode frequency scales and measuring frequency stability |
CN103234624A (en) * | 2012-08-20 | 2013-08-07 | 苏州大学 | High-precision frequency estimation system |
CN104569581A (en) * | 2015-01-30 | 2015-04-29 | 湖南大学 | Multi-level set and single-cycle estimation method of power grid frequency measuring |
US20150115175A1 (en) * | 2013-10-28 | 2015-04-30 | Nathan T. Baltz | Apparatus and method for system identification |
JP5717912B1 (en) * | 2014-09-26 | 2015-05-13 | 雅紀 廣石 | Frequency measuring device, frequency measuring method, and oscillation type sensor |
CN107271772A (en) * | 2017-06-12 | 2017-10-20 | 南京工程学院 | A kind of mains frequency quick determination method of high accuracy and anti-noise jamming |
CN108759809A (en) * | 2018-08-10 | 2018-11-06 | 维沃移动通信有限公司 | A kind of gyroscope detection circuit and terminal |
-
2019
- 2019-06-28 CN CN201910578565.0A patent/CN112151065B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1130577A2 (en) * | 2000-03-02 | 2001-09-05 | Volkswagen Aktiengesellschaft | Method for the reconstruction of low speech frequencies from mid-range frequencies |
CN101031814A (en) * | 2004-09-28 | 2007-09-05 | 秦内蒂克有限公司 | Frequency modulated continuous wave (FMCW) radar having improved frequency sweep linearity |
US20080240301A1 (en) * | 2006-12-19 | 2008-10-02 | Soumyajit Mandal | Architectures for universal or software radio |
US20110103573A1 (en) * | 2008-06-30 | 2011-05-05 | Freescale Semiconductor Inc. | Multi-frequency tone detector |
CN202177705U (en) * | 2011-05-03 | 2012-03-28 | 上海通用卫星导航有限公司 | High-precision caesium optical pump signal processor |
CN102435842A (en) * | 2011-09-20 | 2012-05-02 | 武汉理工大学 | Device for comparing double-mode frequency scales and measuring frequency stability |
CN103234624A (en) * | 2012-08-20 | 2013-08-07 | 苏州大学 | High-precision frequency estimation system |
US20150115175A1 (en) * | 2013-10-28 | 2015-04-30 | Nathan T. Baltz | Apparatus and method for system identification |
JP5717912B1 (en) * | 2014-09-26 | 2015-05-13 | 雅紀 廣石 | Frequency measuring device, frequency measuring method, and oscillation type sensor |
CN104569581A (en) * | 2015-01-30 | 2015-04-29 | 湖南大学 | Multi-level set and single-cycle estimation method of power grid frequency measuring |
CN107271772A (en) * | 2017-06-12 | 2017-10-20 | 南京工程学院 | A kind of mains frequency quick determination method of high accuracy and anti-noise jamming |
CN108759809A (en) * | 2018-08-10 | 2018-11-06 | 维沃移动通信有限公司 | A kind of gyroscope detection circuit and terminal |
Non-Patent Citations (2)
Title |
---|
YAQING TU等: "Method for CMF Signal Processing Based on the Recursive DTFT Algorithm With Negative Frequency Contribution", IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT, pages 2647 - 2648 * |
汪少锋: "基于FPGA带频率自适应跟踪的矢量型数字锁相放大器设计", 中国优秀硕士学位论文全文数据库, pages 19 - 35 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023197933A1 (en) * | 2022-04-11 | 2023-10-19 | 维沃移动通信有限公司 | Signal processing method and apparatus, and electronic device and medium |
Also Published As
Publication number | Publication date |
---|---|
CN112151065B (en) | 2024-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109495410B (en) | High dynamic PCM/FM signal carrier frequency accurate estimation method | |
US8619909B2 (en) | Signal detector using matched filter for training signal detection | |
CN1640082B (en) | Detecting and measuring interference contained within a digital carrier | |
EP3621259B1 (en) | Method and device for fsk/gfsk demodulation | |
CN110089038B (en) | Noise suppression device, noise suppression method, and reception device and reception method using the same | |
CN107231225B (en) | Preamble sequence detection method and device | |
CN112151065A (en) | Method, device and equipment for detecting single tone signal frequency and computer storage medium | |
JP2002537728A (en) | Method for detecting pilot signal | |
JP4543283B2 (en) | Radio identification method and apparatus | |
KR20150035819A (en) | Method for searching for a useful signal in a multiplexing band | |
CN214675207U (en) | Device for recognizing pulse amplitude modulation signal based on level control | |
US20220255829A1 (en) | Received-signal rate detection | |
CN113848391A (en) | Pulse signal detection and extraction method | |
CN114495972A (en) | Signal correction method, device, equipment, storage medium and computer program product | |
CN112865882A (en) | Anti-frequency-deviation DMR interphone signal rapid identification method | |
CN115615467B (en) | Hall signal demodulation method and device | |
JP2000295641A (en) | Tone signal detection method and tone signal detector | |
EP0262860A2 (en) | Digital signal analysis by transformation to vector domain | |
US6449321B1 (en) | Radio receiver | |
JP2007093490A (en) | Am/fm mixed signal processor | |
CN117407818A (en) | Abnormality detection method, abnormality detection device, abnormality detection apparatus, and computer-readable storage medium | |
CN115037320A (en) | Method and device for estimating equivalent baseband direct current offset and transmitter | |
CN113206812A (en) | Device for identifying pulse amplitude modulation signal based on level control and application method | |
CN113378699A (en) | Signal processing method, device, equipment and storage medium | |
KR101912446B1 (en) | Method and apparatus for detecting intruder |
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 |