一种基于统计特性分析的IMDD‑OFDM‑PON定点FFT字长优化的方法与流程

本发明涉及一种ofdm-pon(opticalorthogonalfrequencydivisionmultiplexing-passiveopticalnetwork正交频分复用无源光网络)的定点fft字长优化方法。

背景技术：

光正交频分复用(oofdm)技术是一种新的光传输技术，其中快速傅里叶变换(fft，fastfouriertransform)是oofdm系统中的核心技术之一。光ofdm系统收发机的速率已经达到了几十gbits/s，这就要求系统中的fft处理单元具有非常高的吞吐率。对于oofdm系统，其dsp算法所占用的逻辑资源，尤其是其中fft的计算已经占据了现有fpga资源的80％左右。因此，fft的设计也成为了实现实时系统的巨大的挑战之一。

为了满足高速实时系统的要求，以及实现高吞吐率的ofdm系统，fft需要采用并行流水线结构，同时为了降低fft算法实现所需要的资源，现有fft的实现都采用定点的fft，其资源优化主要集中于定点fft所需要的字长优化。对于fft的字长优化，国内外学者提出了不同的研究策略以及最优字长的搜索算法。文献1[bouzianer,koutsoyannisr,milderp,etal.optimizingfftprecisioninopticalofdmtransceivers[j].ieeephotonicstechnologyletters,2011,23(20):1550-1552.]主要分析了针对不同ad/da分辨率下，fft输入输出字长对系统evm性能的影响，而未对中间级字长进行分析。文献2[zhangjj,yuanwy,giddingsrp,etal.ifftstage-dependentminimumbitresolutionmapsforreal-timeopticalofdmtransceivers[c].opticalfibercommunicationsconferenceandexhibition(ofc),2014:1-3.]对基于dit-2结构的32点fft/ifft进行了联合优化设计，但未提及字长优化策略。文献3[wangcy,kuocb,joujy.hybridwordlengthoptimizationmethodsofpipelinedfftprocessors[j].ieeetransactionsoncomputers,2007,56(8):1105-1118.]利用统计分析提出了搜索fft中间级字长的方案，但是搜索复杂度偏高。

技术实现要素：

本发明的目的在于针对已有fft字长优化方案存在的不足之处，提供了一种ofdm-pon的定点fft字长优化方法，能够实现基于radix-2(基2)结构的定长fft中间级输出字长以及旋转因子字长优化，并将fft最优字长的搜索复杂度降至o(1)。

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

一种ofdm-pon定点fft字长优化的方法，fft的点数设置为n，采用radix-2dit(decimation-in-time，时域抽取)结构的fft共有log2n级，每级定点数的输出采用minf，所述m为定点数的整数部分字长，i表示整数，n为定点数的小数部分字长，f表示小数，每级fft的输出字长为m+n，每级m的字长比前一级的字长增加1位，fft的输入字长为强度调制直接检测ofdm-pon系统中adc的量化字长q，所述定点fft字长优化包括中间级输出字长优化以及参与fft蝶形计算的旋转因子字长优化，所述中间级输出字长以及旋转因子字长都独立可调。

具体操作步骤如下：

1)ofdm-pon接收端对ofdm-pon信号进行数据采集和离线分析，用浮点fft对ofdm-pon信号进行解调并计算ofdm-pon系统的理想逆误差向量幅度ievmideal_db，其单位采用db；

2)预设系统可接受的定点fft计算的逆误差向量幅度ievmtotal_db，其单位采用db，所述ievmtotal_db<ievmideal_db；

3)进行旋转因子字长优化，每级旋转因子的字长都采用同样的数值，并把fft每级输出设置为浮点数，并采用旋转因子字长递增的方法，计算ofdm-pon系统的旋转因子逆误差向量幅度ievmtwiddle_db，其单位采用db来表示，当ievmtwiddle_db≥(1+α)×ievmtotal_db时，获取fft计算所需要的旋转因子字长ltwiddle，所述α为预设的由旋转因子引起的逆误差向量幅度调整因子；

4)进行中间级输出字长优化，所述中间级输出字长优化包括以下步骤：

4-1)根据预设误差ievmtotal_db，并设定fft中间级输出字长对ofdm-pon系统逆误差向量幅度的影响为加性高斯非相关，得出每级输出字长逆误差向量幅度应满足的指标：所述ofdm-pon为强度调制直接检测方式时，fft中间级每级输出字长逆误差向量幅度所述

