双向大气湍流光信道的信号起伏相关系数测量系统及方法与流程

本发明属于大气信道光信号传输与通信技术领域，涉及一种双向大气湍流光信道的信号起伏相关系数测量系统及方法。

背景技术：

对于双向大气湍流光信道来说，如果两个方向上的光信号随机起伏具有良好的相关性，则收发端机就可以利用这种相关性直接获得瞬时信道状态信息(无需从接收端向发射端反馈瞬时信道状态信息)，从而实现发射自适应技术(后面称为基于发端获取信道状态信息的发射自适应技术)。具体地说，如果收发端机a和收发端机b之间有一条双向传输信道，收发端机a通过测量从收发端机b发来的光信号的瞬时起伏就可以获得从收发端机a发送到收发端机b的光信号瞬时起伏信息。该方法能有效工作的前提是，从收发端机a发给收发端机b的光信号瞬时起伏与从收发端机b发给收发端机a的光信号瞬时起伏之间具有良好的相关性。实际测量双向大气湍流光信道的两个传输方向上的光信号起伏相关系数对研究和分析基于发端获取信道状态信息的发射自适应技术的适用性具有重要意义。在测量双向大气湍流光信道的两个传输方向上的光信号起伏相关系数时，可以同时在收发端机a和收发端机b记录接收到的光信号瞬时强度采样值，然后再对收发端机a和收发端机b记录的光信号瞬时强度采样值进行统计相关分析，以得到相关系数。然而，如果测量系统使用的激光器不是价格高昂的高端激光器，激光器的输出光强可能具有较大波动。收发端机a和收发端机b的激光器的输出光强波动通常不相关，这会导致双向大气湍流光信道的两个传输方向上的光信号起伏相关系数的测量结果包含激光器输出光强波动的影响，不利于准确研究和分析大气湍流造成的两个传输方向上的光信号起伏相关性。因此，在测量双向大气湍流光信道的两个传输方向上的光信号起伏相关系数时，在使用廉价激光器的条件下排除激光器输出光强波动影响非常必要。本发明提供一种双向大气湍流光信道的信号起伏相关系数测量系统及方法，可以排除激光器输出光强波动对双向大气湍流光信道的两个传输方向上的光信号起伏相关系数的测量结果的影响。

技术实现要素：

本发明的目的在于，提供一种双向大气湍流光信道的信号起伏相关系数测量系统及方法，其可以排除激光器输出光强波动对测量结果的影响。

本方法的技术方案是这样实现的：一种双向大气湍流光信道的信号起伏相关系数测量系统及方法，其特征在于，所需的硬件系统和执行步骤如下：

如图1所示，硬件系统包括收发端机a和收发端机b，收发端机a包括激光器a(101)、光纤分路器a(102)、光电探测器a1(103)、光纤环形器a(104)、光电探测器a2(105)、计算机a(106)和光学收发系统a(107)，收发端机b包括激光器b(201)、光纤分路器b(202)、光电探测器b1(203)、光纤环形器b(204)、光电探测器b2(205)、计算机b(206)和光学收发系统b(207)。

激光器a(101)发出的激光信号c001被耦合进入光纤分路器a(102)的端口fsa1并分成两路分别从光纤分路器a(102)的端口fsa2和光纤分路器a(102)的端口fsa3输出，从光纤分路器a(102)的端口fsa2输出的激光信号c001入射到光电探测器a1(103)上，从光纤分路器a(102)的端口fsa3输出的激光信号c001入射进入光纤环形器a(104)的端口fca1并从光纤环形器a(104)的端口fca2输出进入光学收发系统a(107)；光学收发系统a(107)把通过光纤环形器a(104)的端口fca2输入的激光信号c001发射到大气湍流信道中。

激光器b(201)发出的激光信号d001被耦合进入光纤分路器b(202)的端口fsb1并分成两路分别从光纤分路器b(202)的端口fsb2和光纤分路器b(202)的端口fsb3输出，从光纤分路器b(202)的端口fsb2输出的激光信号d001入射到光电探测器b1(203)上，从光纤分路器b(202)的端口fsb3输出的激光信号d001入射进入光纤环形器b(204)的端口fcb1并从光纤环形器b(204)的端口fcb2输出进入光学收发系统b(207)；光学收发系统b(207)把通过光纤环形器b(204)的端口fcb2输入的激光信号d001发射到大气湍流信道中。

光学收发系统a(107)收到的来自收发端机b的激光信号d001进入光纤环形器a(104)的端口fca2并通过光纤环形器a(104)的端口fca3入射到光电探测器a2(105)上，光电探测器a2(105)输出的电信号传送到计算机a(106)的信号采集卡1上。光电探测器a1(103)输出的电信号传送到计算机a(106)的信号采集卡2上。

光学收发系统b(207)收到的来自收发端机a的激光信号c001进入光纤环形器b(204)的端口fcb2并通过光纤环形器b(204)的端口fcb3入射到光电探测器b2(205)上，光电探测器b2(205)输出的电信号传送到计算机b(206)的信号采集卡1上。光电探测器b1(203)输出的电信号传送到计算机b(206)的信号采集卡2上。

1)本方法的第一部分使收发端机a和收发端机b正常工作，具体操作包括：

使收发端机a的光学收发系统a(107)和收发端机b的光学收发系统b(207)相互对准；使收发端机a和收发端机b包含的各个组成部分都正常工作。

2)本方法的第二部分在收发端机a中执行采样测量操作，具体步骤如下：

步骤step201：在计算机a(106)的内存中创建一个列表lista1，列表lista1的元素用于存储计算机a(106)的信号采集卡1采集到的入射激光信号强度采样测量值；在计算机a(106)的内存中创建一个列表lista2，列表lista2的元素用于存储计算机a(106)的信号采集卡2采集到的入射激光信号强度采样测量值；令列表lista1为空；令列表lista2为空；令ta＝0；令ia＝1；

步骤step202：在时刻ta，计算机a(106)通过信号采集卡1获得一个激光信号强度采样测量值e001，计算机a(106)通过信号采集卡2获得一个激光信号强度采样测量值e002；在列表lista1的末尾添加一个新元素ea1，新元素ea1是列表lista1的第ia个元素，把列表lista1的第ia个元素赋值为采样测量值e001的值；在列表lista2的末尾添加一个新元素ea2，新元素ea2是列表lista2的第ia个元素，把列表lista2的第ia个元素赋值为采样测量值e002的值；令ia＝ia+1；令ta＝ta+δs；δs表示采样时间间隔；

步骤step203：如果ia≤num，num表示采集的采样测量值总数，则转步骤step202，否则转步骤step204；

步骤step204：收发端机a的采样测量操作结束。

3)本方法的第三部分在收发端机b中执行采样测量操作，具体步骤如下：

步骤step301：在计算机b(206)的内存中创建一个列表listb1，列表listb1的元素用于存储计算机b(206)的信号采集卡1采集到的入射激光信号强度采样测量数据；在计算机b(206)的内存中创建一个列表listb2，列表listb2的元素用于存储计算机b(206)的信号采集卡2采集到的入射激光信号强度采样测量数据；令列表listb1为空；令列表listb2为空；令tb＝0；令ib＝1；

步骤step302：在时刻tb，计算机b(206)通过信号采集卡1获得一个激光信号强度采样测量值f001，计算机b(206)通过信号采集卡2获得一个激光信号强度采样测量值f002；在列表listb1的末尾添加一个新元素eb1，新元素eb1是列表listb1的第ib个元素，把列表listb1的第ib个元素赋值为采样测量值f001的值；在列表listb2的末尾添加一个新元素eb2，新元素eb2是列表listb2的第ib个元素，把列表listb2的第ib个元素赋值为采样测量值f002的值；令ib＝ib+1；令tb＝tb+δs；δs表示采样时间间隔；

步骤step303：如果ib≤num，num表示采集的采样测量值总数，则转步骤step302，否则转步骤step304；

步骤step304：收发端机b的采样测量操作结束。

4)本方法的第四部分计算从收发端机a到收发端机b的信道的光信号起伏与从收发端机b到收发端机a的信道的光信号起伏的相关系数，具体步骤如下：

步骤step401：分别针对i＝1,2,…,num，令列表lista1的第i个元素的值等于列表lista1的第i个元素的值除以列表lista2的第i个元素的值的结果；分别针对i＝1,2,…,num，令列表listb1的第i个元素的值等于列表listb1的第i个元素的值除以列表listb2的第i个元素的值的结果；

步骤step402：从收发端机a到收发端机b的信道的光信号起伏与从收发端机b到收发端机a的信道的光信号起伏的相关系数μ为：

其中sa[i]表示列表lista1的第i个元素的值，sb[i]表示列表listb1的第i个元素的值，表示列表lista1的所有元素的值的平均值，表示列表listb1的所有元素的值的平均值。

在实施本方法时，首先执行本方法的第一部分，然后同时开始执行本方法的第二部分和第三部分，最后执行本方法的第四部分。

本发明的积极效果是，通过同时在发射端测量激光器输出光强波动和在接收端测量经大气湍流信道传输后的激光信号起伏，可以在计算从收发端机a到收发端机b的信道的光信号起伏与从收发端机b到收发端机a的信道的光信号起伏的相关系数时，排除激光器输出光强波动的影响，从而有助于更加准确地研究大气湍流导致的从收发端机a到收发端机b的信道的光信号起伏与从收发端机b到收发端机a的信道的光信号起伏的相关性。

附图说明

图1为双向大气湍流光信道的信号起伏相关系数测量系统示意图。

具体实施方式

为了使本方法的特征和优点更加清楚明白，下面结合具体实施例对本方法作进一步的描述。在本实施例中，收发端机a和收发端机b分别位于两栋高楼的屋顶，激光器a(101)选择波长为808nm的半导体激光器，激光器b(201)选择波长为793nm的半导体激光器，这两个波长都属于800nm波段。光电探测器a1(103)、光电探测器a2(105)、光电探测器b1(203)、光电探测器b2(205)都为pin光电探测器。

本方法的技术方案是这样实现的：一种双向大气湍流光信道的信号起伏相关系数测量系统及方法，其特征在于，所需的硬件系统和执行步骤如下：

如图1所示，硬件系统包括收发端机a和收发端机b，收发端机a包括激光器a(101)、光纤分路器a(102)、光电探测器a1(103)、光纤环形器a(104)、光电探测器a2(105)、计算机a(106)和光学收发系统a(107)，收发端机b包括激光器b(201)、光纤分路器b(202)、光电探测器b1(203)、光纤环形器b(204)、光电探测器b2(205)、计算机b(206)和光学收发系统b(207)。

激光器a(101)发出的激光信号c001被耦合进入光纤分路器a(102)的端口fsa1并分成两路分别从光纤分路器a(102)的端口fsa2和光纤分路器a(102)的端口fsa3输出，从光纤分路器a(102)的端口fsa2输出的激光信号c001入射到光电探测器a1(103)上，从光纤分路器a(102)的端口fsa3输出的激光信号c001入射进入光纤环形器a(104)的端口fca1并从光纤环形器a(104)的端口fca2输出进入光学收发系统a(107)；光学收发系统a(107)把通过光纤环形器a(104)的端口fca2输入的激光信号c001发射到大气湍流信道中。

激光器b(201)发出的激光信号d001被耦合进入光纤分路器b(202)的端口fsb1并分成两路分别从光纤分路器b(202)的端口fsb2和光纤分路器b(202)的端口fsb3输出，从光纤分路器b(202)的端口fsb2输出的激光信号d001入射到光电探测器b1(203)上，从光纤分路器b(202)的端口fsb3输出的激光信号d001入射进入光纤环形器b(204)的端口fcb1并从光纤环形器b(204)的端口fcb2输出进入光学收发系统b(207)；光学收发系统b(207)把通过光纤环形器b(204)的端口fcb2输入的激光信号d001发射到大气湍流信道中。

光学收发系统a(107)收到的来自收发端机b的激光信号d001进入光纤环形器a(104)的端口fca2并通过光纤环形器a(104)的端口fca3入射到光电探测器a2(105)上，光电探测器a2(105)输出的电信号传送到计算机a(106)的信号采集卡1上。光电探测器a1(103)输出的电信号传送到计算机a(106)的信号采集卡2上。

