CN114665928A - 一种基于mimo-vlc系统中的电功率分配与qr-osic预编码方法 - Google Patents

一种基于mimo-vlc系统中的电功率分配与qr-osic预编码方法 Download PDF

Info

Publication number
CN114665928A
CN114665928A CN202210247771.5A CN202210247771A CN114665928A CN 114665928 A CN114665928 A CN 114665928A CN 202210247771 A CN202210247771 A CN 202210247771A CN 114665928 A CN114665928 A CN 114665928A
Authority
CN
China
Prior art keywords
matrix
electric power
power distribution
mimo
precoding
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
Application number
CN202210247771.5A
Other languages
English (en)
Other versions
CN114665928B (zh
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.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
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 Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN202210247771.5A priority Critical patent/CN114665928B/zh
Publication of CN114665928A publication Critical patent/CN114665928A/zh
Application granted granted Critical
Publication of CN114665928B publication Critical patent/CN114665928B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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/0426—Power distribution
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明公开了一种基于MIMO‑VLC系统中的电功率分配与QR‑OSIC预编码方法，属于无线通信领域。具体为：首先构建MIMO‑VLC系统模型，发射端基站发送信号，经过光信道到达接收端的PD，构建直射链路信道增益矩阵H，并采用QR分解法得到预编码矩阵F，消除接收信号的部分干扰；通过MMSE‑OSIC检测接收机，分步消除剩余干扰；然后，对消除干扰的接收信号构建最大化合速率的电功率分配优化模型，并将非凸的电功率分配优化模型分为两个子问题，分步迭代求解相关参数，CVX工具箱利用梯度下降法循环迭代参数直至收敛，进而得到最大化系统可达合速率下的电功率分配最佳值。本发明有效地兼顾了MIMO‑VLC通信系统中的通信与照明功能，能够有效地提高系统的复用增益，增强系统的误码率性能。

Description

一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码 方法
技术领域
本发明属于无线通信领域，涉及可见光MIMO的预编码，具体是一种基于MIMO-VLC(Multiple Input Multiple Output-Visible Light Communication，多输入多输出-可见光通信)系统中的电功率分配与QR-OSIC(QR decomposition-Ordered SuccessiveInterference Cancellation，QR分解-排序串行干扰消除)预编码方法。
背景技术
无线通信技术的迅速发展支撑了智享生活、智赋生产和智慧社会等各类新兴应用场景的快速普及。面向未来2030+，6G应用场景中各类智能通信业务对带宽、通信速率、信道容量、时延和移动性等指标需求的增强，可见光通信(Visible Light Communication，VLC)是满足未来高峰值速率指标需求的潜在使能技术[1]。由于VLC是采用发光二极管(LightEmitting Diodes，LED)作为光源，且可见光波段具有近400THz无需授权的频谱资源。
故VLC可在LED照明的基础上实现高速短距离通信，并能有效缓解无线通信频段资源紧张的问题[2]。此外，可见光通信具有安全、节能环保以及无辐射等特性，由于室内场景存在多个LED，MIMO技术与VLC相结合成为当下研究热点之一。
现有MIMO-VLC系统中，在同时发送不同数据时，存在多用户间干扰或单用户内不同数据流间的干扰预编码，其中典型的预编码算法包括：迫零预编码(Zero Forcing，ZF)[3]、最小均方误差预编码(Minimum Mean Square Error，MMSE)[4]、奇异值分解预编码(Singular Value Decomposition，SVD)[5]、块对角化(Block Diagonalization，BD)[6]、脏纸编码(Minimum Mean Square Error，DPC)[7]以及非线性模代数预编码(Tomlinson-Harashima Precoding，THP)[8]等。
ZF预编码的优势在于：可完全消除用户之间以及用户内部各数据流之间的干扰，劣势在于会使得加性高斯噪声被加权放大；
MMSE预编码的优势在于：综合考虑了用户干扰和信道噪声的问题，尽可能的降低了收发信号之间的误差，劣势在于信噪比较低时，可获得较好的系统性能增益，但在信噪比较高时，性能增益提升没有低信噪比时明显，且复杂度高于ZF预编码；
SVD预编码的优势在于：对收发天线的数量没有要求；
BD预编码的优势在于：抑制了多用户之间的干扰；
DPC预编码被认为是最优的预编码，在发送端对部分干扰进行预消除，在不增加发送功率的情况下可达到系统容量的理论上界，劣势在于复杂度极高；
THP预编码是在DPC基础上的改进，可获得近似MIMO系统容量的上限，劣势在于只是DPC预编码中的一种次优方案。
目前，MIMO-VLC系统中的幅度自适应调光控制存在许多挑战：由于LED具有非线性特性，发送信号必须是非负实数；多LED场景下MIMO-VLC系统信道的高相关性，对系统性能存在一定影响；同时兼容通信与照明要求；随着LED数目的增加，在空间复用条件下，系统的复用效率问题不容忽视。
因此，需要研究一种适用于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，以解决信道相关性较高时提升系统通信速率的问题。
[1]刘光毅,金婧,王启星等.6G愿景与需求:数字孪生,智能泛在[J].移动通信,2020(6).
[2]A.Sevincer,A.Bhattarai,M.Bilgi,M.Yuksel,and N.Pala,“LIGHTNETs:Smart LIGHTing and mobile optical wireless NETworks—A survey,”IEEECommun.Surveys Tuts.,vol.15,no.4,pp.1620–1641,4th Quart.2013.
[3]Joung J,Sun S.Energy Efficient Power Control for DistributedTransmitters with ZF-Based Multiuser MIMO Precoding[J].IEEE CommunicationsLetters,2013,17(9):1766-1769.
[4]Dabbagh A D,Love D J.Multiple antenna MMSE based downlinkprecoding with quantized feedback or channel mismatch[J].IEEE Transactions onCommunications,2008,56(11):1859-1868.
[5]Busche H,Vanaev A,Rohling H.SVD based MIMO Precoding andEqualization Schemes for Realistic Channel Estimation Procedures[J].Frequenz,61(7-8):146-151.
[6]Lee D,Kim K.Reliability comparison of opportunistic scheduling andBD-precoding in downlink MIMO systems with multiple users[J].IEEE Press,2010.
[7]Sumangala.G,Gowdar C V.A Novel Precoder Design for LTE MulticellMultiuser MIMO Systems Using BD,BD_DPC&ZF Precoding&It's Validation by NashEquilibrium.
[8]Souihli O,Hatakawa Y,Konishi S.An LQ-Less Implementation ofTomlinson Harashima Precoding(THP)[J].Technical Report of Ieice Rcs,2012,112(443):13-18.
发明内容
针对上述问题，本发明提出了一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，为了在发送端实现同时发送不同的信号，首先，在给定的调光水平下，确定发送信号的直流偏执；再利用QR分解对发送信号中的部分干扰进行消除；接着，根据LED的非线性约束和总电功率约束等条件构建最大化合速率的优化问题，从中求得电功率分配的次优化方案；最后，在接收端利用OSIC检测接收机消除剩余的干扰。本发明能自适应的分配电功率，提升系统复用增益，有效地兼顾VLC的高速通信和照明的需求。
所述的基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，具体的步骤如下：
步骤一、构建包括发射端基站和接收端的MIMO-VLC系统模型，初始化收发机数目以及布局空间位置；
发射端包括NT个LED，接收端包括NR个光探测器(Photo Detector,PD)；NT个LED均匀分布在室内天花板上，PD与LED的相对位置在初始化阶段根据实际需要人为调控。
步骤二、发射端基站发送信号，经过光信道到达接收端的PD，构建不同LED和PD之间的直射链路信道增益矩阵H；
信道增益矩阵H表示为：
Figure BDA0003545730710000031
其中，hj,i表示第i(1≤i≤NT)个LED到第j(1≤j≤NR)个PD视距(Line of Sight,LoS)链路的直射链路增益。
步骤三、采用QR分解法将信道增益矩阵H分解得到预编码矩阵F，消除接收信号的部分干扰；
计算公式为：
Figure BDA0003545730710000032
其中，矩阵
Figure BDA0003545730710000033
均是副对角线全为1的单位矩阵，
Figure BDA0003545730710000034
Figure BDA0003545730710000035
是上三角矩阵，
Figure BDA0003545730710000036
是正交矩阵。
由此得到预编码矩阵F：
Figure BDA0003545730710000037
步骤四、通过MMSE-OSIC检测接收机，分步消除接收信号中的剩余干扰；
具体步骤为：
步骤401、MMSE-OSIC检测接收机由一组并列的接收机组成，第一个接收机利用NR×NT的信道增益矩阵H和预编码矩阵F，计算等效信道矩阵Heff；
Figure BDA0003545730710000041
步骤402、利用等效信道矩阵Heff计算MMSE加权矩阵，并给出判决式；
wj,MMSE为MMSE加权矩阵的第j行，1≤j≤NR，计算公式为：
Figure BDA0003545730710000042
其中，
Figure BDA0003545730710000043
是单位矩阵，
Figure BDA0003545730710000044
为噪声的功率；
判决式表示为：
Figure BDA0003545730710000045
其中，
Figure BDA0003545730710000046
为发送信号的平均电功率，
Figure BDA0003545730710000047
为等效信道增益矩阵Heff的第j列向量；
步骤403、当所有NR个接收信号经过判决式后，选出最大的SINR对应的索引数据进行检测后剔除估计信号；
步骤404、将剩余的NR-1个接收信号传输到第二个接收机，第一次检测的数据流为第k路数据，此时的等效信道增益矩阵大小变为NR×(NR-1)，更新等效信道矩阵
Figure BDA0003545730710000048
Figure BDA0003545730710000049
步骤405、返回步骤402利用等效信道矩阵
Figure BDA00035457307100000410
计算MMSE加权矩阵
Figure BDA00035457307100000411
经判决式选出最大的SINR对应的索引数据进行检测，并剔除这一阶段检测出来的估计信号；
步骤405、继续对后续剩余NR-2个接收信号进行重复检测，直至最后一个接收机完成信号检测。
步骤五、对接收信号构建最大化合速率的电功率分配优化模型：
Figure BDA00035457307100000412
s.t.
Figure BDA00035457307100000413
Figure BDA00035457307100000414
Figure BDA00035457307100000415
其中，
Figure BDA00035457307100000416
优化目标是系统的可达合速率；Pj代表第j个数据流分配的电功率；
Figure BDA00035457307100000417
为等效信道增益矩阵Heff的第j行第j列的元素；
约束条件C1表示所有数据流分配的电功率之和不能超过最大值Pmax；
约束条件C2表示电功率与信号斜率的等式关系，λj是第j路数据流的发送信号斜率；
约束条件C3代表发送信号限定在LED线性区间内，F(i,j)代表矩阵F的第i行，第j列元素，ΔI是发送信号的阈值，ΔI＝min(IH-BL,BL-IL)；
步骤六、将非凸的电功率分配优化模型分为两个子问题，分步迭代求解相关参数γ、θ和λ；
第一个子问题为固定矢量P*，引入中间变量γ，将原问题转换为凸优化问题，具体表示为：
Figure BDA0003545730710000051
s.t.
C1、C2、C3
Figure BDA0003545730710000052
γj代表中间变量γ的第j个变量，
Figure BDA0003545730710000053
代表已知的固定矢量P*的第j个电功率值。
利用拉格朗日法求解第一个子问题，得到中间变量γ；
拉格朗日表达式：
Figure BDA0003545730710000054
其中，
Figure BDA0003545730710000057
代表已知电功率P*的条件下，关于中间变量γ的拉格朗日函数。
固定第一个子问题得到的中间变量γ*，将第一个子问题转为关于矢量P的表达式：
Figure BDA0003545730710000055
其中，
Figure BDA0003545730710000056
代表已知中间变量γ*的条件下，关于电功率P的拉格朗日函数。
得到第二个子问题为：
Figure BDA0003545730710000061
s.t.C1、C2、C3
进一步优化比率和，得到：
Figure BDA0003545730710000062
s.t.C1、C2、C3
C5:θj≥0
θj代表辅助变量θ的第j个中间值。
利用拉格朗日法求解第二个子问题，得到：
Figure BDA0003545730710000063
其中，β是约束条件C1的拉格朗日因子，矢量
Figure BDA0003545730710000064
是约束条件C3的拉格朗日因子。
步骤七、初始化迭代求解的参数γ(0)、θ(0)和P(0)，CVX工具箱利用梯度下降法循环迭代参数γ*、θ*和P*直至收敛，进而得到最大化系统可达合速率下的电功率分配最佳值。
本发明的优点在于：
1、一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，有效地提高了系统的复用增益；
2、一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，实现了调光控制与电功率分配联合优化；
3、一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，有效地兼顾了MIMO-VLC系统的通信和照明功能；
4、一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，部分干扰消除后可获得更高的通信速率。
附图说明
图1是本发明典型室内场景中的MIMO-VLC系统框图；
图2是本发明一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法流程图；
图3是本发明典型室内场景中的MIMO-VLC通信仿真模型；
图4是本发明中有向视距链路的信道模型；
图5是本发明与多种传统预编码的电功率分配下，系统可达速率随室内LED总数变化的对比图；
图6是本发明与多种传统预编码的电功率分配下，系统可达速率随信噪比变化的对比图；
图7是本发明与多种传统预编码的电功率分配下，系统误码率随信噪比变化的对比图。
具体实施方式
为了使本领域的技术人员更加清楚地理解和实施本发明，下面将结合附图以具体实施例的方式，详细阐述本发明实施例的技术方案。
本发明一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，其系统实现框图如图1所示，数据比特流经过M-QAM调制、Hermitian变换和IFFT变换，并经过线性变换、预编码矩阵和LED直流偏置后，将LED阵列上的电信号转化为光信号，通过光信道传输，接收端的PD将光信号转化为电信号，然后利用MMSE-OSIC检测器检测O-OFDM符号，最后对检测得到的信号进行解映射、解调和FFT变换。相比于现有线性预编码算法，本发明具有更好的抗信道间干扰性能，同时可提升系统合速率与系统复用增益。
如图2所示，本发明的具体步骤为：
步骤一、构建包括发射端基站和接收端的MIMO-VLC系统模型，初始化收发机数目以及布局空间位置；
如图3所示，发射端包括NT个LED，接收端包括NR个光探测器(Photo Detector,PD)；NT个LED均匀分布在室内天花板上，PD与LED的相对位置在初始化阶段根据实际需要人为调控。
步骤二、发射端基站发送信号，经过光信道到达接收端的PD，考虑到不同LED和PD之间的直射链路增益不同，构建基于MIMO-VLC系统的信道增益矩阵H；
设计发送信号的功率分配形式，构建室内MIMO-VLC系统的调光控制模块；具体过程如下：
步骤201、将数据流中的比特信息调制为M-QAM符号，对其进行厄尔米特(Hermitian)对称变换和IFFT变换；
Hermitian对称变换表示为：
Figure BDA0003545730710000071
其中，j＝1,…,NR，l＝1,2,…,N/2-1，N代表子载波的数目。
IFFT变换表示为：
Figure BDA0003545730710000081
步骤202、将变换后的信号进行自适应电功率分配：
Figure BDA0003545730710000082
λj是第j路数据流的发送信号斜率。
步骤203、根据MIMO-VLC系统给定的调光水平η，计算直流偏置BL；
表示为：
BL＝η(IH-IL)+IL
其中，0＜η＜1，IH代表LED最高电平，IL代表LED最低电平。
步骤204、利用QR分解的方法得到预编码矩阵F，并根据LED非线性特性得到发送信号表达式：
Figure BDA0003545730710000083
其中，
Figure BDA0003545730710000084
Λ＝diag(λ1,λ2,...,λNR)，
Figure BDA0003545730710000085
Figure BDA0003545730710000086
为实数集。
发送信号必须在LED线性范围内：|F(i,:)(Λx)|≤ΔI，其中，ΔI＝min(IH-BL,BL-IL)。
本发明中建立的信道矩阵H表示为：
Figure BDA0003545730710000087
其中，hj,i表示第i(1≤i≤NT)个LED到第j(1≤j≤NR)个PD视距(Line of Sight,LoS)链路的直射链路增益。
如图4所示，可将第i个LED和第j个PD构成链路的直流增益表示为：
Figure BDA0003545730710000088
其中，m＝-ln(2)/ln(cos(Φ1/2))表示朗泊辐射阶数，Φ1/2为LED的半功率辐射角，A表示每个PD的有效检测面积，dj,i表示第j个PD到发送端第i个LED的距离，φ表示相对于发送端平面轴线的出射角，g(ψj,i)代表接收端光学接收机的增益，可以用数学公式表示为：
Figure BDA0003545730710000091
其中，β'表示PD的折射系数，Ψc代表PD的视场角(Field of View,FoV)。
步骤三、在接收端，接收光信号首先通过PD被转换为电信号，利用信道矩阵H构造预编码矩阵F；
首先，接收信号矩阵Y可由下式得到：Y＝HF(Λx)+H(BL1)+W，
其中，
Figure BDA0003545730710000092
为噪声矩阵，矩阵W中各元素相互独立且服从高斯分布
Figure BDA0003545730710000093
Figure BDA0003545730710000094
为噪声的功率，
Figure BDA0003545730710000095
为接收信号，1为NT×1大小的全1矢量。
利用QR分解可以将H矩阵表示为：
Figure BDA0003545730710000096
其中，矩阵
Figure BDA0003545730710000097
均是副对角线全为1的单位矩阵，
Figure BDA0003545730710000098
Figure BDA0003545730710000099
是上三角矩阵，
Figure BDA00035457307100000910
是正交矩阵。
由此得到预编码矩阵F：
Figure BDA00035457307100000911
步骤四、通过MMSE-OSIC(Minimum Mean Square Error-Ordered SuccessiveInterference Cancellation，最小均方误差-连续干扰消除)检测接收机判决接收O-OFDM符号，并进一步对判决得到的向量及符号进行解调，分步消除接收信号中的剩余干扰；
具体步骤为：
步骤401、MMSE-OSIC检测接收机由一组并列的接收机组成，第一个接收机利用NR×NT的信道增益矩阵H和预编码矩阵F，计算等效信道矩阵Heff；
Figure BDA00035457307100000912
步骤402、利用等效信道矩阵Heff计算MMSE加权矩阵，并给出判决式；
wj,MMSE为MMSE加权矩阵的第j行，1≤j≤NR，计算公式为：
Figure BDA00035457307100000913
其中，
Figure BDA00035457307100000914
是单位矩阵，
Figure BDA00035457307100000915
为噪声的功率；
判决式表示为：
Figure BDA0003545730710000101
其中，
Figure BDA0003545730710000102
为发送信号的平均电功率，
Figure BDA0003545730710000103
为等效信道增益矩阵Heff的第j列向量；
步骤403、当所有NR个接收信号经过判决式后，选出最大的SINR对应的索引数据进行检测后剔除估计信号；
经过MMSE检测，均方误差达到最小，而且检测后的SINR达到最大。
步骤404、将剩余的NR-1个接收信号传输到第二个接收机，第一次检测的数据流为第k路数据，此时的等效信道增益矩阵大小变为NR×(NR-1)，更新等效信道矩阵
Figure BDA0003545730710000104
Figure BDA0003545730710000105
步骤405、返回步骤402利用等效信道矩阵
Figure BDA0003545730710000106
计算MMSE加权矩阵
Figure BDA0003545730710000107
经判决式选出最大的SINR对应的索引数据进行检测，并剔除这一阶段检测出来的估计信号；
其中，
Figure BDA0003545730710000108
代替了MMSE加权矩阵中的H，重新计算MMSE加权矩阵，并只需要计算NR-1次SINR检测排序计算式，在下一个具有最大SINR的检测符号被剔除后，对剩余的信号重复同样的过程。MMSE-OSIC信号检测的顺序会对检测信号的准确性造成一定的影响，只有将干扰最小的信号进行优先检测并剔除，循环往复，直到完成所有信号的检测。
步骤405、继续对后续剩余NR-2个接收信号进行重复检测，直至最后一个接收机完成信号检测。
接着，对按次序检测到的符号
Figure BDA0003545730710000109
重新排列为原来的次序得到符号为
Figure BDA00035457307100001010
步骤406、对检测后的符号进行逆变换：
Figure BDA00035457307100001011
接着，对变化后的信号进行FFT变换：
Figure BDA00035457307100001012
最后，对FFT变换后的频域信号进行M-QAM解调。
步骤五、对接收信号构建最大化合速率的电功率分配优化模型：
Figure BDA0003545730710000111
s.t.
Figure BDA0003545730710000112
Figure BDA0003545730710000113
Figure BDA0003545730710000114
其中，
Figure BDA0003545730710000115
优化目标是系统的可达合速率；Pj代表第j个数据流分配的电功率；
Figure BDA0003545730710000116
为等效信道增益矩阵Heff的第j行第j列的元素；
约束条件C1表示所有数据流分配的电功率之和不能超过最大值Pmax；
约束条件C2表示电功率与信号斜率的等式关系，λj是第j路数据流的发送信号斜率；
约束条件C3代表发送信号限定在LED线性区间内，F(i,j)代表矩阵F的第i行，第j列元素，ΔI是发送信号的阈值，ΔI＝min(IH-BL,BL-IL)；
步骤六、将非凸的电功率分配优化模型分为两个子问题，分步迭代求解相关参数γ、θ和λ；
第一个子问题为固定矢量P*，引入中间变量γ，将原问题转换为凸优化问题，具体表示为：
Figure BDA0003545730710000117
s.t.
C1、C2、C3
Figure BDA0003545730710000118
γj代表中间变量γ的第j个变量，
Figure BDA0003545730710000119
代表已知的固定矢量P*的第j个电功率值。
利用拉格朗日法求解第一个子问题，得到中间变量γ；
可利用拉格朗日法求解：
Figure BDA00035457307100001110
μj为中间矢量μ的第j个值。
对
Figure BDA0003545730710000121
进行求偏导：
Figure BDA0003545730710000122
可得到：
Figure BDA0003545730710000123
将
Figure BDA0003545730710000124
带入原式中可得到新的表达形式：
Figure BDA0003545730710000125
其中，
Figure BDA0003545730710000126
代表已知电功率P*的条件下，关于中间变量γ的拉格朗日函数。
固定第一个子问题得到的中间变量γ*，将第一个子问题转为关于矢量P的表达式：
Figure BDA0003545730710000127
其中，
Figure BDA0003545730710000128
代表已知中间变量γ*的条件下，关于电功率P的拉格朗日函数。
得到第二个子问题为：
Figure BDA0003545730710000129
s.t.C1、C2、C3
由于优化目标仍然不是凸优化问题，需要将比率和的问题进行进一步的化解得到如下优化问题：
Figure BDA0003545730710000131
s.t.C1、C2、C3
C5:θj≥0
θj代表辅助变量θ的第j个中间值。
利用拉格朗日法求解第二个子问题，得到：
Figure BDA0003545730710000132
其中，β是约束条件C1的拉格朗日因子，矢量
Figure BDA0003545730710000133
是约束条件C3的拉格朗日因子。
步骤七、初始化迭代求解的参数γ(0)、θ(0)和P(0)，CVX工具箱利用梯度下降法循环迭代参数γ*、θ*和P*直至收敛，得到最大化系统可达合速率下的电功率分配最佳值。
本发明与传统预编码在空间复用下电功率分配方法，系统可达合速率和误码率随信噪比、LED数目的变化趋势如图5-图7所示。如图5所示，当系统信噪比固定为30dB，本发明所提方法与传统预编码下电功率分配方法在系统可达速率随LED数目变化上有明显优势；如图6所示，当室内LED总数固定为36时，本发明所提方法与传统预编码下电功率分配方法在系统可达速率随信噪比变化上有明显优势。如图7所示，当室内LED总数固定为36时，本发明所提方法的误码率随信噪比变化的性能曲线，这意味着本发明所提方法在系统信道相关性强时可获得性能增益。
综上所述，通过实施本发明一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，不仅可实现调光控制，而且可实现高复用率、高可达合速率的通信；相比于现有的ZF预编码、MMSE预编码方案，本发明有效地避免用户内数据流之间的干扰以及信道高相关性对系统性能的影响。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。

Claims (5)

1.一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，其特征在于，具体步骤为：
首先、构建包括发射端基站和接收端检测接收机的MIMO-VLC系统模型，初始化收发机数目以及布局空间位置；
发射端基站发送信号，经过光信道到达接收端的PD，构建不同LED和PD之间的直射链路信道增益矩阵H；
信道增益矩阵H表示为：
Figure FDA0003545730700000011
其中，hj,i表示第i个LED到第j个PD视距链路的直射链路增益；
然后、采用QR分解法将信道增益矩阵H分解得到预编码矩阵F，消除接收信号的部分干扰；并通过MMSE-OSIC检测接收机，分步消除接收信号中的剩余干扰；
接着，对消除了全部干扰的接收信号构建最大化合速率的电功率分配优化模型：
Figure FDA0003545730700000012
s.t.
C1:
Figure FDA0003545730700000013
C2:
Figure FDA0003545730700000014
C3:
Figure FDA0003545730700000015
其中，
Figure FDA0003545730700000016
优化目标是系统的可达合速率；Pj代表第j个数据流分配的电功率；
Figure FDA0003545730700000017
为等效信道增益矩阵Heff的第j行第j列的元素；
Figure FDA0003545730700000018
为噪声的功率；
约束条件C1表示所有数据流分配的电功率之和不能超过最大值Pmax；
约束条件C2表示电功率与信号斜率的等式关系，λj是第j路数据流的发送信号斜率；
约束条件C3代表发送信号限定在LED线性区间内，F(i,j)代表矩阵F的第i行，第j列元素，ΔI是发送信号的阈值，ΔI＝min(IH-BL,BL-IL)；
最后、将非凸的电功率分配优化模型分为两个子问题，分步迭代求解相关参数，利用梯度下降法循环迭代参数直至收敛，进而得到最大化系统可达合速率下的电功率分配最佳值。
2.如权利要求1所述的一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，其特征在于，所述发射端包括NT个LED，接收端包括NR个光探测器PD；NT个LED均匀分布在室内天花板上，PD与LED的相对位置在初始化阶段根据实际需要人为调控。
3.如权利要求1所述的一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，其特征在于，所述采用QR分解法得到预编码矩阵F的计算公式为：
Figure FDA0003545730700000021
由此得到预编码矩阵F：
Figure FDA0003545730700000022
其中，矩阵
Figure FDA0003545730700000023
均是副对角线全为1的单位矩阵，
Figure FDA0003545730700000024
Figure FDA0003545730700000025
是上三角矩阵，
Figure FDA0003545730700000026
是正交矩阵。
4.如权利要求1所述的一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，其特征在于，所述MMSE-OSIC检测接收机分步消除接收信号中的剩余干扰，具体步骤为：
步骤401、MMSE-OSIC检测接收机由一组并列的接收机组成，第一个接收机利用NR×NT的信道增益矩阵H和预编码矩阵F，计算等效信道矩阵Heff；
Figure FDA0003545730700000027
步骤402、利用等效信道矩阵Heff计算MMSE加权矩阵，并给出判决式；
wj,MMSE为MMSE加权矩阵的第j行，1≤j≤NR，计算公式为：
Figure FDA0003545730700000028
其中，
Figure FDA0003545730700000029
是单位矩阵；
判决式表示为：
Figure FDA00035457307000000210
其中，
Figure FDA00035457307000000211
为发送信号的平均电功率，
Figure FDA00035457307000000212
为等效信道增益矩阵Heff的第j列向量；
步骤403、当所有NR个接收信号经过判决式后，选出最大的SINR对应的索引数据进行检测后剔除估计信号；
步骤404、将剩余的NR-1个接收信号传输到第二个接收机，第一次检测的数据流为第k路数据，此时的等效信道增益矩阵大小变为NR×(NR-1)，更新等效信道矩阵
Figure FDA0003545730700000031
Figure FDA0003545730700000032
步骤405、返回步骤402利用等效信道矩阵
Figure FDA0003545730700000033
计算MMSE加权矩阵
Figure FDA0003545730700000034
经判决式选出最大的SINR对应的索引数据进行检测，并剔除这一阶段检测出来的估计信号；
步骤405、继续对后续剩余NR-2个接收信号进行重复检测，直至最后一个接收机完成信号检测。
5.如权利要求1所述的一种基于MIMO-VLC系统中的电功率分配与QR-OSIC预编码方法，其特征在于，所述第一个子问题为固定矢量P*，引入中间变量γ，将原问题转换为凸优化问题，具体表示为：
Figure FDA0003545730700000035
s.t.
C1、C2、C3
C4:
Figure FDA0003545730700000036
γj代表中间变量γ的第j个变量，
Figure FDA0003545730700000037
代表已知的固定矢量P*的第j个电功率值；
利用拉格朗日法求解第一个子问题，得到中间变量γ；
拉格朗日表达式：
Figure FDA0003545730700000038
其中，
Figure FDA0003545730700000039
代表已知电功率P*的条件下，关于中间变量γ的拉格朗日函数；
固定第一个子问题得到的中间变量γ*，将第一个子问题转为关于矢量P的表达式：
Figure FDA00035457307000000310
其中，
Figure FDA00035457307000000311
代表已知中间变量γ*的条件下，关于电功率P的拉格朗日函数；
得到第二个子问题为：
Figure FDA0003545730700000041
s.t.C1、C2、C3
进一步优化比率和，得到：
Figure FDA0003545730700000042
s.t.C1、C2、C3
C5:θj≥0
θj代表辅助变量θ的第j个中间值；
利用拉格朗日法求解第二个子问题，得到：
Figure FDA0003545730700000043
其中，β是约束条件C1的拉格朗日因子，矢量
Figure FDA0003545730700000044
是约束条件C3的拉格朗日因子。
CN202210247771.5A 2022-03-14 2022-03-14 一种基于mimo-vlc系统中的电功率分配与qr-osic预编码方法 Active CN114665928B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210247771.5A CN114665928B (zh) 2022-03-14 2022-03-14 一种基于mimo-vlc系统中的电功率分配与qr-osic预编码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210247771.5A CN114665928B (zh) 2022-03-14 2022-03-14 一种基于mimo-vlc系统中的电功率分配与qr-osic预编码方法

Publications (2)

Publication Number Publication Date
CN114665928A true CN114665928A (zh) 2022-06-24
CN114665928B CN114665928B (zh) 2023-06-20

Family

ID=82029195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210247771.5A Active CN114665928B (zh) 2022-03-14 2022-03-14 一种基于mimo-vlc系统中的电功率分配与qr-osic预编码方法

Country Status (1)

Country Link
CN (1) CN114665928B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115967423A (zh) * 2022-12-28 2023-04-14 郑州大学 Mimo-vlc系统中基于几何均值分解容量最大化的预编码设计方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120250780A1 (en) * 2009-12-18 2012-10-04 Huawei Technologies Co., Ltd. Method and device for pre-coding, and method and device for decoding
CA2761461A1 (en) * 2011-04-21 2012-10-21 Ntt Docomo, Inc. Method and apparatus for determining a precoding vector for precoding data to be transmitted to a wireless device in a wireless communication system
US20150326293A1 (en) * 2014-05-07 2015-11-12 Qualcomm Incorporated Hybrid virtual antenna mapping for multiple-input multiple-output system
CN106059640A (zh) * 2016-06-30 2016-10-26 东南大学 一种基于QoS的VLC保密通信系统发射端设计方法
CN107070517A (zh) * 2017-04-13 2017-08-18 东南大学 一种可见光通信非正交多址系统预编码矩阵优化方法
CN107453795A (zh) * 2017-08-21 2017-12-08 东南大学 多用户毫米波通信系统的波束分配方法及其装置和系统
US20180270019A1 (en) * 2017-03-17 2018-09-20 Uurmi Systems Pvt Ltd System and method for detecting spatially multiplexed and space time block coded multiple-input multiple-output (mimo) signals
CN108736971A (zh) * 2018-05-28 2018-11-02 东南大学 一种可见光mimo系统中多用户收发机设计方法
CN108832978A (zh) * 2018-05-25 2018-11-16 哈尔滨工程大学 一种包含直传链路的多用户mimo中继系统的联合预编码方法
CN110266382A (zh) * 2019-05-29 2019-09-20 北京邮电大学 一种基于可见光通信mu-mimo-ofdm系统的多维度混合调光方法
CN112583485A (zh) * 2020-12-04 2021-03-30 东南大学 一种基于速率匹配的预编码与均衡联合优化方法
WO2021238634A1 (zh) * 2020-05-27 2021-12-02 华为技术有限公司 一种下行预编码方法、装置及基站

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120250780A1 (en) * 2009-12-18 2012-10-04 Huawei Technologies Co., Ltd. Method and device for pre-coding, and method and device for decoding
CA2761461A1 (en) * 2011-04-21 2012-10-21 Ntt Docomo, Inc. Method and apparatus for determining a precoding vector for precoding data to be transmitted to a wireless device in a wireless communication system
US20150326293A1 (en) * 2014-05-07 2015-11-12 Qualcomm Incorporated Hybrid virtual antenna mapping for multiple-input multiple-output system
CN106059640A (zh) * 2016-06-30 2016-10-26 东南大学 一种基于QoS的VLC保密通信系统发射端设计方法
US20180270019A1 (en) * 2017-03-17 2018-09-20 Uurmi Systems Pvt Ltd System and method for detecting spatially multiplexed and space time block coded multiple-input multiple-output (mimo) signals
CN107070517A (zh) * 2017-04-13 2017-08-18 东南大学 一种可见光通信非正交多址系统预编码矩阵优化方法
CN107453795A (zh) * 2017-08-21 2017-12-08 东南大学 多用户毫米波通信系统的波束分配方法及其装置和系统
CN108832978A (zh) * 2018-05-25 2018-11-16 哈尔滨工程大学 一种包含直传链路的多用户mimo中继系统的联合预编码方法
CN108736971A (zh) * 2018-05-28 2018-11-02 东南大学 一种可见光mimo系统中多用户收发机设计方法
CN110266382A (zh) * 2019-05-29 2019-09-20 北京邮电大学 一种基于可见光通信mu-mimo-ofdm系统的多维度混合调光方法
WO2021238634A1 (zh) * 2020-05-27 2021-12-02 华为技术有限公司 一种下行预编码方法、装置及基站
CN112583485A (zh) * 2020-12-04 2021-03-30 东南大学 一种基于速率匹配的预编码与均衡联合优化方法

Non-Patent Citations (14)

* Cited by examiner, † Cited by third party
Title
A.SEVINCER: "LIGHTNETs:Smart LIGHTing and mobile optical wireless NETworks—A survey", 《IEEE COMMUN.SURVEYS TUTS》 *
BUSCHE H: "SVD based MIMO Precoding and Equalization Schemes for Realistic Channel Estimation Procedures", 《FREQUENZ》 *
DABBAGH A D: "Multiple antenna MMSE based downlink precoding with quantized feedback or channel mismatch", 《IEEE TRANSACTIONS ON COMMUNICATIONS》 *
JOUNG J: ".Energy Efficient Power Control for Distributed Transmitters with ZF-Based Multiuser MIMO Precoding", 《IEEE COMMUNICATIONS LETTERS》 *
LEE D: "Reliability comparison of opportunistic scheduling and BD-precoding in downlink MIMO systems with multiple users", 《IEEE PRES》 *
SOUIHLI O: "An LQ-Less Implementation of Tomlinson Harashima Precoding(THP)[", 《TECHNICAL REPORT OF IEICE RCS》 *
SUMANGALA.G: "A Novel Precoder Design for LTE Multicell Multiuser MIMO Systems Using BD,BD_DPC&ZF Precoding&It\'s Validation by Nash Equilibrium", 《INTERNATIONAL JOURNAL OF SCIENTIFIC ENGINEERING AND RESEARCH (IJSER) 》 *
YANGYANG: "Joint LED Selection and Precoding Optimization for Multiple-User Multiple-Cell VLC Systems", 《IEEE INTERNET OF THINGS JOURNAL》 *
丁俊杰: "实用室内可见光通信连续干扰的消除", 《光通信研究》 *
刘光毅: "数字孪生,智能泛在", 《6G愿景与需求》 *
张颖等: "预编码室内MIMO可见光通信系统空间相关性分析", 《光电工程》 *
张颖等: "预编码室内MIMO可见光通信系统空间相关性分析", 《光电工程》, no. 03, 15 March 2020 (2020-03-15) *
杨博瑞等: "自适应功率分配可见光MIMO信道容量研究", 《应用光学》 *
杨博瑞等: "自适应功率分配可见光MIMO信道容量研究", 《应用光学》, no. 03, 15 May 2020 (2020-05-15) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115967423A (zh) * 2022-12-28 2023-04-14 郑州大学 Mimo-vlc系统中基于几何均值分解容量最大化的预编码设计方法
CN115967423B (zh) * 2022-12-28 2023-09-12 郑州大学 Mimo-vlc系统中基于几何均值分解容量最大化的预编码设计方法

Also Published As

Publication number Publication date
CN114665928B (zh) 2023-06-20

Similar Documents

Publication Publication Date Title
Liu et al. Optimized uplink transmission in multi-antenna C-RAN with spatial compression and forward
CN110266382B (zh) 一种基于可见光通信mu-mimo-ofdm系统的多维度混合调光方法
Naser et al. Rate-splitting multiple access: Unifying NOMA and SDMA in MISO VLC channels
Shen et al. Rate maximization for downlink multiuser visible light communications
CN107547130B (zh) 波束域光无线通信方法和系统
Ma et al. Coordinated beamforming for downlink visible light communication networks
Jiang et al. Energy-efficient transmission for MIMO interference channels
Ma et al. Robust MMSE linear precoding for visible light communication broadcasting systems
Sun et al. Beam domain massive MIMO for optical wireless communications with transmit lens
Utschick et al. Hybrid LISA precoding for multiuser millimeter-wave communications
Mitra et al. Precoding technique for ill-conditioned massive MIMO-VLC system
Yang et al. Joint LED selection and precoding optimization for multiple-user multiple-cell VLC systems
Naser et al. Rate-splitting multiple access for indoor visible light communication networks
CN114665928B (zh) 一种基于mimo-vlc系统中的电功率分配与qr-osic预编码方法
Zhai et al. Capacity maximized linear precoder design for spatial-multiplexing MIMO VLC systems
CN113285903B (zh) 一种大规模mimo-ofdm光无线通信系统及其低峰均比通信方法
Sifaou et al. Optimal linear precoding for indoor visible light communication system
Liu et al. Quantization and feedback of spatial covariance matrix for massive MIMO systems with cascaded precoding
Liu et al. Power allocation algorithm of optical MIMO NOMA visible light communications
Sun et al. Networked optical massive mimo communications
Li et al. Block-Level Interference Exploitation Precoding without Symbol-by-Symbol Optimization
Thakur et al. Performance analysis of energy-efficient multi-cell massive MIMO system
Feng et al. A novel hybrid dimming scheme for MU-MIMO-OFDM VLC system
Amara et al. Optimal MU-MIMO precoder with miso decomposition approach
Mostafa et al. Downlink optimization in cloud radio access networks with hybrid RF/FSO fronthaul

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