基于AmBC的双选信道LMMSE估计及ML检测方法

基于ambc的双选信道lmmse估计及ml检测方法
技术领域
1.本发明涉及基于ambc的双选信道lmmse估计及ml检测方法，属于无线通信中的信道估计与检测技术领域。


背景技术：

2.环境反向散射通信(ambient backscatter communication，ambc)是一种利用无源标签吸收环境中的射频源能量，并以环境中的信号为载波调制自身信息发射到接收端的技术。由于ambc不占用原通信系统以外的频谱资源，硬件设备结构简单，不需要额外的能量供应，灵活度较高，因而被广泛应用于射频识别、物联网通信等场景。
3.正交频分复用(orthogonal frequency division multiplexing,ofdm)是现代通信系统中,被广泛使用的一种技术。其原理是将输入信号由串行转为n组并行信号并对其中第k组信号采用频率为的子载波进行调制，其中对于任意m≠n,都有和之积在一个周期内积分恒为0，故任意两个不同频率的子载波互相正交。为避免符号间串扰，通常将每个ofdm信号的末端复制到其前端构成循环前缀，避免了子载波之间的相互干扰。最后将混频器输出的各路子信号相加即可得到ofdm信号。在接收端将ofdm信号由串转并后分别与各路子载波相乘，利用子载波之间的互相正交的关系即可积分得到符号映射后的估计值，从而恢复出发送信号。实际上，ofdm信号使用子载波调制和解调的过程分别相当于一组离散傅立叶逆变换(inverse discrete fourier transform，idft)和离散傅立叶变换(discrete fourier transform，dft)。因而，ofdm技术能被广泛应用，得益于快速傅里叶逆变换(inverse fast fourier transform，ifft)和快速傅里叶变换(fast fourier transform，fft)算法的提出带来的巨大简化，凸显出其可以对抗频率选择性衰落或窄带干扰和信道利用率高的优点。
4.目前，环境反向散射通信下的信道估计通常是假设平坦衰落信道。为了更好的模拟实际情况，本发明是基于较复杂的时间-频率双选信道模型。在ambc技术下，读写器收到的信号是直接由发射机经信道到达到读写器(称为直接链路)和从发射机先传输到反射设备，再从反射设备到达读写器(称为反射链路)的叠加，并且考虑到时间-频率双选信道模型对通信性能的负面影响，进行信号检测和信道估计就面临着更大的难度和挑战。
5.对于双选信道而言，信道的未知参数个数远大于大于导频符号个数，导致其无法被完全估计。为了兼顾运算复杂度和对ici的消除效果，本发明采用被广泛运用于ofdm系统时变信道建模的基扩展模型(basis expansion model,bem)，将信道冲激响应变换到一个由基向量张成的低维空间中从而进行信道估计。根据信号模型设计检测方法，再恢复出射频源符号和反射设备符号。


技术实现要素：

