system: [SYSTEM]

user: 2. 系统方案设计与实现 2.1. 系统原理 由于生命信号的频率很小，如果要分辨出生命信号，就要采用高频电磁波。目前对穿墙生命信号探测最有利的电磁波段为毫米波和厘米波，一般选择在L或者S波段。Michigan州立大学和希腊雅典国家技术大学的科学家通过不断地研制、试验，最后确定探测生命信号的最佳频率为S波段2.45 GHz [ 9 ] [ 10 ] 。 在生命体情绪平稳的时候，心跳和呼吸的频率维持在一个稳定的范围，呈周期性的变化，可以从时域和频域的积累上观察到这一特性。由多普勒效应公式可知，当雷达采用不同的发射频率时，所得生命信号的频率变化范围不同。俄罗斯莫斯科科学技术大学的科学家们利用他们的单频探地雷达进行了生命信号探测的实验，雷达发射频率为1~10 GHz的电磁波，得到心跳和呼吸的频率变化范围是：心跳0.8~2.5 Hz，呼吸0.2~0.5 Hz。当发射电磁波频率为1.6 GHz时，得到静止人的稳定心跳频率为1.6 Hz，呼吸频率为0.03 Hz。对于最佳探测频率2.45 GHz，人的生命信号的频率变化范围是：0.1~3 Hz。 雷达(Radar， Radio Detection and ranging)发送电磁波信号，并感应目标引起的回波，从而获得这些目标的信息。一个典型的雷达系统包括发射机、天线、接收机、信号处理软硬件。发射机产生射频波形，并放大到足够的功率。雷达天线收发电磁波，将能量集中于目标的方向，并测定目标。接收机将射频回波信号转换为中频或基带信号，并分离信号和噪声、干扰，最后将信号放大到足够的幅度以满足后续数字电路的需求。 在非接触生命探测过程中，雷达发射的电磁波照射人体，根据多普勒原理，人体反射的回波信号被人体生命活动(如呼吸运动等)引起的体表微动所调制，使得这些回波信号的某些参数(如频率、相位)发生改变 [ 11 ] ，经过后期信号处理，提取出人体的生命参数。人体主要是依靠心跳和呼吸时的胸腔运动和体动对雷达发出的电磁波产生多普勒效应。 该系统就是模拟人体的胸壁起伏，系统的设计思路也是由此而来。由心肺活动频率相当的电信号产生信号源，通过声卡由扬声器末端鼓膜的往复运动产生，本系统就是运用以上原理来模拟人体的心肺运动，得到定量稳定的目标信号，从而进行非接触生命探测实验。 2.2. 雷达监测模拟系统设计 根据非接触式生物雷达监测系统的工作原理，进行整个系统设计。本项目采用2.5 GHz的连续波雷达来发射信号，通过天线(收发双置天线)接收回波；经低噪放大、射频滤波、混频等一系列处理，将射频信号转换到中频信号送入数字下变频模块；其后对信号进行带通滤波，将信号下变频到零中频信号，并通过抽取滤波处理来实现采样率的转换，数字下变频模块输出低速的I/Q两路基带信号送至数字信号处理模块；数字信号处理先对I/Q两路基带信号进行反正切解调，提取体征信号的相位信息，再通过自适应滤波算法来进行呼吸和心跳信号的分离；最终在控显终端获得被测人体的具体体征信息，来对人体的健康状况进行推断。 生物雷达监测系统主要包括雷达探测前端、信号预处理模块、数字下变频模块、数字信号处理模块和控显终端五个部分。 雷达探测前端主要包括产生发射连续波信号的本地振荡器和将发射信号发送给目标物体的发射天线。振荡器产生的信号经定向偶合器，一路通过环行器由天线发射出去，另一路去混频器。当发射信号遇到人体等目标后反射回来被接受天线接受，并送入混频器，该信号再经过放大、预处理、A/D转换，通过后面的自适应滤波器在频域上将呼吸、心跳信号分开，送到PC机进行处理、显示。 预处理模块由课题组前期自行研制开发，主要组成部分分为前置输入级、主放大级、滤波器组、增益控制级、50 Hz陷波级、信号输出级。前置级主要是滤波和保护电路；主放大级设置60 dB的增益放大模块实现对弱信号的放大，对强信号的限幅；滤波器组分为0.5 Hz高通滤波器组和3.3 Hz的低通滤波器组，经行对信号的选频，同时滤除直流和高频的干扰；增益控制电路设置调节模块实现对输出信号增益倍数的调节；陷波电路滤除50 Hz的工频干扰。预处理的目的就是提高混频后的低频信号的信噪比，以利于后面A/D转换和数字信号处理。 数字下变频部分包含A/D采样和数字下变频两个过程。A/D采样是在中频进行带通采样，中频频率为75 MHz，采样频率为100 MHz；数字下变频从高速中频数字信号中提取有用信号，将75 MHz中频信号下变频至零中频信号，下变频采用数字正交解调。先对中频信号 进行采样，采样的数字信号进入数字正交解调环节，通过改变采样频率来进行选取我们所需要解调的信号，通过奇偶抽取与符号校正，把中频信号中的载频(不含有用信息)下变频成零载频，中频信号变成I/Q两路正交的数字基带(零中频)信号，并采用时延滤波器将这两路信号采样点在时间上对齐。然后再将这两路信号送入DSP中进行后续处理。数字下变频的主要作用：其一是把中频信号下变频至零中频信号；其二是降低采样率。 数字信号处理部分主要包括反正切解调和自适应滤波算法。从数字下变频模块输出的I/Q两路基带信号送入数字信号处理部分，经过反正切解调提取I/Q两路基带信号的信息，根据心跳和呼吸信号的频域不同，利用自适应滤波技术将这两信号进行分离，最后用傅立叶变换进行这两种信号的频率提取，将结果送至PC显示，最终获得被测者实际体征情况，数字信号处理部分由DSP来实现。 控显终端采用的是一款基于Windows CE操作系统下专为嵌入式控制而设计的控制平台LJD-eWin4S。该产品具有实时性好、移植性强、处理速度快、多任务切换等优势。系统框图如图1所示。 2.3. 基于雷达监测系统的生命模拟系统组成 模拟系统包括信号源、功率放大器和电–振动换能器三部分组成。信号源产生与心肺活动频率相当的电信号，在该实验系统中，实验模拟信号由计算机中的音频播放软件通过声卡发出。音频源文件的产生是基于单次的心跳和呼吸信号模板，它来自典型的人体呼吸和心跳信号振动信号，由医学数据库中的真实人体心肺相关数据记录转换而来，或通过实验测定来获取。然后将心跳与呼吸信号的模板按二者幅度的比例合成，并根据需要仿真的心率和呼吸频率，进行插值和重采样，形成一段心肺运动波形数据，存储为音频文件。音频文件由Adobe Audition等软件播放，即可从声卡输出心肺振动的波形。 由于心肺模拟信号频率相对音频要低很多，尤其是呼吸的低频段可能低于0.2 Hz，普通音频输出装置并不支持20 Hz以下的低频信号输出，因而在声卡的设计上需要做到支持低频和直流信号的传输。具体实现方法是用TI的USB音频芯片PCM2902实现了数模(DA)转换，然后用运算放大器以直流稱合方式进行电平校正和放大。功率放大器将信号放大至足以推动换能器产生与人体振动幅度相仿的功率。同样出于频率范围的考虑，功率放大器也需要支持频率低至直流的信号。电–振动换能器用低音扬声器实现。考虑到人呼吸时胸壁运动幅度较大，选用线性行程较长的低音扬声器，例如HiViResearch的SP10 (线性位移15.5 mm)、seas的D1001 (线性行程28 mm)等等。另外，在扬声器表面覆盖一层金属平面，以获得较强的电磁波反射能力，更有效的模拟体表反射。生物雷达监测系统实验设计的实物图如图2所示。 这样实验系统使音响末端鼓膜产生3.00 cm位移、频率为2 Hz的周期正弦信号。生物雷达监测系统的发射接收模块平行于音响末端并保持一定距离。这样雷达发射接受模块以50 Hz的频率采集反射回来的雷达信号。每个距离采集两组数据：有模拟源数据和无模拟源数据。其中受检距离为1.6 cm~600.0 cm。 图1. 生物雷达监测系统框图 图2. 生物雷达监测系统实物图 并完成了这些所采集数据的处理，包括数据的筛选、解调、加窗、FFT变换等等。生物雷达监测系统实验设计的示意图如图3所示。 2.4. 系统组成 本章对生命信号检测试验系统、对试验数据处理(包括去尖峰、解卷绕、滤除基线、去噪、加窗和FFT)过程进行了设计，最后根据仿真结果给出结论。其中，采用滑动滤波去除的基线漂移，基于FIR数字滤波器的线性去噪，汉宁窗的选择，尤其是可变点数的FFT变换等，对实验仿真结果给出了有力证据。 实验数据采集及预处理过程如下：首先，获取数据，一定距离采集两组数据包括：有模拟源的截取的6000个数据，没有模拟源的截取的3000个数据。第二，对所有数据借助MATLAB去尖峰处理，以便处理掉采集过程中可能出现的异点和误码。第三，对所有数据加窗处理，由于需要做频域分析，故采用汉宁窗，窗长为数据的长度，为使检测结果明显故需要增大比例，这里的系数扩大100倍。第四，对所有数据FFT处理，FFT处理后的结果做平方处理再次增大结果比例。第五，借助MATLAB中的SUM函数对所有数据在频域内求它们的功率，得到4组功率数据：去声噪I路信号功率、去声噪Q路信号功率、I路声噪功率、Q路声噪功率。第六，求去声噪的有用总功率、信噪比及基带信号的信噪比。 本实验数据获取部分主要是将各个模块进行协调和进行数据交互。首先完成初始化工作，初始化后进入循环处理，在循环过程中主要处理获得采集的数据，并将得到的数据发送到上位机。采集到数据后，首先删除明显的异常数据，由计算机软件通过设置阈值来自动完成。针对心跳的节律性，可近似认为其为准周期信号，对其进行相关处理，其自相关函数也具有周期性，而噪声和短时随机干扰在一定程度上被抑制。并且采样数据以文本格式(.txt)存储在指定位置。 实验采样频率为50 Hz，故每一个距离采集2 min有源数据6000点以及采集1 min有源数据3000点，此时数据未经任何处理。实时数据采集如图4所示。 获得数据后，由于误码或者突然人为或者环境改变造成数据的急剧变化，这些突出的数据时无效数据需要去尖峰处理。图5是I路去尖峰前后数据对比。 图3. 生物雷达监测系统示意图 图4. 实时数据 图5. I路去尖峰前后数据对比 对于生命探测，我们希望在距离上将多个生命体目标分辨开来。为提高距离分辨力，选择汉宁窗对离散差拍信号进行加权处理。由于频域加权处理对应时域是周期加权处理，因此在信号处理中，对每个周期的离散差拍信号，都要进行汉宁窗加权处理。线性调频连续波信号的第一级FFT输出的频谱与距离对应，用频域加权匹配处理可有效抑制距离旁瓣，提高多目标的分辨能力，但无法避免主瓣变宽，并使线性调频连续波信号的匹配滤波输出有信噪比损失。若采用汉宁加权，信噪比损失为3.52 dB [ 12 ] 。图6是生命雷达信号的加窗处理。 本文所提及的信噪比(SNR)其计算公式为： 其中， 及 分别表示信号及噪声在频域里的幅值。 基于FFT的频域积累借助快速傅里叶变换将信号由时域转换到频域后进行积累，这里采用可变点数FFT积累即不断增加参与FFT运算的时间点数，达到频域积累目的。图7、图8为该仿真信号常规的1024点FFT频域和时域仿真图，为使结果明显数据的频率都扩大了10。图中横坐标表示信号不同的频率成分，纵坐标表示各频率分量在频域内的幅度，此时可以清楚的看到心率在0.3 Hz左右，呼吸在1 Hz左右。