
Citation: | LIU Cheng, ZHAO Zhan, DU Lidong, FANG Zhen. Method of Measuring Two-dimensional Wind Based on Diametrical Pressure Differences Developed by Flow Around Cylinder[J]. Journal of Electronics & Information Technology, 2017, 39(3): 737-742. doi: 10.11999/JEIT160468 |
弹道导弹已成为现代战争中的主导武器,而导弹飞行中段往往会伴有诱饵,弹头和诱饵等目标微动特性因其质量分布不同而存在明显差异,因此估计目标的微动参数成为导弹防御系统目标识别重要技术[1]。随着高分辨雷达技术和现代信号处理技术的不断发展,空间目标的微动参数估计已由基于RCS序列的目标微动参数估计逐渐发展到基于微多普勒序列[2,3]、高分辨距离像序列的目标参数估计等技术。结合相位测距等精细化处理技术,利用宽带雷达高分辨序列距离像的参数估计,近年来获得了广泛关注[4,5]。
目标微动是在3维空间进行的,关键散射点瞬时位置蕴含丰富的目标运动几何信息,为3维干涉测量实现目标微动参数估计提供基础[6—9]。本文考虑分布式相参雷达体制多天线接收系统的空间分集[10],利用接收波程相位差对目标瞬时空间位置进行测量,进而获得目标3维运动几何参数估计。相较于高分辨距离像序列的参数估计,宽带3维干涉测量技术可为微动目标分析在宽带测距的基础上提供干涉测角信息,为目标微动几何参数估计提供更多自由度和更高稳健性。文献[11]利用目标微动在宽带系统下回波信号的调制特征差异来区分各散射点的回波信号,通过干涉处理重构微动目标的3维坐标。文献[12]采用对目标距离-慢时间像进行干涉处理的方法重构出弹道目标的真实3维坐标并有效估计出目标的微动参数与结构参数。考虑到实际应用中,目标会不可避免地偏离电轴,即在斜视情况下[13,14],上述文献算法重构目标3维坐标将发生坐标扭曲,对目标几何运动参数估计应用受限。
针对以上问题,本文提出了一种斜视校正干涉测量微动几何参数估计方法。首先分别对L型分布式相参雷达体制下各天线接收距离-慢时间回波进行干涉处理,通过求解2元2次非线性方程组并利用坐标变换来校正斜视引起的误差,重构出目标真实的3维运动轨迹;接着利用
L型分布式相参雷达与弹道目标几何关系如图1所示,
对目标各散射点回波进行解线频调处理,并消除包络斜置、残余视频项并距离维脉冲压缩后,天线
SAd(f,tm)=σkTpsinc{Tp[f+2μcRΔAk(tm)]}⋅exp[−j4πλRΔAk(tm)] | (1) |
其中,
ΔφAB=angle(S∗Ad(f,tm)⋅SBd(f,tm))=2πλ(RΔAk−RΔBk) | (2) |
ΔφAC=angle(S∗Ad(f,tm)⋅SCd(f,tm))=2πλ(RΔAk−RΔCk) | (3) |
实际3维干涉测量技术应用中,目标通常会偏离电轴,即斜视测量。当斜视达到一定程度时,两对天线测出的坐标信息不能直接对应散射点真实坐标值,并且散射点到发射天线的斜距也不能近似为其
以天线对
RAk−RBk=2L(X+xk)−L2RAk+RBk | (4) |
RAO−RBO=2LX−L2RAO+RBO | (5) |
进一步可推导出
ΔRAB=RΔAk−RΔBk=(RAk−RAO)−(RBk−RBO)=(RAk−RBk)−(RAO−RBO)=2Lxk(RAO+RBO)(RAk+RBk)(RAO+RBO)+2LX(RAO+RBO−RAk−RBk)(RAk+RBk)(RAO+RBO)−L2(RAO+RBO−RAk−RBk)(RAk+RBk)(RAO+RBO) | (6) |
其中,
根据上文分析,斜视时由于误差相位项不可忽略,下面介绍通过求解一个2元2次非线性方程组校正斜视引起的误差。
类似式(6),有
ΔRAC=2Lxk(RAO+RCO)(RAk+RCk)(RAO+RCO)+2LX(RAO+RCO−RAk−RCk)(RAk+RCk)(RAO+RCO)−L2(RAO+RCO−RAk−RCk)(RAk+RCk)(RAO+RCO) | (7) |
式(6)和式(7)中可做如下近似:
f1=RAO+RBO−RAk−RBk=LykR0−1R0(2xkX+2ykY+2zkZ+x2k+y2k+z2k)f2=RAO+RCO−RAk−RCk=LxkR0−1R0(2xkX+2ykY+2zkZ+x2k+y2k+z2k)} | (8) |
将式(6)—式(8)代入式(2)和式(3)中并作近似后可得
ΔφAB=2πLλR0yk+2Y−L2λR20Lπf1ΔφAC=2πLλR0xk+2X−L2λR20Lπf2} | (9) |
式(8)中
斜视时雷达视线与目标纵轴存在较大夹角,散射点到天线的斜距不能直接近似为
首先以
[wuv]=M1⋅M2[xyz] | (10) |
其中
z=(v−ysinβy)(1/cosβ) | (11) |
根据式(11)求得斜视状态下散射点
由于噪声和距离分辨率的影响,上述步骤获得的散射点重构3维运动轨迹曲线会出现波动和阶跃现象,与真实运动轨迹偏差过大,若利用其直接进行微动参数估计,会导致较大的参数估计误差。因此需要对重构3维运动轨迹进行平滑滤波,以减少误差带来的不利影响。微动目标散射中心运动模型不易建立,本文用CV模型进行近似。本文目标的角速度和角加速度变化较小,利用
由上节重构目标3维运动轨迹可以看出,目标的微动几何参数与3维运动轨迹有直接关系。本节在获取目标真实3维坐标曲线的基础上首先估计进动频率,然后结合重构运动轨迹3维图像的几何关系估计进动角和几何参数。
对于进动目标而言,锥顶散射点的微动变化曲线为正弦形式,频域为单频信息,频点位置由进动频率确定,可直接对目标进动频率进行估计。
对上一节估计所得的锥顶散射点坐标曲线
利用3维干涉测量技术对弹道目标回波进行处理后,可得到在观测时间内与目标真实空间位置一致的各散射点3维运动轨迹图像,利用该图像可以提取目标进动角及几何参数[12]。
根据图2的几何关系,设向量
\theta = \arccos \frac{{{n} \cdot {OK}}}{{\left\| {n} \right\| \cdot \left\| {{OK}} \right\|}} | (12) |
半锥角
\eta = \arccos \frac{{{KP} \cdot {KO}}}{{\left\| {{KP}} \right\| \cdot \left\| {{KO}} \right\|}} | (13) |
根据几何关系,锥体高度
H = \left\| {{KP}} \right\| \cdot \cos \eta | (14) |
根据勾股定理即可得出底面半径
r = \sqrt {{{\left\| {{KP}} \right\|}^2} - {H^2}} | (15) |
至此完成斜视测量下的宽带雷达3维干涉测量目标微动参数估计。为了更清楚地说明本文算法,具体流程如图4所示。
本实验所估计弹道目标形状如图2,主要目标参数如表1所示,主要雷达参数如表2所示。
参数 | 数值 |
锥体高度H (m) | 0.96 |
质心到锥顶距离d (m) | 0.64 |
底面半径r (m) | 0.25 |
自旋频率fs (Hz) | 2 |
锥旋频率fz (Hz) | 3 |
进动角θ (°) | 10 |
参数 | 数值 |
载频f0 (GHz) | 10 |
带宽B (GHz) | 1 |
脉冲宽度tp (μs) | 10 |
脉冲重复周期prf (Hz) | 100 |
驻留时间T (s) | 1 |
基线长L (m) | 200 |
初始时刻目标中心在雷达坐标系下的坐标为(2 km, 2 km, 10 km)。为对比本文所提方法重构3维运动轨迹的精确性,在本文设置的实验条件下,使用参考文献[12]中的估计方法,直接通过干涉测量求解各散射点的
从图6可以看出,由于斜视角的存在,利用文献[12]中所提方法直接通过干涉测量技术重构的3维坐标与真实值有较大偏差。在相同的实验条件下,利用本文算法,在干涉测量的基础上通过求解式(9)和式(11)完成斜视校正,得到的3维重构坐标如图7所示。可以看出,图7中的重构坐标更接近理论值。
由于图7中重构的3维坐标存在波动和阶跃现象,不利于后续微动参数估计。故采用
结合瞬时3维坐标估计值
A = \left( {1 - \frac{{\displaystyle\sum\limits_{{t_m}} {\left| {{F_{\rm r}}({t_m}) - {F_{\rm e}}({t_m})} \right|} }}{{\displaystyle\sum\limits_{{t_m}} {\left| {{F_{\rm r}}({t_m})} \right|} }}} \right) \times 100\% | (16) |
实验是在信噪比为20 dB条件下进行,所得实验结果相似度记录在表3中。
锥顶散射点x | 锥顶散射点y | 锥顶散射点z | 锥底散射点x | 锥底散射点y | 锥底散射点z | |
文献[12]干涉测量 | 68.26 | 60.42 | 88.72 | 62.93 | 64.91 | 82.73 |
斜视校正干涉测量 | 88.56 | 90.41 | 96.59 | 90.84 | 91.24 | 87.28 |
从表3和图8中可以看出,相较于文献[12]中所提算法,本文算法能够有效校正斜视带来的重构坐标误差,有较高的曲线估计相似度。
对锥顶散射点的坐标曲线进行傅里叶变换即可求得目标进动频率
进动频率{f_c}\left( {\rm Hz} \right) | 进动角\theta \left( {^ \circ } \right) | 半锥角\eta \left( {^ \circ } \right) | 锥体高度H\left( {\rm m} \right) | 底面半径r\left( {\rm m} \right) | |
真实值 | 3 | 10.00 | 14.59 | 0.96 | 0.25 |
平滑处理估计值 | 3 | 9.60 | 15.25 | 0.95 | 0.26 |
未平滑处理估计值 | 3 | 12.49 | 18.17 | 0.87 | 0.31 |
上述实验在SNR为20 dB的实验条件下得出,为验证方法的稳定性,设置SNR=–15~20 dB进行噪声稳健性分析,并以式(17)定义估计相对误差
Q = \frac{{\left| {a - \bar a} \right|}}{a} \times 100\% | (17) |
其中,
从表4可以看出,各参数估计值与实验模型参数设置值都较为接近,从图10可以看出,SNR较低时,由于
本文以L型分布式相参雷达作用在锥体弹道目标为模型,利用斜视3维干涉测量技术处理回波信号,通过求解一个2元2次非线性方程组和坐标变换重构斜视情况下散射点3维运动轨迹,并在此基础上通过滤波和优化求解实现目标微动几何参数估计。实验结果表明该方法通过斜视3维干涉测量可有效重构散射点3维空间位置,并且具有较高的稳健性。但目前本文仅针对单个光滑锥体目标在进动模式下的参数估计,而实际情况中目标的个数、形状和运动方式会更加复杂,重构运动轨迹和参数的估计也会更加困难,后续工作中可对上述问题进行进一步研究。
孙学金, 王晓蕾, 李浩, 等. 大气探测学[M]. 北京: 气象出版社, 2009: 172-247.
|
SUN Xuejin, WANG Xiaolei, LI Hao, et al. Atmospheric Observation[M]. Beijing: China Meteorological Press, 2009: 172247.
|
BRUSCHI P, DEI M, and PIOTTO M. A low-power 2-D wind sensor based on integrated flow meters[J]. IEEE Sensors Journal, 2009, 9(12): 16881696. doi: 10.1109/JSEN.2009. 2030652.
|
ZHU Yanqing, CHEN Bei, QIN Ming, et al. 2-D micromachined thermal wind sensorsA review[J]. IEEE Internet of Things Journal, 2014, 1(3): 216232. doi: 10.1109/ JIOT.2014.2319296.
|
VERESHCHAGINA E, TIGGELAAR R M, SANDERS R G P, et al. Low power micro-calorimetric sensors for analysis of gaseous samples[J]. Sensors and Actuators B: Chemical, 2015, 206: 772787. doi: 10.1016/j.snb.2014.08.077.
|
LIU Haobing, LIN Nay, PAN Shanshan, et al. High sensitivity, miniature, full 2-D anemometer based on MEMS hot-film sensors[J]. IEEE Sensors Journal, 2013, 13(5): 19141920. doi: 10.1109/JSEN.2012.2236014.
|
KOWALSKI L, JIMENEZ V, DOMINGUEZ-PUMAR M, et al. Low pressure spherical thermal anemometer for space missions[C]. 2013 IEEE Sensors, Baltimore, Maryland, USA, 2013: 14. doi: 10.1109/ICSENS.2013.6688493.
|
KOWALSKI L, ATIENZA M T, GORRETA S, et al. Spherical wind sensor for the atmosphere of Mars[J]. IEEE Sensors Journal, 2016, 16(7): 18871897. doi: 10.1109/JSEN. 2015.2509168.
|
DONG Ziqiang, CHEN Jingjing, QIN Yukun, et al. Fabrication of a micromachined two-dimensional wind sensor by Au-Au wafer bonding technology[J]. Journal of Microelectromechanical Systems, 2012, 21(2): 467475. doi: 10.1109/JMEMS.2011.2179014.
|
ZHU Yanqing, CHEN Bei, GAO Di, et al. A robust and low-power 2-D thermal wind sensor based on a glass-in-silicon reflow process[J]. Microsystem Technologies, 2015, 22(1): 151162. doi: 10.1007/s00542-015-2423-9.
|
GEORGIOU D P and KYPROS F M. Fabrication and calibration of a sub-miniature 5-hole probe with embedded pressure sensors for use in extremely confined and complex flow areas in turbomachinery research facilities[J]. Flow Measurement and Instrumentation, 2014, 39: 5463. doi: 10.1016/j.flowmeasinst.2014.07.005.
|
PIOTTO M, DEI M, PENNELLI G, et al. A miniaturized 2D solid state anemometer based on thermal flow sensors[J]. Procedia Chemistry, 2009, 1(1): 14631466. doi: 10.1016/j. proche.2009.07.365.
|
PIOTTO M, PENNELLI G, and BRUSCHI P. Fabrication and characterization of a directional anemometer based on a single chip MEMS flow sensor[J]. Microelectronic Engineering, 2011, 88(8): 22142217. doi: 10.1016/j.mee.2010.11.009.
|
张兆顺, 崔桂香. 流体力学[M]. 北京: 清华大学出版社, 1999: 167168.
|
ZHANG Zhaoshun and CUI Guixiang. Hydromechanics[M]. Beijing: Tsinghua University Press, 1999: 167168.
|
SCHLICHTING H. Boundary-Layer Theory[M]. 7th Edition, New York: McGraw-Hill, 1979: 1223.
|
ZURELL Cory. Aerodynamics of a circular cylinder inclined to airflow and wind-induced vibrations of dry, inclined cables at high wind speeds[D]. [Ph.D. dissertation], University of Ottawa, 2004: 2225.
|
CARE I and ARENAS M. On the impact of anemometer size on the velocity field in a closed wind tunnel[J]. Flow Measurement and Instrumentation, 2015, 44: 210. doi: 10.1016/j.flowmeasinst.2014.11.007.
|
CHEN Wenli, LI Hui, and HU Hui. An experimental study on a suction flow control method to reduce the unsteadiness of the wind loads acting on a circular cylinder[J]. Experiments in Fluids, 2014, 55(4): 14. doi: 10.1007/s00348-014-1707-7.
|
LYSENKO D A, ERTESVG I S, and RIAN K E. Large- eddy simulation of the flow over a circular cylinder at Reynolds number 2104[J]. Flow, Turbulence and Combustion, 2014, 92(3): 673698. doi: 10.1007/s10494-013- 9509-1.
|
GIACOMO P. Equation for the determination of the density of moist air (1981)[J]. Metrologia, 1982, 18(3): 3340. doi: 10.1088/0026-1394/18/3/011.
|
LIU Cheng, DU Lidong, and ZHAO Zhan. A directional cylindrical anemometer with four sets of differential pressure sensors[J]. Review of Scientific Instruments, 2016, 87(3): 035105. doi: 10.1063/1.4943222.
|
1. | 冯存前,李江,黄大荣,胡晓伟,韩立珣. 弹道中段不同平动多目标的平动参数估计方法. 电子与信息学报. 2021(03): 564-571 . ![]() | |
2. | 王红敏,宁生科,董涛. 基于红外触发技术的外弹道目标微动参数估计研究. 兵器装备工程学报. 2021(07): 92-97 . ![]() |
参数 | 数值 |
锥体高度H (m) | 0.96 |
质心到锥顶距离d (m) | 0.64 |
底面半径r (m) | 0.25 |
自旋频率fs (Hz) | 2 |
锥旋频率fz (Hz) | 3 |
进动角θ (°) | 10 |
参数 | 数值 |
载频f0 (GHz) | 10 |
带宽B (GHz) | 1 |
脉冲宽度tp (μs) | 10 |
脉冲重复周期prf (Hz) | 100 |
驻留时间T (s) | 1 |
基线长L (m) | 200 |
锥顶散射点x | 锥顶散射点y | 锥顶散射点z | 锥底散射点x | 锥底散射点y | 锥底散射点z | |
文献[12]干涉测量 | 68.26 | 60.42 | 88.72 | 62.93 | 64.91 | 82.73 |
斜视校正干涉测量 | 88.56 | 90.41 | 96.59 | 90.84 | 91.24 | 87.28 |
进动频率{f_c}\left( {\rm Hz} \right) | 进动角\theta \left( {^ \circ } \right) | 半锥角\eta \left( {^ \circ } \right) | 锥体高度H\left( {\rm m} \right) | 底面半径r\left( {\rm m} \right) | |
真实值 | 3 | 10.00 | 14.59 | 0.96 | 0.25 |
平滑处理估计值 | 3 | 9.60 | 15.25 | 0.95 | 0.26 |
未平滑处理估计值 | 3 | 12.49 | 18.17 | 0.87 | 0.31 |
参数 | 数值 |
锥体高度H (m) | 0.96 |
质心到锥顶距离d (m) | 0.64 |
底面半径r (m) | 0.25 |
自旋频率fs (Hz) | 2 |
锥旋频率fz (Hz) | 3 |
进动角θ (°) | 10 |
参数 | 数值 |
载频f0 (GHz) | 10 |
带宽B (GHz) | 1 |
脉冲宽度tp (μs) | 10 |
脉冲重复周期prf (Hz) | 100 |
驻留时间T (s) | 1 |
基线长L (m) | 200 |
锥顶散射点x | 锥顶散射点y | 锥顶散射点z | 锥底散射点x | 锥底散射点y | 锥底散射点z | |
文献[12]干涉测量 | 68.26 | 60.42 | 88.72 | 62.93 | 64.91 | 82.73 |
斜视校正干涉测量 | 88.56 | 90.41 | 96.59 | 90.84 | 91.24 | 87.28 |
进动频率{f_c}\left( {\rm Hz} \right) | 进动角\theta \left( {^ \circ } \right) | 半锥角\eta \left( {^ \circ } \right) | 锥体高度H\left( {\rm m} \right) | 底面半径r\left( {\rm m} \right) | |
真实值 | 3 | 10.00 | 14.59 | 0.96 | 0.25 |
平滑处理估计值 | 3 | 9.60 | 15.25 | 0.95 | 0.26 |
未平滑处理估计值 | 3 | 12.49 | 18.17 | 0.87 | 0.31 |