CN113162689A - 一种多输入多输出可见光通信特征信道功率分配方法 - Google Patents

Abstract

本发明公开了一种多输入多输出可见光通信特征信道功率分配方法，包括如下步骤：步骤1，通过可见光通信系统的发送端获取信道矩阵，并对其进行奇异值分解；步骤2，对信道容量进行KKT条件分析，得到功率分配数值；步骤3，确保N_t个子信道的功率参数逐级递减或不增加，并去除在总功率限制条件下无法求解出实数μ的情况；步骤4，在步骤3基础上进一步去除

时，φ_k≠0的情况，确保所选用方案均满足步骤2中的功率分配数值，并计算这些方案下的信道容量；步骤5，取步骤4所得信道容量的最大值作为最终的信道容量。本发明在保留较好的可见光通信系统性能同时，又具有低复杂度特点，在计算时间上有明显改善。

Description

一种多输入多输出可见光通信特征信道功率分配方法

技术领域

本发明属于无线通信技术领域，尤其涉及一种多输入多输出可见光通信特征信道功率分配方法。

背景技术

可见光通信(Visible Light Communication，VLC)是一种新型的宽带无线通信技术，是传统光纤通信和无线射频通信相结合的产物。它以发光二极管(Light EmittingDiode，LED)作为光源，通过发出高速明暗闪烁的可见光信号来传输信息，在接收端利用光电二极管(Photo-Diode，PD)完成光电转换，然后进行电信号的接收、再生、解调来实现信息的传递。

与频谱受到严格管制的无线射频通信相比，VLC利用光波频段进行通信，不存在频谱分配问题，无需申请频段使用执照，而且具有很大的带宽。发展VLC，在缓解无线频谱压力的同时，还符合未来无线通信的高速率发展要求。VLC可以实现高速的数据传输，目前现有的实验系统结合多输入多输出技术，已达到每秒数吉比特甚至数十吉比特的传输速率。

可达传输速率(或信道容量)是评估通讯系统性能的重要衡量标准。然而由于可见光通讯的特殊性，例如信号的非负性，使得多输入多输出的VLC的容量分析与传统无线射频通信的容量分析有本质差异。在传统的射频通讯领域中，已经有大量的文献研究预编码矩阵的优化问题以求获得最大的可达传输速率。由于注水算法可以使性能较好的子信道分配更多的功率，从而实现在功率限制情况下信道容量最大化，所以，在传统射频通信中常常采用注水算法进行特征信道功率分配。但在可见光通信系统基于奇异值(Singular ValueDecomposition，SVD)分解的预编码设计中，一个明显的现状是，信道容量是关于各子信道分配的功率参数平方而非本身的函数，这就导致与射频相比更加复杂困难，急需一种简单的方法，快速最大化信道容量。

发明内容

发明目的：针对技术背景中所提到的可见光通信系统中无法直接采用注水算法这一局限性，本发明依据注水算法思想，提出一种多输入多输出可见光通信特征信道功率分配方法。在功率分配问题中，本发明能够快速有效完成各特征信道功率分配，实现信道容量最大化。

技术方案：为实现上述目的，本发明采用的技术方案为：

一种多输入多输出可见光通信特征信道功率分配方法，包括如下步骤：

步骤1，通过可见光通信系统的发送端获取信道矩阵，并对其进行奇异值分解；

步骤2，对信道容量进行KKT条件分析，得到功率分配数值；

输入信号近似服从高斯分布，所得信道容量下界为

其中，k表示第k个子信道，N_t为子信道的个数，e为自然数，

为高斯白噪声方差，λ_k为相应于第k个子信道的特征值，φ_k为第k个子信道的功率分配数值；

通过KKT条件分析得到的功率分配数值

其中，

v_jk为步骤1对信道矩阵奇异值分解所得酉矩阵V中第(j，k)个元素；μ＝2ηln2，η为拉格朗日参数；

步骤3，确保N_t个子信道的功率参数逐级递减或不增加，并去除在总功率限制条件下无法求解出实数μ的情况；

步骤4，在步骤3基础上进一步去除

