阅读说明：本技术 正交时频空系统的信号接收方法、接收装置及存储介质 (Signal receiving method, receiving device and storage medium of orthogonal time-frequency space system ) 是由 别志松 靳宸茜 许文俊 高晖 林雪红 于 2019-09-03 设计创作，主要内容包括：本发明公开了一种正交时频空系统的信号接收方法,包括：对接收信号进行正交频分复用(OFDM)解调,得到时频域接收信号；对所述时频域的接收信号进行时频域线性均衡；其中,在时频域线性均衡时使用带状矩阵代替频域信道响应矩阵；将经过时频域均衡的时频域接收信号变换至时延-多普勒域；对时延-多普勒域的接收信号进行时延-多普勒域均衡；以及输出经过时延-多普勒域均衡的时延-多普勒域的接收信号。本发明还公开了接收装置以及计算机可读存储介质。(The invention discloses a signal receiving method of an orthogonal time-frequency space system, which comprises the following steps: carrying out Orthogonal Frequency Division Multiplexing (OFDM) demodulation on the received signal to obtain a time-frequency domain received signal; performing time-frequency domain linear equalization on the received signals of the time-frequency domain; wherein, a banded matrix is used for replacing a frequency domain channel response matrix in time-frequency domain linear equalization; transforming the time-frequency domain received signal after the time-frequency domain equalization to a time delay-Doppler domain; carrying out delay-Doppler domain equalization on a received signal in a delay-Doppler domain; and outputting the delay-Doppler domain equalized delay-Doppler domain received signal. The invention also discloses a receiving device and a computer readable storage medium.)

1. A signal receiving method of an orthogonal time-frequency space system is characterized by comprising the following steps:

carrying out Orthogonal Frequency Division Multiplexing (OFDM) demodulation on the received signal to obtain a time-frequency domain received signal;

performing time-frequency domain linear equalization on the received signals of the time-frequency domain; wherein, a banded matrix is used for replacing a frequency domain channel response matrix in time-frequency domain linear equalization;

transforming the time-frequency domain received signal after the time-frequency domain equalization to a time delay-Doppler domain;

carrying out delay-Doppler domain equalization on a received signal in a delay-Doppler domain; and

and outputting the received signal of the delay-Doppler domain after the delay-Doppler domain equalization.

2. The signal receiving method as claimed in claim 1, wherein the time-frequency domain linear equalization comprises:

taking N M-dimensional non-zero sub-blocks on the main diagonal of the frequency domain channel response matrix as N M-dimensional first sub-matrices, wherein N is the number of OFDM symbols; m is the number of sub-carriers contained in one OFDM symbol;

replacing the first sub-matrix with a strip matrix; and

and performing linear equalization on the time-frequency domain received signal based on the strip matrix.

3. The signal receiving method according to claim 2, wherein the strip matrix retains a strip region with a main diagonal width 2Q +1 of the first sub-matrix and elements in a triangle with a right angle side length Q at upper left and lower right corners, and the elements at the rest positions are all 0 values, where Q is a preset width parameter.

4. The signal receiving method according to claim 2, wherein the time-frequency domain linear equalization further comprises:

splitting the banded matrix into M (2Q +1) × (4Q +1) second sub-matrices; wherein Q is a preset width parameter; and

the performing linear equalization on the time-frequency domain received signal comprises: and performing linear equalization on the time-frequency domain received signal based on the second sub-matrix.

5. The signal receiving method according to claim 2 or 4, wherein the linear equalization comprises a least mean square error equalization or a zero-forcing equalization.

6. The signal receiving method according to claim 4, wherein the splitting comprises: splitting the strip matrix according to the following formula:

wherein the content of the first and second substances,is the element of the n row and m column position of the strip matrix.

7. A signal receiving apparatus, characterized in that the apparatus comprises:

the OFDM demodulator is used for carrying out OFDM demodulation on the received signal to obtain a time-frequency domain received signal;

the time-frequency domain equalizer is used for performing time-frequency domain linear equalization on the received signals of the time-frequency domain; wherein, a banded matrix is used for replacing a frequency domain channel response matrix in time-frequency domain linear equalization;

the SFFT receiving window is used for converting the time-frequency domain receiving signals subjected to time-frequency domain equalization into a time delay-Doppler domain; and

and the delay-Doppler domain equalizer is used for performing delay-Doppler domain equalization on the received signal in the delay-Doppler domain and outputting the received signal in the delay-Doppler domain subjected to delay-Doppler domain equalization.

8. The receiving apparatus as claimed in claim 7, wherein the time-frequency domain equalizer comprises:

a splitting module, configured to split the frequency domain channel response matrix into N M × M first sub-matrices, where M is the number of subcarriers and N is the number of OFDM symbols;

a replacement module to replace the first sub-matrix with a strip matrix; and

and the equalization module is used for performing linear equalization on the time-frequency domain received signals based on the banded matrix.

9. The receiving apparatus as claimed in claim 8, wherein the time-frequency domain equalizer further comprises:

a second splitting module, configured to split the strip matrix into M (2Q +1) × (4Q +1) second sub-matrices; wherein Q is a preset width parameter; and

the equalization module is used for performing linear equalization on the time-frequency domain received signal based on the second sub-matrix.

10. A computer-readable medium, on which a computer program is stored which, when being executed by a processor, carries out the signal receiving method as claimed in any one of claims 1 to 4 or claim 6.

Technical Field

The present invention relates to mobile communication technologies, and in particular, to a signal receiving method and device for an orthogonal time-frequency-space system, and a computer-readable storage medium.

Background

Orthogonal Frequency Division Multiplexing (OFDM) is one of the most widely used communication technologies, and is mainly used for resisting multipath effects causing inter-symbol interference and simultaneously realizing high-rate data transmission.

However, future wireless communication networks, such as the new generation broadband wireless communication (5G/B5G), will face highly dynamic communication channel environments, such as in high mobility scenarios (e.g., high-speed rail) and millimeter wave (mmWave) communication. The highly dynamic channel exhibits double dispersion properties including time dispersion due to multipath and frequency dispersion due to doppler broadening. OFDM modulation used by current communication systems can be used to combat Inter Symbol Interference (ISI) due to time dispersion. However, Inter-carrier interference (ICI) caused by frequency dispersion can significantly impair the performance of the OFDM system.

In this case, an Orthogonal Time Frequency Space (OTFS) modulation technique works. Specifically, the OTFS system carries transmission information in a delay-doppler domain, and at a transmitting end, spreads each delay-doppler domain information symbol to a whole time-frequency domain within a certain range through Inverse Symplectic Finite Fourier Transform (ISFFT), thereby ensuring that each symbol in an OTFS frame experiences a relatively stable channel, that is, a double-dispersion channel is converted into a channel almost free of frequency/time dispersion. At the receiving end, the OTFS system converts the information of the time-frequency domain to the delay-doppler domain by using Symplectic Finite Fourier Transform (ISFFT) to perform operations such as equalization demodulation. Therefore, the OTFS, as a novel multi-carrier modulation technique, can effectively combat a highly dynamic communication channel environment, and exhibits strong robustness to high doppler spread.

For the OTFS system, the complexity of the equalizer at the receiving end is very important for its practical application. The existing OTFS system equalization techniques are mainly classified into two types: the first is a nonlinear equalization method, which has good bit error rate performance but high complexity, so the practicability is low; the second category is linear equalization, but conventional linear equalization techniques generally involve matrix inversion operations, and the complexity of matrix inversion is unacceptable for large data dimensions in real-world communications.

Disclosure of Invention

In view of this, an embodiment of the present invention provides a signal receiving method for an orthogonal time-frequency space system. The method can comprise the following steps:

carrying out OFDM demodulation on the received signal to obtain a time-frequency domain received signal;

performing time-frequency domain linear equalization on the received signals of the time-frequency domain; wherein, a banded matrix is used for replacing a frequency domain channel response matrix in time-frequency domain linear equalization;

transforming the time-frequency domain received signal after the time-frequency domain equalization to a time delay-Doppler domain;

carrying out delay-Doppler domain equalization on a received signal in a delay-Doppler domain; and

and outputting the received signal of the delay-Doppler domain after the delay-Doppler domain equalization.

Wherein, the time-frequency domain linear equalization includes: taking N M-dimensional non-zero sub-blocks on the main diagonal of the frequency domain channel response matrix as N M-dimensional first sub-matrices, wherein N is the number of OFDM symbols; m is the number of sub-carriers contained in one OFDM symbol; replacing the first sub-matrix with a strip matrix; and performing linear equalization on the time-frequency domain received signal based on the strip matrix.

The strip matrix reserves a strip area with the width of the main diagonal of the first sub-matrix being 2Q +1 and triangle inner elements with the length of a right angle being Q at the upper left corner and the lower right corner, and the elements at the rest positions are all 0 values, wherein Q is a preset width parameter.

Wherein the time-frequency domain linear equalization further comprises: splitting the banded matrix into M (2Q +1) × (4Q +1) second sub-matrices; wherein Q is a preset width parameter; and said linearly equalizing said time-frequency domain received signal comprises: and performing linear equalization on the time-frequency domain received signal based on the second sub-matrix.

The linear equalization includes a minimum mean square error equalization or a zero-forcing equalization.

The splitting comprises the following steps: splitting the strip matrix according to the following formula:

wherein the content of the first and second substances,

is the element of the n row and m column position of the strip matrix.

The embodiment of the invention also provides a receiving device of the orthogonal time-frequency space system, which comprises:

the OFDM demodulator is used for carrying out OFDM demodulation on the received signal to obtain a time-frequency domain received signal;

the time-frequency domain equalizer is used for performing time-frequency domain linear equalization on the received signals of the time-frequency domain; wherein, a banded matrix is used for replacing a frequency domain channel response matrix in time-frequency domain linear equalization;

the SFFT receiving window is used for converting the time-frequency domain receiving signals subjected to time-frequency domain equalization into a time delay-Doppler domain; and

