CN109818885B - 一种信道估计方法、装置、设备及计算机可读存储介质 - Google Patents
一种信道估计方法、装置、设备及计算机可读存储介质 Download PDFInfo
- Publication number
- CN109818885B CN109818885B CN201711171479.5A CN201711171479A CN109818885B CN 109818885 B CN109818885 B CN 109818885B CN 201711171479 A CN201711171479 A CN 201711171479A CN 109818885 B CN109818885 B CN 109818885B
- Authority
- CN
- China
- Prior art keywords
- pdp
- ifft
- fft
- current system
- points
- 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
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000001228 spectrum Methods 0.000 claims abstract description 18
- 238000004891 communication Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims description 130
- 238000009499 grossing Methods 0.000 claims description 72
- 238000001914 filtration Methods 0.000 claims description 58
- 230000009467 reduction Effects 0.000 claims description 50
- 238000012545 processing Methods 0.000 claims description 47
- 238000012937 correction Methods 0.000 claims description 27
- 238000005259 measurement Methods 0.000 claims description 25
- 230000008569 process Effects 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 16
- 230000007704 transition Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000006866 deterioration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Landscapes
- Complex Calculations (AREA)
Abstract
本发明提供一种信道估计方法、装置、设备及计算机可读存储介质，涉及通信技术领域，用以提高系统性能。本发明的信道估计方法，包括：获取当前系统带宽；根据所述当前系统带宽确定IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；根据所述IFFT点数和所述当前系统带宽确定FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。本发明可提高系统性能。
Description
技术领域
本发明涉及通信技术领域，尤其涉及一种信道估计方法、装置、设备及计算机可读存储介质。
背景技术
在LTE(Long Term Evolution，长期演进)系统中，信道估计的一般步骤为：首先基于导频进行LS(Least Square，最小二乘)初估计，然后采用二维MMSE(minimum meansquare error，最小均方差)滤波得到所有RE(Resource Element，资源元素)对应位置的信道估计值。
由于二维MMSE滤波的滤波系数复杂度较高，一般分解成频域和时域两个一维的MMSE滤波。由无线通信理论可知，信道的时延功率谱和频域相关性是一对傅里叶变换对。频域的MMSE滤波系数需要通过频域相关性来获得。一般地，频域相关性通过估计信道的PDP(Power Delay Profile，功率时延谱)做FFT(Fast Fourier Transformation，快速傅里叶变换)得到。PDP的估计一般通过将LS估计的结果变化到时域获得。但是PDP的采样率决定了频域相关性的间隔，从而决定了信道估计的性能。就工程常用系统而言，频域相关性的频率间隔和子载波间隔并没有完全匹配，从而导致频域MMSE滤波的系数不够准确，造成性能下降。
发明内容
有鉴于此，本发明提供一种信道估计方法、装置、设备及计算机可读存储介质，用以提高系统性能。
为解决上述技术问题，第一方面，本发明实施例提供一种信道估计方法，包括：
获取当前系统带宽；
根据所述当前系统带宽确定IFFT(Inverse Fast Fourier Transform，快速傅里叶逆变换)点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；
根据所述IFFT点数和所述当前系统带宽确定FFT(Fast FourierTransformation，快速傅里叶变换)点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；
根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，在所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH之后，所述方法还包括：
计算所述PDP的采样率和系统采样率；
判断所述PDP的采样率是否和所述系统采样率相等；
若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正。
其中，所述根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计，包括：
根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，在所述根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP之后，所述方法还包括：
对所述PDP进行平滑滤波和减噪声处理；
所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP，包括：
判断所述当前系统带宽是否大于第一带宽门限值；
若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；
利用IFFT点数进行IFFT运算，估计PDP。
其中，所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
判断所述当前系统带宽是否小于第二带宽门限值；
若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；
若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；
若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；
利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
将所述PDP进行补零操作；
利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述计算所述PDP的采样率和系统采样率，包括：
确定所述PDP的采样率为IFFT的长度和第四系数的乘积；
确定所述系统采样率为OFDM(Orthogonal Frequency Division Multiplexing，正交频分复用)调制时的FFT的长度和第五系数的乘积。
其中，所述若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正，包括：
计算所述PDP的采样率和所述系统采样率的比值；
根据所述比值对所述RHH进行线性插值。
其中，所述对所述PDP进行平滑滤波和减噪声处理，包括：
进行噪声测量，获得噪声测量结果；
将所述PDP减去所述噪声测量结果，获得减噪声后的PDP；
利用减噪声后的PDP和前一平滑滤波处理后的PDP平滑值，获得当前PDP平滑值。
其中，根据以下公式，利用减噪声后的PDP和前一平滑滤波和减噪声处理后的PDP平滑值，获得当前PDP平滑值，包括：
PDP_Smooth(n)＝(1-α)PDP_Smooth(n-1)+αPDP_Current(n)
其中，PDP_Smooth(n)表示当前PDP平滑值，PDP_Smooth(n-1)表示前一平滑滤波处理后的PDP平滑值，PDP_Current(n)表示当前子帧减噪声后的PDP，α表示平滑因子。
第二方面，本发明实施例提供一种信道估计装置，包括：
第一获取模块，用于获取当前系统带宽；
第一处理模块，用于根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；
第二处理模块，用于根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；
第二获取模块，用于根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，所述装置还包括：
计算模块，用于计算所述PDP的采样率和系统采样率；
判断模块，用于判断所述PDP的采样率是否和所述系统采样率相等；
插值修正模块，用于若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正；
所述第二获取模块具体用于，根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，所述装置还包括：
滤波和减噪声模块，用于对所述PDP进行平滑滤波和减噪声处理；
所述第二处理模块具体用于，根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述第一处理模块包括：
判断子模块，用于判断所述当前系统带宽是否大于第一带宽门限值；
确定子模块，用于若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；
处理子模块，用于利用IFFT点数进行IFFT运算，估计PDP。
其中，所述第二处理模块包括：
第一判断子模块，用于判断所述当前系统带宽是否小于第二带宽门限值；
第一确定子模块，用于若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；
第二判断子模块，用于若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；
第二确定子模块，用于若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；
处理子模块，用于利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
第三方面，本发明实施例提供一种通信设备，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序；所述处理器执行所述程序时实现以下步骤：
获取当前系统带宽；
根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；
根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
计算所述PDP的采样率和系统采样率；
判断所述PDP的采样率是否和所述系统采样率相等；
若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
对所述PDP进行平滑滤波和减噪声处理；
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
判断所述当前系统带宽是否大于第一带宽门限值；
若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；
利用IFFT点数进行IFFT运算，估计PDP。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
判断所述当前系统带宽是否小于第二带宽门限值；
若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；
若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；
若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；
利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
将所述PDP进行补零操作；
利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
确定所述PDP的采样率为IFFT的长度和第四系数的乘积；
确定所述系统采样率为正交频分复用OFDM调制时的FFT的长度和第五系数的乘积。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
计算所述PDP的采样率和所述系统采样率的比值；
根据所述比值对所述RHH进行线性插值。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
进行噪声测量，获得噪声测量结果；
将所述PDP减去所述噪声测量结果，获得减噪声后的PDP；
利用减噪声后的PDP和前一平滑滤波处理后的PDP平滑值，获得当前PDP平滑值。
其中，所述处理器还用于读取存储器中的程序，执行下列过程：
根据以下公式，利用减噪声后的PDP和前一平滑滤波和减噪声处理后的PDP平滑值，获得当前PDP平滑值，包括：
PDP_Smooth(n)＝(1-α)PDP_Smooth(n-1)+αPDP_Current(n)
其中，PDP_Smooth(n)表示当前PDP平滑值，PDP_Smooth(n-1)表示前一平滑滤波处理后的PDP平滑值，PDP_Current(n)表示当前子帧减噪声后的PDP，α表示平滑因子。
第四方面，本发明实施例提供一种计算机可读存储介质，用于存储计算机程序，所述计算机程序被处理器执行时实现如第一方面所述的方法中的步骤。
本发明的上述技术方案的有益效果如下：
在本发明实施例中，根据当前系统带宽确定的IFFT点数、根据IFFT点数和当前系统带宽确定的FFT点数，可使得在获得频域相关性项链RHH时的间隔更为准确，使得频域MMSE滤波的系数更为准确，提高了系统性能。
附图说明
图1为本发明实施例的信道估计方法的流程图；
图2为本发明实施例的信道估计方法的流程图；
图3为步骤202的流程图；
图4为步骤203的流程图；
图5为步骤204的流程图；
图6为线性插值示意图；
图7为本发明实施例的信道估计装置的示意图；
图8为本发明实施例的信道估计装置的第一结构图；
图9为本发明实施例的信道估计装置的第二结构图；
图10为本发明实施例的通信设备的示意图。
具体实施方式
下面将结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
如图1所示，本发明实施例的信道估计方法，包括：
步骤101、获取当前系统带宽。
步骤102、根据所述当前系统带宽确定IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP。
具体的，在此步骤中，判断所述当前系统带宽是否大于第一带宽门限值。若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积，利用IFFT点数进行IFFT运算，估计PDP。
其中，所述第一带宽门限值、第一系数可根据经验设置。
步骤103、根据所述IFFT点数和所述当前系统带宽确定FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
具体的，判断所述当前系统带宽是否小于第二带宽门限值。若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积。若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值。若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积。最后，利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
在利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH时，首先将所述PDP进行补零操作，利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述第二、第三带宽门限值、第一系数和第二系数可根据经验设置。
步骤104、根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
在本发明实施例中，根据当前系统带宽确定的IFFT点数、根据IFFT点数和当前系统带宽确定的FFT点数，可使得在获得频域相关性项链RHH时的间隔更为准确，使得频域MMSE滤波的系数更为准确，提高了系统性能。
如图2所示，本发明实施例的信道估计方法，包括：
步骤201、获取当前系统带宽。
步骤202、根据所述当前系统带宽确定IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP。
在本发明实施例中，使用IFFT来估计PDP，避免IDFT(Inverse Discrete FourierTransform，离散傅里叶逆变换)带来的复杂度。同时，结合频域相关性计算时的FFT以及后续的插值校正，使得在快速简单实现的同时，保证了频域相关性间隔的精确性。
具体的，在此步骤中，如图3所示，此步骤包括：
步骤301、判断所述当前系统带宽是否大于第一带宽门限值ThB。
步骤302、若所述当前系统带宽大于所述第一带宽门限值ThB，则IFFT点数为实际信道估计向量长度的最小2的整数幂次。
步骤303、若所述当前系统带宽小于或等于所述第一带宽门限值ThB，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数(αIFFT)的乘积。
步骤304、利用IFFT点数进行IFFT运算，估计PDP。
其中，所述第一带宽门限值、第一系数可根据经验设置。
通过以上描述可以看出，在本发明实施例中，对小于ThB系统带宽时的IFFT点数进行了特殊处理，将IFFT点数乘以一可配系数αIFFT进行放大，从而提高PDP估计的准确性。
步骤203、对所述PDP进行平滑滤波和减噪声处理。
结合图4所示，本步骤包括：
步骤401、进行噪声测量，获得噪声测量结果。
其中，噪声的测量方法如下：为了避免能量泄露导致的影响，对每根接收天线只通过中间两段数据来算噪声，并且取min值作为最后结果。公式如下：
其中，
和
分别表示第r根接收天线第一段和第二段的噪声功率。lIFFT表示IFFT点数，
表示第r根接收天线估计PDP的第k′根径的功率。
为最终得到的第r根接收天线的PDP的平均噪声功率。
步骤402、将所述PDP减去所述噪声测量结果，获得减噪声后的PDP。
通过上式得到了噪声的粗估计，每次平滑滤波前的PDP向量PDP_Current'(n)减去
使得PDP向量去除噪声的影响，然后再进行连续的alpha滤波，使得最终的PDP更加精确。即，
其中，PDP_Current(n)表示当前子帧减噪声后的PDP。
步骤403、利用减噪声后的PDP和前一平滑滤波处理后的PDP平滑值，获得当前PDP平滑值。
在本发明实施例中，通过alpha滤波获得当前PDP平滑值。
具体的，根据以下公式，利用减噪声后的PDP和前一平滑滤波和减噪声处理后的PDP平滑值，获得当前PDP平滑值，包括：
PDP_Smooth(n)＝(1-α)PDP_Smooth(n-1)+αPDP_Current(n)
其中，PDP_Smooth(n)表示当前PDP平滑值，PDP_Smooth(n-1)表示前一平滑滤波处理后的PDP平滑值，PDP_Current(n)表示当前子帧减噪声后的PDP，α表示平滑因子。对不同的估计值，α取值可能不同。
通过此过程获得的结果除了可用于下一步的计算，还可进行其他测量项测量如定时、频偏等。
步骤204、根据所述IFFT点数和所述当前系统带宽确定FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
在本发明实施例中，与现有技术不同的是，FFT点数并不完全取OFDM系统的FFT长度，而是联合系统带宽和IFFT点数确定FFT点数。在此，设置两档带宽门限，分别为第二、三带宽门限ThB2，ThB3(ThB2＜ThB3)用以判断，设置三个可配因子λ1,λ2,λ3，具体方案如图5所示：
步骤501、判断所述当前系统带宽是否小于第二带宽门限值。
步骤502、若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积。
步骤503、若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值。
步骤504、若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；
步骤505、若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积。
步骤506、利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
在利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH时，首先将所述PDP进行补零操作，利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，其中λ1,λ2,λ3可配，主要结合当前系统的采样率确定。
步骤205、计算所述PDP的采样率和系统采样率。
确定所述PDP的采样率为IFFT的长度和第四系数的乘积，确定所述系统采样率为OFDM调制时的FFT的长度和第五系数的乘积。
具体的，对PDP的采样率fspdp和系统采样率fs进行计算。对于LTE系统的基于CRS的估计，fspdp＝lIFFT×90kHz，fs＝lTFT×15kHz，lIFFT为PDP估计时的IFFT的长度，lTFT为OFDM调制时的FFT的长度。
步骤206、判断所述PDP的采样率是否和所述系统采样率相等。
判断PDP的采样率fspdp和系统采样率fs是否相等，以确定是否需要进行频域相关性间隔的校正。
步骤207、若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正。若相等，则频域MMSE滤波结束。
在此步骤中，计算所述PDP的采样率和所述系统采样率的比值，根据所述比值对所述RHH进行线性插值。
具体的，(1)计算PDP的采样率fspdp和系统采样率fs的比值。RHH_Δfcur表示当前的频域相关性向量的间隔，RHH_Δfnormal表示当前OFDM系统的子载波间隔，LTE中为15KHz。其中，
(2)根据
进行插值，化简
(n,m为正整数)。
在此，进行简单线性插值，每m个样点插值成n个样点，插值系数记为ε1,ε2,…εm，部分频域相关性向量RHH1,RHH2,…RHHn中，每个元素由原部分向量RHH1,RHH2,…RHHm和ε1,ε2,…εm加权求和组成，ε1,ε2,…εm需要做归一化处理。
步骤208、根据插值修正后的RHH，获取频域滤波系数。
之后，再利用获取的频域滤波系数进行后续的信道估计。
利用本发明实施例的方案，可使得RHH的间隔精确达到了15KHz，解决了现有技术中RHH间隔不准确造成频域相关值不准确进而导致性能恶化的问题；使得PDP估计中的IDFT的点数得以缩减，降低了复杂度；提升了PDP对于噪声的抑制效果，改善了低信噪比下的性能。
在具体应用中，将频域相关值记为RHH，RHH间隔表示的是频域相关性的间隔，因为做IDFT时导频的RE紧密排列，所以CRS RE频率间隔为90kHz(6×15KHz)。RHH间隔Δf计算公式如下：
PDP的采样率计算公式如下：
fs＝lIDFT×90k
那么，在现有技术中，
1)如果系统带宽为20MHz，则lIDFT＝324，lDFT＝2048，PDP的采样率为29.16MHz，而基带采样率为30.72MHz，RHH间隔14.238kHz；
2)如果系统带宽为15MHz，则lIDFT＝256，lDFT＝1536，PDP的采样率为23.04MHz，而基带采样率为23.04MHz，RHH间隔15kHz；
3)如果系统带宽为10MHz，则lIDFT＝180，lDFT＝1024，PDP的采样率为16.2MHz，而基带采样率为15.36MHz，RHH间隔15.82kHz；
4)如果系统带宽为5MHz，则lIDFT＝96，lDFT＝512，PDP的采样率为8.64MHz，而基带采样率为7.68MHz，RHH间隔16.875kHz；
5)如果系统带宽为3MHz，则lIDFT＝48，lDFT＝256，PDP的采样率为4.32MHz，而基带采样率为3.84MHz，RHH间隔16.875kHz；
6)如果带宽为1.4MHz，则lIDFT＝24，lDFT＝128，PDP的采样率为2.16MHz，而基带采样率为1.92MHz，RHH间隔16.875kHz。
具体如表1所示。
表1
如果，需要的频域相关性的频域间隔为15kHz，但是现在算法得到的间隔与15kHz最小差了0.82kHz，最大的5M以下相差了1.875kHz。此偏差会导致性能恶化。
而在本发明方案中，考虑IFFT远比IDFT速率快，复杂度低，所以点数选取2的整数次幂。将原方案的IDFT使用IFFT替代，并且缩减了IDFT(IFFT)的点数，进一步降低了复杂度。并调整lDFT的点数，并采用本专利的方法进行处理，从而保持RHH的间隔Δf仍然为15KHz。
1)如果系统带宽为20MHz，则lIFFT＝256，子载波间隔90kHz，lFFT＝1536，PDP采样率23.04MHz，基带采样率为23.04MHz，RHH间隔Δf15kHz；
如果此时lFFT＝2048，则PDP采样率23.04MHz，基带采样率为30.72MHz，RHH间隔Δf11.25kHz；
2)如果系统带宽为15MHz，则lIFFT＝256，子载波间隔90kHz，lFFT＝1536，PDP采样率23.04MHz，基带采样率为23.04MHz，RHH间隔Δf15kHz；
如果lFFT＝2048，PDP采样率23.04MHz，基带采样率为30.72MHz，RHH间隔Δf11.52kHz；
3)如果系统带宽为10MHz，则lIFFT＝128，子载波间隔90kHz，lFFT＝1024，PDP采样率11.52MHz，基带采样率为15.36MHz，RHH间隔Δf11.25kHz；
4)如果系统带宽为5MHz，则lIFFT＝64，子载波间隔90kHz，lFFT＝512，PDP采样率5.76MHz，基带采样率为7.68MHz，RHH间隔Δf11.25kHz；
5)如果系统带宽为3MHz，则lIFFT＝128，子载波间隔30kHz(IFFT输入，每2个RS RE之间插入1个0。频域插0，时域周期重复，PDP估计只有1个周期可用)，lFFT＝256，PDP采样率3.84MHz，基带采样率为3.84MHz，RHH间隔Δf15kHz；
6)如果系统带宽为1.4MHz，则lIFFT＝64，子载波间隔30kHz(IFFT输入，每2个RSRE之间插入1个0。频域插0，时域周期重复，PDP估计只有1个周期可用)，lFFT＝128，PDP采样率1.92MHz，基带采样率为1.92MHz，RHH间隔Δf15kHz。
具体如表2所示。
表2
有些系统带宽下，RHH间隔只有11.25kHz，需要插值到15kHz的频域间隔RHH’。11.25kHz为15kHz的3/4，每4个RHH可以插值出3个RHH’，所以可以用图6所示的方式，进行简单的线性插值得到。
通过以上可以看到，即使部分带宽下增加了线性插值的复杂度，但是总体而言，IFFT代替IDFT使得复杂度大为降低，同时IFFT点数相比于IDFT进行了缩减，进一步带来复杂度的降低。同时，RHH的间隔不再是近似，而是直接等于15KHz，进而消除了频域相关值不准确带来的性能恶化。
如图7所示，本发明实施例的信道估计装置，包括：
第一获取模块701，用于获取当前系统带宽；
第一处理模块702，用于根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；
第二处理模块703，用于根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；
第二获取模块704，用于根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
如图8所示，为了进一步提高RRH间隔的准确性，所述装置还可包括：
计算模块705，用于计算所述PDP的采样率和系统采样率；
判断模块706，用于判断所述PDP的采样率是否和所述系统采样率相等；
插值修正模块707，用于若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正；
所述第二获取模块704具体用于，根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，计算模块705具体用于，确定所述PDP的采样率为IFFT的长度和第四系数的乘积，确定所述系统采样率为正交频分复用OFDM调制时的FFT的长度和第五系数的乘积。插值修正模块707具体用于，计算所述PDP的采样率和所述系统采样率的比值，根据所述比值对所述RHH进行线性插值。
如图9所示，为了进一步提高RRH间隔的准确性，所述装置还可包括：
滤波和减噪声模块708，用于对所述PDP进行平滑滤波和减噪声处理；
所述第二处理模块703具体用于，根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述滤波和减噪声模块708，具体用于进行噪声测量，获得噪声测量结果，将所述PDP减去所述噪声测量结果，获得减噪声后的PDP，利用减噪声后的PDP和前一平滑滤波处理后的PDP平滑值，获得当前PDP平滑值。
具体的，根据以下公式，利用减噪声后的PDP和前一平滑滤波和减噪声处理后的PDP平滑值，获得当前PDP平滑值，包括：
PDP_Smooth(n)＝(1-α)PDP_Smooth(n-1)+αPDP_Current(n)
其中，PDP_Smooth(n)表示当前PDP平滑值，PDP_Smooth(n-1)表示前一平滑滤波处理后的PDP平滑值，PDP_Current(n)表示当前子帧减噪声后的PDP，α表示平滑因子。
其中，所述第一处理模块702包括：判断子模块，用于判断所述当前系统带宽是否大于第一带宽门限值；确定子模块，用于若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；处理子模块，用于利用IFFT点数进行IFFT运算，估计PDP。
其中，所述第二处理模块703包括：
第一判断子模块，用于判断所述当前系统带宽是否小于第二带宽门限值；第一确定子模块，用于若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；第二判断子模块，用于若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；第二确定子模块，用于若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；处理子模块，用于利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
具体的，所述处理子模块，将所述PDP进行补零操作，利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
本发明所述装置的工作原理可参照前述方法实施例的描述。
在本发明实施例中，根据当前系统带宽确定的IFFT点数、根据IFFT点数和当前系统带宽确定的FFT点数，可使得在获得频域相关性项链RHH时的间隔更为准确，使得频域MMSE滤波的系数更为准确，提高了系统性能。
如图10所示，本发明实施例的通信设备，包括：处理器1000，用于读取存储器1020中的程序，执行下列过程：
获取当前系统带宽；根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计；
收发机1010，用于在处理器1000的控制下接收和发送数据。
其中，在图10中，总线架构可以包括任意数量的互联的总线和桥，具体由处理器1000代表的一个或多个处理器和存储器1020代表的存储器的各种电路链接在一起。总线架构还可以将诸如外围设备、稳压器和功率管理电路等之类的各种其他电路链接在一起，这些都是本领域所公知的，因此，本文不再对其进行进一步描述。总线接口提供接口。收发机1010可以是多个元件，即包括发送机和收发机，提供用于在传输介质上与各种其他装置通信的单元。处理器1000负责管理总线架构和通常的处理，存储器1020可以存储处理器1000在执行操作时所使用的数据。
处理器1000负责管理总线架构和通常的处理，存储器1020可以存储处理器1000在执行操作时所使用的数据。
处理器1000还用于读取所述计算机程序，执行如下步骤:
计算所述PDP的采样率和系统采样率；
判断所述PDP的采样率是否和所述系统采样率相等；
若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正。
处理器1000还用于读取所述计算机程序，执行如下步骤:
根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
处理器1000还用于读取所述计算机程序，执行如下步骤:
对所述PDP进行平滑滤波和减噪声处理；
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
处理器1000还用于读取所述计算机程序，执行如下步骤:
判断所述当前系统带宽是否大于第一带宽门限值；
若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；
利用IFFT点数进行IFFT运算，估计PDP。
处理器1000还用于读取所述计算机程序，执行如下步骤:
判断所述当前系统带宽是否小于第二带宽门限值；
若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；
若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；
若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；
利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
处理器1000还用于读取所述计算机程序，执行如下步骤:
将所述PDP进行补零操作；
利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
处理器1000还用于读取所述计算机程序，执行如下步骤:
确定所述PDP的采样率为IFFT的长度和第四系数的乘积；
确定所述系统采样率为正交频分复用OFDM调制时的FFT的长度和第五系数的乘积。
处理器1000还用于读取所述计算机程序，执行如下步骤:
计算所述PDP的采样率和所述系统采样率的比值；
根据所述比值对所述RHH进行线性插值。
处理器1000还用于读取所述计算机程序，执行如下步骤:
进行噪声测量，获得噪声测量结果；
将所述PDP减去所述噪声测量结果，获得减噪声后的PDP；
利用减噪声后的PDP和前一平滑滤波处理后的PDP平滑值，获得当前PDP平滑值。
处理器1000还用于读取所述计算机程序，执行如下步骤:
根据以下公式，利用减噪声后的PDP和前一平滑滤波和减噪声处理后的PDP平滑值，获得当前PDP平滑值，包括：
PDP_Smooth(n)＝(1-α)PDP_Smooth(n-1)+αPDP_Current(n)
其中，PDP_Smooth(n)表示当前PDP平滑值，PDP_Smooth(n-1)表示前一平滑滤波处理后的PDP平滑值，PDP_Current(n)表示当前子帧减噪声后的PDP，α表示平滑因子。
此外，本发明实施例的计算机可读存储介质，用于存储计算机程序，所述计算机程序可被处理器执行实现以下步骤：
获取当前系统带宽；
根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；
根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，在所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH之后，所述方法还包括：
计算所述PDP的采样率和系统采样率；
判断所述PDP的采样率是否和所述系统采样率相等；
若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正。
其中，所述根据所述频域相关性向量RHH，，并利用所述频域滤波系数进行信道估计，包括：
根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
其中，在所述根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP之后，所述方法还包括：
对所述PDP进行平滑滤波和减噪声处理；
所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP，包括：
判断所述当前系统带宽是否大于第一带宽门限值；
若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；
利用IFFT点数进行IFFT运算，估计PDP。
其中，所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
判断所述当前系统带宽是否小于第二带宽门限值；
若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；
若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；
若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；
利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
将所述PDP进行补零操作；
利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
其中，所述计算所述PDP的采样率和系统采样率，包括：
确定所述PDP的采样率为IFFT的长度和第四系数的乘积；
确定所述系统采样率为正交频分复用OFDM调制时的FFT的长度和第五系数的乘积。
其中，所述若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正，包括：
计算所述PDP的采样率和所述系统采样率的比值；
根据所述比值对所述RHH进行线性插值。
其中，所述对所述PDP进行平滑滤波和减噪声处理，包括：
进行噪声测量，获得噪声测量结果；
将所述PDP减去所述噪声测量结果，获得减噪声后的PDP；
利用减噪声后的PDP和前一平滑滤波处理后的PDP平滑值，获得当前PDP平滑值。
其中，根据以下公式，利用减噪声后的PDP和前一平滑滤波和减噪声处理后的PDP平滑值，获得当前PDP平滑值，包括：
PDP_Smooth(n)＝(1-α)PDP_Smooth(n-1)+αPDP_Current(n)
其中，PDP_Smooth(n)表示当前PDP平滑值，PDP_Smooth(n-1)表示前一平滑滤波处理后的PDP平滑值，PDP_Current(n)表示当前子帧减噪声后的PDP，α表示平滑因子。
在本申请所提供的几个实施例中，应该理解到，所揭露方法和装置，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理包括，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元，可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述收发方法的部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-Only Memory，简称ROM)、随机存取存储器(Random Access Memory，简称RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (25)
1.一种信道估计方法，其特征在于，包括：
获取当前系统带宽；
根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；
根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计；
在所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH之后，所述方法还包括：
计算所述PDP的采样率和系统采样率；
判断所述PDP的采样率是否和所述系统采样率相等；
若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正。
2.根据权利要求1所述的方法，其特征在于，所述根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计包括：
根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
3.根据权利要求1-2任一项所述的方法，其特征在于，在所述根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP之后，所述方法还包括：
对所述PDP进行平滑滤波和减噪声处理；
所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
4.根据权利要求1所述的方法，其特征在于，所述根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP，包括：
判断所述当前系统带宽是否大于第一带宽门限值；
若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；
利用IFFT点数进行IFFT运算，估计PDP。
5.根据权利要求1所述的方法，其特征在于，所述根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
判断所述当前系统带宽是否小于第二带宽门限值；
若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；
若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；
若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；
利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
6.根据权利要求5所述的方法，其特征在于，所述利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH，包括：
将所述PDP进行补零操作；
利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
7.根据权利要求1所述的方法，其特征在于，所述计算所述PDP的采样率和系统采样率，包括：
确定所述PDP的采样率为IFFT的长度和第四系数的乘积；
确定所述系统采样率为正交频分复用OFDM调制时的FFT的长度和第五系数的乘积。
8.根据权利要求1所述的方法，其特征在于，所述若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正，包括：
计算所述PDP的采样率和所述系统采样率的比值；
根据所述比值对所述RHH进行线性插值。
9.根据权利要求3所述的方法，其特征在于，所述对所述PDP进行平滑滤波和减噪声处理，包括：
进行噪声测量，获得噪声测量结果；
将所述PDP减去所述噪声测量结果，获得减噪声后的PDP；
利用减噪声后的PDP和前一平滑滤波处理后的PDP平滑值，获得当前PDP平滑值。
10.根据权利要求9所述的方法，其特征在于，根据以下公式，利用减噪声后的PDP和前一平滑滤波和减噪声处理后的PDP平滑值，获得当前PDP平滑值，包括：
PDP_Smooth(n)＝(1-α)PDP_Smooth(n-1)+αPDP_Current(n)
其中，PDP_Smooth(n)表示当前PDP平滑值，PDP_Smooth(n-1)表示前一平滑滤波处理后的PDP平滑值，PDP_Current(n)表示当前子帧减噪声后的PDP，α表示平滑因子。
11.一种信道估计装置，其特征在于，包括：
第一获取模块，用于获取当前系统带宽；
第一处理模块，用于根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；
第二处理模块，用于根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；
第二获取模块，用于根据所述频域相关性向量RHH，获取频域滤波系数，利用所述频域滤波系数进行信道估计；
所述装置还包括：
计算模块，用于计算所述PDP的采样率和系统采样率；
判断模块，用于判断所述PDP的采样率是否和所述系统采样率相等；
插值修正模块，用于若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正；
所述第二获取模块具体用于，根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
12.根据权利要求11所述的装置，其特征在于，所述装置还包括：
滤波和减噪声模块，用于对所述PDP进行平滑滤波和减噪声处理；
所述第二处理模块具体用于，根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
13.根据权利要求11所述的装置，其特征在于，所述第一处理模块包括：
判断子模块，用于判断所述当前系统带宽是否大于第一带宽门限值；
确定子模块，用于若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；
处理子模块，用于利用IFFT点数进行IFFT运算，估计PDP。
14.根据权利要求11所述的装置，其特征在于，所述第二处理模块包括：
第一判断子模块，用于判断所述当前系统带宽是否小于第二带宽门限值；
第一确定子模块，用于若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；
第二判断子模块，用于若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；
第二确定子模块，用于若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；
处理子模块，用于利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
15.一种通信设备，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序；其特征在于，所述处理器执行所述程序时实现以下步骤：
获取当前系统带宽；
根据所述当前系统带宽确定快速傅里叶逆变换IFFT点数，并利用IFFT点数进行IFFT运算，估计功率时延谱PDP；
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH；
根据所述频域相关性向量RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计；
所述处理器还用于读取存储器中的程序，执行下列过程：
计算所述PDP的采样率和系统采样率；
判断所述PDP的采样率是否和所述系统采样率相等；
若PDP的采样率和所述系统采样率不相等，则对所述RHH进行插值修正。
16.根据权利要求15所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
根据插值修正后的RHH，获取频域滤波系数，并利用所述频域滤波系数进行信道估计。
17.根据权利要求15-16任一项所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
对所述PDP进行平滑滤波和减噪声处理；
根据所述IFFT点数和所述当前系统带宽确定快速傅里叶变换FFT点数，并利用FFT点数对进行平滑滤波和减噪声处理后的PDP进行FFT运算，获得频域相关性向量RHH。
18.根据权利要求15所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
判断所述当前系统带宽是否大于第一带宽门限值；
若所述当前系统带宽大于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次；若所述当前系统带宽小于或等于所述第一带宽门限值，则IFFT点数为实际信道估计向量长度的最小2的整数幂次与第一系数的乘积；
利用IFFT点数进行IFFT运算，估计PDP。
19.根据权利要求15所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
判断所述当前系统带宽是否小于第二带宽门限值；
若所述当前系统带宽小于所述第二带宽门限值，则FFT点数为所述IFFT点数和第三系数的乘积；
若所述当前系统带宽大于或等于所述第二带宽门限值，判断所述当前系统带宽是否大于第三带宽门限值，其中，所述第二带宽门限值小于所述第三带宽门限值；
若所述当前系统带宽大于或等于所述第三带宽门限值，则FFT点数为所述IFFT点数和第一系数的乘积；若所述当前系统带宽小于所述第三带宽门限值，则FFT点数为所述IFFT点数和第二系数的乘积；
利用FFT点数对所述PDP进行FFT运算，获得频域相关性向量RHH。
20.根据权利要求19所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
将所述PDP进行补零操作；
利用FFT点数对进行补零操作后的PDP进行FFT运算，获得频域相关性向量RHH。
21.根据权利要求15所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
确定所述PDP的采样率为IFFT的长度和第四系数的乘积；
确定所述系统采样率为正交频分复用OFDM调制时的FFT的长度和第五系数的乘积。
22.根据权利要求15所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
计算所述PDP的采样率和所述系统采样率的比值；
根据所述比值对所述RHH进行线性插值。
23.根据权利要求17所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
进行噪声测量，获得噪声测量结果；
将所述PDP减去所述噪声测量结果，获得减噪声后的PDP；
利用减噪声后的PDP和前一平滑滤波处理后的PDP平滑值，获得当前PDP平滑值。
24.根据权利要求23所述的设备，其特征在于，所述处理器还用于读取存储器中的程序，执行下列过程：
根据以下公式，利用减噪声后的PDP和前一平滑滤波和减噪声处理后的PDP平滑值，获得当前PDP平滑值，包括：
PDP_Smooth(n)＝(1-α)PDP_Smooth(n-1)+αPDP_Current(n)
其中，PDP_Smooth(n)表示当前PDP平滑值，PDP_Smooth(n-1)表示前一平滑滤波处理后的PDP平滑值，PDP_Current(n)表示当前子帧减噪声后的PDP，α表示平滑因子。
25.一种计算机可读存储介质，用于存储计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至10中任一项所述的方法中的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711171479.5A CN109818885B (zh) | 2017-11-22 | 2017-11-22 | 一种信道估计方法、装置、设备及计算机可读存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711171479.5A CN109818885B (zh) | 2017-11-22 | 2017-11-22 | 一种信道估计方法、装置、设备及计算机可读存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109818885A CN109818885A (zh) | 2019-05-28 |
| CN109818885B true CN109818885B (zh) | 2020-07-31 |
Family
ID=66601030
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711171479.5A Active CN109818885B (zh) | 2017-11-22 | 2017-11-22 | 一种信道估计方法、装置、设备及计算机可读存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109818885B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116055263B (zh) * | 2023-03-06 | 2023-06-13 | 南京创芯慧联技术有限公司 | 信道估计方法、装置、通信设备及存储介质 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101682588A (zh) * | 2007-03-27 | 2010-03-24 | 摩托罗拉公司 | 用于ofdm系统的信道估计器 |
| CN101702696A (zh) * | 2009-11-25 | 2010-05-05 | 北京天碁科技有限公司 | 信道估计的实现方法和装置 |
| CN102932291A (zh) * | 2012-10-29 | 2013-02-13 | 合肥东芯通信股份有限公司 | 一种lte系统下行链路rs的信道估计方法及装置 |
| CN104052555A (zh) * | 2014-06-19 | 2014-09-17 | 北京交通大学 | 一种ofdm系统下无线信道多径参数估计的方法 |
| CN105024951A (zh) * | 2014-04-28 | 2015-11-04 | 电信科学技术研究院 | 一种功率时延谱pdp估计方法及装置 |
| EP3110095A1 (en) * | 2015-06-26 | 2016-12-28 | Sequans Communications S.A. | Ambiguity resolution in positioning measurements |
| CN106789773A (zh) * | 2017-01-23 | 2017-05-31 | 北京松果电子有限公司 | 高速移动场景下的频偏和多普勒扩展的确定方法及装置 |
-
2017
- 2017-11-22 CN CN201711171479.5A patent/CN109818885B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101682588A (zh) * | 2007-03-27 | 2010-03-24 | 摩托罗拉公司 | 用于ofdm系统的信道估计器 |
| CN101702696A (zh) * | 2009-11-25 | 2010-05-05 | 北京天碁科技有限公司 | 信道估计的实现方法和装置 |
| CN102932291A (zh) * | 2012-10-29 | 2013-02-13 | 合肥东芯通信股份有限公司 | 一种lte系统下行链路rs的信道估计方法及装置 |
| CN105024951A (zh) * | 2014-04-28 | 2015-11-04 | 电信科学技术研究院 | 一种功率时延谱pdp估计方法及装置 |
| CN104052555A (zh) * | 2014-06-19 | 2014-09-17 | 北京交通大学 | 一种ofdm系统下无线信道多径参数估计的方法 |
| EP3110095A1 (en) * | 2015-06-26 | 2016-12-28 | Sequans Communications S.A. | Ambiguity resolution in positioning measurements |
| CN106789773A (zh) * | 2017-01-23 | 2017-05-31 | 北京松果电子有限公司 | 高速移动场景下的频偏和多普勒扩展的确定方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 基于导频的OFDM系统信道估计算法研究;李鑫;《中国优秀硕士学位论文全文数据库 信息科技辑》;20150531;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109818885A (zh) | 2019-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10419137B2 (en) | Estimating the frequency response of multipath channels | |
| TWI422193B (zh) | 通道估測裝置與方法 | |
| CN101909024B (zh) | 最大多普勒频偏的估计方法和装置 | |
| WO2016019657A1 (zh) | 一种频偏估计的方法和装置、存储介质 | |
| JP2017188890A (ja) | アンチエイリアスチャネル推定装置、方法及び受信機 | |
| CN114095103B (zh) | 频谱平坦度的校正方法与装置 | |
| CN109818885B (zh) | 一种信道估计方法、装置、设备及计算机可读存储介质 | |
| CN103581065B (zh) | 一种维纳滤波信道估计方法和装置 | |
| CN104753834B (zh) | 一种信道估计方法和装置 | |
| KR101018140B1 (ko) | Ｏｆｄｍ 시스템의 샘플링 주파수 동기 장치 | |
| CN112583753B (zh) | 一种相位补偿方法及电子设备 | |
| US20110064170A1 (en) | Apparatus and method for estimating a channel in a broadband wireless communication system | |
| EP2941852A1 (en) | Method for estimating and removing dc offset | |
| CN111181886A (zh) | 频偏估计方法和装置 | |
| JP6050177B2 (ja) | ノイズ電力推定方法 | |
| EP2328310B1 (en) | Method for estimating a received signal and corresponding device | |
| CN117880021B (zh) | Ofdm信号的直流估计方法及装置 | |
| CN108075992B (zh) | 一种信道估计方法及装置 | |
| US20180183527A1 (en) | Channel estimation using intra-symbol frequency domain averaging | |
| EP2533441B1 (en) | Channel estimate interpolation circuit and channel estimate interpolation method | |
| CN112804170B (zh) | 直流偏移处理方法及接收机 | |
| EP2685686B1 (en) | Method and apparatus for channel estimation based on estimated auto correlation | |
| TWI607633B (zh) | 無線通訊之雜訊變異量估測電路與方法 | |
| CA2915645A1 (en) | Method of estimating the frequency of response of multipath channels | |
| JP6720595B2 (ja) | マルチキャリア信号のクリッピング歪みの推定装置、補償装置及び受信機 |
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 | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road Patentee after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY |
|
| CP01 | Change in the name or title of a patent holder | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210601 Address after: 100085 1st floor, building 1, yard 5, Shangdi East Road, Haidian District, Beijing Patentee after: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd. Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY |
|
| TR01 | Transfer of patent right |