CN108429581B - 一种用于大气激光通信的混合调制解调方法 - Google Patents

一种用于大气激光通信的混合调制解调方法 Download PDF

Info

Publication number
CN108429581B
CN108429581B CN201810040246.XA CN201810040246A CN108429581B CN 108429581 B CN108429581 B CN 108429581B CN 201810040246 A CN201810040246 A CN 201810040246A CN 108429581 B CN108429581 B CN 108429581B
Authority
CN
China
Prior art keywords
time slot
demodulation
pulse
modulation
polarization
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.)
Expired - Fee Related
Application number
CN201810040246.XA
Other languages
English (en)
Other versions
CN108429581A (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 Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201810040246.XA priority Critical patent/CN108429581B/zh
Publication of CN108429581A publication Critical patent/CN108429581A/zh
Application granted granted Critical
Publication of CN108429581B publication Critical patent/CN108429581B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/11Arrangements specific to free-space transmission, i.e. transmission through air or vacuum
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/50Transmitters
    • H04B10/501Structural aspects
    • H04B10/503Laser transmitters
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/50Transmitters
    • H04B10/516Details of coding or modulation
    • H04B10/5161Combination of different modulation schemes
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/64Heterodyne, i.e. coherent receivers where, after the opto-electronic conversion, an electrical signal at an intermediate frequency [IF] is obtained

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Optics & Photonics (AREA)
  • Optical Communication System (AREA)

Abstract

本发明公开了一种用于大气激光通信的混合调制解调方法,将多脉冲的脉冲位置调制与高阶偏振移位键控调制混合在一起,解调时对光的偏振状态提取出斯托克斯参数,偏振状态经过大气传输后几乎不变,并且利用斯托克斯参数的三维解调与一维解调或二维解调相比具有更大的编码空间和判决空间,因此具有很强的抗干扰性能,而且,采用三维解调扩大了编码空间和判决空间,在保证误码性能的前提下,能将PolSK推向多阶使用;在接收端利用提取出的斯托克斯参数判断各时隙接收信号是否为加性高斯白噪声信号,由此可判定此时隙是否接收到脉冲,最终实现对脉冲位置的解调;该方法避免使用了传统方法中采用信号幅度解调位置,简便了解调过程,还提高了解调效率。

Description

一种用于大气激光通信的混合调制解调方法
技术领域
本发明属于大气激光通信技术领域,具体涉及一种用于大气激光通信的混合调制解调方法。
背景技术
激光通信具有频带宽、速率高、容量大等优势,随着现代社会的信息量日益膨胀,如何进一步提高激光通信系统通信容量、带宽利用率成为新的研究热点。大气激光通信易受到大气湍流的影响,导致误码率提高,直接影响了激光通信系统的通信质量,提高大气激光通信系统的抗干扰性能也成为了重中之重。调制解调技术是激光通信的一项关键技术,研究具有通信容量大、带宽利用率高、抗干扰性能强等优良性能的调制方式是提高激光通信系统性能的有效途径。
脉冲位置调制(PPM,Pulse Position Modulation)是大气激光通信系统中应用广泛的一种调制方式,基于脉冲位置调制,有人提出了差分脉冲位置调制(DPPM,Differential Pulse Position Modulation)、多脉冲的脉冲位置调制(MPPM,MultiplePulse Position Modulation)等,与脉冲位置调制相比,性能有一定程度的提高,但是仍然存在着带宽利用率低,通信容量不够高的不足,难以满足大容量通信的需求。大部分脉冲位置调制的解调都是利用光的强度作为判决依据进行解调,容易受到干扰。光的偏振状态经过大气传输后几乎不发生变化,偏振移位键控(PolSK,Polarization Shift Keying)调制用光的偏振状态进行编码携带信息,是一种很有潜力的调制方式。但是偏振移位键控对光源的平均发射功率的要求很高,是其很大的一个不足。常见的偏振移位键控调制一般为二进制,即采用两个不同的偏振状态进行编码携带信息,导致偏振移位键控调制的带宽利用率和通信容量受到很大的限制。如何在保证解调端误码率保持在一定程度的条件下,寻找简单的调制解调方式把偏振移位键控推向多阶,充分利用光的偏振状态经过大气传输后几乎不发生变化这一优良性能,是本领域技术人员很想解决的问题。
现有的PPM-PolSK调制技术,采用单脉冲,带宽利用率低,一般采用两个偏振态进行编码,通信容量有限。而且有人用偏振状态的不同排列方式来表示所携带的信息,增加了调制解调实现的技术难度。
所以,寻找进一步扩大通信容量、提高带宽利用率、抗干扰性能强的调制解调方法,是本领域技术人员关心的问题。
发明内容
有鉴于此,本发明的目的是提供一种用于大气激光通信的混合调制解调方法,可以增强抗干扰性能,扩大通信容量以及提高带宽利用率。
一种用于大气激光通信的混合调制解调方法,包括:
1)在发送端,将通信数据以帧为单位进行划分,并确定每一帧通信数据的长度:
Figure BDA0001549261440000021
比特;其中,n为每帧采用的脉冲数;
Figure BDA0001549261440000022
为每帧2m个时隙中n个脉冲的组合数,m为正整数,N为偏振移位键控混合调制用于调制的偏振状态的数量,
Figure BDA0001549261440000023
表示向下取整函数;
2)、将通信数据的每一帧的前
Figure BDA0001549261440000024
比特信息用脉冲位置调制,得到携带位置信息的自然光脉冲;
3)、将通信数据的每一帧的后
Figure BDA0001549261440000025
比特信息用偏振移位键控混合调制,由此将步骤2)的自然光脉冲调制为偏振光脉冲,然后发射到大气信道中;
4)、在接收端,采用相干探测外差接收方式接收发送端发射的信号;每一个时隙接收的信号表示为
Figure BDA0001549261440000026
针对每个时隙提取出接收信号,提取出斯托克斯参数(S0,S1,S2,S3);
5)、利用每个时隙接收信号的斯托克斯参数,判断(S1,S2,S3)是否在以(0,0,0)为球心,
Figure BDA0001549261440000027
为半径的球内:
如果在所述球内,则判定该时隙接收的信号为加性高斯白噪声,即为非脉冲信号;如果不在所述球内,则不是加性高斯白噪声,即该时隙接收的为脉冲信号;
每一帧信号各个时隙判断完毕后,即完成该帧的脉冲位置解调;
6)、对于步骤5)判定出的脉冲信号,利用提取出的斯托克斯参数解出各脉冲的偏振状态,由此实现偏振态的解调。
较佳的,所述偏振状态的数量为4、6、8、12或20。
较佳的,所述偏振状态的数量为4;偏振态分别为右旋圆偏振光,左旋-60°椭圆偏振光、左旋0°椭圆偏振光以及左旋60°椭圆偏振光。
较佳的,在步骤5)中,利用提取出的斯托克斯参数解出各脉冲的偏振状态时,斯托克斯参数呈正四面体星座分布。
本发明具有如下有益效果:
(1)本发明的一种用于大气激光通信的混合调制解调方法,可以增强抗干扰性能:将多脉冲的脉冲位置调制(MPPM)与高阶偏振移位键控调制(NPolSK)混合在一起,MPPM调制所发射的脉冲不是自然光脉冲而是偏振光脉冲,解调时对光的偏振状态提取出斯托克斯参数,属于三维解调方法,一方面偏振状态经过大气传输后几乎不变,另一方面,利用斯托克斯参数的三维解调与一维解调或二维解调相比具有更大的编码空间和判决空间,因此MPPM-NPolSK调制具有很强的抗干扰性能,而且,采用三维解调扩大了编码空间和判决空间,在保证误码性能的前提下,能将PolSK推向多阶使用。
(2)在接收端利用提取出的斯托克斯参数判断各时隙接收信号是否为加性高斯白噪声信号,由此可判定此时隙是否接收到脉冲,最终实现对脉冲位置的解调;该方法避免使用了传统方法中采用信号幅度解调位置,简便了解调过程,还提高了解调效率。
(3)扩大通信容量。本调制解调方法不止利用偏振光传输MPPM的脉冲,而且偏振光本身也携带信息,这样起到了扩大通信容量的效果。
(4)扩大通信容量的同时,不占用额外带宽,带宽利用率很高。
附图说明
图1为MPPM-NPolSK调制解调系统发送端示意图;
图2为4PolSK的斯托克斯参数的正方形星座分布图;
图3为4PolSK的斯托克斯参数的正四面体星座分布图;
图4为6PolSK的斯托克斯参数的正八面体星座分布图;
图5为8PolSK的斯托克斯参数的正方体星座分布图;
图6为几种调制方式的通信容量曲线;
图7为几种调制方式的带宽利用率曲线;
图8为MPPM-NPolSK调制解调系统接收端示意图;
图9为球坐标系
Figure BDA0001549261440000041
示意图;
图10为正四面体4PolSK某一顶点处判决区域的划分;
图11为正八面体6PolSK某一顶点处判决区域的划分;
图12为正方体8PolSK某一顶点处判决区域的划分;
图13为正四面体4PolSK某一顶点处误码区域的划分方式;
图14为解调流程图;
图15为不考虑大气湍流时几种调制方式的误码率曲线;
图16为在不考虑大气湍流、弱湍流、中湍流、强湍流条件下(M,M/4)MPPM-4PolSK误码率曲线。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
一种用于大气激光通信的混合调制解调方法,激光通信系统包括发送端和接收端两个部分。
发送端主要实现以下过程:1)将输入的通信数据以帧为单位进行划分并调制,以
Figure BDA0001549261440000051
为一帧对激光载波进行多脉冲的脉冲位置调制(MPPM)和偏振移位键控混合(NPolSK)调制,其中
Figure BDA0001549261440000052
为MPPM调制每帧2m个时隙中n个脉冲的组合数,m为正整数,N为NPolSK中用于调制的偏振状态的数量,
Figure BDA0001549261440000053
表示向下取整函数。2)每一帧的前
Figure BDA0001549261440000054
信息用MPPM调制。3)每一帧的后
Figure BDA0001549261440000055
信息用NPolSK,把2)中调制后的自然光脉冲调制为偏振光脉冲。4)调制完成后把调制的光信号发射到大气信道中。
接收端采用相干探测外差接收方式,主要实现以下过程:1)接收到的光信号经过光学带通滤波器后,用偏振分束器分为x方向和y方向的线偏振光。2)本振光为45°线偏振光,经过偏振分束分为x方向和y方向的线偏振光分别与信号光的x方向和y方向的线偏振分量进行叠加。3)两个方向的信号分别经过光电探测器和带通滤波器进行滤波后,经过斯托克斯参数提取单元提取出斯托克斯参数。4)对提取出的斯托克斯参数进行两步判断:①根据判决标准1,判定提取出的斯托克斯参数如果为加性高斯白噪声,则解调出非信号时隙;如果不是加性高斯白噪声,进行第二步解调;②根据判决标准2,判断是何偏振状态,解调出信号时隙的偏振状态。
一种发送端的具体实施方式的示意图如图1所示,主要进行以下步骤:1)将输入数据以帧为单位进行划分并调制,以
Figure BDA0001549261440000056
为一帧对激光载波进行多脉冲的脉冲位置调制(MPPM)和偏振移位键控混合(NPolSK)调制,其中
Figure BDA0001549261440000057
为MPPM调制每帧2m个时隙中n个脉冲的组合数,m为正整数,N为NPolSK中用于调制的偏振状态的数量,
Figure BDA0001549261440000058
表示向下取整函数。2)在脉冲位置调制单元,将每一帧的前
Figure BDA0001549261440000059
信息首先转换为MPPM序列,根据MPPM序列对激光进行调制,使之携带脉冲位置信息。3)将MPPM调制的信号时隙所发送的光脉冲,经过偏振分束器分为偏振方向分别x方向和y方向的线偏振光。
4)光的偏振状态可以用斯托克斯参数来描述,若光的电场矢量的两个正交分量Ex,Ey分别为:
Figure BDA0001549261440000061
则斯托克斯参数为:
Figure BDA0001549261440000062
其中,ax表示Ex方向电场矢量的振幅,ay表示Ey方向电场矢量的振幅,
Figure BDA0001549261440000063
Figure BDA0001549261440000064
为Ex方向光场的相位,
Figure BDA0001549261440000065
为Ey方向光场的相位。斯托克斯参数之间满足:S0 2=S1 2+S2 2+S3 2,若把(S1,S2,S3)作为笛卡尔坐标系的三个坐标(x,y,z),可以得到以(0,0,0)为球心,半径为S0的球面,这就是邦加球。球面上的点表示不同偏振状态的偏振光,球面以内或球面以外的点表示部分偏振光,球心表示自然光。球面赤道上的点表示线偏振光,南北两级表示右旋或左旋圆偏振光,其他点则表示椭圆偏振光。
对于高阶偏振移位键控调制,采用的偏振态的数量可以为4,6,8,12,20及以上。图2为4PolSK的斯托克斯参数的正方形星座分布图,即采用四种偏振方向不同的线偏振光携带信息。图3为4PolSK的斯托克斯参数的正四面体星座分布图,即采用右旋圆偏振光和三种不同状态(左旋-60°椭圆偏振光、左旋0°椭圆偏振光以及左旋60°椭圆偏振光)的椭圆偏振光携带信息。图4为6PolSK的斯托克斯参数的正八面体星座分布图,采用右旋、左旋圆偏振光和四种偏振方向不同的线偏振光携带信息。图5为8PolSK的斯托克斯参数的正方体星座分布图,采用八种不同偏振状态的椭圆偏振光携带信息。已有成果证明,采用正四面体星座分布比采用正方形星座分布具有更好的性能。
5)偏振态调制以4PolSK的斯托克斯参数正四面体星座分布和8PolSK的斯托克斯参数正方体星座分布为具体实施实例。对于4PolSK,一种可能的星座点分布组合为(0,0,S0),
Figure BDA0001549261440000071
对于8PolSK,一种可能的星座点分布组合为
Figure BDA0001549261440000072
Figure BDA0001549261440000073
Figure BDA0001549261440000074
相应的调制参数(ax,ay,δ)可通过联立、求解方程组得到。
6)按照4)5)6)中的调制参数,用振幅调制器和相位调制器对3)中得到的x方向和y方向的线偏振光进行调制。
7)发送端进行了以上步骤之后,调制后的光信号发射至大气信道中。所发送的信号时隙的光波函数为:
Figure BDA0001549261440000075
其中P,ωs
Figure BDA0001549261440000076
分别为每帧发送光载波的平均功率,发送光载波的角频率和相位噪声,
Figure BDA0001549261440000077
Figure BDA0001549261440000078
分别为沿x方向和y方向偏振的单位矢量。δ=α-β,只需满足α与β的差为δ即可,对α和β的具体值没有要求。
8)由于在相同条件下不同调制方式传输速率的比值即为每帧携带的信息位数之比,本实施方式采用的速率比较方法为:以M=2m个时隙为一帧的单脉冲PPM调制每帧携带的信息位数为分母,(M,M/4)MPPM,(M,M/2)MPPM,(M,M/4)MPPM-4PolSK,(M,M/4)MPPM-8PolSK,(M,M/2)MPPM-4PolSK,(M,M/2)MPPM-8PolSK,每帧M=2m个时隙所携带的信息数分别做为分子,用比值来比较不同调制方式的传输速率。不同调制方式每帧所携带的信息数如下表所示:
Figure BDA0001549261440000079
如图6所示,横坐标为m,纵坐标为上文中提出的不同调制方式传输速率的比值,为简便,图例用调制方式的名称表示该调制方式与M时隙单脉冲PPM调制的每帧所携带的信息位数之比。由图6可见,对于MPPM,每帧采用的时隙数越多,每帧中的脉冲数越多,MPPM的传输容量越大。8PolSK比4PolSK每个脉冲多携带1bit信息,MPPM-8PolSK比MPPM-4PolSK每一帧多携带几倍的信息量,MPPM-NPolSK的相对于MPPM,传输容量扩大了几倍,MPPM-NPolSK比传统的单脉冲PPM的的传输容量大了几十倍。
11)各调制方式的带宽利用率为用每帧所携带的信息数与每帧的时隙数M的比值,各调制方式的带宽利用率如下表所示:
Figure BDA0001549261440000081
本实施实例对单脉冲PPM,(M,M/4)MPPM,(M,M/2)MPPM,(M,M/4)MPPM-4PolSK,(M,M/4)MPPM-8PolSK,(M,M/2)MPPM-4PolSK,(M,M/2)MPPM-8PolSK几种调制方式的带宽利用率进行了比较。如图7所示,横坐标为m,纵坐标为带宽利用率。由图7可见,MPPM-NPolSK与PPM相比带宽利用率提高了十倍左右,与MPPM相比提高了几倍。(M,M/2)MPPM-4PolSK,(M,M/2)MPPM-8PolSK的带宽利用率达到了200%以上。
(二)一种解调的具体实施方式的示意图如图8所示,主要分为以下步骤:1)接收到的光首先经过光学带通滤波器,其带宽取决于激光器的线宽,其作用是滤除大气环境中的背景光噪声的干扰使光信号通过,接收端接收到光信号表示为
Figure BDA0001549261440000082
本振光采用与发射光源相同的激光光源,其电场分量为
Figure BDA0001549261440000083
接收光和本振光经过偏振分束其x偏振分量和y偏振分量分别混合后得到
Figure BDA0001549261440000084
Figure BDA0001549261440000091
用两个相同的光电探测器分别对x方向和y方向的混合光分别进行探测,输出的电信号经过带通滤波器(其带宽和中心频率分别为2Rs和ωIF,Rs为时隙速率)滤出中频信号,得到x(t)和y(t)。5)在发送端,信号时隙发送偏振脉冲,非信号时隙不发送脉冲。因此,在接收端,对应信号时隙和非信号时隙,所接收到的信号有不同的表达式,如下表所示:
Figure BDA0001549261440000092
Figure BDA0001549261440000101
5)对于信号时隙,不考虑加性高斯白噪声时,经过光电探测器和带通滤波器后,产生差频信号,经过斯托克斯提取单元,解调出(S10,S20,S30),令
Figure BDA0001549261440000102
Figure BDA0001549261440000103
Figure BDA0001549261440000104
(S10,S20,S30)在以S1,S2,S3为坐标轴的球心为(0,0,0)、半径为S00(S00=S10 2+S20 2+S30 2)的球面上。经过接收后,产生的加性高斯白噪声,表示为(nx,ny,nz),其均值为0,方差为σ2,附加有高斯白噪声的信号表示为(S10+nx,S20+ny,S30+nz)。
6)对于非信号时隙,
Figure BDA0001549261440000105
Figure BDA0001549261440000106
只包含本振光的偏振分量,经过光电探测器之后转化成为的电信号无法通过带通滤波器,因此经过斯托克斯参数提取后,非信号时隙得到的是加性高斯白噪声,即5)中的(nx,ny,nz)。
7)判决标准1:根据已有研究结果,偏振光的偏振度和偏振状态经过传输后几乎不发生变化。接收到的斯托克斯参数(S10,S20,S30)在以S1,S2,S3为坐标轴的球心为(0,0,0)、半径为S00的球面上。考虑(S10,S20,S30)也会存在高斯白噪声的影响,取以(0,0,0)为球心,
Figure BDA0001549261440000107
为半径的球内为判决区域,如果(nx,ny,nz)处于判决区域内,则将(nx,ny,nz)判决为加性高斯白噪声。
8)判决标准2:偏振态的编码方式如图2-图5所示,对于不同的N-PolSK编码方式,使误码率最低的判决区域为:从球心向正多边形或正多面体的每一条棱做垂线,再以球心为圆心,经过这条垂线,做垂直于棱的圆,依次方法就可以将球面划分为N个全等的球面多边形,对应各调制方式理想的斯托克斯参数分别位于这些球面多边形的中心,这些球面多边形就是相应的偏振状态的判决区域,这样的判决方式根据不同的NPolSK充分利用了整个球面作为判决区,可以使误码率最低。
图10、图11、图12是从某一顶点看过去的俯视图,分别表示正四面体4PolSK、正八面体6PolSK,正方体8PolSK的判决区域划分,图10、11和12中弧Q1Q2Q3、弧Q1Q2Q3Q4和弧Q1Q2Q3分别为以上三种编码方式的对该顶点对应的偏振状态的判决区域。N取值不同的PolSK寻找判决区域的方式相同,判决区域也都为球面正多边形,区别是球面正多边形的边数不同。
以图13中的4PolSK的正四面体编码方式为例,对于如图9所示的球坐标系
Figure BDA0001549261440000111
对于以D点所在直线为竖直的球坐标系,圆C1为弧Q1Q2Q3的内切圆,所在坐标为θ1,圆C2为弧Q1Q2Q3的外接圆所在坐标为θ2,定义
Figure BDA0001549261440000112
判决产生误码的区域为弧Q1Q2Q3以外的区域,分为两部分,一部分为圆C2以外的部分(θ∈[θ2,π],φ∈[0,2π]),一部分为圆C2以内弧Q1Q2Q3以外的部分(θ∈[θ12],φ的大小为3α(θ))。
不同的NPolSK的判决区域的参数如下表所示:
Figure BDA0001549261440000113
9)在解调单元进行两步解调:首先根据7)中的判决标准1判断是否为加性高斯白噪声,如果是,则该信号为非信号时隙,如果不是,则该信号对应信号时隙;第二步,根据8)中的判决标准2,解调出第一步得到的信号时隙的偏振状态。解调流程图如图14所示。
10)经过以上步骤的解调,产生误码的情况为:经过判决标准1把加性高斯白噪声判决为Stokes参数,对应的误码率为Pe1;在判决标准1不发生错误的前提下,经过判决标准2,错判脉冲的偏振状态,对应的误码率为Pe2。所以,总的误码率为:
Pe=Pe1+(1-Pe1)Pe2
11)定义信噪比为:SNR=μ2/2σ2,μ2=R2PPlo(ax 2+ay 2),其中P为每帧的平均发射功率。
12)把(nx,ny,nz)正确判决为加性高斯白噪声的概率为:
Figure BDA0001549261440000121
其中
Figure BDA0001549261440000122
表示在正确判决范围内积分,
Figure BDA0001549261440000123
为(nx,ny,nz)的联合概率密度分布。所以Pe1=1-Pc1
Figure BDA0001549261440000124
12)激光载波经过大气传输,受弱、中、强大气湍流影响而产生的光强起伏概率分布模型均符合Gamma-Gamma分布,其概率密度函数为:
Figure BDA0001549261440000125
其中,Kn(·)是n阶第二类修正贝塞尔函数,Γ(·)是伽马函数,α和β分别代表大规模和小规模的光强起伏,表示为:
Figure BDA0001549261440000131
Figure BDA0001549261440000132
其中,
Figure BDA0001549261440000133
为无量Rytov方差,可用来判定光强起伏强度。δR 2<1时,大气湍流为弱湍流;δR 2≈1时,大气湍流为中湍流;δR 2>1时,大气湍流为强湍流。Cn 2为大气折射率结构常数,k=2π/λ为波数,L为传播距离。
13)(S1,S2,S3)的联合概率密度分布为
Figure BDA0001549261440000134
ρ>0,θ∈[0,π],φ∈[0,2π],其中,I0是0阶修正的第一类贝塞尔函数,ρ和θ与
Figure BDA0001549261440000135
是独立的,
Figure BDA0001549261440000136
在[0,2π]上均匀分布。使用判决准则2判决的误码率为
Figure BDA0001549261440000137
对于N-PolSK,每个偏振状态发送的概率相等,则对任意一个偏振状态解调与其他偏振状态解调时的正确或错误概率都相等,所以
Figure BDA0001549261440000138
Figure BDA0001549261440000139
其中,
Figure BDA00015492614400001310
Figure BDA00015492614400001311
在大气湍流的影响下,
Figure BDA00015492614400001312
所以,MPPM-NPolSK的总的误码率表达式为:
Figure BDA0001549261440000141
14)不考虑大气湍流时,(M,M/4)MPPM-4PolSK,(M,M/4)MPPM-8PolSK,(M,M/2)MPPM-4PolSK,(M,M/2)MPPM-8PolSK四种调制方式的误码率曲线如图14所示,其中4PolSK的星座图分布为正四面体,8PolSK的星座图分布为正方体。在相同条件下,四种调制方式的误码性能按(M,M/2)MPPM-8PolSK,(M,M/2)MPPM-4PolSK(M,M/4)MPPM-8PolSK,,(M,M/4)MPPM-4PolSK的顺序依次变优,对于(M,M/4)MPPM-4PolSK,在信噪比只需25dB时,误码率即可达到10-10
15)分别在不考虑大气湍流,弱湍流:α=14.55,β=13.00,中湍流:α=6.05,β=4.47强湍流:α=4.52,β=1.40的条件下,(M,M/4)MPPM-4PolSK误码率曲线如图15所示。由图可见弱湍流时,该调制方式几乎不受影响,在中湍流和强湍流条件下误码率会升高,提高信噪比可以降低误码率。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种用于大气激光通信的混合调制解调方法,其特征在于,包括:
1)在发送端,将通信数据以帧为单位进行划分,并确定每一帧通信数据的长度:
Figure FDA0002418329140000011
比特;其中,n为每帧采用的脉冲数;
Figure FDA0002418329140000012
为每帧2m个时隙中n个脉冲的组合数,m为正整数,N为偏振移位键控混合调制用于调制的偏振状态的数量,
Figure FDA0002418329140000013
表示向下取整函数;
2)、将通信数据的每一帧的前
Figure FDA0002418329140000014
比特信息用脉冲位置调制,得到携带位置信息的自然光脉冲;
3)、将通信数据的每一帧的后
Figure FDA0002418329140000015
比特信息用偏振移位键控混合调制,由此将步骤2)的自然光脉冲调制为偏振光脉冲,然后发射到大气信道中,具体为:
将步骤2)得到的自然光脉冲,经过偏振分束器分为偏振方向分别x方向和y方向的线偏振光;
光的偏振状态用斯托克斯参数来描述,若光的电场矢量的两个正交分量Ex,Ey分别为:
Figure FDA0002418329140000016
则斯托克斯参数为:
Figure FDA0002418329140000017
其中,ax表示Ex方向电场矢量的振幅,ay表示Ey方向电场矢量的振幅,
Figure FDA0002418329140000018
Figure FDA0002418329140000019
为Ex方向光场的相位,
Figure FDA00024183291400000110
为Ey方向光场的相位,斯托克斯参数之间满足:S0 2=S1 2+S2 2+S3 2
对于4PolSK,一种可能的星座点分布组合为(0,0,S0),
Figure FDA00024183291400000111
Figure FDA0002418329140000021
对于8PolSK,一种可能的星座点分布组合为
Figure FDA0002418329140000022
Figure FDA0002418329140000023
相应的调制参数(ax,ay,δ)通过联立、求解方程组得到;
按照以上获得的调制参数,用振幅调制器和相位调制器对3)中得到的x方向和y方向的线偏振光进行调制;
4)、在接收端,采用相干探测外差接收方式接收发送端发射的信号;每一个时隙接收的信号表示为
Figure FDA0002418329140000024
针对每个时隙提取出接收信号,提取出斯托克斯参数(S0,S1,S2,S3);
5)、利用每个时隙接收信号的斯托克斯参数,判断(S1,S2,S3)是否在以(0,0,0)为球心,
Figure FDA0002418329140000025
为半径的球内:
如果在所述球内,则判定该时隙接收的信号为加性高斯白噪声,即为非脉冲信号;如果不在所述球内,则不是加性高斯白噪声,即该时隙接收的为脉冲信号;
每一帧信号各个时隙判断完毕后,即完成该帧的脉冲位置解调;
6)、对于步骤5)判定出的脉冲信号,利用提取出的斯托克斯参数解出各脉冲的偏振状态,由此实现偏振态的解调。
2.如权利要求1所述的用于大气激光通信的混合调制解调方法,其特征在于,所述偏振状态的数量为4、6、8、12或20。
3.如权利要求2所述的用于大气激光通信的混合调制解调方法,其特征在于,所述偏振状态的数量为4;偏振态分别为右旋圆偏振光,左旋-60°椭圆偏振光、左旋0°椭圆偏振光以及左旋60°椭圆偏振光。
4.如权利要求3所述的用于大气激光通信的混合调制解调方法,其特征在于,在步骤5)中,利用提取出的斯托克斯参数解出各脉冲的偏振状态时,斯托克斯参数呈正四面体星座分布。
CN201810040246.XA 2018-01-16 2018-01-16 一种用于大气激光通信的混合调制解调方法 Expired - Fee Related CN108429581B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810040246.XA CN108429581B (zh) 2018-01-16 2018-01-16 一种用于大气激光通信的混合调制解调方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810040246.XA CN108429581B (zh) 2018-01-16 2018-01-16 一种用于大气激光通信的混合调制解调方法