and the delay-Doppler domain equalizer is used for performing delay-Doppler domain equalization on the received signal in the delay-Doppler domain and outputting the received signal in the delay-Doppler domain subjected to delay-Doppler domain equalization.

Wherein, the time-frequency domain equalizer comprises:

a splitting module, configured to split the frequency domain channel response matrix into N M × M first sub-matrices, where M is the number of subcarriers and N is the number of OFDM symbols;

a replacement module to replace the first sub-matrix with a strip matrix; and

and the equalization module is used for performing linear equalization on the time-frequency domain received signals based on the banded matrix.

The time-frequency domain equalizer further comprises:

a second splitting module, configured to split the strip matrix into M (2Q +1) × (4Q +1) second sub-matrices; wherein Q is a preset width parameter; and

the equalization module is used for performing linear equalization on the time-frequency domain received signal based on the second sub-matrix.

An embodiment of the present invention proposes a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the above-described signal receiving method.

The embodiment of the invention adopts a two-stage equalization technology to equalize the received signals in a time-frequency domain and a time delay-Doppler domain respectively, thereby obtaining additional diversity gain. Specifically, in the time-frequency domain equalization, a simplified low-complexity linear equalization algorithm is adopted, and the complexity required by equalization inversion is reduced. And the equalization technology adopted in the time delay-Doppler domain can eliminate ISI caused by the equalization error of the previous stage, and simultaneously effectively utilizes diversity gain by means of the structure similar to a DFE equalizer, thereby greatly improving the performance of the receiver of the OTFS system.

Drawings

FIG. 1 is a schematic diagram of the internal structure of an OTFS system 100 according to some embodiments of the present invention;

fig. 2 is a schematic diagram of an internal structure of an OTFS receiving end 200 according to some embodiments of the present invention;

FIG. 3 shows the Doppler frequency f at normalized_dFrequency domain channel response matrix in case of 0.1

A power profile of;

FIG. 4 shows a schematic view of a ribbon matrix according to an embodiment of the invention;

fig. 5 is a schematic flow chart of a signal receiving method of an orthogonal time-frequency space system according to some embodiments of the present invention;

fig. 6 is a schematic diagram illustrating performance comparison between two-stage equalization and a conventional linear equalization method adopted in the signal receiving method according to the embodiment of the present invention; and

fig. 7 is a schematic diagram of an internal structure of a computing device for implementing signal reception 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 described in further detail below with reference to the accompanying drawings in conjunction with specific embodiments.

As described above, OTFS, as a novel multi-carrier modulation technique, can effectively combat a highly dynamic communication channel environment, and exhibits strong robustness to high doppler spread.

Fig. 1 is a schematic diagram of an internal structure of an OTFS system 100 according to some embodiments of the present invention. As shown in fig. 1, the OTFS system 100 includes: an OTFS sender 101 and an OTFS receiver 102. Wherein, signals are transmitted between the OTFS sending end 101 and the OTFS receiving end 102 through the channel 103. It should be noted that, in the embodiment of the present invention, the channel 103 may be a linear time-varying (LTV) channel.

As shown in fig. 1, the OTFS sender 101 may include: an ISFFT transmission window 1011 and an OFDM modulator 1012. The inside of the OTFS receiving end 102 may include: OFDM demodulator 1021, SFFT receive window 1022, and equalizer 1033.

As can be seen from fig. 1, the model of the OTFS system 100 is similar to the conventional OFDM system model, and may be regarded as that, on the basis of OFDM, a conversion pre-processing from the time-domain to the time-domain is newly added at the transmitting end, and a conversion post-processing from the time-domain to the time-domain and a signal equalization post-processing are newly added at the receiving end.

Specifically, in some embodiments of the present invention, the OTFS sending end 101 may divide the delay-doppler domain in advance. In particular, the delay-doppler domain may be divided into M × N delay-doppler bins Γ, wherein,

where NT is the duration and Δ f is the subcarrier spacing. Next, the OTFS transmitter 101 may transmit a modulation sequence S [ l, k ] at the lattice point Γ]It is used. Here, N may be the number of OFDM symbols included in the OTFS system; m may be the number of subcarriers included in one OFDM symbol.

Specifically, the ISFFT transmission window 1011 of the OTFS transmitting end 101 may convert the modulation sequence S [ l, k ] from the delay-doppler domain to the time-frequency domain through ISFFT transformation, which is denoted as X [ M, N ], where M is 0,1, …, M-1, N is 0,1, …, N-1. In some embodiments of the present invention, the time-frequency domain modulation sequence X [ m, n ] can be obtained by the following formula (1).

As can be seen from the above equation (1), SFFT^-1The Transform is actually an M-point Discrete Fourier Transform (DFT) performed on the matrix S, and an N-point Inverse DFT (IDFT) performed after the Transform. For ease of understanding, equation (1) may be rewritten as the form of equation (2) by matrix multiplication as follows.

Wherein, F_MIs a DFT transformation matrix of M points and is a unitary matrix, namely satisfies

