信噪比估计方法、装置、通信设备和存储介质与流程

1.本技术实施例属于通信技术领域，特别是涉及一种信噪比估计方法、装置、通信设备和存储介质。


背景技术：

2.信噪比(signal noise ratio，snr或s/n)是指电子设备、电子系统或信道中信号与噪声的比例，常常用分贝(decibel，db)数表示。一般来说，信噪比越大，表示混在信号里的噪声越小，声音回放的质量越高；否则，则相反。
3.在应用正交频分复用技术(orthogonal frequency division multiplexing，ofdm)的通信系统中，一般可以基于信道的功率时延分布(power
‑
delay profile，pdp)来进行信噪比估计。
4.但是，现有的基于信道的pdp来计算信噪比的方法误差较大。


技术实现要素：

5.有鉴于此，本技术实施例提供了一种信噪比估计方法、装置、通信设备和存储介质，用以提高信噪比估计的准确性。
6.本技术实施例的第一方面提供了一种信噪比估计方法，包括：
7.获取信道的功率时延分布；
8.将所述功率时延分布平均划分为多个分布段；
9.循环计算所述多个分布段的滑动平均值，得到滑动平均值集合，所述滑动平均值集合中的最小值用于表征噪声功率；
10.确定所述功率时延分布中用于统计信号功率的多个目标位置点；
11.基于所述多个目标位置点计算信号功率，所述信号功率等于所述多个目标位置点对应的功率之和减去保留径数与所述噪声功率之积所得的值；
12.基于所述信号功率与所述噪声功率，计算所述信道的信噪比。
13.本技术实施例的第二方面提供了一种信噪比估计装置，包括：
14.功率时延分布获取模块，用于获取信道的功率时延分布；
15.功率时延分布划分模块，用于将所述功率时延分布平均划分为多个分布段；
16.滑动平均值计算模块，用于循环计算所述多个分布段的滑动平均值，得到滑动平均值集合，所述滑动平均值集合中的最小值用于表征噪声功率；
17.目标位置点确定模块，用于确定所述功率时延分布中用于统计信号功率的多个目标位置点；
18.信号功率计算模块，用于基于所述多个目标位置点计算信号功率，所述信号功率等于所述多个目标位置点对应的功率之和减去保留径数与所述噪声功率之积所得的值；
19.信噪比计算模块，用于基于所述信号功率与所述噪声功率，计算所述信道的信噪比。
20.本技术实施例的第三方面提供了一种通信设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述第一方面所述的信噪比估计方法。
21.本技术实施例的第四方面提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上述第一方面所述的信噪比估计方法。
22.本技术实施例的第五方面提供了一种计算机程序产品，当所述计算机程序产品在通信设备上运行时，使得所述通信设备执行上述第一方面所述的信噪比估计方法。
23.与现有技术相比，本技术实施例包括以下优点：
24.本技术实施例，使用功率时延分布来估计噪声与信号功率。在实际处理过程中，通过使用分段求平均，然后再进行循环求平均的方式，减少了处理过程中的计算量。其次，对于获取到的功率时延分布中的最大值及其相邻的两个值，在抑噪时予以保留，保证了信道的主径不会丢失。第三，本技术实施例在处理过程中通过区分噪声、循环前缀和时偏范围，相较于现有技术中只是单纯地根据阈值进行抑噪效果更好。第四，在计算信号功率时，本技术实施例并非单纯地将多径相加，而是减去了每根径上的噪声，理论上的计算结果也更为准确。这样，最终计算得到的信噪比也更准确。
附图说明
25.为了更清楚地说明本技术实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍。显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
26.图1是本技术实施例提供的一种信噪比估计方法的步骤流程示意图；
27.图2是本技术实施例提供的一种信噪比估计方法的步骤s101的一种实现方式的示意图；
28.图3是本技术实施例提供的一种信噪比估计方法的步骤s103的一种实现方式的示意图；
29.图4是本技术实施例提供的一种分段求平均的示意图；
30.图5是本技术实施例提供的一种滑动求平均的示意图；
31.图6是本技术实施例提供的一种信噪比估计方法的步骤s104的一种实现方式的示意图；
32.图7是本技术实施例提供的一种信噪比估计装置的示意图；
33.图8是本技术实施例提供的一种通信设备的示意图。
具体实施方式
34.以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本技术实施例。然而，本领域技术人员应当清楚，在没有这些具体细节的其他实施例中也可以实现本技术。在其他情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本技术的描述。
35.基于信道的pdp来进行信噪比估计时，ofdm理想的信道冲击响应(channel impulse response，cir)可以如下式(1)表示：
[0036][0037]
其中，l为信号传播过程中所经过的多个路径的路径数量，a
l
为每个路径的幅值，n
l
为信号经每个路径的到达时间。
[0038]
由于系统中存在噪声，且cir的长度在大多情况下不会超过循环前缀(cyclic prefix，cp)的长度，cir也可以使用下式(2)表示：
[0039][0040]
其中，w(n)为噪声，cp内表示信号位于循环前缀内，相应地，cp外表示信号位于循环前缀外。
[0041]
因此，上述信道的pdp可以使用下式(3)表示：
[0042][0043]
由于信号能量一般比噪声大，通过将pdp分布中位于cp外的部分视为噪声，cp内的部分视为信道能量，可以计算得到该信道的信噪比。但是，采用这种方式计算得到的信噪比，往往误差较大。
[0044]
针对上述问题，本技术实施例提供了一种信噪比估计方法、装置、通信设备和存储介质，本方法同样基于pdp来进行信噪比的估计，但相较于现有技术，并非直接将位于cp外的部分视为噪声，cp内的部分视为信道能量，而是通过区分噪声、循环前缀和时偏范围来进行抑噪，保留pdp中的最大值及其相邻的两个值，保证了信道的主径不会丢失；功率计算时也并非单纯地将多径相加，而是减去了每根径上的噪声，使计算结果更准确。
[0045]
下面通过具体实施例来说明本技术的技术方案。
[0046]
参照图1，示出了本技术实施例提供的一种信噪比估计方法的步骤流程示意图，具体可以包括如下步骤：
[0047]
s101、获取信道的功率时延分布。
[0048]
需要说明的是，本方法可以应用于ofdm系统。具体地，本方法的执行主体可以是ofdm系统中具备信噪比估计功能的通信设备，本技术实施例对通信设备的具体类型不作限定。
[0049]
在本技术实施例中，通信设备可以使用pdp来估计噪声功率与信号功率。pdp又被称作功率时延谱，描述的是信道在时间上的色散，表示某一时延处接收信号功率的期望。
[0050]
在本技术实施例的一种可能的实现方式中，如图2所示，通信设备获取信道的pdp可以包括如下子步骤s1011
‑
s1013：
[0051]
s1011、根据最小二乘准则(least
‑
square，ls)，对所述信道进行信道估计。
[0052]
信道估计，就是从接收数据中将假定的某个信道模型的模型参数估计出来的过程。在实际应用中，通信设备可以采用多种方法进行信道估计，如非盲信道估计、盲信道估
计等等，本技术实施例对此不作限定。
[0053]
在本技术实施例中，通信设备可以采用ls信道估计的方式来进行信道估计，得到相应的结果。
[0054]
s1012、对所述信道估计的结果作快速傅里叶逆变换(inverse fast fourier transform，ifft)，得到变换结果。
[0055]
快速傅里叶变换(fast fourier transform，fft)是一种高效的离散傅里叶变换(discrete fourier transform，dft)算法，在通信领域具有广泛的应用。ifft即是上述fft算法的逆运算。
[0056]
s1013、依次对所述变换结果进行求模和平方运算，得到所述信道的功率时延分布。
[0057]
在本技术实施例中，根据时频域转换原理，通信设备可以对信道估计得到的结果进行ifft变换，然后对变换结果依次进行求模和平方运算，可以得到信道的pdp。
[0058]
作为一种示例，可以采用采用公式(4)、(5)表示上述获取pdp的过程：
[0059][0060][0061]
其中，h
ls
表示采用ls信道估计得到的结果，ifft(h
ls
)为ifft变换结果。
[0062]
需要说明的是，若信道的长度为k，则pdp的长度可以如下式(6)表示：
[0063][0064]
在本技术实施例中，在获取到pdp后，通信设备可以记录该pdp中的最大值pdp
max
及其位置点n
max
：
[0065]
pdp
max
＝max(pdp(n))
……
(7)
[0066]
n
max
＝argmax(pdp(n))
……
(8)
[0067]
通过记录上述pdp中的最大值及其位置点，通信设备在后续进行抑制噪声处理时，可以保留该最大值，保证其不会被去除。
[0068]
在本技术实施例的一种可能的实现方式中，除记录上述pdp中的最大值及其位置点外，通信设备还可以记录与该最大值相邻的位置点，例如，该最大值的前一个位置点n
max,l
和后一个位置点n
max,r
：
[0069]
n
max,l
＝(n
max
‑
1)&(n
pdp
‑
1)
……
(9)
[0070]
n
max,r
＝(n
max
+1)&(n
pdp
‑
1)
……
(10)
[0071]
其中，n
pdp
‑
1表示该位置点位于n
pdp
的范围之内。
[0072]
在后续进行抑制噪声处理时，通信设备可以保留上述pdp中的最大值的位置点及其相邻的位置点。这样，可以保证信道的主径不会丢失。
[0073]
s102、将所述功率时延分布平均划分为多个分布段。
[0074]
在本技术实施例，为了减少后续的计算量，通信设备可以将pdp平均划分为多个分布段。这样，每个分布段中包含的位置点数量相等。
[0075]
在具体实现中，pdp的长度为n
pdp
，可以将其平均划分为g
pdp
段，则每个分布段内包含的位置点的数量可表示为n
g
＝n
pdp
/g
pdp
。
[0076]
上述g
pdp
为分段常数，其具体大小可以根据实际需要确定。一般地，g
pdp
的取值应满足为2的n次方。示例性的，在n＝5时，g
pdp
＝32，即表示将pdp平均划分为32段。
[0077]
s103、循环计算所述多个分布段的滑动平均值，得到滑动平均值集合，所述滑动平均值集合中的最小值用于表征噪声功率。
[0078]
通常，pdp包含有噪声和多径。其中，多径分布在cp的长度范围内。可以设定cp所占的段数为g
cp
，该g
cp
的取值应当满足需要包含多径的位置范围。因此，噪声所占的段数可以相应地表示为g
noise
＝g
pdp
‑
g
cp
。
[0079]
除此之外，由于时偏的影响，pdp可能会产生偏移，可以设置偏移段数为g
ta
，该g
ta
的取值应当满足包含时偏的位置范围。
[0080]
在本技术实施例中，为了尝试将所有可能的噪声值相加，可以对已经平均划分的多个分布段进行滑动求平均计算，得到包含多个滑动平均值的滑动平均值集合。
[0081]
在本技术实施例的一种可能的实现方式中，如图3所示，通信设备通过循环计算多个分布段的滑动平均值，得到滑动平均值集合，可以包括如下子步骤s1031
‑
s1033：
[0082]
s1031、分别计算每个分布段的算术平均值。
[0083]
在本技术实施例中，为了方便滑动平均值的计算，可以首先计算每个分布段的算术平均值。计算每个分布段的算术平均值也就是分段求平均，针对每个分布段中n
g
个位置点，计算这些位置点的平均值。pdp被划分为多少个分布段，计算得到的算术平均值就有多少个。即，pdp被划分为g
pdp
段，则算术平均值也就有g
pdp
个。
[0084]
如图4所示，是本技术实施例提供的一种分段求平均的示意图。图4中示出了将pdp平均划分为g
pdp
＝32段(0,1,2,
…
,31)时的分段求平均的过程。pdp的长度为n
pdp
，每段内包含n
g
个位置点，最终得到的算术平均值包括32个。
[0085]
具体地，可以采用如下公式(11)计算每个分布段的算术平均值：
[0086][0087]
其中，s1(g1)为g
pdp
段的算术平均值集合，g
pdp
为分段常数，也就是划分后得到的多个分布段的总段数，g1为对应分段值，n
g
为每个分布段内包含的位置点的数量。
[0088]
s1032、确定滑动起始段和滑动结束段。
[0089]
确定滑动起始段和滑动结束段是指确定滑动求平均值时，应当从哪一段开始，到哪一段结束。
[0090]
在本技术实施例中，滑动起始段可以表示为g
s
＝mod(
‑
g
ta
,g
pdp
)，滑动结束段可以表示为g
e
＝g
cp
+g
ta
‑
1。其中，mod为求余预算，g
ta
为偏移段数，g
pdp
为多个分布段的总段数，g
cp
为循环前缀所占段数。
[0091]
s1033、基于所述每个分布段的算术平均值，以预定段数为长度，计算在所述滑动起始段至所述滑动结束段范围内的滑动平均值，得到所述滑动平均值集合。
[0092]
在本技术实施例中，预定段数可以与噪声所占的段数g
noise
相等，也就等于多个分布段的总段数减去循环前缀所占段数。
[0093]
如图5所示，是本技术实施例提供的一种滑动求平均的示意图。按照图5所示，滑动求平均的过程，也就是从滑动起始段g
s
开始，每次计算长度为g
noise
的多个分布段的平均值，
每次计算完成后，向后滑动一段，再次计算长度为g
noise
的多个分布段的平均值，直到计算到从滑动结束段g
e
开始的长度为g
noise
的多个分布段的平均值。
[0094]
具体地，可以采用如下公式(12)计算滑动平均值s2(g2)：
[0095][0096]
上述s2(g2)为滑动平均值的集合。
[0097]
s104、确定所述功率时延分布中用于统计信号功率的多个目标位置点。
[0098]
用于统计信号功率的多个目标位置点可以看作是计算信号功率时不会被去除的位置点，而需要被去除的位置点可以认为是pdp上的噪声。本技术实施例中的多个目标位置点可以包括pdp中属于cp的位置点和该pdp中最大值的位置点及其相邻的位置点。
[0099]
在本技术实施例的一种可能的实现方式中，如图6所示，通信设备确定pdp中用于统计信号功率的多个目标位置点，可以包括如下子步骤s1041
‑
s1042：
[0100]
s1041、根据所述最小值在所述滑动平均值集合中的位置，确定所述功率时延分布中属于所述循环前缀的位置索引集合。
[0101]
在本技术实施例中，可以认为最小的滑动求平均结果中只包含噪声。因此，上述s2(g2)集合中的最小值可以用于表征噪声功率：
[0102]
p
noise
＝min(s2(g2))
……
(13)
[0103]
那么，其位置g
min
可以表示为：
[0104][0105]
这样，就可以根据噪声段的位置，去反推cp所在的位置。
[0106]
在具体实现中，通信设备可以根据最小值在滑动平均值集合中的位置，确定噪声段的起始位置g
start
。
[0107]
一般地，若g
min
＞g
ta
，则噪声段的起始位置可以表示为g
start
＝g
min
‑
g
ta
‑
1；若g
min
≤g
ta
，则噪声段的起始位置可以表示为g
start
＝g
pdp
+g
min
‑
g
ta
‑
1；
[0108]
其中，g
min
为滑动平均值集合中最小值的位置，g
ta
为偏移段数，g
pdp
为多个分布段的总段数。
[0109]
需要说明的是，由于时偏的影响，pdp分布可能会产生偏移。因此，若g
min
＞g
ta
，则表示实际噪声段的起始位置g
start
在pdp的前半段；若g
min
≤g
ta
，则表示实际噪声段的起始位置g
start
在pdp的后半段。
[0110]
然后，通信设备可以根据噪声段的起始位置g
start
，确定pdp中属于cp的位置索引集合idx
cp
为：
[0111]
idx
cp
＝(n
g
·
g
start
‑
n)&(n
pdp
‑
1)，n＝0,1,2,
…
,n
g
·
g
cp
‑1……
(15)
[0112]
其中，n
g
·
g
start
表示pdp上噪声的起始位置，n
g
·
g
cp
表示pdp上属于cp段的位置点的个数，n
pdp
‑
1表示位置索引在n
pdp
的范围之内。
[0113]
s1042、将属于所述位置索引集合中的多个位置点以及所述功率时延分布中最大值的位置点及其相邻的位置点，确定为用于统计信号功率的多个目标位置点。
[0114]
在本技术实施例中，可以将属于idx
cp
集合中的多个位置点以及pdp中最大值的位
置点n
max
及其相邻的位置点n
max,l
和n
max,r
，确定为用于统计信号功率的多个目标位置点。
[0115]
在具体实现中，可以将cp外的位置点置零，并保留上述最大值的位置点n
max
及其相邻的位置点n
max,l
和n
max,r
。
[0116]
在本技术实施例的一种可能的实现方式中，对于上述目标位置点，还可以将idx
cp
集合中的多个位置点中功率小于阈值的位置点删除，上述阈值等于噪声功率p
noise
的预设倍数。
[0117]
示例性的，可以根据p
noise
确定如下阈值thr：
[0118]
thr＝r
·
p
noise
……
(16)
[0119]
其中，上述r可以为经验值。
[0120]
在具体实现中，对于上述功率小于阈值的位置点，可以将其置为零。
[0121]
经过上述处理，pdp中位于cp外和小于阈值的位置点均被置零，保留了cp内且大于阈值的位置点n
max
及其相邻的位置点n
max,l
和n
max,r
。通信设备可以根据被保留的目标位置点计算信号功率。
[0122]
上述处理过程可以以如下代码段表示：
[0123][0124]
其中，n
path
为保留径数。
[0125]
s105、基于所述多个目标位置点计算信号功率，所述信号功率等于所述多个目标位置点对应的功率之和减去保留径数与所述噪声功率之积所得的值。
[0126]
在本技术实施例中，在计算信号功率时，可以在将多个路径的信号功率相加的基础上，减去每个路径上的噪声。这样，计算得到的信号相较于现有技术中单纯地只是将多径上的信号功率相加，更为准确。
[0127]
在具体实现中，n
path
与多个目标位置点的数量相等。因此，将多个路径的信号功率相加可以是指将多个目标位置点对应的功率相加；需要被减去的每个路径上的噪声也就是n
path
与噪声功率p
noise
的乘积。
[0128]
因此，可以采用如下公式(17)计算pdp上的信号功率：
[0129][0130]
s106、采用所述信号功率与所述噪声功率，计算所述信道的信噪比。
[0131]
计算得到的信噪比可以表示为：
[0132][0133]
在本技术实施例中，使用功率时延分布来估计噪声与信号功率。并且，在实际处理过程中，通过使用分段求平均，然后再进行循环求平均的方式，减少了处理过程中的计算量。其次，对于获取到的功率时延分布中的最大值及其相邻的两个值，在抑噪时予以保留，保证了信道的主径不会丢失。第三，本技术实施例在处理过程中通过区分噪声、循环前缀和时偏范围，相较于现有技术中只是单纯地根据阈值进行抑噪效果更好。第四，在计算信号功率时，本技术实施例并非单纯地将多径相加，而是减去了每根径上的噪声，理论上的计算结果也更为准确。
[0134]
需要说明的是，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本技术实施例的实施过程构成任何限定。
[0135]
参照图7，示出了本技术实施例提供的一种信噪比估计装置的示意图，具体可以包括功率时延分布获取模块701、功率时延分布划分模块702、滑动平均值计算模块703、目标位置点确定模块704、信号功率计算模块705和信噪比计算模块706，其中：
[0136]
功率时延分布获取模块701，用于获取信道的功率时延分布；
[0137]
功率时延分布划分模块702，用于将所述功率时延分布平均划分为多个分布段；
[0138]
滑动平均值计算模块703，用于循环计算所述多个分布段的滑动平均值，得到滑动平均值集合，所述滑动平均值集合中的最小值用于表征噪声功率；
[0139]
目标位置点确定模块704，用于确定所述功率时延分布中用于统计信号功率的多个目标位置点；
[0140]
信号功率计算模块705，用于基于所述多个目标位置点计算信号功率，所述信号功率等于所述多个目标位置点对应的功率之和减去保留径数与所述噪声功率之积所得的值；
[0141]
信噪比计算模块706，用于基于所述信号功率与所述噪声功率，计算所述信道的信噪比。
[0142]
在本技术实施例中，保留径数与多个目标位置点的数量相等。
[0143]
在本技术实施例中，功率时延分布获取模块701具体用于：根据最小二乘准则，对所述信道进行信道估计；对所述信道估计的结果作快速傅里叶逆变换，得到变换结果；依次对所述变换结果进行求模和平方运算，得到所述信道的功率时延分布。
[0144]
在本技术实施例中，滑动平均值计算模块703具体用于：分别计算每个分布段的算术平均值；确定滑动起始段和滑动结束段；基于所述每个分布段的算术平均值，以预定段数为长度，计算在所述滑动起始段至所述滑动结束段范围内的滑动平均值，得到所述滑动平均值集合。
[0145]
在本技术实施例中，滑动起始段可以为g
s
＝mod(
‑
g
ta
,g
pdp
)，滑动结束段可以为g
e
＝g
cp
+g
ta
‑
1，预定段数与噪声所占段数相等，噪声所占段数等于多个分布段的总段数减去循环前缀所占段数；其中，mod为求余预算，g
ta
为偏移段数，g
pdp
为多个分布段的总段数，g
cp
为循环前缀所占段数。
[0146]
在本技术实施例中，目标位置点确定模块704具体用于：根据所述最小值在所述滑
动平均值集合中的位置，确定所述功率时延分布中属于所述循环前缀的位置索引集合；将属于所述位置索引集合中的多个位置点以及所述功率时延分布中最大值的位置点及其相邻的位置点，确定为用于统计信号功率的多个目标位置点。
[0147]
在本技术实施例中，目标位置点确定模块704还用于：根据所述最小值在所述滑动平均值集合中的位置，确定噪声段的起始位置；根据所述噪声段的起始位置，确定所述功率时延分布中属于所述循环前缀的位置索引集合。
[0148]
在本技术实施例中，若g
min
＞g
ta
，则所述噪声段的起始位置g
start
＝g
min
‑
g
ta
‑
1；若g
min
≤g
ta
，则所述噪声段的起始位置g
start
＝g
pdp
+g
min
‑
g
ta
‑
1，循环前缀的位置索引集合为：idx
cp
＝(n
g
·
g
start
‑
n)&(n
pdp
‑
1)，n＝0,1,2,
…
,n
g
·
g
cp
‑
1；其中，g
min
为所述最小值在所述滑动平均值集合中的位置，g
ta
为偏移段数，g
pdp
为所述多个分布段的总段数。
[0149]
在本技术实施例中，目标位置点确定模块704还用于：将所述位置索引集合中的多个位置点中功率小于阈值的位置点删除，所述阈值等于所述噪声功率的预设倍数。
[0150]
对于装置实施例而言，由于其与方法实施例基本相似，所以描述得比较简单，相关之处参见方法实施例部分的说明即可。
[0151]
参照图8，示出了本技术实施例提供的一种通信设备的示意图。如图8所示，该通信设备800包括：处理器810、存储器820以及存储在所述存储器820中并可在所述处理器810上运行的计算机程序821。所述处理器810执行所述计算机程序821时实现上述信噪比估计方法各个实施例中的步骤，例如图1所示的步骤s101至s106。或者，所述处理器810执行所述计算机程序821时实现上述各装置实施例中各模块/单元的功能，例如图7所示模块701至706的功能。
[0152]
示例性的，所述计算机程序821可以被分割成一个或多个模块/单元，所述一个或者多个模块/单元被存储在所述存储器820中，并由所述处理器810执行，以完成本技术。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段可以用于描述所述计算机程序821在所述通信设备800中的执行过程。例如，所述计算机程序821可以被分割成功率时延分布获取模块、功率时延分布划分模块、滑动平均值计算模块、目标位置点确定模块、信号功率计算模块和信噪比计算模块，各模块具体功能如下：
[0153]
功率时延分布获取模块，用于获取信道的功率时延分布；
[0154]
功率时延分布划分模块，用于将所述功率时延分布平均划分为多个分布段；
[0155]
滑动平均值计算模块，用于循环计算所述多个分布段的滑动平均值，得到滑动平均值集合，所述滑动平均值集合中的最小值用于表征噪声功率；
[0156]
目标位置点确定模块，用于确定所述功率时延分布中用于统计信号功率的多个目标位置点；
[0157]
信号功率计算模块，用于基于所述多个目标位置点计算信号功率，所述信号功率等于所述多个目标位置点对应的功率之和减去保留径数与所述噪声功率之积所得的值；
[0158]
信噪比计算模块，用于基于所述信号功率与所述噪声功率，计算所述信道的信噪比。
[0159]
所述通信设备800可以是基站或终端。其中终端包括手机、平板、个人计算机。所述通信设备800可包括，但不仅限于，处理器810、存储器820。本领域技术人员可以理解，图8仅仅是通信设备800的一种示例，并不构成对通信设备800的限定，可以包括比图示更多或更
少的部件，或者组合某些部件，或者不同的部件，例如所述通信设备800还可以包括输入输出设备、网络接入设备、总线等。
[0160]
所述处理器810可以是中央处理单元(central processing unit，cpu)，还可以是其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现成可编程门阵列(field
‑
programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
[0161]
所述存储器820可以是所述通信设备800的内部存储单元，例如通信设备800的硬盘或内存。所述存储器820也可以是所述通信设备800的外部存储设备，例如所述通信设备800上配备的插接式硬盘，智能存储卡(smart media card，smc)，安全数字(secure digital，sd)卡，闪存卡(flash card)等等。进一步地，所述存储器820还可以既包括所述通信设备800的内部存储单元也包括外部存储设备。所述存储器820用于存储所述计算机程序821以及所述通信设备800所需的其他程序和数据。所述存储器820还可以用于暂时地存储已经输出或者将要输出的数据。
[0162]
本技术实施例还公开了一种通信设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如前述各个实施例所述的信噪比估计方法。
[0163]
本技术实施例还公开了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如前述各个实施例所述的信噪比估计方法。
[0164]
本技术实施例还公开了一种计算机程序产品，当所述计算机程序产品在通信设备上运行时，使得所述通信设备执行前述各个实施例所述的信噪比估计方法。
[0165]
以上所述实施例仅用以说明本技术的技术方案，而非对其限制。尽管参照前述实施例对本技术进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本技术各实施例技术方案的精神和范围，均应包含在本技术的保护范围之内。