CN102253376B - 一种基于二维微波成像的低散射共形天线rcs测试方法 - Google Patents

一种基于二维微波成像的低散射共形天线rcs测试方法 Download PDF

Info

Publication number
CN102253376B
CN102253376B CN2011100944616A CN201110094461A CN102253376B CN 102253376 B CN102253376 B CN 102253376B CN 2011100944616 A CN2011100944616 A CN 2011100944616A CN 201110094461 A CN201110094461 A CN 201110094461A CN 102253376 B CN102253376 B CN 102253376B
Authority
CN
China
Prior art keywords
rcs
dimensional
mrow
conformal antenna
theta
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.)
Expired - Fee Related
Application number
CN2011100944616A
Other languages
English (en)
Other versions
CN102253376A (zh
Inventor
李南京
郭淑霞
胡楚锋
刘宁
刘琦
李瑛�
陈卫军
杨博
张麟兮
王保平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN2011100944616A priority Critical patent/CN102253376B/zh
Publication of CN102253376A publication Critical patent/CN102253376A/zh
Application granted granted Critical
Publication of CN102253376B publication Critical patent/CN102253376B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种基于二维微波成像的低散射共形天线RCS测试方法,其特征在于步骤如下:对装机状态下的金属蒙皮及共形天线进行RCS测试,得到二维微波像;对新的二维微波像进行二维傅立叶变换,得到目标谱域的数据;对一个RCS已知的金属球进行RCS测试,得到二维微波像对二维微波像进行二维傅立叶变换,得到金属球的散射场随频率和角度变化的数据G0(f,θ),最终可得:共形天线的RCS=共形天线数据G1(f,θ)-金属球数据G0(f,θ)+金属球RCS。本发明提出的方法,基于二维微波成像技术获取装机状态下共形天线的RCS测试方法,经过实践,所得的效果良好,具有极大的工程推广价值。

Description

一种基于二维微波成像的低散射共形天线RCS测试方法
技术领域
本发明涉及一种基于二维微波成像的低散射共形天线RCS测试方法,是借助二维微波成像技术,分离和提取出装机状态下共形天线的RCS。这是一种低散射目标RCS测量的新方法,属于微波技术领域。
背景技术
雷达是一种有效的远程探测设备,它通过发射电磁波并且接收目标反射的回波来精确定位,雷达所收到的回波大小跟目标的有效截面积有关,我们把这个有效截面积称为雷达散射截面(简称RCS)。
随着电子战的迅猛发展,隐身武器的大规模使用是现代战争的一个显著特征。例如美国的第三代隐身战斗机F-117A,其雷达散射截面仅为0.001-0.01m2。由于我国的隐身飞行器处于起步阶段,因而低散射目标的测试成为目前RCS测量遇到的新问题,尤其是如何获取装机状态下共形天线所贡献的RCS,成为人们所关注的焦点。
如图一所示,在机体表面,天线通过一块介质板与机体形成共形,由于入射电场E入射角小,天线安装在机体内侧,天线大部分结构(特别是天线结构的外侧)均未暴露给入射电场。因此,天线的RCS贡献主要由安装结构内腔产生。因此,如何界定嵌入机体区域天线的RCS是我们必须解决的问题。
目前国内外通常采用的方法如图二所示,天线发射步进频率宽带信号,对回波的频率响应作逆傅立叶变换后得到目标的一维距离像,即目标的散射中心随距离的分布情况,不同的强散射点会在不同的距离上体现出峰值,通过时域加门截取目标区的回波,定标后,最终得到目标中心频点的RCS。
该方法是通过成一维像后提取RCS,但在一维像合成的过程中,实质上是将整个横向区域内的散射信息全部叠加到一起,如图一中的一维像,除了包含目标的散射信息,还包含着B1和B2区域的散射,如果B1和B2区域的散射大于目标区域的散射,那么通过一维像提取的目标RCS具较大的误差。
因此,低散射天线装机状态的RCS特性分离和提取存在以下困难:
1)天线嵌入机体结构内,给常规的RCS测试带来了新问题,难以区分;
2)低散射天线的RCS值与机体结构的RCS值在同一量级,必须区分出来,才能获得共形天线可信的RCS。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于二维微波成像的低散射共形天线RCS测试方法。
技术方案
一种基于二维微波成像的低散射共形天线RCS测试方法,其特征在于步骤如下:
步骤1:对装机状态下的金属蒙皮及共形天线进行RCS测试,得到二维微波像
Figure BSA00000474155400021
其中,
Figure BSA00000474155400022
就是估计的x-y平面二维图像。k是波数,k=4πf/c,与频率f相关,kmin与kmax对应扫频的最小频率和最大频率,θ是转台转角,θmin与θmax分别对应最小和最大转角,G(k,θ)是测试频点和转角的目标回波复数数据;
步骤2:在二维微波像上保留共形天线所在几何区域的反射率分布,区域外的反射率设置为零,形成新的二维微波像:
Figure BSA00000474155400023
所述几何区域的半径等于共形天线最大半径a;
步骤3:对新的二维微波像进行二维傅立叶变换,得到目标谱域的数据:
G 1 ( K x , K y ) = ∫ y 1 y 2 ∫ x 1 x 2 g ^ 1 ( x , y ) exp [ - ( j 2 π K x x + j 2 π K y y ) ] dxdy
Kx,Ky分别是谱域的横轴和纵轴,x1、x2,y1、y2分别对应目标像
Figure BSA00000474155400032
所在x-y平面的上下限范围。对G1(Kx,Ky)进行均匀插值,得到随频率和角度变化的谱域数据G1(f,θ),其中,
Figure BSA00000474155400033
θ=tan-1(Ky/Kx);
步骤4:对一个RCS已知的金属球进行RCS测试,得到二维微波像
Figure BSA00000474155400034
对二维微波像进行二维傅立叶变换,得到金属球的散射场随频率和角度变化的数据G0(f,θ),最终可得:共形天线的RCS=共形天线数据G1(f,θ)-金属球数据G0(f,θ)+金属球RCS。
有益效果
本发明提出的一种基于二维微波成像的低散射共形天线RCS测试方法,基于二维微波成像技术获取装机状态下共形天线的RCS测试方法,经过实践,所得的效果良好,具有极大的工程推广价值。
经过试验的测试目标为一个直径为75mm的金属球,采用直径为150mm的金属球定标,测试频段为8-12GHz,转动角度为360度,选取80度的数据进行二维微波成像,比较常规方法和本方法的RCS。
图4给出了通过常规RCS测量获得的结果,从图中可以看出随着角度的变化,金属球RCS的起伏约为0.7dB,并且有较大的毛刺,反映了测试背景对目标的影响。通过二维微波成像后,提取和分离出的金属球RCS如图5所示,金属球的起伏约为0.5dB,并且曲线的毛刺明显减小,说明了该方法的有效性。
附图说明
图1:共形天线嵌入机体安装结构示意图
图2:步进频率法测量RCS示意图
图3:本方法的信号处理与RCS提取示意图
图4:金属球RCS常规方法与二维成像后提取RCS的比较
图5:模拟装机状态下的低散射共形天线
图6:系统组成图
图7:装机状态下共形天线的二维成像
图8:天线反射率分布的波谱图
图9:几种结果的比较
具体实施方式
现结合实施例、附图对本发明作进一步描述:
如图5所示,测量装机状态下的共形天线,所采用的二维微波成像系统如图6所示,系统发射步进频率信号,信号经过功率放大器放大后直接由标准增益天线发射出去,回波信号通过另一个相隔很近的标准增益天线送入矢量网络分析仪接收端,转台旋转,在间隔一定的角度下获取回波数据进行二维成像。
实验测量步骤如下:
1)测量空暗室,进行时域对消,减小杂波影响。
2)测量目标,将其变换到时域,用距离波门截取高于背景的目标区,然后反变换回频域,记录数据。
3)测量定标体,在相同距离波门内记录其频域数据。
测量得到的场与像的关系如式(1)所示:
g ^ ( x , y ) = ∫ θ min θ max ∫ k min k max kG ( k , θ ) exp [ j 2 πk ( y cos θ - x sin θ ) ] dkdθ - - - ( 1 )
式中x-y坐标系是固定于目标上的一组坐标,坐标原点在O点,它随目标变化而变化。
Figure BSA00000474155400051
即所要恢复目标的像。
由于式(1)中积分式的积分限不满足IFFT条件,故在实际的算法实现中需要对k做频移kmin,相当于解调过程,令B′=kmax-kmin,B′是空间频率k的带宽,因此有
P θ ( l ) = ∫ 0 B ′ ( k + k min ) G ( k + k min , θ ) exp ( j 2 πkl ) dk - - - ( 2 )
g ^ ( x , y ) = ∫ θ min θ max P θ ( l ) exp ( j 2 π k min l ) dθ - - - ( 3 )
l=ycosθ-xsinθ
(4)
系统采用的是步进频率信号,所以频率点是离散的,根据采样的频率点数N,可以对k进行离散化,令k=n·B′/N,n=0,1,2,…,N-1。于是式(2)改写成:
P θ ( l ) = Σ k = 0 ( N - 1 ) B ′ / N ( B ′ N · n + k min ) G ( n , θ ) e j 2 π B ′ N · n · l - - - ( 5 )
投影线l按照距离分辨率进行等间距离散化lm=m/B′,m=0,1,…,N-1,令Gθ(n)=(B′·n/N+kmin)G(n,θ),则式(5)可表示为:
P θ ( l m ) = Σ n = 0 N - 1 G θ ( n ) e j 2 π m N n = IFFT [ G θ ( n ) ] - - - ( 6 )
Pθ(lm)是在不同的lm点上所得的投影值,但是积分过程中的投影线l是随θ角变化,对于空间任一点的
Figure BSA00000474155400056
每个θ对应不同的l,因此Pθ(l)是在离散的l处得到的投影值,所以需要通过对Pθ(lm)进行插值来获得Pθ(l),插值公式如下:
Pθ(l)=Pθ(lm-1)+(l-lm-1)tgα     (7)
tgα=[Pθ(lm)-Pθ(lm-1)]/(lm-lm-1)(8)
做完插值后,便可对角度进行积分,公式如下:
g ^ ( x , y ) = ∫ θ min θ max P θ ( l ) exp ( j 2 π k min l ) dθ = Σ θ P θ ( l ) e j 2 π k min l - - - ( 9 )
通过信号空间变换、IFFT、插值、角度积分四步,便完成了成像的处理过程。
二维成像的结果如图7(a)所示,除了天线的反射率分布之外,二维像还包括了载体目标的反射率分布,然后采用一个二维窗,选取出其中天线的反射率分布,如图7(b)所示。
对提取出的天线反射率分布进行傅立叶变换,得到天线随角度和频率变化的散射波谱图,从中提取出9.375GHz频率的曲线,即得到图9中虚线给出的曲线。对天线仰角5度进行测试的效果分析如下:
红线是常规测试得到的RCS(含工装的全域范围),蓝线是从目标的二维像合成的RCS曲线(含工装的全域范围),两者在全角域范围吻合良好;虚线是从二维像中截取天线区域后,再合成获得的天线RCS贡献,它比前者约低3~5dB,通过这种处理方式能进一步降低外部吸波材料工装的RCS影响。因天线外有金属工装过渡,外部吸波材料工装的制作要求可以大大降低。

Claims (1)

1.一种基于二维微波成像的低散射共形天线RCS测试方法,其特征在于步骤如下:
步骤1:对装机状态下的金属蒙皮及共形天线进行RCS测试,得到二维微波像 其中, 
Figure FSB00000890899500012
就是估计的x-y平面二维图像,k是波数,k=4πf/c,与频率f相关,kmin与kmax对应扫频的最小频率和最大频率,θ是转台转角,θmin与θmax分别对应最小和最大转角,G(k,θ)是测试频点和转角的目标回波复数数据;
步骤2:在二维微波像上保留共形天线所在几何区域的反射率分布,区域外的反射率设置为零,形成新的二维微波像:
Figure FSB00000890899500013
所述几何区域的半径等于共形天线最大半径a;
步骤3:对新的二维微波像进行二维傅立叶变换,得到目标谱域的数据:
Figure FSB00000890899500014
Kx,Ky分别是谱域的横轴和纵轴,x1、x2,y1、y2分别对应目标像 
Figure FSB00000890899500015
所在x-y平面的上下限范围,对G1(Kx,Ky)进行均匀插值,得到随频率和角度变化的谱域数据
G1(f1,θ1),其中,
Figure FSB00000890899500016
步骤4:对一个RCS已知的金属球进行RCS测试,得到二维微波像 
Figure FSB00000890899500017
对二维微波像 
Figure FSB00000890899500018
进行二维傅立叶变换,得到金属球的散射场随频率和角度变化的数据G0(f0,θ0),最终可得:共形天线的RCS=共形天线谱域数据G1(f1,θ1)-金属球数据G0(f0,θ0)+金属球RCS。 
CN2011100944616A 2011-04-14 2011-04-14 一种基于二维微波成像的低散射共形天线rcs测试方法 Expired - Fee Related CN102253376B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100944616A CN102253376B (zh) 2011-04-14 2011-04-14 一种基于二维微波成像的低散射共形天线rcs测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100944616A CN102253376B (zh) 2011-04-14 2011-04-14 一种基于二维微波成像的低散射共形天线rcs测试方法

Publications (2)

Publication Number Publication Date
CN102253376A CN102253376A (zh) 2011-11-23
CN102253376B true CN102253376B (zh) 2012-11-21

Family

ID=44980742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100944616A Expired - Fee Related CN102253376B (zh) 2011-04-14 2011-04-14 一种基于二维微波成像的低散射共形天线rcs测试方法

Country Status (1)

Country Link
CN (1) CN102253376B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102590792B (zh) * 2012-02-29 2014-10-22 中国电子科技集团公司第二十二研究所 基于矢量网络分析仪的l波段岸基海杂波观测系统
CN102590795B (zh) * 2012-02-29 2014-03-26 中国电子科技集团公司第二十二研究所 基于矢量网络分析仪的微波散射特性测试系统
CN102608591A (zh) * 2012-03-30 2012-07-25 西北工业大学 一种基于相位修正的rcs外推方法
CN102967855B (zh) * 2012-12-12 2014-04-09 北京航空航天大学 异地连续定标rcs测量中改进背景相减技术的信号处理方法
CN103792520B (zh) * 2013-01-30 2017-05-17 承德石油高等专科学校 太赫兹目标散射特性测试装置及测试方法
CN104199026B (zh) * 2014-08-28 2016-08-24 中国科学院电子学研究所 基于线迹扫描二维近场成像的反向散射截面测量方法
CN104613852B (zh) * 2015-01-09 2017-06-16 北京环境特性研究所 一种用于rcs测量的弧形导轨定位方法及系统
CN107315846B (zh) * 2016-08-29 2020-07-31 南京航空航天大学 一种快速分析宽带rcs的方法
CN107394416B (zh) * 2017-06-01 2020-10-23 西南电子技术研究所(中国电子科技集团公司第十研究所) 自适应改变辐射和散射特性的智能蒙皮天线
CN107995765B (zh) * 2017-10-19 2020-09-04 上海无线电设备研究所 一种等离子腔测试背景消除方法
CN108693506B (zh) * 2018-04-02 2020-03-10 北京环境特性研究所 一种对天线的侧向雷达散射截面积进行测试的载体
CN108647183B (zh) * 2018-04-02 2021-09-24 北京环境特性研究所 基于压缩感知的复rcs数据插值方法
CN110441745B (zh) * 2019-08-16 2021-04-30 北京环境特性研究所 一种基于宽带雷达俯视测量目标rcs的方法和系统
CN110850384B (zh) * 2019-11-04 2020-06-23 北京航空航天大学 一种基于扫频数据产生宽带去斜回波的方法
CN112541279B (zh) * 2020-12-25 2024-03-19 中国航空工业集团公司沈阳飞机设计研究所 一种飞机近场散射特征构造方法

Non-Patent Citations (12)

* Cited by examiner, † Cited by third party
Title
一种新的近距离微波成像系统的建模与仿真;胡楚锋 等;《科学技术与工程》;20070531;第7卷(第10期);2358-2360 *
刘密歌 等.基于矢量网络分析仪的RCS测量系统及应用.《电子测量与仪器学报》.2007,第21卷(第1期),82-85.
刘密歌 等.微波暗室中两种RCS测量系统的比较.《计算机测量与控制》.2007,(第3期),300-301,341.
基于扫频时域法测量的RCS外推技术研究;李南京 等;《微波学报》;20070831;第23卷(第4期);24-27 *
基于矢量网络分析仪的RCS测量系统及应用;刘密歌 等;《电子测量与仪器学报》;20070228;第21卷(第1期);82-85 *
室内UHF波段高精度RCS测试方法及定标技术研究;李南京 等;《系统工程与电子技术》;20090930;第31卷(第9期);2114-2116 *
微波暗室中两种RCS测量系统的比较;刘密歌 等;《计算机测量与控制》;20071231(第3期);300-301,341 *
李南京 等.基于扫频时域法测量的RCS外推技术研究.《微波学报》.2007,第23卷(第4期),24-27.
李南京 等.室内UHF波段高精度RCS测试方法及定标技术研究.《系统工程与电子技术》.2009,第31卷(第9期),2114-2116.
李南京 等.远场RCS的精确测试方法研究.《现代雷达》.2006,第28卷(第8期),70-73.
胡楚锋 等.一种新的近距离微波成像系统的建模与仿真.《科学技术与工程》.2007,第7卷(第10期),2358-2360.
远场RCS的精确测试方法研究;李南京 等;《现代雷达》;20060831;第28卷(第8期);70-73 *

Also Published As

Publication number Publication date
CN102253376A (zh) 2011-11-23

Similar Documents

Publication Publication Date Title
CN102253376B (zh) 一种基于二维微波成像的低散射共形天线rcs测试方法
CN107942330B (zh) 一种基于等离子体近场测试的雷达散射特征数据提取方法及系统
CN108983208B (zh) 一种基于近场稀疏成像外推的目标rcs测量方法
CN107783092B (zh) 基于链条关系式的近场后向rcs测量系统及方法
CN106569191A (zh) 一种利用高分辨率成像获取目标rcs的方法
CN109932719A (zh) 基于sar成像的rcs高精度测量方法
Setsu et al. Super-Resolution Doppler Velocity Estimation by Kernel-Based Range–$\tau $ Point Conversions for UWB Short-Range Radars
Tulgar et al. Improved pencil back-projection method with image segmentation for far-field/near-field SAR imaging and RCS extraction
Li et al. NUFFT-based near-field imaging technique for far-field radar cross section calculation
Piotrowsky et al. Antenna pattern characterization with an industrial robot assisted fmcw radar system
Dürr et al. Highly efficient angular array calibration based on the modal wave expansion technique
Chen et al. Direct wave removal in anechoic chamber range imaging from planar scanned data
CN112782692A (zh) 一种基于弧形mimo线阵扫描的成像系统及成像方法
CN112965061B (zh) 一种基于柱面mimo面阵的成像系统及成像方法
Massaloux et al. Indoor 3D spherical near field RCS measurement facility: Localization of scatterers
Noh et al. RCS feature extraction using discretized point scatterer with compressive sensing
Morales et al. Comparison of echo reduction techniques for one-single cut antenna measurements
CN107765230B (zh) 链条关系式在近场测量系统的近场到远场变换中的应用方法
Sun et al. Antenna planes based wall-clutter mitigation in through-wall-imaging applications
Chen et al. Estimation of pose angle for trihedral in ultrawideband virtual aperture radar
RU2717823C1 (ru) Способ измерения угла места воздушного объекта в метровом диапазоне электромагнитных волн
Wang et al. Comparison of EMC Chamber Debugging Techiniques above 1GHz
CN114415140B (zh) 基于近场平面扫描三维波数域成像的目标rcs测量方法
CN113589280B (zh) 一种频域加窗单视快速雷达成像优化解析方法
Kabourek Objects characterization by means of wideband signals

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121121

Termination date: 20150414

EXPY Termination of patent right or utility model