Channel estimation method, apparatus and channel sounding reference signal receiver

Abstract

This disclosure provides a channel estimation method applicable to the field of wireless communication technology. The method includes: estimating and compensating for the group time offset TO0; and using a cyclic shift matrix W... cs After decoupling the group time offset compensated channel, the channel pre-estimate for each port is obtained. Based on the channel pre-estimate for each port, the coarse residual time offset for each port is estimated. Based on the channel estimation matrix W, the channel estimate H for each port is estimated. i The channel estimation matrix W is based on the coarse residual time offset of each port and the cyclic shift matrix W0. cs Constructed; based on the channel estimation H of each port i The present disclosure also provides a channel estimation apparatus, device, storage medium, and program product, which calculate the precise remaining time offset of each port and perform time offset compensation on the channel estimation of each port based on the precise remaining time offset of each port to generate a target channel estimation result.

Description

Technical Field

[0001] This disclosure relates to the field of wireless communication technology, and more specifically to a channel estimation method, apparatus, channel sounding reference signal receiver, medium, and program product. Background Technology

[0002] In a 3GPP NR system, multiple uplink UE antenna ports transmit SRS signals in a time-frequency multiplexing manner. The arrival times of the SRS signals received by the base station from each port will deviate from each other. When there is a large time deviation between the received SRS signals from each port, the relevant channel estimation method cannot effectively remove the inter-port interference in channel estimation, resulting in large channel interference in the end.

[0003] It should be noted that the information disclosed in the background section above is only used to enhance the understanding of the background of this disclosure, and therefore may include information that does not constitute prior art known to those skilled in the art. Summary of the Invention

[0004] In view of the above problems, this disclosure provides a channel estimation method, apparatus, channel sounding reference signal receiver, medium, and program products for removing inter-port interference.

[0005] According to a first aspect of this disclosure, a channel estimation method is provided, comprising: estimating and compensating for group time offset TO0;

[0006] Through the cyclic shift matrix W cs Decouple the group time offset compensated channel to obtain the channel pre-estimate for each port. in This represents the maximum number of ports that can be reused in the current OFDM symbol for time and frequency.

[0007] Channel pre-estimation based on each port Estimate the coarse residual time bias of each port

[0008] The channel estimate H for each port is estimated based on the channel estimation matrix W. i The channel estimation matrix W is based on the coarse residual time offset of each port. and the cyclic shift matrix W cs Constructed;

[0009] Based on the channel estimation H of each port i Calculate the exact remaining time offset for each port as well as

[0010] Based on the precise remaining time offset of each port Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results.

[0011] According to embodiments of this disclosure, estimating and compensating for group time bias TO0 includes:

[0012] Extract the first frequency domain data y1 of the channel sounding reference signal at the receiving end;

[0013] The base sequence in the first frequency domain data is removed by the least squares method to generate the second frequency domain data y2;

[0014] Based on the second frequency domain data y2, the group time bias TO0 is estimated;

[0015] The second frequency domain data y2 is time-biased and compensated according to the group time bias TO0 to generate the third frequency domain data y3.

[0016] According to embodiments of this disclosure, estimating the group time bias TO0 based on the second frequency domain data y2 includes:

[0017] Calculate the second frequency domain data y2 per interval The average phase difference e between the data points, in For each interval The channel after extracting y2 from each point L is the length of the ZC sequence; and

[0018] The group time offset TO0 is calculated based on the average phase difference e. Where angle is the arctangent function. N FFT It is the FFT length.

[0019] According to embodiments of this disclosure,

[0020] Third frequency domain data Δ is the starting subcarrier index value of the SRS received signal across the entire frequency band.

[0021] According to embodiments of this disclosure, the cyclic shift matrix W... cs Decouple the group time offset compensated channel to obtain the channel pre-estimate for each port. include:

[0022] Construct the cyclic shift matrix W cs = [d0, ..., d i ,...], Among them W cs yes The matrix, d i It is the cyclic shift phase rotation factor. It is the occupied circular shift, N p This is the current number of multiplexed ports;

[0023] The channel pre-estimation for each port is obtained using the following formula.

[0024]

[0025] in (·) H This represents the conjugate transpose of a matrix.

[0026] According to embodiments of this disclosure, the channel pre-estimation based on each port... Estimate the coarse residual time bias of each port include:

[0027] Calculate the phase difference between adjacent subcarriers as well as

[0028] Estimate the coarse residual time bias for each port

[0029] According to embodiments of this disclosure, based on the coarse residual time offset of each port... and the cyclic shift matrix W cs Construct the channel estimation matrix W, including:

[0030] Construct matrix A, A = W cs .*W To W To = [b0, b1, ... b i ], Where .* represents matrix dot product. as well as

[0031] Construct the channel estimation matrix W, W = (A H *A) -1 *A H W is The matrix.

[0032] According to embodiments of this disclosure, the channel estimate H for each port is estimated based on the channel estimation matrix W. i ,include:

[0033] For each of y3 For each subcarrier, the channel estimate H is obtained using the following formula. i H i (m) = W*y3, where

[0034] According to embodiments of this disclosure, the channel estimation H based on each port... i Calculate the exact remaining time offset for each port include:

[0035] Calculate the phase difference between adjacent subcarriers as well as

[0036] Estimate the remaining capacity of each port

[0037] According to embodiments of this disclosure, based on the precise remaining time offset of each port... Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results. include:

[0038] in,

[0039] A second aspect of this disclosure provides a channel estimation apparatus, comprising:

[0040] The group time bias compensation module is used to estimate and compensate for the group time bias TO0;

[0041] The channel pre-estimation module is used to estimate the channel using the cyclic shift matrix W. cs Decouple the group time offset compensated channel to obtain the channel pre-estimate for each port. in This represents the maximum number of ports that can be reused in the current OFDM symbol for time and frequency.

[0042] The coarse residual time offset estimation module is used for channel pre-estimation based on each port. Estimate the coarse residual time bias of each port

[0043] The channel estimation module is used to estimate the channel estimate H for each port based on the channel estimation matrix W. i The channel estimation matrix W is based on the coarse residual time offset of each port. and the cyclic shift matrix W cs Constructed;

[0044] The precise remaining time offset estimation module is used for channel estimation H based on each port. i Calculate the exact remaining time offset for each port

[0045] The time offset compensation module is used to accurately calculate the remaining time offset of each port. Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results.

[0046] According to embodiments of this disclosure, the group time offset compensation module includes:

[0047] The signal receiving submodule is used to extract the first frequency domain data y1 of the channel sounding reference signal at the receiving end;

[0048] The second frequency domain data generation submodule is used to remove the base sequence in the first frequency domain data by the least squares method to generate the second frequency domain data y2.

[0049] The group time offset estimation submodule is used to estimate the group time offset TO0 based on the second frequency domain data y2;

[0050] The first time offset compensation submodule is used to perform time offset compensation on the second frequency domain data y2 according to the group time offset TO0 to generate the third frequency domain data y3.

[0051] According to embodiments of this disclosure, the group time offset estimation submodule includes: a phase difference calculation unit and a group time offset estimation unit.

[0052] The phase difference calculation unit is used to calculate the second frequency domain data y2 per interval. The average phase difference e between the data points, in L is the length of the ZC sequence; and

[0053] The group time offset estimation unit is used to calculate the group time offset TO0 based on the average phase difference e. Where angle is the arctangent function. N FFT It is the FFT length.

[0054] According to embodiments of this disclosure, third frequency domain data Δ is the starting subcarrier index value of the SRS received signal across the entire frequency band.

[0055] According to embodiments of this disclosure, the channel pre-estimation module includes a cyclic shift matrix construction submodule and a channel pre-estimation submodule.

[0056] The cyclic shift matrix construction submodule is used to construct the cyclic shift matrix W. cs = [d0, ..., d i ,...], Among them W cs yes The matrix, d i It is the cyclic shift phase rotation factor. It is the occupied circular shift, Np This is the current number of multiplexed ports;

[0057] The channel pre-estimation submodule is used to obtain the channel pre-estimation for each port using the following formula.

[0058]

[0059] in (·) H This represents the conjugate transpose of a matrix.

[0060] According to embodiments of this disclosure, the coarse residual time offset estimation module includes: a phase difference calculation submodule and a coarse residual time offset estimation submodule.

[0061] The phase difference calculation submodule is used to calculate the phase difference between adjacent subcarriers. as well as

[0062] The coarse residual time-partial estimation submodule is used to estimate the residual for each port.

[0063] According to embodiments of this disclosure, a channel estimation matrix construction module is also included.

[0064] The channel estimation matrix construction module is specifically used to construct matrix A, A = W. cs .*W To W To = [b0, b1, ... b i ], Where .* represents matrix dot product. And construct the channel estimation matrix W, W = (A H *A) -1 *A H W is The matrix.

[0065] According to embodiments of this disclosure, the channel estimation module is specifically used for each of y3 For each subcarrier, the channel estimate H is obtained using the following formula. i H i (m) = W*y3, where

[0066] According to embodiments of this disclosure, the precise residual time offset estimation module is specifically used to calculate the phase difference between adjacent subcarriers. And estimate the remaining capacity of each port

[0067] A third aspect of this disclosure provides a channel sounding reference signal receiver, comprising: one or more processors; and a memory for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors perform the channel estimation method described above.

[0068] A fourth aspect of this disclosure also provides a computer-readable storage medium having executable instructions stored thereon, which, when executed by a processor, cause the processor to perform the channel estimation method described above.

[0069] The fifth aspect of this disclosure also provides a computer program product, including a computer program that, when executed by a processor, implements the channel estimation method described above.