光学收发系统b(207)收到的来自收发端机a的激光信号c001进入光纤环形器b(204)的端口fcb2并通过光纤环形器b(204)的端口fcb3入射到光电探测器b2(205)上，光电探测器b2(205)输出的电信号传送到计算机b(206)的信号采集卡1上。光电探测器b1(203)输出的电信号传送到计算机b(206)的信号采集卡2上。

1)本方法的第一部分使收发端机a和收发端机b正常工作，具体操作包括：

使收发端机a的光学收发系统a(107)和收发端机b的光学收发系统b(207)相互对准；使收发端机a和收发端机b包含的各个组成部分都正常工作。

2)本方法的第二部分在收发端机a中执行采样测量操作，具体步骤如下：

步骤step201：在计算机a(106)的内存中创建一个列表lista1，列表lista1的元素用于存储计算机a(106)的信号采集卡1采集到的入射激光信号强度采样测量值；在计算机a(106)的内存中创建一个列表lista2，列表lista2的元素用于存储计算机a(106)的信号采集卡2采集到的入射激光信号强度采样测量值；令列表lista1为空；令列表lista2为空；令ta＝0；令ia＝1；

步骤step202：在时刻ta，计算机a(106)通过信号采集卡1获得一个激光信号强度采样测量值e001，计算机a(106)通过信号采集卡2获得一个激光信号强度采样测量值e002；在列表lista1的末尾添加一个新元素ea1，新元素ea1是列表lista1的第ia个元素，把列表lista1的第ia个元素赋值为采样测量值e001的值；在列表lista2的末尾添加一个新元素ea2，新元素ea2是列表lista2的第ia个元素，把列表lista2的第ia个元素赋值为采样测量值e002的值；令ia＝ia+1；令ta＝ta+δs；δs表示采样时间间隔；

步骤step203：如果ia≤num，num表示采集的采样测量值总数，则转步骤step202，否则转步骤step204；

步骤step204：收发端机a的采样测量操作结束。

3)本方法的第三部分在收发端机b中执行采样测量操作，具体步骤如下：

步骤step301：在计算机b(206)的内存中创建一个列表listb1，列表listb1的元素用于存储计算机b(206)的信号采集卡1采集到的入射激光信号强度采样测量数据；在计算机b(206)的内存中创建一个列表listb2，列表listb2的元素用于存储计算机b(206)的信号采集卡2采集到的入射激光信号强度采样测量数据；令列表listb1为空；令列表listb2为空；令tb＝0；令ib＝1；

步骤step302：在时刻tb，计算机b(206)通过信号采集卡1获得一个激光信号强度采样测量值f001，计算机b(206)通过信号采集卡2获得一个激光信号强度采样测量值f002；在列表listb1的末尾添加一个新元素eb1，新元素eb1是列表listb1的第ib个元素，把列表listb1的第ib个元素赋值为采样测量值f001的值；在列表listb2的末尾添加一个新元素eb2，新元素eb2是列表listb2的第ib个元素，把列表listb2的第ib个元素赋值为采样测量值f002的值；令ib＝ib+1；令tb＝tb+δs；δs表示采样时间间隔；

步骤step303：如果ib≤num，num表示采集的采样测量值总数，则转步骤step302，否则转步骤step304；

步骤step304：收发端机b的采样测量操作结束。

4)本方法的第四部分计算从收发端机a到收发端机b的信道的光信号起伏与从收发端机b到收发端机a的信道的光信号起伏的相关系数，具体步骤如下：

步骤step401：分别针对i＝1,2,…,num，令列表lista1的第i个元素的值等于列表lista1的第i个元素的值除以列表lista2的第i个元素的值的结果；分别针对i＝1,2,…,num，令列表listb1的第i个元素的值等于列表listb1的第i个元素的值除以列表listb2的第i个元素的值的结果；

步骤step402：从收发端机a到收发端机b的信道的光信号起伏与从收发端机b到收发端机a的信道的光信号起伏的相关系数μ为：

其中sa[i]表示列表lista1的第i个元素的值，sb[i]表示列表listb1的第i个元素的值，表示列表lista1的所有元素的值的平均值，表示列表listb1的所有元素的值的平均值。

在实施本方法时，首先执行本方法的第一部分，然后同时开始执行本方法的第二部分和第三部分，最后执行本方法的第四部分。

在本实施例中，num＝20000，δs＝1毫秒。