Further, the above formula (2) can be further rewritten as the following formula (3) according to the definition of the Kronecker product.

Wherein S and X are the results of vectorization representation of the delay-Doppler domain symbol S and the time-frequency domain symbol X respectively,

in general, before sending to the OFDM modulator 1012 for processing, the OTFS sender 101 further needs to perform windowing on each column of elements in the time-frequency domain symbol matrix X (the above-mentioned elements are not shown in fig. 1)Windowing). In embodiments of the present invention, the windowed processing result may be used

Indicating, windowing the results

Can be expressed by the following formula (4):

wherein U is a matrix representation of the transmit window function. In some embodiments of the present invention, the transmission window function U may be selected as an ideal window, a rectangular window, or the like.

Up to this point, the preprocessing part of the delay-doppler domain to the time-frequency domain has been completely finished, and the result of the preprocessing can be expressed as the following formula (5) according to the above formulas (3) and (4).

Next, the OFDM modulator 1012 of the OTFS transmitter 101 will perform the preprocessing on the result

And carrying out OFDM modulation. Specifically, the above-mentioned OFDM modulator 1012 will

The OFDM symbols are regarded as N OFDM symbols with M subcarriers, and each symbol is subjected to M-point IDFT to form a time domain OFDM symbol. Next, a Cyclic Prefix (CP) may be added before each symbol to eliminate inter-symbol interference, as processed by the OFDM system. The OFDM symbols are then transmitted over a wireless channel 103.

In particular, in some embodiments of the present invention, the kth sample of the nth OFDM symbol may be represented as

The time domain OFDM symbol can be represented by equation (6) as follows:

wherein the content of the first and second substances,

m sub-carrier representing the nth OFDM symbol, M being the number of DFT points, M_cpThe length of the CP is greater than or equal to the maximum multipath delay of the channel,

in some embodiments of the present invention, assuming that the time domain impulse response of the channel is h [ n, l ], the signal obtained after removing the CP by the OFDM demodulator 1021 of the OTFS receiving end 102 is as shown in the following formula (7):

where L is the number of multipaths, h^k[n,l]Time-varying channel time-domain impulse response, w [ n, k ], for the kth sample of an OFDM symbol]Is a mean value of 0 and a variance of σ²Gaussian white noise sample.

Further using A_cpExpressed as a matrix of symbols plus CP, R_cpRepresenting the matrix with CP removed, H represents the time domain response matrix of the channel, and thus, the above formula (7) can be expressed in the form of a matrix vector as shown in the following formula (8):

next, OFDM demodulator 1021 may pair

The time-frequency domain signal can be obtained by DFT conversionThat is, the input/output relation of the OFDM demodulator 1021 can be expressed in a matrix form as shown in the following equation (9):

wherein the content of the first and second substances,

is a frequency domain channel response matrix of the channel, which is a block diagonal matrix, and can be specifically represented by the following formula (10):

wherein the content of the first and second substances,

a frequency response matrix for the nth OFDM symbol.

After the OFDM demodulator 1021 performs demodulation processing, the SFFT reception window 1022 further performs SFFT processing on the received signal. Specifically, in some embodiments of the present invention, the SFFT processing may specifically include a reception window and an SFFT transform processing. This process is equivalent to the inverse of the processing of the ISFFT transmission window 1011.

In some embodiments of the present invention, the received signal in the delay-doppler domain is denoted as r, and the received signal r can be represented by the following equation (11):

where V is a matrix representation of the receive window function.

Further, according to the above equations (4), (9) and (11), it can be derived that the input-output model of the OTFS system will be shown by the following equation (12):

wherein the content of the first and second substances,

representing a gaussian white noise vector in the delay-doppler domain.

The above equation (12) can be expressed in a matrix form as the following equation (13):

thus, a channel response matrix H of the OTFS system in the delay-Doppler domain is obtained_2DWherein, in the step (A),

the calculation derivation process lays a foundation for the subsequent equalization recovery of the transmitted signal.

In consideration of practical feasibility, the transmission window U and the reception window V may be selected as rectangular windows, i.e., another U-V-I_NMAt this time, the channel response matrix of the delay-doppler domain can be further simplified into the form of the following equation (14):

since ICI caused by high-speed moving scenes reduces the signal-to-noise ratio of the received signal, an equalizer 1033 is added at the receiving end to equalize the received signal.

The following describes a conventional method for performing equalization compensation on a received signal of an OTFS system by taking Minimum Mean Square Error (MMSE) linear equalization as an example. In general, the equalizer 1033 needs to assume that the channel state information is known at present when performing equalization compensation.

If MMSE linear equalization is employed, assuming that the coefficient matrix of MMSE equalization is G, the equalization result of equalizer 1033 can be represented by the following equation (15) according to the above equations (13) and (14):

according to the requirements of MMSE algorithm, it is necessary to minimize the mean square error between the transmitted signals s and Gr, i.e. to make E { ee }^H}＝E{(Gr-s)(Gr-s)^HIs smallest.

In this case, assuming that the received signal r is uncorrelated with the error e, the above assumption can be converted to the requirement that the following equation (16) holds, thereby satisfying the MMSE requirement.

E{er^H}＝0 (16)

Substituting the above formula (16) into the formula (15) gives the following formula (17)

G·E{rr^H}-E{sr^H}＝0 (17)

Further, the following equations (18) and (19) can be obtained from equation (13).

Thus, the coefficient matrix G of the delay-doppler domain MMSE linear equalization can be represented by the following equation (20):

also, the result of MMSE linear equalization can be represented by the following equation (21):

as can be seen from the above equation (21), the main complexity of MMSE linear equalizer in the delay-doppler domain of OTFS system is determined by the matrix inversion calculation of NM × NM, and this equalization method is hardly applicable once the dimension of NM is increased.

In order to reduce the complexity of the equalization processing, the embodiment of the invention provides a scheme of secondary equalization, and specifically, a receiving end performs primary equalization in a time-frequency domain at first and then performs primary equalization in a time delay-doppler domain. The complexity of the equalization processing can be effectively reduced through the two equalization processing, and the equalization effect can be improved.

Fig. 2 shows a structure of an OTFS receiving end 200 according to an embodiment of the present invention. As shown in fig. 2, an OTFS receiving end 200 according to an embodiment of the present invention includes: an OFDM demodulator 201, a time-frequency domain equalizer 202, an SFFT receive window 203, and a delay-doppler domain equalizer 204.

The OFDM demodulator 201 may perform OFDM demodulation on a signal received from a channel to obtain a time-frequency domain received signal.

For the specific operation process of the OFDM demodulator 201, reference may be made to the processing process of the OFDM demodulator 1021, which is not described herein again.

The time-frequency domain equalizer 202 may perform time-frequency domain equalization processing on the received signal in the time-frequency domain. The specific processing procedure will be described in detail below.

The SFFT receiving window 203 may perform receiving window and SFFT transform processing on the time-frequency domain received signal subjected to time-frequency domain equalization to obtain a delay-doppler domain received signal.

For the specific operation process of the SFFT receive window 203, reference may be made to the processing process of the SFFT receive window 1022, which is not described herein again.

The delay-doppler domain equalizer 204 may perform a delay-doppler domain equalization process on the delay-doppler domain received signal. The specific processing procedure will also be described in detail below.

The following describes in detail the procedure of the time-frequency domain equalization processing performed by the time-frequency domain equalizer 202 according to the embodiment of the present invention.

In an embodiment of the present invention, the frequency domain channel response matrix H of the channel shown according to the above equation (10)_fStructure of (1), H_fIs a block diagonal matrix, and because there are N nonzero sub-blocks with M × M dimensions on the main diagonal, each nonzero sub-block can be used as a small matrixN small matrixes H with M multiplied by M dimensions can be obtained_f,nWherein N is more than or equal to 0 and less than or equal to N-1. In this case, the above equation (9), that is, the input/output relational expression of the OFDM demodulator 201, may be rewritten into the following equation (22), and the following equation (23) may be further obtained:

y_n＝H_f,nx_n+w_n,0≤n≤N-1 (23)

y in the above formula (23)_n＝[y_N×n,y_N×n+1,…,y_N×n+M-1]^T，x_n＝[x_N×n,x_N×n+1,…,x_N×n+M-1]^T. Therefore, in the embodiment of the present invention, the above equation (23) can be regarded as an input-output model of the nth OFDM symbol in the time-frequency domain.

Next, the time-frequency domain equalization algorithm for the nth OFDM symbol is derived.

For the convenience of derivation, in the embodiment of the present invention, the frequency domain channel response matrix H of the nth OFDM symbol is used_f,nReferred to as the first sub-matrix, denotedWill y_nB, denoting x_nIs denoted as a, w_nDenoted by v, the above formula (23) can be expressed as the following formula (24):

FIG. 3 shows the Doppler frequency f at normalized_dFrequency domain channel response matrix in case of 0.1

Where N-M-32 is assumed.

As can be seen from fig. 3, the matrix is due to the effect of ICI caused by the fast time-varying channelIt will no longer be a diagonal matrix and therefore a conventional single tap equalizer does not apply to this case. The ICI of the channel has the characteristic of sparse distribution, and the influence on the subcarrier signal is mainly caused by a limited number of adjacent subcarriers in the channel. In light of this characteristic, in embodiments of the present invention, matrices having a banded structure (referred to as banded matrices for short) may be used to approximate the actual frequency domain channel response matrixTo simplify the complexity of equalization.

Therefore, in the embodiment of the present invention, the channel response matrix H is generated in the frequency domain of the channel_fSplitting into N MxM dimensional first sub-matrices

Thereafter, the time-frequency domain equalizer 202 may further use a strip matrix

Instead of the first sub-matrix described above

Fig. 4 shows a schematic view of a ribbon matrix according to an embodiment of the invention. As can be seen from FIG. 4, the approximated band matrix

And

the power distribution profiles are very similar, in the form of a strip matrix

The first sub-matrix is reserved

Main diagonal width 2Q +1 banded regionAnd elements in the fields and the triangular regions with the right-angle side length of Q at the positions of the lower left corner and the upper right corner, and the elements at the other positions are all 0 values. Wherein Q is a preset width parameter, and the larger the Q value is, the closer the strip matrix is to the first sub-matrix

Next, based on the above-mentioned band matrix

The time-frequency domain equalizer 202 may equalize the time-frequency domain received signal output by the OFDM demodulator 201 by a linear equalization method.

Specifically, taking MMSE linear equalization as an example, the input-output relation of OFDM demodulator 201 shown in the following equation (25) can be obtained according to the above equation (24):

based on the above equation (25), it can be obtained that the result of MMSE equalization on the nth OFDM symbol is expressed by the following equation (26):

based on the above equation (25), it can be obtained that the result of zero-forcing equalization on the nth OFDM symbol is expressed by the following equation (27):

therefore, compared with the MMSE equalization in the delay-Doppler domain, the frequency domain channel response matrix of the nth OFDM symbol is replaced by the strip matrix in the time-frequency domain equalization process, so that the calculation complexity of channel matrix inversion can be greatly reduced, and the equalization processing complexity is reduced.

In the embodiment of the present invention, in orderTo further simplify the complexity of the time-frequency domain equalization, the time-frequency domain equalizer 202 may further apply the banded matrixSplitting into M (2Q +1) × (4Q +1) second sub-matrices

In the embodiment of the present invention, the splitting may be specifically obtained by the following formula (28):

wherein the content of the first and second substances,

is the element of the n row and m column position of the strip matrix.

The second sub-matrix is constructed according to the formula (24)

Associated channel information equation (29):

wherein, in the above channel information formula (29), b_k＝[b_{(k-Q)mod M},…,b_{(k+Q)mod M}]^T，a_k＝[a_{(k-2Q)mod M},…,a_{(k+2Q)mod M}]^T。

For example, in the case where Q is 1, and a band matrix is assumed

It can be expressed as a structure shown in the following equation (30):

at this time, if the first subcarrier a is desired₀Is constructed according to the splitting method shown in the above formula (28)

The following formula (31) can be obtained

And can obtain the related channel information formula (32)

Next, based on the second sub-matrix

The time-frequency domain equalizer 202 may equalize the time-frequency domain received signal output by the OFDM demodulator 201 by a linear equalization method.

Specifically, when MMSE equalization is adopted, the equalization result of the k-th subcarrier of the nth OFDM symbol can be represented by the following formula (33):

in the above formula (33), e_2QThe column vector is a column vector with length of 4Q +1, wherein the 2Q element is 1, and the rest elements are 0.

Specifically, when ZF equalization is employed, the equalization result of the k-th subcarrier of the nth OFDM symbol can be represented by the following formula (34):

in the above formula (34), e_2QThe column vector is a column vector with length of 4Q +1, wherein the 2Q element is 1, and the rest elements are 0.

Based on the above description, in the embodiment of the present invention, the time-frequency domain equalizer 202 may include:

the device comprises a splitting module, a frequency domain channel response matrix generating module and a frequency domain channel response matrix generating module, wherein the splitting module is used for splitting the frequency domain channel response matrix into N first sub-matrixes with M multiplied by M dimensions, M is the number of subcarriers, and N is the number of OFDM symbols;

a replacement module for replacing the first sub-matrix with a strip matrix; and

and the equalization module is used for performing linear equalization on the time-frequency domain received signals based on the banded matrix.

The time-frequency domain equalizer 202 may further include: a second splitting module, configured to split the strip matrix into M (2Q +1) × (4Q +1) second sub-matrices. In this case, the equalization module is configured to perform linear equalization on the time-frequency domain received signal based on the second sub-matrix.

It can be seen that, with the above method, only the second sub-matrix of (2Q +1) × (4Q +1) in the above formula (27) needs to be inverted, and the computation complexity is much less than that of directly inverting the channel matrix of NM × NM, and is also less than that of inverting the first sub-matrix of M × M or the band matrix.

In the embodiment of the present invention, because the strip matrix only retains the elements of a part of the original frequency domain channel response matrix, the performance of the algorithm is related to the Q of the data width of the retained matrix, and the larger Q means that the more channel information is retained, the better the performance of the algorithm is, but the complexity is increased. Generally, in practical application, Q ≧ f is selected_d]+1 can approach the performance of MMSE equalization algorithm that directly inverts the NM × NM matrix.