[0070] The channel estimation method provided in this disclosure obtains a common channel estimate for each port by removing the base sequence at the receiving end using the least squares method, and then utilizes the characteristics of SRS itself to estimate the channel estimate for each port. Point extraction is used to calculate the group time offset estimate without needing to calculate and average the time offsets of all ports. Furthermore, compensation only needs to be performed on the common channel, resulting in low computational complexity. The inverse matrix is ​​calculated based on the cyclic shift of each port and the estimated time offset, and then the inverse matrix is ​​used to calculate the continuous... The subcarriers are smoothed and filtered to obtain a channel estimate without interference from other ports. Attached Figure Description

[0071] The foregoing contents, as well as other objects, features, and advantages of this disclosure, will become clearer from the following description of embodiments with reference to the accompanying drawings, in which:

[0072] Figure 1 The illustration schematically depicts application scenarios of channel estimation methods, apparatus, devices, media, and program products according to embodiments of the present disclosure;

[0073] Figure 2 A flowchart illustrating a channel estimation method according to an embodiment of the present disclosure is shown schematically;

[0074] Figure 3 A flowchart illustrating a group time bias estimation compensation method according to an embodiment of the present disclosure is shown schematically.

[0075] Figure 4 A flowchart illustrating a channel estimation matrix construction method according to an embodiment of the present disclosure is shown schematically;

[0076] Figure 5 A schematic block diagram of a channel estimation apparatus according to an embodiment of the present disclosure is shown; and

[0077] Figure 6A block diagram of an electronic device suitable for implementing a channel estimation method according to an embodiment of the present disclosure is shown schematically. Detailed Implementation

[0078] The embodiments of the present disclosure will now be described with reference to the accompanying drawings. However, it should be understood that these descriptions are exemplary only and are not intended to limit the scope of the disclosure. In the following detailed description, numerous specific details are set forth to provide a thorough understanding of the embodiments of the present disclosure for ease of explanation. However, it will be apparent that one or more embodiments may be practiced without these specific details. Furthermore, descriptions of well-known structures and techniques are omitted in the following description to avoid unnecessarily obscuring the concepts of the present disclosure.

[0079] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit this disclosure. The terms “comprising,” “including,” etc., as used herein indicate the presence of the stated features, steps, operations, and / or components, but do not exclude the presence or addition of one or more other features, steps, operations, or components.

[0080] All terms used herein (including technical and scientific terms) have the meanings commonly understood by those skilled in the art, unless otherwise defined. It should be noted that the terms used herein are to be interpreted in a manner consistent with the context of this specification, and not in an idealized or overly rigid way.

[0081] When using expressions such as "at least one of A, B, and C", they should generally be interpreted in accordance with the meaning that is commonly understood by a person skilled in the art (e.g., "a system having at least one of A, B, and C" should include, but is not limited to, a system having A alone, a system having B alone, a system having C alone, a system having A and B, a system having A and C, a system having B and C, and / or a system having A, B, and C, etc.).

[0082] In a 3GPP NR system employing OFDM technology, the SRS signal transmitted from a set of time-frequency multiplexed antenna ports is generated as follows: using a ZC sequence q of length L. 0 =[q(1) q(2) … q(m) … q(L)] T The base sequence is used; the SRS signal transmitted by each port is the base sequence after passing through... The result of performing a DFT transform after cyclic shifting of points; the sequence obtained by the DFT transform is divided into K-order intervals. TC (Two or four) subcarriers are loaded with a sequence element on the starting subcarrier Δ of the OFDM symbol, where... and All are integers. This indicates the maximum number of ports that can be time-frequency reused in the current OFDM symbol.

[0083] In related technologies, time offset compensation for channel estimation includes the following steps:

[0084] 1) Extract the frequency domain data of each user at the receiving end;

[0085] 2) Generate the transmission sequence for each user according to the 3GPP protocol;

[0086] 3) Using the sequences in 1) and 2), obtain a coarse channel estimate for each user based on the least squares method;

[0087] 4) Coarse channel estimation for each user, every N m Subcarrier averaging yields the intermediate channel estimate for each user;

[0088] 5) Calculate the first time offset TA1 for each user using the intermediate channel estimation for each user;

[0089] 6) Average TA1 for all users to obtain TA. est ;

[0090] 7) Using TA est Time-bias compensation is applied to the coarse channel estimate obtained in step 3) for each user to obtain the compensated channel estimate for each user.

[0091] 8) Perform continuous N-stage compensation channel estimation for each user. m Subcarrier averaging yields a smoothed channel estimate for each user.

[0092] 9) Calculate the second time offset TA2 for each user using the smoothed channel estimation.

[0093] 10) Use TA2 to perform time-bias compensation on the smoothed channel estimate of each user obtained in 8) to obtain the final channel estimate of each user.

[0094] When there is a time deviation between the received SRS signals of each port, the coarse channel estimate of each port obtained by steps 1), 2), and 3) includes the interference components of the signals of other ports, that is, inter-port interference. All subsequent operations cannot effectively remove the inter-port interference of the channel estimate, especially when the time deviation is large, the performance deteriorates sharply.

[0095] Based on the above-mentioned technical problems, embodiments of this disclosure provide a channel estimation method, including: estimating and compensating for group time offset TO0; and using a cyclic shift matrix W cs Decouple the group time offset compensated channel to obtain the channel pre-estimate for each port. in The maximum number of ports that can be reused in the current OFDM symbol; channel pre-estimation based on each port. Estimate the coarse residual time bias of each port The channel estimate H for each port is estimated based on the channel estimation matrix W. i The channel estimation matrix W is based on the coarse residual time offset of each port. and the cyclic shift matrix W cs Constructed; based on the channel estimation H of each port i Calculate the exact remaining time offset for each port And based on the precise remaining time offset of each port Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results.

[0096] Figure 1 The illustration schematically depicts application scenarios of channel estimation methods, apparatus, devices, media, and program products according to embodiments of the present disclosure.

[0097] like Figure 1 As shown, application scenario 100 according to this embodiment may include a scenario of transmitting and receiving signals from a multi-antenna system. Network 103 is used as a medium to provide a communication link between terminal devices 101, 102, base station 104, and server 105. Network 103 may be a wireless communication link, etc.

[0098] Users can use terminal devices 101 and 102 to receive or send messages via network 103. Various communication client applications can be installed on terminal devices 101 and 102, such as shopping applications, web browser applications, search applications, instant messaging tools, email clients, social media platform software, etc. (for example only).

[0099] Terminal devices 101 and 102 can be various electronic devices with displays and web browsing capabilities, including but not limited to smartphones, tablets, laptops, and desktop computers.

[0100] Base station 104 can be a base station providing signal transmission services, and server 105 can be a server providing signal encoding services or a server providing signal processing services. It should be noted that the channel estimation method provided in this embodiment can generally be executed by server 105. Correspondingly, the channel estimation device provided in this embodiment can generally be located in server 105. The channel estimation method provided in this embodiment can also be executed by a server or server cluster that is different from server 105 and can communicate with terminal devices 101, 102 and / or server 105. Correspondingly, the channel estimation device provided in this embodiment can also be located in a server or server cluster that is different from server 105 and can communicate with terminal devices 101, 102 and / or server 105.

[0101] It should be understood that Figure 1 The number of terminal devices, networks, and servers shown is merely illustrative. Depending on implementation needs, any number of terminal devices, networks, and servers can be included.

[0102] Figure 2 A flowchart illustrating a channel estimation method according to an embodiment of the present disclosure is shown schematically.

[0103] like Figure 2 As shown, the channel estimation method of this embodiment includes operations S210 to S260, which can be executed by a server or other computing device.

[0104] When operating S210, the bias TO0 is estimated and compensated for in the group.

[0105] For methods of estimating group time bias, please refer to Figure 3 Operations 211 to 214 are shown. In this embodiment of the present disclosure, after removing the base sequence at the receiving end using the least squares method, a common channel estimate for each port is obtained. Then, taking advantage of the characteristics of SRS itself, a value is extracted for each N points to calculate the group time offset estimate. It is not necessary to calculate the time offset values ​​of all ports and then average them. Furthermore, only the common channel needs to be compensated during compensation, resulting in low computational complexity. Subsequent calculations only require storing one common channel estimate, thus occupying little memory.

[0106] In operation S220, the cyclic shift matrix W is used. cs Decouple the group time offset compensated channel to obtain the channel pre-estimate for each port. in This represents the maximum number of ports that can be reused in the current OFDM symbol for time and frequency.

[0107] According to embodiments of this disclosure, a cyclic shift matrix W is constructed. cs = [d0, ..., d i ,...], Among them W cs yes The matrix, d i It is the cyclic shift phase rotation factor. It is the occupied circular shift, N p This is the number of multiplexed ports; the channel pre-estimate for each port is obtained using the following formula.

[0108]

[0109] in (·) H This represents the conjugate transpose of a matrix.

[0110] During operation S230, channel pre-estimation is performed based on each port. Estimate the coarse residual time bias of each port

[0111] According to embodiments of this disclosure, the phase difference between adjacent subcarriers is calculated. Estimate the coarse residual time bias for each port

[0112] In operation S240, the channel estimate H for each port is estimated based on the channel estimation matrix W. i The channel estimation matrix W is based on the coarse residual time offset of each port. and the cyclic shift matrix W cs Constructed.

[0113] According to embodiments of this disclosure, for each of y3 For each subcarrier, the channel estimate H is obtained using the following formula. i H i (m) = W*y3, where m = 0, 1, ..., M-1,

[0114] In operation S250, based on the channel estimation H of each port... i Calculate the exact remaining time offset for each port

[0115] According to embodiments of this disclosure, the phase difference between adjacent subcarriers is calculated. conj(H i (m-1)); Estimate the exact remaining time offset for each port.

[0116] During operation S260, based on the precise remaining time offset of each port... Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results.

[0117] According to embodiments of this disclosure, in,

[0118] The above steps are the same for multiple SRS symbols and multiple receiving antennas. Using the signal estimation method provided in this disclosure, the group time offset is estimated and compensated, the coarse residual time offset of each port is estimated after decoupling, a channel estimation matrix W is constructed using the coarse residual time offset, W is used for decoupling, the precise residual time offset of each port is calculated using the decoupled channel, and the precise residual time offset is compensated to obtain a channel estimate without inter-port interference.

[0119] Figure 3 A flowchart illustrating a group time bias estimation and compensation method according to an embodiment of this disclosure is shown schematically. Figure 3 As shown, operation S210 includes operations S211 to S214.

[0120] In operation S211, the first frequency domain data y1 of the channel sounding reference signal is extracted at the receiving end.

[0121] In operation S212, the base sequence in the first frequency domain data is removed by the least squares method to generate the second frequency domain data y2.

[0122] In operation S213, the group time bias TO0 is estimated based on the second frequency domain data y2.

[0123] According to embodiments of this disclosure, the second frequency domain data y2 is calculated per interval. The average phase difference e between the data points, in For each interval The channel after y2 extraction at each point, where L is the length of the ZC sequence.

[0124] The group time offset TO0 is calculated based on the average phase difference e. Where angle is the arctangent function. N FFT It is the FFT length.

[0125] In operation S214, time offset compensation is performed on the second frequency domain data y2 according to the group time offset TO0 to generate the third frequency domain data y3.

[0126] Third frequency domain data Δ is the starting subcarrier index value of the SRS received signal across the entire frequency band.

[0127] After removing the base sequence at the receiving end using the least squares method, a common channel estimate for each port is obtained. Then, taking advantage of the characteristics of SRS itself, N points are sampled to calculate the group time offset estimate. It is not necessary to calculate the time offset values ​​of all ports and then average them. Furthermore, only the common channel needs to be compensated during compensation, resulting in low computational complexity. Subsequent calculations only require storing one common channel estimate, thus occupying little memory.

[0128] Figure 4 A flowchart illustrating a channel estimation matrix construction method according to an embodiment of this disclosure is shown schematically. Figure 4 As shown, it includes operation S310 and operation S320.

[0129] In operation S310, construct matrix A.

[0130] Construct matrix A, A = W cs .*W To W To = [b0, b1, ... b i ], Where .* represents matrix dot product.

[0131] In operation S320, the channel estimation matrix W is constructed.

[0132] Construct the channel estimation matrix W, W = (A H *A) -1 *A H W is The matrix.

[0133] When the time offset values ​​of all ports are the same, the time offset effect can be eliminated through common time offset estimation and compensation. In this way, the channel coefficients of different ports are reduced to only the cyclic shift coefficient. The cyclic shift coefficient is continuously... The subcarriers are orthogonal and can be connected in successive waves. The subcarriers are added together to cancel out the channel at other ports. However, when the time offsets of each port are different, the channel after common time offset estimation and compensation still has the influence of relative time offsets of different ports, and the channel coefficients of different ports are no longer just cyclic shift coefficients, but continuous... The channel coefficients on each subcarrier are no longer orthogonal, therefore it is impossible to pass through each The method of averaging subcarriers to eliminate the channels of other ports results in a large amount of final channel interference. This embodiment of the present disclosure calculates the inverse matrix based on the cyclic shift of each port and the estimated time offset, and then uses the inverse matrix to perform continuous... The subcarriers are smoothed and filtered to obtain a channel estimate without interference from other ports.

[0134] Based on the above-described channel estimation method, this disclosure also provides a channel estimation apparatus. The following will be combined with... Figure 5 The device is described in detail.

[0135] Figure 5 A schematic block diagram of a channel estimation apparatus according to an embodiment of the present disclosure is shown.

[0136] like Figure 5 As shown, the channel estimation device 500 in this embodiment includes a group time offset compensation module 410, a channel pre-estimation module 420, a coarse residual time offset estimation module 430, a channel estimation module 440, a precise residual time offset estimation module 450, and a time offset compensation module 460.

[0137] The group time offset compensation module 410 is used to estimate and compensate for the group time offset TO0. In one embodiment, the group time offset compensation module 410 can be used to perform the operation S210 described above, which will not be repeated here.

[0138] The channel pre-estimation module 420 is used to estimate the channel using the cyclic shift matrix W. cs Decouple the group time offset compensated channel to obtain the channel pre-estimate for each port. in This represents the maximum number of ports for time-frequency multiplexing in the current OFDM symbol. In one embodiment, the channel pre-estimation module 420 can be used to perform the operation S220 described above, which will not be repeated here.

[0139] The coarse residual time offset estimation module 430 is used for port-based channel pre-estimation. Estimate the coarse residual time bias of each port In one embodiment, the coarse residual time bias estimation module 430 can be used to perform the operation S230 described above, which will not be repeated here.

[0140] The channel estimation module 440 is used to estimate the channel estimate H for each port based on the channel estimation matrix W. i The channel estimation matrix W is based on the coarse residual time offset of each port. and the cyclic shift matrix W cs Constructed. In one embodiment, the channel estimation module 440 can be used to perform the operation S240 described above, which will not be repeated here.

[0141] The precise remaining time offset estimation module 450 is used for channel estimation H based on each port. i Calculate the exact remaining time offset for each port In one embodiment, the precise remaining time offset estimation module 450 can be used to perform the operation S230 described above, which will not be repeated here.

