CN114397008B - 低空飞行器噪声源识别的解多普勒方法 - Google Patents
低空飞行器噪声源识别的解多普勒方法 Download PDFInfo
- Publication number
- CN114397008B CN114397008B CN202111487257.0A CN202111487257A CN114397008B CN 114397008 B CN114397008 B CN 114397008B CN 202111487257 A CN202111487257 A CN 202111487257A CN 114397008 B CN114397008 B CN 114397008B
- Authority
- CN
- China
- Prior art keywords
- grid
- doppler
- point
- low
- noise source
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种低空飞行器噪声源识别的解多普勒方法,其包括:利用传感器进行扫描,得到对低空飞行器的扫描网格;根据所要计算的区域边界点,确定最大允许的扫描网格边长;在确定的扫描网格边长范围内,计算中心网格点的解多普勒数据。本发明具有原理简单、操作简便、可降低运算量、提高计算效率等优点。
Description
技术领域
本发明主要涉及到飞行器噪声识别技术领域,特指一种低空飞行器噪声源识别的解多普勒方法。
背景技术
目前,低空飞行器的噪声源识别方法多采用波束形成方法,它利用传声器阵列采集运动/静止目标的辐射噪声,对目标所在区域划分网格后采用球面波假设对各网格点进行数据处理,进而获得目标各部位的噪声源强度与频谱特征。
传统技术中,针对低空飞行器的噪声源识别处理流程一般如下:
(1)利用传声器阵列采集飞行器飞行噪声,并基于航迹数据对声信号同步和解多普勒;
(2)基于阵列在不同频率下固有的点扩展函数(PSF, point spread function),采用不同的波束形成类方法对飞行噪声进行声成像;
(3)在波束形成方法基础上采用解卷积类方法或Clean类方法等高分辨方法,对波束形成结果优化,得到高分辨成像结果;
(4)结合机体结构和飞行状态,分析不同倍频程不同机体部位的噪声分布特性,给出当前飞行状态下的飞行器噪声辐射特征。
在现有技术中,低空飞行器与行驶的汽车、高铁等目标类似,其运动造成的多普勒效应都对其噪声源识别造成了一定影响。
针对此问题,传统手段大致分为两类处理方法:
(1)忽略多普勒效应;
此类方法通过提高采样率,截取短时段内的阵列接收信号,假定此时运动目标为静止目标,忽略多普勒效应直接处理。该方法要求在截取时段内,目标运动距离较小,同时为保证在该时段内有足够的采样数据进行处理,需要保持较高的采样率。对于低空飞行器而言,由于其飞行速度一般在100km/h以上,且机体尺寸长宽方向均较大,较难满足上述条件。此时忽略多普勒效应的计算结果仅能定性给出噪声源识别结果,即强噪声源的大致区域。
(2)逐点解多普勒;
低空飞行器飞行速度快,尺寸大,造成了不同部位的多普勒效应差异明显。为此通常对飞行器所在平面划分细密网格,假定每个网格点上存在一个辐射强度未知的声源,通常地面阵列尺寸也较大,可达直径百米左右,造成每个网格点与地面阵列中的每个传声器存在一个多普勒效应。按照亚音速莫尔斯声学理论,为了获得解多普勒后的阵列接收信号,需要声源面上每个网格点对阵列各个传声器进行解多普勒,由此造成解多普勒运算量巨大,尤其是不同网格划分方式下都需要重新解多普勒。
发明内容
本发明要解决的技术问题就在于:针对现有技术存在的技术问题,本发明提供一种原理简单、操作简便、可降低运算量、提高计算效率的低空飞行器噪声源识别的解多普勒方法。
为解决上述技术问题,本发明采用以下技术方案:
一种低空飞行器噪声源识别的解多普勒方法,其包括:
利用传感器进行扫描,得到对低空飞行器的扫描网格;
根据所要计算的区域边界点,确定最大允许的扫描网格边长;
在确定的扫描网格边长范围内,计算中心网格点的解多普勒数据。
作为本发明的进一步改进:当需要在更精细网格下进行解多普勒和波束形成时,就近选取“粗网格点”的解多普勒结果,所述“粗网格点”为确定的扫描网格边长范围内的网格点。
作为本发明的进一步改进:所述确定最大允许的扫描网格边长是通过判定相邻网格点造成的幅度误差是否可以忽略来确定。
作为本发明的进一步改进:相邻网格点进行替换,在替换之后,进行波束形成导向矢量计算时,相位误差方面达到设定阈值,视为确定最大允许的扫描网格的临界节点。
作为本发明的进一步改进:对于网格点A和B,假定它们与地面阵列上第i个阵元的距离为和,那么用A点代替B点进行波束形成导向矢量计算时,引入的幅度和相位误差为:
幅度误差: i=1,…M
相位误差: i=1,…M
此时若AB距离为设定阈值(如1m),造成的幅度误差可以忽略;相位误差方面满足:
i=1,…M。
作为本发明的进一步改进:若是单网格点对单传声器,那么:
对于某个网格点所截取的时段内的辐射信号的时间轴是等间隔时间序列;根据航迹和阵列位置关系,得到其对应的时段内的非等间隔接收时间序列;通过降采样插值并补偿传播损失后,得到解多普勒后信号。
作为本发明的进一步改进:若是单网格点对整个传声器阵列,那么:
对于单个网格点其截取的辐射时段是确定的,传声器阵列中各个传声器接收信号时间序列,其中
j=1,…,
M,
M为阵元数,起止点各不相同。
传声器阵列中各通道应截取相同接收时刻开始的接收信号,即:
同时,可截取的最长时长为:
。
作为本发明的进一步改进:若是多网格点对整个传声器阵列,那么:
在飞行器上划分了多个网格点;假定每个网格点下传感器阵列解多普勒后信号的时间长度为,
i=1,…,
N,
N为网格点数,它们存在一定差异,对各网格点采取相同的时间长度。
作为本发明的进一步改进:根据截取的辐射信号时段[
tea,
teb],分别计算各网格点的时长,取其最小值作为最终采用的辐射时间长度,即:
接下来,进行解多普勒。
与现有技术相比,本发明的优点就在于:
1、本发明的低空飞行器噪声源识别的解多普勒方法,原理简单、操作简便,通过较为稀疏的网格划分来降低解多普勒运算量和数据缓存量,降低运算量的同时可以为后续计算缓存可用的解多普勒结果。
2、本发明的低空飞行器噪声源识别的解多普勒方法,简化了解多普勒的声源面网格区域划分,在不影响波束形成性能的基础上,大幅降低解多普勒的计算量;同时,本发明支持同一种网格划分下解多普勒结果缓存,并应用于其他各种网格划分模式下的噪声源识别。
3、本发明的低空飞行器噪声源识别的解多普勒方法,将相近区域网格点统一采用区域中心点的解多普勒结果,由此造成的幅度补偿误差对波束形成的影响可以忽略,造成的相频误差对波束形成的主瓣影响可以忽略,仅在旁瓣水平略有差异。因此,本发明在不影响后续波束形成精度的前提下,大幅减小解多普勒运算量,同时解多普勒结果可以缓存用于不同网格划分模式下的波束形成运算,加快外场飞行试验数据处理。
附图说明
图1是本发明方法的流程示意图。
图2是本发明在具体应用实例中网格点选取造成的误差示意图。
图3是本发明在具体应用实例中声成像粗网格划分示意图。
具体实施方式
以下将结合说明书附图和具体实施例对本发明做进一步详细说明。
如图1所示,本发明的低空飞行器噪声源识别的解多普勒方法,包括:
利用传感器进行扫描,得到对低空飞行器的扫描网格;
根据所要计算的区域边界点,确定最大允许的扫描网格边长;
在确定的扫描网格边长范围内,计算中心网格点的解多普勒数据。
在具体应用实例中,进一步,当需要在更精细网格下进行解多普勒和波束形成时,可以就近选取“粗网格点”的解多普勒结果,所述“粗网格点”为确定的扫描网格边长范围内的网格点,从而减少解多普勒运算量,并对整体噪声源识别性能影响不大。
通过采用本发明的上述方案,利用确定的扫描网格边长范围内的网格,即较为稀疏的网格划分来降低解多普勒运算量和数据缓存量,在降低了运算量的同时,还可以为后续计算缓存可用的解多普勒结果。
由上可知,本发明在逐点解多普勒方法基础上,简化了解多普勒的声源面网格区域划分,在不影响波束形成性能的基础上,大幅降低解多普勒的计算量;例如,仍对5m×5m的部件区域解多普勒,按照1m网格边长划分有36个网格点,阵列传声器数目若仍为100个,则只需进行3600次解多普勒运算,该结果仍可用于网格边长0.2m的噪声源识别,即计算量降为原来的5%。
本发明在逐点解多普勒方法基础上,能够支持同一种网格划分下解多普勒结果缓存,并应用于其他各种网格划分模式下的噪声源识别。例如,仍对5m×5m的部件区域解多普勒,按照1m网格边长划分,该解多普勒结果可用于该部件区域任意网格边长或非网格划分方式下的噪声源识别。此种情况下,噪声源识别时新的网格点取距离其最近的缓存网格点解多普勒结果即可,由此造成的多普勒和幅度衰减效应对波束形成结果基本无影响。
在具体应用实例中,所述确定最大允许的扫描网格边长是通过判定相邻网格点造成的幅度误差是否可以忽略来确定。也就是说,如果相邻网格点在替换之后,进行波束形成导向矢量计算时,相位误差方面达到设定阈值,即可视为确定最大允许的扫描网格的临界节点。
如图2所示,对于网格点A和B,假定它们与地面阵列上第i个阵元的距离为和,那么用A点代替B点进行波束形成导向矢量计算时,引入的幅度和相位误差为:
幅度误差: i=1,…M
相位误差: i=1,…M
通常低空飞行器的过顶高度在20-100m,此时若AB距离为1m,造成的幅度误差可以忽略。相位误差方面只要满足:
i=1,…M
此时,相位误差对常规波束形成空间谱的主瓣基本没有影响,旁瓣也仅略有波动,对地面阵列的整体噪声源识别性能影响不大。
在具体应用实例中,如图3所示,根据所要计算的区域边界点,可以按照上式确定最大允许的网格边长,仅计算这些中心网格点的解多普勒数据。
在具体应用实例中,若是单网格点对单传声器,那么:
对于某个网格点所截取的时段内的辐射信号的时间轴是等间隔时间序列。根据航迹和阵列位置关系,可得到其对应的时段内的非等间隔接收时间序列。
那么,通过降采样插值并补偿传播损失后,即可得到解多普勒后信号。
在具体应用实例中,若是单网格点对整个传声器阵列,那么:
对于单个网格点其截取的辐射时段是确定的,但是传声器阵列中各个传声器接收信号时间序列,其中
j=1,…,
M,
M为阵元数,起止点各不相同。
为便于后续噪声源识别处理,传声器阵列中各通道应截取相同接收时刻开始的接收信号,即:
同时,可截取的最长时长为:
。
在具体应用实例中,若是多网格点对整个传声器阵列,那么:
在飞行器上划分了多个网格点;
假定每个网格点下传感器阵列解多普勒后信号的时间长度为,
i=1,…,
N,
N为网格点数,它们存在一定差异。这不利于比较各网格点下的声源辐射能量,故要对各网格点采取相同的时间长度。
在具体应用时,实际处理中可根据截取的辐射信号时段,分别计算各网格点的时长,取其最小值作为最终采用的辐射时间长度,即:
。
接下来,解多普勒步骤与单网格点对整个传声器阵列的情形相同,在此就不在赘述。
按照上述步骤,解多普勒后的传感器阵列接收信号和静止声源的接收信号相似,接下来可利用其计算采样协方差矩阵,以估计噪声源的参数。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (6)
1.一种低空飞行器噪声源识别的解多普勒方法,其特征在于,包括:
利用传感器进行扫描,得到对低空飞行器的扫描网格;
根据所要计算的区域边界点,确定最大允许的扫描网格边长;
在确定的扫描网格边长范围内,计算中心网格点的解多普勒数据;
所述确定最大允许的扫描网格边长是通过判定相邻网格点造成的幅度误差是否可以忽略来确定;相邻网格点在替换之后,进行波束形成导向矢量计算时,相位误差方面达到设定阈值,即可视为确定最大允许的扫描网格的临界节点;对于网格点A和B,假定它们与地面阵列上第i个阵元的距离为和,那么用A点代替B点进行波束形成导向矢量计算时,引入的幅度和相位误差为:
幅度误差: i=1,…M
相位误差: i=1,…M
此时若AB距离为在预设阈值内,造成的幅度误差忽略;相位误差方面满足:
i=1,…M。
2.根据权利要求1所述的低空飞行器噪声源识别的解多普勒方法,其特征在于,当需要在更精细网格下进行解多普勒和波束形成时,就近选取“粗网格点”的解多普勒结果,所述“粗网格点”为确定的扫描网格边长范围内的网格点。
3.根据权利要求1或2所述的低空飞行器噪声源识别的解多普勒方法,其特征在于,若是单网格点对单传声器,那么:
对于某个网格点所截取的时段内的辐射信号的时间轴是等间隔时间序列;根据航迹和阵列位置关系,得到其对应的时段内的非等间隔接收时间序列;通过降采样插值并补偿传播损失后,得到解多普勒后信号。
4.根据权利要求1或2所述的低空飞行器噪声源识别的解多普勒方法,其特征在于,若是单网格点对整个传声器阵列,那么:
对于单个网格点其截取的辐射时段是确定的,传声器阵列中各个传声器接收信号时间序列,其中j=1,…,M,M为阵元数,起止点各不相同;
传声器阵列中各通道应截取相同接收时刻开始的接收信号,即:
同时,截取的最长时长为:
。
5.根据权利要求1或2所述的低空飞行器噪声源识别的解多普勒方法,其特征在于,若是多网格点对整个传声器阵列,那么:
在飞行器上划分了多个网格点;假定每个网格点下传感器阵列解多普勒后信号的时间长度为,i=1,…, N,N为网格点数,它们存在一定差异,对各网格点采取相同的时间长度。
6.根据权利要求5所述的低空飞行器噪声源识别的解多普勒方法,其特征在于,根据截取的辐射信号时段,分别计算各网格点的时长,取其最小值作为最终采用的辐射时间长度,即:
接下来,进行解多普勒。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111487257.0A CN114397008B (zh) | 2021-12-07 | 2021-12-07 | 低空飞行器噪声源识别的解多普勒方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111487257.0A CN114397008B (zh) | 2021-12-07 | 2021-12-07 | 低空飞行器噪声源识别的解多普勒方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114397008A CN114397008A (zh) | 2022-04-26 |
CN114397008B true CN114397008B (zh) | 2023-04-28 |
Family
ID=81227793
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111487257.0A Active CN114397008B (zh) | 2021-12-07 | 2021-12-07 | 低空飞行器噪声源识别的解多普勒方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114397008B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116086599B (zh) * | 2023-04-10 | 2023-10-20 | 中国飞行试验研究院 | 航空器低空辐射声场获取方法、装置、设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112180329A (zh) * | 2020-09-07 | 2021-01-05 | 黑龙江工程学院 | 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5864799B1 (ja) * | 2015-03-30 | 2016-02-17 | 株式会社日立パワーソリューションズ | 音源探査装置および音源探査方法 |
CN104777453B (zh) * | 2015-04-23 | 2017-10-10 | 西北工业大学 | 舰船线谱噪声源定位的波束域时频分析方法 |
CN106772248B (zh) * | 2016-12-09 | 2019-04-09 | 中国空气动力研究与发展中心低速空气动力研究所 | 高速列车模型在声学风洞内的噪声源定位方法 |
CN107563551A (zh) * | 2017-08-28 | 2018-01-09 | 国网湖南省电力公司 | 一种电网舞动关键线路辨别方法及其系统 |
US11099268B2 (en) * | 2018-12-19 | 2021-08-24 | GM Global Technology Operations LLC | Doppler ambiguity resolution via high order phase terms |
CN112904302B (zh) * | 2021-04-22 | 2022-06-28 | 中国人民解放军海军航空大学 | 网格化frft域雷达目标检测与多级联合虚警剔除方法 |
-
2021
- 2021-12-07 CN CN202111487257.0A patent/CN114397008B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112180329A (zh) * | 2020-09-07 | 2021-01-05 | 黑龙江工程学院 | 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114397008A (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102955150B (zh) | 基于天线主瓣占优强度约束的sar方位模糊抑制方法 | |
CN103018727A (zh) | 一种基于样本训练的机载雷达非平稳杂波抑制方法 | |
CN114397008B (zh) | 低空飞行器噪声源识别的解多普勒方法 | |
Wang et al. | An experimental study of passive bistatic radar using uncooperative radar as a transmitter | |
US20120173191A1 (en) | Airspeed And Velocity Of Air Measurement | |
CN109669182B (zh) | 无源双基地sar动/静目标联合稀疏成像方法 | |
CN108267723A (zh) | 陆基短距k波段雷达的离线在线式幅相误差校正方法 | |
CN112014806B (zh) | 一种复杂干扰场景下的机载雷达无意干扰抑制方法 | |
JP5025359B2 (ja) | レーダ装置 | |
CN102288950B (zh) | 基于stap的高速空中机动目标检测方法 | |
CN105738887A (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
KR101426226B1 (ko) | 레이더의 신호 처리 방법 | |
CN109557561B (zh) | 基于稳健协方差矩阵的ads-b压制式干扰抑制方法 | |
CN110261837B (zh) | 一种基于航迹信息的复杂目标rcs计算方法 | |
CN114509733A (zh) | 一种基于联合对消的多通道sar干扰抑制方法 | |
CN107783111B (zh) | 一种基于最大熵准则的雷达前视超分辨成像方法 | |
CN113484859A (zh) | 一种基于融合技术的二维超分辨雷达成像方法 | |
CN104569970B (zh) | 一种用于机载雷达前视单脉冲成像的自聚焦方法 | |
JP3850610B2 (ja) | ウィンドプロファイラ装置 | |
CN104898127A (zh) | 一种多发多收合成孔径激光雷达系统及其体制设计方法 | |
CN114910905A (zh) | 相似性约束下geo星机双基sar动目标智能成像方法 | |
CN114325613A (zh) | 一种区域性非均匀蒸发波导下的雷达探测威力预测方法 | |
Li et al. | DME interference suppression algorithm based on signal separation estimation theory for civil aviation system | |
Denny et al. | Anomalous propagation limitations to high-resolution SAR performance | |
Boone et al. | Microphone array technology for aircraft noise measurements |
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 |