Publications (2)

Publication Number Publication Date
CN108429581A CN108429581A (zh) 2018-08-21
CN108429581B true CN108429581B (zh) 2020-08-11

Family

ID=63155906

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810040246.XA Expired - Fee Related CN108429581B (zh) 2018-01-16 2018-01-16 一种用于大气激光通信的混合调制解调方法

Country Status (1)

Country Link
CN (1) CN108429581B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111130646B (zh) * 2019-12-13 2023-07-18 重庆邮电大学 一种抗时延抖动的高速率mppm星座映射方法
CN111431566B (zh) * 2020-03-20 2022-04-08 哈尔滨工业大学 一种应用于大气激光通信空间调制系统的多路光信号快速搜索方法
CN111510221B (zh) * 2020-04-17 2021-05-07 中山大学 一种多偏振光无线通信的调制解调方法及系统
CN112383354B (zh) * 2020-10-29 2022-07-08 深圳大学 基于光编码的空间-偏振混合分集自由空间光通信系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102208943A (zh) * 2011-05-09 2011-10-05 中国科学技术大学 一种用于激光通信的调制解调方法
CN103795472A (zh) * 2014-03-04 2014-05-14 上海交通大学 基于单个双驱动调制器的三维调制超宽带信号的产生方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4841528A (en) * 1988-09-06 1989-06-20 California Institute Of Technology Frequency doubled, cavity dumped feedback laser
CN100574156C (zh) * 2007-06-18 2009-12-23 西安理工大学 无线激光通信ppm偏振调制解调方法
CN102223177B (zh) * 2011-06-21 2014-10-15 中国科学院上海技术物理研究所 一种基于单光子探测的超远距离光通信系统及方法
CN102780531B (zh) * 2012-08-01 2014-09-24 清华大学深圳研究生院 一种用于无线激光通信的调制方法及解调方法
CN104779997B (zh) * 2014-01-09 2017-09-22 中国科学院西安光学精密机械研究所 基于Stokes参量识别的偏振调制空间激光通信方法
WO2017205652A1 (en) * 2016-05-27 2017-11-30 Schafer Aerospace, Inc. System and method for high speed satellite-based free-space laser communications using automatic gain control

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102208943A (zh) * 2011-05-09 2011-10-05 中国科学技术大学 一种用于激光通信的调制解调方法
CN103795472A (zh) * 2014-03-04 2014-05-14 上海交通大学 基于单个双驱动调制器的三维调制超宽带信号的产生方法