时，φ_k≠0的情况，确保所选用方案均满足步骤2中的功率分配数值，并计算这些方案下的信道容量；

步骤5，取步骤4所得信道容量的最大值作为最终的信道容量。

所述步骤1中，获取信道矩阵的方法为以下两种之一：第一种是发射端发送导频，接收端通过接收到的信号和导频计算出信道矩阵，并通过控制链路反馈给发送端；第二种是根据测量具体的发光二极管、光电二极管的空间分布、数量和辐射特性情况来计算信道矩阵。

所述获取信道矩阵的第二种方法中，根据测量具体的发光二极管、光电二极管的空间分布、数量和辐射特性情况来计算信道矩阵的公式为：

其中，h_ij代表信道矩阵第i行第j列的元素值，物理意义是第j个发光二极管和第i个光电二极管之间的直流增益，系数κ与发光二极管半角值

有关，

A是光电二极管的有效光接收面积，d_ij是第j个发光二极管与第i个光电二极管之间的距离，

为发光二极管光线与发射轴之间的夹角，而ψ为入射光线与接收轴之间的夹角，

为视场角半角。

所述步骤1中，奇异值分解的公式为Η＝UΛV^T；其中，

均为信道矩阵Η的特征值，N_t为发光二极管的个数，U和V均为酉矩阵。

所述步骤2中，关于直流偏置约束条件为abs(F)Δ＝d，其中，F是预编码矩阵，abs(·)代表求矩阵所有元素的绝对值，Δ是元素值均为Δ的列向量，Δ为原始信号的最大值，d是直流偏置列向量。

所述步骤2中，关于光功率的约束条件为1^Td＝P_t，其中，1^T是元素都为1的行向量，P_t代表总发送光功率限制，d是直流偏置列向量。

所述步骤3中，由步骤2得到的功率分配数值φ_k的取值可能性，记

G₃＝0，所以，G₁＞G₂＞G₃；由于步骤1对信道矩阵H进行奇异值分解，确保了越靠前的子信道性能越优；于是，前面的子信道有更大的概率取值G₁，后面的子信道性能较差，所以其功率参数φ_k取值不能大于前面的子信道，从而确保N_t个子信道的功率参数逐级递减或不增加；并且，功率分配应该满足总功率限制条件

其中P_t为输入端总功，如果上述等式求解不出实数μ，则意味着此情况无法满足总功率限制，所以应该去除在总功率限制条件下无法求解实数μ的情况。

有益效果：本发明针对可见光通信系统基于SVD分解的预编码设计中无法直接使用注水算法实现信道容量最大化这一限制，提出一种多输入多输出可见光通信特征信道功率分配方法，根据子信道性能好坏分配子信道功率，从而实现信道容量最大化。相比于传统的遍历方法，本发明所涉及的方法在保留可见光通信性能优越性的同时，又具有低复杂度特点，运算速度更快，可以较好地节约时间成本。

附图说明

图1为本发明实施例所涉及的可见光通讯系统模型示意图；

图2为实施例1与传统遍历方法信道容量与信噪比关系图；

图3为实施例2与现有其他预编码方法信道容量与信噪比关系图；

图4为本发明预编码方法的流程示意图。

具体实施方式

下面结合附图对本发明做进一步说明。

如图1所示为本发明的方法所涉及的可见光通讯系统模型，在可见光通信系统基于SVD分解的预编码设计中，无法直接采用注水算法分配子信道功率，从而快速最大化信道容量；而如果采用传统的遍历方法，则复杂度又较高。为解决这一不足，本发明对优化模型进行KKT条件分析，得到功率参数的解析式。然后，确保性能较好的子信道获取较大的功率参数，得到一种多输入多输出可见光通信系统特征信道功率分配方案，进而实现最大化信道容量。

如图4所示，本发明的一种多输入多输出可见光通信特征信道功率分配方法，包括如下步骤：

步骤1，通过可见光通信系统的发送端获取信道矩阵，并对其进行奇异值分解；奇异值分解的公式为Η＝UΛV^T；其中，

均为信道矩阵Η的特征值，N_t为发光二极管的个数，U和V均为酉矩阵。

其中，获取信道矩阵的方法为以下两种之一：第一种是发射端发送导频，接收端通过接收到的信号和导频计算出信道矩阵，并通过控制链路反馈给发送端；第二种是根据测量具体的发光二极管、光电二极管的空间分布、数量和辐射特性情况来计算信道矩阵。

其中，根据测量具体的发光二极管、光电二极管的空间分布、数量和辐射特性情况来计算信道矩阵的公式为：

其中，h_ij代表信道矩阵第i行第j列的元素值，物理意义是第j个发光二极管和第i个光电二极管之间的直流增益，系数κ与发光二极管半角值

有关，

A是光电二极管的有效光接收面积，d_ij是第j个发光二极管与第i个光电二极管之间的距离，

为发光二极管光线与发射轴之间的夹角，而ψ为入射光线与接收轴之间的夹角，

为视场角半角。

步骤2，对信道容量进行KKT条件分析，得到功率分配数值；

输入信号近似服从高斯分布，所得信道容量下界为

其中，k表示第k个子信道，N_t为子信道的个数，e为自然数，

为高斯白噪声方差，λ_k为相应于第k个子信道的特征值，φ_k为第k个子信道的功率分配数值；

通过KKT条件分析得到的功率分配数值

其中，

v_jk为步骤1对信道矩阵奇异值分解所得酉矩阵V中第(j，k)个元素；μ＝2ηln2，η为拉格朗日参数；

其中，关于直流偏置约束条件为abs(F)Δ＝d，其中，F是预编码矩阵，abs(·)代表求矩阵所有元素的绝对值，Δ是元素值均为Δ的列向量，Δ为原始信号的最大值，d是直流偏置列向量。

其中，关于光功率的约束条件为1^Td＝P_t，其中，1^T是元素都为1的行向量，P_t代表总发送光功率限制，d是直流偏置列向量。

步骤3，确保N_t个子信道的功率参数逐级递减或不增加，并去除在总功率限制条件下无法求解出实数μ的情况；

由步骤2得到的功率分配数值φ_k的取值可能性，记

G₃＝0，所以，G₁＞G₂＞G₃；由于步骤1对信道矩阵H进行奇异值分解，确保了越靠前的子信道性能越优；于是，前面的子信道有更大的概率取值G₁，后面的子信道性能较差，所以其功率参数φ_k取值不能大于前面的子信道，从而确保N_t个子信道的功率参数逐级递减或不增加；并且，功率分配应该满足总功率限制条件

其中P_t为输入端总功，如果上述等式求解不出实数μ，则意味着此情况无法满足总功率限制，所以应该去除在总功率限制条件下无法求解实数μ的情况。

步骤4，在步骤3基础上进一步去除

时，φ_k≠0的情况，确保所选用方案均满足步骤2中的功率分配数值，并计算这些方案下的信道容量；

步骤5，取步骤4所得信道容量的最大值作为最终的信道容量。

下面结合具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所限定的范围。

实施例1

本实施例包括如下步骤：

步骤1，通过可见光通信系统的发送端获取信道矩阵，并对其进行奇异值分解；

列出可见光通信系统信号模型的表达式为

y_r＝Hx+n (1)

发送端LED数量为N_t，接收端的PD数量为N_r，则式中：y_r为N_r×1的接收信号矢量；n为N_r×1的高斯噪声矢量，x为N_t×1发送信号矢量,x＝[x₁，x₂，x₃…x_Nt]^T，H为N_r×N_t的信道矩阵；本实例中设置LED和PD的规模N_t＝N_r＝2。

可见光MIMO的信道矩阵H如下：

其中第i行第j和元素h_ij代表第j个LED和第i个PD之间的直流增益，表达式为：

其中，系数κ与LED半角值

有关，

A是PD的有效光接收面积,d_ij是第j个LED与第i个PD之间的距离，

为LED光线与发射轴之间的夹角，而ψ为入射光线与接收轴之间的夹角，

为视场角(FOV)半角。

本实施例所采用的模型如图1所示，设置LED1坐标(-0.2,0,2.5)，LED2坐标(0.2,0,2.5)；两个探测器坐标分别为PD1：(-0.1,0,0.75),PD2：(-0.1,0,0.75)。

A＝1cm²。所以，

将H矩阵做SVD分解得到

Η＝UΛV^T (4)

其中U和V是左右奇异矩阵，分别是N_r×N_r的酉矩阵和N_t×N_t的酉矩阵，Λ是N_r×N_t的对角矩阵，对角线上是奇异值。(·)^T代表矩阵的转置。

步骤2：

发送信号为

x＝d+Fs (5)

d为N_t×1的直流偏置向量，用来保证x的非负性，F为N_t×N_t的预编码矩阵，s为N_t维的独立调制信号，假设s的符号服从均匀分布，范围是[-Δ,Δ]。

带入公式(1)的信道模型，得：

y_r＝HFs+Hd+n (6)

去除直流的影响，得到下式：

y＝y_r-Hd＝HFs+n (7)

为了满足x的非负性，令

abs(F)Δ＝d (8)

其中Δ＝[Δ,Δ，Δ…]^T；

实例中限制发射机的总光功率是P_t，由于s的均值为0，对总光功率的限制为：

1^Td＝P_t (9)

设置预编码矩阵：

F＝VΦ (10)

其中，Φ是N_t×N_t的对角阵；

公式(7)左右两边同时左乘U^T得到

y′＝U^THFs+U^Tn＝ΛΦs+n′ (11)

其中n′＝U^Tn,y′＝U^Ty，

对于第k个子信道，

y′_k＝λ_kφ_ks_k+n′_k,k＝1,2,3…N_t (12)

从公式(12)可以看出，φ_k影响着相应子信道的信噪比，所以，φ_k可以作为子信道功率分配的参数。并且，

φ_k≥0 (13)

本实施例中设置Δ＝1方便计算，总能量限制为1^Tabs(VΦ)1＝P_t

其中1^T是N_t维的元素均为1的行向量,(·)^T代表矩阵的转置。

第k个子信道的可达传输速率下界为

其中，e为自然数，

为高斯白噪声方差；

根据功率约束条件最大化信道容量，建立优化模型如下：

φ_k≥0 (15)

对上述优化模型建立拉格朗日函数：

对该函数进行KKT条件分析，共得到6个约束关系。具体如下：

t_k≥0 k＝1，2，…，N_t， (18)

t_kφ_k＝0 k＝1，2，…，N_t， (19)

φ_k≥0 k＝1，2，…，N_t， (21)

η≠0 (22)

现需要根据上述约束条件求解出φ_k，首先，需要消除参数t_k。

由公式(17)可知：

将其代入公式(18)，(19)得：

由于ε＞0，φ_k≥0，所以

并且v_k＞0，从而η≥0，然后结合公式(22)可得：

η＞0 (25)

随后，可将公式(23)改写为：

对于公式(26)，由于

将其看为关于φ_k的一个开口向上的一元二次曲线，下面将针对所有可能性进行分类讨论。

当

时，即

时，在这种情况下，由于开口向上，所以，

即：

于是，公式(26)恒成立。此时，若想公式(24)成立，必须满足φ_k＝0。

当

时，即

时，公式(26)恒成立。此时，若想公式(24)成立，则需要满足

或0。

当

时，即

时，只有当

或者

公式(26)才成立。此外，在确保公式(26)成立前提下，为满足公式(24)要求，求得：

或0。

所以，综上所述，直接根据信道容量下限进行模型优化，所得功率分配方案如下：

步骤3：

由步骤2已知φ_k的取值可能性，记

所以，G₁＞G₂＞G₃。步骤1对信道矩阵H进行SVD分解，确保了越靠前的子信道性能越优。于是，前面的子信道有更大的概率取值G₁，后面的子信道性能较差，所以其φ_k取值不能大于前面的子信道，从而确保N_t个子信道的功率参数逐级递减或不增加。并且，功率分配应该满足总功率限制条件

(P_t为输入端总功率)，如果上述等式求解不出实数μ，则意味着此情况无法满足总功率限制，所以这类情况应该去除；

步骤4：

在步骤3基础上进一步去除

时，φ_k≠0的情况，确保所选用方案均满足步骤2中的功率取值，并计算这些方案下的信道容量；

步骤5：

选用步骤4中最大的容量值作为本发明所得最终信道容量。

现重新考虑公式(27)所得功率分配方式。

或0共3种取值可能，则针对N_t个子信道，将会有3N_t种可能性。现这些可能性采用遍历方法，寻找其中信道容量最大的一个为最优解。

图2为本实施例中采用本发明提出的方法与3N_t种情况下遍历结果的对比图。从图2可以看出，本发明的方法与原遍历方法结果几乎完全相同，两条曲线几近完全重合。由此，可以说明本发明提出的方法可以较好地保留可见光通信系统性能的优越性，并且，相对于原3^Nt遍历方法，可以较好地提高运算效率，节约了时间成本。

实施例2

本实施例包括以下步骤：

步骤1：

列出可见光通信系统信号模型的表达式为

y_r＝Hx+n (1)

发送端LED数量为N_t，接收端的PD数量为N_r，则式中：y_r为N_r×1的接收信号矢量；n为N_r×1的高斯噪声矢量，x为N_t×1发送信号矢量,x＝[x₁，x₂，x₃…x_Nt]^T，H为N_r×N_t的信道矩阵；本实例中设置LED和PD的规模N_t＝N_r＝4。

可见光MIMO的信道矩阵H如下：

其中第i行第j和元素h_ij代表第j个LED和第i个PD之间的直流增益，表达式为：

其中，系数κ与LED半角值

有关，

A是PD的有效光接收面积,d_ij是第j个LED与第i个PD之间的距离，

为LED光线与发射轴之间的夹角，而ψ为入射光线与接收轴之间的夹角，

为视场角(FOV)半角。

本实施例所采用的模型与图1相似，只不过增加了LED与PD的个数，设置LED1坐标(-0.3,-0.3,2.5)，LED2坐标(-0.3,0.3,2.5)；LED3坐标(0.3,-0.3,2.5)，LED4坐标(0.3,0.3,2.5)；四个探测器坐标分别为PD1：(-0.05,-0.05,0.75),PD2：(-0.05,0.05,0.75)，PD3：(0.05,-0.05,0.75),PD4：(0.05,0.05,0.75)。

A＝1cm²。所以，

将H矩阵做SVD分解得到：

Η＝UΛV^T (4)

其中U和V是左右奇异矩阵，分别是N_r×N_r的酉矩阵和N_t×N_t的酉矩阵，Λ是N_r×N_t的对角矩阵，对角线上是奇异值。(·)^T代表矩阵的转置。

步骤2：

发送信号为

x＝d+Fs (5)

d为N_t×1的直流偏置向量，用来保证x的非负性，F为N_t×N_t的预编码矩阵，s为N_t维的独立调制信号，假设s的符号服从均匀分布，范围是[-Δ，Δ]。

带入公式(1)的信道模型，得：

y_r＝HFs+Hd+n (6)

去除直流的影响，得到下式：

y＝y_r-Hd＝HFs+n (7)

为了满足x的非负性，令

abs(F)Δ＝d (8)

其中Δ＝[Δ,Δ，Δ…]^T；

本实施例中限制发射机的总光功率是P_t，由于s的均值为0，对总光功率的限制为：

1^Td＝P_t (9)设置预编码矩阵：

F＝VΦ (10)

其中，Φ是N_t×N_t的对角阵；

公式(7)左右两边同时左乘U^T得到

y′＝U^THFs+U^Tn＝ΛΦs+n′ (11)

其中n′＝U^Tn,y′＝U^Ty，

对于第k个子信道，

y′_k＝λ_kφ_ks_k+n′_k,k＝1,2,3…N_t (12)

从公式(12)可以看出，φ_k影响着相应子信道的信噪比，所以，φ_k可以作为子信道功率分配的参数。并且，

φ_k≥0 (13)

本实施例中设置Δ＝1方便计算，总能量限制为1^Tabs(VΦ)1＝P_t

其中1^T是N_t维的元素均为1的行向量,(·)^T代表矩阵的转置。

第k个子信道的可达传输速率下界为

其中，e为自然数，

为高斯白噪声方差；

根据功率约束条件最大化信道容量，建立优化模型如下：

φ_k≥0 (15)

对上述优化模型建立拉格朗日函数：

对该函数进行KKT条件分析，共得到6个约束关系。具体如下：

t_k≥0 k＝1，2，…，N_t， (18)

t_kφ_k＝0 k＝1，2，…，N_t， (19)

φ_k≥0 k＝1，2，…，N_t， (21)

η≠0 (22)

现需要根据上述约束条件求解出φ_k，首先，需要消除参数t_k。

由公式(17)可知：

将其代入公式(18)，(19)得：

由于ε＞0，φ_k≥0，所以

并且v_k＞0，从而η≥0，然后结合公式(22)可得：

η＞0 (25)

随后，可将公式(23)改写为：

对于公式(26)，由于

将其看为关于φ_k的一个开口向上的一元二次曲线，下面将针对所有可能性进行分类讨论。

当

时，即

时，在这种情况下，由于开口向上，所以，

即：

于是，公式(26)恒成立。此时，若想公式(24)成立，必须满足φ_k＝0。

当

时，即

时，公式(26)恒成立。此时，若想公式(24)成立，则需要满足

或0。

当

时，即

时，只有当

或者

公式(26)才成立。此外，在确保公式(26)成立前提下，为满足公式(24)要求，求得：

或0。

所以，综上所述，直接根据信道容量下限进行模型优化，所得功率分配方案如下：

步骤3：

由步骤2已知φ_k的取值可能性，记

G₃＝0。所以，G₁＞G₂＞G₃。步骤1对信道矩阵H进行SVD分解，确保了越靠前的子信道性能越优。于是，前面的子信道有更大的概率取值G₁，后面的子信道性能较差，所以其φ_k取值不能大于前面的子信道，从而确保N_t个子信道的功率参数逐级递减或不增加。并且，功率分配应该满足总功率限制条件

(P_t为输入端总功率)，如果上述等式求解不出实数μ，则意味着此情况无法满足总功率限制，这类情况应该不予考虑；

步骤4：

在步骤3基础上进一步去除

时，φ_k≠0的情况，确保所选用方案均满足步骤2中的功率取值，并根据公式(14)计算这些方案下的信道容量。

现采用另外一种预编码设计方案用以对比。首先，本发明是建立在信道矩阵H的SVD分解基础上的，现在考虑基于信道矩阵H的GMD分解的预编码矩阵设计方案。将H分解为：

H＝QR_gmdG^T (28)

并且，预编码矩阵F设计为F_gmd＝εG。此外，还需要满足功率限制条件：1^Tabs(F_gmd)Δ＝P_t。其中，1^T是N_t维的元素均为1的行向量,(·)^T代表矩阵的转置。abs(·)代表取矩阵元素的绝对值，Δ＝[Δ，Δ，Δ…]^T，为了计算方便，选取Δ＝1。

在GMD分解方案下，信道容量可估计为：

其中，K为信道矩阵H的秩，e为自然数，

为高斯白噪声方差；并且，定义

为：

步骤5：

选用步骤4中满足注水算法思想方案下所得最大的容量值作为本发明所得最终信道容量。

而在基于GMD分解的预编码矩阵设计中，由于默认各子信道分配功率一样，仅需要根据功率限制条件1^Tabs(F_gmd)Δ＝P_t求解出未知参数ε，从而可以得到该方案下所估计的信道容量。

最后，比较这两种不同分解方法下的信道容量，如图3所示。可以看出基于SVD分解的预编码矩阵方案下所得信道容量较大，说明该方案下具有更高的可达传输速率。

本发明在基于SVD分解的预编码矩阵设计基础上，提出来一种多输入多输出可见光通信特征信道功率分配方案，可以确保性能较好的子信道分配较大的功率参数，从而最大化信道容量。所以，本发明在保留较好的通信性能同时，还具有低复杂度优点，节约了时间成本。

