CN116016061A - Short wave double-selection channel double-iteration Turbo equalization method for high maneuvering platform - Google Patents
Short wave double-selection channel double-iteration Turbo equalization method for high maneuvering platform Download PDFInfo
- Publication number
- CN116016061A CN116016061A CN202211624805.4A CN202211624805A CN116016061A CN 116016061 A CN116016061 A CN 116016061A CN 202211624805 A CN202211624805 A CN 202211624805A CN 116016061 A CN116016061 A CN 116016061A
- Authority
- CN
- China
- Prior art keywords
- iteration
- turbo
- sequence
- internal
- distribution
- 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
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 58
- 238000005562 fading Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 8
- 239000000654 additive Substances 0.000 claims abstract description 4
- 230000000996 additive effect Effects 0.000 claims abstract description 4
- 238000009826 distribution Methods 0.000 claims description 75
- 230000008569 process Effects 0.000 claims description 35
- 238000013016 damping Methods 0.000 claims description 23
- 230000005540 biological transmission Effects 0.000 claims description 16
- 230000009977 dual effect Effects 0.000 claims description 14
- 238000012804 iterative process Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 108091026890 Coding region Proteins 0.000 claims description 4
- 238000004891 communication Methods 0.000 abstract description 15
- 230000006870 function Effects 0.000 description 16
- 238000001514 detection method Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000005433 ionosphere Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- BTCSSZJGUNDROE-UHFFFAOYSA-N gamma-aminobutyric acid Chemical compound NCCCC(O)=O BTCSSZJGUNDROE-UHFFFAOYSA-N 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Abstract
The invention belongs to the field of communication, and relates to a short wave double-selection channel double-iteration Turbo equalization method of a high maneuvering platform; the method comprises the steps of adopting a transmitter to encode, interleave and modulate an information sequence, and outputting a transmitting symbol sequence; transmitting the sending symbol sequence through a fading channel, and generating a receiving symbol sequence after adding additive Gaussian white noise; processing the received symbol sequence by using a Turbo-BEP receiver, performing internal expected propagation iteration, processing the received symbol sequence and the prior feedback signal output by external Turbo iteration by using a linear minimum mean square error equalization algorithm, and approximately outputting a posterior feedback signal by using a cavity function; through external Turbo iteration, the posterior feedback signal is demapped, deinterleaved and decoded, and a decoded information sequence is output; the information sequence is interleaved and mapped to generate an a priori feedback signal. The invention improves the performance, simultaneously effectively reduces the average iteration times and improves the real-time performance of the double iteration structure.
Description
Technical Field
The invention belongs to the field of communication, and particularly relates to a short wave double-selection channel double-iteration Turbo equalization method of a high-mobility platform.
Background
The short-wave communication is widely applied to emergency rescue scenes such as rescue and relief work by virtue of the remote communication advantages that the short-wave communication is not limited by a network hub and an active relay, and the communication emergency command network is quickly established under extreme conditions such as damage and interruption of facilities such as basic communication, electric power and the like caused by overregional, long-distance and large-range oversized disasters or war.
A short wave channel is a typical channel with a dual-choice fading characteristic, and the main propagation path is sky wave propagation, and after a signal is sent out by an antenna, the signal is reflected back to the ground through an ionosphere. Because the height and density of the ionosphere are easily affected by factors such as day and night, seasons, climate, etc., the short wave communication environment is severe, and intersymbol interference (ISI, inter-symbol interference) caused by time and frequency dispersion exists. In a single carrier system, inter-symbol interference (ISI) due to multipath delay spread may be spread to tens or even hundreds of baseband symbols, and equalization technology is an effective method for coping with channel distortion and eliminating ISI.
Soft channel equalization or probabilistic channel equalization is a technique to mitigate inter-symbol interference caused by the dispersive nature of the channel. It provides the posterior probability of the transmission symbol estimated under given observation conditions. The two tasks of equalization and decoding were initially considered separately, but adding them to a turbo equalization scheme significantly improves performance. In turbo equalization, the equalizer and decoder exchange information according to log likelihood ratios (LLRs, log Likelihood Ratio). After one or more iterations, the channel decoder generates LLRs and sends them back to the equalizer as updated a priori information.
While early Turbo equalization techniques, such as maximum a posteriori probability (MAP) detectors using BCJR estimation, can approach channel capacity under properly designed coding schemes, the computational complexity of the BCJR algorithm and the number of trellis branches, M L The complexity is proportional to the increase of the channel tap number L and the constellation size M, and the memory requirement of each step of BCJR is increased along with the increase of the state number. The algorithm is therefore not applicable to the case of a high number of channel taps and high order modulation. Equalizer based on LMMSE is proposed subsequently [10] This is a suboptimal alternative. To reduce this complexitySome window versions have also been developed. In particular, m.tuchler et al propose a sliding window LMMSE algorithm, l.liu et al improve these results by replacing the sliding window with an extended window. Other technicians have also proposed some approximate windowing solutions to further reduce complexity. From the Bayesian perspective, the LMMSE algorithm replaces the discrete prior distribution of the transmission symbols with Gaussian prior, so that the algorithm can obtain Gaussian posterior distribution results.
The core task of Turbo equalization is to calculate the posterior distribution of the transmitted symbol x given the observation sequence y. However, in an actual short-wave communication system, since the characteristics of the short-wave channel are complex, the posterior distribution form of the transmitted symbol is difficult to calculate, and the true value of the distribution cannot be obtained by a simple mathematical method, an approximation method is required to solve the posterior distribution.
To obtain a more accurate posterior distribution, the prior distribution may be replaced by an approximation of the probability distribution to improve the complexity of the equalization algorithm in handling higher order modulations. I.e. the joint posterior probability of the transmitted bits is approximated using an equalization algorithm based on the expected propagation (EP, expectation Propagation). The algorithm can decompose posterior distribution, and ensure that each factor after decomposition can obey specific distribution as much as possible; or assuming that a certain parameter in the posterior distribution obeys a specific distribution. The desired propagation algorithm approximates the posterior probability with a family of gaussian functions and approximates the true posterior probability by iteratively performing moment matching. In recent years, based on accurate estimation characteristics of EP algorithm, extensive researches have been made in wireless communication systems, including channel estimation, multiple-Input multiple-Output (MIMO) system receiver signal detection technology, and the like. The Expected Propagation (EP) algorithm is a deterministic approximation method, and theoretical studies show that although a deterministic approximation does not yield an accurate value of the posterior distribution, the expected propagation-based approximation method is closer to the true distribution than the LMMSE criterion-based approximation.
The idea of the EP algorithm is combined with the equalization technology and applied to the double-selection fading channel, so that the gain can be effectively improved through the double-iteration structure, but the decoding delay caused by the calculation complexity of the double-iteration structure is ignored by another problem. In the case of excessive number of iterations, although the equalization performance is guaranteed, the improvement of complexity is unacceptable.
Disclosure of Invention
Based on the problems existing in the prior art, the invention provides a short wave double-selection channel double-iteration Turbo equalization method of a high maneuvering platform, which comprises the following steps:
adopting a transmitter to encode, interweave and modulate the information sequence to form a processed transmitting symbol sequence;
the sending symbol sequence is transmitted through a fading channel, and after additive Gaussian white noise is added, a receiving symbol sequence is generated;
processing the received symbol sequence by using a Turbo-BEP receiver, performing internal expected propagation iteration, processing the received symbol sequence and the prior feedback signal output by external Turbo iteration by using a linear minimum mean square error (LMSE) equalization algorithm, and approximately outputting a posterior feedback signal by using a cavity function; after the posterior feedback signal is demapped, deinterleaved and decoded through external Turbo iteration, a decoded information sequence is output; the information sequence is interleaved and mapped to generate an a priori feedback signal.
Further, the transmitter is adopted to encode, interleave and modulate the information sequence, and the formation of the processed transmission symbol sequence comprises the encoding of the information sequence with a certain code rate to obtain a transmission encoding sequence; interleaving the transmission coding sequence to obtain an interleaved data sequence; and carrying out M-order modulation on the interleaved data sequence to obtain a transmitting symbol sequence.
Further, the process of the internal expected propagation iteration includes:
step 1) obtaining an initial parameter pair of a single-variable complex Gaussian product in each round of internal iteration process;
step 2) calculating a second moment of Gaussian index family approximate distribution of posterior probability of a transmitted symbol sequence according to a parameter pair in the current internal iterative process, wherein the second moment comprises a mean value and a variance;
step 3) edge distribution and outer distribution of each symbol in the Gaussian index family approximate distribution are calculated in sequence;
step 4) estimating the second moment of the posterior probability distribution of the transmitted symbol sequence according to the symbol edge distribution and the outer distribution;
step 5) calculating the product distribution of the Gaussian index family of the outer distribution and the Gaussian index family of the next inner iterative process, enabling the product distribution to be equal to the second moment of the posterior probability distribution of the estimated transmitted symbol sequence, and calculating an intermediate parameter pair;
step 6) based on the intermediate parameter pairs, updating the parameter pairs in the current internal iteration process by adopting damping factors, and returning to the step 2) until the internal iteration times are reached.
Further, in the step 6), the calculation formula of the intermediate parameter pair is expressed as:
wherein ,a first intermediate parameter representing the kth symbol during the 1+1th internal iteration,/a>A second intermediate parameter representing the kth symbol during the 1+1th internal iteration,/a second intermediate parameter representing the kth symbol during the 1+1th internal iteration>Representing the variance of the posterior probability distribution of the estimated transmitted symbol sequence during the first internal iteration, ++>Representing the mean value of the posterior probability distribution of the estimated transmitted symbol sequence in the first internal iteration process,/, and>mean value of cavity function>Representing the cavity function variance.
Further, in the step 6), the formula for updating the parameter pair in the current internal iteration process by using the damping factor is expressed as:
wherein ,a first parameter representing the kth symbol during the 1+1st internal iteration, beta represents the damping factor,a first intermediate parameter representing the kth symbol during the 1+1th internal iteration,/a>A first parameter representing a kth symbol in a first internal iteration; />A second parameter representing the kth symbol during the 1+1th internal iteration, +.>Representing the kth symbol in the 1+1st internal iteration processIs>A second parameter representing a kth symbol in the first internal iteration.
Further, the damping factor is dynamically set, expressed as:
wherein β represents a damping factor; t represents the number of outer Turbo iterations.
Further, the external Turbo iterative process includes:
step 11), calculating posterior feedback signals of a transmitted symbol sequence in the current external iteration process by using the updated parameter pairs of the internal iteration process;
step 12) inputting the posterior feedback signal into a demapper, calculating an external log-likelihood ratio and transmitting the external log-likelihood ratio to a channel decoder;
step 13) according to each bit soft output of the channel decoder, recalculating the prior probability distribution of each transmitted symbol in the decoder, and calculating the average value and variance thereof;
step 14) returns to step 11) and enters the next round of external iteration process until the number of external iterations is reached.
The invention has the beneficial effects that:
the invention aims at a short wave fading channel of a high mobility platform, and adopts combined signal equalization and detection based on an LMMSE algorithm to process signal distortion and ISI phenomena in high frequency transmission. The algorithm is applied to a short wave fading channel, and the equalization complexity of high-order modulation is effectively reduced. While equalization performance benefits from a dual iteration structure, another problem is the excessive number of iterations. Aiming at the situation, the invention provides a double-iteration balanced Turbo method based on dynamic improvement, which ensures better performance than the traditional LMMSE algorithm. On the premise of improving the performance, the iteration times of the double-iteration structure are reduced, so that the calculation complexity is reduced. Simulation experiments show that the performance of the double-iteration Turbo equalization method based on the invention is superior to that of the traditional LMMSE algorithm, and when the number of outer layer iterations reaches 5, the performance is approximately converged. When the outer layer iteration times are fixed to be 5 times, the shortwave double-selection channel double-iteration Turbo equalization method provided by the invention can reduce the calculation complexity under the condition of almost no performance loss, and is superior to the common double-iteration equalization method.
Drawings
FIG. 1 is a flow chart of a dual iterative equalization scheme employed in an embodiment of the present invention;
FIG. 2 shows a transmitting end structure according to an embodiment of the present invention;
fig. 3 is a Turbo receiver according to an embodiment of the present invention;
fig. 4 shows an EP equalization receiver according to an embodiment of the present invention;
FIG. 5 is a block diagram of a dual iterative equalization architecture employed in an embodiment of the present invention
FIG. 6 is a simulation diagram of a Turbo-BEP equalization effect employed in an embodiment of the present invention;
FIG. 7 is a diagram of a three-dimensional effect of Turbo-BEP and Turbo-LMMSE equalization employed in an embodiment of the present invention;
FIG. 8 is a simulation diagram of BER curves of different inner iteration times S and outer iteration times T according to the embodiment of the present invention;
fig. 9 is a simulation diagram of a dual iterative equalization BER curve employing a modified EP algorithm employed in an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Fig. 1 is a flow chart of a short-wave dual-selection channel dual-iteration Turbo equalization method of a high mobility platform according to an embodiment of the present invention, as shown in fig. 1, the method includes:
101. adopting a transmitter to encode, interweave and modulate the information sequence to form a processed transmitting symbol sequence;
in the embodiment of the invention, the information sequence is coded with a certain code rate to obtain a transmission coding sequence; interleaving the transmission coding sequence to obtain an interleaved data sequence; and carrying out M-order modulation on the interleaved data sequence to obtain a transmitting symbol sequence.
The structure block diagram of the transmitting end adopted by the embodiment of the invention is shown in fig. 2, and an information sequence a= [ a ] formed by K information bits 1 ,…,a K ] T Encoding with code rate r=k/V as b= [ b ] 1 ,…,b V ] T The total code length after coding is V, and c= [ c ] is obtained after interleaving 1 ,…,c V ] T . M-element constellation modulation is carried out on the interleaved data sequence to obtain a transmission symbol sequence x= [ x ] 1 ,…,x N ] T Where N represents the number of transmitted symbols, n= [ V/log 2 M]。
102. The sending symbol sequence is transmitted through a fading channel, and after additive Gaussian white noise is added, a receiving symbol sequence is generated;
in the embodiment of the present invention, the transmission symbol sequence passes through a channel h= [ h ] in a fourier transform form x=r (x) +ji (x) of a data frame 1 ,…,h L ]Transmission, L denotes the channel length, each transmitted symbol being denoted as wherein />Represents->A set of symbols of a rank constellation.
103. Processing the received symbol sequence by using a Turbo-BEP receiver, performing internal expected propagation iteration, processing the received symbol sequence and the prior feedback signal output by external Turbo iteration by using a linear minimum mean square error (LMSE) equalization algorithm, and approximately outputting a posterior feedback signal by using a cavity function; after the posterior feedback signal is demapped, deinterleaved and decoded through external Turbo iteration, a decoded information sequence is output; the information sequence is interleaved and mapped to generate an a priori feedback signal.
wherein ,yk Represents the kth received symbol, h l Represents the first transmission channel, x k-l+1 Represents the k-l+1 th transmission symbol, w k Representative variance isThe average energy per symbol and energy per bit of the transmitted signal are respectively E s and Eb And (3) representing. />X when k < 1 or k > N k =0。
The corresponding matrix form of the formula (1) is y=hx+w, and H represents a channel matrix; w represents a noise matrix, and is specifically expressed as follows:
the key of Turbo equalization at the receiving end is that the equalizer and decoder perform iterative exchange of information on the same set of received symbols. The decoder inputs information LLRL calculated for the equalizer E (b t Y), i.e., information output by the equalizer during iteration, for iteration, abbreviated as L E (b t ) The decoder calculates an estimate of the information bits after one or more iterationsExternal LLR of coded bits:
these are with a priori probabilities p D (x k ) The represented LLRs are remapped as shown in fig. 3 (pi and pi -1 Representing the interleaving map and its inverse) and returned to the equalizer as an updated prior probability. This process is performed a given maximum number of iterations T until convergence.
On the basis of a traditional Turbo receiver, the posterior probability of transmitting a symbol vector x is expressed as:
if the prior information is unknown to the channel decoder, then the transmitted symbols are assumed to satisfy the uniform distribution, and the prior probability thereof satisfies:
wherein ,δ(xk -x) represents an impulse function.
This assumption is used to provide the equalizer with a priori information before the Turbo receiver iterates.
In the case of a priori known channel information, the LMMSE equalizer provides a solution that approximates the posterior probability with a mean valueSum of variances->The discrete prior probability p (x) in (3) is approximated by an independent gaussian function, the approximate distribution being expressed as: />
μ MMSE and ΣMMSE The mean value and covariance matrix of the symbol sequence are as follows:
In the iterative process, assuming equal probability of the symbol of the kth prior information, the out-of-symbol information calculated by (5) is passed to the channel decoder, equalizer feedback statistics and />They are derived from updated prior probabilities p D (x k ) The calculation method is as follows:
in the technology, a desired propagation algorithm is adopted to approximate the joint posterior probability of transmission bits, and the desired propagation (EP) is a technology in Bayesian machine learning, is essentially an accurate estimation algorithm based on message passing, and approximates complex probability distribution through exponential family distribution. Given a statistical distribution containing a latent variable x and an observed variable y, a complex posterior probability distribution consisting of I non-negative factor products is:
wherein Index family->Representing a set of known approximate probability functions, the set having sufficient statistics Φ (x) and obeying an exponential distribution. Because (11) contains a substance other than +.>Is not a negative factor t of (2) i (x) Therefore, it is not feasible to directly calculate (11). Adopting EP algorithm, according to "moment" matching process, making approximate iterative approximation of objective function not belonging to index family, i.e. using +.>Is->To replace t i (x) In this process, optimize +.>To obtain a precise approximation solution:
the algorithmic process of a conventional expected propagation algorithm is shown in table 1, where l represents the number of algorithm iterations.
Defining i=1, … I, the expected number of propagation iterations l=1, …, S, initializing an approximate distribution factor;
1) Calculating an approximate solution q (x) from (12);
2) For the non-negative factor i=1, … I and the expected number of propagation iterations i=1, …, S, the following steps are taken:
3) Calculating probability distribution and corresponding moment:
4) By "moment" matching:
5) Ending the iterative process;
For practical communication systems, the EP algorithm model in (11) has high complexity, and in order to apply it to the communication field, an approximation of the EP algorithm is required. The use of the approximate EP algorithm in the communication system can effectively improve the accuracy of the algorithm. In case the channel conditions are known a priori, the structure of the true distribution is preserved as much as possible and the unknown factor t is calculated i (x) From the expression of the probability distribution, the following gaussian index family is a reasonable approximation to (3):
wherein the product of the indicator functions is replaced by the product of a univariate complex gaussian, each function being replaced by a parameter pair (gamma k ,Λ k ) K=1, N corresponds to. The mean and variance of q (x) are:
this solution is the Block EP (BEP) algorithm, which shows a structure similar to the LMMSE in (5).
The equalizer structure used by the conventional algorithm using BEP algorithm is shown in fig. 4, and q is calculated in S iterations (l) (x k ) Refinement is performed to improve accuracy.
Based on the above analysis, the BEP equalizer in embodiments of the present invention may be further improved by using a turbo scheme. The equalizer thus produced, denoted Turbo-BEP, only needs to be replaced by the result of BEP algorithm 2 in (13) on the basis of Turbo-LMMSE. The detailed implementation of Turbo-BEP is contained in algorithm 2. The T-BEP may actually be considered as a Turbo equalizer with two loops. First, an internal iteration scheme, the BEP, is run. After S iterations of the BEP iteration process, the extrinsic LLR is provided to the decoder in the outer iterations, repeating T times. The outer distribution is approximated by a cavity function at the end of the EP algorithm. Limitations from channel coding are exploited. The output of the channel decoder is used to initialize the BEP iterative process, and then its output is fed forward to the channel decoder. This is the main difference of the previous EP scheme used in the conventional art, where the estimation is refined using only the output of the decoder, as shown in fig. 5, the dual iterative Turbo structure of the present invention consists of inner and outer iterations, which can make the proposed inner loop nonexistent.
The internal expected propagation iteration process in the embodiment of the invention comprises the following steps:
step 1) obtaining an initial parameter pair of a single-variable complex Gaussian product in each round of internal iteration process;
step 2) calculating a second moment of Gaussian index family approximate distribution of posterior probability of a transmitted symbol sequence according to a parameter pair in the current internal iterative process, wherein the second moment comprises a mean value and a variance;
step 3) edge distribution and outer distribution of each symbol in the Gaussian index family approximate distribution are calculated in sequence;
step 4) estimating the second moment of the posterior probability distribution of the transmitted symbol sequence according to the symbol edge distribution and the outer distribution;
step 5) calculating the product distribution of the Gaussian index family of the outer distribution and the Gaussian index family of the next inner iterative process, enabling the product distribution to be equal to the second moment of the posterior probability distribution of the estimated transmitted symbol sequence, and calculating an intermediate parameter pair;
step 6) based on the intermediate parameter pairs, updating the parameter pairs in the current internal iteration process by adopting damping factors, and returning to the step 2) until the internal iteration times are reached.
The external Turbo iteration process comprises the following steps:
step 11), calculating posterior feedback signals of a transmitted symbol sequence in the current external iteration process by using the updated parameter pairs of the internal iteration process;
step 12) inputting the posterior feedback signal into a demapper, calculating an external log-likelihood ratio and transmitting the external log-likelihood ratio to a channel decoder;
step 13) according to each bit soft output of the channel decoder, recalculating the prior probability distribution of each transmitted symbol in the decoder, and calculating the average value and variance thereof;
step 14) returns to step 11) and enters the next round of external iteration process until the number of external iterations is reached.
To further illustrate the inner desired propagation iteration and outer Turbo iteration processes of the present invention, a similar description will be provided below in connection with algorithm 2.
Algorithm 2Turbo-BEP Algorithm
Defining the outer Turbo iteration times t=1, t, the inner EP iteration times l=1, S, the symbol digits k=1, N and the initialization parameter pairs
1) For the outer Turbo iteration times, t=1, T, executing step 2);
2) Step 3) is executed for the internal EP iteration number satisfying i=1, S;
3) With initial parameter pairsStarting the EP iteration, calculating the approximate solution q in the first internal iteration process in the formula (14) and the formula (15) under the initial condition (l) (x) Mean and variance (second moment) of (a) in the matrix.
4) Step 5) is performed for the number of sign bits to satisfy k=1, …, N;
5) Calculate q (l) (x) The kth symbol edge distribution in (a)And the outer distribution of the kth symbol +.>
wherein :
7) Order theAnd->Equal second moment (mean and variance) of the intermediate parameter pair, a calculation formula of the intermediate parameter pair can be obtained:
wherein ,a complex gaussian function first intermediate parameter representing the kth symbol in the 1 st +1 internal iteration,a second intermediate parameter of the complex Gaussian function representing the kth symbol in the 1+1th internal iteration process, ">Representing the variance of the posterior probability distribution of the estimated kth transmitted symbol sequence during the first internal iteration,/, for>Representing the mean value of the posterior probability distribution of the kth transmitted symbol sequence estimated in the first internal iteration,
8) Updating the parameter pairs with the intermediate parameter pairs, expressed as:
wherein β represents a damping factor.
9) The loop of steps 4 to 8 ends.
10 Ending the loop of steps 2 to 9.
11 Using parameter pairs after the end of an EP iterationThe final distribution q (x) in (14) is calculated. Will feed back p E (x k |y)=q (S+1)\k (x k ) Input demapper to compute external LLRsL E (b t I y) and then transmitted to the channel decoder. Recalculating each symbol p in the decoder from each bit soft output of the channel decoder D (x k ) Probability distribution of (2)And calculate the average +.>Sum of variancesIs given by (9) and (10).
12 Reinitializing the parameter pairs
13 Ending the loop of steps 1 to 12.
Table 3 compares in detail the complexity of the present invention (Turbo-BEP) and the Turbo-LMMSE. The Turbo-BCJR calculation complexity is also included. From a complexity point of view, as the values of M and L increase, the BCJR equalization method is not burdened, which is also consistent with the conclusions above. Where S' =s+1 and α are the complexity of the LDPC encoder.
Table 3 complexity comparison
The improved Turbo-BEP receiver and the corresponding iterative process thereof can bring excessive iterative times for a double-iterative structure, and the excessive iterative times can influence the received decoding delay.
The iteration number set by the iterative process defaults to s=10 times in EP. The damping factor β in EP iterations determines the convergence rate, and a more conservative damping factor value is typically chosen to ensure the accuracy of the algorithm.
In the conventional EP equalization algorithm, the damping factor is usually set to 0.1 and the convergence number s=10 is set to ensure accurate convergence at the EP iteration, but the convergence speed is slow. It is contemplated that the damping factor may be dynamically set to achieve lower decoding delay with reduced iteration times while ensuring some accuracy. Experiments show that in the Turbo iteration process of the previous three times, the average of the Turbo iterations is equal to the average of the Turbo iterationsThe balance performance is obviously improved, at the moment, a high damping factor can be adopted to accelerate convergence, the damping factor is dynamically reduced after three times to obtain accurate convergence, and a new damping factor value is obtained according to the thoughtThe dynamic setting of the damping factor can reduce S to 5 times within the allowable error range compared to the way the damping factor is fixed. This improved EP algorithm is called dynamic improved double iterative equalization (DI-DIE). />
The channel characteristics of radio waves in short wave communication include reflection, scattering and diffraction, which cause multipath fading problems, considering the influence of topography, atmosphere and the like. The Rayleigh model is a typical multipath communication channel model. The channel is modeled by a two-path Rayleigh model, the multipath delay is 2ms, and the doppler shift is 1Hz. The channel parameter settings are in accordance with the description of the short wave channel in ITU-R, all simulations are performed with complete knowledge of the channel state information, the specific parameters being shown in table 4.
Table 4 simulation parameters
Fig. 6 is a BEP equalizer versus LMMSE equalizer BER curve for the inner EP iteration s=10, standard equalization (t=0) versus the number of iterations t=1, 3,5, 10. In the case of separation of equalization and decoding (t=0), the BEP performance of the decoder without feedback is much lower than LMMSE. After 8 iterations, at ber=10 -2.5 In the case of (2), the signal to noise ratio of the BEP algorithm is improved by about 1.7dB compared to LMMSE.
Fig. 7 shows a three-dimensional graph of the inner iteration number, the bit error rate and the signal to noise ratio, and it can be seen from fig. 7 that the trend of the bit error rate decreases gradually with the increase of the outer iteration number under the condition of different inner EP iteration parameters S, and reaches approximate convergence when t=5.
Fig. 8 shows BER curves for different inner and outer iteration numbers S, T for better observation of convergence. It can be observed that the equalization performance can be considered to be approximately converged when the number of internal and external iterations is 5, and the subsequent iteration performance improvement has a larger gap compared with the previous 5 iterations. As can be seen from fig. 8, by increasing the number of iterations T when s=3, there is a relatively significant gain relative to the first iteration s=1. After s=3, the gain gradually decreases. When s=10, the performance does not significantly improve over s=5 as the number of loops increases, and BEP equalization can be found to converge around the 5 th iteration. Different values of the outer iteration parameter show this trend.
If the fixed inner and outer iteration number and T are preset when performing Turbo-BEP equalization, redundant computation and unnecessary decoding delay may be caused. The method can accelerate the convergence speed of the internal EP on the premise of ensuring the performance by using the improved dynamic damping factor on the basis of the initial preset iteration times, thereby reducing the iteration times on the premise of meeting the balance performance. As shown in fig. 9, the DI-DIE algorithm has little loss of equalization performance when convergence is reached, i.e., the outer iteration number is 5 times, compared with the Turbo-BEP algorithm, and the inner iteration number S is reduced from 10 times to 5 times compared with the Turbo-BEP algorithm of the original fixed damping factor.
The invention aims at a short wave fading channel of a high mobility platform, and adopts combined signal equalization and detection based on an LMMSE algorithm to process signal distortion and ISI phenomena in high frequency transmission. The algorithm is applied to a short wave channel of a maneuvering platform, and the equalization complexity of high-order modulation is effectively reduced. While equalization performance benefits from a dual iteration structure, another problem is the excessive number of iterations. Aiming at the situation, the invention also provides a joint equalization and signal detection algorithm based on dynamic improved double iterative equalization (DI-DIE), which ensures better performance than the LMMSE algorithm. On the premise of improving the performance, the iteration times of the double-iteration structure are reduced, so that the calculation complexity is reduced. Simulation experiments show that the performance of the secondary iterative equilibrium detection algorithm based on the Turbo-BEP algorithm is superior to that of the LMMSE algorithm, and the performance approximately converges when the outer layer iteration number reaches 5 times. When the outer layer iteration number is fixed to be 5 times, the detection algorithm based on the DI-DIE algorithm provided by the invention can reduce the computational complexity by improving the inner EP iteration number under the condition of almost no performance loss, and is superior to the common double-iteration balanced detection algorithm.
Those of ordinary skill in the art will appreciate that all or part of the steps in the various methods of the above embodiments may be implemented by a program to instruct related hardware, the program may be stored in a computer readable storage medium, and the storage medium may include: ROM, RAM, magnetic or optical disks, etc.
Although embodiments of the present invention have been shown and described, it will be understood by those skilled in the art that various changes, modifications, substitutions and alterations can be made therein without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (7)
1. The utility model provides a high maneuver platform shortwave double-selection channel double-iteration Turbo equalization method which is characterized in that the method comprises the following steps:
adopting a transmitter to encode, interweave and modulate the information sequence to form a processed transmitting symbol sequence;
the sending symbol sequence is transmitted through a fading channel, and after additive Gaussian white noise is added, a receiving symbol sequence is generated;
processing the received symbol sequence by using a Turbo-BEP receiver, performing internal expected propagation iteration, processing the received symbol sequence and the prior feedback signal output by external Turbo iteration by using a linear minimum mean square error (LMSE) equalization algorithm, and approximately outputting a posterior feedback signal by using a cavity function; after the posterior feedback signal is demapped, deinterleaved and decoded through external Turbo iteration, a decoded information sequence is output; the information sequence is interleaved and mapped to generate an a priori feedback signal.
2. The method for dual-iteration Turbo equalization of short-wave dual-selection channel of high mobility platform as recited in claim 1 wherein said using a transmitter to encode, interleave and modulate information sequences, forming processed transmit symbol sequences includes encoding said information sequences at a code rate to obtain transmit code sequences; interleaving the transmission coding sequence to obtain an interleaved data sequence; and carrying out M-order modulation on the interleaved data sequence to obtain a transmitting symbol sequence.
3. The high mobility platform short wave dual selection channel dual iteration Turbo equalization method of claim 1, wherein said process of internal desired propagation iteration comprises:
step 1) obtaining an initial parameter pair of a single-variable complex Gaussian product in each round of internal iteration process;
step 2) calculating a second moment of Gaussian index family approximate distribution of posterior probability of a transmitted symbol sequence according to a parameter pair in the current internal iterative process, wherein the second moment comprises a mean value and a variance;
step 3) edge distribution and outer distribution of each symbol in the Gaussian index family approximate distribution are calculated in sequence;
step 4) estimating the second moment of the posterior probability distribution of the transmitted symbol sequence according to the symbol edge distribution and the outer distribution;
step 5) calculating the product distribution of the Gaussian index family of the outer distribution and the Gaussian index family of the next inner iterative process, enabling the product distribution to be equal to the second moment of the posterior probability distribution of the estimated transmitted symbol sequence, and calculating an intermediate parameter pair;
step 6) based on the intermediate parameter pairs, updating the parameter pairs in the current internal iteration process by adopting damping factors, and returning to the step 2) until the internal iteration times are reached.
4. The method for high mobility platform short wave dual selection channel dual iterative Turbo equalization according to claim 3, wherein in said step 6), the calculation formula of the intermediate parameter pair is expressed as:
wherein ,a first intermediate parameter representing the kth symbol during the 1+1th internal iteration,/a>A second intermediate parameter representing the kth symbol during the 1+1th internal iteration,/a second intermediate parameter representing the kth symbol during the 1+1th internal iteration>Representing the variance of the posterior probability distribution of the estimated transmitted symbol sequence during the first internal iteration, ++>Representing the mean value of the posterior probability distribution of the estimated transmitted symbol sequence in the first internal iteration process,/, and>mean value of cavity function>Representing the cavity function variance.
5. The method for dual-iteration Turbo equalization of short-wave dual-selection channel of high mobility platform according to claim 3, wherein in said step 6), the formula for updating the parameter pair in the current internal iteration process by using the damping factor is expressed as:
wherein ,a first parameter representing the kth symbol during the 1+1th internal iteration, β representing the damping factor,/v>A first intermediate parameter representing the kth symbol during the 1+1th internal iteration,/a>A first parameter representing a kth symbol in a first internal iteration; />A second parameter representing the kth symbol during the 1+1th internal iteration, +.>A second intermediate parameter representing the kth symbol during the 1+1th internal iteration,/a second intermediate parameter representing the kth symbol during the 1+1th internal iteration>A second parameter representing a kth symbol in the first internal iteration.
7. A high mobility platform short wave dual selection channel dual iteration Turbo equalization method according to claim 1 or 3, wherein the process of external Turbo iteration comprises:
step 11), calculating posterior feedback signals of a transmitted symbol sequence in the current external iteration process by using the updated parameter pairs of the internal iteration process;
step 12) inputting the posterior feedback signal into a demapper, calculating an external log-likelihood ratio and transmitting the external log-likelihood ratio to a channel decoder;
step 13) according to each bit soft output of the channel decoder, recalculating the prior probability distribution of each transmitted symbol in the decoder, and calculating the average value and variance thereof;
step 14) returns to step 11) and enters the next round of external iteration process until the number of external iterations is reached.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211624805.4A CN116016061A (en) | 2022-12-16 | 2022-12-16 | Short wave double-selection channel double-iteration Turbo equalization method for high maneuvering platform |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211624805.4A CN116016061A (en) | 2022-12-16 | 2022-12-16 | Short wave double-selection channel double-iteration Turbo equalization method for high maneuvering platform |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116016061A true CN116016061A (en) | 2023-04-25 |
Family
ID=86037355
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211624805.4A Pending CN116016061A (en) | 2022-12-16 | 2022-12-16 | Short wave double-selection channel double-iteration Turbo equalization method for high maneuvering platform |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116016061A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116723069A (en) * | 2023-08-08 | 2023-09-08 | 华侨大学 | Multi-module iterative Turbo equalization method, device, equipment, server and medium |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003092170A1 (en) * | 2002-04-26 | 2003-11-06 | Kongsberg Defence Communications As | Method and apparatus for the reception of digital communication signals |
CN1860690A (en) * | 2003-07-24 | 2006-11-08 | 科达无线私人有限公司 | Method and system for communication in a multiple access network |
US20150372785A1 (en) * | 2013-02-13 | 2015-12-24 | Orange | A method and a device for predicting the performance of a communication system over a transmission channel |
CN106301517A (en) * | 2016-08-10 | 2017-01-04 | 清华大学 | The satellite multi-beam joint-detection propagated based on expectation and interpretation method and system |
US20190068294A1 (en) * | 2017-08-25 | 2019-02-28 | Yahong Rosa Zheng | Turbo receivers for multiple-input multiple-output underwater acoustic communications |
CN111682924A (en) * | 2020-04-07 | 2020-09-18 | 杭州电子科技大学 | Bidirectional frequency domain Turbo equalization method adopting expected propagation |
EP3800813A1 (en) * | 2019-10-03 | 2021-04-07 | Thales | Method and device for predicting the performance of a receiver in a communication system |
CN113242190A (en) * | 2021-04-13 | 2021-08-10 | 华南理工大学 | Multichannel communication minimum bit error rate Turbo equalization method based on posterior soft symbol |
CN115208480A (en) * | 2022-06-30 | 2022-10-18 | 哈尔滨工程大学 | Under-ice underwater acoustic communication method based on joint message transfer |
-
2022
- 2022-12-16 CN CN202211624805.4A patent/CN116016061A/en active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003092170A1 (en) * | 2002-04-26 | 2003-11-06 | Kongsberg Defence Communications As | Method and apparatus for the reception of digital communication signals |
CN1860690A (en) * | 2003-07-24 | 2006-11-08 | 科达无线私人有限公司 | Method and system for communication in a multiple access network |
US20150372785A1 (en) * | 2013-02-13 | 2015-12-24 | Orange | A method and a device for predicting the performance of a communication system over a transmission channel |
CN106301517A (en) * | 2016-08-10 | 2017-01-04 | 清华大学 | The satellite multi-beam joint-detection propagated based on expectation and interpretation method and system |
US20190068294A1 (en) * | 2017-08-25 | 2019-02-28 | Yahong Rosa Zheng | Turbo receivers for multiple-input multiple-output underwater acoustic communications |
EP3800813A1 (en) * | 2019-10-03 | 2021-04-07 | Thales | Method and device for predicting the performance of a receiver in a communication system |
CN111682924A (en) * | 2020-04-07 | 2020-09-18 | 杭州电子科技大学 | Bidirectional frequency domain Turbo equalization method adopting expected propagation |
CN113242190A (en) * | 2021-04-13 | 2021-08-10 | 华南理工大学 | Multichannel communication minimum bit error rate Turbo equalization method based on posterior soft symbol |
CN115208480A (en) * | 2022-06-30 | 2022-10-18 | 哈尔滨工程大学 | Under-ice underwater acoustic communication method based on joint message transfer |
Non-Patent Citations (3)
Title |
---|
CONGJI YIN; WENJIANG FENG; JUNBING LI; GUOJUN LI: "Serially Connected Bidirectional Double EP-Based DFE for Turbo Equalization", IEEE, 26 July 2022 (2022-07-26) * |
李国军;龙锟;叶昌荣;梁佳文: "高速移动环境下低复杂度OTSM迭代rake均衡方法", 通信学报, 25 October 2022 (2022-10-25) * |
王晓春;陈佳怡;董超;: "低复杂度频域迭代均衡技术", 无线电工程, no. 02, 5 February 2020 (2020-02-05) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116723069A (en) * | 2023-08-08 | 2023-09-08 | 华侨大学 | Multi-module iterative Turbo equalization method, device, equipment, server and medium |
CN116723069B (en) * | 2023-08-08 | 2023-12-05 | 华侨大学 | Multi-module iterative Turbo equalization method, device, equipment, server and medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113242189B (en) | Adaptive equalization soft information iteration receiving method combined with channel estimation | |
CN109246039B (en) | Soft information iteration receiving method based on bidirectional time domain equalization | |
US20060256888A1 (en) | Multi input multi output wireless communication reception method and apparatus | |
CN104202271B (en) | Based on the iteration equalizing method handled by survivor path in Direct Sequence Spread Spectrum Communication | |
CN113242190B (en) | Multichannel communication minimum bit error rate Turbo equalization method based on posterior soft symbol | |
CN109981501B (en) | Underwater sound direct self-adaptive MIMO communication method | |
CN112866151A (en) | Underwater sound MPSK signal blind Turbo equalization method based on channel blind estimation | |
CN116016061A (en) | Short wave double-selection channel double-iteration Turbo equalization method for high maneuvering platform | |
CN104410593B (en) | Numerical chracter nonlinearity erron amendment equalization methods based on decision-feedback model | |
Xu et al. | Spatial and time-reversal diversity aided least-symbol-error-rate turbo receiver for underwater acoustic communications | |
CN101119177A (en) | Bit-symbol signal processing method for coherent communication machine | |
CN105553903A (en) | Adaptive turbo equalization method, equalizer and underwater acoustic communication system | |
CN108111446B (en) | Receiver equalization module and equalization method | |
Santos et al. | A double EP-based proposal for turbo equalization | |
CN111682924B (en) | Bidirectional frequency domain Turbo equalization method adopting expected propagation | |
CN112039809B (en) | Block iterative equalizer based on mixed soft information and bidirectional block iterative equalizer | |
Santos et al. | Block expectation propagation equalization for ISI channels | |
Abdulkader et al. | Neural networks-based turbo equalization of a satellite communication channel | |
CN107659523B (en) | BPSK modulation equalization system and method in wireless mobile communication | |
CN111901262A (en) | High-order modulation Turbo time domain equalization algorithm suitable for short-wave communication | |
Rahman et al. | Iterative soft decision based complex K-best MIMO decoder | |
Pan et al. | Equalization Techniques for Unmanned Aerial Vehicles Communication Based on Early Termination of Iteration | |
Kwon et al. | SVR-based blind equalization on HF channels with a Doppler spread | |
Fki et al. | New criteria for blind equalization based on pdf fitting | |
CN102355434B (en) | Orthogonal wavelet transform constant modulus blind equalization algorithm based on chaos and steepest descent joint optimization |
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 |