4-2)推出每级输出字长逆误差向量幅度与字长的关系：ievmoutput_db(k)＝6.02×loutput(k)-27.7-(k-3)×3+β,k∈[3,log2n]，所述loutput(k)为fft中间级k级输出字长，k为中间级索引号，β为对应的系统修正因子，所述系统修正因子与fft所采用的点数以及fft发送端的限幅参数相关。

4-3)第一级与第二级的输出字长设置为q+1与q+2，中间级第k级的输出字长k∈[3,log2n]，所述ceil为向上取整。

本发明与现有技术相比，具有如下显而易见的突出实质性特点和显著的进步：

1.本发明提出的radix-2dit的fft字长优化方法不依赖于fft点数，使用范围广；

2.本发明提出的fft字行优化方法不依赖于强度调制直接检测光链路的性能，能够根据物理链路进行自适应的调整，具有较强的鲁棒性；

3.本发明提出的fft优化方法能够对fft中间级的旋转因子字长以及输出字长进行联合优化，极大地降低了ofdm系统的逻辑资源开销；

4.本发明提出的fft优化方法中fft中间级的字长输出的算法复杂度为o(1)，与其他fft字长搜索算法相比，算法复杂度低。

附图说明

图1是本发明的ofdm-pon定点fft字长优化方法的程序框图。

图2是本发明中基于强度调制直接检测的ofdm-pon仿真平台流程图。

图3是本发明实施例中ofdm-pon物理帧。

图4是本发明实施例中光链路vpi仿真结构图。

图5是本发明实施例中1024定点fft中间级整数部分增长图。

图6是本发明实施例中1024定点fft旋转因子字长搜索结果图。

图7是本发明实施例中1024定点fft每级输出字长映射图。

具体实施方式

本发明的优选实施例结合附图说明如下：

实施例一：

1.参见图1，本ofdm-pon定点fft字长优化方法，fft的点数设置为n，采用radix-2dit结构的fft共有log2n级，每级定点数的输出采用minf，所述m为定点数的整数部分字长，i表示整数，n为定点数的小数部分字长，f表示小数，每级fft的输出字长为m+n，每级m的字长比前一级的字长增加1位，fft的输入字长为强度调制直接检测ofdm-pon系统中adc的量化字长q；所述定点fft字长优化包括中间级输出字长优化以及参与fft蝶形计算的旋转因子字长优化，所述中间级输出字长以及旋转因子字长都独立可调；具体操作步骤如下：

1)ofdm-pon接收端对ofdm-pon信号进行数据采集和离线分析，用浮点fft对ofdm-pon信号进行解调并计算ofdm-pon系统的理想逆误差向量幅度ievmideal_db，其单位采用db；

2)预设系统可接受的定点fft计算的逆误差向量幅度ievmtotal_db，其单位采用db，所述ievmtotal_db<ievmideal_db；

3)进行旋转因子字长优化，每级旋转因子的字长都采用同样的数值，并把fft每级输出设置为浮点数，并采用旋转因子字长递增的方法，计算ofdm-pon系统的旋转因子逆误差向量幅度ievmtwiddle_db，其单位采用db来表示，当ievmtwiddle_db≥(1+α)×ievmtotal_db时，获取fft计算所需要的旋转因子字长ltwiddle，所述α为预设的由旋转因子引起的逆误差向量幅度调整因子；

4)进行中间级输出字长优化。

实施例二：

本实施例与实施例一基本相同，特别之处如下：

所述步骤4)进行中间级输出字长优化包括以下步骤：

4-1)根据预设误差ievmtotal_db，并设定fft中间级输出字长对ofdm-pon系统逆误差向量幅度的影响为加性高斯非相关，得出每级输出字长逆误差向量幅度应满足的指标：所述ofdm-pon为强度调制直接检测方式时，fft中间级每级输出字长逆误差向量幅度所述

4-2)推出每级输出字长逆误差向量幅度与字长的关系：ievmoutput_db(k)＝6.02×loutput(k)-27.7-(k-3)×3+β,k∈[3,log2n]，所述loutput(k)为fft中间级k级输出字长，k为中间级索引号，β为对应的系统修正因子，所述系统修正因子与fft所采用的点数以及fft发送端的限幅参数相关；

4-3)第一级与第二级的输出字长设置为q+1与q+2，中间级第k级的输出字长k∈[3,log2n]，所述ceil为向上取整。

实施例三：

为了研究变长fft的字长优化方案，本实施例参见图2-图7所示ofdm系统的收发结构，同时基于vpi和matlab仿真软件搭建了imdd-oofdm-pon系统仿真平台。ofdm的收发系统的数字信号调制与解调模块通过matlab的离线分析来完成，其参数如表1所示。光纤传输信道通过vpi软件仿真完成。在该仿真平台中，结合本发明实时实例，fft/ifft的点数n设置为1024点。

