湍流大气信道中的激光脉冲传输时间展宽及信道传递函数测量方法
【专利摘要】本发明公开一种湍流大气信道中的激光脉冲传输时间展宽及信道传递函数测量方法，属于无线激光通信【技术领域】。本方法将窄带高斯光脉冲序列发射到湍流大气信道中传输，在接收端对光脉冲序列进行探测和波形采样，并将采样数据保存为文件；在对脉冲波形采样数据求系综平均的基础上，用数据拟合方法解算出系综平均高斯脉冲波形参数，进而求出湍流大气信道传递函数。利用本方法可得到湍流大气信道传递函数的解析表达式，这使得基于实测数据对湍流大气信道特性进行理论分析变得更加容易。
【专利说明】湍流大气信道中的激光脉冲传输时间展宽及信道传递函数测量方法
【技术领域】
[0001]本发明属于无线激光通信【技术领域】，涉及一种湍流大气信道中的激光脉冲传输时间展宽及信道传递函数测量方法。
【背景技术】
[0002]无线激光通信具有保密性好、部署机动灵活、无需频带申请等优点，近年来受到人们的普遍关注。当无线激光通信系统工作在地球大气环境中时，地球大气会造成激光脉冲在时间上出现展宽。对于使用强度调制/直接探测方式的无线激光通信系统，通信激光脉冲时间展宽会引起码间串扰，造成通信系统性能降低。因此对实际湍流大气信道中的激光脉冲时间展宽进行测量，以此为依据进而获得湍流大气信道传递函数，对大气环境中的无线激光通信系统性能分析及优化设计具有重要意义。
[0003]发表在《Applied Optics》37卷33期7655-7660页和SPIE会议论文集3266卷231-240页的两篇论文通过理论建模研究指出:窄带高斯光脉冲信号经弱湍流大气信道传输后，在近场和远场条件下接收机收到的系综平均光脉冲信号在时间上都保持高斯函数形状，只是系综平均光脉冲信号的时间宽度变得更宽，即发生了脉冲时间展宽。发表在《Optics Express)) 20卷7期7749-7757页的论文发现,窄带高斯光脉冲信号经强瑞流大气信道传输后，接收机收到的系综平均光脉冲信号在时间上也保持高斯函数形状。对于实际无线激光通信系统来说，接收机收到的系综平均光脉冲时间宽度与发射机发射的光脉冲时间宽度之比，能够从统计上定量地描述光脉冲经湍流大气信道传输后发生的时间展宽程度，是一个重要的物理量。对于高斯光脉冲来说，可以定义高斯光脉冲经湍流大气信道传输后的系综平均时间展宽比为接收机收到的高斯光脉冲系综平均时间宽度与发射机发射的高斯光脉冲时间宽度之比。在数学上，发射机发射的光脉冲可表示为Pi(t)=Vi(t)exp(-iω0t), Vi (t)表示脉冲时间形状，ωo表示光载波角频率，接收机收到的系综平均光脉冲可表示为po(t)=Vo(t)exp(-iωt),Vo(t)表示系综平均脉冲时间形状。从线性系统的角度，可得Po(ω)=H(ω)Pi(ω), Po(ω)为Po(t)的傅里叶变换，Pi(ω)为Pi (t)的傅里叶变换，Η(ω)表示湍流大气信道的传递函数。由于发射机发射的光脉冲的峰值和接收机收到的系综平均光脉冲的峰值通常是不同的，因此Η(ω)包含一个不依赖于角频率ω的比例系数Cko比例系数Ck只表征了光脉冲经湍流大气信道传输所受的衰减情况，其受激光束散角、光束半径、传输距离等多种因素影响。实际上，在分析光脉冲湍流大气传输时间展宽以及脉冲形状变化时，令比例系数Ck为1也不会影响分析结果。因此，在测量湍流大气信道的传递函数时，可不测比例系数ck。从前面的分析可知，只要知道湍流大气信道的传递函数Η(ω)和发射机发射的光脉冲信号Pi (t)，即可求出poω)，进而用逆傅里叶变换得到po(t)，在此基础上就可分析湍流大气信道对通信光脉冲传输产生的时间展宽和形状变化影响。所以测量实际湍流大气信道的传递函数非常有意义。由于高斯函数的傅里叶变换仍然是高斯函数，即vo(t)kexp(-t2/b2)的傅里叶变换为
【权利要求】
1.一种湍流大气信道中的激光脉冲传输时间展宽及信道传递函数测量方法，其特征在于，所需的设备及实现步骤如下: 目的在于提供一种湍流大气信道中的激光脉冲传输时间展宽及信道传递函数测量方法；本方法需要使用激光信号发射装置和激光信号接收与探测装置；激光信号发射装置包括激光器101、信号发生器102、电光调制器103、光学发射系统104 ;激光信号接收与探测装置包括光学接收系统201、光电探测器202、示波器203 ;从激光器101发出的激光信号AOOl入射到电光调制器103的光信号输入端；信号发生器102产生的电脉冲信号A002传送到电光调制器103的电信号输入端；电光调制器103将电脉冲信号A002调制到激光信号AOOl上，调制后的激光信号A003再经光学发射系统104发射到湍流大气信道中；调制后的激光信号A003经湍流大气信道传输后，通过光学接收系统201入射到光电探测器202上，光电探测器202的输出电信号再传送到示波器203的输入端； 本方法的第一部分探测并记录激光脉冲经湍流大气信道传输后的波形采样数据，具体步骤如下: StepOOl:将激光信号发射装置和激光信号接收与探测装置分别置于湍流大气信道的两端；开通激光器101,使其输出连续激光信号；设置信号发生器102,使其输出一电脉冲序列A004 ;在电脉冲序列A004中，各脉冲都为形状相同的高斯脉冲，相邻两个脉冲的峰值的间隔时间为T，每个脉冲的时间宽度为τ，且相邻两个脉冲之间互不重叠；电光调制器103输出的光脉冲序列Α005经光学发射系统104进入瑞流大气信道传输； St印002:使光学接收系统201和光学发射系统104的光轴对准，开通示波器203，设置示波器203的采样时间间隔At=T/M，M为一个正奇数；并用示波器203记录光电探测器202探测到的脉冲序列波形采样数据A005 ； 通过不断重复StepOOl和Step002,可以记录下大量的脉冲序列波形采样数据A005 ;每执行一次StepOOl和Step002可得到一组脉冲序列波形采样数据A005，并将其保存为一个单独的波形采样文件；` 本方法的第二部分根据保存的所有波形采样文件中的脉冲序列波形采样数据A005，计算高斯脉冲经湍流大气信道传输后的平均时间展宽比，具体步骤如下: SteplOl:创建一个具有M个元素的数组INTP，并将数组INTP的每个元素都设置为O ;令NUM=O ;对本方法的第一部分中保存的所有波形采样文件中的每个波形采样文件FWjA行如下步骤: SteplOl-1:读取波形采样文件FW，将其中的脉冲序列波形采样数据A005的采样序列赋值给一个数组W，数组W的第i个元素表示为W[i]，i表示数组元素的编号，I < i ^Ν,Ν为采样序列的采样个数；数组W中的每个元素与一个特定的时刻对应； St印101-2:令脉冲序列波形采样数据Α005的采样序列的第一个采样对应的时刻为`0，则W[i]对应的采样对应的时刻为(1-1) At ;计算脉冲序列波形采样数据A005中的第一个脉冲对应的波形采样序列中的幅度值最大的采样对应的数组W的元素的编号j ;如果j-(M-1)/2 < 1，则令]_=」+皿;
St印 101-2-1:令 idx=j，n=l, NUM=NUM+1 ；
St印 101-2-2:令 INTP [n] =INTP [n] +W[ j- (M-1) /2+n-l]，INTP [n]表示数组 INTP 的第 η个元素；St印101-2-3:令n=n+l ;如果η < M，则转步骤St印101-2-2，否则转步骤St印101-2-4 ； St印101-2-4:如果idx+M+(M-l)/2 ( N，则令j=j+M，转步骤St印101-2-5，否则转步骤Stepl01-3 ；
St印 101-2-5:令 n=l，NUM=NUM+1 ；
St印101-2-6:令 INTP[n]=W[j-(M-1)/2+n-l]+INTP[η]； St印101-2-7:令n=n+l ;如果n ^ Μ,转步骤St印101-2-6，否则转步骤St印101-2-8 ； SteplOl-2-8:令 idx=idx+M,转步骤 SteplOl-2-4: St印101-3 ;针对波形采样文件FW所需执行的操作结束；
Stepl02:令 n=l ；
St印 102-1:令 INTP [n] = INTP [n]/NUM ； St印102-2:令η=η+1 ;如果η < Μ，转步骤St印102-1，否则转步骤St印102-3 ； St印102-3:令数组INTP的第一个元素对应的时刻为0，则数组INTP的第i个元素对应的时刻为(1-1) At，I≤i≤M ;根据数组INTP的各个元素值及其对应的时刻，可以得到一个数据序列(Xpyi) ,Xi=Q-1) Δ t, Yi=INTPti], I ≤ i ≤ M ;用函数
【文档编号】H04B10/07GK103618568SQ201310585108
【公开日】2014年3月5日 申请日期:2013年11月16日 优先权日:2013年11月16日
【发明者】陈纯毅, 杨华民, 蒋振刚, 李华, 范静涛, 娄岩 申请人:长春理工大学