数字信号传输性能分析仪设计

1、简易数字信号传输性能分析仪摘要：本文以FPGA为核心，C8051f330作为控制部分，设计了一种数字信号传输性能分析装置。它由数字信号发生器、伪随机信号发生器、模拟低通滤波器、加法电路以及数字信号分析模块等构成。数字信号发生模块和伪随机信号发生模块由FPGA产生，伪随机信号发生器和低通滤波器模拟传输信道，数字信号分析模块由滤波和整形电路构成。，此外本系统还能实现COD等功能，经测试，各项指标均能满足设计要求。信号经过四阶巴特沃斯低通滤波后，满足滤波器带外衰减不少于40dB/十倍频程的要求，两者通过反向加法电路，相互叠加后输入数字信号分析电路，此电路中有滤波，整形，产生方波后输入到FPGA中测量

2、其脉宽，以便提取同步信号，显示出它的眼图，还能产生曼彻斯特编码，提取其信号的同步时钟，与之产生眼图，达到了数字信号传输分析的功能及其各个指标。关键字：m序列数字信号 低通滤波器 同步信号 眼图一、方案设计与比较1.1 数字信号发生器的设计： 方案一：使用中规模集成电路芯片实现指定数字信号及伪随机信号的发生，然后用门电路处理后输出信号；方案二：使用大规模集成电路（如FPGA）中的移位寄存器及其门阵列实现指定信号的发生，直接输出TTL电平；方案一成本较低，但使用的芯片数量较多，硬件设计较复杂，而且进阶要求中曼彻斯特编码会增加硬件设计的负担；方案二使用现成的FPGA开发板，可较方便且较快实现指定信号

3、的设计输出。而且可较容易实现对信号的曼彻斯特编码。所以采用方案二。1.2 低通滤波器的设计：方案一：搭建RC无源滤波器。电路比较简单，可靠性高；但通带内的信号有能量损耗，负载效应比较明显。对于500K的截止频率，额外增加噪声，不适合高频滤波，给予排除。方案二：四阶巴特沃斯低通滤波器，通带内没有能量损耗，负载效应不明显，多级级联相互影响小。四阶巴特沃斯滤波满足带外衰减大于40dB/十倍频程，通带平坦，满足设计要求。较方案一，方案二能量损耗小，负载效应不明显，切合题目要求，所以选择方案二。1.3 同步信号提取的设计方案一：先测曼切斯特信号的脉宽，得到较大波的脉宽，根据产生曼切斯特信号的原理可得，该

4、脉宽恰是其原始时钟的周期长度，根据此得到原始时钟的频率，进而得到同步信号。方案二：原序列一个循环周期中曼切斯特上升沿和下降沿的总和是固定的，在同一段时间内，可以把上升沿和下降沿的总数跟频率对应起来，即可以得到原来同步时钟的频率，频率输出由计数器实现，通过对计数器的异步清零与原序列进行相位同步，这样就可以得到同步信号。第二个方案测量的频率简单，易行，频率比原始的多或少，就在此基础上减掉几个或增加几个，对于不稳的频率或有过冲，更加准确，所以选择方案二。二、系统总体设计通过分析知，本检测装置主要有数字信号、伪随机信号发生模块、低通滤波模块、反相加法模块、数字分析模块。系统整体框图如下图一所示。图1

5、系统总体框图数字信号和伪随机信号的m序列都在FPGA中产生，数字信号经过四阶巴特沃斯滤波后，相应频带的信号与伪随机信号叠加，进入数字信号分析电路，滤去噪声，整形成方波，输入到FPGA中，测量其脉宽，得到同步信号，曼彻斯特编码产生的信号就能与提取的同步信号产生眼图了。三、理论分析与计算3.1 低通滤波器的设计巴特沃思滤波器具有通带最大平坦幅度特性，式（3-1）是n阶巴特沃思低通滤波器的幅频响应表达式 。 （3-1）截止角频率 : ; 截止频率：；电压增益：;由于设计要求每个滤波器带外衰减不少于40dB/十倍频程，采用四阶电路，理论可达80dB/十倍频程。四阶巴特沃思低通滤波器的传输函数

6、为；A(S)=A1(S)A2(S) （3-2） 100kHz低通滤波电路如下：图2 低通滤波电路参数计算：R1=R2=15K； C9=C2=100PF; R3=R4=6.2K ； C12=500PF;截止频率： ；电压增益： ；同样200kHz低通滤波电路参数值如下：R1=R2=6.8K； C9=C2=100PF; R3=R4=3.3K ； C12=500PF;截止频率：电压增益： ；500kHz低通滤波电路参数值如下：R1=R2=3K； C9=C2=100PF; R3=R4=1.2K ； C12=500PF;截止频率：；电压增益： ；由于设计要求滤波器通带增益在0.24.0内可调，在滤波电路

