CN113556305B - FBMC iterative channel equalization method and system suitable for high-frequency selective fading - Google Patents

FBMC iterative channel equalization method and system suitable for high-frequency selective fading Download PDF

Info

Publication number
CN113556305B
CN113556305B CN202110682086.0A CN202110682086A CN113556305B CN 113556305 B CN113556305 B CN 113556305B CN 202110682086 A CN202110682086 A CN 202110682086A CN 113556305 B CN113556305 B CN 113556305B
Authority
CN
China
Prior art keywords
symbol
equalization
fbmc
oqam
symbols
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110682086.0A
Other languages
Chinese (zh)
Other versions
CN113556305A (en
Inventor
李靖
吕宝均
李慧芳
任德锋
葛建华
闫伟平
武思同
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Zexin Technology Co.,Ltd.
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN202110682086.0A priority Critical patent/CN113556305B/en
Publication of CN113556305A publication Critical patent/CN113556305A/en
Application granted granted Critical
Publication of CN113556305B publication Critical patent/CN113556305B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2689Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2695Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3911Fading models or fading generators
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3912Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03878Line equalisers; line build-out devices
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2649Demodulators
    • H04L27/26534Pulse-shaped multi-carrier, i.e. not using rectangular window
    • H04L27/2654Filtering per subcarrier, e.g. filterbank multicarrier [FBMC]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2697Multicarrier modulation systems in combination with other modulation techniques
    • H04L27/2698Multicarrier modulation systems in combination with other modulation techniques double density OFDM/OQAM system, e.g. OFDM/OQAM-IOTA system
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • H04L2027/0024Carrier regulation at the receiver end
    • H04L2027/0026Correction of carrier offset
    • H04L2027/0038Correction of carrier offset using an equaliser

Abstract

The invention belongs to the technical field of multi-carrier communication, and discloses an FBMC iterative channel equalization method and system suitable for high-frequency selective fading, wherein the FBMC iterative channel equalization method suitable for high-frequency selective fading comprises the following steps: generating a transmission signal and a reception signal, respectively; calculating the corresponding symbol a p,q The analysis filter bank output signal of (a); initializing an iterative equalization algorithm, wherein q =2; carrying out first equalization; reconstructing the symbol for the first time; carrying out second equalization; reconstructing the symbol for the second time; let q = q +1, repeat the first equalization, the first reconstructed symbol, the second equalization and the second reconstructed symbol if q < N-1; otherwise, all equalization information is obtained, and the channel equalization is finished. The analysis filter output signal model used by the invention is an accurate model, does not use strict slow fading or flat fading approximation in the traditional method, is suitable for any fading FBMC/OQAM communication system, eliminates residual interference under a high-frequency selective channel, and improves the system performance.

Description

FBMC iterative channel equalization method and system suitable for high-frequency selective fading
Technical Field
The invention belongs to the technical field of multi-carrier communication, and particularly relates to an FBMC iterative channel equalization method and system suitable for high-frequency selective fading.
Background
At present, compared with an Orthogonal Frequency Division Multiplexing (OFDM) technique, a Filter Bank Multi-Carrier (FBMC) technique adopts a prototype Filter with good time-Frequency localization characteristics, so that the Filter Bank FBMC technique can well resist interference without a cyclic prefix, and has the advantages of high spectral efficiency, low out-of-band leakage and the like, and therefore receives wide attention. The Offset Quadrature Amplitude Modulation (OQAM) technique adopted by the FBMC system can only satisfy the orthogonality of the real number domain between the subcarriers, and there is inherent imaginary part interference in the complex number domain, so how to reduce the equalization error caused by the inherent interference becomes a hot point of research in recent years.
At present, the research method for channel equalization of FBMC/OQAM system mainly still uses Zero Forcing (ZF) equalization and Minimum Mean Square Error (MMSE) equalization methods commonly used in OFDM systems, wherein Baltar LG and Waldhauser D S, etc. design a fractional-interval Decision Feedback Equalizer (DFE) in "MMSE sub-channel decision feedback equalization for filter bank based multicarrier systems", which can reduce mean square error (MMSE) and has better equalization performance, but this scheme can only be used in flat fading, and once the frequency selectivity is improved, the equalization performance of this scheme will be greatly reduced. In recent years, FBMC/OQAM channel estimation and equalization techniques combined with Deep learning have become a research hotspot, wherein x.cheng and d.liu et al in "Deep learning-based channel estimation and equalization scheme for FBMC/OQAM systems" and "arescan-DNN based channel estimation and equalization scheme in FBMC/OQAM systems" have performed Deep learning on channel estimation and equalization schemes, wherein channel state information and constellation demapping methods are learned through a Deep neural network model, and then a distorted frequency domain sequence is implicitly equalized to directly obtain binary bits. Although equalization techniques have been extensively studied in recent years, most of them are still studied under severe slow fading or flat fading channels, and in high frequency selective channels, the traditional equalization scheme will cause residual interference to remain and affect the system performance. Therefore, an FBMC iterative channel equalization method suitable for high frequency selective fading is needed.
Through the above analysis, the problems and defects of the prior art are as follows:
(1) The existing equalization method has the problem of poor equalization performance under high-frequency selective fading.
(2) Most of the existing equalization technologies are studied under strict slow fading or flat fading channels, and the traditional equalization scheme can cause residual interference to remain and affect the system performance under high-frequency selective channels.
The difficulty in solving the above problems and defects is:
the existing equalization scheme has insufficient accuracy on interference calculation during equalization under a high-frequency selective fading channel, and the performance of the system is seriously influenced.
The significance of solving the problems and the defects is as follows:
through symbol reconstruction after equalization, pseudo pilot frequency calculation and iterative equalization, inherent interference can be greatly inhibited, and the equalization performance can be effectively improved.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides an FBMC iterative channel equalization method and system suitable for high-frequency selective fading, and aims to solve the problem that the existing equalization method has poor equalization performance under high-frequency selective fading.
The invention is realized in such a way that an FBMC iterative channel equalization method suitable for high-frequency selective fading comprises the following steps:
step one, generating a sending signal and a receiving signal respectively;
step two, calculating corresponding symbol a p,q The analysis filter bank output signal of (1);
initializing an iterative equalization algorithm, wherein q =2;
step four, carrying out first equalization;
step five, reconstructing the symbol for the first time;
step six, performing second equalization;
step seven, reconstructing the symbol for the second time;
step eight, enabling q = q +1, and if q is smaller than N-1, repeating the step four to the step seven; otherwise, all equalization information is obtained, and the channel equalization is finished.
Further, in the first step, the generating the transmission signal and the reception signal respectively includes:
(1) Constructing pilot symbols a composed of FBMC/OQAM symbols 0 All zero symbols a of FBMC/OQAM 1 FBMC/OQAM, data symbol a 2 ,a 3 ,…,a N-1 A frame structure of the composition; wherein N is more than or equal to 2 and is a positive integer, and represents the total number of FBMC/OQAM symbols in a frame; a is n =[a 0,n ,a 1,n ,…,a M-1,n ] T Wherein M =2 z Z is more than or equal to 2, z is a positive integer, and M represents the number of subcarriers in one FBMC/OQAM symbol;
(2) Calculating a baseband transmission signal s (l) consisting of M subcarriers and N FBMC/OQAM symbols:
Figure BDA0003123171430000031
wherein l =0,1, \8230, (N-1) M/2+ L g -1,
Figure BDA0003123171430000032
g m,n (l) Denotes a m,n G (i) adopts a PHYDYAS prototype filter with the length of L g KM +1, k representing the overlap factor, take a typical value of 4;
Figure BDA0003123171430000033
Figure BDA0003123171430000034
G 0 =1,G 1 =0.97196,
Figure BDA0003123171430000035
G 3 =0.235147;
(3) And (2) passing the baseband transmission signal s (l) through a multipath fading channel to obtain a baseband receiving signal y (l):
Figure BDA0003123171430000036
wherein [ h (0, n), h (1, n), \ 8230;, h (L) h -1,n)]Denotes the firstImpulse response of multipath channel during transmission of n FBMC/OQAM symbols, L h Denotes maximum, w (l) denotes mean zero and variance σ 2 Complex white gaussian noise.
Further, in step (1), the constructing of the pilot and data frame structure includes:
1) The first FBMC/OQAM symbol is a pilot symbol a 0 =[a 0,0 ,a 1,0 ,…,a M-1,0 ] T Wherein a is m,0 Represents a pseudo-randomly generated real number with a value of 1 or-1,m =0,1, \8230;, M-1;
2) The second FBMC/OQAM symbol is pilot symbol a 1 =[0,0,…,0] T For isolating pilot symbols from data symbols;
3) The remaining N-2 FBMC/OQAM symbols a 2 ,a 3 ,…,a N-1 For transmitting data, wherein a m,n Is a real-valued OQAM symbol, M =0,1, \ 8230;, M-1, N =4,1, \ 8230;, N-1.
Further, in step two, the calculation of the corresponding symbol a p,q Comprises:
output signal y of symbol obtained after baseband receiving signal passes through analysis filter bank p,q Comprises the following steps:
Figure BDA0003123171430000041
wherein, time frequency points (m, n) and (p, q) respectively represent OQAM symbols positioned on the nth (q) time symbol and the mth (p) subcarrier; Δ m = m-p, Δ n = n-q representing the difference in indices; h (Δm,Δn) (p + Δ m, q + Δ n) belongs to
Figure BDA0003123171430000042
A channel frequency domain response vector H corresponding to M sub-carriers on the nth FBMC/OQAM symbol n =[H(0,n),H(1,n),...,H(M-1,n)] T And interference factor moment E (Δm,Δn) Multiplication of the arrays:
Figure BDA0003123171430000043
further, in step three, the iterative equalization algorithm is initialized, and q =2, including:
symbol index q =2 for initial equalization, symbol after first equalization
Figure BDA0003123171430000044
Symbol after second equalization
Figure BDA0003123171430000045
Further, in step four, the first equalizing includes:
considering the first-order left-adjacent domain of a symbol
Figure BDA0003123171430000046
Interference analysis filterbank output signal y using channel response matrix and q-th symbol q Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q First time estimate of
Figure BDA0003123171430000051
(1) The inherent interference of the FBMC/OQAM system mainly comes from first-order interference, and according to the sequence of equalization on time symbols, the first equalization only considers the inherent interference brought by a first-order left neighborhood, and is expressed as:
Figure BDA0003123171430000052
(2) Let y q =[y 0,q ,y 1,q ,...,y M-1,q ] T The output signal y of the analysis filterbank for the qth symbol vector q Expressed as:
Figure BDA0003123171430000053
(3) Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q First time estimate of
Figure BDA0003123171430000054
Wherein the expression of the zero-forcing equalization is:
Figure BDA0003123171430000055
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003123171430000056
the operation of taking the real part of the pair.
Further, in step five, the reconstructing the symbol for the first time includes:
demodulated data symbol estimates
Figure BDA0003123171430000057
Recovering the bit data stream of the transmitting end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q First OQAM reconstructed symbol vector
Figure BDA0003123171430000058
Further, in the sixth step, the second equalizing includes:
considering a first order neighborhood of symbols
Figure BDA0003123171430000059
Interference analysis filterbank output signal y using channel response matrix and q-1 symbol q- 1, obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q-1 Second order estimate of
Figure BDA00031231714300000510
(1) According to the data symbol after the first equalization reconstruction, the inherent interference brought by the first-order neighborhood is considered:
Figure BDA00031231714300000511
(2) Analysis filterbank output signal y for the q-1 th symbol q-1 Expressed as:
Figure BDA0003123171430000061
(3) Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q-1 Second order estimate of
Figure BDA0003123171430000062
Wherein the expression of the zero-forcing equalization is:
Figure BDA0003123171430000063
further, in step seven, the reconstructing the symbol for the second time includes:
demodulated data symbol estimates
Figure BDA0003123171430000064
Recovering the bit data stream of the sending end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q-1 Second OQAM reconstruction of the symbol vector
Figure BDA0003123171430000065
Another object of the present invention is to provide an FBMC iterative channel equalization system suitable for high frequency selective fading, which applies the FBMC iterative channel equalization method suitable for high frequency selective fading, and the FBMC iterative channel equalization system suitable for high frequency selective fading comprises:
a signal generating module for generating a transmitting signal s (l) and a receiving signal y (l) respectively;
an output signal calculation module for calculating the corresponding symbol a p,q Analysis filter bank output signal y p,q
The initialization module is used for initializing an iterative equalization algorithm, and q =2;
a first-order equalization module for considering a first-order left-adjacent domain of symbols
Figure BDA0003123171430000066
Interference analysis filterbank output signal y using channel response matrix and qth symbol q Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q First time estimate of
Figure BDA0003123171430000067
A first symbol reconstruction module for demodulating the data symbol estimation value
Figure BDA0003123171430000068
Recovering the bit data stream of the transmitting end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q First OQAM reconstruction of symbol vectors
Figure BDA0003123171430000069
A second equalization module for considering a first order neighborhood of symbols
Figure BDA0003123171430000071
Interference analysis filterbank output signal y using channel response matrix and q-1 symbol q-1 Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q-1 Second order estimate of
Figure BDA0003123171430000072
A second symbol reconstruction module for demodulating the data symbol estimate
Figure BDA0003123171430000073
Recovering the bit data stream of the transmitting end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q-1 Second OQAM reconstructing the symbol vector
Figure BDA0003123171430000074
The judging module is used for enabling q = q +1, and if q is smaller than N-1, the first equalizing module, the first symbol reconstructing module, the second equalizing module and the second symbol reconstructing module are repeatedly carried out; otherwise, all equalization information is obtained, and the channel equalization is finished.
Another object of the present invention is to provide an information data processing terminal for implementing the FBMC iterative channel equalization method suitable for high-frequency selective fading.
By combining all the technical schemes, the invention has the advantages and positive effects that: the FBMC iterative channel equalization method suitable for high-frequency selective fading provided by the invention utilizes the symbols reconstructed twice as known data information to calculate first-order neighborhood interference, can eliminate residual interference under a high-frequency selective channel, obviously improves the equalization performance of a system, and can be used for a channel equalizer in an FBMC communication system.
The analysis filter output signal model used by the invention is an accurate model, does not use strict slow fading or flat fading approximation in the traditional method, and can be suitable for any fading FBMC/OQAM communication system.
The method takes the neighborhood interference into account during equalization, eliminates the inherent interference among data which cannot be considered by channel estimation, and further improves the equalization performance of the system through iteration.
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 it is obvious that the drawings described below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a flowchart of an FBMC iterative channel equalization method suitable for high frequency selective fading according to an embodiment of the present invention.
Fig. 2 is a schematic diagram of an FBMC iterative channel equalization method suitable for high frequency selective fading according to an embodiment of the present invention.
Fig. 3 is a block diagram of an FBMC iterative channel equalization system suitable for high frequency selective fading according to an embodiment of the present invention;
in the figure: 1. a signal generation module; 2. an output signal calculation module; 3. initializing a module; 4. a first equalization module; 5. a first symbol reconstruction module; 6. a second equalization module; 7. a second symbol reconstruction module; 8. and a judging module.
Fig. 4 is a frame structure diagram according to an embodiment of the present invention.
Fig. 5 is a comparison graph of system error performance of the method of the present invention and the conventional method according to the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and do not limit the invention.
In view of the problems in the prior art, the present invention provides an FBMC iterative channel equalization method and system suitable for high frequency selective fading, and the following describes the present invention in detail with reference to the accompanying drawings.
As shown in fig. 1, an FBMC iterative channel equalization method suitable for high frequency selective fading provided by the embodiment of the present invention includes the following steps:
s101, respectively generating a sending signal and a receiving signal;
s102, calculating corresponding symbol a p,q The analysis filter bank output signal of (1);
s103, initializing an iterative equalization algorithm, wherein q =2;
s104, carrying out first equalization;
s105, reconstructing the symbol for the first time;
s106, performing second equalization;
s107, reconstructing the symbol for the second time;
s108, making q = q +1, if q is less than N-1, repeating S104-S107; otherwise, all equalization information is obtained, and the channel equalization is finished.
A schematic diagram of an FBMC iterative channel equalization method suitable for high frequency selective fading according to an embodiment of the present invention is shown in fig. 2.
As shown in fig. 3, an FBMC iterative channel equalization system suitable for high frequency selective fading according to an embodiment of the present invention includes:
a signal generating module 1, configured to generate a transmission signal s (l) and a reception signal y (l), respectively;
an output signal calculation module 2 for calculating the corresponding symbol a p,q Analysis filter bank output signal y p,q
The initialization module 3 is used for initializing an iterative equalization algorithm, and q =2;
a first equalization block 4 for considering a first order left neighborhood of the symbols
Figure BDA0003123171430000091
Interference analysis filterbank output signal y using channel response matrix and q-th symbol q Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q First time estimate of
Figure BDA0003123171430000092
A first symbol reconstruction module 5 for demodulating the data symbol estimation value
Figure BDA0003123171430000093
Recovering the bit data stream of the sending end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q First OQAM reconstructed symbol vector
Figure BDA0003123171430000094
A second equalization module 6 for considering the first order neighbourhood of the symbol
Figure BDA0003123171430000095
Interference analysis filterbank output signal y using channel response matrix and q-1 symbol q-1 Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q-1 Second order estimate of
Figure BDA0003123171430000096
A second symbol reconstruction module 7 for demodulating the data symbol estimation value
Figure BDA0003123171430000097
Recovering the bit data stream of the transmitting end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q-1 Second OQAM reconstructing the symbol vector
Figure BDA0003123171430000098
A judging module 8, configured to make q = q +1, and if q is less than N-1, repeat the first equalizing module, the first symbol reconstructing module, the second equalizing module, and the second symbol reconstructing module; otherwise, all equalization information is obtained, and the channel equalization is finished.
The technical solution of the present invention will be further described with reference to the following examples.
The invention discloses an FBMC iterative channel equalization method suitable for high-frequency selective fading, which mainly solves the problem that the existing equalization method is poor in equalization performance under high-frequency selective fading. The technical scheme comprises the following steps: generating a transmission signal; generating a received signal; calculating an analysis filter bank output signal of the corresponding symbol; balancing for the first time; reconstructing the symbol for the first time; reconstructing the symbol for the second time; second equalization; and judging whether the equalization of the last symbol is reached, if not, starting the equalization of the next symbol, otherwise, obtaining the equalization signals of all the symbols, and ending the iterative equalization. The invention uses the symbols reconstructed twice as the known data information to calculate the first-order neighborhood interference, obviously improves the equalization performance and can be used for a channel equalizer in an FBMC communication system.
The implementation of the invention is further described with reference to fig. 2.
Step 1, generating a transmission signal
1a) Constructing pilot symbols a composed of FBMC/OQAM symbols 0 All zero symbols a of FBMC/OQAM 1 FBMC/OQAM data symbols a 2 ,a 3 ,…,a N-1 The frame structure is formed, wherein N is more than or equal to 2 and is a positive integer, and the N represents the total number of FBMC/OQAM symbols in a frame. a is a n =[a 0,n ,a 1,n ,…,a M-1,n ] T Wherein M =2 z Z is more than or equal to 2, z is a positive integer, and M represents the number of subcarriers in one FBMC/OQAM symbol.
Referring to fig. 4, the specific structure of the pilot and data frame structure is:
first, the first FBMC/OQAM symbol is a pilot symbol a 0 =[a 0,0 ,a 1,0 ,…,a M-1,0 ] T Wherein a is m,0 Representing a pseudo-randomly generated real number with a value of 1 or-1,m = -0, 1, \8230, M-1;
secondly, the second FBMC/OQAM symbol is pilot symbol a 1 =[0,0,…,0] T The method is used for isolating the pilot frequency symbol and the data symbol, so that the inherent interference of the data to the pilot frequency symbol can be reduced, and the estimation performance is improved;
finally, N-2 FBMC/OQAM symbols a remain 2 ,a 3 ,…,a N-1 For transmitting data, wherein a m,n Is a real-valued OQAM symbol, M =0,1, \ 8230;, M-1, N =4,1, \ 8230;, N-1.
Calculating a baseband transmission signal s (l) consisting of M subcarriers and N FBMC/OQAM symbols:
Figure BDA0003123171430000111
wherein l =0,1, \8230; (N-1) M/2+ L g -1,
Figure BDA0003123171430000112
g m,n (l) Denotes a m,n G (i) filtering with a PHYDYAS prototypeDevice with length L g K = KM +1, representing an overlap factor, inappropriately general, K being a typical value of 4;
Figure BDA0003123171430000113
Figure BDA0003123171430000114
G0=1,G1=0.97196,
Figure BDA0003123171430000115
G3=0.235147。
step 2, generating a receiving signal
And (2) passing the baseband transmission signal s (l) generated in the step (1) through a multipath fading channel to obtain a baseband receiving signal:
Figure BDA0003123171430000116
wherein [ h (0, n), h (1, n), \ 8230;, h (L) h -1,n)]Represents the impulse response of the multipath channel during the transmission of the nth FBMC/OQAM symbol, L h Denotes maximum, w (l) denotes mean zero and variance σ 2 Complex white gaussian noise.
Step 3, calculating corresponding symbol a p,q Analysis filter bank output signal
Output signal y of symbol obtained after baseband receiving signal passes through analysis filter bank p,q Comprises the following steps:
Figure BDA0003123171430000117
wherein, time-frequency points (m, n) and (p, q) respectively represent OQAM symbols positioned on the nth (q) time symbol and the mth (p) subcarrier, Δ m = m-p, Δ n = n-q represent difference of indexes, H (Δm,Δn) (p + Δ m, q + Δ n) belongs to
Figure BDA0003123171430000118
The channel frequency domain response vector corresponding to M sub-carriers on the nth FBMC/OQAM symbolH n =[H(0,n),H(1,n),...,H(M-1,n)] T And interference factor moment E (Δm,Δn) Multiplication of the arrays:
Figure BDA0003123171430000121
step 4, initializing iterative equalization algorithm
Symbol index q =2 for initial equalization, symbol after first equalization
Figure BDA0003123171430000122
Symbol after second equalization
Figure BDA0003123171430000123
Step 5, first equalization
(5a) The inherent interference of the FBMC/OQAM system mainly comes from first-order interference, and according to the sequence of equalization on the time symbol, the first equalization can only consider the inherent interference brought by the first-order left neighborhood, which can be expressed as:
Figure BDA0003123171430000124
(5b) Let y q =[y 0,q ,y 1,q ,...,y M-1,q ] T Then the output signal of the analysis filter for the qth symbol vector is:
Figure BDA0003123171430000125
Figure BDA0003123171430000126
(5c) Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q First time estimate of
Figure BDA0003123171430000127
Figure BDA0003123171430000128
Wherein the content of the first and second substances,
Figure BDA0003123171430000129
the operation of taking the real part is carried out.
Step 6. First reconstructing the symbol
Demodulated data symbol estimates
Figure BDA00031231714300001210
Recovering the bit data stream of the sending end; then, OQAM modulation is carried out on the recovered bit data stream to obtain a corresponding symbol vector a q First OQAM reconstruction of symbol vectors
Figure BDA00031231714300001211
Step 7. Second equalization
(7a) According to the data symbol after the first equalization reconstruction, the inherent interference caused by the first-order neighborhood is considered, which can be expressed as:
Figure BDA0003123171430000131
(7b) The analysis filter output signal for the q-1 th symbol vector is then:
Figure BDA0003123171430000132
Figure BDA0003123171430000133
Figure BDA0003123171430000134
(7c) Obtaining symbols by adopting a zero forcing equalization algorithmVector a q-1 Second order estimate of (2)
Figure BDA0003123171430000135
Figure BDA0003123171430000136
Step 8, reconstructing the symbol for the second time
Demodulated data symbol estimates
Figure BDA0003123171430000137
Recovering the bit data stream of the sending end; then, OQAM modulation is carried out on the recovered bit data stream to obtain a corresponding symbol vector a q- Second OQAM reconstructed symbol vector of 1
Figure BDA0003123171430000138
Step 9. Let q = q +1, if q < N-1, repeat steps (5) - (8); otherwise, equalization information is obtained, and channel equalization is finished.
The technical effects of the present invention will be described in detail with reference to simulation experiments.
1. Simulation conditions
The effect of the present invention is illustrated by simulation in two frequency selective fading scenarios, and two frequency selective fading scenario simulation parameters are shown in table 1.
TABLE 1 simulation parameters
Figure BDA0003123171430000139
Figure BDA0003123171430000141
2. Simulation content and results
The performance of the channel estimation normalized mean square error NMSE of the method and the traditional AFB output signal model, which changes along with the signal-to-noise ratio SNR, is simulated and compared, and the simulation result is shown in figure 5.
As can be seen from fig. 5: 1) Under two frequency selective fading scenes, the NMSE of the channel estimation provided by the invention has better performance than that of the traditional method; 2) The method provided by the invention has the advantages that the performance is obviously improved in a high-frequency selectivity scene, and particularly, the gain of 6-9 dB can be obtained under a high signal-to-noise ratio.
In the above embodiments, the implementation may be wholly or partially realized by software, hardware, firmware, or any combination thereof. When used in whole or in part, can be implemented in a computer program product that includes one or more computer instructions. When loaded or executed on a computer, cause the flow or functions according to embodiments of the invention to occur, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a network of computers, or other programmable device. The computer instructions may be stored in a computer readable storage medium or transmitted from one computer readable storage medium to another computer readable storage medium, for example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wire (e.g., coaxial cable, fiber optic, digital Subscriber Line (DSL), or wireless (e.g., infrared, wireless, microwave, etc.)). The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device, such as a server, a data center, etc., that includes one or more of the available media. The usable medium may be a magnetic medium (e.g., floppy Disk, hard Disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., solid State Disk (SSD)), among others.
The above description is only for the purpose of illustrating the embodiments of the present invention, and the scope of the present invention should not be limited thereto, and any modifications, equivalents and improvements made by those skilled in the art within the technical scope of the present invention as disclosed in the present invention should be covered by the scope of the present invention.