6.本发明目的在于针对高速移动场景下的信道复杂，导致多普勒频移和多径干扰，
进而使得环境反向散射系统接收信号误码率高的技术缺陷，提出了基于ambc的双选信道lmmse估计及ml检测方法。
7.本发明的核心思想为：1)将环境射频源到标签的信道以及环境射频源到接收端的信道建模为时间-频率双选信道，标签到接收端的信道建模为平坦衰落信道；2)针对环境信号和标签信号设计了块状导频簇；3)采用bem模型拟合信道状态，将估计信道单位冲激响应(channel impulse response，cir)转换为估计bem基的系数向量；4)用线性最小均方误差(linear mean square error，lmmse)估计分别对直接链路和反射链路信道进行估计；5)在信道估计后采用最小均方误差(minimum mean square error，mmse)检测算法检测环境信号；6)采用最大似然检测算法检测标签信号。
8.为了达到上述目的，本发明采取如下技术方案：
9.所述双选信道lmmse估计及ml检测方法，依托的系统包括环境射频源、标签和接收端；从环境射频源到接收端的直接链路信道和从环境射频源到标签的反射链路信道为双选信道；
10.其中，环境射频源简称射频源，射频源发射的ofdm信号简称为环境信号；标签发射ask信号，该ask信号简称为标签信号；双选信道具有多径效应和快时变特性；将接收端与标签间的信道建模为平坦衰落信道；原因在于接收端与标签间的距离比较短、距离恒定致使增益恒定。
11.所述双选lmmse信道估计及检测方法，包括如下步骤：
12.s1、基于块状导频簇结构分别构建环境信号和标签信号；
13.s2、环境信号经直接链路信道到达接收端得到直接链路信号；环境信号经过反射链路信道到达标签，标签吸收环境信号的能量，根据自身携带的比特信息选择对射频源信号进行调制，得到调制后符号，再将调制后符号经过平坦衰落信道到达接收端得到反射链路信号；接收端同时收到直接链路信号、反射链路信号和高斯白噪声得到总的接收信号；
14.s3、采用bem模型对直接链路信道和反射链路信道进行建模，然后分别构建直接链路观测方程和反射链路观测方程；
15.s3中，bem模型为时域信道矩阵，该矩阵为bem基矩阵与bem基的系数向量相乘的形式；
16.s4、分别对s3直接链路观测方程和反射链路观测方程中的直接链路bem基的系数向量及反射链路bem基的系数向量进行估计，然后再把bem基的系数向量恢复为直接链路时域信道矩阵和反射链路时域信道矩阵，具体包括如下子步骤：
17.s41、根据s3计算得到的直接链路观测方程进行直接链路信道估计，得到直接链路时域信道矩阵；
18.s42、根据s3计算得到的反射链路观测方程进行反射链路信道估计，得到反射链路时域信道矩阵；
19.s5、得到s4的直接链路时域信道矩阵后，对射频源信号进行mmse检测，得到ofdm符号估计值；
20.s6、采用最大似然检测法对标签信号进行检测，得到标签符号的检测结果，具体如下：
21.s61、分别对s4的直接链路时域信道矩阵和反射链路时域信道矩阵左乘一个傅里
叶矩阵后，再右乘一个“傅里叶矩阵的共轭转置”，得到直接链路频域信道矩阵和反射链路频域信道矩阵；
22.s62、将反射链路频域信道矩阵与反射系数相乘后再与直接链路频域信道矩阵相加得到复合链路频域信道矩阵；
23.s63、分别计算第i个数据符号时刻的直接链路频域信道矩阵和复合链路频域信道矩阵与s5得到的ofdm符号估计值相乘得到y0和y1；
24.s64、对总的接收信号做傅里叶变换得到频域接收信号，其中第i个数据符号时刻的频域接收信号为yi；
25.s65、分别计算yi与y0和y1的欧几里得距离d0和d1，若d0小于等于d1，则输出当前标签符号的检测结果为“0”；若d0大于d1，则输出当前标签符号的检测结果为“1”。
26.所述s1中块状导频簇的一组为两个符号，前一个符号持续时间内标签信号为导频“0”，后一个符号持续时间内标签信号为导频“1”；构建环境信号和标签信号，具体为：每间隔数个数据符号后插入一组块状导频簇；
27.s3具体包括如下子步骤：
28.s31、将bem基向量分别与维度为直接链路信道多径数与反射链路信道多径数的单位阵进行克罗内克求积，得到直接链路的bem基矩阵与反射链路的bem基矩阵；
29.s32、将第i个符号周期第l条路径上的环境信号用对角化表示为s
li
；直接链路信道多径数为ld，l的取值范围为0到l
d-1，反射链路信道多径数为lr，l的取值范围为0到l
r-1；
30.s33、将s32中直接链路所有路径上的环境信号用对角化表示s
li
合并为一个矩阵s
di
；
31.s34、对s33中的s
di
做离散傅里叶变换后与s31直接链路的bem基矩阵相乘得到直接链路观测矩阵a
di
，观测方程为y
di
＝a
di
×cdi
；
32.其中，y
di
为直接链路的频域接收信号，c
di
为直接链路的bem基的系数向量；
33.s35、将s32中反射链路所有路径上的环境信号用对角化表示s
li
合并为一个矩阵s
ri
；
34.s36、对s35中的s
ri
做离散傅里叶变换后与s31中的反射链路的bem基矩阵相乘得到反射链路观测矩阵a
ri
，观测方程为y
ri
＝a
ri
×cri
；
35.其中，y
ri
为反射链路的频域接收信号，c
ri
为反射链路的bem基的系数向量。
36.s41、具体为：
37.s411、在标签发送导频符号“0”的时刻用lmmse估计出直接链路的bem基的系数向量；
38.s412、将s3计算得到的bem基矩阵与s411计算得到的直接链路的bem基的系数向量相乘得到“0”时刻的直接链路时域信道矩阵；
39.s413、在标签发送导频符号“1”以及数据符号的时刻用信道插值计算出“1”时刻及数据符号时刻的直接链路时域信道矩阵。
40.s42、具体为：
41.s421、在标签发送导频符号“1”的时刻用估计出的直接链路时域信道矩阵与ofdm导频符号相乘恢复出“1”时刻的直接链路接收信号；
42.s422、在标签发送导频符号“1”的时刻用总的接收信号减去s421恢复出的直接链
路接收信号从而获得剩余信号；
43.s423、在标签发送导频符号“1”的时刻将s422计算得到的剩余信号做离散傅里叶变换得到反射链路频域接收信号y
ri
；
44.s424、在标签发送导频符号“1”的时刻用lmmse估计出反射链路的bem基的系数向量；
45.s425、将s3计算得到的bem基矩阵b与s424计算得到的反射链路的bem基的系数向量相乘得到1时刻的反射链路时域信道矩阵；在标签发送导频符号“0”以及数据符号的时刻用信道插值估计出“0”时刻及数据符号时刻的反射链路时域信道矩阵。
46.有益效果
47.所述基于ambc的双选信道lmmse估计及ml检测方法，与现有技术相比，具有如下有益效果：
48.1、所述方法通过设计基于块状导频簇的通信协议，如s1所述从而实现对信道的时变特性进行跟踪，使信道估计及信号检测算法的准确度显著提升；
49.2、所述方法采用bem模型对信道单位冲激响应进行建模，如s3所述在模型不丢失双选信道性质的情况下减少了信道待估计参数的数量，显著降低了信道估计方法的复杂度；
50.3、所述方法通过设计针对环境反向通信系统的双选信道估计，如s4所示实现了对信道时变特性的有效跟踪，使得反射链路和直射链路的信道估计性能显著提升；
51.4、所述方法通过mmse检测算法对时变信道下的射频源ofdm符号进行检测，对带状信道频域矩阵反卷积并利用了噪声的统计特性，如s5所示，使得环境信号的检测性能显著提升；
52.5、所述方法通过最大似然检测，如s6所示，直接使用总的接收信号作为判决门限，降低复杂度的同时，使得标签符号的检测性能及显著提升。
附图说明
53.图1是本发明基于ambc的双选信道lmmse估计及ml检测方法依托的双选信道下ambc系统模型图；
54.图2是本发明基于ambc的双选信道lmmse估计及ml检测方法的流程图；
55.图3是本发明基于ambc的双选信道lmmse估计及ml检测方法与现有的lmmse信道估计方法进行信道估计后的nmse(normalized mean squared error,归一化均方误差)对比图；
56.图4是本发明基于ambc的双选信道lmmse估计及ml检测方法与现有的lmmse信道估计及最大似然检测方法进行环境信号检测后的误码率对比图；
57.图5是本发明基于ambc的双选信道lmmse估计及ml检测方法与现有的lmmse信道估计及能量检测方法进行标签检测后的误码率对比图；
58.图6是本发明基于ambc的双选信道lmmse估计及ml检测方法与现有方法进行标签符号检测的误码率对比图；
59.图7是本发明基于ambc的双选信道lmmse估计及ml检测平台的示意图。
具体实施方式
60.下面结合附图和实施例对本发明基于ambc的双选信道lmmse估计及ml检测方法做进一步说明和详细描述。
61.实施例1
62.本方法可用于将无人机或者高速移动台等高速移动的信号源作为标签的能量来源的环境反射散射系统。首先针对环境信号和标签信号设计了块状导频簇，然后用bem模型对直接链路信道和反射链路信道建模并计算观测方程；根据导频结构及观测方程，首先用lmmse估计直接链路信道，然后将直链信号消除后估计反射链路信道；在信道估计的基础上，采用mmse检测算法检测环境信号；最后采用最大似然算法检测标签信号。所述方法的准确度显著提升，显著降低了信道估计方法的复杂度，实现了对信道时变特性的有效跟踪，信道估计及检测性能显著提升。
63.设定环境射频源ofdm的子载波数n为512，子载波调制方式为qpsk(quadrature phase shift keying，正交相移键控)；从环境射频源到接收端的信道为5径双选信道，多径功率服从指数分布，多普勒频移为120hz；从环境射频源到标签的信道为4径双选信道，多径功率服从指数分布，多普勒频移为120hz；从标签到接收端的信道为平坦衰落信道，信道增益为1；bem模型采用ce-bem(complex exponential bem，复指数bem信道模型)；标签反射系数α为0.1；具体实施步骤如下：
64.s1、基于块状导频簇结构分别构建环境信号和标签信号，具体为：构建以2个符号为一组的块状导频簇，并基于块状导频簇结构分别构建环境信号和标签信号，每间隔4个数据符号后插入一组块状导频簇；对于一组块状导频簇，前一个符号持续时间内标签信号为导频“0”，后一个符号持续时间内标签信号为导频“1”；即通过设计基于块状导频簇的通信协议：每间隔4个数据符号后插入一组块状导频簇，实现了对信道的时变特性进行跟踪，这使得信道估计及信号检测算法的准确度显著提升；
65.s2、环境信号经直接链路信道到达接收端得到直接链路信号；环境信号经过反射链路信道到达标签，标签吸收环境信号的能量，并根据自身携带的比特信息选择对射频源信号反射或者不反射，从而发送2ask符号经过平坦衰落信道到达接收端得到反射链路信号；接收端同时收到直接链路信号、反射链路信号和高斯白噪声得到总的接收信号；
66.即标签根据自身携带的比特信息选择对射频源信号进行调制，得到调制后符号，再将调制后符号经过平坦衰落信道到达接收端得到反射链路信号；接收端同时收到直接链路信号、反射链路信号和高斯白噪声得到总的接收信号；
67.其中，当标签准备发送符号“0”时，标签吸收射频源信号的能量后不对射频源信号进行反射；当标签准备发送符号“1”时，标签吸收射频源信号的能量后对射频源信号进行反射。
68.其中，总的接收信号包括环境信号经过直接链路信道到达接收端得到直接链路信号和经过反射链路信道到达标签后在通过平坦衰落信道到达接收端得到的反射链路信号以及高斯白噪声wi；
69.其中，是直接链路时域信道矩阵，代表第i个ofdm符号持续时间内第l条路径上第n个采样点的直接链路信道增益，n＝0,1,l,n-1；l＝0,1,l,l
d-1，ld为直接链路信道的多径数；是反射链路时域信道矩阵，代表第i个ofdm符号持续时间内第l条路径上第n个采样点的反射链路信道增益，n＝0,1,l,n-1；l＝0,1,l,l
r-1，lr为直接链路信道的多径数；α为反射链路能量衰减因子；f为平坦衰落信道的信道增益；ci为第i个标签符号；
70.s3、采用bem模型对直接链路信道和反射链路信道进行建模，然后分别构建直接链路观测方程和反射链路观测方程；
71.s3中，bem模型为时域信道矩阵，该矩阵为bem基矩阵与bem基的系数向量相乘的形式；
72.s3包括如下子步骤：
73.s31、采用bem模型对直接链路信道的单位冲激响应进行建模，具体为：
74.s311、将直接链路信道第i个符号持续时间内第l条路径上的第n个采样点用bem模型表示为
75.其中，bk＝[b
k,0
,k,b
k,q-1
]
t
为bem基向量，q为基向量的维度，为第i个符号持续时间内直接链路信道第l条路径的bem基的系数向量；
[0076]
s312、将s311中的扩展到n个采样点和ld条路径后有其中表示第i个符号上第l条路径上的直接链路cir向量；cir向量；为维度为ld×
ld的单位矩阵，为克罗克内积；为第i个符号持续时间内直接链路信道的bem基的系数矩阵。
[0077]
s32、在bem模型的基础上构建直接链路信道的观测方程，具体如下：
[0078]
s321、将第i个符号周期第l条路径上的环境射频源ofdm时域信号向量用对角化表示
[0079]
s322、将s321中直接链路信道所对应的合并为
[0080]
s323、对s322中的做离散傅里叶变换后与s312中的b相乘得到直接链路的观测矩阵观测方程表示为其中为直接链路的频域接收信号。
[0081]
s33、采用bem模型对反射链路信道的单位冲激响应进行建模，具体为：
[0082]
s331、将反射链路信道第i个符号持续时间内第l条路径上的第n个采样点用bem模型表示为其中，bk＝[b
k,0
,k,b
k,q-1
]
t
为bem基向量，q为基向量的维度，为第i个符号持续时间内直接链路信道第l条路径的bem基的系数向量；
[0083]
s332、将s331中的扩展到n个采样点和lr条路径后有其中表示第i个符号上第l条路径上的反射链路cir向量；cir向量；为维度为lr×
lr的单位矩阵，为克罗克内积；为第i个符号持续时间内反射链路信道的bem基的系数矩阵。
[0084]
s34、在bem模型的基础上构建反射链路信道的观测方程，具体如下：
[0085]
s341、将第i个符号周期第l条路径上的环境射频源ofdm时域信号向量用对角化表示
[0086]
s342、将s341中从第1条径到第lr条径所对应的并入一个矩阵
[0087]
s343、对s342中的做离散傅里叶变换后与s332中的b相乘得到反射链路的观测矩阵观测方程表示为其中为反射链路的频域接收信号。
[0088]
因此，s3采用bem模型对信道单位冲激响应进行建模，在模型不丢失双选信道性质的情况下减少了信道待估计参数的数量，显著降低了信道估计方法的复杂度；
[0089]
s4、在s3计算得到的直接链路观测方程和反射链路观测方程的基础上进行信道估计，得到直接链路时域信道矩阵和反射链路时域信道矩阵，具体为：分别对s3直接链路观测方程和反射链路观测方程中的直接链路bem基的系数向量及反射链路bem基的系数向量进行估计，然后再把bem基的系数向量恢复为直接链路时域信道矩阵和反射链路时域信道矩阵，包括如下子步骤：
[0090]
s41、根据s3计算得到的直接链路观测方程进行直接链路信道估计，得到直接链路时域信道矩阵，具体为：
[0091]
s411、在标签发送导频符号“0”的时刻用lmmse估计直接链路时域信道矩阵，具体如下：
[0092]
s4111、首先根据s3的直接链路观测矩阵和直接链路频域接收信号估计出
bem基的最小二乘系数向量估计值
[0093]
s4112、基于s312中的b与信道时域相关矩阵r
hh
计算出bem基的系数向量相关矩阵r
cc
＝[(bh×
b)-1
×bh
]
×rhh
×
[b
×
(bh×
b)-1
]；
[0094]
s4113、计算出映射到bem基空间后的噪声协方差矩阵其中为噪声功率，i
nl
为维度为n乘l的单位阵；
[0095]
s4114、计算出bem基的lmmse系数向量估计值
[0096]
s4115、将s332中的b与s414估计出的bem基的lmmse系数向量相乘得到直接链路时域信道矩阵；
[0097]
s412、对标签发送导频符号“1”以及数据符号的时刻用信道插值计算出直接链路时域信道矩阵；
[0098]
s42、根据s3计算得到的反射链路观测方程进行反射链路信道估计，得到反射链路时域信道矩阵，具体为：
[0099]
s421、在标签发送导频符号“1”的时刻用估计出的直接链路时域信道矩阵与ofdm导频符号相乘恢复出直接链路接收信号；
[0100]
s422、在标签发送导频符号“1”的时刻用总的接收信号减去s43恢复出的直接链路接收信号从而获得剩余信号；
[0101]
s423、在标签发送导频符号“1”的时刻将s44计算得到的剩余信号做离散傅里叶变换得到反射链路频域接收信号
[0102]
s424、在标签发送导频符号“1”的时刻用lmmse估计反射链路时域信道矩阵，具体如下：
[0103]
s4241、首先根据s343的反射链路观测矩阵和s423中的反射链路频域接收信号估计出bem基的最小二乘系数向量估计值：
[0104][0105]
s4242、基于s312中的b与信道时域相关矩阵r
hh
计算出bem基的系数向量相关矩阵r
cc
＝[(bh×
b)-1
×bh
]
×rhh
×
[b
×
(bh×
b)-1
]；
[0106]
s4243、计算出映射到bem基空间后的噪声协方差矩阵其中为噪声功率，i
nl
为维度为n
×
l的单位阵；
[0107]
s4244、计算出bem基的lmmse系数向量估计值
[0108]
s4245、将s332中的b与s4244估计出的bem基的lmmse系数向量相乘得到发射链路时域信道矩阵；
[0109]
s425、对标签发送导频符号“0”以及数据符号的时刻用信道插值估计出反射链路时域信道矩阵；
[0110]
s4对环境反向通信系统的双选信道进行估计，实现了对信道时变特性的有效跟踪，使得反射链路和直射链路的信道估计性能显著提升，
[0111]
具体包括如下子步骤：
[0112]
s41、根据s3计算得到的直接链路观测方程进行直接链路信道估计，得到直接链路时域信道矩阵；
[0113]
s42、根据s3计算得到的反射链路观测方程进行反射链路信道估计，得到反射链路时域信道矩阵；
[0114]
s41、具体为：
[0115]
s411、在标签发送导频符号“0”的时刻用lmmse估计出直接链路的bem基的系数向量；
[0116]
s412、将s3计算得到的bem基矩阵与s411计算得到的直接链路的bem基的系数向量相乘得到“0”时刻的直接链路时域信道矩阵；
[0117]
s413、在标签发送导频符号“1”以及数据符号的时刻用信道插值计算出“1”时刻及数据符号时刻的直接链路时域信道矩阵。
[0118]
s42、具体为：
[0119]
s421、在标签发送导频符号“1”的时刻用估计出的直接链路时域信道矩阵与ofdm导频符号相乘恢复出“1”时刻的直接链路接收信号；
[0120]
s422、在标签发送导频符号“1”的时刻用总的接收信号减去s421恢复出的直接链路接收信号从而获得剩余信号；
[0121]
s423、在标签发送导频符号“1”的时刻将s422计算得到的剩余信号做离散傅里叶变换得到反射链路频域接收信号y
ri
；
[0122]
s424、在标签发送导频符号“1”的时刻用lmmse估计出反射链路的bem基的系数向量；
[0123]
s425、将s3计算得到的bem基矩阵b与s424计算得到的反射链路的bem基的系数向量相乘得到1时刻的反射链路时域信道矩阵；在标签发送导频符号“0”以及数据符号的时刻用信道插值估计出“0”时刻及数据符号时刻的反射链路时域信道矩阵。
[0124]
s5、得到s4的直接链路时域信道矩阵后，对射频源信号进行mmse检测，得到ofdm符号估计值，具体为：
[0125]
s51、对总的时域接收信号做离散傅里叶变换后得到频域接收信号yi；
[0126]
s52、将反射链路信号视为干扰，对环境信号ofdm符号进行mmse均衡得到其中为直接链路频域响应矩阵；
[0127]
s 53、对s52计算得到的均衡结果进行星座图映射从而解调得到ofdm符号检测结果；
[0128]
s5通过mmse检测算法对时变信道下的射频源ofdm符号进行检测，对带状信道频域矩阵反卷积并利用了噪声的统计特性，使得环境信号的检测性能显著提升。
[0129]
s6、采用最大似然检测法对标签信号进行检测，得到标签符号的检测结果，具体如下：
[0130]
s61、分别对s4的直接链路时域信道矩阵和反射链路时域信道矩阵左乘一个傅里叶矩阵后，再右乘一个“傅里叶矩阵的共轭转置”，得到直接链路频域信道矩阵和反射链路
频域信道矩阵；
[0131]
s62、将反射链路频域信道矩阵与反射系数相乘后再与直接链路频域信道矩阵相加得到复合链路频域信道矩阵；
[0132]
s63、分别计算第i个数据符号时刻的直接链路频域信道矩阵和复合链路频域信道矩阵与s5得到的ofdm符号估计值相乘得到y0和y1；
[0133]
s64、对总的接收信号做傅里叶变换得到频域接收信号，其中第i个数据符号时刻的频域接收信号为yi；
[0134]
s65、分别计算yi与y0和y1的欧几里得距离d0和d1，若d0小于等于d1，则输出当前标签符号的检测结果为“0”；若d0大于d1，则输出当前标签符号的检测结果为“1”。
[0135]
s6直接使用总的接收信号作为判决门限，降低复杂度的同时，使得标签符号的检测性能及显著提升。至此，双选信道下ambc的信道估计及信号检测完成，然后对信道估计得到的两条链路的时域信道矩阵与实际的时域信道矩阵进行nmse的计算，具体为:“估计得到矩阵与真实矩阵的各个位置元素之差的平方和”除以“真实矩阵各个位置元素的平方和”。
[0136]
图3是采用本发明所述信道估计方法进行直接链路信道估计得到的时域信道矩阵与实际时域直接链路信道矩阵的nmse和采用现有方法进行直接链路信道估计得到的时域信道矩阵与实际时域直接链路信道矩阵的nmse的结果对比图；其中，现有方法在双选信道条件下nmse维持在0.05左右；采用本发明所述方法在双选信道条件下nmse随着直接链路信噪比的提升也逐渐降低，在直接链路信噪比为35db时，nmse达到了4
×
10-5
。
[0137]
图4是采用本发明所述信道估计方法进行反射链路信道估计得到的时域信道矩阵与实际时域反射链路信道矩阵的nmse和采用现有方法进行反射链路信道估计得到的时域信道矩阵与实际时域反射链路信道矩阵的nmse的结果对比图；其中，现有方法在双选信道条件下nmse随着直接链路信噪比的提升稍微降低，在直接链路信噪比为35db时，nmse为10-1
；采用本发明所述方法在双选信道条件下nmse随着直接链路信噪比的提升也逐渐降低，在直接链路信噪比为35db时，nmse达到了2
×
10-3
。
[0138]
图5是采用所述方法与现有方法进行射频源ofdm信号检测的误码率对比图；采用所述方法的误码率性能相较于现有方法提升了1db左右，在直接链路信噪比为35db时，所述方法的ofdm符号误码率达到了2
×
10-2
。
[0139]
图6是采用所述方法与现有方法进行标签符号检测的误码率对比图；其中采用所述方法的误码率性能相较于现有方法提升了10db左右。在直接链路信噪比为35db时，现有方法标签符号误码率仅为9
×
10-1
，所述方法的ofdm符号误码率达到了9
×
10-2
。
[0140]
实施例2
[0141]
如图7所示，本发明第2实施方式涉及基于ambc的双选信道lmmse估计及ml检测平台，包括至少一个处理器，以及与至少一个处理器通信连接的存储器；其中，存储器存储有可被至少一个处理器执行的指令，指令被至少一个处理器执行，以使至少一个处理器能够执行上述基于ambc的双选信道lmmse估计及ml检测方法。
[0142]
其中，存储器和处理器采用总线方式连接，总线可以包括任意数量的互联的总线和桥，总线将一个或多个处理器和存储器的各种电路连接在一起。总线还可以通过接口将诸如外围设备、稳压器和功率管理电路等之类的各种其他电路连接在一起，这些都是本领域所公知的。接口在总线和收发机之间提供接口，例如通信接口、用户接口。收发机可以是
一个元件，也可以是多个元件，比如多个接收器和发送器，提供用于在传输介质上与各种其他装置通信的单元。经处理器处理的数据通过天线在无线介质上进行传输，进一步，天线还接收数据并将数据传送给处理器。处理器负责管理总线和通常的处理，还可以提供各种功能，包括定时，外围接口，电压调节、电源管理以及其他控制功能。而存储器可以被用于存储处理器在执行操作时所使用的数据。
[0143]
实施例3
[0144]
本发明第3实施例涉及一种计算机可读存储介质，存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例。
[0145]
本领域技术人员通过上述说明可以理解，实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序存储在一个存储介质中，包括若干指令用以使得一个设备(可以是单片机，芯片等)或处理器(processor)执行本技术各个实施例所述方法的全部或部分步骤。而前述的存储介质包括但不限于u盘、移动硬盘、磁性存储器、光学存储器等各种可以存储程序代码的介质。
[0146]
通过上述说明可知，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。以上所述为本发明的较佳实施例而已，本发明不应该局限于该实施例和附图所公开的内容。凡是不脱离本发明所公开的精神下完成的等效或修改，都落入本发明保护的范围。