CN120321069B - A channel estimation device based on power delay spectrum estimation - Google Patents
A channel estimation device based on power delay spectrum estimationInfo
- Publication number
- CN120321069B CN120321069B CN202510241720.5A CN202510241720A CN120321069B CN 120321069 B CN120321069 B CN 120321069B CN 202510241720 A CN202510241720 A CN 202510241720A CN 120321069 B CN120321069 B CN 120321069B
- Authority
- CN
- China
- Prior art keywords
- power
- data
- channel estimation
- frequency domain
- channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W72/00—Local resource management
- H04W72/04—Wireless resource allocation
- H04W72/044—Wireless resource allocation based on the type of the allocated resource
- H04W72/0446—Resources in time domain, e.g. slots or frames
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W72/00—Local resource management
- H04W72/04—Wireless resource allocation
- H04W72/044—Wireless resource allocation based on the type of the allocated resource
- H04W72/0473—Wireless resource allocation based on the type of the allocated resource the resource being transmission power
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W72/00—Local resource management
- H04W72/12—Wireless traffic scheduling
- H04W72/1263—Mapping of traffic onto schedule, e.g. scheduled allocation or multiplexing of flows
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
The invention discloses a channel estimation device based on power delay spectrum estimation, which is characterized in that firstly, resource particles of pilot frequency are extracted from frequency domain resource grids aiming at different physical channels, LS channel estimation is carried out with generated local pilot frequency, then, according to different physical channels, two strategies of segmentation calculation and inversion edge data are adopted, data are changed into power point numbers of 2, channel estimation values at the pilot frequency are converted into a delay domain through a base 2IFFT fast algorithm, a power delay spectrum window is calculated in the delay domain, noise reduction treatment is carried out on the channel estimation values at the pilot frequency of the delay domain, finally, interpolation is carried out through a base 2FFT fast algorithm, and channel estimation values of the whole frequency domain are obtained, and data are integrated and output according to different physical channels. The invention can give consideration to performance and implementation complexity, process time delay and hardware resources, improve the universality, compatibility, portability and expandability of the system and effectively solve the problem of channel estimation realized based on FPGA or ASIC.
Description
Technical Field
The invention belongs to the wireless field, and particularly relates to a channel estimation device based on power delay spectrum estimation.
Background
Wireless communication environments are quite complex and diverse and are often affected by terrain, buildings, weather, and the like. These factors present challenges such as multipath effects, channel fading, time-variability, frequency selectivity, etc., making the signal susceptible to interference and attenuation in transmission. In such an environment, channel estimation plays a vital role. The method is characterized in that the influence of the signal in the transmission process is accurately known by estimating the state of a transmission channel. The channel estimation helps the system effectively compensate signal distortion, track dynamic changes of the channel in real time, and adjust communication parameters according to current channel conditions. The capacity can ensure that the system can still keep lower error rate and strong anti-interference capacity under a complex environment, and ensures stable communication quality.
Currently, there are relatively few implementations of minimum mean square error (Minimum Mean Square Error, MMSE) channel estimation algorithms. The conventional MMSE channel estimation algorithm needs to be implemented based on the result of Least Square (LS) channel estimation, and needs prior information such as a cross-correlation matrix of a known frequency domain pilot frequency band and a frequency domain full-band channel, and an autocorrelation matrix of the frequency domain pilot frequency band.
Because the MMSE channel estimation algorithm involves matrix inversion and matrix multiplication, which is not beneficial to the hardware implementation of a field programmable gate array (Field Programmable GATE ARRAY, FPGA) or an Application SPECIFIC INTEGRATED Circuit (ASIC), the MMSE channel estimation needs to be approximately simplified and deformed to meet the requirements of performance and implementation complexity, computation delay, FPGA or ASIC hardware resource consumption, and the like. Furthermore, in order to use the radix-2 fast fourier transform (Fast Fourier Transform, FFT), the number of data to be processed needs to be a power of 2. At present, the mainstream algorithm is to zero-fill data of non-2 powers to two powers, however, spectrum leakage caused by the strategy is more serious, and the performance of the MMSE algorithm is obviously reduced.
Disclosure of Invention
The invention aims to solve the problems that the calculation of the current MMSE channel estimation algorithm is complex and is not beneficial to the hardware realization, after zero padding is carried out on non-2 power data, the frequency spectrum leakage is serious when a base 2FFT fast algorithm is used, the algorithm performance is seriously reduced, and the like.
The invention discloses a channel estimation device based on power delay spectrum estimation, which is used for processing a signal of a receiving end of a wireless channel and specifically comprises the following steps:
the resource element extraction module is used for extracting resource elements for placing pilot frequency in a frequency domain resource grid according to a 3GPP protocol aiming at different physical channels, wherein the different physical channels comprise a physical shared channel (Physical Downlink SHARED CHANNEL, PDSCH), a physical broadcast channel (Physical Broadcast Channel, PBCH) and a physical control channel (Physical Control Channel, PDCCH);
the LS channel estimation module is used for generating a local pilot frequency sequence and carrying out LS channel estimation with the received pilot frequency;
the IFFT operation module is used for transforming the data into power point numbers of 2 by adopting two strategies of segmentation calculation and inversion edge data according to different physical channels, and converting an LS channel estimated value at a pilot frequency into a time delay domain through base 2 inverse fast Fourier transform (INVERSE FAST Fourier Transform, IFFT);
the power delay spectrum estimation and noise reduction module is used for calculating a power delay spectrum window in a delay domain and carrying out noise reduction treatment on a channel estimation value at a pilot frequency of the delay domain;
the FFT operation module is used for interpolating through a base 2FFT fast algorithm and obtaining a channel estimation value of a full frequency domain;
and the integration output module is used for integrating and outputting the data according to different physical channels.
Furthermore, the device is deployed and realized on an FPGA or ASIC hardware platform, the unified and fixed-point data bit width is realized, the top layer module, the internal sub-module and the calling IP all use the interface protocol of the unified AXI4-Stream, the internal IP of the module adopts a pipeline structure, and the data processing adopts a parallel processing mode.
Further, the local pilot sequence is generated according to a protocol, specifically, a gold pseudo-random sequence is modulated by Quadrature phase shift coding (Quadrature PHASE SHIFT KEYING, QPSK). In the physical shared channel, the generation of the gold sequence is related to factors such as the sequence number of the time slot in the radio frame, the scrambling ID, the number of orthogonal frequency division multiplexing (Orthogonal Frequency Division Multiplexing, OFDM) symbols contained in one time slot, and the OFDM symbol index where the DMRS is placed in one time slot. In the physical broadcast channel, the gold sequence is generated in relation to factors such as SS/PBCH block index, cell ID, field indication, etc. In the physical control channel, the generation of the gold sequence is related to factors such as the sequence number of the time slot in the radio frame, the scrambling ID, the number of OFDM symbols contained in one time slot, the OFDM symbol index where the DMRS is placed in one time slot, and the like.
Furthermore, in the physical shared channel, configuration of different partial Bandwidths (BWP) is supported, so that the number of input frequency domain data is calculated according to configuration information, such as the size of the input RB, and then the maximum number of power points smaller than 2 of the number of input frequency domain data is obtained by table lookup, and the input data is segmented according to the table lookup result. The method is characterized in that data are divided into two sections, wherein the data of the power point number of the first 2 are the first half section, the data of the power point number of the second 2 are the second half section, the first half section is fetched from the top of the complete frequency domain data, the second half section is fetched from the tail of the complete frequency domain data, and the data at the two ends are partially overlapped in the middle. If the size of the input frequency domain data is just the power of 2, the two pieces of data are the size of the input frequency domain data, namely the original data are not segmented, so that the universality of the whole module can be improved. And then performing IFFT fast operation on the two sections of data respectively.
Further, in the physical broadcast channel, the pilot signals on the complete OFDM symbols of columns 2 and 4 in the physical broadcast block are used for estimation. Since the OFDM symbols of the 2 nd and 4 th columns of the synchronous broadcast block fixedly occupy 240 subcarriers in the frequency domain, and the pilots are four times spread, 60 pilot data are fixed on one column. And reversing the last 4 data, supplementing the total data to 64 points, and performing 64-point IFFT fast operation.
Furthermore, in the physical control channel, due to the resource mapping characteristic, pilot symbols are distributed on 1-3 OFDM symbols, and a four-time sparse distribution placing mode is presented on a frequency domain. Therefore, it is necessary to calculate the number of input frequency domain data according to configuration information, such as the size of the input RB, then look up a table to obtain the maximum number of power points smaller than 2 of the input data, and segment the input data according to the result of the table look up. The method is characterized in that data are divided into two sections, wherein the data of the power point number of the first 2 are the first half section, the data of the power point number of the second 2 are the second half section, the first half section is fetched from the top of the complete frequency domain data, the second half section is fetched from the tail of the complete frequency domain data, and the data at the two ends are partially overlapped in the middle. If the size of the input frequency domain data is just the power of 2, the two pieces of data are the size of the input frequency domain data, namely the original data are not segmented, so that the universality of the whole module can be improved. And then performing IFFT fast operation on the two sections of data respectively.
Further, the power delay profile estimation includes the following steps:
(a) Intercepting the LS channel estimation value of the time delay domain, carrying out first noise reduction, and configuring the noise reduction window length from the outside of the module;
(b) Calculating the power of an LS channel estimated value at a pilot frequency position on a time delay domain;
(c) Calculating the average value of the LS channel estimation value power value on the same subcarrier in the same time slot or in a plurality of time slots;
(d) Calculating a time delay domain power spectrum soft window by combining noise power, wherein the noise power is input from the outside of a module and is compared with a time delay domain power average value obtained on the same subcarrier, so that a numerator and a denominator of a power time delay spectrum window are obtained, if the time delay domain power is smaller than the noise power and the noise power is 0, the denominator of the power time delay spectrum window is 1, if the time delay domain power is smaller than the noise power and the noise power is not 0, the denominator of the power time delay spectrum window is the noise power, otherwise, the denominator of the power time delay spectrum window is the time delay domain power, and if the time delay domain power is not smaller than the noise power, the numerator of the power time delay spectrum window is the time delay domain power minus the noise power, otherwise, the denominator is 0;
(e) And adding a soft window to the LS channel estimation value at the pilot frequency in the time delay domain to perform secondary noise reduction.
Furthermore, in the physical shared channel, because the pilot frequency is placed in a double sparse way, the LS channel estimation value after the noise reduction of the time delay domain is subjected to zero padding firstly, so that the data volume after the zero padding reaches twice the pilot frequency data volume after segmentation, then the FFT fast operation of the corresponding point number is carried out, the frequency domain interpolation is completed, and finally the overlapping parts of the frequency domain channel estimation values obtained by the two segments of interpolation are combined, and the complete frequency domain full channel estimation value is output.
Furthermore, in the physical broadcast channel, because the pilot frequency is placed in a four-time sparse way, the LS channel estimation value after the noise reduction of the time delay domain is subjected to zero filling firstly, so that the data volume after the zero filling reaches 4 times of the pilot frequency data volume after the segmentation, then the FFT fast operation of the corresponding point number is carried out, the frequency domain interpolation is completed, the data is finally truncated to 240 points, and the channel estimation value on the 3 rd OFDM symbol of the SSB block can be obtained by averaging or interpolating the channel estimation value on the 2 nd OFDM symbol and the 4 th OFDM symbol.
Furthermore, in the physical control channel, because the pilot frequency is placed in a four-fold sparse way, the LS channel estimation value after the time delay domain noise reduction is subjected to zero padding firstly, so that the data volume after the zero padding reaches 4 times of the pilot frequency data volume after segmentation, then the FFT fast operation of the corresponding point number is carried out, the frequency domain interpolation is completed, and finally the overlapped parts of the frequency domain channel estimation values obtained by the two segments of interpolation are combined, and the complete frequency domain full channel estimation value is output.
Compared with the prior art, the invention has the following advantages:
1. The architecture can support the channel estimation processing of a receiving end of a plurality of physical channels including a physical shared channel, a physical broadcast channel, a physical control channel and the like, has good channel estimation performance under the channel conditions of a Gaussian Noise channel (ADDITIVE WHITE AWGN), a multipath channel and the like, and improves the universality of the module;
2. The method has the advantages that the data quantity is controlled to be the power point number of 2 by using the segmentation and inversion edge data, the hardware implementation is more convenient by using the FFT fast algorithm, the frequency spectrum leakage caused by zero padding is optimized, and the problem of serious performance degradation caused by the frequency spectrum leakage is solved;
3. The data processing adopts a pipelining parallel processing mode, so that the data processing time delay is effectively reduced;
4. For the physical control channel, channel estimation caused by interleaving when the resource grid is discontinuous in the frequency domain can be supported, and the compatibility of the system is improved;
5. the reasonable design of the calculation structure unifies the bit width of the fixed-point data, thereby saving the hardware realization resources;
6. The module portability is improved by using the interface protocol of unified AXI 4-Stream;
7. the hardware implementation logic of the data processing module is optimized, and the expandability of the module is improved, namely, when the power value of the delay domain signal is averaged, multiple groups of additional pilots in one time slot or multiple groups of pilots and additional pilots in multiple time slots can be averaged, so that the overall performance is improved.
Drawings
Fig. 1 is a schematic diagram of a channel estimation device based on power delay spectrum estimation according to an embodiment of the present invention.
Fig. 2 is a schematic diagram of a PDSCH channel frequency domain resource grid structure.
Fig. 3 is a schematic diagram of an implementation structure of an LS channel estimation module according to an embodiment of the present invention.
Fig. 4 is a schematic diagram of a data segment implementation structure according to an embodiment of the present invention.
Fig. 5 is a schematic diagram of a power delay spectrum window implementation structure according to an embodiment of the present invention.
Fig. 6 is a schematic diagram of an implementation structure of a power accumulation average module according to an embodiment of the invention.
Fig. 7 is a flowchart of power delay spectrum window calculation according to an embodiment of the present invention.
Fig. 8 is a schematic diagram of a zero padding module implementation structure according to an embodiment of the present invention.
Detailed Description
The technical scheme of the invention is further described below with reference to the attached drawings and specific embodiments.
The embodiment of the invention discloses a channel estimation device structure based on power delay spectrum estimation, which can be deployed and realized on an FPGA or ASIC hardware platform.
As shown in fig. 1, the channel estimation device based on power delay spectrum estimation includes a resource element extraction module, an LS channel estimation module, an IFFT operation module, a power delay spectrum estimation and noise reduction module, an FFT operation module and an integration output module. The system comprises a base station side or a user terminal, a resource element extraction module, an IFFT operation module, a power delay spectrum estimation and noise reduction module, an FFT operation module and an integration output module, wherein the base station side or the user terminal is used for extracting resource elements for placing pilots in a frequency domain resource grid according to different physical channels, the different physical channels comprise a PDSCH channel, a PBCH channel and a PDCCH channel, the LS channel estimation module is used for generating local pilot sequences through the protocols and carrying out LS channel estimation with received pilots, the IFFT operation module is used for changing data into power point numbers of 2 according to the different physical channels, the LS channel estimation value at the pilots is converted into a time delay domain through a base 2IFFT fast algorithm, the power delay spectrum estimation and noise reduction module is used for calculating a power delay spectrum window in the time delay domain and carrying out noise reduction processing on the channel estimation value at the pilots of the time delay domain, the FFT operation module is used for carrying out interpolation through the base 2FFT fast algorithm and obtaining the channel estimation value of the whole frequency domain, and the integration output module is used for integrating and outputting the data according to the different physical channels.
The working flow of the channel estimation device based on the power time delay spectrum estimation is that at a receiving end, frequency domain data, noise variance and system configuration information, such as RB size to be processed, frame timing information including subcarrier index, OFDM symbol index, time slot number index and frame number index, are input by an upper module, and after the frequency domain data is processed by the flow, frequency domain channel coefficient estimation values are output by taking time slots as units.
In this embodiment, in the PDSCH, the time-frequency resource position of the Demodulation reference signal (DMRS) is related to factors such as the duration of the physical channel in one slot, the DMRS mapping method, and the additional pilot configuration. In the PBCH, the time-frequency resource position of the DMRS is placed on the 2,3,4 th symbol of the synchronous broadcast block (Synchronization Signal/PBCH, SSB). In each Resource Block (RB), a DMRS specific mapping location is related to a cell ID configuration. In the PDCCH, the position of the DMRS on the resource grid is the 1 st, 5 th and 9 th resource elements in each resource block.
Fig. 2 is a schematic diagram of a PDSCH channel frequency domain resource grid structure. Taking PDSCH channel as an example, in the resource element extraction module, pilot symbols are extracted according to the frequency domain position of the physical shared channel indicated in the physical control channel and the pilot mapping mode, and then subsequent channel estimation is performed. In the physical shared channel, there are various pilot mapping modes, as shown in fig. 2, the pilot is mapped on the 3 rd and 10 th OFDM symbols in the time domain, and the frequency domain adopts a mapping mode of 2 times sparse division.
For the PBCH channel, it is necessary to first find the resource grid where the SSB is located, and then extract the complete pilot symbols in columns 2 and 4 in the PBCH to perform subsequent channel estimation.
For the PDCCH channel, blind detection is required, and then pilot symbols in a continuous frequency domain on a control resource set (Control Resource Set, CORESET) are extracted to perform subsequent channel estimation.
Fig. 3 is a schematic diagram of an implementation structure of the LS channel estimation module. The LS channel estimation module mainly realizes the function of carrying out LS channel estimation on the extracted received pilot frequency and a locally generated pilot frequency sequence. Because of the specificity of the local pilot sequence, its modulus is 1, and there is only a phase difference between different pilot symbols. The local pilot sequence is generated according to a protocol, specifically, the gold pseudo-random sequence is subjected to quadrature phase shift coded QPSK modulation. In the PDSCH channel, the generation of the gold sequence is related to factors such as the sequence number of the time slot in the radio frame, the scrambling ID, the number of the OFDM symbols contained in one time slot, and the OFDM symbol index where the DMRS is placed in one time slot. In the PBCH channel, the gold sequence is generated in relation to factors such as SS/PBCH block index, cell ID, field indication, etc. In the PDCCH, the gold sequence is generated according to the sequence number of the time slot in the radio frame, the scrambling ID, the number of OFDM symbols contained in one time slot, the OFDM symbol index placed in one time slot by the DMRS, and other factors.
Specifically, the specific generation mode of the local pilot frequency in the PDSCH channel is that gold pseudo-random sequences are QPSK modulated, namely:
Wherein r represents a generated demodulation reference signal, m represents a symbol sequence number after modulation, c represents a pseudo random sequence, and the specific generation mode is as follows:
c(n)=[x1(n+Nc)+x2(n+Nc)]mod2
Where N c = 1600, the sequences x 1 (N) and x 2 (N) are defined as:
x1(n+31)=[x1(n+3)+x1(n)]mod2
x2(n+31)=[x2(n+3)+x2(n+2)+x2(n+1)+x2(n)]mod2
the initial values of x 1 (n) and x 2 (n) are respectively:
Wherein, the Wherein, the For the number of OFDM symbols in one slot,For the slot number in a frame, l is the OFDM symbol index in which the DMRS is placed in a slot, n SCID is the scrambling ID, n SCID ε {0,1} is given by the DMRS initialization field,Representing an identifier associated with the reference signal,Given by the relevant parameters in the DMRS-DownlinkConfig IE, if no DMRS-DownlinkConfig IE is providedWherein the method comprises the steps ofIs the physical cell ID of the serving cell.
In the PBCH channel, the reference sequence and the pseudo random sequence are generated in the same manner as the PDSCH, and the initialization value is generated in the following manner:
where n hf is the field indication of the PBCH, i SSB is the index of the SS/PBCH block, Is the physical cell ID of the serving cell.
In the PDCCH, the reference sequence and the pseudo random sequence are generated in the same manner as the PDSCH, and the initialization value is generated in the following manner:
Wherein, the For the number of OFDM symbols in one slot,For a slot number within a frame, i is an OFDM symbol index within a slot, N ID is a scrambling ID, N ID e {0,1} is given by the DMRS initialization field, N ID e {0, 1..65535 } is given by a related parameter in pdcch-DMRS-ScramblingID, if pdcch-DMRS-ScramblingID is not providedWherein the method comprises the steps ofIs the physical cell ID of the serving cell.
Because of the specificity of the local pilot sequence, the hardware implementation of the LS channel estimation algorithm can convert division into phase rotation of the local pilot by means of the CORDIC algorithm, so that the hardware resource consumption of the FPGA or the ASIC can be reduced, and the processing time delay is reduced. The module outputs LS channel estimates at the pilot symbols.
The IFFT operation module is used for carrying out fast IFFT operation after mainly segmenting the LS channel estimation value of the frequency domain into powers of 2 and converting the LS channel estimation value into the LS channel estimation value of the time delay domain. In the physical shared channel, different partial bandwidth configurations are supported, so that the number of input frequency domain data is calculated according to configuration information, such as the size of an input RB, then the maximum number of power points which is less than 2 of the number of input frequency domain data is obtained through table lookup, and the input data is segmented according to the table lookup result. The method is characterized in that data are divided into two sections, wherein the data of the power point number of the first 2 are the first half section, the data of the power point number of the second 2 are the second half section, the first half section is fetched from the top of the complete frequency domain data, the second half section is fetched from the tail of the complete frequency domain data, and the data at the two ends are partially overlapped in the middle. If the size of the input frequency domain data is just the power of 2, the two pieces of data are the size of the input frequency domain data, namely the original data are not segmented, so that the universality of the whole module can be improved. And then performing IFFT fast operation on the two sections of data respectively.
In the physical broadcast channel, the pilot signals on the complete OFDM symbols of columns 2 and 4 in the physical broadcast block are used for estimation. Since the OFDM symbols of the 2 nd and 4 th columns of the synchronous broadcast block fixedly occupy 240 subcarriers in the frequency domain, and pilots are four times spread, 60 pilot data are fixed on one column. And reversing the last 4 data, supplementing the total data to 64 points, and performing 64-point IFFT fast operation.
In the physical control channel, due to the resource mapping characteristic, pilot symbols are distributed on 1-3 OFDM symbols, and a four-time sparse placement mode is presented on a frequency domain. Therefore, it is necessary to calculate the number of input frequency domain data according to configuration information, such as the size of the input RB, then look up a table to obtain the maximum number of power points smaller than 2 of the input data, and segment the input data according to the result of the table look up. The method is characterized in that data are divided into two sections, wherein the data of the power point number of the first 2 are the first half section, the data of the power point number of the second 2 are the second half section, the first half section is fetched from the top of the complete frequency domain data, the second half section is fetched from the tail of the complete frequency domain data, and the data at the two ends are partially overlapped in the middle. If the size of the input frequency domain data is just the power of 2, the two pieces of data are the size of the input frequency domain data, namely the original data are not segmented, so that the universality of the whole module can be improved. And then performing IFFT fast operation on the two sections of data respectively.
Illustratively, FIG. 4 is a schematic diagram of a data segment implementation. As shown in fig. 4 (a), for LS channel estimation values of non-2 power points in PDCCH and PDSCH, data is divided into two pieces of power data of 2 before and after, and the two pieces of data overlap partially. If the data quantity of the LS channel estimation value is just the power of 2, the two segments of data are both original data, so that the universality of the module is improved. The length of the intercepted power point number of the front section 2 and the rear section 2 can be obtained through table lookup according to the input system configuration information, such as configuration information of RB size and the like.
As shown in fig. 4 (b), for PBCH, the LS channel estimation value is fixed to 60 data, so that the end data is inverted to 4 bits, i.e., the 56 th to 59 th data are directly inverted and spliced to the 60 th data, and then 64 bit data are directly formed, which is the smallest two power data larger than 60, and thus the algorithm performance is better than the segmentation into 232 point data.
As shown in fig. 4 (c), the hardware implementation of the IFFT operation is to use the FFT IP core in the FPGA or ASIC hardware platform to complete the power point serial IFFT operation corresponding to 2. The input time delay domain channel estimation value is divided into the data of the power point number of the front section and the rear section 2 through the counting of the frame timing information, and then, for the PDSCH and the PDCCH channels, a FIFO IP is used for buffering, for the PBCH channels, a stack function module is used for reversing the data, and finally, the two sections of data are serially input into an FFT IP core. And configuring the information such as the operation point number of the FFT IP, the flag bit of the last data, the valid flag bit of the data and the like according to externally input system configuration information such as the size of a resource block and the like so as to enable the FFT IP to perform dynamic FFT operation point number configuration.
Since the receiving end data may be affected by noise, the first noise reduction operation is needed to be performed on the delay domain data, namely, the noise reduction window length is configured outside the module, and a window can be configured before and after the whole data so as to cope with frequency domain leakage. And cutting off output data of the FFT IP according to the window length, reserving data in front and rear windows, and carrying out zero setting processing on data outside the windows. The window length can be configured from outside, so that the universality and the flexibility of the system are improved.
The power delay spectrum estimation and noise reduction module is mainly used for completing the calculation of the power delay spectrum window in the delay domain and reducing the noise of the channel estimation value at the pilot frequency on the delay domain. The power delay spectrum estimation comprises the following steps of (a) cutting off a delay domain LS channel estimation value, carrying out first noise reduction, (b) calculating the power of the LS channel estimation value at a pilot frequency position on the delay domain, (c) calculating the average value of LS channel estimation value power values on the same subcarrier in the same time slot or multiple time slots, (d) calculating a delay domain power spectrum soft window by combining noise power, and (e) adding the soft window to the LS channel estimation value at the pilot frequency position on the delay domain, and carrying out second noise reduction.
Fig. 5 is a schematic diagram of a power delay spectrum window implementation structure, and the module mainly comprises a power accumulation average module, a power delay spectrum window module, a secondary windowing module and a zero padding module. The module obtains the power of the input complex-valued signal data and obtains the average power on the same subcarrier of a plurality of OFDM symbols in the same time slot or on the same subcarrier of a plurality of OFDM symbols in a plurality of time slots. The module supports window length calculation of single time slot and multiple time slots, and expandability and universality of the module are improved.
The implementation structure of the power accumulation average module is shown in fig. 6, and mainly consists of two DSPIP. The time delay domain data after the first noise reduction is in a complex form, the real part and the imaginary part of the data are separated, and the real part square is calculated by using a function of DSPIP A, B and C. Wherein, the A and B ports are both input with real part of data, and the C port is selected to input with output data of the subsequent stage DSPIP or 0. The module simultaneously realizes the accumulation function, so that when the input data is the data on the first column of the accumulated OFDM symbols, the selector selects 0 for input, and if the input data is the data on other columns of OFDM symbols in the same time slot or the data on OFDM symbols in the subsequent time slot, the second DSPIP is selected for output. The second DSPIP selects the function of a x B + C for the calculation of the square of the imaginary part and sums the square of the real part. Wherein, the A and B ports are both input with the imaginary part of the data, and the C port is selected to input the output result of the first DSPIP. The output result of the second stage DSPIP needs to be selectively processed, that is, if the accumulated values on the same subcarrier on a plurality of OFDM symbols in the same time slot are calculated, or the accumulated values on the same subcarrier on a plurality of OFDM symbols in a plurality of time slots are calculated, the power values are averaged and output, and if the accumulated calculation is not completed, the calculation result is stored in the FIFO, and the data in the next group of OFDM symbols is waited for input.
The power time delay spectrum window generating module is mainly composed of DIV IP, and the generating mode of the power time delay spectrum window numerator and denominator is shown as a power time delay spectrum window calculation flow chart in figure 7. The noise power is input from the outside and then compared with the average value of the delay domain power obtained on the same subcarrier, so as to obtain the numerator and denominator of the power delay spectrum window. The method comprises the following steps that if the delay domain power is smaller than the noise power and the noise power is 0, the denominator of the power delay spectrum window is 1, if the delay domain power is smaller than the noise power and the noise power is not 0, the denominator of the power delay spectrum window is the noise power, otherwise, the denominator of the power delay spectrum window is the delay domain power, if the delay domain power is not smaller than the noise power, the numerator of the power delay spectrum window is the delay domain power minus the noise power, otherwise, the numerator is 0. And finally, division operation is carried out by using a DIV IP core in the FPGA or ASIC hardware platform, and a power delay spectrum window is output. The reasonable design of the calculation structure is realized, and the magnitude of the numerator and the denominator are close in the designed soft window, so that the data can be scaled according to a unified format, the problem that the bit width of the fixed-point data is required to be increased due to overlarge or undersize divider results is avoided, and the hardware realization resources are saved.
In the secondary windowing module, the time delay domain LS channel estimation data after the first noise reduction is multiplied with the generated soft window result, the multiplication can be realized through MULT IP in FPGA or ASIC, and the module completes the secondary noise reduction in the time delay domain.
As shown in the implementation structure diagram of the zero padding module in fig. 8, in the zero padding module, zero padding output is required to be performed on data after the secondary windowing. The specific method is that for PDSCH channel, the length of power data of each segment 2 is assumed to be N, because the pilot frequency is double sparse in frequency domain, N0 are needed to be added in the middle of effective data of front and back windows to make the total length of data be 2N, for PBCH and PDCCH channel, the length of power data of each segment 2 is assumed to be N, because the pilot frequency is four times sparse in frequency domain, 3N 0 are needed to be added in the middle of effective data of front and back windows to make the total length of data be 4N. The module outputs LS channel estimation data after zero padding of the time delay domain after the secondary noise reduction.
And inputting the time delay domain data after the secondary noise reduction into an FFT operation module, and correspondingly configuring the information such as the operation point number of the FFT IP, the flag bit of the last data, the effective flag bit of the data and the like. The FFT operation point number can be obtained by looking up a table according to the input system configuration information, such as the configuration information of the size of the resource block.
The main function of the integration output module is to integrate and output the segmented frequency domain data finally. And for PDSCH and PDCCH channels, averaging overlapped parts of the segmented full-frequency-domain channel estimation and outputting. For the PBCH channel, the first 240 points of the frequency domain channel estimation value are intercepted and output as the channel estimation value of the PBCH frequency domain channel, and for the channel estimation value of the PBCH on the 3 rd incomplete OFDM symbol in the SSB, the channel estimation value of the PBCH on the 2 rd and 4 th OFDM symbols can be obtained through average or interpolation.
In the channel estimation device based on power delay spectrum estimation, the top layer module, the inner sub-module and the calling IP all use the interface protocol of the unified AXI4-Stream, so that the portability of the module is improved. And the IP inside the module adopts a pipeline structure, and the data processing adopts a parallel processing mode, so that the data processing time delay is effectively reduced.
Claims (8)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202510241720.5A CN120321069B (en) | 2025-03-03 | 2025-03-03 | A channel estimation device based on power delay spectrum estimation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202510241720.5A CN120321069B (en) | 2025-03-03 | 2025-03-03 | A channel estimation device based on power delay spectrum estimation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN120321069A CN120321069A (en) | 2025-07-15 |
| CN120321069B true CN120321069B (en) | 2026-04-21 |
Family
ID=96326730
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202510241720.5A Active CN120321069B (en) | 2025-03-03 | 2025-03-03 | A channel estimation device based on power delay spectrum estimation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN120321069B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117579437A (en) * | 2024-01-15 | 2024-02-20 | 北京智联安科技有限公司 | Method and device for determining frequency domain correlation coefficient and computer equipment |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008042575A (en) * | 2006-08-07 | 2008-02-21 | Matsushita Electric Ind Co Ltd | Receiver |
| KR20080020934A (en) * | 2006-09-01 | 2008-03-06 | 한국전자통신연구원 | Uplink Signal Transmission Method, Transmission Device, Generation Method, and Generation Device in Communication System |
| US8243836B2 (en) * | 2008-12-17 | 2012-08-14 | Electronics And Telecommunications Research Institute | Channel estimation apparatus and method in mobile communication system having dispersed pilot |
| CN101808053A (en) * | 2010-02-26 | 2010-08-18 | 大唐联诚信息系统技术有限公司 | Channel estimation method and device based on OFDM (Orthogonal Frequency Division Multiplexing) |
| CN107483373B (en) * | 2017-08-03 | 2020-12-15 | 北京睿信丰科技有限公司 | An anti-multipath iterative weighted LMMSE channel estimation method and device |
| CN109818885B (en) * | 2017-11-22 | 2020-07-31 | 电信科学技术研究院 | Channel estimation method, device, equipment and computer readable storage medium |
| CN111342919B (en) * | 2018-12-18 | 2022-04-05 | 大唐移动通信设备有限公司 | Method and equipment for estimating frequency domain channel correlation value of channel |
| CN111010359B (en) * | 2019-12-19 | 2022-06-03 | 紫光展锐(重庆)科技有限公司 | Timing and frequency offset estimation method and device based on PBCH (physical broadcast channel), and terminal equipment |
| CN115733591A (en) * | 2021-08-30 | 2023-03-03 | 华为技术有限公司 | Downlink transmission method and related device based on multi-station cooperation |
| CN114978822B (en) * | 2022-05-20 | 2024-05-10 | Oppo广东移动通信有限公司 | Signal processing method, device, chip and storage medium |
| CN116232812B (en) * | 2023-02-27 | 2025-06-03 | 山东大学 | Least Squares Channel Estimation Method for OTFS System Based on Deep Neural Network |
-
2025
- 2025-03-03 CN CN202510241720.5A patent/CN120321069B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117579437A (en) * | 2024-01-15 | 2024-02-20 | 北京智联安科技有限公司 | Method and device for determining frequency domain correlation coefficient and computer equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN120321069A (en) | 2025-07-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101242383B (en) | A Channel Estimation Method | |
| KR100923892B1 (en) | Fast fourier transform twiddle multiplication | |
| JP2010016913A (en) | Transmitter and transmitting method, and device and method for calculating delay times | |
| CN101669312B (en) | OFDM reception device, OFDM reception method, OFDM reception circuit, integrated circuit, and program | |
| KR100958231B1 (en) | Fast fourier transform processing in an ofdm system | |
| KR20040032683A (en) | Fast fourier transform apparatus for wireless lan system | |
| CN100479445C (en) | Orthogonal frequency division multiplexing communication method using prefix and suffix signals | |
| CN101729479A (en) | Blind channel estimation method based on OFDM signal cyclostationarity | |
| JP3691357B2 (en) | Carrier arrangement method, transmission apparatus, and reception apparatus in orthogonal frequency division multiplexing transmission system | |
| CN120321069B (en) | A channel estimation device based on power delay spectrum estimation | |
| CN1909528B (en) | Channel estimation method and device for crossing frequency division multiplexing system | |
| CN101258488A (en) | Fast Fourier Transform Processing in OFDM System | |
| CN100574305C (en) | The signal-noise ratio estimation method and the device that shine upon based on planisphere in the multicarrier system | |
| CN115913420A (en) | SRS-based SNR estimation method and device in 5G small base station system | |
| CN107566311B (en) | Transmission method of RB F-OFDM system based on resource block filtering | |
| CN106953822B (en) | novel generalized multi-carrier communication method suitable for time-frequency double-selective fading channel | |
| US20080219360A1 (en) | System and method for performing channel estimation using interpolation | |
| CN103188192B (en) | It is applied to the base band processing device of video sensor | |
| Tsai et al. | A low-power multicarrier-CDMA downlink baseband receiver for future cellular communication systems | |
| Chang et al. | FPGA-based Design of a Pulsed-OFDM System | |
| CN101557366B (en) | Channel estimation method based on orthogonal frequency division multiplexing system | |
| KR101275693B1 (en) | Method and apparatus of estimating frequency offset based on blind approach in wireless communication system | |
| CN201107851Y (en) | Signal-noise ratio estimating device based on planisphere mapping in multi-carrier system | |
| CN101039288B (en) | Channel Estimation Method and Device Against Timing Deviation in OFDM System | |
| Nirmala et al. | An area‐efficient FFT processor for the OFDMA transceiver communication system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |