CN109283529B - 一种基于完全随机波束照射的可定制微波成像系统及算法 - Google Patents
一种基于完全随机波束照射的可定制微波成像系统及算法 Download PDFInfo
- Publication number
- CN109283529B CN109283529B CN201811073445.7A CN201811073445A CN109283529B CN 109283529 B CN109283529 B CN 109283529B CN 201811073445 A CN201811073445 A CN 201811073445A CN 109283529 B CN109283529 B CN 109283529B
- Authority
- CN
- China
- Prior art keywords
- imaging
- random
- transmitting
- transmitting antenna
- transmitting antennas
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 107
- 238000004422 calculation algorithm Methods 0.000 title abstract description 8
- 238000003491 array Methods 0.000 claims abstract description 7
- 230000005284 excitation Effects 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 26
- 230000005540 biological transmission Effects 0.000 claims description 14
- 238000005286 illumination Methods 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 abstract description 5
- 238000005259 measurement Methods 0.000 description 13
- 239000013598 vector Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004613 tight binding model Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于完全随机波束照射的可定制微波成像系统及算法。每个发射天线经移相器连接到发射端口,接收射天线连接到接收端口,多个发射天线和接收天线布置于待测目标的一侧,多个发射天线阵列排布形成发射天线阵列;通过各个移相器随机调节各个发射天线的激励电流使得所有发射天线的发出随机电磁波束照射到待测目标所在的成像区域;成像后通过不断调整成像分辨率、成像区域的大小、成像距离、发射天线阵列中的发射天线总数、发射天线之间的间距和随机波束照射的次数,获得最优的奇异值分布曲线,最优的奇异值分布曲线上的所有奇异值之和最大。本发明成像系统具有结构紧凑,实现简单、随机性能优良,能量效率高,频谱效率高等优点。
Description
技术领域
本发明涉及了一种微波成像系统及算法,尤其是涉及了一种基于完全随机波束照射的可定制微波成像系统及算法。
背景技术
近些年来,微波成像的理论和技术得到了广泛的研究和发展。微波是频率在300MHz-300GHz,对应波长为1m-1mm的电磁波。与可见光波不同,在微波频段具有对非金属物质的穿透能力。利用这一独特性质,微波成像能够对获取目标的位置、形状、材料属性等多种信息量。微波成像具有安全、非接触、成本低等特点,因此广泛应用于安全检查、隔墙监视、医学成像等军用、民用领域。
不同于传统的微波成像方法,基于随机照射的微波成像方法引起了国内外研究人员的关注。参考光学中的ghost成像理论,将随时间变化的随机辐射场照射在物体上,以随机照射的方式获取多组目标散射测量值,通过入射场与散射场的反演计算就能够提取散射体的轮廓和表面信息。不同于传统电磁成像中采用结构复杂的天线阵进行波束扫描的相控阵雷达成像或者移动天线进行空间扫描的合成孔径雷达成像,随机照射方式的微波成像方法在有限空间和无法进行相位扫描的场合有着突出应用价值。从天线的角度来看,这种随机照射可以通过频率扫描天线或者相控阵天线来实现。其中频率扫描天线往往依赖于大带宽,其频谱利用效率低下,极大地浪费了宝贵的频谱资源。而工作于单一频率的成像系统具有频谱效率高的优点。
发明内容
为了解决背景技术中存在的问题,本发明是提供了一种基于完全随机波束照射的可定制微波成像系统及算法。
本发明解决其技术问题所采用的技术方案是:
一、一种基于随机波束照射的可定制微波成像系统:
微波成像系统包括多个发射天线和一个接收天线,每个发射天线经移相器连接到发射端口,接收射天线连接到接收端口,多个发射天线和接收天线布置于待测目标的一侧,多个发射天线阵列排布形成发射天线阵列;通过各个移相器随机调节各个发射天线的激励电流相位使得所有发射天线的发出随机电磁波束照射到待测目标所在的成像区域。
所述的多个发射天线阵列排列方式呈方形或者圆形。
所有发射天线朝向待测目标发出电磁波作为入射场,入射场与待测目标相互作用形成散射场,接收天线对散射场进行接收,收发天线工作于单一频率。
二、一种基于随机波束照射的可定制微波成像算法,方法采用以下步骤:
1)搭建上述微波成像系统进行随机波束照射,采用1比特随机相位调制;
2)经过多次随机波束照射,由接收天线采集的接收信号构建传输矩阵;
3)预设成像分辨率、成像区域的大小、成像距离、发射天线阵列中的发射天线总数、发射天线之间的间距,多次随机波束照射后计算传输矩阵的奇异值分布曲线,通过不断调整成像分辨率、成像区域的大小、成像距离、发射天线阵列中的发射天线总数、发射天线之间的间距和随机波束照射的次数,获得最优的奇异值分布曲线,最优的奇异值分布曲线上的所有奇异值之和最大;
成像分辨率是指成像区域的网格大小。成像距离是指天线和成像区域之间的距离。
4)取最优的奇异值分布曲线下对应的成像分辨率、成像区域的大小、成像距离、发射天线阵列中的发射天线总数、发射天线之间的间距、随机波束照射的次数设置微波成像系统进行成像。
采用1比特随机相位调制的发射天线阵列实现完全随机的波束照射。
所述1比特随机相位调制具体是从0度至360度的相位值中选取任意两个不同的相位值进行调制。
所述1比特随机相位调制通过采用射频开关切换不同长度的传输线实现。
所述的随机波束照射的最少次数与发射天线阵列中的发射天线总数量相同。
本发明具有结构紧凑,随机性能优良,能量效率高,频谱效率高等优点。
附图说明
图1是本发明基于随机波束照射的微波成像示意图。
图2是本发明基于随机波束照射的微波成像系统框图。
图3是本发明实施例归一化奇异值曲线图。
图4是本发明实施例随测量次数和成像距离变化的奇异值分布曲线图。
图5是本发明实施例的图4在测量次数为25处的二维剖面图。
图6是本发明实施例不同天线组阵参数和成像距离和分辨率的关系曲线图。
图7是本发明实施例不同天线单元个数N和成像距分辨率之间的关系曲线图。
图8是本发明实施例的不同成像距离的仿真成像结果图。
图9是本发明实施例的NRMSE随成像距离D的变化曲线图。
图10是本发明实施例的实验重建图像结果图。
具体实施方式
下面结合本发明实施例中的附图,详细描述本发明的实施过程。
具体实施如示意图1所示,x-y平面内的发射天线阵列的总数量为N,并标记第i个单元的空间位置为i=1,2,…,N。假设每个发射天线阵列的激励电流为其中I0是一个常数。为了实现随机照射,发射天线阵列中的每个发射天线的相位由移相器实现1比特移相。
其中,j表示虚数,k和η0分别为自由空间中的波数和波阻抗,Φ表示所有N个天线单元随机调相值矩阵,N表示天线单元的总数,表示第M次测量下的N个天线单元的随机调相向量,表示第M次测量下的成像区域处的入射场;G为每个发射天线与成像区域上所有点之间的格林函数矩阵;
其中,R(Φ)表示随机调相值矩阵Φ的秩,R(G)表示格林函数矩阵G的秩;
如图2所示,成像过程中将成像区域划分成N′个网格,其中i′表示网格的序号,i′=1,2,…,N′。那么对于第m次随机照射,m=1,2,…,M,对应的传输矩阵H的元素,表达为:
其中,和分别是发射天线、接收天线和成像区域第i′个网格上点之间的传输系数;表示第m次测量下传输矩阵H中成像区域第i′个网格上点对应的元素;表示发射天线阵列中第i个发射天线由移相器控制的随机相位,N表示发射天线的总数。
最终得到成像模型,表达为:
g=Hf
其中,g是测量得到的散射场,H是传输矩阵,f为成像区域的反射率分布。
在M=N=N′时,传输矩阵H为一个方阵,通过利用矩阵求逆f=H-1g重建图像。通过奇异值来评价传输矩阵H的非相关特性,奇异值越大说明矩阵越不相关,所包含的有效信息越多,重建图像的成像结果越好。因此本发明能够通过选取最优的奇异值分布曲线,对于给定的成像分辨率和成像区域大小,获得最优的成像距离、发射天线阵列中的发射天线总数、发射天线之间的间距、随机波束照射的次数的配置,实现最优可定制的成像效果。
不失一般性,假定平面阵列中N=10×10,组阵间距Δr=1λ,成像距离D=5λ,工作频率5.8GHz。为了观察不同随机相位对随机性能的影响,选取了四种不同的1比特随机相位,分别为(0,π),(0,π/2),(0,π/4)和(0,π/8)的高斯随机数。成像区域划分网格数量为N′=20×20,其间距Δr′=0.5λ。随机调相次数M取值为140。根据上述参数,H矩阵一个大小M×N′的矩阵。图3为仿真得到的归一化奇异值曲线,当测量次数大于N时奇异值曲线急剧变化到接近零,表明此时新增的行向量之间具有很强的相关性。因此,选择测量次数M等于天线单元个数N。对于不同的调相数值,0和π的相位切换所对应的奇异值曲线最高,是最佳的1比特调相数值。
在成像系统中,作为指标的成像分辨率Δr′和成像距离D是给定的,待设计的参数包括天线单元数量N和组阵间距Δr。为了验证参数之间的相互关系,图 4为随测量次数和成像距离变化的奇异值分布曲线图,其中N=5×5,Δr=1λ,Δr′=2λ,作为系统设置参数。根据奇异值曲线,对于不同的测量次数,存在最优成像距离。图5为图4在测量次数为25处的二维剖面图。因此,可以依据不同的组阵参数得到相应的成像距离和分辨率的关系曲线,见图6。保持组阵间距为典型的半波长,可以得到不同天线单元个数N和成像距分辨率之间的关系曲线,见图7。
本发明的具体实施例如下:
为了验证所设计的成像系统性能,选取和图6相同的系统设置参数,其中H 矩阵是一个25×25的方阵,因此对成像模型公式(9)采用直接矩阵求逆的方法重建图像。成像结果如图8所示,其中,图8(a)为“T”形的原始图像,有目标散射体的位置反射率为1,空气背景的反射率为0。图8(b)(c)(d)分别为距离5λ、10λ和20λ处成像反演的归一化结果,在近距离(5λ)和远距离 (20λ)时,反演结果模糊,而在最佳成像距离10λ处的成像效果最好。
为了得到成像质量随距离D的变化关系,对其进行参数扫描仿真:成像距离D变化范围5λ-70λ,变化间隔1λ。采用归一化均方根误差NRMSE来衡量反演结果的质量:
其中,fi分别为第i个离散网格处估计值和原始值,N为离散网络的数量,fmax和fmin分别是原始图像的最大值和最小值。图9为不同成像距离D处的 NRMSE曲线,显示在D=10λ处具有最小的误差值,与相关性分析的优化结果一致。
在硬件电路实现中,1比特随机相位调制可以通过采用射频开关切换不同长度的传输线实现。为了实现0°和180°两种相位的切换,两段不同长度的传输线l1和l2的可以分别为λ/2和λ长度。在电路集成化设计中,电路板顶层为5×5 微带贴片天线,底层主要为功分电路、调相电路和控制电路三大功能模块,并且采用矢量网络分析仪得到H矩阵以及散射场的幅度和相位,利用功率放大器来提高信噪比。
图10为不同形状的目标图案的重建结果,显示均能够对原始目标的形状和位置实现正确成像。
由此上述实施例可见,本发明公开了一种基于完全随机波束照射的可定制微波成像系统及算法。包括采用1比特随机调相的阵列天线实现随机照射,并利用矩阵求逆的方法得到重新图像。本发明工作于单一频率,对于给定的成像分辨率和成像距离,依据奇异值分布曲线,给出天线阵列的单元数量和组阵间距以及测量次数,整个成像系统具有结构紧凑,随机性能优良,能量效率高,频谱效率高等优点。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或者替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (5)
1.一种基于随机波束照射的可定制微波成像方法,其特征在于:方法采用以下步骤:
1)搭建微波成像系统进行随机波束照射,采用1比特随机相位调制;
微波成像系统包括多个发射天线和一个接收天线,每个发射天线经移相器连接到发射端口,接收射天线连接到接收端口,多个发射天线和接收天线布置于待测目标的一侧,多个发射天线阵列排布形成发射天线阵列;通过各个移相器随机调节各个发射天线的激励电流相位使得所有发射天线的发出随机电磁波束照射到待测目标所在的成像区域;
2)经过多次随机波束照射,由接收天线采集的接收信号构建传输矩阵;
3)预设成像分辨率、成像区域的大小、成像距离、发射天线阵列中的发射天线总数、发射天线之间的间距,多次随机波束照射后计算传输矩阵的奇异值分布曲线,通过不断调整成像分辨率、成像区域的大小、成像距离、发射天线阵列中的发射天线总数、发射天线之间的间距和随机波束照射的次数,获得最优的奇异值分布曲线,最优的奇异值分布曲线上的所有奇异值之和最大;
4)取最优的奇异值分布曲线下对应的成像分辨率、成像区域的大小、成像距离、发射天线阵列中的发射天线总数、发射天线之间的间距、随机波束照射的次数设置微波成像系统进行成像。
2.根据权利要求1所述的一种基于随机波束照射的可定制微波成像方法,其特征在于:采用1比特随机相位调制的发射天线阵列实现完全随机的波束照射。
3.根据权利要求1所述的一种基于随机波束照射的可定制微波成像方法,其特征在于:所述1比特随机相位调制具体是从0度至360度的相位值中选取任意两个不同的相位值进行调制。
4.根据权利要求1所述的一种基于随机波束照射的可定制微波成像方法,其特征在于:所述1比特随机相位调制通过采用射频开关切换不同长度的传输线实现。
5.根据权利要求1所述的一种基于随机波束照射的可定制微波成像方法,其特征在于:所述的随机波束照射的最少次数与发射天线阵列中的发射天线总数量相同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811073445.7A CN109283529B (zh) | 2018-09-14 | 2018-09-14 | 一种基于完全随机波束照射的可定制微波成像系统及算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811073445.7A CN109283529B (zh) | 2018-09-14 | 2018-09-14 | 一种基于完全随机波束照射的可定制微波成像系统及算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109283529A CN109283529A (zh) | 2019-01-29 |
CN109283529B true CN109283529B (zh) | 2020-07-10 |
Family
ID=65181236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811073445.7A Active CN109283529B (zh) | 2018-09-14 | 2018-09-14 | 一种基于完全随机波束照射的可定制微波成像系统及算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109283529B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114594426B (zh) * | 2022-03-04 | 2024-06-18 | 中国科学技术大学 | 随机辐射天线阵列系统及微波凝视关联成像方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103744078A (zh) * | 2013-12-30 | 2014-04-23 | 中国科学技术大学 | 一种基于不同码速随机跳频的微波凝视关联成像装置 |
CN104360345A (zh) * | 2014-09-29 | 2015-02-18 | 林子怀 | 基于随机天线阵列和微波关联成像原理的高分辨率穿墙成像系统和方法 |
DE102016100367A1 (de) * | 2015-01-12 | 2016-07-14 | Siemens Medical Solutions Usa, Inc. | Spärliche Verfolgung in Schallstrahlintensitätsimpuls-Bildgebung |
CN107678028A (zh) * | 2017-11-22 | 2018-02-09 | 中国科学技术大学 | 低信噪比条件下的微波凝视关联成像方法 |
CN108008380A (zh) * | 2017-11-29 | 2018-05-08 | 中国科学技术大学 | 一种基于卫星编队的微波凝视关联成像方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017105566A1 (en) * | 2015-12-17 | 2017-06-22 | Massachusetts Institute Of Technology | Methods and systems for near-field microwave imaging |
-
2018
- 2018-09-14 CN CN201811073445.7A patent/CN109283529B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103744078A (zh) * | 2013-12-30 | 2014-04-23 | 中国科学技术大学 | 一种基于不同码速随机跳频的微波凝视关联成像装置 |
CN104360345A (zh) * | 2014-09-29 | 2015-02-18 | 林子怀 | 基于随机天线阵列和微波关联成像原理的高分辨率穿墙成像系统和方法 |
DE102016100367A1 (de) * | 2015-01-12 | 2016-07-14 | Siemens Medical Solutions Usa, Inc. | Spärliche Verfolgung in Schallstrahlintensitätsimpuls-Bildgebung |
CN107678028A (zh) * | 2017-11-22 | 2018-02-09 | 中国科学技术大学 | 低信噪比条件下的微波凝视关联成像方法 |
CN108008380A (zh) * | 2017-11-29 | 2018-05-08 | 中国科学技术大学 | 一种基于卫星编队的微波凝视关联成像方法及系统 |
Non-Patent Citations (1)
Title |
---|
太赫兹孔径编码成像研究综述;邓斌 等;《红外与毫米波学报》;20170630;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109283529A (zh) | 2019-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5681196B2 (ja) | 稀薄分布のエレメントアレイを有するアンテナ | |
CN113300119B (zh) | 一种用于圆极化波束赋形的透射型超表面及设计方法 | |
CN109768385B (zh) | 一种超宽带大容量太赫兹超材料随机辐射天线 | |
CN113708077A (zh) | 一种基于幅度和相位调制的波束赋形超表面结构及方法 | |
CN113131225A (zh) | 一种基于子阵阵因子增益最大化的不规则天线阵优化方法 | |
CN105024166B (zh) | 基于子阵的平面阵列天线方向图综合方法 | |
CN109283529B (zh) | 一种基于完全随机波束照射的可定制微波成像系统及算法 | |
Li et al. | Pattern synthesis of conformal arrays by a modified particle swarm optimization | |
Dahl et al. | Comparison of virtual arrays for MIMO radar applications based on hexagonal configurations | |
CN104063426B (zh) | 一种面向辐射和散射的有源相控阵天线结构公差的快速确定方法 | |
CN112968299A (zh) | 一种基本单元、超表面及基于信息超表面的到达角估计方法 | |
CN105158736B (zh) | 一种mimo雷达发射方向图与波形设计方法 | |
CN115332762A (zh) | 一种低互相干的非均匀随机辐射阵列口径天线及方法 | |
CN215834727U (zh) | 一种多入多出天线子阵及稀疏化天线 | |
CN106338656B (zh) | 大型地面天线罩空间肋阵的散射云图获取方法 | |
CN111142164B (zh) | 一种柱面雷达成像系统 | |
Zhou et al. | Efficient optimization and realization of a shaped-beam planar array for very large array application | |
Li et al. | Binary array beam forming with simulated annealing method | |
Kim et al. | A Hadamard matrix feed network for a dual-beam forming array antenna | |
Bourges et al. | High-frequency surface wave radar based on a sea floating antenna concept | |
CN113206388B (zh) | 基于相位调制有源频率选择表面的成像系统及其成像方法 | |
CN112526510A (zh) | 基于方向图分集的单通道角度超分辨方法及系统 | |
Azzam et al. | Phased-only pattern synthesis for sunflower and conformal arrays with shaped radiation patterns | |
CN218586361U (zh) | 一种可扩展子阵级非周期稀布阵列天线 | |
Liu et al. | Recent development in nonuniformly spaced array synthesis methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |