基于FPGA的生命探測儀算法研究與系統(tǒng)設(shè)計

摘要：給出了一種基于FPGA的生命探測信號處理系統(tǒng)的設(shè)計方法。從理論上研究了生命探測儀的算法及其軟硬件系統(tǒng)。其中在FPGA軟件設(shè)計中利用模塊化的思想方法分別設(shè)計了FIR濾波器、異步FIFO、UART、電池監(jiān)控、功能控制等功能模塊。最后完成人體特征信號和體動信號的分析與提取，實現(xiàn)了非接觸情況下生命探測與發(fā)現(xiàn)。相對于傳統(tǒng)的生命探測儀，該設(shè)備具有體積小，功耗低，操作簡單，攜帶方便等優(yōu)點，特別適用于野外和戰(zhàn)場生命探測等應(yīng)用場合。
關(guān)鍵詞：生命探測儀；FPGA；FIR濾波器；FIFO；UART

0 引言
雷達式非接觸生命探測技術(shù)是近年來發(fā)展起來的一種新技術(shù)，是融合雷達技術(shù)、生物醫(yī)學(xué)工程技術(shù)于一體，可穿透非金屬介質(zhì)(磚墻、廢墟等)，不需要任何電極或傳感器接觸生命體，可在較遠的距離內(nèi)探測到生命體的生命信號(呼吸、體動)的一種特殊電子裝置。該探測儀克服了基于激光、紅外生命探測儀受溫度影響嚴重、遇物體阻擋失效的問題，也克服了超聲探測空間傳播衰減大、受環(huán)境雜物反射干擾、水、冰、泥土阻擋失效的問題，因此近年來備受國內(nèi)外學(xué)者的關(guān)注。其基本原理是利用雷達天線發(fā)射的電磁波穿透一定障礙物照射到人體時，反射的回波信號受到人體生理運動(如心跳、呼吸)引起表面微動的多普勒調(diào)制，人體表面微動信號就加載到了反射波中，這種人體的微動與回波幅度相位之間具有相關(guān)性。對回波信號進行A／D轉(zhuǎn)換、濾波等處理就可以提取到人體的生命特征信息。
現(xiàn)有設(shè)備多是基于單片機與PC機顯控系統(tǒng)，具有系統(tǒng)體積較大，靈活性差，實時性差等缺點。本文主要研究生命探測儀算法和基于FPGA的信號處理系統(tǒng)設(shè)計，上位機采用基于嵌入式系統(tǒng)的顯示控制單元。該系統(tǒng)體積小、成本低、攜帶方便，具有很高的實用價值，可廣泛應(yīng)用于災(zāi)害救援(地震、塌方傷員的探尋)、反恐斗爭(隔墻監(jiān)控罪犯、解救人質(zhì))等場合。

1 生命探測儀算法研究
1．1 生命探測雷達工作原理
生命探測雷達的基本原理是多普勒效應(yīng)。當(dāng)發(fā)射源與接收者之間有相對徑向運動時，接收到的信號頻率將發(fā)生變化，其頻率差別與兩者的相對運動速度向量有關(guān)，這種現(xiàn)象被稱為多普勒頻移。多普勒雷達發(fā)送連續(xù)的電磁波信號到被探測對象，返回的信號被調(diào)制而具有被探測對象運動的信息，因此只要解調(diào)出返回信號就可以獲得所要結(jié)果。
如果忽略幅度的變化，由單頻連續(xù)波雷達發(fā)射的信號可以表示為：
s(t)=Acos(2πfct+φ) (1)
式(1)中，fc是雷達的發(fā)射波頻率，φ為初相，A為振幅。
雷達接收機接收到的目標(biāo)回波信號sr(t)為：
sr(t)=Ks(t-tr)=KAcos[2πfc(t-tr)+φ] (2)
式(2)中，tr=2R／c為回波信號滯后于發(fā)射信號的時間：R為目標(biāo)和雷達之間的距離：c為電磁波傳播速度，在自由空間傳播時，c等于光速；K為回波的衰減系數(shù)。
當(dāng)目標(biāo)和雷達之間有相對運動時，距離R隨時間變化。設(shè)目標(biāo)以勻速相對雷達運動，則在t時刻目標(biāo)與雷達的距離R月(t)為：
R(t)=R0-vrt (3)
式(3)中，R0為t=0時的距離，vr為目標(biāo)相對雷達的徑向運動速度。
由于通常雷達和目標(biāo)間的相對速度vr遠小于電磁波速度c，故時延tr可近似寫為：