[0142] The time offset compensation module 460 is used to calculate the accurate remaining time offset of each port. Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results. In one embodiment, the time offset compensation module 460 can be used to perform the operation S230 described above, which will not be repeated here.

[0143] According to embodiments of this disclosure, the group time offset compensation module includes:

[0144] The signal receiving submodule is used to extract the first frequency domain data y1 of the channel sounding reference signal at the receiving end;

[0145] The second frequency domain data generation submodule is used to remove the base sequence in the first frequency domain data by the least squares method to generate the second frequency domain data y2.

[0146] The group time offset estimation submodule is used to estimate the group time offset TO0 based on the second frequency domain data y2;

[0147] The first time offset compensation submodule is used to perform time offset compensation on the second frequency domain data y2 according to the group time offset TO0 to generate the third frequency domain data y3.

[0148] According to embodiments of this disclosure, the group time offset estimation submodule includes: a phase difference calculation unit and a group time offset estimation unit.

[0149] The phase difference calculation unit is used to calculate the second frequency domain data y2 per interval. The average phase difference e between the data points, in L is the length of the ZC sequence; and

[0150] The group time offset estimation unit is used to calculate the group time offset TO0 based on the average phase difference e. Where angle is the arctangent function. N FFT It is the FFT length.

[0151] According to embodiments of this disclosure, third frequency domain data Δ is the starting subcarrier index value of the SRS received signal across the entire frequency band.

[0152] According to embodiments of this disclosure, the channel pre-estimation module includes a cyclic shift matrix construction submodule and a channel pre-estimation submodule.

[0153] The cyclic shift matrix construction submodule is used to construct the cyclic shift matrix W. cs = [d0, ..., di ,...], Among them W cs yes The matrix, d i It is the cyclic shift phase rotation factor. It is the occupied circular shift, N p It is the number of multiplexed ports;

[0154] The channel pre-estimation submodule is used to obtain the channel pre-estimation for each port using the following formula.

[0155]

[0156] in (·) H This represents the conjugate transpose of a matrix.

[0157] According to embodiments of this disclosure, the coarse residual time offset estimation module includes: a phase difference calculation submodule and a coarse residual time offset estimation submodule.

[0158] The phase difference calculation submodule is used to calculate the phase difference between adjacent subcarriers. as well as

[0159] The coarse residual time-partial estimation submodule is used to estimate the residual for each port.

[0160] According to embodiments of this disclosure, a channel estimation matrix construction module is also included.

[0161] The channel estimation matrix construction module is specifically used to construct matrix A, A = W. cs .*W To W To = [b0, b1, ... b i ], Where .* represents matrix dot product. And construct the channel estimation matrix W, W = (A H *A) -1 *A H W is The matrix.

[0162] According to embodiments of this disclosure, the channel estimation module is specifically used for each of y3 For each subcarrier, the channel estimate H is obtained using the following formula. i H i (m) = W*y3, where

[0163] According to embodiments of this disclosure, the precise residual time offset estimation module is specifically used to calculate the phase difference between adjacent subcarriers. And estimate the remaining capacity of each port

[0164] According to embodiments of this disclosure, any multiple modules among the group time offset compensation module 410, channel pre-estimation module 420, coarse residual time offset estimation module 430, channel estimation module 440, precise residual time offset estimation module 450, and time offset compensation module 460 can be combined into one module, or any one of these modules can be split into multiple modules. Alternatively, at least some of the functions of one or more of these modules can be combined with at least some of the functions of other modules and implemented in one module. According to embodiments of this disclosure, at least one of the group time offset compensation module 410, channel pre-estimation module 420, coarse residual time offset estimation module 430, channel estimation module 440, precise residual time offset estimation module 450, and time offset compensation module 460 can be at least partially implemented as hardware circuits, such as field-programmable gate arrays (FPGAs), programmable logic arrays (PLAs), systems-on-a-chip, systems-on-a-substrate, systems-on-package, application-specific integrated circuits (ASICs), or any other reasonable means of integrating or packaging circuits, or implemented in hardware or firmware, or in any one of software, hardware, and firmware implementations, or in a suitable combination of any of these. Alternatively, at least one of the group time offset compensation module 410, channel pre-estimation module 420, coarse residual time offset estimation module 430, channel estimation module 440, precise residual time offset estimation module 450, and time offset compensation module 460 can be at least partially implemented as computer program modules, which can perform corresponding functions when the computer program module is run.

[0165] Figure 6 A block diagram of an electronic device suitable for implementing a channel estimation method according to an embodiment of the present disclosure is shown schematically.

[0166] like Figure 6 As shown, an electronic device 500 according to an embodiment of the present disclosure includes a processor 501, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 502 or a program loaded from a storage portion 508 into a random access memory (RAM) 503. The processor 501 may include, for example, a general-purpose microprocessor (e.g., a CPU), an instruction set processor and / or an associated chipset and / or a special-purpose microprocessor (e.g., an application-specific integrated circuit (ASIC)), etc. The processor 501 may also include onboard memory for caching purposes. The processor 501 may include a single processing unit or multiple processing units for performing different actions of the method flow according to an embodiment of the present disclosure.