7、后接一级可调放大电路以满足设计需要。电压放大倍数为010倍内可调。3.2 m序列数字信号m序列数字信号由线性移位寄存器产生，如图3所示。主要由移位寄存器和反馈函数构成。反馈函数的输入端通过系数与移位寄存器的各级状态相连，通过反馈线作为x1的输入。移位寄存器在时钟的作用下把反馈函数的输出存入x1，在下一个时钟周期又把新的反馈函数的输出存入x1而把原x1的内容移入x2，依次循环下去,xn不断输出。由题要求得，m序列数字信号由线性移位寄存器产生，则为，, 的模2加n。 （3-3） （3-4）图3 m序列数字信号产生过程框图数字信号V1：f1(x) =的m 序列,由式子可用右移8位的寄存器，再由4输入

8、异或门，或门及8输入与门实现。伪随机信号V3：的m序列由式子可用右移12位的寄存器，再由4输入异或门，或门及12输入与门实现。3.3 同步信号提取首先要把模拟信号变成数字信号，滤掉高频噪声信号，由于曼彻斯特编码与原序列的关系是曼彻斯特码的上升沿表示原序列的1，下降沿表示原序列的0，可以由原序列与同步时钟异或得到，那么原序列一个循环周期中曼切斯特上升沿和下降沿的总和是固定的，在同一段时间内，可以把上升沿和下降沿的总数跟频率对应起来，即可以得到原来同步时钟的频率，频率输出由计数器实现，通过对计数器的异步清零与原序列进行相位同步，这样就可以得到同步信号。3.4 眼图显示方法观测眼图可以对信号质量作出

9、定量和定性分析，眼图的眼张开越大越好,线条越清晰越好,眼图稳定比抖动好,单线比多线好,希望交叉点在幅值的中间,图形对称。只要系统有同步信号，就可以把同步时钟和信号输到示波器中，由同步时钟外部触发显示眼图。在发挥部分中，输入端同步信号不可用，则需要从曼彻斯特码中提取同步信号时钟，就可以显示眼图。四、单元电路分析与设计41 低通滤波设计由于设计要求每个滤波器带外衰减不少于40dB/十倍频程，采用四阶巴特沃斯电路。当截止频率为100KHZ时，电路如下图所示：图4 100KHZ低通滤波电路截止频率为200KHZ时，电路如下图：图5 200KHZ低通滤波电路截止频率为500khz时，电路如下：图6 50

10、0KHZ低通滤波电路4.2 反向加法电路数字信号与伪随机信号两者相互叠加，因滤波电路使输出反向，于是采用反向求和电路。伪随机信号输入的频率达到10M，使用270M高速宽频的THS4001作为运放。图7反向加法电路4.3 数字分析电路电路 4.3.1 滤波电路 为了更好的显示波形，频率测量要准确，产生的方波要能准确表达信号的频率，要滤掉加入的噪声，因为输入数字信号最大为100KHZ，于是使用前面截止频率为500KHZ的四阶巴特沃斯电路。 4.3.2 整形电路 经过滤波隔直后的信号有正有负，使用LM311作为比较器，进行过零比较，为使输出的信号能接入FPGA中进行同步信号的提取，输出端直接接3.3

11、V电源，以便输给FPGA。图8 整形电路五、系统程序设计5.1 系统软件总体流程图图9 软件流程图信号发生部分：单片机发出控制频率的信号，改变步进来控制数字信号的数据率，数据率由计数器对50M系统时钟的分频得到，然后通过线性移位寄存器得到一个m序列，通过单片机的控制信号切换m序列输出和曼彻斯特编码输出这两种模式。同步信号提取部分：叠加信号输入后进行滤波，滤去噪声，整形成方波，输入到FPGA中，测量其频率，频率输出由计数器实现，通过对计数器的异步清零与原序列进行相位同步，这样就可以得到同步信号。六、系统测试6.1测试仪器 GOS-2064型数字示波器GOS-620型模拟示波器LPS-305电源F

12、10A型数字合成函数信号发生器6.2测试数据6.2.1数据率测试测试方法：由单片机控制FPGA输出伪随机信号的数据率，改变数据率。FPGA输出端口通过探针与示波器相连,比较设定值与实测值误差。表一 数据率测试表设定信号频率（Hz）实际信号频率（Hz）误差%1000010000.20.0022000020000.30.00153000029940.60.1984000040000.70.001755000050000.80.00166000060242.10.40357000070423.70.8000079366.60.791759000089287.30.1000001000020.002测

13、试结果分析：设定频率与实测频率误差最大为0.%，满足误差绝对值不大于1%的要求。6.2.2 通带增益测试 测试方法：先选定一个数据传输通道，输入3V的正弦波，改变放大器部分的滑动变阻器电阻，进而改变放大倍数，测量输出电压。表二为截止频率为100KHZ的增益数据,表三为截止频率为200KHZ的增益数据，表四为截止频率为500KHZ的增益数据。表二 通带增益测试表一输入频率(HZ)Au1Au2Au3Au4Au51000.1880.6866673.1246.5333331K0.1880.6866673.1246.53333310K0.1866670.6866673.1246.53333390K0.1