In the embodiment of the present invention, the equalization result of the time-frequency domain equalization algorithm is recorded as

The SFFT receive window 203 will equalize the result

Mapping to the delay-doppler domain, a delay-doppler domain signal can be obtained as shown in the following formula (35):

due to the effect of the SFFT transform, the equalization error in the time-frequency domain will be transformed into ISI in the delay-doppler domain, and therefore an equalization algorithm needs to be applied in the delay-doppler domain to remove the ISI.

The processing of the delay-doppler domain equalizer 204 according to the embodiment of the present invention will be described in detail below.

In an embodiment of the present invention, the delay-doppler domain equalizer 204 may further pre-process the delay-doppler domain signal output by the SFFT receive window 203 before performing the delay-doppler domain equalization. Specifically, the preprocessing step may include:

first, an autocorrelation matrix is obtained for the delay-doppler domain channel response matrix obtained by the above equation (14)

Then, the simplified autocorrelation matrix is obtained by simplifying the autocorrelation matrix

In particular, the simplification process may include removing the autocorrelation matrix

The main diagonal and the value where the power is smaller. Here, less may mean a value of power less than a predetermined threshold, e.g., 1 e-4.

After the preprocessing, the delay-doppler domain equalization algorithm can be expressed by the following equation (36):

wherein the content of the first and second substances,

