本发明属于无线通信技术领域，涉及一种盲信道估计模糊度消除方法，具体涉及一种用于非合作mimo通信基于信源特性的盲信道估计模糊度消除方法，可应用于信号识别等场景。
背景技术：
mimo(multiple-inputmultiple-output，多输入多输出)通信系统中，发射天线向接收天线发送信号时产生信道，对信道的信道估计是指接收端获得信道状态信息(channelstateinformation，csi)的过程和方法。信道估计方法分为：利用导频信息的非盲信道估计方法、不使用导频信息的盲信道估计方法以及使用部分导频信息的半盲信道估计方法。mimo通信系统又可分为合作通信场景和非合作通信场景。在合作通信场景下，接收端需要利用导频信息进行信道估计来获得信道状态信息，然后使用信道状态信息来恢复出发射端的信号。在非合作mimo通信场景下，接收端在对信道进行信道估计时，无法获取合作收发端的导频信息，因此无法使用非盲信道估计方法及半盲信道估计方法，只能采用盲信道估计方法。但是，在使用盲信道估计方法得到盲信道估计值时，存在天线排序模糊度问题和符号模糊度问题，天线模糊度问题指的是天线可能存在某些天线排列顺序错误，符号模糊度问题指的是每个符号可能存在错误的相位旋转。如果不能将模糊度去除，则严重影响信道的估计值，导致系统误码率的提高，恢复正确发射信号的概率降低，系统通信的可靠性较低，其中，系统误码率指的是恢复出的发射信号的出错数据量占总数据量的比例。
例如申请公布号为cn103763222a，名称为“一种mimo信号盲检测中的信道模糊度去除方法”的专利申请，公开了一种非合作mimo信号盲检测中的信道模糊度去除方法，该方法首先通过将mimo系统的信道划分为多个simo系统信道，消除天线排序所产生的天线排序模糊度影响，然后对各个simo系统信道进行盲信道估计，再通过均衡和估计后将这simo系统信道参数组成mimo信道矩阵，获得mimo信道估计值。该方法主要用来降低系统误码率，提高恢复出正确发射信号的概率，提升通信可靠性，但是其存在的不足之处在于，该方法只对天线排序模糊度进行了消除，缺少对符号模糊度的消除，导致系统误码率的降低不够理想。
技术实现要素：
本发明的目的在于针对上述现有技术存在的不足，提出了一种用于非合作mimo通信基于信源特性的盲信道估计模糊度消除方法，用于解决现有技术中存在的系统误码率较高的技术问题。
为实现上述目的，本发明采取的技术方案包括如下步骤：
(1)获取非合作mimo系统发射天线序号排列的顺序序列：
对包括n个发射天线和m个接收天线的非合作mimo系统中的每个发射天线进行标注，并对n个标注后的发射天线进行全排列，得到n！个发射天线序号排列的顺序序列，其中，n≥2，m≥2；
(2)构建发射天线排序初等矩阵集合hs：
(2a)构建排列单位行向量集合i＝{i1,i2,…，ii,…in}，其中，ii表示第i个单位行向量，ii＝[c1,c2,…,cj,…cn]，cj表示第j个元素的值，当j＝i时cj＝1，当j≠i时cj＝0；
(2b)按发射天线序号排列的顺序序列对i进行排列，得到发射天线排序初等矩阵集合
(3)构建模糊度遍历矩阵集合h：
(3a)设符号遍历系数为α，符号遍历列向量集合为m，并令α＝0，m＝{}；
(3b)对α个-1和n-α个1进行全排列，得到子符号遍历列向量集合
(3c)令α 1，并判断α＝n 1是否成立，若是，得到符号遍历列向量集合m，
(3d)以mr中的n个元素为主对角线、以0为其余元素构建符号遍历对角矩阵
(3e)计算hs与hc的笛卡尔积，得到过渡矩阵组集合h′，
(4)构建发射信号估计矩阵
(4a)将像素大小为x×z的原始灰度图像中的每个像素点转换为二进制像素点，转换后二进制像素点的位数若不满8位，在其左侧补零，得到像素大小为x×(z×8)的灰度图像，并将其作为每个发射天线的发射信号，发送至每个接收天线，其中，x≥28,z≥28；
(4b)对m个接收天线对接收到的n个发射信号进行组合，得到大小为m×(x×z×8/n)的接收信号矩阵y；
(4c)对y进行盲信道估计，得到大小为m×n的盲信道估计矩阵a，并通过a的取逆结果a-1和接收信号矩阵y，计算发射信号估计矩阵
(5)构建带模糊度的发射信号估计矩阵集合p：
对h中的每一个模糊度遍历矩阵hg与发射信号估计矩阵
其中，pg表示第g个带模糊度的发射信号估计矩阵，
(6)计算奇异值矩阵集合s：
(6a)设像素估计矩阵集合u，并令g＝1，u＝{}；
(6b)按照从前到后的顺序对pg中x×z×8/n个发射信号估计列向量进行拼接，得到长度为x×z×8的估计列向量
(6c)将
(6d)以z为单位长度对vg′按从上至下的顺序进行切分，得到包括x个像素点估计子向量的像素点估计子列向量集合
(6e)将vg″′中由前到后排列的x个转置子向量转化为自上而下排列的像素估计矩阵
(6f)对u内每个像素估计矩阵ug进行奇异值分解，得到奇异值矩阵集合
(7)基于信源特性计算消除模糊度后的盲信道估计矩阵adone：
(7a)设累积速度值为ρ，奇异值矩阵序号值为
(7b)构造长度均为q的奇异值向量eg和累积量向量tg：
其中，
(7c)通过tg以及eg中所有元素的和sg计算信源特性向量
其中，
(7d)令d＝1；
(7e)判断
(7f)判断d＜ρ是否成立，若是，令ρ＝d，
(7g)令d＝d 1，执行步骤(7e)；
(7h)判断g＝2n×n！是否成立，若是，将p中的第
(7i)通过pdone的取逆结果(pdone)-1与步骤(4b)中的接收信号矩阵y计算消除模糊度后的盲信道估计矩阵adone＝y×(pdone)-1。
本发明与现有技术相比，具有以下优点：
本发明组合了每种可能存在的天线排序模糊度和每种可能存在的符号模糊度得到像素估计矩阵集合，再将每个像素估计矩阵的奇异值累积量向量与其奇异值的和的比值作为信源特性向量，通过比较每个信源特性向量，找到信源特性最快到达0.7的信源特性向量，通过该信源特性向量计算消除模糊度后的盲信道估计矩阵。避免了现有技术仅能消除天线排序模糊度导致的系统误码率较高的缺陷，有效提高了恢复出正确发射信号的概率，系统通信的可靠性较高。
附图说明
图1是本发明的实现流程图；
图2是本发明中非合作mimo系统的结构示意图。
具体实施方式
下面结合附图和具体实施例，对本发明作进一步详细说明。
参照图1，本发明包括如下步骤：
步骤1)获取非合作mimo系统发射天线序号排列的顺序序列：
构建如图2所示的非合作mimo系统，该系统包括n个发射天线和m个接收天线，p1,p2,…,pn分别对应每个发射天线发送的发射信号，y1,y2,…,yn分别对应每个接收天线接收的n个发射信号，非合作是指接收方无法对其进行半盲信道估计或非盲信道估计，只能对接收信号矩阵进行盲信道估计得到盲信道估计矩阵；每个发射天线分别发送发射信号至每个接收天线，每个接收天线接收到的n个发射信号进行组合得到接收信号矩阵，将非合作mimo系统中的每个发射天线进行标注，并对n个标注后的发射天线进行全排列，得到n！个发射天线序号排列的顺序序列，其中，n≥2，m≥2，在本实施例中，n＝2，m＝2；
步骤2)构建发射天线排序初等矩阵集合hs：
步骤2a)构建排列单位行向量集合i＝{i1,i2,…，ii,…in}，其中，ii表示第i个单位行向量，ii＝[c1,c2,…,cj,…cn]，cj表示第j个元素的值，当j＝i时cj＝1，当j≠i时cj＝0；
步骤2b)针对天线排序模糊度的可能情况，按发射天线序号排列的顺序序列对i进行排列，得到发射天线排序初等矩阵集合
步骤3)构建模糊度遍历矩阵集合h：
步骤3a)为解决符号模糊度问题，对所有可能的符号变化情况构造得到对应的符号遍历列向量集合，通过每个符号遍历列向量构造成符号遍历矩阵，每个符号遍历矩阵对应一种符号模糊度的可能性；设符号遍历系数为α，符号遍历列向量集合为m，并令α＝0，m＝{}；
步骤3b)对α个-1和n-α个1进行全排列，得到子符号遍历列向量集合
步骤3c)令α 1，并判断α＝n 1是否成立，若是，得到符号遍历列向量集合m，
步骤3d)以mr中的n个元素为主对角线、以0为其余元素构建符号遍历对角矩阵
步骤3e)综合考虑天线排序模糊度及符号模糊度的影响，根据表达式h′＝{(η,γ)|η∈hs∩γ∈hc}，计算hs与hc的笛卡尔积，得到过渡矩阵组集合h′，
步骤4)构建发射信号估计矩阵
步骤4a)将像素大小为x×z的原始灰度图像中的每个像素点转换为二进制像素点，转换后二进制像素点的位数若不满8位，在其左侧补零，得到像素大小为x×(z×8)的灰度图像，并将其作为每个发射天线的发射信号，发送至每个接收天线，其中，x≥28,z≥28，原始灰度图像也可以为彩色图像或音频，实际情况中，接收端通过信号检测判断发射信号的类型和大小后就可继续进行步骤(4b)；在本实施例中采用像素大小为512×512的原始灰度图像，原始灰度图像的每个像素点的值为0～255的十进制数，将每个像素点做二进制转换后得到的每个二进制像素点对应的最大位数为8，将所有二进制像素点以8位补齐，确保信号结构的统一性；
步骤4b)对m个接收天线对接收到的n个发射信号进行组合，得到大小为m×(x×z×8/n)的接收信号矩阵y；
步骤4c)对y进行盲信道估计，得到大小为m×n的盲信道估计矩阵a，并通过a的取逆结果a-1和接收信号矩阵y，计算发射信号估计矩阵
步骤5)构建带模糊度的发射信号估计矩阵集合p：
对h中的每一个模糊度遍历矩阵hg与发射信号估计矩阵
将处理盲信道估计矩阵模糊度的问题转化为处理带模糊度的发射信号中模糊度的问题，能够通过发射信号将带模糊度的发射信号构造为结构相同的矩阵形式，可直接按位对比带模糊度的发射信号与原始图像信号；pg表示第g个带模糊度的发射信号估计矩阵，
步骤6)计算奇异值矩阵集合s：
步骤6a)设像素估计矩阵集合u，并令g＝1，u＝{}；
步骤6b)每个接收天线均得到二进制数据，需要将其按原始灰度图像处理中的长度8重新切分，并经过十进制转换得到与原始灰度图像数据量一致的十进制数据；按照从前到后的顺序对pg中x×z×8/n个发射信号估计列向量进行拼接，得到长度为x×z×8的估计列向量
步骤6c)将
步骤6d)以z为单位长度对vg′按从上至下的顺序进行切分，得到包括x个像素点估计子向量的像素点估计子列向量集合
步骤6e)将得到的所有像素点按照原始灰度图像的尺寸，按先行后列的顺序重新排列为像素估计矩阵；将vg″′中由前到后排列的x个转置子向量转化为自上而下排列的像素估计矩阵
步骤6f)对u内每个像素估计矩阵ug进行奇异值分解，得到奇异值矩阵集合
步骤7)基于信源特性计算消除模糊度后的盲信道估计矩阵adone：
步骤7a)设累积速度值为ρ，奇异值矩阵序号值为
步骤7b)构造长度均为q的奇异值向量eg和累积量向量tg：
其中，
步骤7c)通过tg以及eg中所有元素的和sg计算信源特性向量
其中，
步骤7d)令d＝1；
步骤7e)判断
步骤7f)判断d＜ρ是否成立，若是，令ρ＝d，
步骤7g)令d＝d 1，执行步骤(7e)；
步骤7h)判断g＝2n×n！是否成立，若是，通过对所有信源特性向量进行比较后得到的
步骤7i)通过pdone的取逆结果(pdone)-1与接收信号矩阵y计算消除模糊度后的盲信道估计矩阵adone＝y×(pdone)-1。
以上是对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求范围内做出各种变形或修改，这并不影响本发明的实质内容。
1.一种用于非合作mimo通信的基于信源特性的盲信道估计模糊度消除方法，其特征在于，包括如下步骤：
(1)获取非合作mimo系统发射天线序号排列的顺序序列：
对包括n个发射天线和m个接收天线的非合作mimo系统中的每个发射天线进行标注，并对n个标注后的发射天线进行全排列，得到n！个发射天线序号排列的顺序序列，其中，n≥2，m≥2；
(2)构建发射天线排序初等矩阵集合hs：
(2a)构建排列单位行向量集合i＝{i1,i2,…，ii,…in}，其中，ii表示第i个单位行向量，ii＝[c1,c2,…,cj,…cn]，cj表示第j个元素的值，当j＝i时cj＝1，当j≠i时cj＝0；
(2b)按发射天线序号排列的顺序序列对i进行排列，得到发射天线排序初等矩阵集合
(3)构建模糊度遍历矩阵集合h：
(3a)设符号遍历系数为α，符号遍历列向量集合为m，并令α＝0，m＝{}；
(3b)对α个-1和n-α个1进行全排列，得到子符号遍历列向量集合
(3c)令α 1，并判断α＝n 1是否成立，若是，得到符号遍历列向量集合m，
(3d)以mr中的n个元素为主对角线、以0为其余元素构建符号遍历对角矩阵
(3e)计算hs与hc的笛卡尔积，得到过渡矩阵组集合h′，
(4)构建发射信号估计矩阵
(4a)将像素大小为x×z的原始灰度图像中的每个像素点转换为二进制像素点，转换后二进制像素点的位数若不满8位，在其左侧补零，得到像素大小为x×(z×8)的灰度图像，并将其作为每个发射天线的发射信号，发送至每个接收天线，其中，x≥28,z≥28；
(4b)对m个接收天线对接收到的n个发射信号进行组合，得到大小为m×(x×z×8/n)的接收信号矩阵y；
(4c)对y进行盲信道估计，得到大小为m×n的盲信道估计矩阵a，并通过a的取逆结果a-1和接收信号矩阵y，计算发射信号估计矩阵
(5)构建带模糊度的发射信号估计矩阵集合p：
对h中的每一个模糊度遍历矩阵hg与发射信号估计矩阵
其中，pg表示第g个带模糊度的发射信号估计矩阵，
(6)计算奇异值矩阵集合s：
(6a)设像素估计矩阵集合u，并令g＝1，u＝{}；
(6b)按照从前到后的顺序对pg中x×z×8/n个发射信号估计列向量进行拼接，得到长度为x×z×8的估计列向量
(6c)将
(6d)以z为单位长度对vg′按从上至下的顺序进行切分，得到包括x个像素点估计子向量的像素点估计子列向量集合
(6e)将vg″′中由前到后排列的x个转置子向量转化为自上而下排列的像素估计矩阵
(6f)对u内每个像素估计矩阵ug进行奇异值分解，得到奇异值矩阵集合
(7)基于信源特性计算消除模糊度后的盲信道估计矩阵adone：
(7a)设累积速度值为ρ，奇异值矩阵序号值为
(7b)构造长度均为q的奇异值向量eg和累积量向量tg：
其中，
(7c)通过tg以及eg中所有元素的和sg计算信源特性向量
其中，
(7d)令d＝1；
(7e)判断
(7f)判断d＜ρ是否成立，若是，令ρ＝d，
(7g)令d＝d 1，执行步骤(7e)；
(7h)判断g＝2n×n！是否成立，若是，将p中的第
(7i)通过pdone的取逆结果(pdone)-1与步骤(4b)中的接收信号矩阵y计算消除模糊度后的盲信道估计矩阵adone＝y×(pdone)-1。
2.根据权利要求1所述的用于非合作mimo通信的基于信源特性的盲信道估计模糊度消除方法，其特征在于，步骤(3e)中所述的计算hs与hc的笛卡尔积，表达式为：
h′＝{(η,γ)|η∈hs∩γ∈hc}
其中，∩表示逻辑“且”。
技术总结