一种基于LDPC码的自由空间光通信方法及通信系统技术方案

本发明专利技术提供了一种基于LDPC码的自由空间光通信方法及通信系统，针对在激光通信中大气湍流引起的闪烁、光束漂移、扩展以及大气色散等减小接收端信号信噪比，从而影响系统误码率和通信距离、通信带宽的问题；本发明专利技术首先根据大气湍流信道的特点，提出弱湍流下的对数正态分布信道模型和强湍流下的K分布信道模型；在提出的信道模型基础上，构造一种不规则的QC‑LDPC码；推导出LLR概率密度函数，运用LLR‑BP算法译码；经置信译码BP结果显示，随着信噪比的增加，系统误码率明显下降，经实验证实，本发明专利技术所提出的基于LDPC码的自由空间光通信方法及通信系统能够有效地改善FSO系统的性能。

【技术实现步骤摘要】
一种基于LDPC码的自由空间光通信方法及通信系统
本专利技术涉及一种自由空间光通信方法及通信系统，属于光通信

技术介绍
自由空间光通信(Free-Space-Opticalcommunication，简称FSO)又称无线光通信(Wireless-Optical-communication，简称WOC)，是以激光作为信息的传输载体，大气作为传输媒介，对数据、图像、语音等信息进行双向传送的一种通信技术。一方面，FSO系统采用了光纤通信过程中光的高传输速率的特点和光电转换技术，具有高带宽、高调制速率且频谱资源占用量少的特点；另一方面，FSO系统结合了无线电波通信技术的优点，大大降低了通信的成本、线路建设周期短、系统架设也更为方便快捷。而在近地激光通信系统信号传输中，对于大气对激光通信信号干扰的分析，研究主要集中在大气的吸收和散射，直到近年来才开始对大气湍流引起的闪烁、光束漂移、扩展以及大气色散等问题进行研究，而这些因素都会影响接收端信号的信噪比，从而影响系统的误码率和通信距离、通信带宽。因此，选择何种信道模型，采用何种差错控制编码将会直接影响到光通信系统的性能。LDPC码是一类可以用非常稀疏的校验矩阵H(Parity-CheckMatrixH)或二分图(BipartiteGraph)来描述的线性分组纠错码，LDPC码的最小汉明距离会随着码长的增加而线性增加，进行后验概率迭代译码时具有随码字长度增加而比特误码率降低的特征；并且，进行迭代译码算法的LDPC码具有逼近香农限的性能。而现今的LDPC编码技术主要集中在高斯信道和瑞利信道上，这两种信道下具有更低的误码率；但AWGN信道和瑞利信道属于小尺度的衰落效应，满足不了强湍流效应下的系统通信性能需求。
技术实现思路
本专利技术的目的在于克服现有技术的不足，并提出一种基于LDPC码的自由空间光通信方法及通信系统，针对大气湍流影响下的干扰误差，选择合适的信道模型并应用更优化的信道编码从而降低接收信号的误码率。实现本专利技术目的所采用的技术方案为，一种基于LDPC码的自由空间光通信方法，包括如下步骤：(1)建立接收信号的数据模型，提出针对不同大气湍流强度的信道模型：(1-1)建立接收信号的数据模型r(t)，r(t)＝y(t)·n1(t)+n0(t)其中t为时间，y(t)为发射信号，n1(t)是大气湍流引起的乘性噪声，n0(t)是包括热噪声、放大噪声在内的加性高斯噪声，n0(t)均值为0、方差为σ2；(1-2)将大气信道分为强湍流信道和弱湍流信道，提出针对不同大气湍流强度的信道模型；a)弱湍流信道下采用对数正态分布进行模型仿真，光强V在弱湍流条件下的概率密度为其中，σx2为对数振幅方差，V0为平均光强；b)强湍流信道下采用K分布进行模型仿真，光强V在强湍流效应下的分布函数为其中V为光强，Γ(α)为Gamma函数，α为k分布阶数，Kp(x)为修正的第二类贝塞尔函数；(1-3)针对步骤(1-2)中得到的信道模型，在强弱湍流条件下得到不同的对数似然比；c)在弱湍流条件下，设发送信号y＝0时为状态D0，发送信号y＝1时为状态D1；c-a)处于状态D0时，接收信号的条件概率密度函数为：c-b)处于状态D1时，接收信号的条件概率密度函数为：d)在强湍流条件下，设发送信号y＝0时为状态D0，发送信号y＝1时为状态D1；d-a)处于状态D0时，r(t)＝n(t)，湍流以外的加性高斯噪声为0，其接收信号的条件概率密度函数为：d-b)处于状态D1时，接收信号的条件概率密度函数为：(1-4)根据步骤(1-3)中得到的不同信道下的条件概率密度函数，计算转移概率密度fiy，计算公式为：(2)编码，采用不规则的准循环LDPC码构造方法，具体内容如下；(2-1)N、M和q初始化：设定码长为N，信息位为M，码率R＝(N-M)/N，单位矩阵维数为q＝N/6；(2-2)构建校验矩阵H，校验矩阵H由子矩阵A和子矩阵B构成，H＝[AB]，子矩阵B由次级子矩阵B1和次级子矩阵B2构成，B＝[B1B2]，其中，子矩阵A的元素为单位循环子矩阵Ix，定义Ix由q阶单位矩阵I的每行向右循环移动x位得到，ak-1bk-1≥x≥1，a和b均为小于q的素数，k≥2；次级子矩阵B1的元素P中有且仅有3个非零元素且Pk-1为其中一个非零元素，其余元素均为零，定义非零元素为单位循环子矩阵Id，单位循环子矩阵Id由q阶单位矩阵I的每行向右循环移动d位得到，d为随机自然数；(2-3)检验子矩阵A，若子矩阵A为奇异矩阵，则返回步骤(3-1)重新设定N；若A为非奇异矩阵，则生成矩阵G以及长度为N的码cn，G＝[B-1·AIN-M]；cn＝un·G式中，un＝[un，0un，1…un，M-1]，un是一个随机选取的信息比特向量；(3)译码，具体步骤如下：(3-1)初始化：根据转移概率密度fiy设定信道传递给变量节点的初始概率Ln，以及对应的似然比消息数据λ(Ln)，然后设定每个变量节点n传向与其相邻的校验节点m，得到的变量节点传递给校验节点的初始信息为：λ(0)(qn,m)＝λ(Ln)；(3-2)横向处理，又称校验节点更新：即对n∈M(m)以及m∈N(n)的第l次迭代过程中变量节点传递到校验节点的信息进行更新，更新公式如下：(3-3)纵向处理，又称变量节点更新：即对n∈M(m)以及m∈N(n)的第l次迭代过程中校验节点传递到变量节点的信息进行更新，更新公式如下：(3-4)比特判决：(3-4-1)计算伪后验概率，计算公式如下：设变量当λ(l)(qn)＞0时，则否则，(3-4-2)译码输出条件判决，若或迭代次数超过规定的最大迭代次数，则结束，作为译码的有效输出值；否则回到步骤(3-1)继续迭代。本专利技术还对应提供的用于上述方法的通信系统，包括发射端和接收端，所述发射端包括电性连接的信源存储单元、LDPC编码单元、调制器和激光器，所述接收端包括电性连接的光电探测器、解调器、BP译码器和信宿存储单元。所述激光器与光电探测器通过光学发射天线和光学接收天线实现激光信号的传输。由上述技术方案可知，本专利技术提供的基于LDPC码的自由空间光通信方法，包括三个步骤：步骤一，建立接收信号的数据模型，提出针对不同大气湍流强度的信道模型；步骤二，编码，采用不规则的准循环LDPC码构造方法；步骤三，译码，根据步骤一中提供的信道模型，推导出LLR概率密度函数，运用LLR-BP算法译码。与现有技术相比，本专利技术具有以下益处：1、本专利技术根据大气湍流信道的特点，提出弱湍流下的对数正态分布信道模型和强湍流下的K分布信道模型，根据大气湍流强度大小，在强、弱湍流下建立对应信道模型，弱湍流下的小尺度衰落效应选择对数正态分布，而强湍流下多重散射效应更严重，与同概率分布模型相比，K分布能够更好地符合多重散射机制的要求，在此信道模型基础上，本专利技术提出一种基于K分布的LDPC编码，现有技术中的LDPC编码都是在高斯信道的基础上进行的，本专利技术中基于K分布进行的LDPC编码设计更加贴近大气湍流的实际情况，所设计的码型具有更优越的性能；2、由于四环的存在会降低编码的性能，增大误码率，本专利技术在编码方法中获得准循环矩阵B的步骤中，引入由3个单位循环子矩阵Ix和0元素构成的数列，单位循环子矩阵Ix在准循本文档来自技高网...

【技术保护点】
一种基于LDPC码的自由空间光通信方法，其特征在于，包括如下步骤：(1)建立接收信号的数据模型，提出针对不同大气湍流强度的信道模型：(1‑1)建立接收信号的数据模型r(t)，r(t)＝y(t)·n

【技术特征摘要】
1.一种基于LDPC码的自由空间光通信方法，其特征在于，包括如下步骤：(1)建立接收信号的数据模型，提出针对不同大气湍流强度的信道模型：(1-1)建立接收信号的数据模型r(t)，r(t)＝y(t)·n1(t)+n0(t)其中t为时间，y(t)为发射信号，n1(t)是大气湍流引起的乘性噪声，n0(t)是包括热噪声、放大噪声在内的加性高斯噪声，n0(t)均值为0、方差为σ2；(1-2)将大气信道分为强湍流信道和弱湍流信道，提出针对不同大气湍流强度的信道模型；a)弱湍流信道下采用对数正态分布进行模型仿真，光强V在弱湍流条件下的概率密度为其中，σx2为对数振幅方差，V0为平均光强；b)强湍流信道下采用K分布进行模型仿真，光强V在强湍流效应下的分布函数为其中V为光强，Γ(α)为Gamma函数，α为k分布阶数，Kp(x)为修正的第二类贝塞尔函数；(1-3)针对步骤(1-2)中得到的信道模型，在强弱湍流条件下得到不同的对数似然比；c)在弱湍流条件下，设发送信号y＝0时为状态D0，发送信号y＝1时为状态D1；c-a)处于状态D0时，接收信号的条件概率密度函数为：c-b)处于状态D1时，接收信号的条件概率密度函数为：d)在强湍流条件下，设发送信号y＝0时为状态D0，发送信号y＝1时为状态D1；d-a)处于状态D0时，r(t)＝n(t)，湍流以外的加性高斯噪声为0，其接收信号的条件概率密度函数为：d-b)处于状态D1时，接收信号的条件概率密度函数为：(1-4)根据步骤(1-3)中得到的不同信道下的条件概率密度函数，计算转移概率密度fiy，计算公式为：(2)编码，采用不规则的准循环LDPC码构造方法，具体内容如下；(2-1)N、M和q初始化：设定码长为N，信息位为M，码率R＝(N-M)/N，单位矩阵维数为q＝N/6；(2-2)构建校验矩阵H，校验矩阵H由子矩阵A和子矩阵B构成，H＝[AB]，子矩阵B由次级子矩阵B1和次级子矩阵B2构成，B＝[B1B2]，

【专利技术属性】
技术研发人员：王瑾，张琦，张晓锋，吴让仲，
申请(专利权)人：中国地质大学武汉，
类型：发明
国别省市：湖北,42