represent will

And (5) performing a hard decision result. The hard decision specifically includes: finger holderMapping to the constellation point closest thereto.

According to the above equation (36), it is assumed that the hard decision symbols are all correct, i.e., that

The delay-doppler domain equalizer output can be represented by the following equation (37):

as can be seen from the above equation (37), ISI is completely removed by the above equalization method, and only useful signal and additive noise are present in the output signal.

As such, the delay-doppler domain equalizer 204 may include:

an autocorrelation module for calculating an autocorrelation matrix for a channel response matrix of the delay-Doppler domain

A simplification module for simplifying the autocorrelation matrix to obtain a simplified autocorrelation matrix

In particular, the simplification process may include removing the autocorrelation matrix

The main diagonal and the value where the power is smaller.

A delay-Doppler domain equalization module for equalizing the signal based on the simplified autocorrelation matrix

And determining a delay-Doppler domain equalization result.

Specifically, the delay-doppler domain equalization method can be represented by the above equation (36).

Looking at equation (37) above, the result is similar to the output result of applying Decision Feedback Equalization (DFE), so the delay-doppler domain equalization can also be considered as a DFE equalizer. The use of the equalizer described above results in improved performance.

Based on the above research results, an embodiment of the present invention provides a signal receiving method for an orthogonal time-frequency space system, where a specific implementation process of the method is shown in fig. 5, and mainly includes:

step 501: and carrying out OFDM demodulation on the received signal to obtain a time-frequency domain received signal.

In the embodiment of the present invention, the OFDM demodulation may include removing a CP from a received signal, and performing DFT on the received signal.

Step 502: and performing time-frequency domain linear equalization on the received signals of the time-frequency domain. Wherein, a banded matrix is used to replace a frequency domain channel response matrix in time-frequency domain linear equalization.

In an embodiment of the present invention, the time-frequency domain linear equalization specifically may include:

step 5021: the frequency domain channel response matrix H_fTaking N M-dimensional non-zero sub-blocks on the main diagonal as N M-dimensional first sub-matrixes

Step 5022: using a ribbon matrix

Instead of the first sub-matrix described above

Wherein the strip matrixThe first sub-matrix is reserved

The main diagonal width is 2Q +1, and the right-angle side length at the upper left and lower right corners is Q, and the other elements are all approximate to 0.

Step 5023: based on the strip matrix

And performing linear equalization on the time-frequency domain received signal.

In an embodiment of the present invention, the linear equalization may be MMSE equalization or zero-forcing (ZF) equalization.

In the embodiment of the present invention, in order to further reduce the complexity of time-frequency domain equalization, before performing step 5023, the following steps may be further performed: the strip matrix is formed

Splitting into M (2Q +1) × (4Q +1) second sub-matricesIn this case, in step 5023, it will be based on the second sub-matrix described above

And equalizing the time-frequency domain received signal in a linear equalization mode.

Specifically, in the embodiment of the present invention, the band matrix may be formed in a manner shown by the above-described formula (28)

Splitting into M (2Q +1) × (4Q +1) second sub-matrices

When the MMSE equalization method is employed, the equalization result can refer to the above equation (33). When ZF equalization method is adopted, the equalization result can be obtained by referring to the above equation (34)

Step 503: and transforming the time-frequency domain received signal subjected to the time-frequency domain equalization into a time delay-Doppler domain.

Specifically, in this step, SFFT conversion processing may be performed on the time-frequency domain received signal subjected to time-frequency domain equalization, so as to obtain a delay-doppler domain received signal.

In the embodiment of the present invention, before performing step 503, the time-frequency domain equalized time-frequency domain received signal may be subjected to a receive window process. Wherein, the receiving window may be a rectangular window.

Step 504: and carrying out delay-Doppler domain equalization on the received signals in the delay-Doppler domain.

In an embodiment of the present invention, the delay-doppler domain linear equalization specifically may include:

step 5041, obtaining autocorrelation matrix of the channel response matrix in the time delay-Doppler domain

Step 5042, simplifying the autocorrelation matrix to obtain a simplified autocorrelation matrix

In particular, the simplification process may include removing the autocorrelation matrix

The main diagonal and the value where the power is smaller.

Step 5043, simplifying the autocorrelation matrix based on the above