[0167] RAM 503 stores various programs and data required for the operation of electronic device 500. Processor 501, ROM 502, and RAM 503 are interconnected via bus 504. Processor 501 performs various operations of the method flow according to embodiments of the present disclosure by executing programs in ROM 502 and / or RAM 503. It should be noted that the programs may also be stored in one or more memories other than ROM 502 and RAM 503. Processor 501 may also perform various operations of the method flow according to embodiments of the present disclosure by executing programs stored in said one or more memories.

[0168] According to embodiments of this disclosure, the electronic device 500 may further include an input / output (I / O) interface 505, which is also connected to a bus 504. The electronic device 500 may also include one or more of the following components connected to the I / O interface 505: an input section 506 including a keyboard, mouse, etc.; an output section 507 including a cathode ray tube (CRT), liquid crystal display (LCD), etc., and a speaker, etc.; a storage section 508 including a hard disk, etc.; and a communication section 509 including a network interface card such as a LAI card, modem, etc. The communication section 509 performs communication processing via a network such as the Internet. A drive 510 is also connected to the I / O interface 505 as needed. A removable medium 511, such as a disk, optical disk, magneto-optical disk, semiconductor memory, etc., is installed on the drive 510 as needed so that computer programs read from it can be installed into the storage section 508 as needed.

[0169] This disclosure also provides a computer-readable storage medium, which may be included in the device / apparatus / system described in the above embodiments; or it may exist independently and not assembled into the device / apparatus / system. The computer-readable storage medium carries one or more programs, which, when executed, implement the channel estimation method according to the embodiments of this disclosure.

[0170] According to embodiments of this disclosure, the computer-readable storage medium may be a non-volatile computer-readable storage medium, such as including, but not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In this disclosure, the computer-readable storage medium may be any tangible medium that contains or stores a program that can be used by or in conjunction with an instruction execution system, apparatus, or device. For example, according to embodiments of this disclosure, the computer-readable storage medium may include ROM 502 and / or RAM 503 and / or one or more memories other than ROM 502 and RAM 503 described above.

[0171] Embodiments of this disclosure also include a computer program product comprising a computer program containing program code for performing the methods shown in the flowchart. When the computer program product is run on a computer system, the program code is used to cause the computer system to implement the channel estimation method provided in the embodiments of this disclosure.

[0172] When the computer program is executed by the processor 501, it performs the functions defined in the system / apparatus of this disclosure embodiments. According to embodiments of this disclosure, the systems, apparatuses, modules, units, etc., described above can be implemented by computer program modules.

[0173] In one embodiment, the computer program may rely on a tangible storage medium such as an optical storage device or a magnetic storage device. In another embodiment, the computer program may also be transmitted and distributed in the form of signals over a network medium, and may be downloaded and installed via the communication section 509, and / or installed from a removable medium 511. The program code contained in the computer program can be transmitted using any suitable network medium, including but not limited to: wireless, wired, etc., or any suitable combination thereof.

[0174] In such an embodiment, the computer program can be downloaded and installed from a network via communication section 509, and / or installed from removable medium 511. When the computer program is executed by processor 501, it performs the functions defined in the system of this disclosure embodiment. According to embodiments of this disclosure, the systems, devices, apparatuses, modules, units, etc., described above can be implemented by computer program modules.

[0175] According to embodiments of this disclosure, program code for executing the computer programs provided in embodiments of this disclosure can be written in any combination of one or more programming languages. Specifically, these computational programs can be implemented using high-level procedural and / or object-oriented programming languages, and / or assembly / machine languages. Programming languages ​​include, but are not limited to, languages ​​such as Java, C++, Python, "C", or similar programming languages. The program code can execute entirely on the user's computing device, partially on the user's device, partially on a remote computing device, or entirely on a remote computing device or server. In cases involving remote computing devices, the remote computing device can be connected to the user's computing device via any type of network, including a local area network (LAI) or a wide area network (WAI), or it can be connected to an external computing device (e.g., via the Internet using an Internet service provider).

[0176] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this disclosure. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram or flowchart, and combinations of blocks in a block diagram or flowchart, may be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0177] Those skilled in the art will understand that the features described in the various embodiments and / or claims of this disclosure can be combined or combined in various ways, even if such combinations or combinations are not explicitly described in this disclosure. In particular, the features described in the various embodiments and / or claims of this disclosure can be combined or combined in various ways without departing from the spirit and teachings of this disclosure. All such combinations and / or combinations fall within the scope of this disclosure.

[0178] The embodiments of this disclosure have been described above. However, these embodiments are for illustrative purposes only and are not intended to limit the scope of this disclosure. Although various embodiments have been described above, this does not mean that the measures in the various embodiments cannot be used advantageously in combination. The scope of this disclosure is defined by the appended claims and their equivalents. Various substitutions and modifications can be made by those skilled in the art without departing from the scope of this disclosure, and all such substitutions and modifications should fall within the scope of this disclosure.

Claims

1. A method of channel estimation, characterized by, Includes the following steps: Estimate and compensate for group time bias TO0; By cyclically shifting the matrix W cs Decoupling the channel after group delay compensation to obtain channel estimates of each port wherein is the maximum number of ports that can be multiplexed in time and frequency in the current OFDM symbol; Channel pre-estimation based on each port Estimate the coarse residual time offset for each port The channel estimate H for each port is estimated based on the channel estimation matrix W. i The channel estimation matrix W is based on the coarse residual time offset of each port. and the cyclic shift matrix W cs Constructed; Based on the channel estimation H of each port i Calculate the exact remaining time offset for each port as well as Based on the precise remaining time offset of each port Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results.

2. The method according to claim 1, characterized in that, The estimation and compensation of the group time bias TO0 includes: Extract the first frequency domain data y1 of the channel sounding reference signal at the receiving end; The base sequence in the first frequency domain data is removed by the least squares method to generate the second frequency domain data y2; Based on the second frequency domain data y2, the group time bias TO0 is estimated; The second frequency domain data y2 is time-biased and compensated according to the group time bias TO0 to generate the third frequency domain data y3.

3. The method according to claim 2, characterized in that, Estimating the group time bias TO0 based on the second frequency domain data y2 includes: Calculate the second frequency domain data y2 per interval The average phase difference e between the data points, in For each interval The channel after extracting y2 from each point L is the length of the ZC sequence; and The group time offset TO0 is calculated based on the average phase difference e. Where angle is the arctangent function. N FFT It is the FFT length.

4. The method according to claim 3, characterized in that, Third frequency domain data Δ is the starting subcarrier index value of the SRS received signal across the entire frequency band.

5. The method according to claim 3, characterized in that, Through the cyclic shift matrix W cs Decouple the third frequency domain data y3 to obtain the channel pre-estimate for each port. include: Construct the cyclic shift matrix W cs = [d0, ..., d i ,...], Among them W cs yes The matrix, d i It is the cyclic shift phase rotation factor. It is the occupied circular shift, N p This is the current number of multiplexed ports; The channel pre-estimation for each port is obtained using the following formula. Where m = 0, 1, ..., M-1, (·) H This represents the conjugate transpose of a matrix.

6. The method according to claim 5, characterized in that, The channel pre-estimation based on each port Estimate the coarse residual time offset for each port include: Calculate the phase difference between adjacent subcarriers as well as Estimate the coarse residual time bias for each port 7. The method according to claim 6, characterized in that, Based on the coarse residual time offset of each port and the cyclic shift matrix W cs Construct the channel estimation matrix W, including: Construct matrix A, A = W cs .*W To W To = [b0, b1, ... b i ], Where .* represents matrix dot product. as well as Construct the channel estimation matrix W, W = (A H *A) -1 *A H W is The matrix.

8. The method according to claim 7, characterized in that, said estimating the channel estimate H for each port from the channel estimate matrix W i comprising: For each of y3 For each subcarrier, the channel estimate H is obtained using the following formula. i H i (m) = W*y3, where m = 0, 1, ..., M-1, 9. The method according to claim 8, characterized in that, Based on the channel estimation H of each port i Calculate the exact remaining time offset for each port include: Calculate the phase difference between adjacent subcarriers as well as Estimate the exact remaining time offset for each port 10. The method according to claim 9, characterized in that, Based on the precise remaining time offset of each port Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results. include: Where m = 0, 1, ..., M-1, 11. A channel estimation device, characterized in that, include: The group time bias compensation module is used to estimate and compensate for the group time bias TO0; The channel pre-estimation module is used to estimate the channel using the cyclic shift matrix W. cs Decouple the group time offset compensated channel to obtain the channel pre-estimate for each port. in This represents the maximum number of ports that can be reused in the current OFDM symbol for time and frequency. The coarse residual time offset estimation module is used for channel pre-estimation based on each port. Estimate the coarse residual time offset for each port The channel estimation module is used to estimate the channel estimate H for each port based on the channel estimation matrix W. i The channel estimation matrix W is based on the coarse residual time offset of each port. and the cyclic shift matrix W cs Constructed; The precise remaining time offset estimation module is used for channel estimation H based on each port. i Calculate the exact remaining time offset for each port The time offset compensation module is used to accurately calculate the remaining time offset of each port. Time offset compensation is performed on the channel estimates for each port to generate the target channel estimation results.

12. A channel sounding reference signal receiver, comprising: One or more processors; Storage device for storing one or more programs. When the one or more programs are executed by the one or more processors, the one or more processors perform the method according to any one of claims 1 to 10.

13. A computer-readable storage medium having executable instructions stored thereon, which, when executed by a processor, cause the processor to perform the method according to any one of claims 1 to 10.

14. A computer program product comprising a computer program that, when executed by a processor, implements the method according to any one of claims 1 to 10.

Images