CN116156518A - 基于智能反射面辅助的非正交多址接入通信系统及方法 - Google Patents

基于智能反射面辅助的非正交多址接入通信系统及方法 Download PDF

Info

Publication number
CN116156518A
CN116156518A CN202310043881.4A CN202310043881A CN116156518A CN 116156518 A CN116156518 A CN 116156518A CN 202310043881 A CN202310043881 A CN 202310043881A CN 116156518 A CN116156518 A CN 116156518A
Authority
CN
China
Prior art keywords
signal
expressed
formula
channel
ris
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202310043881.4A
Other languages
English (en)
Inventor
李美玲
高文华
董增寿
程志豪
韩贺永
宁永宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Taiyuan University of Science and Technology
Original Assignee
Taiyuan University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Taiyuan University of Science and Technology filed Critical Taiyuan University of Science and Technology
Priority to CN202310043881.4A priority Critical patent/CN116156518A/zh
Publication of CN116156518A publication Critical patent/CN116156518A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/22Traffic simulation tools or models
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/309Measuring or estimating channel quality parameters
    • H04B17/336Signal-to-interference ratio [SIR] or carrier-to-interference ratio [CIR]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W72/00Local resource management
    • H04W72/04Wireless resource allocation
    • H04W72/044Wireless resource allocation based on the type of the allocated resource
    • H04W72/0473Wireless resource allocation based on the type of the allocated resource the resource being transmission power
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Quality & Reliability (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明涉及无线网络通信领域,公开了一种基于智能反射面辅助的非正交多址接入通信系统及方法。系统包括一个井下基站BS、一个井下用的小型智能反射面RIS、两个井下通信设备Dm和Dn以及两个井下工业WiFi接入点W1和W2;本发明首次考虑将基于智能反射面辅助的非正交多址接入技术应用于井下环境中,考虑实际环境中空间的限制、WiFi网络的干扰以及不同用户在不同巷道的情况,分析了井下用户中断性能,推导了中断概率的闭合表达式,并通过仿真进行了验证,证实了在综合考虑井下环境的实际情况下,本发明系统及方法的可靠性和正确性。可为实际井下通信系统运用智能反射面辅助通信方案提供参考,并对使用智能反射单元个数提供依据。

Description

基于智能反射面辅助的非正交多址接入通信系统及方法
技术领域
本发明涉及无线网络通信领域,特别是涉及一种基于智能反射面辅助的非正交多址接入通信系统及方法。
背景技术
非正交多址接入技术有着极高的频谱效率,将非正交多址接入技术(NOMA)应用于井下巷道环境可以有效满足井下智慧矿山建设对高带宽、广连接和低时延的传输需求。现有的技术中,虽然对井下无线通信采用扩频通信或者中继协作通信技术的可靠性传输性能进行了研究。然而智慧矿山的建设中对频谱带宽的要求很大,并且井下环境中空间的局限性,设备的运行以及粉尘,还有巷道的分岔情况会严重降低通信传输性能。因此,如何在复杂情况下满足井下智慧矿山的通信要求是值得深入研究的课题。
发明内容
现有的技术中尚没有关于基于智能反射面辅助非正交多址接入网络在矿井环境下的可靠传输性能研究。因此,矿井环境条件下基于智能反射面辅助的非正交多址接入网络的可靠传输,成为井下通信领域技术人员亟需解决的技术问题。本发明的目的是提供一种井下环境条件下基于智能反射面的非正交多址接入通信系统,能够基于更实际井下环境条件提升无线通信系统的可靠性。
为了达到上述目的,本发明采用了下列技术方案:
本发明提供一种基于智能反射面辅助的非正交多址接入通信系统,所述系统适用于下行链路的井下环境条件,包括:一个井下基站BS、一个井下用的小型智能反射面RIS、两个井下通信设备Dm和Dn以及两个井下工业WiFi接入点W1和W2;其中:井下基站BS部署在主巷道的侧壁上;小型智能反射面RIS配备N个RUs,部署在主巷道与支巷道岔路口的侧壁上;井下通信设备Dm位于支巷道上,Dn位于主巷道上,Dn为近用户,Dm为远用户,基站BS与Dm和Dn不在一个视距链路上,依靠智能反射面RIS的多角度反射特性使信号在BS与Dm和Dn之间传递;工业WiFi接入点W1部署在支巷道上,W2部署在主巷道上。
本发明还提供一种基于智能反射面辅助的非正交多址接入通信方法,所述方法包括以下步骤:
步骤S1、基站BS向智能反射面RIS发送叠加的混合信号;
步骤S2、智能反射面RIS向非正交多址接入用户转发信息;
步骤S3、基于非正交多址接入准则,通信设备根据用户信号功率从高到低依次进行解码;
步骤S4、根据井下巷道特征,分析瑞利衰落信道的参数设置。
进一步,所述步骤S1中基站BS向智能反射面RIS发送叠加的混合信号为:
Figure BDA0004051667660000021
/>
其中xk表示用户Dk的信息,
Figure BDA0004051667660000022
Ps为BS发射功率,αk代表Dk的功率分配系数,αn<αm,αmn=1。
进一步,所述步骤S2的具体过程为:
智能反射面向非正交多址接入用户转发信息,此时用户会收到来自WiFi接入点W1和W2的干扰信号造成信号干扰;假设所有节点均配备单天线,基站BS对第i个RU的通道系数、第i个RU对Dk的通道系数分别表示为hsi和gik,接入点Wq(q∈(1,2))对Dk的通道系数表示为hwk;则设备Dk接收到的信息可以表示为:
Figure BDA0004051667660000023
其中,ri=βiexp(jφi)代表i个反射单元的响应,j为虚数,βi和φi分别代表第i个反射元的振幅和相移反射系数,在不丧失一般性的情况下,假设βi=1,
Figure BDA0004051667660000031
代表Dk处的AWGN,
Figure BDA0004051667660000032
在本文假设/>
Figure BDA0004051667660000033
Figure BDA0004051667660000034
是接入点Wq(q∈(1,2))的发射功率,xW是干扰源发射的信息;
通过微控制器设置相应的电压,可以独立实现RIS元件的不同相移,假设RIS完全了解BS→RIS信道hsi的相位θi与RIS→Dk信道gik的相位ψi,则:
φi=-(θii) (3)
则ri可以进一步写成:
ri=exp(-j(θii)) (4)
然后将公式(4)带入到公式(2)得:
Figure BDA0004051667660000035
其中,
Figure BDA0004051667660000036
再将公式(1)带入公式(5)中得:
Figure BDA0004051667660000037
进一步,所述步骤S3的具体过程为:
SIC的基本原理是根据用户信号功率从高到低依次进行解码,在解码当前信号时,要将其他用户信号视为干扰;根据非正交多址接入准则,分配更多的功率给距离基站较远的用户Dm,因此Dn首先解调Dm的消息,由公式可得Dn检测xm时的信号干扰噪声比SINR:
Figure BDA0004051667660000038
其中,
Figure BDA0004051667660000039
表示终端设备的信噪比,/>
Figure BDA00040516676600000310
表示干扰源W所发出的信噪比,
Figure BDA00040516676600000311
是BS→RIS→Dn的等价信道;/>
在Dn接收到的信号中,Dn先减去Dm的信号,然后再解码自身信号,因此,Dn解码自身信息时的SINR为:
Figure BDA0004051667660000041
当Dm解码自身信号时,此时信号强度较弱的Dn的信号将被视为干扰噪声,则Dm解码自身信号时的SINR可表示为:
Figure BDA0004051667660000042
其中,
Figure BDA0004051667660000043
是BS→RIS→Dm的等价信道。
进一步,所述步骤S4的具体过程为:
井下环境存在着矿车,矿工,大型机电设备等,在空间受限的巷道中,障碍物占据了相当的空间,阻碍着信号的传递,因其特殊性,我们考虑任意通信链路的信道均服从瑞利衰落信道,由于井下电磁波传播的特征,不能用固定的衰落参数描述瑞利衰落特征,本发明根据井下巷道特征,分析瑞利衰落信道的参数设置方法;
瑞利信道的参数由以下方法推导而出,后面的分析过程中,先用t表示发射节点,用r表示接收节点,htr表示从发射节点到接收节点链路的信道增益;
将隧道横截面看成是一个宽为2a,高为2b的矩形,εa、εh、εv分别表示隧道内空气的、水平墙壁的和垂直墙壁的相对介电常数,σa、σh、σv则分别是隧道内空气的、水平墙壁的和垂直墙壁的电导率,μ0表示为电磁波渗透系数,f0表示为电磁波的频率,m和n表示隧道内垂直分量和水平分量的波模阶数,由于隧道内不同反射点处波模阶数不同,m和n的取值组合个数表示为P,(m,n)的任意一种模态阶表示为p阶,相应的m和n对应表示为mp和np;假设巷道内发射端到接收端之间存在M个移动物体,Er,p,{p∈(1,…,P)}表示接收端r处p阶波模的场分布,与m、n、隧道的宽a和高b有关,表达式为公式(14);Ct,p,{p∈(1,…,P)}表示发射端t处p阶波模的强度,与m、n、隧道的宽和高、以及发射端场强E0有关,表达式为公式(15);Γp=αp+jβp,{p∈(1,…,P)}表示p阶模式下电磁波的衰减系数,αp和βp分别为p阶波模的衰减系数和相移系数,可分别由f0、m、n、a、b、εa、εh、εv、σa、σh、σv、μ0这些物理量所表示,表达式为公式(16)和公式(17),
Figure BDA0004051667660000051
表示第k个物体运动时引起的额外损耗,是均值为μp,方差为/>
Figure BDA0004051667660000052
的独立的随机变量;dr为发射天线到接收天线之间的距离。故信道增益可以表示为:
Figure BDA0004051667660000053
其中,Gt和Gr分别为发射端和接收端的天线增益,Er是在接收端位置的本征向量,
Figure BDA0004051667660000054
是发射端到接收端之间的传播矩阵,Ct是在发射端位置的模式强度向量,他们分别表示为:
Er=[Er,1,Er,2,…,Er,P] (11)
Figure BDA0004051667660000055
Ct=[Ct,1,Ct,2,…,Ct,P] (13)
Figure BDA0004051667660000056
Figure BDA0004051667660000057
Figure BDA0004051667660000061
Figure BDA00040516676600000615
式中,x和y表示横截面处以横截面中心为原点建立坐标系的位置坐标,
Figure BDA0004051667660000062
和/>
Figure BDA0004051667660000063
的值由mp和np决定;具体的,当mp为偶数时,/>
Figure BDA0004051667660000064
等于0,当mp为奇数时,/>
Figure BDA00040516676600000612
等于π/2;当np为奇数时,
Figure BDA00040516676600000613
等于0,当np为偶数时,/>
Figure BDA00040516676600000614
等于π/2;k为电磁波的传播强度,与f0、ε0、σ0、μ0这些物理量有关,其表达式为公式(18),kv和kh分别为隧道水平和竖直方向上的相对电气参数,可分别由f0、εa、εh、εv、σa、σh、σv、ε0这些物理量表示,其表达式为公式(19)和公式(20):
Figure BDA0004051667660000067
Figure BDA0004051667660000068
Figure BDA0004051667660000069
为了得到随机变量|htr|2的均值,我们令
Figure BDA00040516676600000610
则|htr|2可以表示为:
Figure BDA00040516676600000611
其中,
Figure BDA00040516676600000616
X1的均值和方差可以表示为:
Figure BDA0004051667660000071
Figure BDA0004051667660000072
式(22)(23)中P(M=m)是移动物体数量M=m的概率,其值为:
Figure BDA0004051667660000073
其中v为移动物体平均速率,λ为移动物体到达的平均速度(个/秒),将(24)代入(22)可得:
Figure BDA0004051667660000074
式(23)中的最后子式可以进一步计算为:
Figure BDA0004051667660000075
将(24)和(26)代入到公式(23)得:
Figure BDA0004051667660000076
X2的均值和方差与X1的均值和方差获得过程相同,最后可得:
Figure BDA0004051667660000081
已知瑞利信道增益的平方服从指数分布,其概率密度函数可表示为:
Figure BDA0004051667660000082
其均值为
Figure BDA0004051667660000083
是瑞利衰落系数,将其与公式(28)联立可求出/>
Figure BDA0004051667660000084
表达式为:/>
Figure BDA0004051667660000085
基于同一发明构思,本申请还提供了一种上述通信系统的中断性能验证装置,包括:模型建立模块,用于确定所述井下通信系统的通信信道;第一构建模块,用于构建目标通信用户在该通信信道下的干扰信号模型和概率密度模型;第二构建模块,适于基于所述干扰信号模型和概率密度模型,构建目标通信用户的中断性能模型。
基于同一发明构思,本申请还提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现上述基于智能反射面辅助的非正交多址接入通信方法。
基于同一发明构思,本申请还提供了一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令用于使计算机执行上述基于智能反射面辅助的非正交多址接入通信方法。
与现有技术相比本发明具有以下优点:
本发明首次考虑将基于智能反射面辅助的非正交多址接入技术应用于井下环境中,考虑实际环境中空间的限制、WiFi网络的干扰以及不同用户在不同巷道的情况,分析了井下用户中断性能,推导了中断概率的闭合表达式,并通过仿真进行了验证,证实了在综合考虑井下环境的实际情况下,本发明系统及方法的可靠性和正确性。可为实际井下通信系统运用智能反射面辅助通信方案提供参考,并对使用智能反射单元个数提供依据。
本申请将6G智能反射面合理的部署在井下巷道中,可以在不增设中继基站的情况下有效提高井下环境中无线通信系统的抗干扰情况,进而提高传输性能。
附图说明
图1为矿井巷道环境下下行IRS辅助的非正交多址接入NOMA通信系统模型。
图2为NOMA用户的中断概率随发射端信噪比的变化情况图。
图3为NOMA用户的中断概率随发射端信噪比和干扰端信噪比的变化情况图。
图4为RIS-NOMA系统和NOMA系统随发射端信噪比的变化情况图。
图5 RIS-NOMA系统和RIS-OMA系统的系统吞吐量随发射端信噪比的变化情况图。
具体实施方式
下面结合本发明实施例和附图,对本发明的技术方案进行具体、详细的说明。应当指出,对于本领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干变型和改进,这些也应视为属于本发明的保护范围。
一种基于智能反射面辅助的非正交多址接入通信系统,如图1所示,包括井下基站BS,井下用的小型智能反射面RIS,两个井下通信设备Dm和Dn,构成了矿井内下行链路通信系统,W1和W2为两个井下WiFi接入端;其中:井下基站BS部署在主巷道的侧壁上;小型智能反射面RIS配备N个RUs,部署在主巷道与支巷道岔路口的侧壁上;井下通信设备Dm位于支巷道上,Dn位于主巷道上,Dn为近用户,Dm为远用户,基站BS与Dm和Dn不在一个视距链路上,依靠智能反射面RIS的多角度反射特性使信号在BS与Dm和Dn之间传递;工业WiFi接入点W1部署在支巷道上,W2部署在主巷道上。在BS经过RIS向Dm和Dn通信时,接入端W1会干扰Dm,接入端W2会干扰Dn。BS发射的信号通过智能反射面相位抵消和补充,以使接收信号干扰噪声比(SINR)最大化。h代表两用户之间的信道增益。
我们假设所有节点均配备单天线。基站BS对第i个RU的通道系数、第i个RU对Dk{k∈(m,n)}的通道系数分别表示为hsi和gik,接入点(Wq(q∈(1,2)))对Dk的通道系数表示为hwk
在所提出的系统中,在第一阶段基站BS向智能反射面RIS发送叠加的混合信号,发送的信号为:
Figure BDA0004051667660000101
其中xk表示用户Dk的信息,
Figure BDA0004051667660000102
Ps为BS发射功率,αk代表Dk的功率分配系数,αn<αm,αmn=1。
在第二阶段,智能反射面RIS向非正交多址接入用户转发信息,此时WiFi接入点W1和W2会主动向用户发送干扰信息造成信号干扰;则设备Dk接收到的信息可以表示为:
Figure BDA0004051667660000103
其中,ri=βiexp(jφi)代表i个反射单元的响应,j为虚数,βi和φi分别代表第i个反射元的振幅和相移反射系数,在不丧失一般性的情况下,假设βi=1,
Figure BDA0004051667660000111
代表Dk处的AWGN,
Figure BDA0004051667660000112
在本文假设/>
Figure BDA0004051667660000113
Figure BDA0004051667660000114
是接入点Wq(q∈(1,2))的发射功率,xW是干扰源发射的信息;
通过微控制器设置相应的电压,可以独立实现RIS元件的不同相移,假设RIS完全了解BS→RIS信道hsi的相位θi与RIS→Dk信道gik的相位ψi,则:
φi=-(θii) (3)
则ri可以进一步写成:
ri=exp(-j(θii)) (4)
然后将公式(4)带入到公式(2)得:
Figure BDA0004051667660000115
其中,
Figure BDA0004051667660000116
再将公式(1)带入公式(5)中得:
Figure BDA0004051667660000117
第三阶段,SIC的基本原理是根据用户信号功率从高到低依次进行解码,在解码当前信号时,要将其他用户信号视为干扰;根据非正交多址接入准则,分配更多的功率给距离基站较远的用户Dm,因此Dn首先解调Dm的消息,由公式可得Dn检测xm时的信号干扰噪声比SINR:
Figure BDA0004051667660000118
其中,
Figure BDA0004051667660000119
表示终端设备的信噪比,/>
Figure BDA00040516676600001110
表示干扰源W所发出的信噪比,
Figure BDA00040516676600001111
是BS→RIS→Dn的等价信道;
在Dn接收到的信号中,Dn先减去Dm的信号,然后再解码自身信号,因此,Dn解码自身信息时的SINR为:
Figure BDA0004051667660000121
当Dm解码自身信号时,此时信号强度较弱的Dn的信号将被视为干扰噪声,则Dm解码自身信号时的SINR可表示为:
Figure BDA0004051667660000122
其中,
Figure BDA0004051667660000123
是BS→RIS→Dm的等价信道。
第四阶段,井下环境存在着矿车,矿工,大型机电设备等,在空间受限的巷道中,障碍物占据了相当的空间,阻碍着信号的传递,因其特殊性,我们考虑任意通信链路的信道均服从瑞利衰落信道。由于井下电磁波传播的特征,不能用固定的衰落参数描述瑞利衰落特征,本发明根据井下巷道特征,分析瑞利衰落信道的参数设置方法;
瑞利信道的参数由以下方法推导而出。后面的分析过程中,先用t表示发射节点,用r表示接收节点,htr表示从发射节点到接收节点链路的信道增益;
将隧道横截面看成是一个宽为2a,高为2b的矩形,εa、εh、εv分别表示隧道内空气的、水平墙壁的和垂直墙壁的相对介电常数,σa、σh、σv则分别是隧道内空气的、水平墙壁的和垂直墙壁的电导率,μ0表示为电磁波渗透系数,f0表示为电磁波的频率,m和n表示隧道内垂直分量和水平分量的波模阶数,由于隧道内不同反射点处波模阶数不同,m和n的取值组合个数表示为P,(m,n)的任意一种模态阶表示为p阶,相应的m和n对应表示为mp和np;假设巷道内发射端到接收端之间存在M个移动物体,Er,p,{p∈(1,…,P)}表示接收端r处p阶波模的场分布,与m、n、隧道的宽a和高b有关,表达式为公式(14);Ct,p,{p∈(1,…,P)}表示发射端t处p阶波模的强度,与m、n、隧道的宽和高、以及发射端场强E0有关,表达式为公式(15);Γp=αp+jβp,{p∈(1,…,P)}表示p阶模式下电磁波的衰减系数,αp和βp分别为p阶波模的衰减系数和相移系数,可分别由f0、m、n、a、b、εa、εh、εv、σa、σh、σv、μ0这些物理量所表示,表达式为公式(16)和公式(17),
Figure BDA0004051667660000131
表示第k个物体运动时引起的额外损耗,是均值为μp、方差为/>
Figure BDA0004051667660000132
的独立的随机变量;dr为发射天线到接收天线之间的距离。故信道增益可以表示为:
Figure BDA0004051667660000133
其中,Gt和Gr分别为发射端和接收端的天线增益,Er是在接收端位置的本征向量,
Figure BDA0004051667660000134
是发射端到接收端之间的传播矩阵,Ct是在发射端位置的模式强度向量,他们分别表示为:
Er=[Er,1,Er,2,…,Er,P] (11)
Figure BDA0004051667660000135
Ct=[Ct,1,Ct,2,…,Ct,P] (13)
Figure BDA0004051667660000136
Figure BDA0004051667660000137
Figure BDA0004051667660000138
Figure BDA0004051667660000139
式中,x和y表示横截面处以横截面中心为原点建立坐标系的位置坐标,
Figure BDA0004051667660000141
和/>
Figure BDA0004051667660000142
的值由mp和np决定;具体的,当mp为偶数时,/>
Figure BDA0004051667660000143
等于0,当mp为奇数时,/>
Figure BDA0004051667660000144
等于π/2;当np为奇数时,/>
Figure BDA0004051667660000145
等于0,当np为偶数时,/>
Figure BDA0004051667660000146
等于π/2;k为电磁波的传播强度,与f0、ε0、σ0、μ0这些物理量有关,其表达式为公式(18),kv和kh分别为隧道水平和竖直方向上的相对电气参数,可分别由f0、εa、εh、εv、σa、σh、σv、ε0这些物理量表示,其表达式为公式(19)和公式(20):
Figure BDA0004051667660000147
Figure BDA0004051667660000148
Figure BDA0004051667660000149
为了得到随机变量|htr|2的均值,我们令
Figure BDA00040516676600001410
则|htr|2可以表示为:/>
Figure BDA00040516676600001411
其中,
Figure BDA00040516676600001412
X1的均值和方差可以表示为:
Figure BDA00040516676600001413
Figure BDA0004051667660000151
式(22)、(23)中P(M=m)是移动物体数量M=m的概率,其值为:
Figure BDA0004051667660000152
其中v为移动物体平均速率,λ为移动物体到达的平均速度(个/秒),将(24)代入(22)可得:
Figure BDA0004051667660000153
式(23)中的最后子式可以进一步计算为:
Figure BDA0004051667660000154
将(24)和(26)代入到公式(23)得:
Figure BDA0004051667660000155
X2的均值和方差与X1的均值和方差获得过程相同,最后可得:
Figure BDA0004051667660000156
已知瑞利信道增益的平方服从指数分布,其概率密度函数可表示为:
Figure BDA0004051667660000161
其均值为
Figure BDA0004051667660000162
是瑞利衰落系数,将其与公式(28)联立可求出/>
Figure BDA0004051667660000163
表达式为:
Figure BDA0004051667660000164
研究方法的可靠性,分析系统的中断性能:为了验证本申请系统及方法和可靠性、准确性,对上述方法进行理论分析和仿真验证。
(1)Dn端中断概率分析
当Dn端未成功解码信息xm和信息xn时Dn端将发生中断事件。所以Dn端中断概率可以表示为:
Figure BDA0004051667660000165
其中
Figure BDA0004051667660000166
Rk是解码xk的目标速率阈值。当/>
Figure BDA0004051667660000167
Figure BDA0004051667660000168
时,原式为:
Figure BDA0004051667660000171
已知信道增益系数hsi,gin均为瑞利衰落增益,则增益系数|hsi||gin|服从双瑞利分布,信道系数
Figure BDA0004051667660000172
可以看成是N个独立且相同的双瑞利过程之和,|An|2的CDF可以写为:/>
Figure BDA0004051667660000173
其中
Figure BDA0004051667660000174
所以:
Figure BDA0004051667660000175
其中γ(·)和Γ()分别表示下不完全Gamma函数和Gamma函数,
Figure BDA0004051667660000176
为瑞利衰落信道的衰落系数,/>
Figure BDA0004051667660000177
将公式(34)带入到公式(32)得:
Figure BDA0004051667660000178
式中
Figure BDA0004051667660000179
为W2的发射信噪比,当/>
Figure BDA00040516676600001710
取值非常大时,原式可写为:
Figure BDA00040516676600001711
Figure BDA0004051667660000181
所以原式为:
Figure BDA0004051667660000182
式中Dp(z)为抛物线圆柱函数。
(2)Dm端中断概率分析
跟分析Dn的中断概率过程一样,当Dm端解码信息xm失败时,Dm端将发生中断事件,当
Figure BDA0004051667660000183
即/>
Figure BDA0004051667660000184
时,Dm端的中断概率可以表示为:/>
Figure BDA0004051667660000185
其中
Figure BDA0004051667660000186
式中γ(·)和Γ()分别表示不完全Gamma函数和Gamma函数,/>
Figure BDA0004051667660000187
为瑞利衰落信道的衰落系数,
Figure BDA0004051667660000188
将其带入公式(38)得:
Figure BDA0004051667660000189
式中
Figure BDA00040516676600001810
为W1的发射信噪比,当/>
Figure BDA00040516676600001811
取值非常大时原式可写为:
Figure BDA00040516676600001812
Figure BDA0004051667660000191
则原式为:
Figure BDA0004051667660000192
式中Dp(z)为抛物线圆柱函数。
(3)系统验证:
图2为NOMA用户的中断概率随发射端信噪比的变化情况。从图2可以看出,在该通信系统中用户Dm的通信性能是要高于用户Dn的;从整体上看,用户的中断概率是随着发射端信噪比的增加而降低的;当发射端信噪比为-10dB时,随着RIS中参与响应元件数目的增加,用户的中断概率也随之而降低。因此可以通过增加发射端信噪比和增加RIS中参与响应元件的数目来改善系统的通信性能。
图3 NOMA用户的中断概率随发射端信噪比和干扰端信噪比的变化情况。在图3中可以明显看到在其他条件不变时,干扰端信噪比越大,用户的中断概率就越大。特殊的当发射端信噪比为-20dB且其他条件不变时,N=8时的中断概率变化情况要大于N=2时的中断概率变化情况,说明N值的增加确实可以减少干扰端的影响进而改善系统的性能。
图4 RIS-NOMA系统和NOMA系统随发射端信噪比的变化情况。从图4可以看出,在井下环境中当N=1时,RIS辅助的NOMA系统性能在低信噪比情况下要比NOMA系统好,但在高信噪比情况下是比不上的。而随着N值的增加,可以很明显的看出RIS辅助的NOMA系统要全面的优于NOMA系统。
图5 RIS-NOMA系统和RIS-OMA系统的系统吞吐量随发射端信噪比的变化情况。从图5可以看出,无论发射端信噪比增加到多大,其系统的吞吐量最终会确定为某一个定值,而N值的大小只会改变起点的大小而对最终值不会产生任何影响。还可以看出,RIS-NOMA系统要比RIS-OMA系统的吞吐量要好。
上述实验结果共同证明了本发明系统和技术方案的可靠性、正确性。
需要说明的是,上述对本申请的一些实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于上述实施例中的顺序来执行并且仍然可以实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
基于同一发明构思,与上述实施例方法相对应的,本申请还提供了一种上述通信系统的中断性能验证装置,包括:模型建立模块,用于确定所述井下通信系统的通信信道;第一构建模块,用于构建目标通信用户在该通信信道下的干扰信号模型和概率密度模型;第二构建模块,适于基于所述干扰信号模型和概率密度模型,构建目标通信用户的中断性能模型。上述针对本申请通信系统提出的中断性能验证装置,进一步保证了通信系统的有效性与可靠性,能够满足用户的服务质量需求。
基于同一发明构思,与上述实施例方法相对应的,本申请还提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现上述基于智能反射面辅助的非正交多址接入通信方法。
基于同一发明构思,与上述实施例方法相对应的,本申请还提供了一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令用于使计算机执行上述基于智能反射面辅助的非正交多址接入通信方法。
尽管上述实施例对已对本申请进行了详细的描述,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (9)

1.一种基于智能反射面辅助的非正交多址接入通信系统,其特征在于:所述系统适用于下行链路的井下环境条件,包括:一个井下基站BS、一个井下用的小型智能反射面RIS、两个井下通信设备Dm和Dn以及两个井下工业WiFi接入点W1和W2;其中:井下基站BS部署在主巷道的侧壁上;小型智能反射面RIS配备N个RUs,部署在主巷道与支巷道岔路口的侧壁上;井下通信设备Dm位于支巷道上,Dn位于主巷道上,Dn为近用户,Dm为远用户,基站BS与Dm和Dn不在一个视距链路上,依靠智能反射面RIS的多角度反射特性使信号在BS与Dm和Dn之间传递;工业WiFi接入点W1部署在支巷道上,W2部署在主巷道上。
2.一种基于智能反射面辅助的非正交多址接入通信方法,其特征在于:所述方法包括以下步骤:
步骤S1、基站BS向智能反射面RIS发送叠加的混合信号;
步骤S2、智能反射面RIS向非正交多址接入用户转发信息;
步骤S3、基于非正交多址接入准则,通信设备根据用户信号功率从高到低依次进行解码;
步骤S4、根据井下巷道特征,分析瑞利衰落信道的参数设置。
3.根据权利要求2所述的一种基于智能反射面辅助的非正交多址接入通信方法,其特征在于:所述步骤S1中基站BS向智能反射面RIS发送叠加的混合信号为:
Figure FDA0004051667650000011
其中xk表示用户Dk的信息,
Figure FDA0004051667650000012
Ps为BS发射功率,αk代表Dk的功率分配系数,αn<αm,αmn=1。
4.根据权利要求2所述的一种基于智能反射面辅助的非正交多址接入通信方法,其特征在于,所述步骤S2的具体过程为:
智能反射面RIS向非正交多址接入用户转发信息,此时用户会收到来自WiFi接入点W1和W2的干扰信号造成信号干扰;假设所有节点均配备单天线,基站BS对第i个RU的通道系数、第i个RU对Dk的通道系数分别表示为hsi和gik,接入点Wq(q∈(1,2))对Dk的通道系数表示为hwk;则设备Dk接收到的信息可以表示为:
Figure FDA0004051667650000021
其中,ri=βiexp(jφi)代表i个反射单元的响应,j为虚数,βi和φi分别代表第i个反射元的振幅和相移反射系数,在不丧失一般性的情况下,假设βi=1,
Figure FDA0004051667650000022
代表Dk处的AWGN,
Figure FDA0004051667650000023
在本文假设/>
Figure FDA0004051667650000024
是接入点Wq(q∈(1,2))的发射功率,xW是干扰源发射的信息;
通过微控制器设置相应的电压,可以独立实现智能反射面RIS元件的不同相移,假设智能反射面RIS完全了解BS→RIS信道hsi的相位θi与RIS→Dk信道gik的相位ψi,则:
φi=-(θii)(3)
则ri可以进一步写成:
ri=exp(-j(θii))(4)
然后将公式(4)带入到公式(2)得:
Figure FDA0004051667650000025
其中,
Figure FDA0004051667650000026
再将公式(1)带入公式(5)中得:
Figure FDA0004051667650000027
5.根据权利要求2所述的一种基于智能反射面辅助的非正交多址接入通信方法,其特征在于,所述步骤S3的具体过程为:
SIC的基本原理是根据用户信号功率从高到低依次进行解码,在解码当前信号时,要将其他用户信号视为干扰;根据非正交多址接入准则,分配更多的功率给距离基站较远的用户Dm,因此Dn首先解调Dm的消息,由公式可得Dn检测xm时的信号干扰噪声比SINR:
Figure FDA0004051667650000031
其中,
Figure FDA0004051667650000032
表示终端设备的信噪比,/>
Figure FDA0004051667650000033
表示干扰源W所发出的信噪比,
Figure FDA0004051667650000034
是BS→RIS→Dn的等价信道;
在Dn接收到的信号中,Dn先减去Dm的信号,然后再解码自身信号,因此,Dn解码自身信息时的SINR为:
Figure FDA0004051667650000035
当Dm解码自身信号时,此时信号强度较弱的Dn的信号将被视为干扰噪声,则Dm解码自身信号时的SINR可表示为:
Figure FDA0004051667650000036
其中,
Figure FDA0004051667650000037
是BS→RIS→Dm的等价信道。
6.根据权利要求2所述的一种基于智能反射面辅助的非正交多址接入通信方法,其特征在于,所述步骤S4的具体过程为:
井下环境存在着矿车,矿工,大型机电设备,在空间受限的巷道中,障碍物占据了相当的空间,阻碍着信号的传递,因其特殊性,我们考虑任意通信链路的信道均服从瑞利衰落信道,由于井下电磁波传播的特征,不能用固定的衰落参数描述瑞利衰落特征,本发明根据井下巷道特征,分析瑞利衰落信道的参数设置方法;
瑞利信道的参数由以下方法推导而出,后面的分析过程中,先用t表示发射节点,用r表示接收节点,htr表示从发射节点到接收节点链路的信道增益;
将隧道横截面看成是一个宽为2a,高为2b的矩形,εa、εh、εv分别表示隧道内空气的、水平墙壁的和垂直墙壁的相对介电常数,σa、σh、σv则分别是隧道内空气的、水平墙壁的和垂直墙壁的电导率,μ0表示为电磁波渗透系数,f0表示为电磁波的频率,m和n表示隧道内垂直分量和水平分量的波模阶数,由于隧道内不同反射点处波模阶数不同,m和n的取值组合个数表示为P,(m,n)的任意一种模态阶表示为p阶,相应的m和n对应表示为mp和np;假设巷道内发射端到接收端之间存在M个移动物体,Er,p,{p∈(1,…,P)}表示接收端r处p阶波模的场分布,与m、n、隧道的宽a和高b有关,表达式为公式(14);Ct,p,{p∈(1,…,P)}表示发射端t处p阶波模的强度,与m、n、隧道的宽和高、以及发射端场强E0有关,表达式为公式(15);Γp=αp+jβp,{p∈(1,…,P)}表示p阶模式下电磁波的衰减系数,αp和βp分别为p阶波模的衰减系数和相移系数,可分别由f0、m、n、a、b、εa、εh、εv、σa、σh、σv、μ0这些物理量所表示,表达式为公式(16)和公式(17),
Figure FDA0004051667650000041
表示第k个物体运动时引起的额外损耗,是均值为μp,方差为/>
Figure FDA0004051667650000042
的独立的随机变量;dr为发射天线到接收天线之间的距离,故信道增益可以表示为:
Figure FDA0004051667650000043
其中,Gt和Gr分别为发射端和接收端的天线增益,Er是在接收端位置的本征向量,
Figure FDA0004051667650000051
是发射端到接收端之间的传播矩阵,Ct是在发射端位置的模式强度向量,他们分别表示为:
Er=[Er,1,Er,2,…,Er,P] (11)
Figure FDA0004051667650000052
Figure FDA00040516676500000515
Figure FDA0004051667650000053
Figure FDA0004051667650000054
Figure FDA0004051667650000055
Figure FDA0004051667650000056
式中,x和y表示横截面处以横截面中心为原点建立坐标系的位置坐标,
Figure FDA0004051667650000057
和/>
Figure FDA0004051667650000058
的值由mp和np决定;具体的,当mp为偶数时,/>
Figure FDA0004051667650000059
等于0,当mp为奇数时,/>
Figure FDA00040516676500000510
等于π/2;当np为奇数时,/>
Figure FDA00040516676500000511
等于0,当np为偶数时,/>
Figure FDA00040516676500000512
等于π/2;k为电磁波的传播强度,与f0、ε0、σ0、μ0这些物理量有关,其表达式为公式(18),kv和kh分别为隧道水平和竖直方向上的相对电气参数,可分别由f0、εa、εh、εv、σa、σh、σv、ε0这些物理量表示,其表达式为公式(19)和公式(20):
Figure FDA00040516676500000513
Figure FDA00040516676500000514
Figure FDA0004051667650000061
为了得到随机变量htr 2的均值,我们令
Figure FDA0004051667650000062
则htr 2可以表示为:
Figure FDA0004051667650000063
其中,
Figure FDA0004051667650000064
X1的均值和方差可以表示为:
Figure FDA0004051667650000065
Figure FDA0004051667650000066
式(22)、(23)中P(M=m)是移动物体数量M=m的概率,其值为:
Figure FDA0004051667650000067
其中v为移动物体平均速率,λ为移动物体到达的平均速度,将(24)代入(22)可得:
Figure FDA0004051667650000068
式(23)中的最后子式可以进一步计算为:
Figure FDA0004051667650000071
将(24)和(26)代入到公式(23)得:
Figure FDA0004051667650000072
X2的均值和方差与X1的均值和方差获得过程相同,最后可得:
Figure FDA0004051667650000073
已知瑞利信道增益的平方服从指数分布,其概率密度函数可表示为:
Figure FDA0004051667650000074
其均值为
Figure FDA0004051667650000075
是瑞利衰落系数,将其与公式(28)联立可求出/>
Figure FDA0004051667650000076
表达式为:
Figure FDA0004051667650000077
7.一种权利要求1所述通信系统的中断性能验证装置,其特征在于,包括:
模型建立模块,用于确定所述井下通信系统的通信信道;
第一构建模块,用于构建目标通信用户在该通信信道下的干扰信号模型和概率密度模型;
第二构建模块,适于基于所述干扰信号模型和概率密度模型,构建目标通信用户的中断性能模型。
8.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于:所述处理器执行所述程序时实现如权利要求2-6任意一项所述的方法。
9.一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,其特征在于:所述计算机指令用于使计算机执行权利要求2-6任一项所述的方法。
CN202310043881.4A 2023-01-29 2023-01-29 基于智能反射面辅助的非正交多址接入通信系统及方法 Pending CN116156518A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310043881.4A CN116156518A (zh) 2023-01-29 2023-01-29 基于智能反射面辅助的非正交多址接入通信系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310043881.4A CN116156518A (zh) 2023-01-29 2023-01-29 基于智能反射面辅助的非正交多址接入通信系统及方法

Publications (1)

Publication Number Publication Date
CN116156518A true CN116156518A (zh) 2023-05-23

Family

ID=86350213

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310043881.4A Pending CN116156518A (zh) 2023-01-29 2023-01-29 基于智能反射面辅助的非正交多址接入通信系统及方法

Country Status (1)

Country Link
CN (1) CN116156518A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117439673A (zh) * 2023-12-21 2024-01-23 厦门大学 基于非正交多址的智能超表面辅助水声数据传输方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117439673A (zh) * 2023-12-21 2024-01-23 厦门大学 基于非正交多址的智能超表面辅助水声数据传输方法
CN117439673B (zh) * 2023-12-21 2024-03-12 厦门大学 基于非正交多址的智能超表面辅助水声数据传输方法

Similar Documents

Publication Publication Date Title
KR20160113869A (ko) 빔 포밍 방식을 지원하는 통신 시스템에서 풀 듀플렉스 방식을 운용하는 장치 및 방법
CN116156518A (zh) 基于智能反射面辅助的非正交多址接入通信系统及方法
Reddy et al. Spectrum cartography techniques, challenges, opportunities, and applications: A survey
CN113965881B (zh) 遮挡效应下的毫米波一体化通信与感知方法
CN220732796U (zh) 一种基于智能反射面辅助的非正交多址接入通信系统
Pitsiladis et al. Improving connectivity in indoor millimeter wave wireless networks using diversity reception
US20230412429A1 (en) Active User Detection and Channel Estimation Method and Device, Using Deep Neural Network
CN110856193B (zh) 基于空间信息网络的中继通信系统和通信方法
Hwang et al. SINR maximizing beamforming schemes for the full duplex amplify-and-forward relay channel
CN113315607B (zh) 一种基于harq辅助多ris系统的通信方法
Teng Outage probability and ergodic capacity analysis of satellite–terrestrial NOMA system with mixed RF/mmWave relaying
KR101564462B1 (ko) 분산 안테나 시스템에서 선택적 신호 처리 장치 및 방법
KR101954527B1 (ko) 협력 통신 시스템에서 중계 단말 선택 방법
US11012177B2 (en) Message and rate based user grouping in non-orthogonal multiple access (NOMA) networks
KR102110493B1 (ko) 무선 네트워크에서 밀리미터파 기반 통신 방법 및 장치들
Lee et al. Opportunities of MIM capture in IEEE 802.11 WLANs: analytic study
Sunil Performance analysis of moving multi-antenna relay cooperation with hybrid relaying scheme in cooperative wireless networks
CN111246578B (zh) 一种定向天线场景下的noma下行链路通信方法
CN113395092B (zh) 改善Massive MIMO系统检测能力的方法和系统
CN115314135B (zh) 一种无人机协作的通信感知一体化波形设计方法
CN102573069B (zh) 探测序列的分配方法和基站
CN118214503A (zh) 一种ngso星座系统动态时变信道频谱感知方法及系统
Stańczak et al. On in-network computation via wireless multiple-access channels with applications
Aprilya et al. The UAV Assisted Wireless Powered on D2D Communication Hybrid AF/DF Multi Relay Based
Xue et al. Performance of Spectrum Sharing in Hybrid Satellite Terrestrial Network with Opportunistic Relay Selection

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination