基于DSP的全数字PGC解调技术研究..doc

基于DSP的全数字PGC 解调技术研究 摘要：基于DSP TMS320VC5509a 平台实现的数字解调算法，利用查表法产生载波，克服了传统电路生成载波与A/D采样得到的干涉信号混频时存在的延迟问题，充分利用了DSP并行处理信号的优势。 实现了在低频（20Hz-1.1KHz）范围内的PGC解调，可用于地震波检测及石油井下勘测及长距离分布式测量。 引言 相位生成载波（phase generated carrier，PGC）[1]调制解调技术是一种广泛应用于干涉型光纤传感器的零差解调技术，可以消除相位漂移引起的干涉输出不稳，具有灵敏度高、动态范围大、线性度好等优点，被广泛应用于光纤水听器、光纤加速度计、光纤振动传感器等光纤传感器的信号解调中。随着数字信号处理技术的发展，全数字PGC解调技术已经成为当前研究的重点。 在PGC解调过程中，一般将干涉输出直接与基频、二倍频信号进行混频相乘，然后进行后续处理。然而，由于载波产生方式和光纤传输等原因，干涉输出中的载波与混频相乘时的基频、二倍频信号之间存在一定延迟，会对PGC解调产生不良影响，严重可导致解调失败[2]。 本文中载波在DSP解调系统中通过查表法产生，查表法产生基频、二倍频数据存放在DSP内存中，通过占用系统的内存空间可以克服干涉信号与载频信号不同步问题，即以空间换取时间的同步。而当前一般通用的54x、55x 系列DSP芯片内存是可以满足这样的要求的。 PGC 调制解调原理 PGC调制分为内调制和外调制，干涉仪结构分为Mach-Zehnder、Michelson、Sagnac等不同结构，但无论采取哪种调制，哪种干涉结构，干涉仪输出光信号均可表示为： （1） 其中，A为直流量、B为干涉项幅度，C为载波幅值，也称为相位调制度，为载波角频率， ，是待测信号，D为被测信号幅值，为被测信号角频率，为相位漂移。 式（1）按Bessel函数展开得： （2） 同样可得和的Bessel函数展开： （3） （4） 图1为采用DCM算法的PGC解调原理框图，图1中，干涉信号与基频和二倍频混频相乘后经过低通滤波后得到： 图1 基于DCM 算法的PGC解调原理框图 （5） （6） 经过微分交叉相乘相减之后得到： （7） 对上式积分得： （8） 再进行高通得到： （9） 即被测信号被完整的解调出来，只是幅值上有一些变化，同时消除了随机相位漂移的影响。 2.相关参数的确定 由式（2）、（3）、（4）可知，干涉信号I的频谱在理论上是无限带宽，信号和也是无限带宽，理论上需要的载波频率无穷大，低通滤波器的截止频率无穷大，这都是无法实现的。由于信号调制只是进行功率的重新分配，即调制前后总功率不变，只是将为调制信号功率的一部分分配给各个边频分量，而分配的原则与信号的幅度即调制深度有关。尽管干涉信号I的频谱是无限宽，但根据Bessel函数的衰减性质，如图2所示。I的频谱的主要成分是集中在一定范围内的。 工程上常采用的原则是：当边频信号的功率小于未调制信号功率的1%，即边频信号的幅度小于未调制信号时幅度10%可以忽略不计[]。 2.1 调制深度C的选择 由DCM算法解调输出知，增益与有关；随着环境的变化，伴生调幅等因素影响，C也会有波动，导致增益的变化，使得原先定标在一定程度上失效，系统解调的时间稳定性及可靠性下降。 C的选取，应遵循以下3条原则：其一，尽可能大。其二，导数的绝对值尽可能小。最后，C尽可能小。其中，前两条保证取得极值。 对于DCM算法，C=2.37为最佳值[]。 Bessel函数的衰减特性如图2所示，调制度C对1、2阶Bessel函数及其导数的影响如如图3所示： 图2.Bessel函数的衰减特性 图3.C对 及其导数的影响 2.2 采样频率的确定 由式（2）、（3）、（4）知 信号I的频谱成分中含有载波的基频以及它们的各次谐波成分，在每一个的两侧携带有的边频信号，k为（2）式中的Bessel函数阶次，n 为式（3）、（4）中的贝塞尔函数阶次。