本发明属于无线通信技术领域,具体涉及一种毫米波系统动态用户混合预编码方法。
背景技术:
作为第五代移动通信系统的关键技术之一,毫米波可以在有效解决全球性频段紧张问题同时极大地提升系统传输速度和频谱效率。但是,毫米波系统中稀疏信道和窄波束使得其性能极易受系统中用户状态的影响。例如,当毫米波波束被人体阻挡,将损失约35db的波束增益。除此之外,当用户运动时,毫米波波束增益也会遭受不同程度的衰减。
针对毫米波部分连接天线阵系统中预编码问题,目前很多学者已进行了深入的研究,例如,利用模拟预编码矩阵块对角化结构和半正定松弛算法设计模拟预编码矩阵以降低混合预编码矩阵和全数字预编码矩阵的欧氏距离;通过在波束空间内穷搜索以得到最优的传输方向并将该方向对应的阵列导向矢量作为模拟预编码矩阵的相应列以简化混合预编码矩阵的设计复杂度;通过将传统模拟预编码矩阵的优化问题转化为幅度和相位分别优化的子问题以降低优化过程中的计算复杂度。
然而,传统混合预编码算法仅考虑优化静态用户的和速率性能以及简化静态用户混合预编码矩阵设计过程中的计算复杂度。当系统中存在动态用户时,传统混合预编码算法的和速率性能将极易受到系统中用户状态影响。因此上述传统算法不适用于动态用户场景,或很难维持其预期性能。
技术实现要素:
为了解决现有技术中存在的上述问题,本发明提供了一种毫米波系统动态用户混合预编码方法。本发明要解决的技术问题通过以下技术方案实现,本发明提供了一种毫米波系统动态用户混合预编码方法,用于毫米波部分连接天线阵系统,包括:
获取每个用户的信道状态信息,并将每个用户相邻时隙的所述信道状态信息存储于存储单元中;
根据存储的所述信道状态信息,得到每个用户的最优传播径;
判断所述最优传播径是否为直射径,得到每个用户的判断结果;
根据所述判断结果得到每个用户的状态;
根据所述用户的状态,设计用户的模拟预编码矩阵;
根据所述模拟预编码矩阵,设计用户的数字预编码矩阵;
根据所述模拟预编码矩阵和所述数字预编码矩阵,进行数据传输。
在本发明的一个实施例中,获取每个用户的信道状态信息,并将每个用户相邻时隙的所述信道状态信息存储于存储单元中,包括:
基站获取每个用户在波束训练时隙内扫描得到的所述信道状态信息hm,
其中,hm表示第m个用户的信道状态信息矩阵,m表示用户个数,
将每个用户的第t-1个时隙和第t个时隙的所述信道状态信息分别存储在第一存储单元和第二存储单元,其中t≥1。
在本发明的一个实施例中,根据存储的所述信道状态信息,得到每个用户的最优传播径,包括:
根据以下公式得到所述第一存储单元和所述第二存储单元中每个用户的最优传播径,
其中,g表示存储单元,
在本发明的一个实施例中,判断所述最优传播径是否为直射径,得到每个用户的判断结果,包括:
若所述最优传播径满足
其中,ε表示直射径和非直射径的门限值,0≤ε≤μas,μas表示传播波束的角度扩展的平均值。
在本发明的一个实施例中,根据所述判断结果得到每个用户的状态,包括:
若用户在所述第一存储单元和所述第二存储单元中的最优传播径均为直射径,则执行以下步骤:
若用户在所述第一存储单元和所述第二存储单元中的最优传播径的出发角满足
其中,δ表示用户运动和静止状态的判断阈值,
若用户在所述第一存储单元中的最优传播径为直射径,在所述第二存储单元中的最优传播径为非直射径,则执行以下步骤:
根据以下公式匹配用户在所述第一存储单元和所述第二存储单元中最相近的非直射径,得到
若
若用户在所述第一存储单元中的最优传播径为非直射径,在所述第二存储单元中的最优传播径为直射径,则执行以下步骤:
根据以下公式匹配用户在所述第一存储单元和所述第二存储单元中最相近的非直射径,得到
若
若用户在所述第一存储单元和所述第二存储单元中的最优传播径均为非直射径,则执行以下步骤:
根据以下公式匹配用户在所述第一存储单元和所述第二存储单元中最相近的非直射径,得到
若
在本发明的一个实施例中,根据所述用户的状态,设计用户的模拟预编码矩阵,包括:
对所述用户的状态进行优先级排序,其中,
静止>阻塞>运动>阻塞+运动;
根据所述优先级排序,设计第一个用户的模拟预编码矩阵;
设计其他用户的模拟预编码矩阵,得到所有用户的模拟预编码矩阵。
在本发明的一个实施例中,根据所述优先级排序,设计第一个用户的模拟预编码矩阵,包括:
i.令m=1,wm为全零向量,模拟预编码矩阵frf为全零矩阵,对第一个用户的状态进行判断,其中,所述模拟预编码矩阵frf为,
其中,模拟预编码矩阵frf的维度为
ii.若所述第一个用户的状态为静止或阻塞状态,则
f1=an(φs),
wm(1:n)=an(θs),
n=1,
其中,φs表示用户的传播波束的出发角,θs表示用户的传播波束的到达角,用户的传播波束为路径增益最大径对应的出发角和到达角,wm(1:n)表示第m个用户接收合并向量的第1个到第n个元素,n表示基站中被使用的发送子阵个数,
iii.若所述第一个用户的状态为运动或运动+阻塞状态,则通过基站中相邻nc个发送子阵的协作,以及用户的nc个接收子阵的协作以进行波束合并,且相邻的所述发送子阵和相邻的所述接收子阵的传输波束角度差为
若所述第一个用户运动方向的切线方向为逆时针方向,那么,
n=nc,
若所述第一个用户运动方向的切线方向为顺时针方向,那么,
n=nc;
其中,wm(n+1:2n)表示第m个用户接收合并向量的第n+1个到第2n个元素,wm[(nc-1)n+1:ncn]表示第m个用户接收合并向量的第(nc-1)n+1个到第ncn个元素;
iv.更新所述模拟预编码矩阵frf的值;
在本发明的一个实施例中,设计其他用户的模拟预编码矩阵,得到所有用户的模拟预编码矩阵,包括:
i.为第m个用户选择传播波束,确定第m个用户的传播波束的出发角和到达角,
其中,||||2表示取2的范数,min表示取最小值;
ii.对第m个用户的状态进行判断,若所述第m个用户的状态为静止或阻塞状态,则,
fn+1=an(φs),
wm(1:n)=an(θs),
n=n+1,
iii.若所述第m个用户状态为运动或运动+阻塞状态,则通过基站中相邻nc个所述发送子阵的协作,以及用户的nc个接收子阵的协作以进行波束合并,且相邻的所述发送子阵和相邻的所述接收子阵的传输波束角度差为
若所述第m个用户运动方向的切线方向为逆时针方向,那么,
n=n+nc,
若所述第m个用户运动方向的切线方向为顺时针方向,那么,
n=n+nc;
iv.更新所述模拟预编码矩阵frf的值,令m=m+1,重复步骤i-iv,迭代m-1次得到用户的所述模拟预编码矩阵frf。
在本发明的一个实施例中,根据所述模拟预编码矩阵,设计用户的数字预编码矩阵,包括:
构造第m个用户的组合干扰矩阵
其中,
对所述第m个用户的组合干扰矩阵
其中,um表示组合干扰矩阵
根据所述右奇异矩阵vm,得到所述第m个用户的数字预编码向量fbb(:,m),
fbb(:,m)=vm(:,m),
其中,fbb(:,m)表示数字预编码矩阵fbb的第m列,vm(:,m)表示所述右奇异矩阵vm的第m列;
令m=m+1,重复上述步骤,迭代m-1次得到用户的所述模拟预编码矩阵fbb。
在本发明的一个实施例中,根据所述模拟预编码矩阵和所述数字预编码矩阵,进行数据传输,包括:
基站根据所述信道状态信息、所述模拟预编码矩阵、所述数字预编码矩阵以及所述用户的接收合并向量,进行数据传输,第m个用户的接收信号ym为,
其中,nm表示第m个用户的接收噪声向量,x表示传输的信号向量。
与现有技术相比,本发明的有益效果在于:
1、本发明的毫米波系统动态用户混合预编码方法中,在不借助用户定位系统或陀螺仪的情况下,只通过对用户相邻时隙的信道状态信息进行比较判断即可得到每个用户的状态,提高了用户状态检测的效率以及降低了检测的复杂度;
2、本发明的毫米波系统动态用户混合预编码方法,通过利用相邻天线子阵的协作以及使得相邻子阵传播波束出发角和到达角间隔为
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂,以下特举较佳实施例,并配合附图,详细说明如下。
附图说明
图1是本发明实施例提供的一种系统模型图;
图2是本发明实施例提供的一种毫米波系统动态用户混合预编码方法流程图;
图3是本发明实施例提供的一种波束模式对比图;
图4是本发明实施例提供的另一种波束模式对比图;
图5是本发明实施例提供的一种用户运动时和速率随时间变化仿真图;
图6是本发明实施例提供的另一种用户运动时和速率随时间变化仿真图;
图7是本发明实施例提供的又一种用户运动时和速率随时间变化仿真图;
图8是本发明实施例提供的再一种用户运动时和速率随时间变化仿真图;
图9是本发明实施例提供的一种用户状态随机时和速率随时间变化仿真图;
图10是本发明实施例提供的另一种用户状态随机时和速率随时间变化仿真图;
图11是本发明实施例提供的又一种用户状态随机时和速率随时间变化仿真图;
图12是本发明实施例提供的再一种用户状态随机时和速率随时间变化仿真图。
具体实施方式
为了进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及具体实施方式,对依据本发明提出的一种毫米波系统动态用户混合预编码方法进行详细说明。
有关本发明的前述及其他技术内容、特点及功效,在以下配合附图的具体实施方式详细说明中即可清楚地呈现。通过具体实施方式的说明,可对本发明为达成预定目的所采取的技术手段及功效进行更加深入且具体地了解,然而所附附图仅是提供参考与说明之用,并非用来对本发明的技术方案加以限制。
实施例一
请参见图1,图1是本发明实施例提供的一种系统模型图,如图所示,本实施例使用的毫米波通信系统由毫米波部分连接天线阵基站、发送信道和用户构成。基站有
请参见图2,图2是本发明实施例提供的一种毫米波系统动态用户混合预编码方法流程图,如图所示,本实施例的方法包括:
s1:获取每个用户的信道状态信息,并将每个用户相邻时隙的所述信道状态信息存储于存储单元中;
具体地,包括:
s11:基站获取每个用户在波束训练时隙内扫描得到的所述信道状态信息hm,
其中,hm表示第m个用户的信道状态信息矩阵,m表示用户个数,
在本实施例中,所述信道状态信息包括每个用户与均匀线性天线阵基站之间的出发角、到达角和路径增益,
其中,
s12:将每个用户的第t-1个时隙和第t个时隙的所述信道状态信息分别存储在第一存储单元和第二存储单元,其中t≥1。
在本实施例中,第一存储单元和第二存储单元只存储每个用户的相邻时隙的信道状态信息,在需要存储下一个时隙的信道状态信息时,之前时隙的信道状态信息被覆盖或被删除。
s2:根据存储的所述信道状态信息,得到每个用户的最优传播径;
具体地,包括:
根据以下公式得到所述第一存储单元和所述第二存储单元中每个用户的最优传播径,
其中,g表示存储单元,
s3:判断所述最优传播径是否为直射径,得到每个用户的判断结果;
具体地,包括:
若所述最优传播径满足
若所述最优传播径满足
其中,ε表示直射径和非直射径的门限值,0≤ε≤μas,μas表示传播波束的角度扩展的平均值,在本实施例中,所述门限值根据系统的角度扩展值和系统需求设定。如上式所示,若用户在第g个存储单元中所述最优传播径的出发角和接收角之和的绝对值小于ε,则所述最优传播径为直射径,否则为非直射径。
s4:根据所述判断结果得到每个用户的状态;
在本实施例中,用户的状态包括静止、运动、阻塞和阻塞+运动四种状态,其中,静止状态表示用户没有阻塞且没有运动,运动状态表示用户处于移动的状态,阻塞状态表示毫米波波束被障碍物阻挡,阻塞+运动状态表示用户同时存在运动和阻塞的情况。
具体地,包括:
s41:若用户在所述第一存储单元和所述第二存储单元中的最优传播径均为直射径,则执行以下步骤:
若用户在所述第一存储单元和所述第二存储单元中的最优传播径的出发角满足
若用户在所述第一存储单元和所述第二存储单元中的最优传播径的出发角满足
其中,δ表示用户运动和静止状态的判断阈值,其取值根据系统要求而定,在本实施例中,
s42:若用户在所述第一存储单元中的最优传播径为直射径,在所述第二存储单元中的最优传播径为非直射径,则执行以下步骤:
根据以下公式匹配用户在所述第一存储单元和所述第二存储单元中最相近的非直射径,得到
若
s43:若用户在所述第一存储单元中的最优传播径为非直射径,在所述第二存储单元中的最优传播径为直射径,则执行以下步骤:
根据以下公式匹配用户在所述第一存储单元和所述第二存储单元中最相近的非直射径,得到
若
s44:若用户在所述第一存储单元和所述第二存储单元中的最优传播径均为非直射径,则执行以下步骤:
根据以下公式匹配用户在所述第一存储单元和所述第二存储单元中最相近的非直射径,得到
若
本实施例的毫米波系统动态用户混合预编码方法中,在不借助用户定位系统或陀螺仪的情况下,只通过对用户相邻时隙的信道状态信息进行比较判断即可得到每个用户的状态,提高了用户状态检测的效率以及降低了检测的复杂度。
s5:根据所述用户的状态,设计用户的模拟预编码矩阵;
具体地,包括:
s51:对所述用户的状态进行优先级排序,其中,
静止>阻塞>运动>阻塞+运动;
在本实施例中,静止状态用户的处理优先级最高,首先进行数据传送处理,阻塞+运动状态的用户处理优先级最低,最后进行数据传送处理则,若同一状态中存在多个用户,多个用户之间的处理优先级不分先后,可以选择其中任一个用户开始进行数据传送处理。
s52:根据所述优先级排序,设计第一个用户的模拟预编码矩阵
具体地,包括以下步骤:
i.令m=1,wm为全零向量,模拟预编码矩阵frf为全零矩阵,对第一个用户的状态进行判断,其中,所述模拟预编码矩阵frf为,
其中,模拟预编码矩阵frf的维度为
ii.若所述第一个用户的状态为静止或阻塞状态,则
f1=an(φs),
wm(1:n)=an(θs),
n=1,
其中,φs表示用户的传播波束的出发角,θs表示用户的传播波束的到达角,用户的传播波束为路径增益最大径对应的出发角和到达角,wm(1:n)表示第m个用户接收合并向量的第1个到第n个元素,n表示基站中被使用的发送子阵个数,
iii.若所述第一个用户的状态为运动或运动+阻塞状态,则通过基站中相邻nc个发送子阵的协作,以及用户的nc个接收子阵的协作以进行波束合并,用以增加传输波束的波束宽度,而且相邻的所述发送子阵和相邻的所述接收子阵的传输波束角度差为
若所述第一个用户运动方向的切线方向为逆时针方向,那么,
n=nc,
若所述第一个用户运动方向的切线方向为顺时针方向,那么,
n=nc;
其中,wm(n+1:2n)表示第m个用户接收合并向量的第n+1个到第2n个元素,wm[(nc-1)n+1:ncn]表示第m个用户接收合并向量的第(nc-1)n+1个到第ncn个元素;
在本实施例中,协作的相邻子阵nc的个数根据系统的需求而设定,基站中被使用的发送子阵的个数n在基站给用户传输数据的过程中累加计算,基站共有
iv.更新所述模拟预编码矩阵frf的值;
具体地,将步骤ii或步骤iii中得到第一个用户的f1或f1到fnc的值更新至所述模拟预编码矩阵frf中。
s53:设计其他用户的模拟预编码矩阵,得到所有用户的模拟预编码矩阵。
具体地,包括以下步骤:
i.为第m个用户选择传播波束,确定第m个用户的传播波束的出发角和到达角,
其中,||||2表示取2的范数,min表示取最小值;
ii.对第m个用户的状态进行判断,若所述第m个用户的状态为静止或阻塞状态,则,
fn+1=an(φs),
wm(1:n)=an(θs),
n=n+1,
iii.若所述第m个用户状态为运动或运动+阻塞状态,则通过基站中相邻nc个所述发送子阵的协作,以及用户的nc个接收子阵的协作以进行波束合并,且相邻的所述发送子阵和相邻的所述接收子阵的传输波束角度差为
若所述第m个用户运动方向的切线方向为逆时针方向,那么,
n=n+nc,
若所述第m个用户运动方向的切线方向为顺时针方向,那么,
n=n+nc;
iv.更新所述模拟预编码矩阵frf的值,令m=m+1,重复步骤i-iv,迭代m-1次得到所有用户的所述模拟预编码矩阵frf。
s6:根据所述模拟预编码矩阵,设计用户的数字预编码矩阵;
具体地,包括:
s61:构造第m个用户的组合干扰矩阵
其中,
s62:对所述第m个用户的组合干扰矩阵
其中,um表示组合干扰矩阵
s63:根据所述右奇异矩阵vm,得到所述第m个用户的数字预编码向量fbb(:,m),
fbb(:,m)=vm(:,m),
其中,fbb(:,m)表示数字预编码矩阵fbb的第m列,vm(:,m)表示所述右奇异矩阵vm的第m列;
s64:令m=m+1,重复s61-s64,迭代m-1次得到用户的所述模拟预编码矩阵fbb。
s7:根据所述模拟预编码矩阵和所述数字预编码矩阵,进行数据传输。
具体地,包括:
基站根据所述信道状态信息、所述模拟预编码矩阵、所述数字预编码矩阵以及所述用户的接收合并向量,进行数据传输,第m个用户的接收信号ym为,
其中,nm表示第m个用户的接收噪声向量,x表示传输的信号向量。
本实施例的毫米波系统动态用户混合预编码方法,通过利用相邻天线子阵的协作以及使得相邻子阵传播波束出发角和到达角间隔为
实施例二
本实施例是对实施例一中的毫米波系统动态用户混合预编码方法的仿真实验。本实施例的仿真条件为:基站有15个发送子阵,每个发送子阵有20根天线,系统中有5个用户,每个用户有3个接收子阵,每个接收子阵有20根天线,莱斯因子k为7db,每个信号的发送功率为单位功率。
请参见图3和图4,其中,图3为采用本发明的方法且相邻协作子阵为nc=2以及未采用本发明方法的波束模式图的对比图;图4为采用本发明的方法相邻协作子阵为nc=2和nc=3的波束模式图对比图。如图3所示,采用nc=2的本发明方法的波束远远宽与未采用本发明方法的波束,如图4所示采用nc=3的本发明方法比采用nc=2的的本发明方法波束更宽且顶端更平稳,通过图3和图4对比,可以看出本发明方法在相邻子阵波束角度差为
请参见图5-图8,是本发明实施例提供的用户在不同运动状态下采用本发明方法与不采用本发明方法的和速率随时间变化仿真图。其中,系统中5个用户都处于运动状态,且每个用户仅有直射径,系统中用户离基站的距离分别为5米、15米和35米,图5中用户的移动速度为2米/秒,图6中用户的移动速度为5米/秒,图7中用户的移动速度为15米/秒,图8中用户的移动速度为35米/秒。在图5中,当用户的移动速度为2米/秒,且离基站15米时,nc=3的波束合并算法(3个子阵的波束合并)可维持和速率性能在1000毫秒内保持不变,nc=2的波束合并算法(2个子阵的波束合并)在1000毫秒后和速率下降到13比特/秒/赫兹,而未采用波束合并算法的和速率性能在750毫秒后下降到5比特/秒/赫兹。如图6所示,利用nc=3的本发明方法,在距离15米、速度5米/秒的情况下,可以维持和速率在500毫秒内保持不变。未采用本发明方法的,在距离15米、速度5米/秒的情况下,和速率在300毫秒后衰减到5比特/秒/赫兹。如图7所示,当用户距离为35米,速度为15米/秒时,nc=2时本发明方法的和速率在200毫秒后衰减到5比特/秒/赫兹,nc=3时本发明方法的和速率在300毫秒之后衰减到5比特/秒/赫兹,而未采用本发明方法的和速率在100毫秒之后衰减到5比特/秒/赫兹。如图8所示,当用户速度为35米/秒时,采用nc=2和nc=3的本发明方法在距离为5米、15米和35米时均比未采用本发明方法的和速率性能好,特别地,当距离为35米,时间为500毫秒时,nc=2的本发明方法的和速率为4.5比特/秒/赫兹,nc=3的本发明方法的和速率为15比特/秒/赫兹和,而未采用本发明方法的和速率为1.5比特/秒/赫兹。从图5、图6、图7和图8可以看出,与未采用本发明方法的曲线相比,本发明方法可以有效地维持和速率性能,并延缓其和速率的下降速度。
请参见图9-图12,是本发明实施例提供的用户在随机状态下采用本发明方法与不采用本发明方法的和速率随时间变化仿真图。其中,系统中所有用户状态随机,其状态包括静止、运动、阻塞和运动+阻塞,每个用户有5条传播径,用户离基站的距离分别为5米、15米和35米,图9中用户的移动速度为2米/秒,图10中用户的移动速度为5米/秒,图11中用户的移动速度为15米/秒,图12中用户的移动速度为35米/秒。第i个传播径从发射端到接收端的传播时延服从指数延迟分布,其公式为,
τ′i=-γτdsln(xi),
其中,γτ=3表示延迟缩放系数,ds表示延迟扩展,lgds的均值为-0.28log10(1+fc)-7.173,fc=28ghz为载波频率。然后,通过减去最小时延以归一化时延,并将时延按降序排列可得到第i径的时延,
τi=sort{τ′i-min(τ′i)},
其中,sort表示降序排列函数,min表示取最小值。
从图9、图10、图11和图12可以看出,本发明的方法可以提高动态用户的和速率性能,相比于nc=2时的本发明方法,nc=3时的本发明方法能将和速率性能保持更长时间。如图9所示,当用户移动速度为2米/秒,且距基站15米时,nc=2和nc=3的自适应混合预编码算法可在1000毫秒内将系统和速率维持在20比特/秒/赫兹之上,而未采用自适应混合预编码算法的和速率在1000毫秒后下降到16比特/秒/赫兹。如图10所示,nc=3时本发明方法可将距离为5米,速度为5米/秒的用户的和速率性能维持500毫秒不变,但是,未采用本发明方法的和速率迅速降低。如图11所示,当距离为15米时,未采用本发明方法、采用nc=2时的本发明方法以及采用nc=3时的本发明的方法,其和速率分别在200毫秒、325毫秒和500毫秒后衰减到14比特/秒/赫兹,可以看出本发明方法可以有效地延缓动态用户的和速率性能衰减,提升系统稳定性。如图12所示,当距离为15米时,未采用本发明方法、采用nc=2时的本发明方法以及采用nc=3时的本发明的方法,其和速率分别在200毫秒、300毫秒和500毫秒后衰减到14比特/秒/赫兹,可以看出本发明方法亦可以有效地缓解当用户速度为35米/秒,距离为15米时的和速率性能地衰落。因此本发明方法可以有效地延缓动态用户的和速率性能,提升系统稳定性。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。