CN111682899A - 海洋无线光通信系统的预编码和均衡矩阵设计方法 - Google Patents
海洋无线光通信系统的预编码和均衡矩阵设计方法 Download PDFInfo
- Publication number
- CN111682899A CN111682899A CN202010477672.7A CN202010477672A CN111682899A CN 111682899 A CN111682899 A CN 111682899A CN 202010477672 A CN202010477672 A CN 202010477672A CN 111682899 A CN111682899 A CN 111682899A
- Authority
- CN
- China
- Prior art keywords
- matrix
- precoding
- equalization
- turbulence
- mse
- 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.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 132
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000004891 communication Methods 0.000 title claims abstract description 55
- 230000003287 optical effect Effects 0.000 title claims abstract description 35
- 238000013461 design Methods 0.000 title abstract description 28
- 230000014509 gene expression Effects 0.000 claims abstract description 30
- 230000000694 effects Effects 0.000 claims abstract description 25
- 238000005457 optimization Methods 0.000 claims description 30
- 239000013598 vector Substances 0.000 claims description 29
- 230000005540 biological transmission Effects 0.000 claims description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 23
- 230000003068 static effect Effects 0.000 claims description 16
- 238000005562 fading Methods 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000009795 derivation Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 8
- 238000010521 absorption reaction Methods 0.000 description 7
- 239000000047 product Substances 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000013307 optical fiber Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 241001168730 Simo Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/11—Arrangements specific to free-space transmission, i.e. transmission through air or vacuum
- H04B10/114—Indoor or close-range type systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/50—Transmitters
- H04B10/501—Structural aspects
- H04B10/503—Laser transmitters
- H04B10/505—Laser transmitters using external modulation
- H04B10/5055—Laser transmitters using external modulation using a pre-coder
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/50—Transmitters
- H04B10/516—Details of coding or modulation
- H04B10/524—Pulse modulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/50—Transmitters
- H04B10/516—Details of coding or modulation
- H04B10/54—Intensity modulation
- H04B10/541—Digital intensity or amplitude modulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/60—Receivers
- H04B10/66—Non-coherent receivers, e.g. using direct detection
- H04B10/69—Electrical arrangements in the receiver
- H04B10/697—Arrangements for reducing noise and distortion
- H04B10/6971—Arrangements for reducing noise and distortion using equalisation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Optics & Photonics (AREA)
- Optical Communication System (AREA)
Abstract
本发明提供的海洋无线光通信系统的预编码和均衡矩阵设计方法,针对采用零均值脉冲幅度调制的MIMO UWOC通信系统,根据湍流效应的统计特性对水下信道矩阵进行建模,推导了在湍流环境下,系统等效接收符号与发送符号之间的MSE表达式,并基于该表达式,提出了一种基于TS‑MMSE的预编码和均衡矩阵设计方法,有效提高MIMO UWOC系统在湍流环境下的通信可靠性。
Description
技术领域
本发明面向水下无线光通信技术领域,更具体的,涉及一种海洋无线光通信系统的预编码和均衡矩阵设计方法。
背景技术
近年来,水下无线光通信(Underwater Wireless Optical Communication,UWOC)逐渐成为了短程水下无线通信的主流方案,受到了越来越多的学者关注。UWOC系统采用大功率激光二极管(Laser Diode,LD)或发光二极管(Light-Emitting Diode,LED)作为信号源,且一般采用在水下衰减较低的蓝绿波段光作为信号载体来传输信息。相比于传统的水下无线通信技术,即水声通信和射频通信,UWOC具有更高的带宽、更低的时延以及更好的安全性等优点[1]Z.Zeng,S.Fu,H.Zhang,Y.Dong and J.Cheng,A Survey of UnderwaterOptical Wireless Communications[J],IEEE Communications Surveys&Tutorials,vol.19,no.1,pp.204-238,Firstquarter 2017.[2]M.V.Jamali,J.A.Salehi,andF.Akhoundi,Performance Studies of Underwater Wireless Optical CommunicationSystems With Spatial Diversity:MIMO Scheme[J],IEEE Transactions onCommunications,vol.65,no.3,pp.1176-1192,March 2017.因此已经被应用在各种水下通信业务中,如实时视频传输、高吞吐传感器网络以及高速短距水下通信等[2]。
然而,水下光无线信道固有特性将很大程度上影响UWOC的系统性能。这些固有特性主要是吸收、散射以及湍流效应,分别会造成信号衰减、符号间干扰(Inter SymbolInterference,ISI)以及接收光强度的随机抖动,从而降低系统的有效数据传输率[3]J.G.Proakis,Digital Communications[M],McGraw-Hill,New York,1995.[4]T.S.Rappaport,Wireless communications:principles and practice[M],PrenticeHall,New Jersey,2nd ed.,2002。近年来,研究者们对考虑了以上水下光无线信道特性的UWOC信道模型进行了全面的研究。其中,使用蒙特卡洛仿真[5]J.W.C.Cox,Simulation,Modeling and Design of Underwater Optical Communication Systems[D],Ph.D.dissertation,North Carolina State University,Raleigh,NC,USA,2012.以及随机信道模型[6]H.Zhang and Y.Dong,General Stochastic Channel Model andPerformance Evaluation for Underwater Wireless Optical Links[J],IEEETransactions on Wireless Communications,vol.15,no.2,pp.1162-1173,Feb.2016.[7]H.Zhang,J.Cheng,and Z.Wang,On Integrated Stochastic Channel Model forUnderwater Optical Wireless Communications[C],in IEEE InternationalConference on Communications(ICC),Kansas City,MO,USA,MAY,2018,pp.1-6.可以获得在考虑吸收和散射效应时,静态水体(无湍流时)的信道状态信息(Channel StateInformation,CSI),如光强度衰减以及信道冲击响应。另一方面,由湍流效应造成的衰落可以被建模为服从对数正态分布的随机乘性因子[2][8]M.V.Jamali,P.Nabavi,andJ.A.Salehi,MIMO Underwater Visible Light Communications:Comprehensive ChannelStudy,Performance Analysis,and Multiple-Symbol Detection[J],IEEE Transactionson Vehicular Technology,vol.67,no.9,pp.8823-8237,Sept.2018.[9]H.Zhang,J.Cheng,Z.Wang,and Y.Dong,On the Capacity of Buoy-Based MIMO Systems forUnderwater Optical Wireless Links with Turbulence[C],in IEEE InternationalConference on Communications(ICC),Kansas City,MO,USA,May,2018,pp.1-6.[10]Y.Dong and J.Liu,On BER performance of Underwater Wireless Optical MISO LinksUnder Weak Turbulence[C],In OCEANS,Shanghai,China,Apr.2016,pp.1-4。
另一方面,已经被应用于室内可见光通信(Visible Light Communication,VLC)的光MIMO技术由于利用了空间分集特性,可以成倍地提高无线信道容量与频谱效率。同时为提高室内MIMO VLC系统的有效性及可靠性,针对室内MIMO VLC系统的预编码器和均衡器设计也在近年来得到了广泛的研究。在文献[11]K.Park,Y.Ko,and M.Alouini,On thePower and Offset Allocation for Rate Adaptation of Spatial Multiplexing inOptical Wireless MIMO Channels[J],IEEE Transactions on Communications,vol.4,no.61,pp.1535-1543,Apr.2013.中,作者提出了一种自适应MIMO预编码和均衡方法,以最大化无调光控制时的系统传输速率。而在考虑具备调光功能以及光信号非负约束的条件下,一种采用零均值调制的MIMO预编码和均衡方法在文献[12]K.Ying,H.Qian,R.J.Baxley,and S.Yao,Joint Optimization of Precoder and Equalizer in MIMO VLCSystems[J],IEEE Journal on Selected Areas in Communications,vol.33,no.9,pp.1949-1958,Sept.2015.中被提出。在此基础上,文献[13]K.Ying,H.Qian,R.J.Baxley,and G.T.Zhou,MIMO Transceiver Design in Dynamic-Range-Limited VLC Systems[J],IEEE Photonics Technology Letters,vol.28,no.22,pp.2593-2596,Nov.2016.进一步将LED的动态范围考虑进优化问题中,而文献[14]B.Li,R.Zhang,W.Xu,C.Zhao,and L.Hanzo,Joint Dimming Control and Transceiver Design for MIMO-Aided Visible LightCommunication[J],IEEE Communications Letters,vol.20,no.11,pp.2193-2196,Nov.2016.则将预编码器和均衡器的设计扩展至采用非零均值调制的MIMO系统中。在文献[15]Q.Gao,C.Cong,and Z.Xu,Joint Transceiver and Offset Design for VisibleLight Communications With Input-Dependent Shot Noise[J],IEEE Transactions onWireless Communications,vol.16,no.5,pp.2736-2747,May 2017.中,作者则考虑了与输入信号相关的散粒噪声,提出了MIMO多颜色VLC系统的联合预编码与均衡矩阵以及直流(Direct Current,DC)偏置设计方法,其中DC偏置被联合优化以最小化发送与恢复数据之间的均方误差(Mean Squared Error,MSE)。
一般室内VLC MIMO信道状态在光源和光电检测器(Photo Detector,PD)的位置相对固定的情况下基本保持稳定。而与室内VLC MIMO信道不同,即使收发机的位置相对固定,UWOC在湍流环境下的MIMO信道还是表现出较强的随机和时变性。这种因湍流造成的信道的随机与时变性将会增大对信道估计结果可靠性的要求,从而降低MIMO系统的性能[1]。虽然已有文献[8][9]指出了利用空间分集的MIMO通信方案可以有效对抗湍流效应,但在这些方案中,所有发射机都发送相同的信号,并未利用到多重信道的特性来提升系统传输速率。此外,在考虑水下湍流衰落下的MIMO UWOC系统的预编码以及均衡矩阵设计还尚未研究。
发明内容
本发明为克服现有UWOC系统在对抗湍流效应过程中,未利用多重信道的特性对系统传输速率进行提升,存在可靠性不高的技术缺陷,提供一种海洋无线光通信系统的预编码和均衡矩阵设计方法。
为解决上述技术问题,本发明的技术方案如下:
海洋无线光通信系统的预编码和均衡矩阵设计方法,包括以下步骤:
S1:构建具有Nr×Nt条通道的MIMO UWOC系统、预编码矩阵F和均衡矩阵G;
S2:将MIMO UWOC系统的发送信号采用直流偏置向量进行处理,得到对预编码矩阵的功率约束;
S3:MIMO UWOC系统接收端接收到发送信号后,经过均衡矩阵和直流去除操作后,得到等效接收符号表达式;
S4:基于静态信道信息以及湍流效应的统计特性,获取等效接收符号与发送信号之间的MSE表达式;
S5:构建在步骤S2所得功率约束条件下最小化MSE的优化问题,其中优化变量为F和G;
S6:推导在分别给定F和G时,另一个变量在优化问题的最优解;
S7:基于步骤S6的推导,采用交替优化算法求解所构建的优化问题,对F和G进行联合优化,得到最小化MSE的最优预编码矩阵和均衡矩阵。
其中,所述步骤S1具体为:
构建一个Nr×Nt的MIMO UWOC系统由Nt个激光发射器和Nr个PD组成,发送数据源向量为s=[s1,…,sK]T,其中K≤min{Nt,Nr}表示数据总流数,[·]T表示向量转置操作符;向量s中的每个元素都取自独立的零均值调制星座符号,且满足-uk≤sk≤uk,其中uk>0是每个符号的上界;同时,构建预编码矩阵和均衡矩阵,将预编码矩阵表示为将均衡矩阵表示为
其中,所述步骤S2具体为:
不失一般性,假设线性的电光转换过程以及单位发送增益,将发送信号x表示为x=Fs+p,其中表示直流偏置向量,用来保证传输信号的非负性;在给定直流偏置向量p的条件下,根据发送信号x的非负约束,得到对预编码矩阵的功率约束,表示为abs(F)u≤p,其中u=[u1`,...,uK]T,abs(·)表示对矩阵各元素取绝对值运算。
其中,所述步骤S3具体为:
在接收端,光信号经PD通过光电转换输出为电信号,其中接收信号的向量表达式为y=Hx+n=H(Fs+p)+n,其中表示水下信道矩阵,表示零均值的加性高斯白噪声(Additive White Gaussian Noise,AWGN)向量,其协方差矩阵为Rn;
当检测器之间的横向距离大于1厘米时,所有子信道之间湍流衰落满足相互独立,这在实际应用中是可以满足的;因此,假设矩阵中的元素是相互独立的随机变量;而在湍流效应建模中,湍流衰落的强度由闪烁指数来评估,其计算表达式为:
其中Ii,j表示从第j个发射机到第i个PD的瞬时光强,Ε(·)表示统计期望;此外,弱海洋湍流衰落系数一般由对数正态分布进行描述,其概率密度函数(ProbabilityDensity Function,PDF)为:
其中和分别是服从高斯分布的对数幅度的均值和方差;此外,对水下湍流效应建模时,湍流衰落系数的期望通常会进行归一化来避免平均功率的变化,从而有成立,因此得出和除此之外,得到闪烁指数与对数幅度的方差之间的关系:
其中,在所述步骤S4中,所述等效接收符号与发送信号之间的MSE表达式具体为:
和
其中,在所述步骤S5中,所述优化问题具体表示为:
s.t.abs(F)u≤p
其中,优化变量为F和G。
其中,所述步骤S6具体为:
当给定F时最优的G,记为G*,由下式获得:
而给定G时最优的F通过求解以下线性约束二次规划凸问题获得:
其中,所述步骤S6的具体推导过程为:
1)给定F更新G
当给定预编码矩阵F时,最优的均衡矩阵G*通过求解下式获得:
因此可得
式中,两个期望项分别表示为:
和
其中X[i,j]表示矩阵X的第(i,j)个元素,Diag(a)是关于向量a的对角矩阵;
2)给定G更新F
类似地,MSE的表达式第三项改写为:
改写后的式子是关于vec(FH)的二次型;另外,易证MSE的表达式的其他项都是关于vec(FH)的线性项或常数项,因为有:
和
其中利用恒等式Tr(PQ)=vec(PH)Hvec(Q)和上式中的矩阵都为实数矩阵的事实;另一方面,同时对F的功率约束也是线性约束,因此当给定均衡矩阵G时,最优F为一个典型的线性约束二次规划凸问题的解,该凸问题用公式表示为:
s.t.abs(F)u≤p
其中,所述步骤S7具体为:交替优化算法会初始化满足优化问题约束的F和G,接着对两个变量进行迭代更新,其中每次更新时算法将其中一个变量进行固定,另一个变量进行更新,以使得MSE逐渐减小至收敛,最终获得最小化MSE的最优预编码矩阵和均衡矩阵。
上述方案中,本发明在现有的传统室内MIMO VLC预编码和均衡矩阵设计的基础上,针对湍流环境下采用零均值调制的MIMO UWOC系统,提出了一种可有效抑制水下湍流效应的预编码和均衡矩阵设计方法,即基于湍流抑制的最小均方误差准则(TurbulenceSuppressed Minimum Mean Squared Error,TS-MMSE)预编码和均衡矩阵设计方法,可有效提高MIMO UWOC系统在湍流环境下的通信可靠性。
与现有技术相比,本发明技术方案的有益效果是:
本发明提供的海洋无线光通信系统的预编码和均衡矩阵设计方法,针对采用零均值脉冲幅度调制(Pulse Amplitude Modulation,PAM)的MIMO UWOC通信系统,根据湍流效应的统计特性对水下信道矩阵进行建模,推导了在湍流环境下,系统等效接收符号与发送符号之间的MSE表达式,并基于该表达式,提出了一种基于TS-MMSE的预编码和均衡矩阵设计方法,有效提高MIMO UWOC系统在湍流环境下的通信可靠性。
附图说明
图1为本发明所述方法流程图;
图2为采用M进制PAM的MIMO UWOC系统框图;
图3为采用TS-MMSE、MMSE、ZF和SVD四种预编码和均衡矩阵设计方法时的系统BER性能对比图;
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示,海洋无线光通信系统的预编码和均衡矩阵设计方法,包括以下步骤:
S1:构建具有Nr×Nt条通道的MIMO UWOC系统、预编码矩阵F和均衡矩阵G;
S2:将MIMO UWOC系统的发送信号采用直流偏置向量进行处理,得到对预编码矩阵的功率约束;
S3:MIMO UWOC系统接收端接收到发送信号后,经过均衡矩阵和直流去除操作后,得到等效接收符号表达式;
S4:基于静态信道信息以及湍流效应的统计特性,获取等效接收符号与发送信号之间的MSE表达式;
S5:构建在步骤S2所得功率约束条件下最小化MSE的优化问题,其中优化变量为F和G;
S6:推导在分别给定F和G时,另一个变量在优化问题的最优解;
S7:基于步骤S6的推导,采用交替优化算法求解所构建的优化问题,对F和G进行联合优化,得到最小化MSE的最优预编码矩阵和均衡矩阵。
在具体实施过程中,本发明提供的海洋无线光通信系统的预编码和均衡矩阵设计方法,针对采用零均值PAM的MIMO UWOC通信系统,根据湍流效应的统计特性对水下信道矩阵进行建模,推导了在湍流环境下,系统等效接收符号与发送符号之间的MSE表达式,并基于该表达式,提出了一种基于TS-MMSE的预编码和均衡矩阵设计方法,有效提高MIMO UWOC系统在湍流环境下的通信可靠性。
更具体的,本发明所提出的基于TS-MMSE的预编码和均衡矩阵设计方法将应用于一个由Nt个激光发射器和Nr个PD组成,具有Nr×Nt条通道的MIMO UWOC系统,如图2所示。
令s=[s1,…,sK]T表示实数数据源向量,其中K≤min{Nt,Nr}表示数据总流数,[·]T表示向量转置操作符。向量s中的每个元素都取自独立的零均值M进制PAM(M-ary PAM,M-PAM)符号,即有Ε(sk)=0,k=1,…,K,其中Ε(·)表示统计期望。同时每个子数据通道上的符号sk满足
-uk≤sk≤uk (1)
其中,uk>0是PAM符号的上界。假设发送符号功率归一化时,uk满足
数据源向量s产生后,将先经过一个预编码矩阵对其进行预编码。根据LD光源的传输特性,预编码后的信号在被发送到信道之前,需加上DC偏置来保证传输信号的非负性。不失一般性,假设线性的电光转换过程以及单位发送增益[13],发送信号x可以被表示为
x=Fs+p (3)
根据式(3),每个发射器上的平均光功率可计算为
或等效地表示为
abs(F)u≤p (6)
其中记u=[u1`,…,uK]T,abs(·)表示对矩阵各元素取绝对值运算。
在接收端,光信号经PD通过光电转换输出为电信号。接收信号向量可表示为
y=Hx+n=H(Fs+p)+n (7)
考虑到由湍流引起的乘性衰落,水下信道增益矩阵可以被分解为
根据文献[16]W.Liu,Z.Xu,and L.Yang,SIMO Detection for UnderwaterOptical Wireless Communication under Turbulence,Photonics Research,vol.3,no.3,pp.48-53,Jun.2015.当检测器之间的横向距离大于1厘米时,所有子信道之间湍流衰落满足相互独立,这在实际应用中通常是可以满足的。因此,本发明假设矩阵中的元素为相互独立的随机变量。
湍流衰落的强度可由闪烁指数来评估,其计算表达式为[17]L.C.Andrews,R.L.Phillips,and C.Y.Hopen,Laser Beam Scintillation With Applications[M],Bellingham,WA,USA:SPIE,2001.
其中和分别是服从高斯分布的对数幅度的均值和方差。此外,对水下湍流效应建模时,湍流衰落系数的期望通常会进行归一化来避免平均功率的变化[8],从而有成立。据此可以推导得[18]S.M.Navidpour,M.Uysal,and M.Kavehrad,BER Performance of Free-Space Optical Transmission withSpatial Diversity[J],IEEE Transactions on Wireless Communications,vol.6,no.8,pp.2813–2819,Aug.2007.
和
除此之外,还能得到闪烁指数与对数幅度的方差之间的关系[17]:
获得接收信号y后,接收端将采用基于均衡检测的方法来估计发送数据,该方法比最大似然检测的复杂度小很多[13]。具体而言,接收信号将经过线性均衡器得到均衡后的符号d=Gy。然后在移除信号中的DC分量后,得到等效接收符号r,表示为
注意到由于湍流效应造成的随机扰动,系统可能无法获取准确的瞬时CSI,仅能获得无湍流影响下的静态CSI,即它可以通过对足够多次信道估计的结果取期望获得。因此在式(14)中的DC偏置项中的信道矩阵是而不是H。这意味着DC分量在接收端将无法被理想地消除,从而导致接收机在信号恢复时的性能下降。另一方面,接收端处的平均接收功率可以由下式计算
考虑在下行通信场景中,接收端可能无法提前获得发射功率的信息,根据式(16),发射功率可以在接收端处通过下式进行估计
针对在湍流环境下MIMO UWOC中信道估计结果不准确的问题,本发明提出基于TS-MMSE的预编码和信道均衡矩阵设计方法,目标是根据静态CSI以及湍流强度的统计信息,在发射功率约束下联合优化预编码矩阵F和均衡矩阵G,以最小化等效接收符号与发送符合之间的MSE。该优化问题可表述为
其中MSE可具体表示为
其中(·)H和Tr(·)分别表示共轭转置和矩阵求迹,Rs=Ε(ssH)是数据源的协方差矩阵,IK是K维单位矩阵,同时定义
和
问题(18)是一个双变量优化问题,本发明采用文献[12]使用的交替优化算法进行求解。具体而言,该算法首先会初始化满足问题(18)约束的F和G,接下来对两个变量进行迭代更新,其中每次更新时算法将固定其中一个变量,而对另一个变量进行更新,以使得MSE逐渐减小至最终收敛。以下说明在F和G分别给定时,另一个变量的更新计算,该部分的推导也是本发明的创新点。
1)给定F更新G
当给定预编码矩阵F时,最优的均衡矩阵G*可通过求解下式获得:
因此可得
式(25)中的两个期望项分别可表示为
和
其中X[i,j]表示矩阵X的第(i,j)个元素,Diag(a)是关于向量a的对角矩阵。
2)给定G更新F
Horn and C.Johnson,Topics in Matrix Analysis[M],New York,NY,USA:Cambridge Univ.Press,1991.,其中vec(·)和分别表示矩阵向量化以及克罗内克积,式(19)的MSE的表达式中第一项可以被改写为
类似地,式(19)中第三项也可改写为
经过改写后,可发现式(28)和式(30)是关于vec(FH)的二次型。另外,易证式(19)的其他项都是关于vec(FH)的线性项或常数项,因为有
和
其中利用了恒等式[19]Tr(PQ)=vec(PH)Hvec(Q)[19]和式(32)中的矩阵都为实数矩阵的事实。另一方面,同时对F的功率约束式(6)也是线性约束。因此当给定均衡矩阵G时,由(18)定义的问题是以下关于F的一个典型的线性约束二次规划凸问题,用公式表示为:
其中是关于F的常数项。该凸问题的最优解F*可以通过CVX软件工具包[20]M.Grant and S.Boyd,CVX:Matlab Software forDisciplined Convex Programming,version 2.0beta[OL],http://cvxr.com/cvx,Sept.2013.进行有效求解。
综上,结合式(25)与问题(33)的解,可应用交替优化算法,最终获得最小化MSE的最优预编码和均衡矩阵。
实施例2
在实施例1的基础上,为更充分地阐述本发明所具有的有益效果,以下结合仿真分析及结果,进一步对本发明的有效性和先进性予以说明。
仿真使用随机信道模型[6]计算水下静态信道矩阵,其中关键仿真参数以如表1所示。此外,仿真中考虑近岸海域水体,它的吸收、散射和消光系数的数值为[a,b,c]T=[0.179,0.219,0.398]T[21]C.D.Mobley,Light and water:radiative transfer innatural waters[M],New York,NY,USA:Academic,1994。本发明提出的预编码和均衡矩阵设计方法也适用于其他湍流环境下的水下信道以及收发机配置参数。
表1关键仿真参数表
具体地,仿真所得水下静态信道矩阵结果为
以下将选用在MIMO VLC系统中常用的与三种预编码和均衡方法与所提出的TS-MMSE方法进行性能对比,这三种方法分别是传统基于MMSE的,基于迫零(Zero Forcing,ZF)的和基于特征值分解(Singular Value Decomposition,SVD)的预编码和均衡矩阵设计方法。
对于基于ZF的方法,其预编码和均衡矩阵分别为
其中[X][k,:]表示矩阵X的第k行向量。
其中αSVD与αZF作用一致,并由下式计算:
图3给出当系统采用PAM调制阶数为2时,四种方法在不同闪烁指数和平均发射功率下的误比特率(Bit Error Rate,BER)性能。
由图3中可见,对于进行对比的四种方法,系统的BER性能都将随着闪烁指数的增大而下降,体现出湍流效应对MIMO UOWC系统的不利影响。但相比而言,本发明所提出的TS-MMSE方法在湍流存在,即时在不同条件下均能获得最低的BER,说明本发明提出的TS-MMSE方法能在湍流环境下有效提升MIMO UOWC的传输可靠性。另外图4给出了TS-MMSE与MMSE方法时,等效接收符号与发送符号间MSE在不同闪烁指数和平均发射功率下的性能比较。
更具体的,如图4可见,在湍流存在的情况下,由于TS-MMSE方法考虑了湍流的统计特性,相比传统未考虑湍流效应的MMSE方法能够实现更低的MSE。且随着湍流强度的上升,MSE的相对下降比例更大。另外根据MSE表达式(19)的第二项可知,由于湍流效应造成的信道估计偏差与发射功率有关,因此当平均发射功率较大时,该项对MSE的影响不可忽略,导致采用传统MMSE方法时的MSE性能在平均发射功率达到某个阈值时出现恶化,这也是在图3中,采用MMSE方法的系统BER性能在当平均发射功率高于31dBm时出现恶化的原因。与之相反,本发明提出的TS-MMSE方法在优化MSE时考虑了湍流效应的影响,因此BER和MSE性能不会出现随着平均发射功率增大而恶化的情况,体现出本发明所提出的TS-MMSE方法在提升UWOC系统在对抗湍流效应方面的优越性。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (10)
1.海洋无线光通信系统的预编码和均衡矩阵设计方法,其特征在于,包括以下步骤:
S1:构建具有多条通道的MIMO UWOC系统、预编码矩阵F和均衡矩阵G;
S2:将MIMO UWOC系统的发送信号采用直流偏置向量进行处理,得到对预编码矩阵的功率约束;
S3:MIMO UWOC系统接收端接收到发送信号后,经过均衡矩阵和直流去除操作后,得到等效接收符号表达式;
S4:基于静态信道信息以及湍流效应的统计特性,获取等效接收符号与发送信号之间的MSE表达式;
S5:构建在步骤S2所得功率约束条件下最小化MSE的优化问题,其中优化变量为F和G;
S6:推导在分别给定F和G时,另一个变量在优化问题的最优解;
S7:基于步骤S6的推导,采用交替优化算法求解所构建的优化问题,对F和G进行联合优化,得到最小化MSE的最优预编码矩阵和均衡矩阵。
当检测器之间的横向距离大于1厘米时,所有子信道之间湍流衰落满足相互独立,这在实际应用中是可以满足的;因此,假设矩阵中的元素是相互独立的随机变量;而在湍流效应建模中,湍流衰落的强度由闪烁指数来评估,其计算表达式为:
其中和分别是服从高斯分布的对数幅度的均值和方差;此外,对水下湍流效应建模时,湍流衰落系数的期望通常会进行归一化来避免平均功率的变化,从而有成立,因此得出和除此之外,得到闪烁指数与对数幅度的方差之间的关系:
9.根据权利要求8所述的海洋无线光通信系统的预编码和均衡矩阵设计方法,其特征在于,所述步骤S6的具体推导过程为:
1)给定F更新G
当给定预编码矩阵F时,最优的均衡矩阵G*通过求解下式获得:
因此可得
式中,两个期望项分别表示为:
和
其中X[i,j]表示矩阵X的第(i,j)个元素,Diag(a)是关于向量a的对角矩阵;
2)给定G更新F
类似地,MSE的表达式第三项改写为:
改写后的式子是关于vec(FH)的二次型;另外,易证MSE的表达式的其他项都是关于vec(FH)的线性项或常数项,因为有:
和
其中利用恒等式Tr(PQ)=vec(PH)Hvec(Q)和上式中的矩阵都为实数矩阵的事实;另一方面,同时对F的功率约束也是线性约束,因此当给定均衡矩阵G时,最优F为一个典型的线性约束二次规划凸问题的解,该凸问题用公式表示为:
s.t.abs(F)u≤p
10.根据权利要求8所述的海洋无线光通信系统的预编码和均衡矩阵设计方法,其特征在于,所述步骤S7具体为:交替优化算法会初始化满足优化问题约束的F和G,接着对两个变量进行迭代更新,其中每次更新时算法将其中一个变量进行固定,另一个变量进行更新,以使得MSE逐渐减小至收敛,最终获得最小化MSE的最优预编码矩阵和均衡矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010477672.7A CN111682899B (zh) | 2020-05-29 | 2020-05-29 | 海洋无线光通信系统的预编码和均衡矩阵设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010477672.7A CN111682899B (zh) | 2020-05-29 | 2020-05-29 | 海洋无线光通信系统的预编码和均衡矩阵设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111682899A true CN111682899A (zh) | 2020-09-18 |
CN111682899B CN111682899B (zh) | 2021-05-04 |
Family
ID=72452790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010477672.7A Active CN111682899B (zh) | 2020-05-29 | 2020-05-29 | 海洋无线光通信系统的预编码和均衡矩阵设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111682899B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113141214A (zh) * | 2021-04-06 | 2021-07-20 | 中山大学 | 基于深度学习的水下光通信失准鲁棒盲接收机设计方法 |
CN114301525A (zh) * | 2021-12-29 | 2022-04-08 | 北京邮电大学 | 一种针对具有子连接结构的mimo vlc系统的sic预编码方法 |
CN115913840A (zh) * | 2022-11-11 | 2023-04-04 | 中国船舶重工集团公司第七一五研究所 | 一种基于矢量阵接收的mimo高速水声通信方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104038321A (zh) * | 2014-05-16 | 2014-09-10 | 桂林电子科技大学 | M-ppm光通信系统中的正交空时分组编译码方法及装置 |
CN106487446A (zh) * | 2016-06-23 | 2017-03-08 | 兰州理工大学 | 适合于强度调制/直接检测式的光空时网格码编码方法 |
US20170201321A1 (en) * | 2016-01-07 | 2017-07-13 | Ozyegin Universitesi | Adaptive multiple input multiple output (mimo) optical orthogonal frequency division multiplexing (o-ofdm) based visible light communication |
KR101786575B1 (ko) * | 2016-11-25 | 2017-11-15 | 서울과학기술대학교 산학협력단 | 다중 입출력 시스템을 기반으로 하는 프리코딩 장치 |
US20180062766A1 (en) * | 2016-08-25 | 2018-03-01 | King Abdullah University Of Science And Technology | Systems and methods for underwater illumination, survey, and wireless optical communications |
CN107872256A (zh) * | 2017-12-19 | 2018-04-03 | 郑州大学 | 无线光通信中盲检测及调制星座优化方法、存储介质 |
CN108809382A (zh) * | 2018-04-26 | 2018-11-13 | 中山大学 | 基于选择调光和预编码的多用户mimo可见光通信方法 |
CN109617576A (zh) * | 2018-12-17 | 2019-04-12 | 兰州理工大学 | 联合效应影响下大气光mimo系统平均容量的近似计算方法 |
-
2020
- 2020-05-29 CN CN202010477672.7A patent/CN111682899B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104038321A (zh) * | 2014-05-16 | 2014-09-10 | 桂林电子科技大学 | M-ppm光通信系统中的正交空时分组编译码方法及装置 |
US20170201321A1 (en) * | 2016-01-07 | 2017-07-13 | Ozyegin Universitesi | Adaptive multiple input multiple output (mimo) optical orthogonal frequency division multiplexing (o-ofdm) based visible light communication |
CN106487446A (zh) * | 2016-06-23 | 2017-03-08 | 兰州理工大学 | 适合于强度调制/直接检测式的光空时网格码编码方法 |
US20180062766A1 (en) * | 2016-08-25 | 2018-03-01 | King Abdullah University Of Science And Technology | Systems and methods for underwater illumination, survey, and wireless optical communications |
KR101786575B1 (ko) * | 2016-11-25 | 2017-11-15 | 서울과학기술대학교 산학협력단 | 다중 입출력 시스템을 기반으로 하는 프리코딩 장치 |
CN107872256A (zh) * | 2017-12-19 | 2018-04-03 | 郑州大学 | 无线光通信中盲检测及调制星座优化方法、存储介质 |
CN108809382A (zh) * | 2018-04-26 | 2018-11-13 | 中山大学 | 基于选择调光和预编码的多用户mimo可见光通信方法 |
CN109617576A (zh) * | 2018-12-17 | 2019-04-12 | 兰州理工大学 | 联合效应影响下大气光mimo系统平均容量的近似计算方法 |
Non-Patent Citations (4)
Title |
---|
S. ALIREZA NEZAMALHOSSEINI;LAWRENCE R. CHEN: "Optimal Power Allocation for MIMO Underwater Wireless Optical Communication Systems Using Channel State Information at the Transmitter", 《IEEE JOURNAL OF OCEANIC ENGINEERING》 * |
YUHANG SONG等: "Experimental demonstration of MIMO-OFDM underwater wireless optical communication", 《OPTICS COMMUNICATIONS》 * |
傅玉青;黄诚惕;杜永兆: "MIMO水下无线光通信系统性能解析研究", 《光电子·激光》 * |
寿钦等: "Gamma-Gamma大气湍流下自由空间光通信MIMO-OFDM系统误比特率分析", 《激光与光电子学进展》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113141214A (zh) * | 2021-04-06 | 2021-07-20 | 中山大学 | 基于深度学习的水下光通信失准鲁棒盲接收机设计方法 |
CN113141214B (zh) * | 2021-04-06 | 2022-03-15 | 中山大学 | 基于深度学习的水下光通信失准鲁棒盲接收机设计方法 |
CN114301525A (zh) * | 2021-12-29 | 2022-04-08 | 北京邮电大学 | 一种针对具有子连接结构的mimo vlc系统的sic预编码方法 |
CN115913840A (zh) * | 2022-11-11 | 2023-04-04 | 中国船舶重工集团公司第七一五研究所 | 一种基于矢量阵接收的mimo高速水声通信方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111682899B (zh) | 2021-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111682899B (zh) | 海洋无线光通信系统的预编码和均衡矩阵设计方法 | |
CN111682900B (zh) | 一种海洋无线光通信mimo系统的判决阈值设计方法 | |
El Ayach et al. | Interference alignment with analog channel state feedback | |
Jamali et al. | On the BER of multiple-input multiple-output underwater wireless optical communication systems | |
CN111786724B (zh) | 基于深度学习的多波长led水下可见光通信调制方法 | |
CN114567397B (zh) | 一种无线通信系统安全符号级预编码方法 | |
Guo et al. | Robust MSE equalizer design for MIMO communication systems in the presence of model uncertainties | |
Saeed et al. | Performance evaluation of MIMO FSO communication with gamma-gamma turbulence channel using diversity techniques | |
Khan et al. | Multiuser detection for MIMO-OFDM system in underwater communication using a hybrid bionic binary spotted hyena optimizer | |
Gursoy et al. | Molecular index modulation with space-time equalization | |
Yao et al. | Minimum bit error rate multiuser transmission designs using particle swarm optimisation | |
Mowlavi et al. | Bit error rate and outage analysis of MIMO-FSO communications over K-distributed atmospheric channels with imperfect feedback | |
Ntontin et al. | Performance analysis of antenna subset selection in space shift keying systems | |
Yeşilyurt et al. | Hybrid ml-mmse adaptive multiuser detection based on joint channel estimation in sdma-ofdm systems | |
Amara et al. | Optimal MU-MIMO precoder with miso decomposition approach | |
Kaur et al. | Analysis of MIMO FSO over different Modulation Techniques. | |
CN111431620B (zh) | 一种基于ppm调制的差分空间调制系统的构建方法 | |
Huang et al. | Reliable rotation transform aided MIMO scheme for underwater visible light communication | |
Alshamali et al. | Performance analysis of spatial modulation over Weibull fading channels | |
Porto et al. | Iterative Detection and Decoding for RIS-Assisted Multiuser Multiple-Antenna Systems | |
Saw et al. | Analysis of Space Shift Keying Modulation over F Turbulence Channel Model in FSO Communication | |
CN110190908A (zh) | 一种最小化isi信道中非相干大规模simo系统误码率的星座图设计方法 | |
Ho et al. | Robust minimax MSE equalizer designs for MIMO wireless communications with time-varying channel uncertainties | |
Vázquez et al. | Space-time ml receiver and its reciprocal transmitter design for interference networks | |
Li et al. | A New Photon-Counting MU-MISO Ultraviolet Communication System With MMSE Precoding in Turbulence Channels |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |