CN112804172A - Method for realizing signal modulation mode identification based on high-order cumulant - Google Patents
Method for realizing signal modulation mode identification based on high-order cumulant Download PDFInfo
- Publication number
- CN112804172A CN112804172A CN202011556395.5A CN202011556395A CN112804172A CN 112804172 A CN112804172 A CN 112804172A CN 202011556395 A CN202011556395 A CN 202011556395A CN 112804172 A CN112804172 A CN 112804172A
- Authority
- CN
- China
- Prior art keywords
- signal
- modulation mode
- order
- order cumulant
- conversion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0012—Modulated-carrier systems arrangements for identifying the type of modulation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
The invention provides a method for realizing signal modulation mode identification based on high-order cumulant, which comprises the following steps: 1. data preprocessing, namely completing analog-to-digital conversion of a received signal, performing FFT (fast Fourier transform) processing and completing rough estimation of a signal carrier frequency; 2. performing down-conversion, namely performing down-conversion processing on the signal through a carrier frequency estimation result of data preprocessing to obtain a complex baseband signal (I, Q paths of data); 3. calculating high-order cumulant of the obtained I, Q paths of data; 4. and comparing the obtained high-order cumulant with a preset threshold value, and judging the modulation mode of the output signal. According to the method, two groups of high-order cumulants of the signals are calculated, and at most three times of comparison and judgment are carried out, so that the identification of a signal modulation mode based on the high-order cumulants is realized, compared with the traditional signal modulation mode identification algorithm based on the high-order cumulants, the number of signal characteristics (high-order cumulants) is less, and a high-precision identification result can be realized through fewer judgment times.
Description
Technical Field
The invention relates to the field of communication, in particular to a method for realizing signal modulation mode identification based on high-order cumulant.
Background
In the electronic reconnaissance activity, the following interference and attack measures can be implemented only by identifying the modulation mode of the intercepted non-partner signal and then demodulating and analyzing the modulation mode. At present, three groups of high-order cumulants are generally required to be selected for classifying 8 modulation signals by adopting a high-order cumulant algorithm, four times of judgment can be realized, and the algorithm is complex and real-time modulation mode identification is difficult to realize. Therefore, it is an urgent problem to develop a method for realizing signal modulation mode identification based on high-order cumulant.
Disclosure of Invention
In view of this, the present invention aims to provide a method for identifying signal modulation modes based on high-order cumulant, which realizes identification of 8 signal modulation modes through at most 3 times of decision.
In order to solve the technical problems, the invention adopts the technical scheme that:
a method for realizing signal modulation mode identification based on high-order cumulant comprises the following steps:
1. data preprocessing, namely completing analog-to-digital conversion of a received signal, performing FFT (fast Fourier transform) processing and completing rough estimation of a signal carrier frequency;
2. performing down-conversion, namely performing down-conversion processing on the signal through a carrier frequency estimation result of data preprocessing to obtain a complex baseband signal (I, Q paths of data);
3. calculating high-order cumulant of the obtained I, Q paths of data;
4. and comparing the obtained high-order cumulant with a preset threshold value, and judging the modulation mode of the output signal.
Preferably, the method for realizing the identification of the signal modulation mode based on the high-order cumulant comprises the following specific steps:
at a signal receiving end, a received signal is preprocessed, and assuming that a carrier frequency, a phase and a timing of the received signal are synchronous, a complex baseband signal after down-conversion processing can be represented as:
X(t)＝Σak√Esq(t-kT)ej(Δθc)+n(t)
where k is 1,2,3, … N, and N is the number of transmission symbol sequencesTo that end, q (t) is the baseband symbol waveform, akIs a sequence of symbols, T is a symbol width, EsTo transmit the energy of the symbol waveform, Δ θcFor the carrier phase difference, n (t) is white gaussian noise with zero mean.
This time, a stationary complex random signal x (t) is obtained, for which the mixing moment of order k is calculated:
Mpq＝E[X(t)(p-q)X*(t)q]
the conversion relationship between the higher order moments and the higher order cumulants of the random signal is as follows:
Cum[X1,X2,…,Xk]＝Σ(-1)q-1(q-1)！∏E{∏Xi}
where Σ (·) denotes X ═ X (X) in all the mutually unconnected ordered partition sets1,X2,…,Xk) Inner summation, q is the number of subsets of the partition, UpA subscript set representing the pth subset element among the q subsets.
The invention selects C41，C80Two high-order cumulants are calculated by the following specific formula:
C41＝M41-3M20M21
C80＝M80-28M20M60-35M40 2+420M20 2M40-630M20 4
comparing and judging the calculated high-order cumulant with a set threshold value, wherein C41_1＝3；C41_2＝1.5；C41_3＝0.5；C41_4＝0.2；C41_5＝0.1；C80_1＝10；C80_220; and finally, outputting the modulation mode of the signal according to the judgment result.
Compared with the prior art, the method for realizing the signal modulation mode identification based on the high-order cumulant has the following advantages: according to the method, two groups of high-order cumulants of the signals are calculated, and at most three times of comparison and judgment are carried out, so that the identification of a signal modulation mode based on the high-order cumulants is realized, compared with the traditional signal modulation mode identification algorithm based on the high-order cumulants, the number of signal characteristics (high-order cumulants) is less, and a high-precision identification result can be realized through fewer judgment times.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate an embodiment of the invention and, together with the description, serve to explain the invention and not to limit the invention. In the drawings:
fig. 1 is a schematic diagram of the present invention.
Fig. 2 is a decision logic proposed by the present invention.
Fig. 3 is a simulation calculation result of the algorithm proposed by the present invention for each signal recognition success rate.
Fig. 4 is a simulation calculation result of the algorithm of the present invention for the average recognition success rate of 8 signals.
Detailed Description
It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.
In the description of the present invention, it is to be understood that the terms "center", "longitudinal", "lateral", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", and the like, indicate orientations or positional relationships based on those shown in the drawings, and are used only for convenience in describing the present invention and for simplicity in description, and do not indicate or imply that the referenced devices or elements must have a particular orientation, be constructed and operated in a particular orientation, and thus, are not to be construed as limiting the present invention. Furthermore, the terms "first", "second", etc. are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first," "second," etc. may explicitly or implicitly include one or more of that feature. In the description of the present invention, "a plurality" means two or more unless otherwise specified.
In the description of the present invention, it should be noted that, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meaning of the above terms in the present invention can be understood by those of ordinary skill in the art through specific situations.
The present invention will be described in detail below with reference to the embodiments with reference to the attached drawings.
As shown in fig. 1 to 4, a method for implementing signal modulation mode identification based on high-order cumulant includes the following steps:
1. data preprocessing, namely completing analog-to-digital conversion of a received signal, performing FFT (fast Fourier transform) processing and completing rough estimation of a signal carrier frequency;
2. performing down-conversion, namely performing down-conversion processing on the signal through a carrier frequency estimation result of data preprocessing to obtain a complex baseband signal (I, Q paths of data);
3. calculating high-order cumulant C41 and C80 of the obtained I, Q paths of data;
4. and comparing the obtained high-order cumulant with a preset threshold value, and judging the modulation mode of the output signal.
The specific embodiment of the method is as follows:
complex representations of digitally modulated signals that are typically disturbed by gaussian noise:
X(t)＝Σak√Esq(t-kT)ej(2πfct+θc)+n(t)
where k is 1,2,3, … N, N is the number of transmitted symbol sequences, q (t) is the baseband symbol waveform, akIs a sequence of symbols, T is a symbol width, EsFor transmitting the energy of the symbol waveform, fcIs the carrier frequency, θcIs the initial phase of the carrier, and n (t) is white Gaussian noise with zero mean.
At a signal receiving end, a received signal is preprocessed, and assuming that a carrier frequency, a phase and a timing of the received signal are synchronous, a complex baseband signal after down-conversion processing can be represented as:
X(t)＝Σak√Esq(t-kT)ej(Δθc)+n(t)
in the formula,. DELTA.theta.cFor the carrier phase difference, n (t) is white gaussian noise with zero mean.
This time, a stationary complex random signal x (t) is obtained, for which the mixing moment of order k is calculated:
Mpq＝E[X(t)(p-q)X*(t)q]
the conversion relationship between the higher order moments and the higher order cumulants of the random signal is as follows:
Cum[X1,X2,…,Xk]＝Σ(-1)q-1(q-1)！∏E{∏Xi}
where Σ (·) denotes X ═ X (X) in all the mutually unconnected ordered partition sets1,X2,…,Xk) Inner summation, q is the number of subsets of the partition, UpA subscript set representing the pth subset element among the q subsets.
The invention selects C41，C80Two high-order cumulants are calculated by the following specific formula:
C41＝M41-3M20M21
C80＝M80-28M20M60-35M40 2+420M20 2M40-630M20 4
comparing and judging the calculated high-order cumulant with a set threshold value, wherein C41_1＝3；C41_2＝1.5；C41_3＝0.5；C41_4＝0.2；C41_5＝0.1；C80_1＝10；C80_220; and finally, outputting the modulation mode of the signal according to the judgment result.
By adopting the scheme, the identification of the modulation mode of the received signal can be realized. Through calculation of two groups of high-order cumulants of the signals, at most three times of comparison and judgment are carried out, and specific judgment logics are shown in fig. 2, the modulation mode identification of 8 modulation signals (2ASK, 2FSK, 2PSK, 4ASK, 4FSK, 4PSK, 8PSK and 16QAM) can be realized, and the carrier frequency is 10MHz, the sampling frequency is 60MHz and the symbol rate is 2.5MHz under simulation conditions. The signal has white Gaussian noise, and the signal-to-noise ratio range is [ -2dB,16dB ]; and each signal is independently subjected to 1000 Monte Carlo experiments, and when the signal-to-noise ratio is higher than 2dB, the identification precision of each signal can reach more than 85%. The average recognition success rate of the algorithm is over 97 percent.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (3)
1. A method for realizing signal modulation mode identification based on high-order cumulant is characterized by comprising the following steps:
1) data preprocessing, namely performing analog-to-digital conversion on a received signal, performing FFT (fast Fourier transform) processing, and completing coarse estimation on a signal carrier frequency;
2) performing down-conversion, namely performing down-conversion processing on the signal through a carrier frequency estimation result of data preprocessing to obtain a complex baseband signal, namely I, Q paths of data;
3) calculating high-order cumulant of the obtained I, Q paths of data;
4) and comparing the obtained high-order cumulant with a preset threshold value, and judging the modulation mode of the output signal.
2. The method for realizing signal modulation mode identification based on the high-order cumulant as claimed in claim 1, wherein:
at a signal receiving end, firstly, preprocessing a received signal, assuming that a carrier frequency, a phase and a timing of the received signal are synchronous, and a complex baseband signal after down-conversion processing is represented as:
X(t)＝Σak√Esq(t-kT)ej(Δθc)+n(t)
where k is 1,2,3, … N, N is the number of transmitted symbol sequences, q (t) is the baseband symbol waveform,akis a sequence of symbols, T is a symbol width, EsTo transmit the energy of the symbol waveform, Δ θcWhite Gaussian noise with zero mean value of n (t) which is the carrier phase difference;
this time, a stationary complex random signal x (t) is obtained, for which the mixing moment of order k is calculated:
Mpq＝E[X(t)(p-q)X*(t)q]
the conversion relationship between the higher order moments and the higher order cumulants of the random signal is as follows:
Cum[X1,X2,…,Xk]＝Σ(-1)q-1(q-1)！∏E{∏Xi}
where Σ (·) denotes X ═ X (X) in all the mutually unconnected ordered partition sets1,X2,…,Xk) Inner summation, q is the number of subsets of the partition, UpA subscript set representing a pth subset element among the q subsets;
selecting C41，C80Two high-order cumulants are calculated by the following specific formula:
C41＝M41-3M20M21
C80＝M80-28M20M60-35M40 2+420M20 2M40-630M20 4
and comparing and judging the calculated high-order cumulant with a set threshold value, and finally outputting a modulation mode of the signal according to a judgment result.
3. The method for realizing signal modulation mode identification based on the high-order cumulant as claimed in claim 2, wherein: c41_1＝3；C41_2＝1.5；C41_3＝0.5；C41_4＝0.2；C41_5＝0.1；C80_1＝10；C80_2＝20。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011556395.5A CN112804172A (en) | 2020-12-24 | 2020-12-24 | Method for realizing signal modulation mode identification based on high-order cumulant |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011556395.5A CN112804172A (en) | 2020-12-24 | 2020-12-24 | Method for realizing signal modulation mode identification based on high-order cumulant |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112804172A true CN112804172A (en) | 2021-05-14 |
Family
ID=75805811
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011556395.5A Pending CN112804172A (en) | 2020-12-24 | 2020-12-24 | Method for realizing signal modulation mode identification based on high-order cumulant |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112804172A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113852384A (en) * | 2021-09-27 | 2021-12-28 | 广西万云科技有限公司 | Method for enhancing signal-to-noise ratio of radio frequency signal |
| CN114422016A (en) * | 2021-12-28 | 2022-04-29 | 成都爱科特科技发展有限公司 | PCMA signal identification method |
| CN114760174A (en) * | 2022-06-14 | 2022-07-15 | 北京邮电大学 | Modulation recognition method and device, electronic equipment and storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101909035A (en) * | 2010-07-16 | 2010-12-08 | 北京航空航天大学 | Method for recognizing subcarrier modulation modes of orthogonal frequency division multiplexing signal in wireless communication |
| CN107948107A (en) * | 2017-11-16 | 2018-04-20 | 成都玖锦科技有限公司 | A kind of digital modulation signals sorting technique based on union feature |
| CN110958199A (en) * | 2019-12-03 | 2020-04-03 | 上海电机学院 | Communication modulation signal identification method based on multi-dimensional feature extraction |
-
2020
- 2020-12-24 CN CN202011556395.5A patent/CN112804172A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101909035A (en) * | 2010-07-16 | 2010-12-08 | 北京航空航天大学 | Method for recognizing subcarrier modulation modes of orthogonal frequency division multiplexing signal in wireless communication |
| CN107948107A (en) * | 2017-11-16 | 2018-04-20 | 成都玖锦科技有限公司 | A kind of digital modulation signals sorting technique based on union feature |
| CN110958199A (en) * | 2019-12-03 | 2020-04-03 | 上海电机学院 | Communication modulation signal identification method based on multi-dimensional feature extraction |
Non-Patent Citations (3)
| Title |
|---|
| 任玉静: "常用数字信号调制样式的自动识别", 《中国优秀硕士学位论文全文库 信息科技辑》 * |
| 任玉静: "常用数字信号调制样式的自动识别，2010年第11期", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 王柳: "宽带无线电信号中的调制识别", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113852384A (en) * | 2021-09-27 | 2021-12-28 | 广西万云科技有限公司 | Method for enhancing signal-to-noise ratio of radio frequency signal |
| CN114422016A (en) * | 2021-12-28 | 2022-04-29 | 成都爱科特科技发展有限公司 | PCMA signal identification method |
| CN114422016B (en) * | 2021-12-28 | 2024-01-02 | 成都爱科特科技发展有限公司 | PCMA signal identification method |
| CN114760174A (en) * | 2022-06-14 | 2022-07-15 | 北京邮电大学 | Modulation recognition method and device, electronic equipment and storage medium |
| CN114760174B (en) * | 2022-06-14 | 2022-09-13 | 北京邮电大学 | Modulation recognition method and device, electronic equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112804172A (en) | Method for realizing signal modulation mode identification based on high-order cumulant | |
| CN101083649B (en) | Method for identifying OFDM modulation system of multi-path Rayleigh fast fading channel | |
| CN101834819B (en) | Analog-digital mixing modulation recognition device and digital modulation recognition device based on parallel judgment | |
| CN108627818B (en) | OFDM-based frequency control array radar communication integrated waveform design method | |
| CN107612867B (en) | MQAM signal modulation order identification method | |
| CN112363132A (en) | FBMC-based radar communication integrated waveform generation method | |
| CN1750523A (en) | Method for estimating channel quality of continuous phase modulation self adaptive frequency-hopping system | |
| CN113612716B (en) | Multi-user communication radar integrated system based on constant envelope orthogonal frequency division multiplexing | |
| CN112751798B (en) | NMW-OFDM system for terahertz communication perception integration | |
| CN106911605A (en) | A kind of radar-communication integration Waveform generation method based on FMCW | |
| CN104158633A (en) | Maximum likelihood modulation recognition method based on Gaussian mixture model | |
| CN109787932A (en) | A kind of radar-communication integration Design of Signal method based on OFDM | |
| CN113612708B (en) | Multi-user communication radar integrated system based on time division multiple access | |
| CN101442512B (en) | Method for recognizing OFDM signal | |
| CN101478524B (en) | Correction method for phase noise of multi-antenna OFDM system receiving terminal | |
| Wu et al. | A non-uniform multi-wideband OFDM system for terahertz joint communication and sensing | |
| CN101909035A (en) | Method for recognizing subcarrier modulation modes of orthogonal frequency division multiplexing signal in wireless communication | |
| CN101032076B (en) | Method for estimating signal magnitude, noise power, and signal-to-noise ratio of received signals | |
| CN111131107B (en) | Self-adaptive soft demodulation method based on 5G downlink shared channel state | |
| CN103546418A (en) | Method and device for recognizing modulation modes on basis of peak-to-average ratios | |
| Alyaoui et al. | The modulation recognition approaches for software radio | |
| CN111313974A (en) | Modulation format identification method for signal amplitude sequencing and ratio calculation | |
| JP2003244263A (en) | Signal processor | |
| JP2003032221A (en) | Measurement equipment for transmission line characteristics | |
| CN103685112A (en) | Phase jump-based recognition method and device for signal modulation modes |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210514 |