发明名称 一种聚焦增强的超宽带成像方法及装置
摘要 本发明公开了一种聚焦增强的超宽带成像方法及装置,控制器FPGA1控制时域信号发射机5发出窄脉冲信号,此脉冲信号宽度纳秒级,具有较高带宽的频谱。窄脉冲信号通过发射天线阵列12辐射,对目标物进行探测。接收天线阵列13将接收到的微弱反射信号送给低噪声放大器11,信号经放大后进入接收机7,接收机7将接收到的高速微弱信号解析成数字信号,通过多路数据采集单元10送入信号处理与成像模块4。同时,控制器FPGA1控制惯性定位传感器2实时测量装置的加速度,数据采集单元3将测量结果送给信号处理与成像模块4,最后由成像算法程序解析出目标物的位置。
申请公布号 CN103176182B 申请公布日期 2015.11.18
申请号 CN201310077279.9 申请日期 2013.03.11
申请人 山东省科学院自动化研究所 发明人 王丰贵;杨秀蔚;王忠民;崔洪亮;成巍;杨传法
分类号 G01S13/89(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 济南圣达知识产权代理有限公司 37221 代理人 张勇
主权项 一种聚焦增强的超宽带成像装置的成像方法,所述聚焦增强的超宽带成像装置,它包括:发射天线阵列,用于将时域信号发射机发射的时域信号辐射,对目标物进行探测;接收天线阵列,用于接收目标物反射的反射信号并将接收到的反射信号送入接收机,接收机将接收到的反射信号解析成数字信号,通过数据采集单元I送入信号处理与成像模块成像;控制器FPGA,控制所述时域信号发射机发射时域信号,同时还控制接收机接收反射信号及控制定位装置实时测量天线装置的加速度,所述定位装置通过数据采集单元II与信号处理与成像模块连接,所述天线装置由发射天线阵列和接收天线阵列组成;其特征是,具体步骤为:步骤一:控制器FPGA控制时域信号发射机发出时域信号,时域信号通过发射天线阵列辐射,对目标物进行探测;步骤二:接收天线阵列将接收到的目标物反射的反射信号经低噪声放大器送入接收机,控制器FPGA控制接收机采用等效采样的方式接收反射信号并将接收到的反射信号解析为数字信号;步骤三:接收机将解析成的数字信号,通过数据采集单元I送入信号处理与成像模块,同时,FPGA控制惯性定位传感器实时测量天线装置的加速度,数据采集单元II将测量结果送给信号处理与成像模块;步骤四:信号处理与成像模块根据接收到的反射信号和由惯性定位传感器测量的加速度,利用成像算法解析出目标物的位置;所述步骤四的具体步骤为:(4‑1)装置具有一路发射天线和多路接收天线,装置紧贴墙体移动,对墙体另一侧的探测区域任意点目标A进行探测,发射信号被目标A反射后再次穿过墙体,由接收天线进行接收,设目标A的坐标是(x,y),x是探测区域横向坐标,y是探测区域纵向坐标,x和y的数值与垂直方向坐标z无关;根据空气和墙的介电常数、电导率、磁导率确定出电磁波在墙体中的传播速度以及电磁波由空气入射到墙和由墙入射到空气时,入射角与折射角的关系;将空气视为理想介质,其介电常数ε<sub>1</sub>=ε<sub>0</sub>,磁导率μ<sub>1</sub>=μ<sub>0</sub>,电导率δ<sub>1</sub>=0,电磁波在理想介质中的传播速度为光速c;在墙壁内,假设墙壁为均匀介质,介电常数为ε<sub>2</sub>=ε<sub>r</sub>ε<sub>0</sub>,磁导率μ<sub>2</sub>=μ<sub>r</sub>μ<sub>0</sub>,电导率为δ<sub>2</sub>,假设电磁波的角频率为ω,则电磁波在墙壁内部的传播速度为:<img file="FDA0000783145180000021.GIF" wi="1285" he="255" />其中,ε<sub>0</sub>真空介电常数,ε<sub>r</sub>相对介电常数,μ<sub>0</sub>真空磁导率,μ<sub>r</sub>相对磁导率;设T是发射天线,坐标为(x<sub>0</sub>,y<sub>0</sub>);设R<sub>1</sub>,R<sub>2</sub>,R<sub>3</sub>,R<sub>4</sub>分别是四路接收天线,对应的坐标分别为(x<sub>1</sub>,y<sub>1</sub>),(x<sub>2</sub>,y<sub>2</sub>),(x<sub>3</sub>,y<sub>3</sub>),(x<sub>4</sub>,y<sub>4</sub>),设R<sub>i</sub>代表接收天线中的第i路,对应的坐标为(x<sub>i</sub>,y<sub>i</sub>),其中i=1,2,3,4;M<sub>0</sub>代表发射信号电磁波穿透墙壁的折射点,M<sub>1</sub>、M<sub>2</sub>、M<sub>3</sub>和M<sub>4</sub>代表四路反射信号到墙壁上的折射点,M<sub>i</sub>代表第i路反射信号到墙壁上的折射点,其中i=1,2,3,4;由于墙壁表面的折射,电磁波从发射天线到A点的路径被分为L<sub>01</sub>和L<sub>02</sub>两段,设L<sub>01</sub>是发射信号由折射点M<sub>0</sub>到目标的距离,L<sub>02</sub>是发射天线发射的信号在墙壁内的距离;由于墙壁表面的折射,第i路反射信号从目标到接收天线R<sub>i</sub>的路径被分为L<sub>i1</sub>,L<sub>i2</sub>两段,其中i=1,2,3,4;设L<sub>i1</sub>是第i路反射信号从目标到折射点M<sub>i</sub>的距离,L<sub>i2</sub>是第i路反射信号从折射点M<sub>i</sub>到接收天线R<sub>i</sub>的距离;由Snell折射定律可得电磁波的入射角α与折射角β的关系式为:<img file="FDA0000783145180000031.GIF" wi="1276" he="136" />式中,<img file="FDA0000783145180000032.GIF" wi="218" he="88" />是复数形式,j是虚数单位;<img file="FDA0000783145180000033.GIF" wi="347" he="80" />,k<sub>1</sub>以rad/m计      (3)<img file="FDA0000783145180000034.GIF" wi="1369" he="225" /><img file="FDA0000783145180000035.GIF" wi="1373" he="212" />将电磁波从发射天线到目标A,然后从目标A反射到任一个接收天线的传播路径分为四段,发射天线到目标A的两段传播路径分别为L<sub>01</sub>,L<sub>02</sub>,目标A到第i路接收天线R<sub>i</sub>的两段传播路径为L<sub>i1</sub>,L<sub>i2</sub>,其中i=1,2,3,4;发射天线发出的任一角频率ω的电磁波经过目标A后反射到第i路天线R<sub>i</sub>的传输时间为:<img file="FDA0000783145180000036.GIF" wi="1373" he="128" />其中,v<sub>ω</sub>是电磁波在墙体中的传播速度;令墙体的磁导率μ<sub>2</sub>=μ<sub>0</sub>,墙体厚度远小于探测区域尺寸,设电导率δ<sub>2</sub>=0,则式(4)和(5)简化为:<img file="FDA0000783145180000041.GIF" wi="1201" he="90" />k″<sub>2</sub>=0       (8) 则电磁波的色散现象消失,墙壁内的电磁波传播速度为:<img file="FDA0000783145180000042.GIF" wi="1243" he="118" />则电磁波由空气入射到墙壁时,入射角与折射角的关系式为<img file="FDA0000783145180000043.GIF" wi="1268" he="123" />设α<sub>0</sub>是电磁波发射时的折射角;β<sub>0</sub>是电磁波发射时的入射角;α<sub>i</sub>是第i路反射的电磁波在达到墙壁时的入射角;β<sub>i</sub>是第i路反射的电磁波在达到墙壁时的折射角;d是折射点与第i路接收天线R<sub>i</sub>在横向x轴上的距离,d=|x<sub>i</sub>‑x<sub>Mi</sub>|;发射信号电磁波穿透墙壁的折射点M<sub>0</sub>的坐标是(x<sub>M0</sub>,y<sub>M0</sub>);第i路反射信号到墙壁上的折射点M<sub>i</sub>的坐标是(x<sub>Mi</sub>,y<sub>Mi</sub>);目标A到第i路接收天线R<sub>i</sub>的路径中,设Δx=|x‑x<sub>i</sub>|,Δy=y‑y<sub>i</sub>‑d<sub>w</sub>,d=|x<sub>i</sub>‑x<sub>Mi</sub>|,其中,i=1,2,3,4;若Δx≠0,则<img file="FDA0000783145180000044.GIF" wi="1339" he="128" /><img file="FDA0000783145180000045.GIF" wi="1333" he="128" />(4‑2)根据电磁波在空气和墙体中的传播速度,以及入射角与折射角的关系,确定关于折射点的方程;由(10)、(11)、(12)式得到关于d的一元四次方程:<img file="FDA0000783145180000046.GIF" wi="1696" he="138" />其中,d<sub>w</sub>代表墙壁厚度;(4‑3)根据公式(13),确定d的值,再利用d值得到第i路反射信号到墙壁上的折射点M<sub>i</sub>的坐标;(4‑4)根据折射点M<sub>i</sub>坐标,得出延迟参数t<sub>i,</sub><sub>ω</sub>;根据折射点坐标,得:<img file="FDA0000783145180000051.GIF" wi="1949" he="470" />(4‑5)利用惯性定位传感器采集的天线装置的加速度数据,得出此时刻天线装置的位移,并确定探测区域的新坐标;装置在移动过程中的加速度a是空间中任意方向上的矢量,与横向x轴的夹角是θ<sub>x</sub>,与纵向y轴的夹角是θ<sub>y</sub>,与垂直方向z轴的夹角是θ<sub>z</sub>,在三维坐标系中,可以分解为横向加速度a<sub>x</sub>,纵向加速度a<sub>y</sub>和垂直加速度a<sub>z</sub>;装置发生移动,则固定在装置上的整个天线阵列发生位移;假如天线阵列移动的加速度为矢量a,是空间中任意方向上的矢量,且是时间的函数,设时间为Δt<sub>s</sub>,则速度矢量为:<img file="FDA0000783145180000052.GIF" wi="1031" he="114" />假如天线阵列的起始位置为P<sub>0</sub>,移动到P<sub>1</sub>位置时所经过的距离为:<img file="FDA0000783145180000053.GIF" wi="1110" he="111" />由于整个天线阵列发生位移时是紧贴墙面运动的,则忽略纵向加速度a<sub>y</sub>,由于探测区域任意点目标A的坐标在横向x轴和纵向y轴方向上的数值与垂直方向z轴无关,则忽略垂直加速度a<sub>z</sub>,则公式(15)和(16)简化为:<img file="FDA0000783145180000061.GIF" wi="911" he="116" /><img file="FDA0000783145180000062.GIF" wi="899" he="111" />将发射天线位置作为当前坐标原点,确定第j次成像时探测区域的新坐标<img file="FDA0000783145180000063.GIF" wi="1050" he="102" /><img file="FDA0000783145180000064.GIF" wi="97" he="101" />表示第j次成像前每个天线相对于上次成像时的位移;(4‑6)根据延迟参数t<sub>i,</sub><sub>ω</sub>和探测区域的新坐标,完成探测区域的单次成像;天线阵列紧贴墙面移动,探测区域的单次成像结果为:<img file="FDA0000783145180000065.GIF" wi="1849" he="286" />s<sub>i</sub>(t<sub>i</sub>)表示第i路接收天线接收到的目标反射信号在t<sub>i</sub>时刻的信号幅度值,其中,i=1,2,3,4;(4‑7)对单次成像的结果进行合成:以最后一次成像时发射天线的位置作为当前坐标原点,对多个单次成像的结果进行坐标变换,然后进行合成;<img file="FDA0000783145180000066.GIF" wi="1603" he="239" />当总成像次数小于N时,对第1~j次单次成像的结果进行合成;当总成像次数大于等于N时,对第j‑N+1~j次单次成像的结果进行合成;若装置要求完成一次合成成像的时间不超过T<sub>syn</sub>,假定单次成像的周期为ΔT<sub>s</sub>,则需合成单次成像的个数为<img file="FDA0000783145180000067.GIF" wi="242" he="125" />
地址 250014 山东省济南市历下区科院路19号