基于多用户SRS信号的SNR估计方法和装置与流程

基于多用户srs信号的snr估计方法和装置【
技术领域：
：】1.本发明涉及移动通信
技术领域：
：，尤其是涉及一种基于多用户srs信号的snr估计方法和装置。
背景技术：
：：2.在移动通信5g系统中，尤其是基站上行系统，需要根据srs(soundingreferencesignal，探测参考信号)进行snr(signal-to-noiseratio，信噪比)的测量。snr是指对接收信号中的信息和噪声功率或能量分别进行估计和测量，并计算信息和噪声之比。信噪比是衡量移动通信5g系统质量的重要指标。一方面，系统使用snr参数衡量信道质量，另一方面可以通过snr参数优化系统中其它算法模块的性能。随着高速通信系统的快速发展，对snr估计的要求也越来越高，要求准确度更高、性能更好、计算更简单且实现更容易。3.目前，5g小站系统计算snr估计方法为最大似然(maximumlikelihood)估计法，该方法是利用接收到的频域信号y，采用最小二乘法准则(leastsquare，ls)算法计算出信道估计响应h，然后利用y-hx(x为本地生成的srs序列信号)计算出信号噪声ni。最后计算出信噪比snr＝hx/(y-hx)。在实际测试中发现，在时偏较大的场景下，传统的ls最小二乘法则算法计算信道估计响应h不够精确，当小基站同时接入多个用户终端(ue)时，存在多ue的srs信号相互叠加干扰现象，从而导致snr测量不准确的问题。因此，信道估计响应计算越精确，snr计算越准确，对整个5g小基站系统性能提升具有至关重要的作用。4.因此，亟需提出一种信道估计更准确的snr估计方法和装置。技术实现要素：5.本发明所要解决的技术问题是提供基于多用户srs信号的snr估计方法和装置，其信道估计更准确。6.为解决上述技术问题，本发明提供一种基于多用户srs信号的snr估计方法，包括以下步骤：7.s1.从接收到的频域数据中提取出多个ue的srs测量信号ysrs(k，l，r)；其中，k为接收到的srs信号的子载波索引，l为ofdm符号，r为接收天线；8.s2.根据3gpp协议生成每个ue的本地发生序列xsrs(k，l，p)；其中，p为发送天线端口索引；9.s3.依据所述多个ue的srs测量信号ysrs(k，l，r)和所述每个ue的本地发生序列xsrs(k，l，p)，基于最小二乘估计算法，计算得到每个ue的粗信道估计10.s4.对所述每个ue的粗信道估计进行时偏预补偿，得到预补偿信道估计11.s5.对所述每个ue的预补偿信道估计进行连续nm子载波平滑去干扰处理，得到每个ue的中间信道估计12.其中，所述其中，所述是srs的端口数，nu是用户数，或8；13.s6.利用所述每个ue的中间信道估计进行时偏估计，得到时偏值ta；14.s7.依据所述时偏值ta对所述每个ue的中间信道估计进行时偏补偿，得到每个ue的补偿信道估计h′srs(k，l，r，p)；15.s8.依据mmse均衡算法对所述每个ue的补偿信道估计h′srs(k，l，r，p)进行插值滤波处理，得到不同时刻不同子载波间协方差矩阵θ(k)及mmse矩阵φ(k′)；16.s9.依据所述不同时刻不同子载波间协方差矩阵θ(k)和所述mmse矩阵φ(k′)计算权值w(k，l；k′，l′)，再依据所述权值w(k，l；k′，l′)对所述每个ue的补偿信道估计h′srs(k，l，r，p)进行re级插值运算，得到每个ue的插值信道估计h"srs(k，l，r，p)；17.s10.对所述每个ue的插值信道估计h"srs(k，l，r，p)进行解时偏补偿处理，得到每个ue的最终信道估计hsrs(k，l，r，p)；18.s11.依据所述每个ue的最终信道估计hsrs(k，l，r，p)和所述每个ue的本地发生序列xsrs(k，l，p)计算频带上的信号功率pu；依据所述srs测量信号ysrs(k，l，r)，所述最终信道估计hsrs(k，l，r，p)和所述本地发生序列xsrs(k，l，p)计算出频带上的噪声功率ni；19.s12.依据所述频带上的信号功率pu和所述频带上的噪声功率ni，基于信噪比计算公式确定中间信噪比snr′；20.s13.根据所述中间信噪比snr′和所述协方差矩阵θ(k)，得到新的mmse矩阵φ(k′)；然后，返回所述步骤s9进行再一次依序逐步计算，直到所述步骤s12依据所述信噪比计算公式得到最终信噪比snrest。21.更进一步的，所述步骤s3进一步包括：22.所述粗信道估计23.所述步骤s4进一步包括：24.所述预补偿信道估计其中，taest,pre＝40ts。25.更进一步的，所述步骤s5进一步包括：26.所述中间信道估计27.更进一步的，所述步骤s6进一步包括：28.所述时偏值29.其中，nfft为4096，若为两梳分，则l＝2×nm，若为四梳分，则l＝4×nm，angle为反正切函数；30.所述步骤s7进一步包括：31.所述补偿信道估计32.更进一步的，所述步骤s8进一步包括：33.所述协方差矩阵34.其中，τmax为信道传播最大时延量，δf为载波空间，k为整个带宽的载波索引值，k′为srs载波索引值；35.所述mmse矩阵其中，snr0为初始信噪比；36.在所述步骤s13中，所述新的mmse矩阵37.更进一步的，所述步骤s9进一步包括：38.所述权值w(k，l；k′，l′)＝θ(k)×φ-1(k′)；39.所述插值信道估计h"srs(k，l，r，p)＝wt(k，l；k′，l′)×h′srs(k′，l，r，p)；其中，t表示矩阵转置；40.所述步骤s10进一步包括：41.所述最终信道估计其中，taest,pre＝40ts。42.更进一步的，所述步骤s11进一步包括：43.所述频带上的信号功率pu＝mean(xe×xeh)；其中，所述xe＝∑pxsrs(k，l，p)×hsrs(k，l，r，p)，所述xeh为xe的转置共轭；44.所述频带上的噪声功率ni＝mean(δ×δh)；其中，δ＝ysrs(k，l，r)-∑pxsrs(k，l，p)×hsrs(k，l，r，p)，所述δh是δ的转置共轭。45.为解决上述技术问题，本发明还提供一种基于多用户srs信号的snr估计装置，包括第一信号单元，第二信号单元，第一运算单元，第二运算单元，时偏单元，信道估计单元和信噪比计算单元；46.所述第一信号单元从接收到的频域数据中提取出多个ue的srs测量信号ysrs(k，l，r)；其中，k为接收到的srs信号子载波索引，l为ofdm符号，r为接收天线；47.所述第二信号单元根据3gpp协议生成每个ue的本地发生序列xsrs(k，l，p)；其中，p为发送天线端口索引；48.所述第一运算单元用于依据所述每个ue的最终信道估计hsrs(k，l，r，p)和所述每个ue的本地发生序列xsrs(k，l，p)计算频带上的信号功率pu；49.所述信道估计单元用于计算所述每个ue的最终信道估计hsrs(k，l，r，p)；所述每个ue的最终信道估计hsrs(k，l，r，p)通过对每个ue的插值信道估计h"srs(k，l，r，p)进行解时偏补偿处理得到；所述每个ue的插值信道估计h"srs(k，l，r，p)依据权值w(k，l；k′，l′)对每个ue的补偿信道估计h′srs(k，l，r，p)进行re级插值运算得到；所述权值w(k，l；k′，l′)依据不同时刻不同子载波间协方差矩阵θ(k)和mmse矩阵φ(k′)计算得到；所述不同时刻不同子载波间协方差矩阵θ(k)和所述mmse矩阵φ(k′)依据mmse均衡算法对每个ue的补偿信道估计h′srs(k，l，r，p)进行插值滤波处理得到；所述每个ue的补偿信道估计h′srs(k，l，r，p)依据时偏值ta对每个ue的中间信道估计进行时偏补偿得到；所述每个ue的中间信道估计通过对每个ue的预补偿信道估计进行连续nm子载波平滑去干扰处理得到；所述每个ue的预补偿信道估计通过对所述每个ue的粗信道估计进行时偏预补偿得到；所述每个ue的粗信道估计依据所述多个ue的srs测量信号ysrs(k，l，r)和所述每个ue的本地发生序列xsrs(k，l，p)，基于最小二乘估计算法计算得到；其中，所述基于最小二乘估计算法计算得到；其中，所述基于最小二乘估计算法计算得到；其中，所述是srs的端口数；nu是用户数，或8；50.所述第二运算单元用于依据所述多个ue的srs测量信号ysrs(k，l，r)，所述最终信道估计hsrs(k，l，r，p)和所述本地发生序列xsrs(k，l，p)计算出频带上的噪声功率ni；51.所述时偏单元用于利用所述每个ue的中间信道估计进行时偏估计，得到所述时偏值ta；52.所述信噪比计算单元用于依据所述频带上的信号功率pu和所述频带上的噪声功率ni，基于信噪比计算公式确定中间信噪比snr′，将中间信噪比snr′输出到所述第一运算单元；再接收所述第一运算单元输出的新的频带上的信号功率pu，和所述第二运算单元输出的新的频带上的噪声功率ni，基于信噪比计算公式得到最终信噪比snrest；53.所述信道估计单元还用于根据所述中间信噪比snr′和所述不同时刻不同子载波间协方差矩阵θ(k)，得到新的mmse矩阵φ(k′)；再依据所述不同时刻不同子载波间协方差矩阵θ(k)和所述新的mmse矩阵φ(k′)计算得到新的权值w(k，l；k′，l′)，再依据新的权值w(k，l；k′，l′)对所述每个ue的补偿信道估计h′srs(k，l，r，p)进行re级插值运算，得到每个ue的插值信道估计h"srs(k，l，r，p)，再对所述每个ue的插值信道估计h"srs(k，l，r，p)进行解时偏补偿处理得到新的最终信道估计hsrs(k，l，r，p)；54.所述第一运算单元还用于依据新的最终信道估计hsrs(k，l，r，p)和所述本地发生序列xsrs(k，l，p)计算得到新的频带上的信号功率pu，再输出给所述信噪比计算单元；55.所述第二运算单元还用于依据新的最终信道估计hsrs(k，l，r，p)，所述多个ue的srs测量信号ysrs(k，l，r)和所述本地发生序列xsrs(k，l，p)，计算出新的频带上的噪声功率ni。56.更进一步的，所述每个ue的粗信道估计更进一步的，所述每个ue的粗信道估计57.所述每个ue的预补偿信道估计其中，taest,pre＝40ts；58.所述每个ue的中间信道估计59.所述协方差矩阵其中，τmax为信道传播最大时延量，δf为载波空间，k为整个带宽的载波索引值，k′为srs载波索引值；60.所述mmse矩阵其中，snr0为初始信噪比；61.所述新的mmse矩阵62.所述权值w(k，l；k′，l′)＝θ(k)×φ-1(k′)；63.所述插值信道估计h"srs(k，l，r，p)＝wt(k，l；k′，l′)×h′srs(k′，l，r，p)；其中，t表示矩阵转置；64.所述最终信道估计65.所述频带上的信号功率pu＝mean(xe×xeh)；其中，所述xe＝∑pxsrs(k，l，p)×hsrs(k，l，r，p)，所述xeh为xe的转置共轭；66.所述频带上的噪声功率ni＝mean(δ×δh)；其中，所述δ＝ysrs(k，l，r)-∑pxsrs(k，l，p)×hsrs(k，l，r，p)，所述δh是δ的转置共轭。67.更进一步的，所述时偏值68.其中，nfft为4096，若为两梳分，则l＝2×nm，若为四梳分，则l＝4×nm，angle为反正切函数；69.所述每个ue的补偿信道估计70.与现有技术相比，本发明具有如下有益效果：本发明采用在5gnr通信协议框架下周期性发送的srs信号计算信号功率与噪声，采用最小二乘法准则ls算法计算出信道估计响应。此外，本发明还对各自用户的信道估计响应进行时偏估计与补偿，使信道估计响应更加精确，从而解决现有技术中snr测量不准确的问题，进而改善5g小基站系统的通信质量。【附图说明】71.图1是本发明实施例的基于多用户srs信号的snr估计方法步骤图；72.图2是本发明实施例的基于多用户srs信号的snr估计装置的结构框图。【具体实施方式】73.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。74.需要说明的是，本技术的说明书和权利要求书中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便按本发明实施例以外的其他顺序实施。75.如图1所示，本发明实施例的基于多用户srs信号的snr估计方法，包括以下步骤：76.s1.从接收到的频域数据中提取出多个ue的srs测量信号ysrs(k，l，r)；其中，k为收到的srs信号的子载波索引，根据协议的规定，k为偶数，可以是0,2,4,8...等，也可以是0，4，8，12...等。1为ofdm(orthogonalfrequencydivisionmultiplexing，正交频分复用)符号，r为接收天线。77.s2.根据3gpp协议生成每个ue的本地发生序列xsrs(k，l，p)；其中，p为发送天线端口索引。78.3gppts38.211协议规定了产生天线端口pi的srs序列xsrs(k，l，p)的生成公式：[0079][0080]其中：[0081]个连续的ofdm符号；[0082][0082]是srs序列长度，i是天线端口的索引值。[0083]msrs,b表示srs在频域占用的rb个数，可参考3gppts38.211协议的table6.4.1.4.3-1表格设置。令b＝bsrs，bsrs∈{0，1，2，3}，csrs∈{0，1，...，63}是srs带宽配置索引。均由上层参数freqhopping设置，以此确定msrs,b的取值。ktc是传输梳的数量，取值2或4，包含在高层参数transmissioncomb中。δ＝log2(ktc)。[0084]天线端口pi的循环移位αi根据如下公式获得：[0085][0086][0087][0088]其中，包含在高层参数transmissioncomb中，协议规定，若ktc＝4，则若ktc＝2，则[0089]为低峰均比系列，由以下公式产生：[0090][0091]其中，ru，v(k)是基序列，是序列长度，是每个rb的载波个数，j为复数，α是循环移位，通过不同的α和δ，可以从单个基序列上产生多个序列。[0092]基序列ru，v(k)分成多个组，其中u∈{0，1，...，29}是组号，v是序列内的基序列号，当一个组只包含一个基序列(v＝0)时，每个基序列的长度时当一个组只包含一个基序列(v＝0)时，每个基序列的长度时其中1/2≤m/2δ≤5。此时用此配置，一个组只包含一个基序列。[0093]基序列ru，v(0)，...，ru，v(mzc-1)的定义依赖于长度mzc。当基序列长度大于等于36，即基序列ru，v(0)，...，ru，v(mzc-1)由以下公式定义：[0094]ru，v(k)＝xq(kmodnzc)[0095][0096]其中：π＝3.14，[0097]长度nzc是满足nzc＜mzc的最大质数。[0098]当基序列长度小于36时，分为以下两种情况：[0099]对于mzc＝30，[0100]对于mzc∈{6，12，18,24}，[0101]由3gppts38.211协议5.2.2节的4个表格定义，分别对应mzc等于6/12/18和24共4种情况，在此不再赘述。[0102]s3.依据多个ue的srs测量信号ysrs(k，l，r)和每个ue的本地发生序列xsrs(k，l，p)，基于最小二乘估计算法，计算得到每个ue的粗信道估计即[0103]s4.对每个ue的粗信道估计进行时偏预补偿，得到预补偿信道估计即其中，taest,pre＝40ts。[0104]s5.对每个ue的预补偿信道估计进行连续nm子载波平滑去干扰处理，得到每个ue的中间信道估计[0105]即，[0106]其中，其中，是srs的端口数，nu是用户数，或8。[0107]s6.利用每个ue的中间信道估计进行时偏估计，得到时偏值ta。[0108]具体的，时偏值[0109]其中，nfft为4096，若为两梳分，则l＝2×nm，若为四梳分，则l＝4×nm，angle为反正切函数。[0110]s7.依据时偏值ta对每个ue的中间信道估计进行时偏补偿，得到每个ue的补偿信道估计h′srs(k，l，r，p)。[0111]具体的，补偿信道估计[0112]s8.依据mmse均衡算法对每个ue的补偿信道估计h′srs(k，l，r，p)进行插值滤波处理，得到不同时刻不同子载波间协方差矩阵θ(k)及mmse矩阵φ(k′)。[0113]对于srs信道估计，由于只需配置单符号，因此只需考虑频域插值，时域插值可以忽略，因此，协方差矩阵mmse矩阵[0114]其中，τmax为信道传播最大时延量，δf为载波空间，可配置为30khz，k为整个带宽的载波索引值，k′为srs载波索引值；snr0为初始信噪比，可设置为30db。[0115]s9.依据不同时刻不同子载波间协方差矩阵θ(k)和mmse矩阵φ(k′)计算权值w(k，l；k′，l′)，再依据权值w(k，l；k′，l′)对每个ue的补偿信道估计h′srs(k，l，r，p)进行re级插值运算，得到每个ue的插值信道估计h"srs(k，l，r，p)。[0116]权值w(k，l；k′，l′)＝θ(k)×φ-1(k′)。[0117]插值信道估计h"srs(k，l，r，p)＝wt(k，l；k′，l′)×h′srs(k′，l，r，p)；其中，t表示矩阵转置。[0118]s10.对每个ue的插值信道估计h"srs(k，l，r，p)进行解时偏补偿处理，得到最终信道估计hsrs(k，l，r，p)。[0119]具体的，最终信道估计其中，taest,pre＝40ts。[0120]s11.依据每个ue的最终信道估计hsrs(k，l，r，p)和每个ue的本地发生序列xsrs(k，l，p)计算频带上的信号功率pu；依据srs测量信号ysrs(k，l，r)，最终信道估计hsrs(k，l，r，p)和本地发生序列xsrs(k，l，p)计算出频带上的噪声功率ni。[0121]具体的，频带上的信号功率pu＝mean(xe×xeh)；其中，xe＝∑pxsrs(k，l，p)×hsrs(k，l，r，p)，xeh为xe的转置共轭。[0122]xe计算出来为复数，xe×xeh表示xe乘以自身的转置共轭，可将其转换为实数，即信号功率。[0123]频带上的噪声功率ni＝mean(δ×δh)；其中，δ＝ysrs(k，l，r)-∑pxsrs(k，l，p)×hrs(k，l，r，p)，δh是δ的转置共轭。[0124]δ表示在srs的ofdm符号上，所有上的噪声值，此噪声值为复数，δ×δh表示δ乘以自身的转置共轭，可将其转换为实数，即噪声功率。[0125]mean是指计算平均值。[0126]s12.依据频带上的信号功率pu和频带上的噪声功率ni，基于信噪比计算公式确定中间信噪比snr′。[0127]s13.根据中间信噪比snr′和协方差矩阵θ(k)，得到新的mmse矩阵φ(k′)；然后，返回步骤s9进行再一次依序逐步计算，直到步骤s12依据信噪比计算公式得到最终信噪比snrest。[0128]具体的，本步骤中新的mmse矩阵[0129]然后，返回步骤9，计算权值w(k，l；k′，l′)，将步骤s11计算得出的新的mmse矩阵φ(k′)输入权值计算公式w(k，l；k′，l′)＝θ(k)×φ-1(k′)，得到新的权值w(k，l；k′，l′)，再将新的权值w(k，l；k′，l′)输入插值信道估计的计算公式h"srs(k，l，r，p)＝wt(k，l；k′，l′)×h′srs(k′，l，r，p)，得到新的插值信道估计h"srs(k，l，r，p)；再执行步骤s10，将新的插值信道估计h"srs(k，l，r，p)再进行解时偏补偿处理，得到新的最终信道估计hsrs(k，l，r，p)；再执行步骤s11，将新的最终信道估计hsrs(k，l，r，p)输入频带上的信号功率的计算公式pu＝mean(xe×xeh)，xe＝∑pxsrs(k，l，p)×hsrs(k，l，r，p)，得到新的频带上的信号功率pu；将新的最终信道估计hsrs(k，l，r，p)输入频带上的噪声功率的计算公式ni＝mean(δ×δh)；δ＝ysrs(k，l，r)-∑pxsrs(k，l，p)×hsrs(k，l，r，p)，得到新的频带上的噪声功率ni；再执行步骤s12，将新的频带上的信号功率pu和新的频带上的噪声功率ni输入信噪比计算公式得到最终信噪比snrest。[0130]如图2所示，本发明实施例基于多用户srs信号的snr估计装置，包括第一信号单元，第二信号单元，第一运算单元，第二运算单元，时偏单元，信道估计单元和信噪比计算单元。[0131]第一信号单元从接收到的频域数据中提取出多个ue的srs测量信号ysrs(k，l，r)；其中，k为接收到的srs信号的子载波索引，l为ofdm符号，r为接收天线。[0132]第二信号单元根据3gpp协议生成每个ue的本地发生序列xsrs(k，l，p)；其中，p为发送天线端口索引。[0133]第一运算单元用于依据每个ue的最终信道估计hsrs(k，l，r，p)和每个ue的本地发生序列xsrs(k，l，p)计算频带上的信号功率pu。[0134]信道估计单元用于计算每个ue的最终信道估计hsrs(k，l，r，p)；每个ue的最终信道估计hsrs(k，l，r，p)通过对每个ue的插值信道估计h"srs(k，l，r，p)进行解时偏补偿处理得到；每个ue的插值信道估计h"srs(k，l，r，p)依据权值w(k，l；k′，l′)对每个ue的补偿信道估计h′srs(k，l，r，p)进行re级插值运算得到；权值w(k，l；k′，l′)依据不同时刻不同子载波间协方差矩阵θ(k)和mmse矩阵φ(k′)计算得到；不同时刻不同子载波间协方差矩阵θ(k)和mmse矩阵φ(k′)依据mmse均衡算法对每个ue的补偿信道估计h′srs(k，l，r，p)进行插值滤波处理得到；每个ue的补偿信道估计h′srs(k，l，r，p)依据时偏值ta对每个ue的中间信道估计进行时偏补偿得到；每个ue的中间信道估计通过对每个ue的预补偿信道估计进行连续nm子载波平滑去干扰处理得到；每个ue的预补偿信道估计通过对每个ue的粗信道估计进行时偏预补偿得到；每个ue的粗信道估计依据多个ue的srs测量信号ysrs(k，l，r)和每个ue的本地发生序列xsrs(k，l，p)，基于最小二乘估计算法计算得到；其中，ue的本地发生序列xsrs(k，l，p)，基于最小二乘估计算法计算得到；其中，是srs的端口数；nu是用户数，或8。[0135]具体的，每个ue的粗信道估计[0136]每个ue的预补偿信道估计其中，taest,pre＝40ts。[0137]每个ue的中间信道估计[0138]协方差矩阵其中，τmax为信道传播最大时延量，δf为载波空间，可配置为30khz，k为整个带宽的载波索引值，k′为srs载波索引值。[0139]mmse矩阵其中，snr0为初始信噪比，可设置为30db。新的mmse矩阵[0140]权值w(k，l；k′，l′)＝θ(k)×φ-1(k′)；插值信道估计h"srs(k，l，r，p)＝wt(k，l；k′，l′)×h′srs(k′，l，r，p)；其中，t表示矩阵转置。[0141]所述插值信道估计h"srs(k，l，r，p)＝wt(k，l；k′，l′)×h′srs(k′，l，r，p)；其中，t表示矩阵转置。[0142]最终信道估计[0143]第二运算单元用于依据多个ue的srs测量信号ysrs(k，l，r)，最终信道估计hsrs(k，l，r，p)和本地发生序列xsrs(k，l，p)计算出频带上的噪声功率ni。[0144]时偏单元用于利用每个ue的中间信道估计进行时偏估计，得到时偏值ta。[0145]具体的，时偏值[0146]其中，nfft为4096，若为两梳分，则l＝2×nm，若为四梳分，则l＝4×nm，angle为反正切函数；[0147]每个ue的补偿信道估计[0148]信噪比计算单元用于依据频带上的信号功率pu和频带上的噪声功率ni，基于信噪比计算公式确定中间信噪比snr′，将中间信噪比snr′输出到第一运算单元；再接收第一运算单元输出的新的频带上的信号功率pu，和第二运算单元输出的新的频带上的噪声功率ni，基于信噪比计算公式得到最终信噪比snrest。[0149]信道估计单元还用于根据中间信噪比snr′和协方差矩阵θ(k)，得到新的mmse矩阵φ(k′)；再依据不同时刻不同子载波间协方差矩阵θ(k)和新的mmse矩阵φ(k′)计算得到新的权值w(k，l；k′，l′)，再依据新的权值w(k，l；k′，l′)对每个ue的补偿信道估计h′srs(k，l，r，p)进行re级插值运算，得到每个ue的插值信道估计h"srs(k，l，r，p)，再对每个ue的插值信道估计h"srs(k，l，r，p)进行解时偏补偿处理得到新的最终信道估计hsrs(k，l，r，p)。[0150]其中，新的mmse矩阵[0151]第一运算单元还用于依据新的最终信道估计hsrs(k，l，r，p)和本地发生序列xsrs(k，l，p)计算得到新的频带上的信号功率pu，再输出给信噪比计算单元。[0152]具体的，频带上的信号功率pu＝mean(xe×xeh)；其中，xe＝∑pxsrs(k，l，p)×hsrs(k，l，r，p)，xeh为xe的转置共轭。[0153]第二运算单元还用于依据新的最终信道估计hsrs(k，l，r，p)，多个ue的srs测量信号ysrs(k，l，r)和本地发生序列xsrs(k，l，p)，计算出新的频带上的噪声功率ni。[0154]具体的，频带上的噪声功率ni＝mean(δ×δh)；其中，δ＝ysrs(k，l，r)-∑pxsrs(k，l，p)×hsrs(k，l，r，p)，δh是δ的转置共轭。[0155]综上所述，本发明采用在5gnr通信协议框架下周期性发送的srs信号计算信号功率与噪声，采用最小二乘法准则ls算法计算出信道估计响应。此外，本发明还对各自用户的信道估计响应进行时偏估计与补偿，使信道估计响应更加精确，从而解决现有技术中snr测量不准确的问题，进而改善5g小基站系统的通信质量。[0156]以上实施例仅表达了本发明的优选实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，如对各个实施例中的不同特征进行组合等，这些都属于本发明的保护范围。当前第1页12当前第1页12
技术特征：
1.一种基于多用户srs信号的snr估计方法，其特征在于，包括以下步骤：s1.从接收到的频域数据中提取出多个ue的srs测量信号ysrs(k,l,r)；其中，k为接收到的srs信号的子载波索引，l为ofdm符号，r为接收天线；s2.根据3gpp协议生成每个ue的本地发生序列xsrs(k,l,p)；其中，p为发送天线端口索引；s3.依据所述多个ue的srs测量信号ysrs(k,l,r)和所述每个ue的本地发生序列xsrs(k,l,p)，基于最小二乘估计算法，计算得到每个ue的粗信道估计s4.对所述每个ue的粗信道估计进行时偏预补偿，得到预补偿信道估计s5.对所述每个ue的预补偿信道估计进行连续nm子载波平滑去干扰处理，得到每个ue的中间信道估计其中，所述是srs的端口数，nu是用户数，或8；s6.利用所述每个ue的中间信道估计进行时偏估计，得到时偏值ta；s7.依据所述时偏值ta对所述每个ue的中间信道估计进行时偏补偿，得到每个ue的补偿信道估计h
′
srs
(k,l,r,p)；s8.依据mmse均衡算法对所述每个ue的补偿信道估计h
′
srs
(k,l,r,p)进行插值滤波处理，得到不同时刻不同子载波间协方差矩阵θ(k)及mmse矩阵φ(k')；s9.依据所述不同时刻不同子载波间协方差矩阵θ(k)和所述mmse矩阵φ(k')计算权值w(k,l；k
′
,l
′
)，再依据所述权值w(k,l；k
′
,l
′
)对所述每个ue的补偿信道估计h
′
srs
(k,l,r,p)进行re级插值运算，得到每个ue的插值信道估计h"
srs
(k,l,r,p)；s10.对所述每个ue的插值信道估计h"
srs
(k，l，r，p)进行解时偏补偿处理，得到每个ue的最终信道估计h
srs
(k，l，r，p)；s11.依据所述每个ue的最终信道估计h
srs
(k，l，r，p)和所述每个ue的本地发生序列xsrs(k，l，p)计算频带上的信号功率pu；依据所述srs测量信号ysrs(k，l，r)，所述最终信道估计h
srs
(k，l，r，p)和所述本地发生序列xsrs(k，l，p)计算出频带上的噪声功率ni；s12.依据所述频带上的信号功率pu和所述频带上的噪声功率ni，基于信噪比计算公式确定中间信噪比snr
′
；s13.根据所述中间信噪比snr
′
和所述协方差矩阵θ(k)，得到新的mmse矩阵φ(k
′
)；然后，返回所述步骤s9进行再一次依序逐步计算，直到所述步骤s12依据所述信噪比计算公式得到最终信噪比snrest。2.如权利要求1所述的基于多用户srs信号的snr估计方法，其特征在于，所述步骤s3进一步包括：所述粗信道估计
所述步骤s4进一步包括：所述预补偿信道估计其中，ta
est,pre
＝40ts。3.如权利要求1所述的基于多用户srs信号的snr估计方法，其特征在于，所述步骤s5进一步包括：所述中间信道估计4.如权利要求1所述的基于多用户srs信号的snr估计方法，其特征在于，所述步骤s6进一步包括：所述时偏值其中，nfft为4096，若为两梳分，则l＝2
×
n
m
，若为四梳分，则l＝4
×
n
m
，angle为反正切函数；所述步骤s7进一步包括：所述补偿信道估计5.如权利要求1所述的基于多用户srs信号的snr估计方法，其特征在于，所述步骤s8进一步包括：所述协方差矩阵其中，τ
max
为信道传播最大时延量，δf为载波空间，k为整个带宽的载波索引值，k
′
为srs载波索引值；所述mmse矩阵其中，snr0为初始信噪比；在所述步骤s13中，所述新的mmse矩阵6.如权利要求1所述的基于多用户srs信号的snr估计方法，其特征在于，所述步骤s9进一步包括：所述权值w(k，l；k
′
，l
′
)＝θ(k)
×
φ-1
(k
′
)；所述插值信道估计h"
srs
(k，l，r，p)＝w
t
(k，l；k
′
，l
′
)
×
h
′
srs
(k
′
，l，r，p)；其中，t表示矩阵转置；所述步骤s10进一步包括：所述最终信道估计其中，ta
est,pre
＝40ts。7.如权利要求1所述的基于多用户srs信号的snr估计方法，其特征在于，所述步骤s11进一步包括：所述频带上的信号功率pu＝mean(x
e
×
x
eh
)；其中，所述x
e
＝∑
p
x
srs
(k，l，p)
×
h
srs
(k，l，r，p)，所述x
eh
为x
e
的转置共轭；所述频带上的噪声功率ni＝mean(δ
×
δ
h
)；其中，δ＝y
srs
(k，l，r)-∑
p
x
srs
(k，l，p)
×
h
srs
(k，l，r，p)，所述δ
h
是δ的转置共轭。8.一种基于多用户srs信号的snr估计装置，其特征在于，包括第一信号单元，第二信号单元，第一运算单元，第二运算单元，时偏单元，信道估计单元和信噪比计算单元；所述第一信号单元从接收到的频域数据中提取出多个ue的srs测量信号ysrs(k，l，r)；其中，k为接收到的srs信号子载波索引，l为ofdm符号，r为接收天线；所述第二信号单元根据3gpp协议生成每个ue的本地发生序列xsrs(k，l，p)；其中，p为发送天线端口索引；所述第一运算单元用于依据每个ue的最终信道估计h
srs
(k，l，r，p)和所述每个ue的本地发生序列xsrs(k，l，p)计算频带上的信号功率pu；所述信道估计单元用于计算所述每个ue的最终信道估计h
srs
(k，l，r，p)；所述最终信道估计h
srs
(k，l，r，p)通过对每个ue的插值信道估计h"
srs
(k，l，r，p)进行解时偏补偿处理得到；所述每个ue的插值信道估计h"
srs
(k，l，r，p)依据权值w(k，l；k
′
，l
′
)对每个ue的补偿信道估计h
′
srs
(k，l，r，p)进行re级插值运算得到；所述权值w(k，l；k
′
，l
′
)依据不同时刻不同子载波间协方差矩阵θ(k)和mmse矩阵φ(k
′
)计算得到；所述不同时刻不同子载波间协方差矩阵θ(k)和所述mmse矩阵φ(k
′
)依据mmse均衡算法对每个ue的补偿信道估计h
′
srs
(k，l，r，p)进行插值滤波处理得到；所述每个ue的补偿信道估计h
′
srs
(k，l，r，p)依据时偏值ta对每个ue的中间信道估计进行时偏补偿得到；所述每个ue的中间信道估计通过对每个ue的预补偿信道估计进行连续nm子载波平滑去干扰处理得到；所述每个ue的预补偿信道估计通过对所述每个ue的粗信道估计进行时偏预补偿得到；所述每个ue的粗信道估计依据所述多个ue的srs测量信号ysrs(k，l，r)和所述每个ue的本地发生序列xsrs(k，l，p)，基于最小二乘估计算法计算得到；其中，所述是srs的端口数；nu是用户数，或8；所述第二运算单元用于依据所述多个ue的srs测量信号ysrs(k，l，r)，所述最终信道估计h
srs
(k，l，r，p)和所述本地发生序列xsrs(k，l，p)计算出频带上的噪声功率ni；所述时偏单元用于利用所述每个ue的中间信道估计进行时偏估计，得到所述时偏值ta；所述信噪比计算单元用于依据所述频带上的信号功率pu和所述频带上的噪声功率ni，基于信噪比计算公式确定中间信噪比snr
′
，将中间信噪比snr
′
输出到所述第一运算单元；再接收所述第一运算单元输出的新的频带上的信号功率pu，和所述第二运算单元输出的新的频带上的噪声功率ni，基于信噪比计算公式得到最终信噪比snrest；所述信道估计单元还用于根据所述中间信噪比snr
′
和所述不同时刻不同子载波间协方差矩阵θ(k)，得到新的mmse矩阵φ(k
′
)；再依据所述不同时刻不同子载波间协方差矩阵θ(k)和所述新的mmse矩阵φ(k
′
)计算得到新的权值w(k，l；k
′
，l
′
)，再依据新的权值w(k，l；
k
′
，l
′
)对所述每个ue的补偿信道估计h
′
srs
(k，l，r，p)进行re级插值运算，得到每个ue的插值信道估计h"
srs
(k，l，r，p)，再对所述每个ue的插值信道估计h"
srs
(k，l，r，p)进行解时偏补偿处理得到新的最终信道估计h
srs
(k，l，r，p)；所述第一运算单元还用于依据新的最终信道估计h
srs
(k，l，r，p)和所述本地发生序列xsrs(k，l，p)计算得到新的频带上的信号功率pu，再输出给所述信噪比计算单元；所述第二运算单元还用于依据新的最终信道估计h
srs
(k，l，r，p)，所述多个ue的srs测量信号ysrs(k，l，r)和所述本地发生序列xsrs(k，l，p)，计算出新的频带上的噪声功率ni。9.如权利要求8所述的基于多用户srs信号的snr估计装置，其特征在于，所述每个ue的粗信道估计所述每个ue的预补偿信道估计其中，ta
est,pre
＝40ts；所述每个ue的中间信道估计所述协方差矩阵其中，τ
max
为信道传播最大时延量，δf为载波空间，k为整个带宽的载波索引值，k
′
为srs载波索引值；所述mmse矩阵其中，snr0为初始信噪比；所述新的mmse矩阵所述权值w(k，l；k
′
，l
′
)＝θ(k)
×
φ-1
(k
′
)；所述插值信道估计h"
srs
(k，l，r，p)＝w
t
(k，l；k
′
，l
′
)
×
h
′
srs
(k
′
，l，r，p)；其中，t表示矩阵转置；所述最终信道估计所述频带上的信号功率pu＝mean(x
e
×
x
eh
)；其中，所述x
e
＝∑
p
x
srs
(k，l，p)
×
h
srs
(k，l，r，p)，所述x
eh
为x
e
的转置共轭；所述频带上的噪声功率ni＝mean(δ
×
δ
h
)；其中，所述δ＝y
srs
(k，l，r)-∑
p
x
srs
(k，l，p)
×
h
srs
(k，l，r，p)，所述δ
h
是δ的转置共轭。10.如权利要求8所述的基于多用户srs信号的snr估计装置，其特征在于，所述时偏值其中，n
fft
为4096，若为两梳分，则l＝2
×
n
m
，若为四梳分，则l＝4
×
n
m
，angle为反正切函数；所述每个ue的补偿信道估计

技术总结
本发明涉及基于多用户SRS信号的SNR估计方法和装置，采用在5G NR通信协议框架下周期性发送的SRS信号计算信号功率与噪声，采用最小二乘法准则LS算法计算出多个用户的信道估计响应。此外，本发明还对各自用户的信道估计响应进行时偏估计与补偿，使信道估计响应更加精确，从而解决现有技术中SNR测量不准确的问题，进而改善5G小基站系统的通信质量。进而改善5G小基站系统的通信质量。进而改善5G小基站系统的通信质量。


技术研发人员：赵强 许秋平 史涛
受保护的技术使用者：深圳国人无线通信有限公司
技术研发日：2023.04.07
技术公布日：2023/7/21