CN110933011B - High-resolution blind scanning method for WFRFT (WFRFT) signals - Google Patents
High-resolution blind scanning method for WFRFT (WFRFT) signals Download PDFInfo
- Publication number
- CN110933011B CN110933011B CN201911263420.8A CN201911263420A CN110933011B CN 110933011 B CN110933011 B CN 110933011B CN 201911263420 A CN201911263420 A CN 201911263420A CN 110933011 B CN110933011 B CN 110933011B
- Authority
- CN
- China
- Prior art keywords
- channel
- modulation order
- wfrft
- calculating
- beta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2626—Arrangements specific to the transmitter only
- H04L27/2627—Modulators
- H04L27/2628—Inverse Fourier transform modulators, e.g. inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2649—Demodulators
- H04L27/265—Fourier transform demodulators, e.g. fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators
Landscapes
- Physics & Mathematics (AREA)
- Discrete Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
A high resolution blind scanning method of WFRFT signal belongs to the technical field of receiving method of WFRFT communication system, the method comprises establishing rule relation of modulation order error and error rate under theoretical condition, determining minimum main value window; setting a left modulation order and a right modulation order in each channel, further performing WFRFT inverse processing of two branches on a received signal by each channel, calculating a peak-to-average ratio to obtain a theoretical boundary value, calculating a mapping channel, and further calculating a demodulation order corresponding to the channel; calculating an estimated modulation order; and performing WFRFT inverse processing by using the estimated modulation order and the received signal, and obtaining demodulation data. The invention solves the scanning problem under the unknown modulation order condition and the fuzzy modulation order condition in the severe environment and the fuzzy resolution problem of the successful channel in the multi-channel scanning.
Description
Technical Field
The invention belongs to the technical field of receiving methods of WFRFT communication systems, and particularly relates to a blind scanning receiving method under the condition of unknown modulation orders.
Background
Weighted Fractional Fourier Transform (WFRFT) is a Transform domain signal processing means in wireless communication systems. The modulation order (parameter) in the general WFRFT system is known only to own users and is unknown to other users. Therefore, considering the uncertain modulation order under unknown conditions or severe environment conditions, the setting of the demodulation order at the receiving end cannot meet the requirement of correctly recovering the original data, so blind scanning processing is required. The current scanning methods are all scanning in one period of modulation order at intervals of 0.01, so the scanning times are 400 times, obviously, for a high-speed processor, 400 times of WFRFT inverse processing are not a problem, but how to distinguish exact successful channels from the demodulation results of 400 channels becomes a difficult problem. If the successful channel cannot be distinguished, the subsequent analysis, analysis and other processing are required to be carried out on the 400 groups of data all the time, and the channel which is real data cannot be distinguished after the analysis. If 400 channels of data processing are required for each cycle, a great complexity and ambiguity resolution are involved for a continuous multi-cycle or real-time data transmission system. Therefore, the blind scanning method with high resolution capability is researched to become the key of the WFRFT signal processing field.
Disclosure of Invention
The invention provides a high-resolution blind scanning method of a WFRFT signal, which discusses the shunt relevant characteristics of a left modulation order and a right modulation order by establishing a rule between a modulation order error and a theoretical bit error rate, gives a rectangular relevant interval, obtains an estimated modulation order by utilizing a rectangular relevant interval boundary and a relevant mapping channel, and further recovers real data to achieve the aim of high-resolution blind scanning.
The technical scheme is as follows:
in a high-resolution blind scanning method of WFRFT signals, firstly, a modulation order error and bit error rate rule relation under a theoretical condition is established, and a minimum principal value window is determined.
Then, a multi-channel demodulation order is established in a main value window, a left modulation order and a right modulation order are set in each channel, each channel further carries out WFRFT inverse processing of two branches on the received signal, correlation operation is carried out by using the inverse processing results of the two branches in each channel, and the peak-to-average ratio is calculated.
Further, a rectangular correlation interval is established by using the theoretical correlation degree, and a theoretical boundary value is obtained.
Then, a mapping channel is calculated by using the peak-to-average ratios of a plurality of channels, and a demodulation order corresponding to the channel is further calculated.
Secondly, the estimated modulation order is calculated by using the calculation result.
And finally, performing WFRFT inverse processing by using the estimated modulation order and the received signal, and obtaining demodulation data.
Further, a method for high resolution blind scanning of WFRFT signal according to claim 1 is specifically as follows:
establishing a regular relation between modulation order error delta beta and an error rate under a theoretical condition, setting a minimum main value window K to be 2, and setting a stepping amount to be 0.01;
in a minimum principal value period, namely 0-2 intervals, carrying out multichannel scanning by taking lambda as a stepping quantity, wherein the scanning upper limit is N:
defining a demodulation order set betaiFormula (7), wherein i ∈ [ 1N ]],
βi=i·λ (7)
Assuming the estimated modulation order is α ', the error between the demodulation order of each channel and α' is:
Δβi=α'-βi (8)
will betaiIs subjected to periodic continuation to obtainThen, the reverse rotation and the translation are carried out to obtainAndin each channel betaiOn the basis, the right modulation order beta is carried outiRLeft modulation order betaiLAnd (3) calculating:
using two branches beta per channeliRAnd betaiLPerforming WFRFT inverse processing on the received signals S (n):
two-branch inverse processing result S in each channeliR(n) and SiL(n) performing correlation operation to obtain correlation result C of each channeli(n), wherein fft (-), ifft (-) are fast fourier transform and inverse transform processing functions;
By Ci(n) calculating the peak-to-average ratioPiWhere max (·) is a maximum computation function;
ΔβRis the right boundary of the rectangular correlation interval, Δ βLIs the left boundary of the rectangular correlation interval;
using peak-to-average ratio P of N channelsiComparing with a threshold G, and further calculating to obtain a continuous sequence R (j), wherein mod (·) is a modular operation function;
and (3) performing minimum value calculation on the continuous sequence R (j) by using a minimum value calculation function min (·), and performing minimum periodic processing to obtain a mapping channel H corresponding to a rectangular correlation interval:
H=mod(min[R(j)],N) (18)
calculating the demodulation order beta corresponding to the channel by using the mapping channel HH:
By using betaHAnd the right boundary Δ β of the rectangular intervalRAnd, in conjunction with equation (8), an estimated modulation order α' can be calculated:
wherein m is an integer and is more than or equal to 0, and the alpha' is the minimum positive value in the range of [04 ] by selecting m;
using the estimated alpha' and the received signal S (n) to perform WFRFT inverse processing:
calculating to obtain data D' (n) after WFRFT inverse processing as a formula (22); if the modulation order alpha of the transmitting end is in a range of 0-2, the demodulation data D '(n) is real data x (n), and if the modulation order alpha of the transmitting end is in a range of 2-4, the demodulation data D' (n) is inverted data x (-n) of the real data;
The data D ' (n) and the data D ' (-n) after the data D ' (n) are inverted are output simultaneously.
The invention has the advantages and effects that:
a high-resolution blind scanning method for WFRFT signals aims to solve the problems of scanning under unknown modulation order conditions and fuzzy modulation order conditions in severe environments and the problem of fuzzy resolution of successful channels in the conventional multi-channel scanning. A high-resolution blind scanning method for WFRFT signals not only can not increase the complexity of a system, but also can accurately estimate the modulation order and correctly output the original data, thereby achieving the purpose of high-resolution blind scanning. The blind scanning method is not only suitable for a 4-WFRFT communication system, but also suitable for a multi-item WFRFT communication system, and is compatible with the conditions of fixed single parameters and variable single parameters.
Drawings
Fig. 1 is a diagram of the relationship between modulation order error and bit error rate rule in the high-resolution blind scanning method of WFRFT signal of the present invention.
Fig. 2 is a diagram of an implementation scheme of a high-resolution blind scanning method for a WFRFT signal according to the present invention.
Detailed Description
The invention is further explained below with reference to the figures and the examples.
Let the data information to be transmitted in the digital communication system be x (n), and the length of data transmitted per period be L, for which n belongs to [1L ]. Defining 4 kinds of state functions of X (n), X (-n) and X (-n) as 0, 1, 2 and 3 times Fourier transform results of X (n), respectively. Since all the multiple WFRFT systems can be converted into corresponding 4 WFRFT systems, the implementation of the 4-WFRFT system can be universally applied to the multiple WFRFT systems.
At the communication transmitting end, alpha is defined as the modulation order (parameter) of WFRFT process, and the weighting coefficient omega can be established by utilizing alphal(α) (l is 0,1,2,3), ω is determined as shown in formula (1)lThe modulation order α of (α) is 4-period as in the fourier transform, for which [ 04) is set as a main period, and all modulation orders can be equivalent to one [ 04) period.
Using established omegal(alpha) subjecting the data x (n) to be transmitted to alpha-order WFRFT processing, as shown in formula (2), wherein Fα(. cndot.) is a WFRFT processing function of order α.
S(n)=Fα(x(n))=ω0(α)x(n)+ω1(α)X(n)+ω2(α)x(-n)+ω3(α)X(-n) (2)
At the receiving end of the communication, the purpose is to recover the data information by performing corresponding processing on the received signal s (n). If the demodulation order of the receiving end is beta, carrying out beta-order WFRFT inverse transformation on the received signal S (n):
setting:
Δβ=α-β (4)
it can be seen that equation (3) can be written as equation (5) only when Δ β → 0, i.e., β ═ α, which indicates that the original data x (n) can be correctly recovered.
SΔ(n)=FΔβ(x(n))≈x(n) (5)
However, for the communication receiving end, the modulation order α is known only to the own user and is unknown to other users. Therefore, considering the uncertain modulation order under unknown conditions or severe environment conditions, the setting of the receiving end β cannot satisfy the condition of β ═ α, and therefore, the modulation order β needs to be scanned blindly, so that Δ β → 0 can receive data x (n) correctly, thereby providing a high-resolution blind scanning method for WFRFT signals.
In a high-resolution blind scanning method of a WFRFT signal, firstly, a modulation order error delta beta and an error rate rule relation under a theoretical condition are established, as shown in FIG. 1, the error rate does not have a linear ascending trend but has a periodic cyclic change along with the increase of the delta beta, and the cyclic period is also 4. In addition, Δ β ═ 2 is a symmetric axis in one period, and the minimum principal value window K is set to be 2. Then, the step amount is set to λ 0.01.
Based on this, in a minimum principal value period, namely 0-2 intervals, with lambda as the stepping quantity, multi-channel scanning is carried out, the scanning upper limit is N:
secondly, define the demodulation order set betaiFormula (7), wherein i ∈ [ 1N ]]There are 201 channels in total.
βi=i·λ (7)
Further, assuming that the estimated modulation order is α ', the error between the demodulation order of each channel and α' is:
Δβi=α'-βi (8)
further, beta is adjusted toiIs subjected to periodic continuation to obtainThen, the reverse rotation and the translation are carried out to obtainAndas can be seen from FIG. 1, the following steps are carried outAnd the change in alpha,andhas the lowest correlation characteristic, for this reason, the error rate curve of the channel has the beta in each channeliOn the basis, the right modulation order beta is carried outiRLeft modulation order betaiLAnd (3) calculating:
in channel 201, two branches beta of each channel are used iRAnd betaiLPerforming WFRFT inverse processing on the received signals S (n):
then, the two-branch inverse processing result S in each channeliR(n) and SiL(n) performing correlation operation to obtain correlation result C of each channeli(n), wherein fft (. cndot.) and ift (. cndot.) are fast Fourier transform and inverseThe processing function is transformed.
Further, use of Ci(n) calculating the Peak-to-average ratio PiWhere max (·) is a maximum calculation function.
As can be seen in fig. 1, within the main value window of α' ═ 0, with βiWhen scanning toIn the interval, two intervals with the maximum shunt correlation degree appear, and the interval presents rectangular distribution and is called as a rectangular correlation interval. The rectangular correlation interval range with respect to Δ β can be obtained by calculation, as shown in equations (15) and (16), Δ βRIs the right boundary of the rectangular correlation interval, Δ βLThe left border of the rectangular correlation interval.
As can be seen from fig. 1, as the main value window is shifted, the variation of the α 'value is equivalent, and no matter what the α' value is, within each shifted main value window, a rectangular correlation interval between Δ β and the peak value appears, and its right boundary Δ βRAnd left boundary Δ βLThe results of equations (15) and (16) are always satisfied.
Then, the peak-to-average ratio P of N channels is usediAnd comparing with a threshold G, and further calculating to obtain a continuous sequence R (j), wherein mod (-) is a modulus operation function.
And then, performing minimum value calculation on the continuous sequence R (j) by using a minimum value calculation function min (-) and performing minimum periodic processing to obtain a mapping channel H corresponding to the rectangular correlation interval:
H=mod(min[R(j)],N) (18)
calculating the demodulation order beta corresponding to the channel by using the mapping channel HH:
By using betaHAnd the right boundary Δ β of the rectangular intervalRAnd, in conjunction with equation (8), an estimated modulation order α' can be calculated:
wherein m is an integer, and m is greater than or equal to 0, and is selected so that alpha' is the minimum positive value within [04 ]. Since the real modulation order α is 4 as a period, the value of α is equivalent to be in the range of [04 ].
Then, WFRFT inverse processing is performed using the estimated α' and the received signal s (n):
finally, the WFRFT inverse data D' (n) is calculated as formula (22). It can be seen that if the modulation order α of the transmitting end is within an interval of 0 to 2, the demodulated data D '(n) is the real data x (n), and if the modulation order α of the transmitting end is within an interval of 2 to 4, the demodulated data D' (n) is the inverted data x (-n) of the real data.
Because it cannot be determined whether the original modulation order is within 0-2 or within 2-4, the data D ' (n) and the data D ' (-n) after D ' (n) inversion are output simultaneously, and the 2-path data analysis process is far more complicated than the 400-path data analysis process.
A high resolution blind scanning method for WFRFT signals carries out modulation order estimation and communication data reception under the condition of unknown modulation order (parameter), so that the same scanning conclusion can be reached no matter the modulation order is in a fixed form or a dynamic change form. In addition, since all the multiple-term WFRFT systems can be converted into corresponding 4-term WFRFT systems, a high-resolution blind scanning method for WFRFT signals is not only applicable to the 4-term WFRFT systems, but also applicable to the multiple-term WFRFT systems.
Claims (1)
1. A method for high resolution blind scanning of WFRFT signals is characterized in that: firstly, establishing a regular relation between a modulation order error and an error rate under a theoretical condition, and determining a minimum main value window; establishing a multi-channel demodulation order in a main value window, setting a left modulation order and a right modulation order in each channel, further performing WFRFT inverse processing of two branches on a received signal by each channel, performing correlation operation by using the inverse processing results of the two branches in each channel, and calculating a peak-to-average ratio; establishing a rectangular correlation interval by using the theoretical correlation degree, and obtaining a theoretical boundary value; calculating a mapping channel by using the peak-to-average ratios of a plurality of channels, and further calculating a demodulation order corresponding to the channel; calculating an estimated modulation order by using the calculation result; performing WFRFT inverse processing by using the estimated modulation order and the received signal, and obtaining demodulation data;
The method specifically comprises the following steps:
establishing a regular relation between modulation order error delta beta and an error rate under a theoretical condition, setting a minimum main value window K to be 2, and setting a stepping amount to be 0.01;
in a minimum principal value period, namely 0-2 intervals, carrying out multichannel scanning by taking lambda as a stepping quantity, wherein the scanning upper limit is N:
defining a demodulation order set betaiFormula (7), wherein i ∈ [1N ]],
βi=i·λ (7)
Assuming the estimated modulation order is α ', the error between the demodulation order of each channel and α' is:
Δβi=α'-βi (8)
will betaiIs subjected to periodic continuation to obtainThen, the reverse rotation and the translation are carried out to obtainAndin each channel betaiOn the basis, the right modulation order beta is carried outiRLeft modulation order betaiLAnd (3) calculating:
using two branches beta per channeliRAnd betaiLPerforming WFRFT inverse processing on the received signals S (n):
two-branch inverse processing result S in each channeliR(n) and SiL(n) performing correlation operation to obtain correlation result C of each channeli(n), wherein fft (-), ifft (-) are fast fourier transform and inverse transform processing functions;
by Ci(n) calculating the Peak-to-average ratio PiWhere max (·) is a maximum calculation function;
ΔβRis the right boundary of the rectangular correlation interval, Δ βLIs the left boundary of the rectangular correlation interval;
using peak-to-average ratio P of N channels iComparing with a threshold G, and further calculating to obtain a continuous sequence R (j), wherein mod (·) is a modular operation function;
and (3) performing minimum value calculation on the continuous sequence R (j) by using a minimum value calculation function min (·), and performing minimum periodic processing to obtain a mapping channel H corresponding to a rectangular correlation interval:
H=mod(min[R(j)],N) (18)
calculating the demodulation order beta corresponding to the channel by using the mapping channel HH:
By using betaHAnd the right boundary Δ β of the rectangular intervalRAnd combining with equation (8), the estimated modulation order α' can be calculated:
wherein m is an integer and is more than or equal to 0, and the alpha' is the minimum positive value in the range of [04 ] by selecting m;
performing WFRFT inverse processing by using the estimated alpha' and a receiving signal S (n):
calculating to obtain data D' (n) after WFRFT inverse processing as a formula (22); if the modulation order alpha of the transmitting end is within the range of 0-2, the demodulated data D '(n) is real data x (n), and if the modulation order alpha of the transmitting end is within the range of 2-4, the demodulated data D' (n) is inverted data x (-n) of the real data;
the data D ' (n) and the data D ' (-n) obtained by inverting the data D ' (n) are outputted simultaneously.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911263420.8A CN110933011B (en) | 2019-12-11 | 2019-12-11 | High-resolution blind scanning method for WFRFT (WFRFT) signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911263420.8A CN110933011B (en) | 2019-12-11 | 2019-12-11 | High-resolution blind scanning method for WFRFT (WFRFT) signals |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110933011A CN110933011A (en) | 2020-03-27 |
CN110933011B true CN110933011B (en) | 2022-06-10 |
Family
ID=69859723
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911263420.8A Active CN110933011B (en) | 2019-12-11 | 2019-12-11 | High-resolution blind scanning method for WFRFT (WFRFT) signals |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110933011B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103441822A (en) * | 2013-09-16 | 2013-12-11 | 哈尔滨工业大学 | Signal modulation method based on weighted-type fractional Fourier transform and energy spreading transform |
WO2017054632A1 (en) * | 2015-09-29 | 2017-04-06 | 华为技术有限公司 | Signal sending method, signal receiving method, transmitting terminal and receiving terminal |
CN107682297A (en) * | 2017-09-06 | 2018-02-09 | 西北工业大学 | A kind of mobile underwater sound communication method |
CN108920419A (en) * | 2018-06-20 | 2018-11-30 | 沈阳理工大学 | Scan method is classified based on one-parameter Fractional Fourier inverse transformation |
-
2019
- 2019-12-11 CN CN201911263420.8A patent/CN110933011B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103441822A (en) * | 2013-09-16 | 2013-12-11 | 哈尔滨工业大学 | Signal modulation method based on weighted-type fractional Fourier transform and energy spreading transform |
WO2017054632A1 (en) * | 2015-09-29 | 2017-04-06 | 华为技术有限公司 | Signal sending method, signal receiving method, transmitting terminal and receiving terminal |
CN107682297A (en) * | 2017-09-06 | 2018-02-09 | 西北工业大学 | A kind of mobile underwater sound communication method |
CN108920419A (en) * | 2018-06-20 | 2018-11-30 | 沈阳理工大学 | Scan method is classified based on one-parameter Fractional Fourier inverse transformation |
Also Published As
Publication number | Publication date |
---|---|
CN110933011A (en) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1221096C (en) | Channel estimation method for orthogonal frequency-division multiplexing communication system | |
US5790598A (en) | Block decision feedback equalizer | |
WO2006079181A1 (en) | Frequency estimation | |
CN1228951C (en) | Multicarrier transmission systems with simple channel response estimation | |
US6421401B1 (en) | Method and apparatus for achieving and maintaining symbol synchronization particularly in an OFDM system | |
CN1110148C (en) | Equalizing method and equalizer for OFDM receiver | |
CN107682297A (en) | A kind of mobile underwater sound communication method | |
CN1161121A (en) | Transmission system with improved symbol processing | |
CN106789813A (en) | The generation method of leading symbol | |
KR100186930B1 (en) | Initial frame synchronization method using ofdm method | |
CN1168576A (en) | Adaptive equalizer for controlling step size in proportion to estimated delay of received signals | |
CN1061186C (en) | Demodulator for manchester-coded FM signals | |
CN101141425A (en) | Time-division pilot based channel estimation method of mobile communication system | |
CN1420647A (en) | Receiver of quadrature frequency division multiplexing system | |
CN107666451B (en) | Channel estimation method for LTE system | |
CN1581740A (en) | Feedback type channel estimating method and device based on PN sequence and pilot frequency in OFDM system | |
CN1181673A (en) | Timing synchronization in receiver employing orthogonal frequency division multiplexing | |
CN108562883A (en) | A kind of maximum likelihood distance estimating algorithm of overloading wave radar system | |
CN111049766A (en) | Estimation method for PDSCH of 5G system | |
EP0913973A2 (en) | Method and system for joint timing recovery and channel estimation for DMT modems | |
CN1913509A (en) | Synchronous method in radio system of orthogonal frequency division multiplex technology | |
CN110933011B (en) | High-resolution blind scanning method for WFRFT (WFRFT) signals | |
CN1329786A (en) | Synchronization tracking method | |
CN1402919A (en) | Bit error estimates from pilot signals | |
CN112383492B (en) | Recursive compressed sensing method and system applied to short-wave OFDM double-selection sky wave channel estimation |
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 |