And determining a delay-Doppler domain equalization result. The specific equalization result can refer to the above equation (36).

Step 505: and outputting the received signal of the delay-Doppler domain after the delay-Doppler domain equalization.

By combining the two-stage equalization operations, it can be seen that the embodiment of the invention can greatly improve the performance of the receiver on the premise of not increasing the complexity of the traditional delay-doppler domain linear MMSE equalization algorithm.

Firstly, the scheme of the second-stage equalization provided by the invention is different from the scheme of performing equalization in a time-frequency domain or a time delay-Doppler domain, firstly performs time-frequency domain equalization to eliminate ICI caused by Doppler spread, and then performs equalization in the time delay-Doppler domain to eliminate ISI caused by the residual ICI of the previous-stage equalization. The time-frequency domain equalization adopts low-complexity linear equalization, and the complexity of a receiver is reduced. The subsequent delay-Doppler domain equalization method is similar to a DFE, so that the ISI is eliminated, and meanwhile, the diversity gain can be effectively utilized, and the performance is improved.

For the second-level equalization algorithm provided by the embodiment of the invention, the performance of the algorithm in a high-speed moving environment can be simulated through the following experiment, and the result shows that the algorithm can obtain extra grading gain and the performance is improved.

The simulation parameter settings used in this experiment are shown in table 1 below:

parameter name	Parameter value	Parameter name	Parameter value
				N	32	Subcarrier spacing	30KHZ
M	32	CP Length	13
				Modulation system	4QAM	Channel model	TU6

TABLE 1

To compare the Bit Error Rate (BER) performance of the new algorithm with the conventional equalization, the normalized Doppler frequency f was simulated_dIn the case of 0.1, the BER of the two-stage equalization method according to the embodiment of the present invention is compared with the BER and the signal-to-noise ratio (SNR) of the conventional MMSE equalization, and accurate channel information is assumed in the simulation, and Q is 2.

Fig. 6 shows simulation results of the experiment according to the embodiment of the present invention. As can be seen from the simulation results shown in fig. 6, the BER performance of the two-stage equalization method of the present invention is greatly improved compared to the conventional one-stage MMSE equalization method. The simplified time-frequency domain equalization sacrifices some performances due to the simplification of the algorithm, and the equalization of the delay-Doppler domain well eliminates the equalization error caused by the former stage of equalization.

As can be seen from the above description, the embodiments of the present invention use the two-stage equalization technique to equalize the received signals in the time-frequency domain and the delay-doppler domain, respectively, so as to obtain additional diversity gain. Furthermore, in the time-frequency domain equalization, a simplified low-complexity linear equalization algorithm is adopted, and the complexity required by equalization inversion is reduced. And the equalization technology adopted in the time delay-Doppler domain can eliminate ISI caused by the equalization error of the previous stage, and simultaneously effectively utilizes diversity gain by relying on the structure similar to a DFE equalizer, thereby greatly improving the performance of the receiver of the OTFS system. Therefore, the two-stage equalization technology provided by the invention has wide application prospect in practical engineering.

An embodiment of the present invention further provides a computing device, an internal structure of which mainly includes, as shown in fig. 7: at least one processor 702, a memory 704, and a bus 706 connecting the above devices. The at least one processor 702 is configured to execute a memory-stored module of machine-readable instructions 708. In an embodiment of the present invention, the one or more processors execute a machine readable instruction module 708 to implement the signal receiving method.

An embodiment of the present invention also provides a computer-readable medium on which a computer program is stored, which, when executed by a processor, implements the above-described signal receiving method.

Those of ordinary skill in the art will understand that: the discussion of any embodiment above is meant to be exemplary only, and is not intended to intimate that the scope of the disclosure, including the claims, is limited to these examples; within the idea of the invention, also features in the above embodiments or in different embodiments may be combined, steps may be implemented in any order, and there are many other variations of the different aspects of the invention as described above, which are not provided in detail for the sake of brevity.

In addition, well known power/ground connections to Integrated Circuit (IC) chips and other components may or may not be shown within the provided figures for simplicity of illustration and discussion, and so as not to obscure the invention. Furthermore, devices may be shown in block diagram form in order to avoid obscuring the invention, and also in view of the fact that specifics with respect to implementation of such block diagram devices are highly dependent upon the platform within which the present invention is to be implemented (i.e., specifics should be well within purview of one skilled in the art). Where specific details (e.g., circuits) are set forth in order to describe example embodiments of the invention, it should be apparent to one skilled in the art that the invention can be practiced without, or with variation of, these specific details. Accordingly, the description is to be regarded as illustrative instead of restrictive.

While the present invention has been described in conjunction with specific embodiments thereof, many alternatives, modifications, and variations of these embodiments will be apparent to those of ordinary skill in the art in light of the foregoing description. For example, other memory architectures (e.g., dynamic ram (dram)) may use the discussed embodiments.

The embodiments of the invention are intended to embrace all such alternatives, modifications and variances that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, substitutions, improvements and the like that may be made without departing from the spirit and principles of the invention are intended to be included within the scope of the invention.