微弱信号程序.doc

Harbin Institute of Technology 课程设计（论文） 课程名称： 微弱信号检测及处理 设计题目： 软正交矢量型LIA在 微弱信号检测中的应用 院 系： 电气工程及自动化学院 班 级： 设 计 者： 学 号： 指导教师： 设计时间： 2013-12-26 摘　要 本报告中利用软件移相技术得到相互正交的参考正弦波信号，通过互相关算法，完成了软正交矢量型LIA相关器的具体实现，利用该方法实现了对微弱信号幅值和相位的检测，有效地抑制了干扰，减少了硬件电路成本，并且通过Matlab仿真，验证了该算法具有一定的优越性，并在实际应用中取得了很好的效果。 关键词：正交矢量；移相；互相关；微弱信号 1 互相关原理 设混有随机噪声的信号： 其中：为有用信号且其重复周期或频率已知，为随机噪声信号； 参考信号为： ， 且和相位差为。 则有： 　　　　　　　（1） 　　　　　　　　　　（2） 　　　　　　　（3） 求解式（1）式（2）和式（3）可以得到： 　　　　　　　　　　 （4） 　　　　　　　　　　　　　 （5） 　　　　　　　　　　　　　　　（6） 2 软正交矢量型LIA的原理及实现 2.1 实现框图 实现框图如图1所示。 图 1 软正交矢量LIA实现框图 由检测线圈感应出来的信号比较微弱、干扰多,通过信号调理电路进入A /D采集卡的1通道,参考信号进入采集卡的通道2,由计算机控制对两路信号进行采样。 2.2 软正交矢量LIA的实现 一般地，实现移相的方法主要有： 直接对模拟信号进行移相，如阻容移相、变压器移相，而采用这种移相器有很多不足之处，输出波形容易受输入波形的影响，移相操作不方便，移相的角度随时间和负载等因素的影响而产生漂移，误差较大。 利用单片机技术形成的数字移相技术，即利用单片机通过2片D/A转换芯片，连续循环地读取存储在ROM中的正弦波数字化的数据表，获得两路正弦信号；相位差由读取数据表的位置不同而不同，但这需要很多的芯片和设计硬件电路。 本文采用软件移相方法，算法的主要实现过程如下： 设定参考信号的3个连续采样值为ik-1，ik，ik+1，信号角频率为ω，采样周期为Ts，要得到滞后ik采样值(角度的值ik(-()。 （7） （8） 由式（7）得： （9） 由式（9）得： （10） 将（10）带入（8）中得： （11） 其中：。 实际应用中，可以通过提高采样频率来降低移相产生的误差。 通过A /D采集，得到检测信号序列和参考信号序列，通过软件移相算法对参考序列进行移相，得到与正交的另外一路参考信号序列；，和可以通过采样序列经过下面的公式得到（N为采样点数）: 　　　　　　　　　　　　（12） 　　　　　　　　　　　　（13） 　　　　　　　　　　　　（14） 根据得到的，和求解式（12）、式（13）和式（14）可以得到信号的幅值A和与参考信号的相位差h。 2.3 软件实现流程框图 软件实现流图如图2所示： 图 2 软件功能实现流图 3　利用软正交矢量LIA方法检测微弱信号的Matlab仿真结果 被测信号频率为1 kHz，幅值A= 0. 5V，N（t）为噪声信号；参考信号的频率也为1 kHz，幅值B=1V，相位差h= 45°；设置采样频率为50 kHz，同时对参考与被测信号采样，将参考波形序列，通过软件移相，得到另一路参考信号。 编写Matlab程序，测试时输入采样频率、方差、采样点数等参数，要求画出被测信号波形、经软件移相产生的两路参考波形、采样被测信号、采样参考信号的波形图，并给出正交计算得到的A、B、h的值。 首先以采样频率为50kHz，采样点数为500，方差为0.5V为例，运行程序得到被测信号波形和参考信号波形分别如图3、图4所示： 图 3 方差为0.5时的被测信号波形 图 4 标准参考信号和移相参考信号 对方差为0.5V不同采样点数和采样频率的所有情况进行测试，测试结果如表1所示： 表 1 方差为0.5V时的幅值相位 噪声方差 采样频率/kHz 50 100 采样点数 500 1000 1000 2000 幅值A/V 0.5284 0.5146 0.5072 0.5016 幅值误差/% 5.68 2.91 1.43 0.31 相位h/° 44