阅读说明：本技术 散射通信系统的均衡方法、系统、计算机设备、处理终端 (Equalization method and system of scattering communication system, computer equipment and processing terminal ) 是由 郭漪 陈镇涛 毛晶晶 孙柏昶 刘刚 于 2021-06-18 设计创作，主要内容包括：本发明属于散射通信系统技术领域,公开了一种散射通信系统的均衡方法、系统、计算机设备、处理终端,给出SC-FDE系统模型,并针对非线性均衡算法对数据帧结构进行设计；对基于MMSE的HDFE均衡算法进行分析；提出一种HDFE-MMSE-RISIC均衡算法,利用HDFE均衡器完成数据符号的判决,提高了均衡算法的性能。本发明提出的HDFE-MMSE-RISIC均衡算法,通过HDFE均衡器的反馈滤波计算,提高了均衡计算过程的判决反馈结果准确性,之后利用该判决结果进行RISI的估计,从而进一步提高残余码间干扰结果的正确性。通过基于SC-FDE系统下的散射信道下的仿真分析,该改进的均衡算法有较为明显的误码性能的提升。(The invention belongs to the technical field of a scattering communication system, and discloses an equalization method, a system, computer equipment and a processing terminal of the scattering communication system, wherein an SC-FDE system model is provided, and a data frame structure is designed aiming at a nonlinear equalization algorithm; analyzing an HDFE equalization algorithm based on MMSE; the HDFE-MMSE-RISIC equalization algorithm is provided, the judgment of data symbols is completed by the HDFE equalizer, and the performance of the equalization algorithm is improved. The HDFE-MMSE-RISIC equalization algorithm provided by the invention improves the accuracy of a decision feedback result in the equalization calculation process through the feedback filtering calculation of the HDFE equalizer, and then utilizes the decision result to estimate RISI, thereby further improving the accuracy of a residual intersymbol interference result. Through simulation analysis under a scattering channel based on an SC-FDE system, the improved equalization algorithm has obvious improvement on error code performance.)

1. An equalization method for a scatterometry communication system, the equalization method comprising:

providing an SC-FDE system model, and designing a data frame structure aiming at a nonlinear equalization algorithm;

analyzing an HDFE equalization algorithm based on MMSE, wherein the algorithm is calculated through feedforward frequency domain equalization and feedback filtering of a time domain;

the HDFE-MMSE-RISIC equalization algorithm utilizes the judgment result of the HDFE equalizer on the data symbol to accurately estimate the RISI interference.

2. The equalization method for a scattercommunications system according to claim 1, wherein said SC-FDE system model comprises: the method comprises the steps that input data are subjected to LDPC coding, constellation mapping and the like to obtain complex transmission data, then a guard interval CP is added for sending, after the CP is removed from the data at a receiving end, a frequency domain receiving signal is obtained through FFT calculation, then channel estimation and frequency domain equalization calculation are carried out, finally the signal is converted into a time domain through IFFT, and signal detection is completed through demapping, channel decoding and the like, so that the design of the whole single carrier frequency domain equalization system is completed.

3. The equalization method for a scattercommunications system according to claim 1, wherein said designing a data frame structure for a non-linear equalization algorithm comprises: data transmission symbols are expanded by using an all-zero sequence, the length of the all-zero sequence is defined as L, and L is ensured to be more than or equal to L_h，L_hIs the channel impulse response CIR order;

the data to be transmitted is represented as x, and after zero sequence expansion with length L, the transmitted data block s can be represented as:

wherein, M represents the number of effective transmission data, and N is the size of the nonlinear equalization transmission data block, that is, N is M + L;

when the data block passes through a multipath channel h, a received signal after linear convolution is obtained:

expressed by using a matrix vector form, the following results are obtained:

wherein the content of the first and second substances,

since the second half of the data vector s is an all-zero sequence L long, the matrix is thenCan be expressed as Λ, formulaExpressed as:

r＝Λs+w；

wherein Λ is a Toepliz matrix:

and (3) converting the frequency domain through FFT calculation to obtain:

R＝HS+W；

wherein R ═ F_Nr，S＝F_Ns，W＝F_Nw，R represents DFT transformation of a received signal, S represents a zero padding expanded frequency domain transformation value of transmitted data x, W represents DFT transformation of a noise vector, and H obtains a calculation result as a diagonal matrix according to Toepliz matrix properties; f_NDenote an N × N DFT matrix, and the elements in the matrix are represented as:

then the formula R ═ HS + W is further simplified to give:

R_k＝H_kS_k+W_k＝H_kX_k+W_k,k＝0,1,...,N-1。

4. the equalization method for a scattered communication system according to claim 1, wherein in step two, the MMSE-based HDFE equalization algorithm comprises:

1) calculating time domain feedback coefficient g by using channel estimation result and signal-to-noise ratio estimation value_FBAnd then obtaining a frequency domain representation G of the feedback coefficient by FFT_FB；

2) Calculating to obtain feedforward frequency domain equalization coefficient value G based on MMSE_FFComplex multiplication calculation of a feedforward filter is completed by using the equilibrium coefficient, then IFFT conversion is carried out on the obtained linear frequency domain calculation result, and the first M effective data lengths are intercepted to obtain an equilibrium calculation value of a time domain;

3) calculating a time domain feedback value, and performing symbol-by-symbol summation feedback on the time domain equalization result obtained in the step 2);

4) carrying out symbol-by-symbol selection judgment on the feedback result obtained in the step 3);

5) and (4) repeating the steps 3) and 4) to perform feedback and judgment symbol by symbol, and performing M times of circulation until the time domain equalization result obtained in the step 2) completely completes feedback filtering calculation to obtain a final time-frequency domain mixed judgment feedback calculation result.

5. The equalizing method for a scattering communication system as claimed in claim 4, wherein in step 4), the feedback result obtained in step 3) is subjected to symbol-by-symbol selection decision to obtain a decided symbol sequence through time-domain feedback calculation, and then the L-long known sequence obtained in step 3) is left-shifted by one bit and removed, and is added and replaced with the symbol after selection decision to re-form the decided symbol sequence with L-length.

6. The equalization method for a scattered communication system according to claim 4, wherein the MMSE-based HDFE equalization algorithm specifically comprises:

receiving frequency domain data after equalization coefficient is G_FFThe feedforward frequency domain filter of (1) is obtained by performing complex multiplication calculation:

Y_k＝R_kG_FF,k,k＝0,1,...N-1；

then, the time domain is converted through IFFT:

frequency selective fading channel H and feed forward filter G_FFIs expressed as U-HG_FFAssuming the length of the equivalent channel response after concatenation to be L_uTurning to the time domain to obtain:

that is, the time domain value of the transmission data after passing through the feedforward filter is:

wherein the content of the first and second substances,

when L is_uWhen L is less than or equal to L, the cyclic convolution is consistent with the linear convolution, and a signal transmission symbol with the length of M is obtained, namely:

the input to the feedback filter is the decided symbol and a known all-zero sequence:

defining the coefficient of the feedback filter as g_FBAnd has a length L_gThe time domain received signal after feedback is expressed as:

then toMaking symbol-by-symbol selection decision to obtain the decided symbol result

According to the minimum mean square error principle, the cost function here is expressed as:

when L is_gL is less than or equal to L, and the feedback coefficient is transformed by the frequency domain:

thus, it isAnd converting the frequency domain to a frequency domain through FFT to obtain:

whereinPower for transmitting data symbols; in pair typePerforming a gradient operation to obtain:

substituted typeComprises the following steps:

defining feedback coefficientsIn pair typePerforming gradient operation and substituting into formulaObtaining a feedback coefficient satisfying:

Ag_FB＝b；

wherein the content of the first and second substances,

wherein A is L_g×L_gOf the Toeplitz matrix, its first behaviorAnd b is L_gColumn vector of length of valueThe IDFT transform of (2) obtains the calculation results of the two by N-point FFT/IFFT.

7. The equalizing method for a scattered communication system according to claim 1, wherein the HDFE-MMSE-RISIC equalizing algorithm comprises:

(1) carrying out HDFE equalization based on MMSE to obtain time domain equalization data with the length of N after selection judgment, and then intercepting the effective result of the previous M;

(2) adding a known all-zero sequence to the obtained HDFE equalization result, thereby expanding the data to N lengths again, and then carrying out FFT (fast Fourier transform) of N points to obtain a frequency domain value X of a judged symbol_k；

(3) Calculating frequency-domain values of residual intersymbol interferenceObtaining the RISI value of the time domain after IFFT transformation, and intercepting the effective data length of the former M;

(4) after the HDFE equalization calculation based on MMSE, obtaining a feedback equalization result with the same length of M, and performing summation compensation operation on the time domain data to finish the calculation of the equalization algorithm of the whole HDFE-MMSE-RISIC;

the HDFE-MMSE-RISIC equalization algorithm specifically comprises the following steps:

in the calculation of the HDFE equalization algorithm based on MMSE, the calculation mainly comprises frequency domain feedforward filtering calculation and time domain feedback filtering calculation, and the calculation process of the feedforward filter is represented as follows:

and performing N-point IFFT on the frequency domain equalization calculation result, and intercepting the effective data with the length of M to obtain:

wherein the content of the first and second substances,the values after the pre-M points are intercepted after the N-point IFFT is performed are shown, that is, the DFT transformation matrix only containing NxNThe first M rows of (C);

in-situ typeWherein the first part represents the feedforward filtered data part and the second part represents the residual intersymbol interference value delta after MMSE equalization_nI.e., the portion that needs to be removed later, and the third portion represents colored noise due to frequency domain equalization of the HDFE

By the formula:

Ag_FB＝b；

calculating to obtain g_FB,nAnd selecting the decided symbolFiltering calculation is carried out to obtain a time domain feedback value

Is directed toResidual intersymbol interference delta of_nCarrying out residual interference estimation on a frequency domain, and when a data block with the length of N is subjected to HDFE (high-level data interchange Format) equalization calculation, equalizing data y_nCarrying out selection judgment to obtain effective transmission data with the length of M; then, expanding the judged sequence to N length by utilizing a known all-zero sequence, carrying out FFT calculation, and then calculating a RISI estimated value of a frequency domain by combining a noise estimated value:

in pair typeIFFT calculation is carried out to obtain a residual intersymbol interference value with the length of N, effective data with the length of M before is intercepted to obtain an RISI estimated value of a time domainIf the selection decision result is completely correct, that is to sayThe equalization results were:

8. a computer device, characterized in that the computer device comprises a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to carry out the steps of:

providing an SC-FDE system model, and designing a data frame structure aiming at a nonlinear equalization algorithm;

analyzing an HDFE equalization algorithm based on MMSE, wherein the algorithm is calculated through feedforward frequency domain equalization and feedback filtering of a time domain;

the HDFE-MMSE-RISIC equalization algorithm utilizes the judgment result of the HDFE equalizer on the data symbol to accurately estimate the RISI interference.

