本发明涉及移动通信中的信道估计领域。本发明的基于序参量非正交变换的通信系统盲信道估计方法，包括如下步骤：步骤一：计算MIMO-OFDM通信系统信道的相关矩阵；步骤二：选取接收信号的序参量，根据接收信号Y(t)={y(t1)，y(t2)，…y(tM)}，选取其二阶矩的N个主要特征值及其对应的特征向量作为序参量；步骤三：估计MIMO-OFDM通信系统非正交信道的特性，即其各个信道的冲激响应，求解非正交基函数；步骤四：将信号按照非正交基函数进行非正交分解，通过盲分离估计出源信号，其分解系数x1，…，xN即是分离出来的源信号向量本发明不需要在发送端加入导频序列，只利用接收端的接收信号来对信道进行盲估计，能有效的节省带宽，并且计算复杂度较低，估计精度较高。
1.一种基于序参量非正交变换的通信系统盲信道估计方法，其特征在于，包括如下步骤：步骤一：计算MIMO-OFDM通信系统信道的相关矩阵，2
将正定矩阵Rsp分解为正定矩阵U的平方，即Rsp=U，其中， 表示发送端天线的相关矩阵， 表示接收端天线的相关矩阵，Rp表示路径相关矩阵；
步骤二：选取接收信号的序参量，根据接收信号Y(t)={y(t1)，y(t2)，…y(tM)}，选取其二阶矩的N个主要特征值及其对应的特征向量作为序参量，使得他们所包含的信息熵为信号总信息熵的90％，包含了源信号的主要信息，剩下的M-N个次要特征值及其对应的特征向量则包含了加性噪声的主要信息，占信号总信息熵的10％；
步骤三：估计MIMO-OFDM通信系统非正交信道的特性，即其各个信道的冲激响应，求解非正交基函数，其中，hi(t)，i=1，2，…，N为非正交基，Rsp是hi(t)，i=1，2，…，N的相关函数矩阵，如果Rsp是单位矩阵I，则相应的基即为单位正交基，说明信道之间不存在相关性，是完全正交信道；如果Rsp不是单位矩阵，则说明信道之间具有相关性；
步骤四：将信号按照非正交基函数 进行非正交分解，通过盲分离估计出源信号，其分解系数x1，…，xN即是分离出来的源信号向量
2.根据权利要求1所述的基于序参量非正交变换的通信系统盲信道估计方法，其特征在于，步骤一所述的计算MIMO-OFDM通信系统信道的相关矩阵，按如下步骤进行：通过Kronecker内积来计算如下公式中的空间、路径相关矩阵，其中， 表示发送端天线的相关矩阵， 表示接收端天线的相关矩阵，Rp表示路径相关矩阵，通过MIMO-OFDM系统散射体模型其中， Rp的计算按如下方法：(1)设信号到达接收天线的第p条路径的角度为 相移 由下式计算：其中，λ为载波波长，p表示第p条路径，r表示到r个阵元， 为相邻阵元之间的距离，则接收天线阵列的相移用向量 表示：(2)发射天线阵列的相移用向量 表示：
其中，式中 是功率延迟截面PDP，γm／γm+1=3dB,(m=0，…，L-1)， p是路径相关因子。
3.根据权利要求1所述的基于序参量非正交变换的通信系统盲信道估计方法，其特征在于，步骤二所述的选取接收信号的序参量，按如下步骤进行：(1)根据接收信号Y(t)={y(t1)，y(t2)，…y(tM)}，令M>N估计其相关矩阵，Y(s)表示接收信号在时刻s点的值；
是所有序参量构成的函数集分别在时刻s、t的表达式；
(3)对λ由大到小排序，选择N个最大λ值作为序特征值，与之相对应有包含N个是序参量，使得他们所包含的信息熵为信号总信息熵的90％，余下10％的信息熵作为加性噪声的信息熵，特征值之间通常有典型的关系λ1≥λ2≥…≥λN≥λN+1≥…≥λM，这里的M>N，选取的N个主特征值及其对应的特征向量包含了源信号的主要信息，剩下的M-N个次要特征值及其对应的特征向量则包含了加性噪声的主要信息；
(4)通过特征值的个数可以估计源信号的个数，通过寻找序参量的方法可以将信号的主信号空间与噪声信号空间分离开来。
4.根据权利要求1所述的基于序参量非正交变换的通信系统盲信道估计方法，其特征在于，步骤三所述的估计MIMO-OFDM通信系统非正交信道的特性，即其各个信道的冲激响应，按如下步骤进行：将 两边左乘U并做如下变形，
即满足方程 为序参量，φ是求序参量运算的算子形式，U是由前2
面求出的信道相关矩阵Rsp计算得到，Rsp＝U，Q＝Rsp＝UU， 即为估计的各个信道的冲激响应。
5.根据权利要求书中所述的基于序参量非正交变换的通信系统盲信道估计方法，其特T征在于，步骤四所述的将信号按照非正交基函数H(t)＝(h1(t)，h2(t)，...，hN(t)) 进行非正交分解，通过盲分离估计出源信号，按如下步骤进行：表示将接收信号分解在各个信道上，每个信道上所包含的信号成分，通过上式消除相关性，求出的X＝(x1，...，xN)即为估计出的源信号。
基于序参量非正交变换的通信系统盲信道估计方法
技术领域
[0001] 本发明涉及移动通信中的信道估计领域，尤其涉及一种应用于多输入多输出-正交频分复用通信系统盲信道估计方法。
背景技术
[0002] 在MIMO-OFDM(MIMO-OFDM通信系统即为多输入多输出-正交频分复用通信系统)系统中不同的天线发送不同的信号，每一个天线、每一个子载波都有多个信道参数，信道估计比较困难。当信道变化较快时，由于通过信道估计所得到的信道状态信息会很快过时，因此为了得到信道的瞬时状态信息需要付出很大的代价，使得系统的传输效率降低，同时也会使得误码率增高。针对目前MIMO-OFDM系统中信道估计存在的困难和问题，急需一种方案来达到更好的信道估计效果。
发明内容
[0003] 本发明的目的是解决在信道估计领域中，特别在MIMO-OFDM系统中，由于天线的空间相关性，信号传输的路径相关性，使得信道特性的估计系统比较复杂以及估计结果不够准确的问题，提供一种基于序参量非正交变换的通信系统盲信道估计方法。
[0004] 为了达到上述目的，本发明所采用的技术方案是：包括如下步骤：
[0005] 步骤一：计算MIMO-OFDM通信系统信道的相关矩阵，
[0006]
[0007] 将正定矩阵Rsp分解为正定矩阵U的平方，即Rsp=U2，其中， 表示发送端天线的相关矩阵， 表示接收端天线的相关矩阵，Rp表示路径相关矩阵；
[0008] 步骤二：选取接收信号的序参量，根据接收信号Y(t)={y(t1)，y(t2)，…y(tM)}，选取其二阶矩的N个主要特征值及其对应的特征向量作为序参量，使得他们所包含的信息熵为信号总信息熵的90％，包含了源信号的主要信息，剩下的M-N个次要特征值及其对应的特征向量包含了加性噪声的主要信息，占信号总信息熵的10％，通过寻找序参量的方法将信号的主信号空间与噪声信号空间分离开来；
[0009] 步骤三：估计MIMO-OFDM通信系统非正交信道的特性，即其各个信道的冲激响应，求解非正交基函数，
[0010]
[0011] 为估计的各个信道的冲激响应，hi(t)，i＝1，2，…，N为非正交基，这些非正交基函数代表了信道的特性，Rsp是hi(t)，i＝1，2，…，N的相关函数矩阵，为正定矩阵，并且Rsp体现了信道函数间的关系，特殊地，如果Rsp是单位矩阵I，则相应的基即为单位正交基，说明信道之间不存在相关性，是完全正交信道，如果Rsp不是单位矩阵，则说明信道之间具有相关性；
[0012] 步骤四：将信号按照非正交基函数 进行非正交分解，通过盲分离估计出源信号，其分解系数x1，…，xN即是分离出来的源信号向量 。
[0013] 其中，步骤一所述的计算MIMO-OFDM通信系统信道的相关矩阵，按如下步骤进行：
[0014] 通过Kronecker内积即克罗内克内积来计算如下公式中的空间、路径相关矩阵，[0015]
[0016] 其中， 表示发送端天线的相关矩阵， 表示接收端天线的相关矩阵，Rp表示路径相关矩阵，Rp的计算既依赖于路径相关，也依赖于功率延迟截面PDP，也就是说，即使在系统中没有路径相关，Rp也不一定是单位矩阵，这是因为Rp对角线上的元素描述了每个径依赖于功率延迟截面的平均能量，其中， Rp的计算按如下方法：
[0017] (1)设信号到达接收天线的第p条路径的角度为 相移 由下式计算：
[0018]
[0019] 其中，λ为载波波长，p表示第p条路径，r表示到r个阵元， 为相邻阵元之间的距离，则接收天线阵列的相移用向量 表示：
[0020]
[0021] (2)发射天线阵列的相移用向量 表示：
[0022]
[0023] 于是得到：
[0024]
[0025]
[0026] (3)就可以得到相关矩阵的表达式为：
[0027]
[0028] (4)通过下式计算路径相关矩阵：
[0029]
[0030] 其中，式中 是功率延迟截面PDP，γm／γm+1=3dB，(m=0，…，L-1)，p是路径相关因子。
[0031] 步骤二所述的选取接收信号的序参量，按如下步骤进行：
[0032] (1)根据接收信号Y(t)={y(t1)，y(t2)，…y(tM)}，令M>N估计其相关矩阵，[0033] 为接收信号的相关矩阵，Y(s)表示接收信号在时刻s点的值；
[0034] (2)解下列积分方程，
[0035]
[0036] 设 则
[0037] 其 中 ， 为 相 应 的 对 角 阵是 其 中 一 个 序 参 量，
是所有序参量构成的函数集分别在时刻s、t的表达式；
[0038] (3)对λ由大到小排序，选择N个最大λ值作为序特征值，与之相对应有包含N个 是序参量，使得他们所包含的信息熵为信号总信息熵的90％，余下10％的信息熵作为加性噪声的信息熵，特征值之间通常有典型的关系λ1≥λ2≥…≥λN≥λN+1≥…≥λM(M＞N)，选取的N个主特征值及其对应的特征向量包含了源信号的主要信息，而剩下的M-N个次要特征值及其对应的特征向量则包含了加性噪声的主要信息；
[0039] (4)通过特征值的个数可以估计源信号的个数，通过寻找序参量的方法可以将信号的主信号空间与噪声信号空间分离开来。
[0040] 步骤三所述的估计MIMO-OFDM通信系统非正交信道的特性，即其各个信道的冲激响应，按如下步骤进行：
[0041] 将 对角变换后两边左乘U并做如下变形
[0042]
[0043] 取
[0044] 即满足方程
[0045] 为序参量，Φ是求序参量运算的算子形式，U是由前面求出的信道相关矩阵2
Rsp计算得到，Rsp=U，Q=Rsp=UU， 即为估计的各个信道的冲激响
应。。
[0046] 步骤四所述的将信号按照非正交基函数 进行非正交分解，通过盲分离估计出源信号，按如下步骤进行：
[0047]
[0048] 表示将接收信号分解在各个信道上，每个信道上所包含的信号成分，通过上式消除相关性，求出的X=(x1，…，xN)即为估计出的源信号。
[0049] 本发明不需要在发送端加入导频序列，只利用接收端的接收信号来对信道进行盲估计，本发明能有效的节省带宽，提高数据传输效率，并且计算复杂度较低，估计精度较高。
附图说明
[0050] 图1为本发明结构示意框图；
[0051] 图2为本发明与EM算法误码率性能比较；
[0052] 图3为本发明MIMO-OFDM系统的一种散射体模型；
[0053] 图4为本发明MIMO-OFDM系统的一种散射体几何模型。
具体实施方式
[0054] 下面结合附图对本发明作进一步的描述：
[0055] 本发明的基于序参量非正交变换的MIMO-OFDM通信系统盲信道估计方法，这里的MIMO-OFDM通信系统即为多输入多输出一正交频分复用通信系统，包括如下步骤：
[0056] 步骤一：计算MIMO-OFDM通信系统信道的相关矩阵，
[0057]
[0058] 将正定矩阵Rsp分解为正定矩阵U的平方，即，其中， 表示发送端天线的相关矩2
阵， 表示接收端天线的相关矩阵，Rp表示路径相关矩阵Rsp＝U ；
[0059] 步骤二：选取接收信号的序参量，根据接收信号Y(t)={y(t1)，y(t2)，…y(tM)}，选取其二阶矩的N个主要特征值及其对应的特征向量作为序参量，使得他们所包含的信息熵为信号总信息熵的90％，包含了源信号的主要信息，剩下的M-N个次要特征值及其对应的特征向量则包含了加性噪声的主要信息，占信号总信息熵的10％，通过寻找序参量的方法将信号的主信号空间与噪声信号空间分离开来；
[0060] 步骤三：估计MIMO-OFDM通信系统非正交信道的特性，即其各个信道的冲激响应，求解非正交基函数，采用如下公式：
[0061]
[0062] 为估计的各个信道的冲激响应，hi(t)，i＝1，2，…，N为非正交基，这些非正交基函数代表了信道的特性，Rsp是hi(t)，i＝1，2，…，N的相关函数矩阵，为正定矩阵，并且Rsp体现了信道函数间的关系，特殊地，如果Rsp是单位矩阵I，则相应的基即为单位正交基，说明信道之间不存在相关性，是完全正交信道，如果Rsp不是单位矩阵，则说明信道之间具有相关性；
[0063] 步骤四：将信号按照非正交基函数 进行非正交分解，通过盲分离估计出源信号，其分解系数x1，…，xN即是分离出来的源信号向量[0064] 其中，步骤一所述的计算MIMO-OFDM通信系统信道的相关矩阵，按如下步骤进行：
[0065] 通过Kronecker内积即克罗内克内积来计算如下公式中的空间、路径相关矩阵，[0066]
[0067] 其中， 表示发送端天线的相关矩阵， 表示接收端天线的相关矩阵，Rp表示路径相关矩阵，Rp的计算既依赖于路径相关，也依赖于功率延迟截面PDP，也就是说，即使在系统中没有路径相关，Rp也不一定是单位矩阵，这是因为Rp对角线上的元素描述了每个径依赖于功率延迟截面的平均能量，其中， Rp的计算按如下方法：
[0068] (1)设信号到达接收天线的第p条路径的角度为 相移 由下式计算：
[0069]
[0070] 其中，λ为载波波长，p表示第p条路径，r表示到r个阵元， 为相邻阵元之间的距离，其中， 为相邻阵元之间的距离，则接收天线阵列的相移用向量 表示：
[0071]
[0072] (2)发射天线阵列的相移用向量 表示：
[0073]
[0074] 于是得到：
[0075]
[0076]
[0077] (3)就可以得到相关矩阵的表达式为：
[0078]
[0079] (4)通过下式计算路径相关矩阵：
[0080]
[0081] 其中，式中 是功率延迟截面PDP，γm/γm+1=3dB，(m=0，…，L-1)，p是路径相关因子。
[0082] 步骤二所述的选取接收信号的序参量，按如下步骤进行：
[0083] (1)根据接收信号Y(t)={y(t1)，y(t2)，…y(tM)}，令M>N估计其相关矩阵，Y(s)表示接收信号在时刻s点的值， 为接收信号的相关矩阵；
[0084] (2)解下列积分方程，
[0085]
[0086] 设 则
[0087] 其 中，Λ 为 相 应 的 对 角 阵 Λ=diag(λ1，λ2， …，λN)，是其中一个序参量， 是所有序参量构成的函数集分别在时刻s、t的表达式；
[0088] (3)对λ由大到小排序，选择N个最大λ值作为序特征值，与之相对应有包含N个 是序参量，使得他们所包含的信息熵为信号总信息熵的90％，余下10％的信息熵作为加性噪声的信息熵，特征值之间通常有典型的关系λ1≥λ2≥…≥λN≥λN+1≥…≥λM(M>N)，选取的N个主特征值及其对应的特征向量包含了源信号的主要信息，而剩下的M-N个次要特征值及其对应的特征向量则包含了加性噪声的主要信息；
[0089] (4)通过特征值的个数可以估计源信号的个数，通过寻找序参量的方法可以将信号的主信号空间与噪声信号空间分离开来。
[0090] 步骤三所述的估计MIMO-OFDM通信系统非正交信道的特性，即其各个信道的冲激响应，按如下步骤进行：
[0091] 将 对角变换后两边左乘U并做如下变形
[0092]
[0093] 取
[0094] 即满足方程
[0095] 为序参量，Φ是求序参量运算的算子形式，U是由前面求出的信道相关矩阵Rsp2
计算得到，Rsp=U，Q=Rsp=UU， 即为估计的各个信道的冲激响应。
[0096] 步骤四所述的将信号按照非正交基函数 进行非正交分解，通过盲分离估计出源信号，按如下步骤进行：
[0097]
[0098] 氛示将接收信号分解在各个信道上，每个信道上所包含的信号成分，通过上式消除相关性，求出的X=(x1，…，xN)即为估计出的源信号。
[0099] 本发明的效果通过仿真进一步说明：
[0100] 源信号是独立同分布的实随机变量，其均值为0、方差为1，每径的噪声是独立的加性零均值高斯白噪声，仿真中使用的采样数据序列长度N=2000，考虑收发天线组数为2、4、6时的MIMO-OFDM系统，SNR是Nt个发送天线的符号总能量与噪声能量的比值，对每个固定的信噪比，检测长为2000个码元的序列，结果取10次独立仿真的平均值。空间相关矩阵根据发送端余接收端的特定物理设置，发送端的天线阵列之间的距离为3λ，接收端的天线阵列之间的距离为0.5λ，λ为波长，路径相关矩阵Rp计算式为
[0101] 如图2所示，比较采用本发明的盲信道估计方法与基于EM算法的半盲信道估计方法，分别进行MIMO-OFDM系统的盲信道估计的误码率BER。从图中可以看出，随着信噪比的提高本发明的方法比基于EM的信道估计算法的误码率更加接近于实际信道的误码率。