14、373330.482.2933332.9333334.626667100K0.1240.4366672.082.6666674.2表三 通带增益测试表二输入频率(HZ)Au1Au2Au3Au4Au51K0.20.6533331.643.2266676.53333310K0.20.6533331.643.2266676.533333100K0.1786670.621.63.0933336.333333180K0.1426670.4933331.2533332.45.133333200K0.1240.4266671.0666672.164.3表四 通带增益测试表三输入频率(HZ)Au1Au2Au3

15、Au4Au510K0.1853331.1133332.2666674.0333336.533333100K0.1853331.1066672.2666674.0333336.533333400K0.1546670.9266671.9466673.3333335.173333450K0.1453330.8666671.683.3166674.666667500K0.1293330.781.442.7066674.366667测试结果分析：在这三个截止频率的通带增益里，最小为0.129333，最大为6.533333，满足增益在0.24.9范围。6.2.3 滤波器测试测试方法：先选定一个数据传输通道

16、，输入3V的正弦波，改变放大器的输入频率，测量输出电压。表五 100KHZ的频率特性表输入频率（HZ）输出电压（V）20log(Au)1003.0002K2.96-0.0911850K2.70-0.7156680K2.34-1.6876790K2.16-2.23135100K1.94-2.961110K1.70-3.85801120K1.44-4.98546140K0.928-7.96984150K0.752-9.39827图10 100KHZ波特图表六 200KHZ的频率特性表输入频率（HZ）输出电压（V）20log(Au)1K3.000100K2.86-0.32462130K2.72-0.

17、66553150K2.60-0.97201180K2.24-1.98433190K2.10-2.4227200K1.96-2.89133210K1.80-3.46977220K1.64-4.10208250K1.18-6.33803图11 200KHZ波特图表七 500KHZ的频率特性表输入频率（HZ）输出电压（V）20log(Au)1K3.000100K2.98-0.04543400K2.42-1.45933450K2.26-1.92395470K2.18-2.16875480K2.12-2.35832490K2.06-2.55333500K2.02-2.68652550K1.74-3.70

18、004600K1.44-4.98546图12 500KHZ波特图测试结果分析：由波特图可知，第一种的滤波器截止频率为95KHZ，第二种的滤波器截止频率为190KHZ，第三种的滤波器截止频率为480KHZ，误差都小于10%，满足设计要求。6.2.4 伪随机信号数据率、峰峰值测试方法：FPGA产生10M数据率，用示波器探头去测量其输出率，再用探头去测量伪随机的输出信号的峰峰值。测试数据：数据率：10.0002MHZ 峰峰值106.8mV测试结果分析：数据率误差绝对值百分比为(10.0002-10)/10*100=0.002% 远小于1%的绝对误差，峰峰值误差为6.8%，小于10%，都满足设计要求。

19、6.2.5 信号眼图显示测试方法：FPGA产生两路伪随机信号，一路信号通过低通滤波器，增益放大电路后与另一路伪随机信号相加作为输出信号。FPGA产生的同步时钟与输出信号同时输出给模拟示波器，观察图形。表八 眼幅度测量表V2（V）0.51.52.833.263.965.968.9610.0012.0013.00眼幅度(V)0.10.270.580.60.751.21.822.42.6下图为电压为8.96V与100mV噪声叠加后的眼图图13 眼图6.2.6 伪随机信号幅度测量伪随机信号峰峰值的范围24.8mV3.14V6.2.7 低信噪比的眼图图14 低噪比的眼图七、结论根据以上设计方案，完成了数

20、字信号传输性能分析仪。它由数字信号发生电路、低通滤波电路、求和电路、数字信号分析电路构成，能产生数字信号和伪随机信号，峰峰值都能到达TTL电平，每个滤波的截止频率误差都小于10%，增益在0.24.0范围内可调，在很低信噪比下能提取产生曼彻斯特编码信号的同步时钟，能与同步信号产生眼图，指标性能较好，完成了设计要求。八、参考文献1、童诗白、华成英 模拟电子技术基础（第四版） 高等教育出版社 2006，052、马忠梅 单片机的C语言应用程序设计 北京航空航天大学出版社 2007，02九、附录附录一 主要程序清单将10MHz时钟分频为100KHz时钟process(clkIn)beginif clkIn'event and clkIn = '1' thencounter <= counter + '1'if counter = "110001" thenclkTmp <= not clkTmp;counter <= "000000"end if;clkOut <= clkTmp;end if;end process;测脉宽process(PulseIn,Gate)beginif Gate = '0' thenCo