CN103166887A - Sideband modulating signal classifying method and device - Google Patents
Sideband modulating signal classifying method and device Download PDFInfo
- Publication number
- CN103166887A CN103166887A CN2013100765460A CN201310076546A CN103166887A CN 103166887 A CN103166887 A CN 103166887A CN 2013100765460 A CN2013100765460 A CN 2013100765460A CN 201310076546 A CN201310076546 A CN 201310076546A CN 103166887 A CN103166887 A CN 103166887A
- Authority
- CN
- China
- Prior art keywords
- signal
- correlation function
- spectrum correlation
- frequency spectrum
- module
- 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
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The invention discloses a sideband modulating signal classifying method and device. The method includes that an enhancement type frequency spectrum correlation function is acquired according to a frequency spectrum correlation function and a conjugate frequency spectrum correlation function of a time domain digital signal x[n], a cyclical characteristic chromatogram of the time domain digital signal x[n] is figured according to the enhancement type frequency spectrum correlation function, a test statistic Tp is acquired by the cyclical characteristic chromatogram, the test statistic Tp is compared with a threshold value gamma, if Tp is smaller than or equal to gamma, the signal is a single sideband modulating signal, and otherwise, the signal is a double-sideband modulating signal. By means of the method and the device, the application range of sideband signal classification detection is increased.
Description
Technical field
The present invention relates to the cognitive radio technology field, relate in particular to a kind of method and device of sideband modulated signal classification.
Background technology
The Classification and Identification of single sideband modulated signal and double-sideband modulation signal is an important step in the modulated-analog signal classification, plays the part of important function in detection of electrons and military communication.
Existing sorting technique is all be based upon precognition signal carrier frequency or estimate on the basis of signal carrier frequency, and still, the carrier frequency information of signal is difficult to know in advance or accurately estimates in most cases, has caused very large difficulty for the classification of sideband modulated signal.Therefore, how when unknown carrier frequency information, sideband modulated signal effectively to be categorized into and to be important subject.
Summary of the invention
The purpose of this invention is to provide a kind of method and device of sideband modulated signal classification, solved the problem that to know in advance carrier frequency information in existing sideband modulated signal sorting technique.
The objective of the invention is to be achieved through the following technical solutions:
A kind of method of sideband modulated signal classification, the method comprises:
According to time-domain digital signal x[n] the frequency spectrum correlation function
With the conjugate spectrum correlation function
Obtain enhancement mode frequency spectrum correlation function
According to described enhancement mode frequency spectrum correlation function
Calculate described time-domain digital signal x[n] the cycle specificity collection of illustrative plates, and obtain test statistics T by this cycle specificity collection of illustrative plates
p
With described test statistics T
pγ compares with threshold value, if T
p≤ γ is single sideband modulated signal, otherwise is the double-sideband modulation signal.
A kind of device of sideband modulated signal classification comprises Anneta module, low noise amplifier, band pass filter, analog to digital converter and time-domain signal pretreatment module, and this device also comprises the signal sort module, and this module comprises:
Enhancement mode frequency spectrum correlation function module is used for according to time-domain digital signal x[n] the frequency spectrum correlation function
With the conjugate spectrum correlation function
Obtain enhancement mode frequency spectrum correlation function
The test statistics computing module is used for according to described enhancement mode frequency spectrum correlation function
Calculate described time-domain digital signal x[n] the cycle specificity collection of illustrative plates, and obtain test statistics T by this cycle specificity collection of illustrative plates
p
Comparison module is used for described test statistics T
pγ compares with threshold value, if T
p≤ γ is single sideband modulated signal, otherwise is the double-sideband modulation signal.
As seen from the above technical solution provided by the invention, by forming enhancement mode frequency spectrum correlation function in conjunction with conjugate spectrum correlation function and frequency spectrum correlation function, obtain the cycle specificity collection of illustrative plates by enhancement mode frequency spectrum correlation function, thereby realize carrying out the classification of sideband modulated signal under unknown carrier frequency information condition, enlarge the usage range of sideband signals classification and Detection, and had advantages of low complex degree.
Description of drawings
In order to be illustrated more clearly in the technical scheme of the embodiment of the present invention, during the below will describe embodiment, the accompanying drawing of required use is done to introduce simply, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite of not paying creative work, can also obtain other accompanying drawings according to these accompanying drawings.
The flow chart of the method for a kind of sideband modulated signal classification that Fig. 1 provides for the embodiment of the present invention one;
The enhancement mode frequency spectrum correlation function mould schematic three dimensional views of the single sideband modulated signal that Fig. 2 provides for the embodiment of the present invention one;
The enhancement mode frequency spectrum correlation function mould schematic three dimensional views of the double-sideband modulation signal that Fig. 3 provides for the embodiment of the present invention one;
The schematic diagram of the cycle specificity figure of the single sideband modulated signal that Fig. 4 provides for the embodiment of the present invention one;
The schematic diagram of the cycle specificity figure of the double-sideband modulation signal that Fig. 5 provides for the embodiment of the present invention one;
The schematic diagram of the device of a kind of sideband modulated signal classification that Fig. 6 provides for the embodiment of the present invention two;
The schematic diagram of the device of the another kind of sideband modulated signal classification that Fig. 7 provides for the embodiment of the present invention two;
The schematic diagram of the device of another sideband modulated signal classification that Fig. 8 provides for the embodiment of the present invention two.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is only the present invention's part embodiment, rather than whole embodiment.Based on embodiments of the invention, those of ordinary skills belong to protection scope of the present invention not making the every other embodiment that obtains under the creative work prerequisite.
Embodiment one
The flow chart of the method for a kind of sideband modulated signal classification that Fig. 1 provides for the embodiment of the present invention one.As shown in Figure 1, mainly comprise the steps:
Described pretreated time-domain digital signal x[n] can assist acquisition by hardware device.For example, at first, utilize Anneta module to receive aerial signal; Then, carry out to the received signal low noise amplification and the signal after amplifying is carried out bandpass filtering; Again filtering signal out is sent to and carries out analog-to-digital conversion in analog to digital converter, at last, digital signal after analog-to-digital conversion is sent into the signal procesing in time domain module, complete in this module processing such as the reception of signal and storages, and then obtain pretreated time-domain digital signal x[n].
This time-domain digital signal x[n] can be that single sideband modulated signal can be also the double-sideband modulation signal.For example and without limitation, we are with x
m[n] expression intermediate frequency modulated-analog signal, can set up following dualism hypothesis:
H
1Single sideband modulated signal:
H
0The double-sideband modulation signal:
Wherein, x
s[n] is by x
m[n] carries out the Hilbert(Hilbert) conversion obtains, f
cCarrier frequency, T
sBe sample frequency,
Be [0,2 π) upper equally distributed phase place.
Obtaining enhancement mode frequency spectrum correlation function
Need to calculate in advance x[n before] the frequency spectrum correlation function
X[n in the embodiment of the present invention] (n=1,2 ..., frequency spectrum correlation function N-1)
Can adopt following formula to calculate obtains:
Wherein, M is positive odd number and M≤N, and α is cycle frequency, and X (f) is x[] discrete time Fourier transform (DTFT) and
The conjugation of expression signal namely for complex signal, obtains after the sign-flip with its imaginary part.
When f ± the α in following formula/2nd, the integral multiple of 1/N, can calculate X (f) by fast Fourier transform (FFT).But the frequency spectrum correlation function that the method calculates can not effectively be showed the difference between single sideband modulated signal and double-sideband modulation signal, so the embodiment of the present invention proposes the conjugate spectrum correlation function
And on this basis in conjunction with frequency spectrum correlation function structure enhancement mode frequency spectrum correlation function
Wherein,
Computing formula as follows:
Wherein, M is positive odd number and M≤N; X (f) is x[] discrete time Fourier transform (DTFT) and
The conjugation of expression signal namely for complex signal, obtains after the sign-flip with its imaginary part.
At last, with x[n] the frequency spectrum correlation function
With the conjugate spectrum correlation function
Combine and obtain enhancement mode frequency spectrum correlation function
Wherein, COEFFICIENT K
1With K
2For greater than 0 real number.
The enhancement mode frequency spectrum correlation function that obtains by the way can effectively be showed the difference between single sideband modulated signal and double-sideband modulation signal, as shown in Fig. 2-3, the graphics of the enhancement mode frequency spectrum correlation function of single sideband modulated signal and double-sideband modulation signal has been described respectively, wherein, signal to noise ratio snr=5dB in Fig. 2-Fig. 3, κ
1=0.1, κ
2=1.Can be from visually seeing both differences, the double-sideband modulation signal in Fig. 3 has an obvious spike in the middle of amplitude larger continuum, and the single sideband modulated signal in Fig. 2 is without this feature.
The cycle specificity collection of illustrative plates can have multiple design, for reducing the data volume of required processing, only gets enhancement mode frequency spectrum correlation function in the embodiment of the present invention
Amplitude peak on each given cycle frequency point, corresponding cycle specificity collection of illustrative plates is defined as:
Wherein,
Be define symbol; ξ is frequency threshold, namely only needs the calculated rate absolute value less than the enhancement mode frequency spectrum correlation function of this value.
As shown in Fig. 4-5, be respectively the cycle specificity collection of illustrative plates of single sideband modulated signal and double-sideband modulation signal, wherein, the signal to noise ratio snr=5dB in Fig. 4-Fig. 5, κ
1=0.1, κ
2=1.It is equivalent to the three-dimensional enhancement mode frequency spectrum correlation function of Fig. 2-Fig. 3 has been done a projection in the cycle frequency territory.
For the ease of embody the difference of single sideband modulated signal and double-sideband modulation signal cycle characteristic spectrum from numerical value, the embodiment of the present invention adopts following test statistics:
Wherein, η is the maximum of cycle specificity figure spectral amplitude, and e is the truth of a matter of natural logrithm, and τ is the symmetry coefficient of cycle specificity collection of illustrative plates.For ease of understanding the implication of τ in above-mentioned formula, can use a
lExpression cycle specificity collection of illustrative plates is used a in the amplitude equalizing value on the amplitude maximum of points left side
rExpression the right amplitude equalizing value:
τ=|a
l-a
r|/(a
l+a
r)。
Wherein, threshold value γ can also can be according to the monte carlo method analog computation according to the actual numerical value analog computation; Described analog computation comprises according to monte carlo method: obtained the statistical distribution characteristic of the test statistics of different signal to noise ratio lower sideband modulation signals by monte carlo method emulation, then ask for threshold value γ under given false alarm rate (referring to that the probability of false-alarm appears in detection system when detecting).
The test statistics T that after calculating threshold value γ, itself and step 14 is calculated
pCompare.Generally, the T of single sideband modulated signal
pLess, the T of double-sideband modulation signal
pRelatively large; Therefore, work as T
pDuring≤γ, judge that current demand signal is single sideband modulated signal, otherwise be the double-sideband modulation signal.
The embodiment of the present invention is carried out the sideband modulated signal sorting technique based on enhancement mode frequency spectrum correlation function, can carry out under different signal to noise ratio snr, in emulation experiment, when SNR=-10dB, still can keep high verification and measurement ratio.
The embodiment of the present invention is by forming enhancement mode frequency spectrum correlation function in conjunction with conjugate spectrum correlation function and frequency spectrum correlation function, obtain the cycle specificity collection of illustrative plates by enhancement mode frequency spectrum correlation function, thereby realize carrying out the classification of sideband modulated signal under unknown carrier frequency information condition, enlarge the usage range of sideband signals classification and Detection, and had advantages of low complex degree.
Through the above description of the embodiments, those skilled in the art can be well understood to above-described embodiment and can realize by software, also can realize by the mode that software adds necessary general hardware platform.Based on such understanding, the technical scheme of above-described embodiment can embody with the form of software product, it (can be CD-ROM that this software product can be stored in a non-volatile memory medium, USB flash disk, portable hard drive etc.) in, comprise some instructions with so that computer equipment (can be personal computer, server, the perhaps network equipment etc.) carry out the described method of each embodiment of the present invention.
Embodiment two
The schematic diagram of the device of a kind of sideband modulated signal classification that Fig. 6 provides for the embodiment of the present invention.As shown in Figure 6, mainly comprise: Anneta module 71, low noise amplifier 72, band pass filter 73, analog to digital converter 74 and time-domain signal pretreatment module 75, it is characterized in that, this device also comprises signal sort module 76, and this module comprises:
Enhancement mode frequency spectrum correlation function module 761 is used for according to time-domain digital signal x[n] the frequency spectrum correlation function
With the conjugate spectrum correlation function
Obtain enhancement mode frequency spectrum correlation function
Test statistics computing module 762 is used for according to described enhancement mode frequency spectrum correlation function
Calculate described time-domain digital signal x[n] the cycle specificity collection of illustrative plates, and obtain test statistics T by this cycle specificity collection of illustrative plates
p
Comparison module 763 is used for described test statistics T
pγ compares with threshold value, if T
p≤ γ is single sideband modulated signal, otherwise is the double-sideband modulation signal.
Described signal sort module 76 also comprises:
Threshold value computing module 764 is used for according to monte carlo method analog computation threshold value or according to actual numerical value analog computation threshold value.
In addition, the signal sort module 76 in the present embodiment can be in the PC(computer) in, referring to Fig. 7, namely Anneta module 71, low noise amplifier 72, band pass filter 73, analog to digital converter 74, time-domain signal pretreatment module 75 are connected successively with PC77.also signal sort module 76 can be placed in the DSP78(digital signal processor), referring to Fig. 8, it is Anneta module 71, low noise amplifier 72, band pass filter 73, analog to digital converter 74, time-domain signal pretreatment module 75 is connected successively with DSP78, this DSP78 is sent to by corresponding interface (as the EMIF mouth) data that time-domain signal pretreatment module 75 sends in data storage 781, DSP78 in real time from memory 781 reading out data complete classification and the Output rusults of sideband modulated signal, result for output can be delivered in the external equipment 782 that is connected with DSP78, also can pass through some interfaces such as USB mouth 783, network interface 784 grades are transferred in other equipment.
Need to prove, have a detailed description in the specific implementation embodiment in front of the function that each processing module that comprises in said apparatus realizes, therefore here repeat no more.
The those skilled in the art can be well understood to, be the convenience described and succinct, only the division with above-mentioned each functional module is illustrated, in practical application, can as required above-mentioned functions be distributed and be completed by different functional modules, the internal structure that is about to device is divided into different functional modules, to complete all or part of function described above.
The above; only for the better embodiment of the present invention, but protection scope of the present invention is not limited to this, anyly is familiar with those skilled in the art in the technical scope that the present invention discloses; the variation that can expect easily or replacement are within all should being encompassed in protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection range of claims.
Claims (9)
1. the method for a sideband modulated signal classification, is characterized in that, the method comprises:
According to time-domain digital signal x[n] the frequency spectrum correlation function
With the conjugate spectrum correlation function
Obtain enhancement mode frequency spectrum correlation function
According to described enhancement mode frequency spectrum correlation function
Calculate described time-domain digital signal x[n] the cycle specificity collection of illustrative plates, and obtain test statistics T by this cycle specificity collection of illustrative plates
p
With described test statistics T
pγ compares with threshold value, if T
p≤ γ is single sideband modulated signal, otherwise is the double-sideband modulation signal.
2. method according to claim 1, is characterized in that, calculates described time-domain digital signal x[n] the frequency spectrum correlation function
Formula be:
3. method according to claim 1, is characterized in that, calculates described time-domain digital signal x[n] the conjugate spectrum correlation function
Formula be:
4. according to claim 1-3 described methods of any one, is characterized in that, calculates described enhancement mode frequency spectrum correlation function
Formula comprise:
Wherein, COEFFICIENT K
1With K
2For greater than 0 real number.
5. method according to claim 1, is characterized in that, and is described according to described enhancement mode frequency spectrum correlation function
Calculate described time-domain digital signal x[n] the formula of cycle specificity collection of illustrative plates be:
Wherein, ξ is frequency threshold.
6. method according to claim 1, is characterized in that, and is described by this cycle specificity collection of illustrative plates acquisition test statistics T
pFormula be:
Wherein, η be cycle specificity figure spectral amplitude; Large value, τ is the symmetry coefficient of cycle specificity collection of illustrative plates, e is the truth of a matter of natural logrithm.
7. method according to claim 1, is characterized in that, the method for calculating described threshold value γ comprises:
Obtained the statistical distribution characteristic of the test statistics of different signal to noise ratio lower sideband modulation signals by monte carlo method emulation, then ask for threshold value under predetermined false alarm rate or according to actual numerical value analog computation threshold value.
8. the device of a sideband modulated signal classification, comprise Anneta module, low noise amplifier, band pass filter, analog to digital converter and time-domain signal pretreatment module, it is characterized in that, this device also comprises the signal sort module, and this module comprises:
Enhancement mode frequency spectrum correlation function module is used for according to time-domain digital signal x[n] the frequency spectrum correlation function
With the conjugate spectrum correlation function
Obtain enhancement mode frequency spectrum correlation function
The test statistics computing module is used for according to described enhancement mode frequency spectrum correlation function
Calculate described time-domain digital signal x[n] the cycle specificity collection of illustrative plates, and obtain test statistics T by this cycle specificity collection of illustrative plates
p
Comparison module is used for described test statistics T
pγ compares with threshold value, if T
p≤ γ is single sideband modulated signal, otherwise is the double-sideband modulation signal.
9. device according to claim 8, is characterized in that, described signal sort module also comprises:
The threshold value computing module is used for according to monte carlo method analog computation threshold value or according to actual numerical value analog computation threshold value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310076546.0A CN103166887B (en) | 2013-03-11 | 2013-03-11 | A kind of method of sideband modulated signal classification and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310076546.0A CN103166887B (en) | 2013-03-11 | 2013-03-11 | A kind of method of sideband modulated signal classification and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103166887A true CN103166887A (en) | 2013-06-19 |
CN103166887B CN103166887B (en) | 2015-11-18 |
Family
ID=48589635
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310076546.0A Expired - Fee Related CN103166887B (en) | 2013-03-11 | 2013-03-11 | A kind of method of sideband modulated signal classification and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103166887B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105812040A (en) * | 2016-03-21 | 2016-07-27 | 中国科学技术大学 | Method for transmitting signals to multiple users by using wireless channel reciprocity |
CN107135176A (en) * | 2017-07-06 | 2017-09-05 | 电子科技大学 | Figure field communication signal modulate method based on fractional lower-order Cyclic Spectrum |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2765198Y (en) * | 2005-01-13 | 2006-03-15 | 阮志彬 | Band-limited pseudocode modulation signal source in RF identification system |
EP2299642A2 (en) * | 2009-09-14 | 2011-03-23 | Siemens Aktiengesellschaft | A method meeting the spectral mask in a radiofrequency identification system and a device therefor |
-
2013
- 2013-03-11 CN CN201310076546.0A patent/CN103166887B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2765198Y (en) * | 2005-01-13 | 2006-03-15 | 阮志彬 | Band-limited pseudocode modulation signal source in RF identification system |
EP2299642A2 (en) * | 2009-09-14 | 2011-03-23 | Siemens Aktiengesellschaft | A method meeting the spectral mask in a radiofrequency identification system and a device therefor |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105812040A (en) * | 2016-03-21 | 2016-07-27 | 中国科学技术大学 | Method for transmitting signals to multiple users by using wireless channel reciprocity |
CN105812040B (en) * | 2016-03-21 | 2018-12-14 | 中国科学技术大学 | A method of using wireless channel reciprocity to multi-user transmission signal |
CN107135176A (en) * | 2017-07-06 | 2017-09-05 | 电子科技大学 | Figure field communication signal modulate method based on fractional lower-order Cyclic Spectrum |
CN107135176B (en) * | 2017-07-06 | 2019-12-27 | 电子科技大学 | Image domain communication signal modulation identification method based on fractional low-order cyclic spectrum |
Also Published As
Publication number | Publication date |
---|---|
CN103166887B (en) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103607361B (en) | The parameter prediction method of time-frequency overlapped signal under Alpha Stable distritation noise | |
CN103457890A (en) | Method for effectively recognizing digital modulating signals in non-Gaussian noise | |
CN104729853A (en) | Rolling bearing performance degradation evaluation device and method | |
CN102830409B (en) | Navigation signal acquiring method based on compressed sensing | |
CN103873170B (en) | A kind of frequency spectrum detecting method of compressed sensing under blind sparse condition | |
CN104535905A (en) | Partial discharge diagnosis method based on naive bayesian classification | |
CN103630916B (en) | Code capturing method based on dual-FFT frequency domain filtering | |
CN110187313B (en) | Radar signal sorting and identifying method and device based on fractional order Fourier transform | |
CN102508031A (en) | Fourier series based measurement method of phase angle of partial discharge pulse | |
CN105515695A (en) | Compressed sampling signal detection method based on modulated wideband converter | |
CN105137297A (en) | Method and device for separating multi-source partial discharge signals of power transmission device | |
CN107037486A (en) | The Time-frequency Spectrum Analysis method and system of earth natural pulses electromagnetic field data processing | |
CN106771456A (en) | Digital oscilloscope triggering method and system based on Power Spectral Entropy | |
CN102914325A (en) | Dissipation synchronization-based detection method of small signal under chaos background | |
CN106526566A (en) | Signal in-pulse characteristic real-time analysis and processing method based on FPGA high-speed preprocessing | |
CN102546499A (en) | Fractional-order channelized receiving method of real linear frequency modulation (LFM) signal | |
CN102158248A (en) | Estimation method of DS (direct sequence)/FH (frequency hopping) spread spectrum signal parameters based on cyclic spectrum theory | |
CN103166887B (en) | A kind of method of sideband modulated signal classification and device | |
CN105721090B (en) | A kind of detection and recognition methods of illegal f-m broadcast station | |
CN103308829A (en) | GIS (gas insulated switchgear) single-time partial discharge signal extraction and trigger moment adjusting method | |
CN105242286A (en) | Capture method based on wavelet domain filtering code | |
CN103929256B (en) | A kind of multiframe compressed sensing signal spectrum detection method | |
CN112881812A (en) | Full-flash real-time positioning method and device based on machine learning coding | |
EP3055704B1 (en) | Analog to information converter | |
Liang et al. | An individual emitter recognition method combining bispectrum with wavelet entropy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151118 |