以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (7)

1.一种多输入多输出可见光通信特征信道功率分配方法，其特征在于：包括如下步骤：

步骤1，通过可见光通信系统的发送端获取信道矩阵，并对其进行奇异值分解；

步骤2，对信道容量进行KKT条件分析，得到功率分配数值；

输入信号近似服从高斯分布，所得信道容量下界为

其中，k表示第k个子信道，N_t为子信道的个数，e为自然数，

为高斯白噪声方差，λ_k为相应于第k个子信道的特征值，φ_k为第k个子信道的功率分配数值；

通过KKT条件分析得到的功率分配数值

其中，

v_jk为步骤1对信道矩阵奇异值分解所得酉矩阵V中第(j，k)个元素；μ＝2ηln2，η为拉格朗日参数；

步骤3，确保N_t个子信道的功率参数逐级递减或不增加，并去除在总功率限制条件下无法求解出实数μ的情况；

步骤4，在步骤3基础上进一步去除

时，φ_k≠0的情况，确保所选用方案均满足步骤2中的功率分配数值，并计算这些方案下的信道容量；

步骤5，取步骤4所得信道容量的最大值作为最终的信道容量。

2.根据权利要求1所述的多输入多输出可见光通信特征信道功率分配方法，其特征在于：所述步骤1中，获取信道矩阵的方法为以下两种之一：第一种是发射端发送导频，接收端通过接收到的信号和导频计算出信道矩阵，并通过控制链路反馈给发送端；第二种是根据测量具体的发光二极管、光电二极管的空间分布、数量和辐射特性情况来计算信道矩阵。

3.根据权利要求2所述的多输入多输出可见光通信特征信道功率分配方法，其特征在于：所述获取信道矩阵的第二种方法中，根据测量具体的发光二极管、光电二极管的空间分布、数量和辐射特性情况来计算信道矩阵的公式为：

其中，h_ij代表信道矩阵第i行第j列的元素值，物理意义是第j个发光二极管和第i个光电二极管之间的直流增益，系数κ与发光二极管半角值

有关，

A是光电二极管的有效光接收面积，d_ij是第j个发光二极管与第i个光电二极管之间的距离，

为发光二极管光线与发射轴之间的夹角，而ψ为入射光线与接收轴之间的夹角，

为视场角半角。

4.根据权利要求1所述的多输入多输出可见光通信特征信道功率分配方法，其特征在于：所述步骤1中，奇异值分解的公式为Η＝UΛV^T；其中，

均为信道矩阵Η的特征值，N_t为发光二极管的个数，U和V均为酉矩阵。

5.根据权利要求1所述的多输入多输出可见光通信特征信道功率分配方法，其特征在于：所述步骤2中，关于直流偏置约束条件为abs(F)Δ＝d，其中，F是预编码矩阵，abs(·)代表求矩阵所有元素的绝对值，Δ是元素值均为Δ的列向量，Δ为原始信号的最大值，d是直流偏置列向量。

6.根据权利要求1所述的多输入多输出可见光通信特征信道功率分配方法，其特征在于：所述步骤2中，关于光功率的约束条件为1^Td＝P_t，其中，1^T是元素都为1的行向量，P_t代表总发送光功率限制，d是直流偏置列向量。

7.根据权利要求1所述的多输入多输出可见光通信特征信道功率分配方法，其特征在于：所述步骤3中，由步骤2得到的功率分配数值φ_k的取值可能性，记

G₃＝0，所以，G₁＞G₂＞G₃；由于步骤1对信道矩阵H进行奇异值分解，确保了越靠前的子信道性能越优；于是，前面的子信道有更大的概率取值G₁，后面的子信道性能较差，所以其功率参数φ_k取值不能大于前面的子信道，从而确保N_t个子信道的功率参数逐级递减或不增加；并且，功率分配应该满足总功率限制条件

其中P_t为输入端总功，如果上述等式求解不出实数μ，则意味着此情况无法满足总功率限制，所以应该去除在总功率限制条件下无法求解实数μ的情况。

Images