表1imdd-ofdm-pon光路的参数设置表

在发送端，将15阶伪随机二进制序列(prbs)作为系统的输入，经过串并转换单元转换成511路并行的数据，接着对并行数据流进行qam映射编码。由于收发系统的逆误差向量幅度独立于系统的调制格式，因此所有511个子载波采用的相同的调制格式即32-qam，同时考虑到实际oofdm-pon系统的发送端信号采用交流输出，因此关闭ifft的0号子载波输入。接着对数据进行共轭对称操作，以保证ifft运算后的输出结果仅包含实部以便在imdd-ofdm系统进行传输，然后在发送的ofdm数据前加上16点的循环前缀以保证子载波的正交性以及克服ofdm符号间的干扰。在本文中为了分析fft字长对系统性能的影响，每一个ofdm帧共包括1000个ofdm符号，接着在发送ofdm数据前增加训练序列以实现符号同步以及信道均衡。上述所示帧结构如图3所示。

然后将上述一帧ofdm信号通过12db的限幅操作(clipping)解决ofdm系统高峰均比从而提升ofdm收发系统的传输性能。最后，通过10位数模转换器(digital-to-analogconverter,dac)将数据转换成模拟信号在信道上传输，其中dac输出的幅度峰峰值为0.04v，使得dac的输出能够通过电压与电流转换模块驱动分布式反馈激光器(dfb)。

光链路的传输系统是通过vpi软件来完成的，将模拟信号调制后通过分布式反馈激光器实现信号的电光转换，从dfb激光器输出的光功率为10dbm，转换后的光信号经25km标准单模光纤(ssmf)传输后，通过衰减器调节使其接收光功率为0dbm，接着经过响应率为1a/w的光电二极管将信号进行光电转换。其仿真结构如图4所示。

在接收端，将光电二极管转换后的电模拟信号进行10位模数转换器(analog-to-digitalconverter,adc)，也即adc的量化字长q＝10。接着根据训练序列实现系统的同步及信道估计。在去除循环前缀之后，通过fft计算将ofdm符号数据转化为qam信号，从而实现数据的解调以及系统逆误差向量幅度的计算。

在本实施例中，把adc量化后的10位有符号数转换为5i5f的定点数，该定点数也即为1024点fft的输入数据。则1024点定长fft中间级定点数的整数部分字长如图5所示，fft中间级每增加一级，则对应的整数部分字长需要增加1位。

为了对定点fft字长进行优化，需把fft设置为浮点fft，计算出上述系统的理想逆误差向量幅度ievmideal_db为30.29db。在本实施例中，由于需要对fft进行字长优化，因此在ofdm系统解调时，必然会带来系统的性能损失。在本次fft字长优化实施例中，允许损失的系统性能为0.5db,也即逆误差向量幅度ievmtotal_db＝29.79db。

1024点定点fft的字长优化包含旋转因子字长优化以及中间级输出字长优化。首先设定中间旋转因子的字长都为相同的数值，然后对中间级的字长进行优化。由于中间级旋转因子字长优化同样会带来系统性能损失，在本次实施中，首先设定α为0.015，也即ievmtwiddle_db≥30.24db时对应的旋转因子字长为满足本实施实例中的旋转因子字长映射关系。利用上述离线仿真得出中间级旋转因子字长与系统性能的映射图如图6所示，从图6可以看出当ltwiddle＝8，系统的逆误差向量幅度满足实时例设定的约束，也即所有旋转因子都设定为8位。

根据实施例中设定的ievmtwiddle_db以及ievmtotal_db，并依据ievmoutput_db的输出满足条件的公式可得出1024点fft中间字长的ievmoutput_db≥48.94db，在本实施例中系统修正因子β＝1.6db。因此根据k∈[3,log2n]，可以得出1024点fft中间级字长的映射表如图7所示。把图7得出的1024点输出字长映射表与旋转因子字长为8，用图2的imdd-ofdm-pon系统进行仿真，得出ofdm-pon接收端的逆误差向量幅度性能为29.89db。因此采用旋转因子与中间级字长映射表对fft进行优化后，其对应的ievm相比于原系统性能减少了0.4db，小于系统的性能恶化容忍程度，因此验证了本发明基于统计特性的fft优化算法的正确性与准确性。

以上实施例仅用于说明本发明的技术方案，并非用于限制本发明的范围，凡在本发明精神和原则之内所做的任何修改、等同替换和改进等，均含于本发明的保护范围之内。