多普勒頻率正比于相對運動的速度，而反比于工作波長λ。當(dāng)目標(biāo)以接近雷達的方向運動時，多普勒頻率為正值，接收信號頻率高于發(fā)射信號頻率；當(dāng)目標(biāo)以背離雷達的方向運動時，多普勒頻率為負值，接收信號的頻率低于發(fā)射信號的頻率。
生命的特征在于運動，如人體呼吸和心跳時胸腔的運動以及人的體動等。當(dāng)電磁波照射到人體時，反射波會產(chǎn)生多普勒頻移。根據(jù)人的呼吸、心跳和體動產(chǎn)生的多普勒頻率的特性，可以探測生命體是否存在。
1．2 生命探測雷達回波信號分析
人體生命體征信號(呼吸、心跳和體動)是一種極微弱的低速目標(biāo)信號。理論上，人體連續(xù)的心臟跳動和呼吸會引起胸腔的起伏運動，進而會產(chǎn)生系列多普勒頻移信號，這些信號可以用連續(xù)波的形式來表示。實際中，由于雷達內(nèi)部運動、人體體動等情況，目標(biāo)回波信號除了包括心跳和呼吸信號外還包括其它連續(xù)波分量。將人體簡化為復(fù)合介電常數(shù)的球體和圓柱體模型。設(shè)人體生命體征運動(呼吸和心跳)是頻率為Ω，幅度為A的簡諧振動：
X(t)=Acos(Ωt+φ) (7)
設(shè)質(zhì)點的振動速度為：
v(t)=dX／dt=-AΩsin(Ωt+φ) (8)
則雷達發(fā)射信號入射到人體表面產(chǎn)生的多普勒頻移為：
fd(t)=2v(t)／λ=-(2AΩ／λ)sin(Ωt+φ) (9)
式(9)中，λ為雷達發(fā)射信號的波長。若設(shè)發(fā)射信號為St(t)=Atcos(ωct+θ0)，則經(jīng)過人體表面反射后的雷達信號為：
Sr(t)=Arcos[ωc(t-tr)]+θ0+kg(t) (10)
式(10)中，為回波信號滯后于發(fā)射信號的時間；k為常數(shù)，用于反映人體振動面對雷達發(fā)射信號的相位調(diào)制度，且k|g(t)|1。分析可知：
Sr(t)≈Arcos(ωct+θ0)-Arkg(t)sin(ωct+θ0) (11)
Sr(t)經(jīng)過放大后與相干本振信號sln(ωct+θ0)進行相干混頻，再經(jīng)過低通濾波，即可得出人體生命信號：
a(t)=-(ArGk／λ)X(t) (12)
由式(12)可見，雷達回波信號經(jīng)過處理后所得的人體生命體征信號a(t)的輸出幅度與雷達的發(fā)射信號Ar、人體振動面雷達相位調(diào)制度k和雷達系統(tǒng)處理增益G成正比，與雷達發(fā)射信號波長λ成反比。
人體的呼吸、心跳信號是一種窄帶、低幅值、準(zhǔn)周期信號，易受噪聲和環(huán)境干擾的影響。呼吸頻率的不均勻和胸腔的多點反射導(dǎo)致信號的頻譜有一定展寬。因為回波信號體現(xiàn)了人體呼吸和心跳產(chǎn)生的多普勒頻移，在頻域上出現(xiàn)不同的諧波，可以將心跳和呼吸的回波信號表示成多個多普勒頻移信號的組合：

式(3)中，ai、fi和φi分別為第i個諧波分量的幅度、歸一化頻率和初始相位，p是諧波個數(shù)。