专利名称：一种用于自适应调制编码的信噪比估算方法
技术领域：
本发明涉及一种精确估算接收信噪比的方法，尤其涉及一种用于广义频分复用系统和单载波系统的自适应调制编码的接收信噪比估算方法，同时，本发明还可直接用于基于循环前缀的码分多址(CP-CDMA)系统。
背景技术：
用于OFDM(正交频分复用)和CP-SC(单载波)系统的自适应调制编码(AMC)在过去的很多文献中都被研究和提出过。在基于帧的自适应调制编码中，一帧数据的调制编码方式(MCS)是不变的，但根据信道的状态信息的不同，帧与帧之间的MCS是改变的。对于OFDM系统，误码率性能主要取决于信噪比最小的那个子载波。而对于CP-SC系统，每个比特的能量都是均匀分布在整个频谱上。因此，由衰落性信道造成的严重的频率衰落对CP-SC系统的影响并不大。这也就是为什么采用了基于帧的自适应调制编码的CP-SC系统的吞吐量可以远远超过同样采用了基于帧AMC的OFDM系统。如果对于OFDM系统的每个子载波采用不同的MCS，系统的吞吐量可以大大提高，但这需要增加大量的信令冗余，非常不实际。
广义频分复用系统(GOFDM)是作为OFDM(正交频分复用)和CP-SC(单载波)系统的折中方案而提出的。在GOFDM中，多个小尺寸的OFDM符号在时域进行级联组成一个FFT长度的GOFDM帧。同样的，循环前缀加在每一个GOFDM帧的前端。在接收端，采用频域均衡器对数据进行均衡。GOFDM系统可以被看成是传统的CP-SC系统当OFDM符号的子载波数为1，也可以被看成OFDM系统当OFDM符号的子载波数等于FFT的长度。
广义频分复用系统(GOFDM)是一个很好的候选方案，基于对复杂度，频域分集，峰均功率比，尤其是基于帧的AMC。
现有技术中，对于广义频分复用系统因为该系统没有被实际应用，因此还没有人提出基于GOFDM的信噪比计算方法。
发明内容
本发明所要解决的技术问题是提供一种用于自适应调制编码的信噪比估算方法，可用于在广义频分复用系统中精确计算接收信噪比。
为了解决上述技术问题上，本发明采用了如下的技术方案本发明的用于自适应调制编码的信噪比估算方法，其特征在于，将接收数据表示为含有UH的函数，其中，U为酉矩阵， 且UHU=UUH=INs×Ns,]]>WNm=1Nm1WNm1-1...WNm1-Nm1WNm2-1...WNm2-Nm............1WNmNm-1...WNmNm·NmNm×Nm,]]>且WNmHWNm=WNmWNmH=INm×Nm;]]>其中，WNm是FFT酉矩阵，也可以是正交Walsh-Hadamad变换矩阵，H表示矩阵的共扼转置。INm×Nm为对角元素为1的对角阵。
将上述含有UH的接收数据用于计算信噪比，获得的信噪比为一个对应于每个星座符号的接收信噪比的向量。
图1是广义频分复用系统的发送装置的结构示意图。
图2是广义频分复用系统的接收装置的结构示意图。
图3是多种调制编码方式对应的系统吞吐量的仿真性能示意图。
图4是采用基于帧的AMC的CP-SC，GOFDM和OFDM的系统吞吐量性能示意图。
具体实施例方式
图1给出了GOFDM系统的发送接收结构框图。在发送端，长为Ns的数据序列被调制到星座点上，进行串并转换将总长为Ns的发送序列d分成K组长度为Nm的短序列dk，Ns＝Nm×k，0≤k≤K-1。将这些短序列分别做Nm点的IFFT，则得到sk，可以表示为sk=WNmHdk---(1)]]>
其中， WNm是FFT酉矩阵，WNm=1Nm1WNm1-1...WNm1-Nm1WNm2-1...WNm2-Nm............1WNmNm-1...WNmNm·NmNm×Nm,]]>且WNmHWNm=WNmWNmH=INm×Nm.]]>H表示矩阵的共扼转置。INm×Nm为对角元素为1的对角阵。然后将sk级联，再并串转换成序列s，s可以表示为s＝UHd(2)其中， 且UHU=UUH=INs×Ns.]]>最后序列s加上循环前缀后进行发送。
如图2所示经过信道后在接收端，首先将接收到的数据去除循环前缀后得到r，r可以表示为r＝Cs+n(3)＝CUHd+n其中，信道矩阵C是循环矩阵表示为 andh＝[h0h1…hl-1]是信道的时域冲击响应，l为信道阶数，n为高斯白噪声。另外信道矩阵C还可以表示为C=WNsHHWNs---(4)]]>其中 为Ns×Ns的对角阵，对角元素为信道的频率响应，其余为零。WNs为Ns×Ns的FFT酉矩阵。将式(4)代入(3)得到r=CUHd+n]]>=WNsHHWNsUHd+n---(5)]]>去除CP后在时域得到的接收信号通过Ns点FFT可以变换到频域，R=WNsr]]>=WNs(WNsHHWNsUHd+n)---(6)]]>=HWNsUHd+WNsn]]>1.1迫零(ZF)均衡采用ZF均衡，即将经过Ns点的FFT变换后的数据序列R乘以频率信道矩阵(HHH)-1HHR^=(HHH)-1HHR]]>=(HHH)-1HH(HWNsUHd+WNsn)---(7)]]>=WNsUHd+H-1WNsn]]>再进行Ns点的IFFT变换得到 ，表示为r^=WNsHR^]]>=WNsH(WNsUHd+H-1WNsn)---(8)]]>=UHd+WNsHH-1WNsn]]>最后，将 进行分组，做Nm点的FFT变换，得到 ，可以表示为d^=Ur^]]>=U(UHd+WNsHH-1WNsn)---(9)]]>=d+UWNsHH-1WNsn]]>因为ZF均衡是没有ISI的，因此接收信噪比可以表示为
SNRZF=σs2·I(UWNsHH-1WNsn)(UWNsHH-1WNsn)H]]>=σs2·IUWNsHH-1WNsn·nHWNsHH-1,HWNsUH---(10)]]>=σs2·IUWNsHσn2(HHH)-1WNsUH]]>其中σn2为噪声方差，σs2为信号能量，假定为1。INs×Ns为对角元素为1的对角阵。由于W和U都是酉矩阵，因此SNRZF是一个对应于Ns个星座符号的接收信噪比的向量。
1.2最小均方误差(MMSE)均衡MMSE均衡是将序列R乘以(HHH+Iσn2)-1HH，得到R^=(HHH+Iσn2)-1HHR]]>=(HHH+Iσn2)-1HH(HWNsUHd+WNsn)---(11)]]>=(HHH+Iσn2)-1HHHWNsUHd+(HHH+Iσn2)-1HHWNsn]]>同样的，再进行Ns点的IFFT变换得到 ，表示为r^=WSsHR^]]>=WSsH((HHH+Iσn2)-1HHHWNsUHd+(HHH+Iσn2)-1HHWNsn)---(12)]]>=WNsH(HHH+Iσn2)-1HHHWNsUHd+WNsH(HHH+Iσn2)-1HHWNsn]]>最后，将 进行分组，做Nm点的FFT变换，得到 ，可以表示为d^=Ur^]]>=U(WNsH(HHH+Iσn2)-1HHHWNsUHd+WNsH(HHH+Iσn2)-1HHWNxn)]]>=UWNsH(HHH+Iσn2)-1HHHWNsUHd+UWNsH(HHH+Iσn2)-1HHWNsn---(13)]]>
其中，第一项是期望的信号，中间一项是残余的符号间干扰ISI，最后一项是加性高斯噪声。残余的符号间干扰和加性噪声之和记为JJ=σs2(UWNsH(HHH+Iσn2)-1HHHWNsUH-INs×Ns)(UWNsH(HHH+Iσn2)-1HHHWNsUH-INs×Ns)H]]>+(UWNsH(HHH+Iσn2)-1HHWNsn)(UWNsH(HHH+Iσn2)-1HHWNsn)H]]>=I-UWNsH(HHH+Iσn2)-1HHHWNsUH---(14)]]>最后接收信噪比可以表示为SNRMMSE=σs2·I-JJ]]>=UWNsH(HHH+Iσn2)-1HHHWNsUHI-UWNsH(HHH+Iσn2)-1HHHWNsUH---(15)]]>由于W和U都是酉矩阵，因此SNRMMSE是一个对应于Ns个星座符号的接收信噪比的向量。这种信噪比的分析方法也可以直接用于其他的系统，例如CP-CDMA。对CP-CDMA而言，唯一的区别就是将U矩阵替换为正交Walsh-Hadamad变换矩阵。
2.用于帧AMC的接收信噪比基于整个帧的AMC中，一帧的数据采用同一种调制编码方式，调制编码的选择是依据所有可选的调制编码方式下，系统在高斯白噪声环境中得到的误码率或者误帧率性能，设定信噪比门限，利用一帧数据的接收信噪比进行调制编码的选择。接收信噪比计算方法如下采用ZF均衡时，假设经过星座点调制的符号的能量σs2=1,]]>根据式(10)一帧数据的接收信噪比SNRZF_all表示为SNRZF_all=trace(σs2·IUWNsHσn2(HHH)-1WNsUH)]]>=NsΣq=0Ns-1|Hqq|2|Hqq|2+σn2---(16)]]>采用MMSE均衡时，根据式(15)，一帧数据的接收信噪比SNRMMSE_all表示为
SNRMMSE_all=trace(UWNsH(HHH+Iσn2)-1HHHWNsUHI-UWNsH(HHH+Iσn2)-1HHHWNsUH)]]>=Σq=0Ns-1|Hqq|2|Hqq|2+σn2Ns-Σq=0Ns-1|Hqq|2|Hqq|2+σn2---(17)]]>在一个GOFDM系统中，用于复用的OFDM符号长度Nm＝16，8个OFDM符号进行复用，k＝8，组成长度为Ns＝128的GOFDM帧，采用MMSE均衡。每一帧的MCS方式根据门限值由式(17)计算出来的SNRMMSE决定。
仿真结果仿真参数如表1所示
调制编码方式如表2所示
调制编码方式对应的信噪比门限如表3所示
信噪比SNR的门限设定式基于在高斯信道下，每一种调制编码方式的系统吞吐量的仿真性能，如图3所示。相邻两条曲线的交点就是门限值。采用基于帧的AMC的CP-SC，GOFDM和OFDM的系统吞吐量性能见图4，MCS的选择是根据式(17)计算出来的接收信噪比。可以明显的看到GOFDM系统的吞吐量介于CP-SC和OFDM之间。
权利要求
1.一种用于自适应调制编码的信噪比估算方法，其特征在于，将接收数据表示为含有UH的函数，U为酉矩阵， 且UHU=UUH=INs×Ns,]]>WNm=1Nm1WNm1.1...WNm1.Nm1WNm2.1...WNm2.Nm............1WNmNm.1...WNmNm.NmNm×Nm]]>，且WNmHWNm=WNmWNmH=INm×Nm]]>；其中，WNm是FFT酉矩阵，也可以是正交Walsh-Hadamad变换矩阵，H表示矩阵的共扼转置。INm×Nm为对角元素为1的对角阵。将上述含有UH的接收数据用于计算信噪比，获得的信噪比为一个对应于每个星座符号的接收信噪比的向量。
2.根据权利要求1所述的用于自适应调制编码的信噪比估算方法，其特征在于，包括步骤在发送端，长为Ns的数据序列被调制到星座点上，进行串并转换将总长为Ns的发送序列d分成K组长度为Nm的短序列dk，NS＝Nm×k，O≤k≤K-1，将所述短序列分别做Nm点的IFFT，则得到sk，可以表示为Sk=WNmHdk.]]>
3.根据权利要求2所述的用于自适应调制编码的信噪比估算方法，其特征在于，包括步骤将序列sk级联，再并串转换成序列s，s表示为s=UHd。
4.根据权利要求3所述的用于自适应调制编码的信噪比估算方法，其特征在于，包括步骤在序列s加上循环前缀后进行发送。
5.根据权利要求4所述的用于自适应调制编码的信噪比估算方法，其特征在于，包括步骤在接收端，首先将接收到的数据去除循环前缀后得到序列r，r表示为r=WNsHHWNsUHd+n,]]>其中， 为Ns×Ns的对角阵，对角元素为信道的频率响应，其余为零，n为高斯白噪声。
6.根据权利要求5所述的用于白适应调制编码的信噪比估算方法，其特征在于，包括步骤将上述去除循环前缀后的序列通过Ns点FFT可以变换到频域，获得数据序列RR=HWNsUHd+WNsn.]]>
7.根据权利要求6所述的用于自适应调制编码的信噪比估算方法，其特征在于，包括步骤对上述数据序列R进行zF均衡，即将数据序列R乘以频率信道矩阵(HHH)-1HH，获得R^=WNsUHd+H-1WNsn;]]>再进行NS点的IFFT变换得到r^=UHd+WNsHH-1WNsn;]]>将 进行分组，做Nm点的FFT变换，得到d^=d+UWNsHH-1WNsn;]]>最后得到接收信噪比为SNRZF=σs2.IUWNsHσn2(HHH)-1WNsUH]]>；其中，σn2为噪声方差，σs2为信号能量，假定为1，INs×Ns为对角元素为1的对角阵。
8.根据权利要求6所述的用于自适应调制编码的信噪比估算方法，其特征在于，包括步骤对上述数据序列R进行MMSE均衡，即将序列R乘以(HHH+Iσn2)-1HH，得到R^=(HHH+Iσn2)-1HHHWNsUHd+(HHH+Iσn2)-1HHWNsn;]]>再进行Ns点的IFFT变换得到 ，表示为r^=WNsH(HHH+Iσn2)-1HHHWNsUHd+WNsH(HHH+Iσn2)-1HHWNsn;]]>将 进行分组，做Nm点的FFT变换，得到 ，可以表示为 其中，第一项是期望的信号，中间一项是残余的符号间干扰ISI，最后一项是加性高斯噪声；残余的符号间干扰和加性噪声之和记为J:J=I-UWNsH(HHH+Iσn2)-1HHHWNsUH;]]>最后，接收信噪比可以表示为SNRMMSE=UWNsH(HHH+Iσn2)-1HHHWNsUHI-UWNsH(HHH+Iσn2)-1HHHWNsUH.]]>
全文摘要
本发明提供一种用于自适应调制编码的信噪比估算方法，将接收数据表示为含有U
文档编号H04L25/02GK101026433SQ20061002412
公开日2007年8月29日 申请日期2006年2月24日 优先权日2006年2月24日
发明者唐琳, 杨秀梅, 李明齐, 卜智勇, 张小东, 王海峰 申请人:上海无线通信研究中心