9. An information data processing terminal, characterized in that the information data processing terminal is used for implementing the equalization method of the scatter communication system according to any one of claims 1 to 7.

10. An equalization system of a scattercommunications system for implementing the equalization method of any one of claims 1 to 7, the equalization system of the scattercommunications system comprising:

the data frame structure design module is used for providing an SC-FDE system model and designing a data frame structure aiming at a nonlinear equalization algorithm;

the HDFE equalization algorithm analysis module is used for analyzing the HDFE equalization algorithm based on MMSE;

and the RISI interference estimation module is used for accurately estimating the RISI interference by the HDFE-MMSE-RISIC equalization algorithm by utilizing the judgment result of the HDFE equalizer on the data symbol.

Technical Field

The invention belongs to the technical field of scattering communication systems, and particularly relates to an equalization method, an equalization system, computer equipment and a processing terminal of a scattering communication system.

Background

At present: troposphere scattering communication is used as a communication mode with beyond visual range, and has the characteristics of strong confidentiality, strong interference resistance and the like, so that the troposphere scattering communication plays an important role in the field of military communication. However, the scattering communication has multipath effect and doppler effect, which often causes dispersion distortion of the received signal in time and frequency domains. And as the communication rate increases, its interference with the transmission signal becomes non-negligible.

The mainstream tropospheric scatter communication systems at present mainly include OFDM and SC-FDE systems. The OFDM carrier modulation scheme divides a fading channel into parallel subcarrier channels, so that the channel fading influence can be effectively resisted, and the OFDM carrier modulation scheme has the advantages of low implementation complexity, high spectrum utilization rate and the like. On the other hand, SC-FDE can achieve complexity and performance similar to OFDM by performing equalization calculation in the frequency domain, and is less affected by nonlinearity of transmitter power amplifier, and thus is receiving increasing attention. However, in the SC-FDE system, for a channel with severe frequency selective fading, the influence of intersymbol interference is more obvious, thereby reducing the performance of the conventional linear frequency domain equalization algorithm.

Through the above analysis, the problems and defects of the prior art are as follows:

in an SC-FDE system, aiming at a channel with serious frequency selective fading, the influence of intersymbol interference is obvious, and the performance of the traditional linear frequency domain equalization algorithm is reduced.

The difficulty in solving the above problems and defects is: in order to reduce ISI interference and improve the error performance of the system, it is necessary to improve the conventional equalization algorithm and design a reasonable equalizer structure, so as to further reduce the influence of intersymbol interference caused by the detected symbols. While the MMSE-based HDFE equalization algorithm completes equalization through feedback calculation, residual intersymbol interference still exists, and in a channel with severe fading, the interference becomes non-negligible, so that improvement on the equalization algorithm is needed. Firstly, an HDFE equalization algorithm based on MMSE needs to be studied in detail and theoretical analysis is combined, so that a residual intersymbol interference value in feedforward linear calculation of the equalization algorithm is obtained, and then a time domain feedback equalization calculation structure needs to be reasonably designed and built, so that the influence of the residual intersymbol interference is reduced through feedback calculation, and the performance of the equalization algorithm is improved.

The significance of solving the problems and the defects is as follows: the HDFE-MMSE-RISIC balanced algorithm improves the symbol selection and decision performance after equalization by utilizing the HDFE equalization algorithm based on MMSE, calculates the residual intersymbol interference after equalization on the basis, and further eliminates the intersymbol interference caused by the detected symbols by time domain feedback calculation. Compared with the traditional equalization algorithm, the improved HDFE-MMSE-RISIC equalization algorithm has the advantage that the equalization performance is further improved under the condition of a scattering channel with more serious frequency selective fading, so that better error code performance can be obtained.

Disclosure of Invention

Aiming at the problems in the prior art, the invention provides an equalization method, an equalization system, computer equipment and a processing terminal of a scattering communication system.

The invention is realized in such a way that an equalization method of a scattering communication system comprises the following steps:

providing an SC-FDE system model, and designing a data frame structure aiming at a nonlinear equalization algorithm; therefore, the performance of the equalization algorithm can be comparatively analyzed by combining the result of the bit error rate curve under the system;

step two, analyzing an HDFE equalization algorithm based on MMSE, wherein the algorithm reduces the interference of ISI through feedforward frequency domain equalization and feedback filtering calculation of a time domain; through theoretical research, residual intersymbol interference exists in the equalization algorithm, so that the RISI interference needs to be further analyzed in order to further improve the performance of the equalization algorithm;

and step three, providing an HDFE-MMSE-RISIC equalization algorithm, accurately estimating RISI interference by using a judgment result of an HDFE equalizer on a data symbol, reducing influence caused by intersymbol interference through summation compensation of a time domain, and improving the performance of the equalization algorithm.

Further, in step one, the SC-FDE system model includes:

the method comprises the steps that input data are subjected to LDPC coding, constellation mapping and the like to obtain complex transmission data, then a protection interval CP is added for sending, after the CP is removed from the data at a receiving end, frequency domain receiving signals are obtained through FFT calculation, then channel estimation and frequency domain equalization calculation are carried out, finally the signals are converted into a time domain through IFFT, demapping, channel decoding and the like are carried out, and signal detection is finished, so that the design of the whole single carrier frequency domain equalization system is finished.

Further, in the first step, the designing the data frame structure for the non-linear equalization algorithm includes:

data transmission symbols are expanded by using an all-zero sequence, the length of the all-zero sequence is defined as L, and L is ensured to be more than or equal to L_h，L_hIs the Channel Impulse Response (CIR) order;

the data to be transmitted is represented as x, and after zero sequence expansion with length L, the transmitted data block s can be represented as:

wherein M represents the number of effective transmission data, and N is the size of a nonlinear equalization transmission data block;

when the data block passes through a multipath channel h, a received signal after linear convolution is obtained:

expressed by using a matrix vector form, the following results are obtained:

wherein the content of the first and second substances,

since the second half of the data vector s is an all-zero sequence L long, the matrix is thenCan be expressed as Λ, equation (3) can be expressed as:

r＝Λs+w (5)

wherein Λ is a Toepliz matrix,

and (3) converting the frequency domain through FFT calculation to obtain:

R＝HS+W (7)

wherein R ═ F_Nr，S＝F_Ns，W＝F_Nw，R represents DFT transformation of a received signal, S represents a zero-filling expanded frequency domain transformation value of transmitted data x, W represents DFT transformation of a noise vector, and H obtains a calculation result as a diagonal matrix according to Toepliz matrix properties; f_NDenote an N × N DFT matrix, and the elements in the matrix are represented as:

then equation (7) is further simplified to yield:

R_k＝H_kS_k+W_k＝H_kX_k+W_k,k＝0,1,...,N-1 (9)

further, in step two, the MMSE-based HDFE equalization algorithm includes:

1) calculating time domain feedback coefficient g by using channel estimation result and signal-to-noise ratio estimation value_FBAnd then obtaining a frequency domain representation G of the feedback coefficient by FFT_FB；

2) Calculating to obtain feedforward frequency domain equalization coefficient value G based on MMSE_FFThe equalization coefficient is used for completing complex multiplication calculation of a feedforward filter, then IFFT conversion is carried out on the obtained linear frequency domain calculation result, and the first M effective data lengths are intercepted to obtain an equalization calculation value of a time domain;

3) calculating a time domain feedback value, and performing symbol-by-symbol summation feedback on the time domain equalization result obtained in the step 2);

4) carrying out symbol-by-symbol selection judgment on the feedback result obtained in the step 3);

5) and (4) repeating the steps 3) and 4) to perform symbol-by-symbol feedback and judgment, and performing M times of circulation until the time domain equalization result obtained in the step 2) completely completes feedback filtering calculation to obtain a final time-frequency domain mixed judgment feedback calculation result.

Further, in step 4), after performing symbol-by-symbol selection decision on the feedback result obtained in step 3), a decided symbol sequence subjected to time domain feedback calculation is obtained, and then the L-long known sequence obtained in step 3) is shifted left by one bit and removed, and is added and replaced with a symbol subjected to selection decision, so as to re-form a decided symbol sequence with L length.

Further, the MMSE-based HDFE equalization algorithm specifically includes:

receiving frequency domain data after equalization coefficient is G_FFThe feedforward frequency-domain filter of (1) is obtained by performing complex multiplication calculation:

Y_k＝R_kG_FF,k,k＝0,1,...N-1 (10)

then, the time domain is converted through IFFT:

frequency selective fading channel H and feed forward filter G_FFIs expressed as U-HG_FFAssuming the length of the equivalent channel response after concatenation to be L_uTurning to the time domain to obtain:

that is, the time domain value of the transmission data after passing through the feedforward filter is:

wherein the content of the first and second substances,

when L is_uWhen L is less than or equal to L, the cyclic convolution is consistent with the linear convolution, and a signal transmission symbol with the length of M effective is obtained, namely:

the input to the feedback filter is the decided symbol and a known all-zero sequence:

defining the coefficient of the feedback filter as g_FBAnd has a length L_gThe time domain received signal after feedback is expressed as:

then toMaking symbol-by-symbol selection decision to obtain the decided symbol result

According to the minimum mean square error principle, the cost function here is expressed as:

when L is_gL is less than or equal to L, and the feedback coefficient is transformed by the frequency domain:

then equation (18) is transformed to the frequency domain by the FFT to obtain:

whereinPower for transmitting data symbols; performing a gradient operation on equation (20) to obtain:

the substitution formula (20) comprises:

defining feedback coefficientsAnd carrying out gradient operation on the formula (22) and substituting the formula (19) to obtain a feedback coefficient satisfying:

Ag_FB＝b (23)

wherein the content of the first and second substances,

wherein A is L_g×L_gOf the Toeplitz matrix, its first behaviorAnd b is L_gColumn vector of length of valueThe IDFT transform of (2) obtains the calculation results of the two by N-point FFT/IFFT.

Further, in step three, the HDFE-MMSE-RISIC equalization algorithm includes:

(1) carrying out HDFE equalization based on MMSE to obtain time domain equalization data with the length of N after selection judgment, and then intercepting the effective result of the previous M;

(2) adding a known all-zero sequence to the obtained HDFE equalization result, thereby expanding the data to N lengths again, and then carrying out FFT (fast Fourier transform) of N points to obtain a frequency domain value X of a judged symbol_k；

(3) Calculating frequency-domain values of residual intersymbol interferenceObtaining the RISI value of the time domain after IFFT transformation, and intercepting the effective data length of the former M;

(4) after the HDFE equalization calculation based on MMSE, a feedback equalization result with the same length of M is obtained, and the time domain data is subjected to summation compensation operation, so that the calculation of the equalization algorithm of the whole HDFE-MMSE-RISIC is completed.

Further, the HDFE-MMSE-RISIC equalization algorithm specifically includes:

in the calculation of the HDFE equalization algorithm based on MMSE, the calculation mainly comprises frequency domain feedforward filtering calculation and time domain feedback filtering calculation, and the calculation process of the feedforward filter is represented as follows:

and performing N-point IFFT on the frequency domain equalization calculation result, and intercepting the effective data with the length of M to obtain:

wherein the content of the first and second substances,the values after the pre-M points are intercepted after the N-point IFFT is performed are shown, that is, the DFT transformation matrix only containing NxNThe first M rows of (C);

in equation (26), the first part represents the feedforward-filtered data part, and the second part represents the residual intersymbol interference value δ after MMSE equalization_nI.e., the portion that needs to be removed later, and the third portion represents colored noise due to frequency domain equalization of the HDFE

Passing through type

Ag_FB＝b (23)

Calculating to obtain g_FB,nAnd selecting the decided symbolFiltering calculation is carried out to obtain a time domain feedback value

Residual intersymbol interference δ for equation (26)_nCarrying out residual interference estimation on a frequency domain, and when a data block with the length of N is subjected to HDFE (high-level data interchange Format) equalization calculation, equalizing data y_nCarrying out selection judgment to obtain effective transmission data with the length of M; then, expanding the judged sequence to N length by utilizing a known all-zero sequence, carrying out FFT calculation, and then calculating a RISI estimated value of a frequency domain by combining a noise estimated value:

IFFT calculation is carried out on the formula (28) to obtain a residual intersymbol interference value with the length of N, and effective data with the length of M before is intercepted to obtain an RISI estimated value of a time domainIf the selection decision result is completely correct, that is to sayThe equalization results were:

it is a further object of the invention to provide a computer device comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of:

providing an SC-FDE system model, and designing a data frame structure aiming at a nonlinear equalization algorithm;

analyzing an HDFE equalization algorithm based on MMSE, wherein the algorithm is calculated through feedforward frequency domain equalization and feedback filtering of a time domain;

the HDFE-MMSE-RISIC equalization algorithm utilizes the judgment result of the HDFE equalizer on the data symbol to accurately estimate the RISI interference.

Another object of the present invention is to provide an information data processing terminal, which is used for implementing the equalization method of the scattering communication system.

Another object of the present invention is to provide an equalization system of a scatterometry system for implementing the equalization method of the scatterometry system, the equalization system of the scatterometry system including:

the data frame structure design module is used for providing an SC-FDE system model and designing a data frame structure aiming at a nonlinear equilibrium method;

the HDFE equalization algorithm analysis module is used for analyzing the HDFE equalization algorithm based on MMSE;

and the RISI interference estimation module is used for accurately estimating the RISI interference by the HDFE-MMSE-RISIC equalization algorithm by utilizing the judgment result of the HDFE equalizer on the data symbol.

By combining all the technical schemes, the invention has the advantages and positive effects that: the invention provides an HDFE-MMSE-RISIC equalization algorithm on the basis of analysis and research of a time-frequency domain Hybrid Decision Feedback Equalization (HDFE) equalization algorithm based on Minimum Mean Square Error (MMSE) under an SC-FDE system, wherein the improved equalization algorithm utilizes the HDFE equalization algorithm to improve the estimation accuracy of residual intersymbol interference (RISI), and compared with the HDFE equalization algorithm based on MMSE, the influence of intersymbol interference is further eliminated. Through simulation comparison, the HDFE-MMSE-RISIC equalization algorithm can obtain better error code performance under a scattering channel.

The HDFE-MMSE-RISIC equalization algorithm provided by the invention further improves the error code performance of the equalization algorithm by eliminating the influence of residual intersymbol interference on the basis of the HDFE equalizer.

In summary, the HDFE-MMSE-RISIC equalization algorithm provided by the invention improves the accuracy of the decision feedback result in the equalization calculation process through the feedback filtering calculation of the HDFE equalizer, and then utilizes the decision result to estimate RISI, thereby further improving the accuracy of the residual intersymbol interference result. Through simulation analysis under a scattering channel based on an SC-FDE system, the improved equalization algorithm has obvious improvement on error code performance.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings used in the embodiments of the present application will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present application, 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 equalization method of a scatter communication system according to an embodiment of the present invention.

Fig. 2 is a schematic structural diagram of an SC-FDE system model according to an embodiment of the present invention.

Fig. 3 is a schematic diagram of a structure of a transmission data frame according to an embodiment of the present invention.

Fig. 4 is a block diagram of a calculation flow of an HDFE equalizer according to an embodiment of the present invention.

Fig. 5 is a block diagram of the calculation flow of the HDFE-MMSE-RISIC equalization algorithm according to an embodiment of the present invention.

Fig. 6 is a schematic diagram comparing constellation diagrams of MMSE-HDFE equalization and HDFE-MMSE-RISIC equalization results provided by the embodiment of the present invention.

Fig. 7 is a schematic diagram illustrating error performance comparison between HDFE equalization and HDFE-MMSE-RISIC equalizer according to an embodiment of the present invention.

Fig. 8 is a frame structure diagram of two transmitting antennas according to an embodiment of the present invention.

Fig. 9 is a schematic structural diagram of a transmitting end of an HDFE-MMSE-HDFE equalizer in an MIMO system according to an embodiment of the present invention.

Fig. 10 is a schematic structural diagram of a receiving end of an HDFE-MMSE-HDFE equalizer in an MIMO system according to an embodiment of the present invention.

Fig. 11 is a schematic design structure diagram of an equalization algorithm under a system according to an 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 are not intended to limit the invention.

In view of the problems in the prior art, the present invention provides an equalization method, system, computer device, and processing terminal for a scatter communication system, and the present invention is described in detail below with reference to the accompanying drawings.

As shown in fig. 1, an equalization method for a scattering communication system according to an embodiment of the present invention includes:

s101, providing an SC-FDE system model, and designing a data frame structure aiming at a nonlinear equalization algorithm;

s102, analyzing an HDFE equalization algorithm based on MMSE, wherein the algorithm reduces ISI interference through feedforward frequency domain equalization and time domain feedback filtering calculation;

s103, an HDFE-MMSE-RISIC equalization algorithm is provided, the algorithm accurately estimates the interference of RISI by using the decision result of the HDFE equalizer on the data symbol, and the influence caused by intersymbol interference is reduced through the summation compensation of a time domain.

The invention is further described with reference to specific examples.

1. System model

The structural block diagram of the SC-FDE system is shown in fig. 2, the input data is subjected to LDPC coding, constellation mapping and the like to obtain complex transmission data, then a guard interval CP is added and the complex transmission data is transmitted, and after the CP is removed from the data at the receiving end, a frequency domain received signal is obtained through FFT calculation, then channel estimation and frequency domain equalization calculation are performed, and finally the complex transmission data is transformed to a time domain through IFFT, and signal detection is completed through demapping, channel decoding and the like, thereby completing the design of the entire single carrier frequency domain equalization system.

In the nonlinear equalization calculation of the invention, the data transmission symbols are expanded by using the all-zero sequence. As shown in the figure3, to avoid intersymbol interference, the sequence length is defined as L, while ensuring that L ≧ L_hHere L_hIs the Channel Impulse Response (CIR) order. The length of the data transmitted effectively is M-N-L. Finally, in the framing of the transmitted data, a zero sequence with the length of L needs to be added in front of the consecutive data subframes.

Here, the data to be transmitted is denoted as x, and after zero sequence expansion with length L, the transmitted data block s can be expressed as:

where M represents the number of effective transmission data and N is the size of the non-linear equalized transmission data block.

Thus, when the data block passes through the multipath channel h, a received signal after linear convolution can be obtained:

with matrix vector form representation, we can get:

wherein the content of the first and second substances,

since here the second half of the data vector s is an all-zero sequence L long, here the matrix is thenCan be expressed as Λ, then equation (3) can be expressed as:

r＝Λs+w (5)

where Λ is the Toepliz matrix,

the FFT calculation is carried out to the frequency domain, and the following results can be obtained:

R＝HS+W (7)

wherein R ═ F_Nr，S＝F_Ns，W＝F_Nw，R represents DFT transformation of a received signal, S represents a zero-filling expanded frequency domain transformation value of transmitted data x, W represents DFT transformation of a noise vector, and H obtains a calculation result as a diagonal matrix according to Toepliz matrix properties. Here, F_NDenote an nxn DFT matrix with the elements in the matrix expressed as:

then equation (7) is further simplified to yield:

R_k＝H_kS_k+W_k＝H_kX_k+W_k,k＝0,1,...,N-1 (9)

2. improved equalization algorithm

2.1 MMSE-based HDFE equalization

The structure of the HDFE equalizer based on MMSE is shown in fig. 4, the equalizer first obtains a data block with length N through serial-to-parallel conversion, the feedforward frequency domain filter performs linear equalization calculation, then the frequency domain equalization result is converted to the time domain through IFFT, and symbol-by-symbol selection decision is performed, so that the feedback filter can perform time domain filtering calculation by using the decided symbol, and perform cancellation compensation on the equalization result, thereby reducing the inter-symbol interference.

Specifically, the received frequency domain data is equalized by a coefficient G_FFThe feedforward frequency-domain filter of (2) needs to perform complex multiplication to obtain:

Y_k＝R_kG_FF,k,k＝0,1,...N-1 (10)

then, the time domain is converted through IFFT:

frequency selective fading channel H and feed forward filter G_FFMay be equivalently expressed as U-HG_FFIt is assumed here that the equivalent channel response length after concatenation is L_uThen, we can go to the time domain to get:

that is, the time domain value of the transmission data after passing through the feedforward filter is:

wherein the content of the first and second substances,

when L is_uWhen L is less than or equal to L, the above-mentioned cyclic convolution is identical to the linear convolution, so that the signal transmission symbol with effective length M can be obtained, namely:

and the remaining L-length data is an all-zero sequence. It is noted that when L is_uAbove L, the subsequent time-domain feedback filter cannot completely eliminate ISI interference.

Then as shown in the lower part of fig. 4, the input to the feedback filter is the decided symbol and the known all-zero sequence:

the coefficient of the feedback filter is here defined as g_FBAnd has a length L_gThen, the time domain received signal after feedback can be expressed as:

then need to be pairedMaking symbol-by-symbol selection decision to obtain the decided symbol result

In an HDFE equalizer, ZF or MMSE based equalization can be implemented according to different equalization calculation criteria. The implementation complexity of the two is similar, but the ZF equalization performance is inferior to the latter, so that only further research on the MMSE-based HDFE equalization algorithm is made here.

According to the minimum mean square error principle, the cost function here can be expressed as:

when L is_gL is less than or equal to L, and the feedback coefficient is transformed by the frequency domain:

then equation (18) can be obtained by FFT transformation to the frequency domain:

whereinIs the power at which the data symbols are transmitted. Performing a gradient operation on the above equation can obtain:

the substitution formula (20) comprises:

defining feedback coefficientsThen, by performing gradient operation on the above formula and substituting equation (19), the feedback coefficient can be obtained to satisfy:

Ag_FB＝b (23)

wherein the content of the first and second substances,

where A is L_g×L_gOf the Toeplitz matrix, its first behaviorAnd b is L_gColumn vector of length of valueThe IDFT transform of (3), so that the calculation results of both can be obtained by N-point FFT/IFFT here.

Combining the above analysis, the basic steps of the MMSE-based HDFE equalization algorithm are as follows:

1) calculating a time domain feedback coefficient g by using a channel estimation result and a signal-to-noise ratio estimation value (which is a noise power value after normalization) through an equation (23) and an equation (24)_FBThen feedback is obtained by FFTFrequency domain representation G of coefficients_FB。

2) Then, through the calculation of the formula (21), a feedforward frequency domain equalization coefficient value G based on MMSE is obtained_FFAnd the equalization coefficient is used for completing the complex multiplication calculation of the feedforward filter, then the IFFT conversion is carried out on the obtained linear frequency domain calculation result, and the first M effective data lengths are intercepted to obtain the equalization calculation value of the time domain.

3) And (2) calculating a time domain feedback value by using the formula (17) and combining the known all-zero (or decided) sequence with the length of L, and performing symbol-by-symbol summation feedback on the time domain equalization result obtained in the step (2).

4) And 4, carrying out symbol-by-symbol selection judgment on the feedback result obtained in the step 3. Then, a judged symbol sequence after time-domain feedback calculation can be obtained, and then the L-length known sequence obtained in step 3 is shifted by one bit to the left and removed, and symbols after selective judgment are added and replaced, so as to reconstruct the judged symbol sequence with the length of L.

5) And (4) repeating the steps 3 and 4 to perform symbol-by-symbol feedback and decision, and performing M times of circulation until the time domain equalization results obtained in the step 2 all complete feedback filtering calculation to obtain final time-frequency domain mixed decision feedback calculation results.

2.2 improved HDFE-MMSE-RISIC equalization Algorithm

According to the time-frequency domain mixed feedback equalization principle, the calculation of the HDFE equalization algorithm based on MMSE mainly comprises frequency domain feedforward filtering calculation and time domain feedback filtering calculation. Specifically, here the calculation process of the feedforward filter can be expressed as:

then, the frequency domain equalization calculation result is subjected to N-point IFFT, and the M-length effective data is intercepted to obtain:

wherein the content of the first and second substances,the values after the pre-M points are intercepted after the N-point IFFT is performed are shown, that is, the DFT transformation matrix only containing NxNThe first M rows of (a).

Thus, in equation (26), the first part represents the feedforward filtered data portion and the second part represents the residual intersymbol interference value δ after MMSE equalization_nThe third part represents the colored noise due to frequency domain equalization of HDFE

Then g is calculated by the formula (23)_FB,nAnd selecting the decided symbolFiltering calculation is carried out to obtain a time domain feedback value

The HDFE-MMSE-RISIC equalization algorithm provided by the invention aims at the residual intersymbol interference delta of the formula (26)_nWhen the data block with length N is subjected to HDFE equalization calculation, as shown in the dashed-line diagram in the lower half of fig. 5, the equalized data y needs to be estimated for residual interference in the frequency domain_nAnd carrying out selection judgment to obtain effective transmission data with the length of M. Then, expanding the judged sequence to N length by utilizing a known all-zero sequence, carrying out FFT calculation, and then calculating a RISI estimated value of a frequency domain by combining a noise estimated value:

IFFT calculation is performed on equation (28) to obtain a residual intersymbol interference value with a length of N, as shown in FIG. 4, where the effective data with the previous length of M is truncated to obtain an estimated RISI value in the time domainAssuming here that the selection decision result is completely correct, i.e.The equalization results can be obtained as:

according to the principle analysis of the HDFE-MMSE-RISIC equilibrium algorithm, the implementation steps of the improved algorithm can be obtained as follows:

1) firstly, carrying out HDFE equalization based on MMSE to obtain time domain equalization data with the length of N after selection judgment, and then intercepting the effective result of M before.

2) And adding a known all-zero sequence to the HDFE equalization result obtained in the previous step, thereby expanding the data to be N-long again. Then FFT conversion of N points is carried out to obtain the frequency domain value X of the judged symbol_k。

3) Calculating a frequency-domain value of residual intersymbol interference according to equation (28)And obtaining the RISI value of the time domain after IFFT transformation, and intercepting the effective data length of the former M.

4) After the MMSE-based HDFE equalization calculation, a feedback equalization result with the same length as M can be obtained, and the time domain data is subjected to summation compensation operation as shown in the formula (29), so that the calculation of the equalization algorithm of the whole HDFE-MMSE-RISIC is completed.

The technical effects of the present invention will be described in detail with reference to simulation experiments.

1. The invention is based on the SC-FDE system and utilizes the scattering channel model described in the table 1 to simulate the error code performance of the improved equalization algorithm. The equalization data block length N is 1024, the guard interval length L is 32, LDPC channel coding with a coding rate of 2/3 is used, and the modulation scheme is QPSK modulation.

TABLE 1 Scattering communication simulation channel model

The following is a detailed development of the simulation and analysis of the performance of the improved equalization algorithm:

the signal-to-noise ratio E is shown in FIG. 6_b/N₀When the signal is 9dB, comparing the equalization result constellation diagram of the HDFE equalization algorithm based on MMSE with the HDFE-MMSE-RISIC equalization algorithm provided by the invention, it can be obviously seen that the improved equalization algorithm reduces the influence of noise on the received signal by further eliminating the residual intersymbol interference.

2. Further, the error code performance of the MMSE-based HDFE equalization algorithm and the HDFE-MMSE-RISIC equalization algorithm under the scattered communication is compared below. The filter coefficients defining decision feedback equalization are all of length L_g＝12(L_h-1), i.e. close to optimal decision feedback performance.

As can be seen from fig. 7, at a bit error rate of 10^-5Compared with the HDFE equalization algorithm, the improved MMSE-RISIC equalization algorithm based on the HDFE has the error code performance improvement of about 1 dB. In addition, compared with a frequency domain linear MMSE equalization mode, the performance is improved by about 3 dB.

3. Scheme expansion for improved equalization algorithms

3.1 HDFE-MMSE-RISIC equalization under STBC based MIMO system

In an STBC based MIMO system, according to the space-time block coding principle, the two transmit antenna frame structure here is as shown in fig. 8:

wherein, the transmission signals of the two sending antennas are x respectively₁(n) and x₂(N), i and i + N respectively denote oddData transmission block of even time slot, where N is data symbol size. Thus, when the transmit data is denoted as s (n), the two transmit antenna data time domain value may be expressed as:

wherein s is_r＝[s₀ s_N-1 s_N-2 .... s₁]。

Based on the coding and decoding principle of STBC, the HDFE-MMSE-HDFE equalizer in MIMO system is designed, and its structural block diagram is shown in fig. 9 and fig. 10:

in STBC decoding calculation at receiving end, firstly, frequency domain estimation value of channel response is utilizedThen, the data of the two receiving antennas are combined and calculated, so as to obtain the frequency domain decoding result divided into odd and even time slots, and the specific calculation process is as follows:

wherein a value of a channel parameter lambda is definedIs composed ofRepresenting the colored noise frequency domain values of the two receive antennas.

Thus, the design implementation of the equalization algorithm in the SC-FDE system is similar to that in the SC-FDE system, wherein the channel parameter λ corresponds to the above R_k＝H_kX_k+W_kChannel value H in_kTherefore, according to the calculation steps of the HDFE-MMSE-RISIC equalization algorithm in the single-antenna system, the design realization of the improved equalization algorithm under the MIMO system can be realized, and the system capacity and the link reliability under the scattering communication system are improved.

3.2 HDFE-MMSE-RISIC equalization based on OFDM

Compared with the SC-DFE system, if the peak-to-average ratio and the carrier synchronization and phase sensitivity are not considered, the design implementation of the equalization algorithm can be realized by utilizing the OFDM technology, and the signal transmission data frame structure is consistent with the SC-FDE system. The most different of the improved equalization algorithm based on the OFDM system is that the transmitting end needs to perform IFFT processing first, so as to obtain the design structure of the equalization algorithm in the system, as shown in fig. 11.

When the guard interval is greater than the maximum time delay of the channel, similar to a single carrier modulation mode, a frequency domain signal transmission model can be obtained as follows:

Y_k＝H_kX_k+W_k,k＝0,1,...,N-1

then, the complex value after constellation mapping in the upper graph is equivalent frequency domain transmission data, and through OFDM modulation, the data symbol is equivalent to pass through N independent transmission channels, that is, the frequency selective fading channel is divided into N parallel flat fading subcarrier channels. And then, at a receiving end, obtaining the transmitted frequency domain transmission data through FFT calculation processing, then, keeping the equalization calculation steps consistent with those of the SC-FDE system, and combining the HDFE-MMSE-RISIC equalization algorithm principle to realize the design and implementation of an improved equalization algorithm under the OFDM system.

Of course, under the OFDM carrier modulation, each data symbol experiences a channel that is approximately flat fading, and thus, the ISI interference on the data transmitted by the receiving end is small at this time, so the performance of the improved HDFE-MMSE-RISIC equalization algorithm is only improved a little theoretically.

In general, the equalization algorithm provided by the present invention is applicable to different communication technical schemes, including MIMO systems based on STBC, and different carrier modulation modes, i.e. OFDM systems, etc. Therefore, by combining different system implementation schemes, the improved equalization algorithm can effectively improve the error code performance, the communication link reliability and the like of a communication system in the environment with serious frequency selective fading such as a troposphere scattering channel and the like, and has stronger practicability and wider applicability.

It should be noted that the embodiments of the present invention can be realized by hardware, software, or a combination of software and hardware. The hardware portion may be implemented using dedicated logic; the software portions may be stored in a memory and executed by a suitable instruction execution system, such as a microprocessor or specially designed hardware. Those skilled in the art will appreciate that the apparatus and methods described above may be implemented using computer executable instructions and/or embodied in processor control code, such code provided on a carrier medium such as a diskette, CD-or DVD-ROM, a programmable memory such as read-only memory (firmware), or a data carrier such as an optical or electronic signal carrier, for example. The apparatus and its modules of the present invention may be implemented by hardware circuits such as very large scale integrated circuits or gate arrays, semiconductors such as logic chips, transistors, or programmable hardware devices such as field programmable gate arrays, programmable logic devices, etc., or by software executed by various types of processors, or by a combination of hardware circuits and software, e.g., firmware.

The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any modification, equivalent replacement, and improvement made by those skilled in the art within the technical scope of the present invention disclosed in the present invention should be covered within the scope of the present invention.