Also Published As

Publication number Publication date
CN108429581A (zh) 2018-08-21

Similar Documents

Publication Publication Date Title
CN108429581B (zh) 一种用于大气激光通信的混合调制解调方法
US10833911B2 (en) Modulation systems and methods in optical transceivers
US6996189B1 (en) Symmetric spherical QAM constellation
CN111884975B (zh) 基于时延-多普勒域的索引调制解调方法和系统
US20200252318A1 (en) Machine learning training system for identification or classification of wireless signals
JP2007506291A (ja) 無線通信用偏波状態技術
EP0527966A1 (en) MULTILEVEL CONSISTENT OPTICAL SYSTEM.
US20120051452A1 (en) Modified coded hybrid subcarrier amplitude phase polarization modulation
CN108173599A (zh) 基于DENCLUE聚类的Stokes空间相干光调制格式识别方法
JPS6324342B2 (zh)
US11218355B2 (en) Multi-dimensional signal encoding
CN109076039A (zh) 使用多维非均匀星座的编码和调制装置
Yen et al. Multiple-access bosonic communications
CN106788626B (zh) 一种能够获得二阶发射分集的改进正交空间调制传输方法
JP2007503734A (ja) 下位互換性のあるマルチキャリア送信システム
CN110958051A (zh) 一种基于参数调整提升noma vlc系统误码率性能的方法
US20170353247A1 (en) Constellation design for use in communication systems
Yu et al. Group-based CSS modulation: A novel enhancement to LoRa physical layer
CN105991508A (zh) 一种8qam调制方法及其系统
CN106878228A (zh) 一种基于四维信号星座图映射与解映射的数字通信方法
US10236989B2 (en) Data transport using pairwise optimized multi-dimensional constellation with clustering
Swetha et al. Energy detection spectrum sensing on DPSK modulation transceiver using GNU radio
CN107040491B (zh) 一种基于16进制三维格子星座图映射的数字通信系统
CN109005139A (zh) 一种正交幅度调制解调方法及其星座图优化方法
EP1469628A2 (en) Constellation design for a multiple access communication systems

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200811

Termination date: 20210116