Claims (4)

1. An FBMC iterative channel equalization method suitable for high-frequency selective fading, which is characterized by comprising the following steps:
step one, generating a sending signal and a receiving signal respectively;
step two, calculating corresponding symbol a p,q The analysis filter bank output signal of (a);
initializing an iterative equalization algorithm, wherein q =2;
step four, carrying out first equalization;
step five, reconstructing the symbol for the first time;
step six, carrying out second equalization;
step seven, reconstructing the symbol for the second time;
step eight, enabling q = q +1, and if q is smaller than N-1, repeating the step four to the step seven; otherwise, all equalization information is obtained, and the channel equalization is finished;
in the first step, the generating the transmission signal and the reception signal respectively includes:
(1) Constructing pilot symbols a composed of FBMC/OQAM symbols 0 All zero symbols a of FBMC/OQAM 1 FBMC/OQAM, data symbol a 2 ,a 3 ,θ,a N-1 A frame structure of components; wherein N is more than or equal to 2 and is a positive integer, and represents the total number of FBMC/OQAM symbols in a frame; a is a n =[a 0,n ,a 1,n ,...,a M-1,n ] T Wherein M =2 z Z is more than or equal to 2, z is a positive integer, and M represents the number of subcarriers in one FBMC/OQAM symbol;
(2) Calculating a baseband transmission signal s (l) consisting of M subcarriers and N FBMC/OQAM symbols:
Figure FDA0003879800370000011
wherein l =0,1, \8230, (N-1) M/2+ L g -1,
Figure FDA0003879800370000012
g m,n (l) Denotes a m,n G (i) adopts a PHYDYAS prototype filter with the length of L g K = KM +1, representing an overlap factor, with a typical value of 4;
Figure FDA0003879800370000021
Figure FDA0003879800370000022
G 0 =1,G 1 =0.97196,
Figure FDA0003879800370000023
G 3 =0.235147;
(3) Passing the baseband transmission signal s (l) through a multipath fading channel to obtain a baseband receiving signal y (l):
Figure FDA0003879800370000024
wherein, h (0, n), h (1, n), \8230;, h (L) h -1,n)]Represents the impulse response of the multipath channel during the transmission of the nth FBMC/OQAM symbol, L h Denotes maximum, w (l) denotes mean zero and variance σ 2 Complex white gaussian noise of (a);
in step (1), the constructing of the pilot frequency and data frame structure includes:
1) The first FBMC/OQAM symbol is a pilot symbol a 0 =[a 0,0 ,a 1,0 ,…,a M-1,0 ] T Wherein a is m,0 Representing a pseudo-randomly generated real number with a value of 1 or-1,m = -0, 1, \8230, M-1;
2) The second FBMC/OQAM symbol is pilot symbol a 1 =[0,0,...,0] T For isolating pilot symbols and data symbols;
3) The remaining N-2 FBMC/OQAM symbols a 2 ,a 3 ,…,a N-1 For transmitting data, wherein a m,n Is a real-valued OQAM symbol, M =0,1, \8230, M-1, N =4,1, \8230, N-1;
in step two, the corresponding symbol a is calculated p,q Comprises:
output signal y of symbol obtained after baseband receiving signal passes through analysis filter bank p,q Comprises the following steps:
Figure FDA0003879800370000025
wherein, time frequency points (m, n) and (p, q) respectively represent OQAM symbols positioned on the nth (q) time symbol and the mth (p) subcarrier; Δ m = m-p, Δ n = n-q representing the difference in indices; h (Δm,Δn) (p + Δ m, q + Δ n) belongs to
Figure FDA0003879800370000026
A channel frequency domain response vector H corresponding to M sub-carriers on the nth FBMC/OQAM symbol n =[H(0,n),H(1,n),...,H(M-1,n)] T And interference factor moment E (Δm,Δn) Product of the arrays:
Figure FDA0003879800370000031
in step four, the first equalizing includes:
considering the first order left neighbourhood of the symbol
Figure FDA0003879800370000032
Interference analysis filterbank output signal y using channel response matrix and q-th symbol q Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q First time estimate of
Figure FDA0003879800370000033
(1) The inherent interference of the FBMC/OQAM system mainly comes from first-order interference, and according to the sequence of equalization on time symbols, the first equalization only considers the inherent interference brought by a first-order left neighborhood, and is expressed as:
Figure FDA0003879800370000034
(2) Let y q =[y 0,q ,y 1,q ,...,y M-1,q ] T The output signal y of the analysis filterbank for the qth symbol vector q Expressed as:
Figure FDA0003879800370000035
(3) Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q First time estimate of
Figure FDA0003879800370000036
Wherein the expression of the zero-forcing equalization is as follows:
Figure FDA0003879800370000037
wherein the content of the first and second substances,
Figure FDA0003879800370000038
the operation of taking the real part is carried out;
in the fifth step, the first symbol reconstruction includes:
demodulated data symbol estimates
Figure FDA0003879800370000039
Recovering the bit data stream of the sending end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q First OQAM reconstructed symbol vector
Figure FDA00038798003700000310
In step six, the second equalization includes:
considering a first order neighborhood of symbols
Figure FDA00038798003700000311
Interference analysis filterbank output signal y using channel response matrix and q-1 symbol q-1 Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q-1 Second order estimate of
Figure FDA0003879800370000041
(1) According to the data symbol after the first equalization reconstruction, the inherent interference brought by the first-order neighborhood is considered:
Figure FDA0003879800370000042
(2) Analysis filterbank output signal y for the q-1 th symbol q-1 Expressed as:
Figure FDA0003879800370000043
(3) Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q-1 Second order estimate of (2)
Figure FDA0003879800370000044
Wherein the expression of the zero-forcing equalization is:
Figure FDA0003879800370000045
in step seven, the second reconstructing the symbol includes:
demodulated data symbol estimates
Figure FDA0003879800370000046
Recovering the bit data stream of the sending end; oqam modulation of recovered bit streamTo obtain the corresponding symbol vector a q-1 Second OQAM reconstruction of the symbol vector
Figure FDA0003879800370000047
2. The FBMC iterative channel equalization method for high frequency selective fading as claimed in claim 1, wherein in step three, said iterative equalization algorithm is initialized, and q =2, comprising:
symbol index q =2 for initial equalization, symbol after first equalization
Figure FDA0003879800370000048
Symbols after second equalization
Figure FDA0003879800370000049
3. An FBMC iterative channel equalization system for high frequency selective fading, which implements the FBMC iterative channel equalization method for high frequency selective fading according to any one of claims 1 to 2, wherein the FBMC iterative channel equalization system for high frequency selective fading comprises:
a signal generating module for generating a transmitting signal s (l) and a receiving signal y (l) respectively;
an output signal calculation module for calculating a corresponding symbol a p,q Analysis filter bank output signal y p,q
The initialization module is used for initializing an iterative equalization algorithm, and q =2;
a first-order equalization module for considering a first-order left-adjacent domain of symbols
Figure FDA0003879800370000051
Interference analysis filterbank output signal y using channel response matrix and q-th symbol q Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q First time estimate of
Figure FDA0003879800370000052
A first time symbol reconstruction module for demodulating the data symbol estimation value
Figure FDA0003879800370000053
Recovering the bit data stream of the transmitting end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q First OQAM reconstructed symbol vector
Figure FDA0003879800370000054
A second equalization module for considering a first order neighborhood of symbols
Figure FDA0003879800370000055
Interference analysis filterbank output signal y using channel response matrix and q-1 symbol q-1 Obtaining a symbol vector a by adopting a zero-forcing equalization algorithm q-1 Second order estimate of
Figure FDA0003879800370000056
A second symbol reconstruction module for demodulating the data symbol estimate
Figure FDA0003879800370000057
Recovering the bit data stream of the sending end; carrying out OQAM modulation on the recovered bit data stream to obtain a corresponding symbol vector a q-1 Second OQAM reconstructing the symbol vector
Figure FDA0003879800370000058
The judging module is used for enabling q = q +1, and if q is smaller than N-1, the first equalizing module, the first symbol reconstructing module, the second equalizing module and the second symbol reconstructing module are repeatedly carried out; otherwise, all equalization information is obtained, and the channel equalization is finished.
4. An information data processing terminal, characterized in that the information data processing terminal is configured to implement the FBMC iterative channel equalization method suitable for high-frequency selective fading according to any one of claims 1 to 2.
CN202110682086.0A 2021-06-19 2021-06-19 FBMC iterative channel equalization method and system suitable for high-frequency selective fading Active CN113556305B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110682086.0A CN113556305B (en) 2021-06-19 2021-06-19 FBMC iterative channel equalization method and system suitable for high-frequency selective fading

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110682086.0A CN113556305B (en) 2021-06-19 2021-06-19 FBMC iterative channel equalization method and system suitable for high-frequency selective fading

Publications (2)

Publication Number Publication Date
CN113556305A CN113556305A (en) 2021-10-26
CN113556305B true CN113556305B (en) 2022-12-13

Family

ID=78130757

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110682086.0A Active CN113556305B (en) 2021-06-19 2021-06-19 FBMC iterative channel equalization method and system suitable for high-frequency selective fading

Country Status (1)

Country Link
CN (1) CN113556305B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116192576B (en) * 2023-03-02 2023-08-11 哈尔滨工程大学 PSRO prototype filter for FBMC system against frequency selective fading

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016117980A1 (en) * 2015-01-23 2016-07-28 Samsung Electronics Co., Ltd. Method and auxiliary method, apparatus, base station and user equipment for interference cancellation
CN109672637A (en) * 2019-02-28 2019-04-23 重庆邮电大学 A kind of FBMC channel estimation methods based on interferential loads iteration
CN112491756A (en) * 2020-11-23 2021-03-12 西安电子科技大学 RLS channel estimation method for joint channel equalization in FBMC system

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2782304B1 (en) * 2013-03-19 2016-06-22 Fundació Centre Tecnologic de Telecomunicacions de Catalunya Method for equalizing filterbank multicarrier (fbmc) modulations
EP2819476A1 (en) * 2013-06-27 2014-12-31 Alcatel Lucent Resource allocation with low signalling overhead
EP2840749B1 (en) * 2013-08-23 2020-09-30 Alcatel Lucent Receiver and receive method for a filtered multicarrier signal
FR3013928B1 (en) * 2013-11-28 2015-12-25 Commissariat Energie Atomique CHANNEL ESTIMATION METHOD FOR FBMC TELECOMMUNICATION SYSTEM
WO2015154799A1 (en) * 2014-04-08 2015-10-15 Huawei Technologies Duesseldorf Gmbh Transmitter, receiver and system for filterbank multicarrier communication
EP3076576A1 (en) * 2015-04-01 2016-10-05 NTT DoCoMo, Inc. Transmit diversity from orthogonal design for fbmc/oqam
CN106302279A (en) * 2016-08-04 2017-01-04 成都极比特通信技术有限公司 FBMC system equalization method based on interference variance statistics
CN106788935B (en) * 2016-12-23 2020-01-10 西安交通大学 Pilot frequency design and channel estimation method for high frequency spectrum efficiency in FBMC system
CN107438038B (en) * 2017-06-07 2020-03-13 西安交通大学 Pilot design and synchronous channel estimation method of FBMC/OQAM

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016117980A1 (en) * 2015-01-23 2016-07-28 Samsung Electronics Co., Ltd. Method and auxiliary method, apparatus, base station and user equipment for interference cancellation
CN109672637A (en) * 2019-02-28 2019-04-23 重庆邮电大学 A kind of FBMC channel estimation methods based on interferential loads iteration
CN112491756A (en) * 2020-11-23 2021-03-12 西安电子科技大学 RLS channel estimation method for joint channel equalization in FBMC system

Also Published As

Publication number Publication date
CN113556305A (en) 2021-10-26

Similar Documents

Publication Publication Date Title
CN109922020B (en) Low-computation-complexity orthogonal time-frequency space modulation balancing method
CN108833311B (en) Transform domain quadratic estimation method combining time domain clustering denoising and equalization judgment
Zhou et al. Finite-alphabet based channel estimation for OFDM and related multicarrier systems
CN106302298B (en) A method of eliminating OFDM underwater sound communication system clipped noise
CN101945066B (en) Channel estimation method of OFDM/OQAM system
CN110266617B (en) Multipath channel estimation method of super-Nyquist system
Han et al. Iterative per-vector equalization for orthogonal signal-division multiplexing over time-varying underwater acoustic channels
CN107483373B (en) Anti-multipath iterative weighting LMMSE channel estimation method and device
CN113452641B (en) FBMC channel estimation method, system, computer equipment and terminal
Luzio et al. SC-FDE for offset modulations: an efficient transmission technique for broadband wireless systems
CN107359906B (en) Method for suppressing impulse noise in low-voltage power line communication system
CN112003804B (en) Multipath multivariate differential chaotic shift keying iterative receiving method
CN113381951A (en) MFTN joint channel estimation and equalization method under time-frequency conversion selective fading channel
CN111865863A (en) RNN neural network-based OFDM signal detection method
JP4352035B2 (en) OFDM demodulator, method and program
CN113852575A (en) Iterative OTFS symbol detection method based on time domain channel equalization assistance
CN115664898A (en) OFDM system channel estimation method and system based on complex convolution neural network
CN107181704B (en) Low-complexity weighted LS soft iteration mobile channel estimation method
CN113556305B (en) FBMC iterative channel equalization method and system suitable for high-frequency selective fading
CN110324271A (en) A kind of F-OFDM system transceiver design method of compressed sensing based clipping
CN112636855A (en) OFDM signal detection method
WO2019075918A1 (en) Method of combining channel estimation to avoid pts transmitting sideband information
Khan et al. Data aided channel estimation for MIMO-OFDM wireless systems using reliable carriers
CN111147157B (en) Soft interference elimination Turbo equalization method for orthogonal signal division multiplexing on underwater acoustic channel
CN111431837B (en) OFDM signal iteration detection method for coping with interference between subcarriers and symbols

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
TR01 Transfer of patent right

Effective date of registration: 20230707

Address after: Room 3202, Building B, Wangzuo International City, No. 1 Tangyan Road, High tech Zone, Xi'an City, Shaanxi Province, 710075

Patentee after: Shaanxi Zexin Technology Co.,Ltd.

Address before: 710071 Xi'an Electronic and Science University, 2 Taibai South Road, Shaanxi, Xi'an

Patentee before: XIDIAN University

TR01 Transfer of patent right