CN104280737A - 加权宽带时反算子分解声成像方法 - Google Patents
加权宽带时反算子分解声成像方法 Download PDFInfo
- Publication number
- CN104280737A CN104280737A CN201410436820.5A CN201410436820A CN104280737A CN 104280737 A CN104280737 A CN 104280737A CN 201410436820 A CN201410436820 A CN 201410436820A CN 104280737 A CN104280737 A CN 104280737A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- omega
- max
- function
- target
- 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.)
- Granted
Links
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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
一种加权宽带时反算子分解声成像方法:将收发合置阵接收到的信号进行快速傅里叶变换,经过窄带滤波器划分为各个子带,将各个子带对应的时反算子奇异值分解,将最大特征值对应的特征向量与基于声传播模型的传递向量相关后获得模糊度函数,提取各子带模糊度函数的最大值及所对应的二维声场空间,同时以最大特征值设计加权系数,加权模糊度函数,并将所有子带的结果相干累加后呈现为三维图像。本发明的优点在于将目标的散射特性通过加权的方式应用于宽带时反算子分解成像方法,充分利用了宽带信息,真实地显示扩展性目标的空间特征,并通过实验验证了该发明方案的有效性和可靠性。
Description
技术领域
本发明属于海洋资源开发和利用的声探测领域,具体的说是一种采用阵列对点目标和扩展性目标进行声成像的方法。
背景技术
由于声波在水中传播的性能优于光波和电磁波,所以,应用于海洋资源的探测基本上是以声学手段为主。目前实际应用于声纳***的空间处理技术仍然是以上世纪五、六十年代兴起并发展起来的以白噪声背景和自由场传播条件下的平面波波束形成。然而,在浅海环境下,传统的基于平面波模型的处理方法不再适用,由于声信号在海洋中的传播呈现多路径/多简正模状态,很容易造成检测过程中的伪目标或虚影现象,影响声纳***的检测性能。时反处理从上世纪末开始应用于水声领域,是目前颇有应用前景的主动目标探测手段,时反算子分解是通过分析目标散射矩阵达到聚焦发射、定位以及成像目标的一种方法。
然而,目前时反算子分解方法是每个频点对应一个目标散射矩阵,每个频带单独处理的方式,存在如下缺点:(1)宽带信号的处理通常采用仅仅考虑中心频率或者整个频段的处理结果非相干累加的方式,没有有效地利用宽带信号的信息;(2)特别是波导环境中的扩展性目标,由于受目标本身的散射特性、波导传播效应以及换能器与频率的关系等因素的影响,最终表现出的目标散射特性是频率的函数,而目前处理过程中并没有考虑和利用这种散射特性。
发明内容
本发明要克服现有技术的上述缺点,提供一种加权宽带时反算子分解声成像方法,在现有的时反算子分解方法的基础上,有效地利用目标散射特性和宽带信息,从而达到消除伪目标以及有效显示扩展性目标空间特征的目的。
本发明提出的加权宽带时反算子分解声成像方法,提取每个子带模糊度函数的最大值以及相应的二维声场空间,根据最大特征值与频率的函数关系设计加权系数,加权子带的模糊度函数,并将加权后的模糊度函数相干累加,将最终结果显示为三维声场图像,确定目标的二维声场空间。
下面对本发明方法做进一步说明。
对环境的声学参数如水深、声速剖面、密度、沉积层参数等有足够了解的情况下,按如下步骤完成目标的三维声学图像及目标空间位置的识别:
(1)建立坐标系。设P个换能器组成的收发合置线阵,垂直放置在水中,以线阵为z轴,水平方向为r轴建立坐标系,水面为坐标原点。第一个换能器距水面的深度为z1,第P个换能器距水面的深度为zP;
(2)将换能器阵列接收到数据快拍排列成列向量,yi(t)=[y1i(t)…ypi(t)…yPi(t)]T,p=1,…,P,i=1,…,T,T表示数据快拍数;
(3)根据预定采样频率对每个数据快拍进行采样,
yi(n)=[y1i(n)…ypi(n)…yPi(n)]T;
(4)对采样后的信号进行短时傅里叶变换
其中w(n-m)为窗口函数序列;
(5)将所有的T个数据快拍排列成矩阵
则时反算子为K(ω)=Y(ω)YH(ω);
(6)将时反算子划分为多个子带,进行奇异值分解:
K(ω)=U(ω)Λ(ω)V*(ω) (3)
(7)输出各子带非相干累加后的特征值与水平距离r的函数曲线,根据最大特征值与水平距离的函数曲线判断包含目标的时间窗,找出所对应的时反算子Ko(ω),特征值Λo(ω)以及特征向量Uo(ω)和Vo*(ω);
(8)根据特征值Λo(ω)设计加权系数,比较所有子带的最大特征值λ1(ω),找出其最大峰值λmax,其余各子带的最大特征值
与峰值λmax的比值即为加权系数:
η(ω)=λ1(ω)/λmax (4)
(9)将感兴趣的空间区域划分为网格,网格的交点即为假定目标所在的位置(r,z),其中r表示假定目标距收发合置阵的水平距离,z表示假定目标的水深;
(10)根据成像方法所面临的水声环境,确定所适用的传播模型,分别得到假定目标与收发合置阵之间的传递向量
g(r,z,ω)=[g1(z1,r,z,ω)…gp(zp,r,z,ω)…gP(zP,r,z,ω)]T,
其中g1(z1,r,z,ω)表示假定目标与第一个换能器之间的传递函数,例如,简正模传播模型的传递函数为:
其中Zl(·)表示第l号传播模对应的特征函数,H0 (2)(·)为汉克尔函数,κl表示第l号传播模对应的波数;
(11)各个子带的模糊度函数为:
I(r,z,ω)=|gH(r,z,ω)u1(ω)|2 (6)
其中u1(ω)表示Λo(ω)的最大特征值λ1(ω)对应的特征向量;
(12)提取各个子带模糊度函数的最大值Imax(r,z,ω)及其对应的声场空间位置(rmax(ω),zmax(ω));
(13)将所有子带的声场空间(rmax(ω),zmax(ω))整理为一个集合Ω(rmax,zmax),并重新定义模糊度函数:
(14)对所有子带的模糊度函数进行加权累加:
Iend(r,z,ω)=Ση(ω)I(r,z,ω) (8)
将最终获得的模糊度函数显示为距离r和深度z的三维图像,可确定目标所对应的二维空间区域。
本发明的优点在于将目标的散射特性通过加权的方式应用于宽带时反算子分解成像方法,真实地显示扩展性目标的空间特征,充分利用了宽带信息,并通过实验验证了该发明方案的有效性和可靠性。
附图说明
图1本发明的实验布置及坐标设置图。
图2本发明的特征值与距离的函数关系图。
图3本发明根据最大特征值设计的加权系数曲线。
图4为本发明目标在二维波导空间的定位图。
图5为本发明目标空间特征放大展示图。
具体实施方案
下面通过具体实施例子对本发明作进一步地描述。在一个深度为1.44m的波导实验室水池中,水池的三面贴消声瓦,水池的底部铺了一层0.22m厚的沙。水池的环境参数如下:水中声速为常数,通过测量水体温度计算得知为c1=1493m/s,水体密度ρ1=1000kg/m3,沉积层参数密度ρ2=1800kg/m3,声速为c2=1650m/s,衰减系数α2=0.67dB/λ,基底的密度为ρ3=1800kg/m3,声速c3=1580m/s,衰减系数α3=0.8dB/λ。
(1)本实施例子中,实验布置及坐标系的建立如图1所示,收发合置阵32个阵元,全场布阵,第一个阵元距水面0.04m,阵元间距0.04m,垂直放置在水池中。一个直径为0.21m,长度为0.51m的圆柱型目标放置在距阵约8m的水池底部。发射信号为脉冲宽度为1ms的中心频率为12kHz,带宽为4kHz的线性调频信号。
(2)将换能器阵列接收到数据快拍排列成列向量,
yi(t)=[y1i(t)…ypi(t)…yPi(t)]T,p=1,…,32,i=1,…,8;
(3)根据采样频率为50kHz对每个数据快拍进行采样,
yi(n)=[y1i(n)…ypi(n)…yPi(n)]T;
(4)选择的短时窗为发射脉宽的两倍,即点数为100的矩形窗,对采样后的信号进行短时傅里叶变换
(5)将所有的8个数据快拍排列成矩阵
则时反算子为K(ω)=Y(ω)YH(ω);
(6)将4kHz的带宽的时反算子划分为401个子带进行奇异值分解:
K(ω)=U(ω)Λ(ω)V*(ω) (3)
(7)输出各子带非相干累加后的特征值与水平距离r的函数曲线,根据最大特征值与水平距离的函数曲线判断包含目标的时间窗,找出所对应的时反算子Ko(ω),特征值Λo(ω)以及特征向量Uo(ω)和Vo*(ω);图2是前5个特征值与距离的函数关系图,以所有子带最大特征值的最大值规一化后以dB的形式显示,由上往下依次为最大、第二、第三、第四和第五特征值的函数关系曲线。由图2可知,最大特征值在目标所在的距离(8m)附近有一个局部峰值,从而可以确定目标所在的时间窗。
(8)根据特征值Λo(ω)设计加权系数,比较所有子带的最大特征值λ1(ω),找出其最大峰值λmax,由图3所见,在中心频率12kHz附近,纵坐标显示为0dB,其余各子带的最大特征值与峰值λmax的比值即为加权系数:
η(ω)=λ1(ω)/λmax (4)
图3是根据最大特征值设计的加权系数,结果以dB形式显示。
(9)将感兴趣的空间区域划分为网格,水平起始距离为0.1m,搜索步长为0.05m,截止距离为15m,垂直起始深度为0m,搜索步长为0.0075m,截止深度为1.5m。网格的交点即为假定目标所在的位置(r,z),其中r表示假定目标距收发合置阵的水平距离,z表示假定目标的水深;
(15)本例采用简正模传播模型,传递向量为
g(r,z,ω)=[g1(z1,r,z,ω)…gp(zp,r,z,ω)…gP(zP,r,z,ω)]T,
其中g1(z1,r,z,ω)表示假定目标与第一个换能器之间的传递函数,表达式为:
其中Zl表示第l号传播模对应的特征函数,H0 (2)(·)为汉克尔函数,κl表示波数;
(10)各个子带的模糊度函数为:
I(r,z,ω)=|gH(r,z,ω)u1(ω)|2 (6)
其中u1(ω)表示Λo(ω)的最大特征值λ1(ω)对应的特征向量;
(11)提取各个子带模糊度函数的最大值Imax(r,z,ω)及其对应的声场空间位置(rmax(ω),zmax(ω));
(12)将所有子带的声场空间(rmax(ω),zmax(ω))整理为一个集合Ω(rmax,zmax),并重新定义模糊度函数:
(13)对所有子带的模糊度函数进行加权累加:
Iend(r,z,ω)=Ση(ω)I(r,z,ω) (8)
图4是目标在二维波导空间的定位图,即最终获得的模糊度函数与距离r和深度z的三维图像,展示了目标在整个水池所处的空间位置,图5目标空间特征放大展示图。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (1)
1.加权宽带时反算子分解声成像方法,包括如下步骤:
(1)建立坐标系。设P个换能器组成的收发合置线阵,垂直放置在水中,以线阵为z轴,水平方向为r轴建立坐标系,水面为坐标原点。第一个换能器距水面的深度为z1,第P个换能器距水面的深度为zP;
(2)将换能器阵列接收到数据快拍排列成列向量,
yi(t)=[y1i(t)…ypi(t)…yPi(t)]T,p=1,…,P,i=1,…,T,T表示数据快拍数;
(3)根据预定采样频率对每个数据快拍进行采样,
yi(n)=[y1i(n)…ypi(n)…yPi(n)]T;
(4)对采样后的信号进行短时傅里叶变换
其中w(n-m)为窗口函数序列;
(5)将所有的T个数据快拍排列成矩阵
则时反算子为K(ω)=Y(ω)YH(ω);
(6)将时反算子划分为多个子带,进行奇异值分解:
K(ω)=U(ω)Λ(ω)V*(ω) (3)
(7)输出各子带非相干累加后的特征值与水平距离r的函数曲线,根据最大特征值与水平距离的函数曲线判断包含目标的时间窗,找出所对应的时反算子Ko(ω),特征值Λo(ω)以及特征向量Uo(ω)和Vo*(ω);
(8)根据特征值Λo(ω)设计加权系数,比较所有子带的最大特征值λ1(ω),找出其最大峰值λmax,其余各子带的最大特征值与峰值λmax的比值即为加权系数:
η(ω)=λ1(ω)/λmax (4)
(9)将感兴趣的空间区域划分为网格,网格的交点即为假定目标所在的位置(r,z),其中r表示假定目标距收发合置阵的水平距离,z表示假定目标的水深;
(10)根据成像方法所面临的水声环境,确定所适用的传播模型,分别得到假定目标与收发合置阵之间的传递向量
g(r,z,ω)=[g1(z1,r,z,ω)…gp(zp,r,z,ω)…gP(zP,r,z,ω)]T,
其中g1(z1,r,z,ω)表示假定目标与第一个换能器之间的传递函数,例如,简正模传播模型的传递函数为:
其中Zl(·)表示第l号传播模对应的特征函数,H0 (2)(·)为汉克尔函数,κl表示第l号传播模对应的波数;
(11)各个子带的模糊度函数为:
I(r,z,ω)=|gH(r,z,ω)u1(ω)|2 (6)
其中u1(ω)表示Λo(ω)的最大特征值λ1(ω)对应的特征向量;
(12)提取各个子带模糊度函数的最大值Imax(r,z,ω)及其对应的声场空间位置(rmax(ω),zmax(ω));
(13)将所有子带的声场空间(rmax(ω),zmax(ω))整理为一个集合Ω(rmax,zmax),并重新定义模糊度函数:
(14)对所有子带的模糊度函数进行加权累加:
Iend(r,z,ω)=Ση(ω)I(r,z,ω) (8)
将最终获得的模糊度函数显示为距离r和深度z的三维图像,可确定目标所对应的二维空间区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410436820.5A CN104280737B (zh) | 2014-08-29 | 2014-08-29 | 加权宽带时反算子分解声成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410436820.5A CN104280737B (zh) | 2014-08-29 | 2014-08-29 | 加权宽带时反算子分解声成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104280737A true CN104280737A (zh) | 2015-01-14 |
CN104280737B CN104280737B (zh) | 2017-01-11 |
Family
ID=52255801
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410436820.5A Active CN104280737B (zh) | 2014-08-29 | 2014-08-29 | 加权宽带时反算子分解声成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104280737B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106019290A (zh) * | 2016-05-26 | 2016-10-12 | 浙江工业大学 | 加权宽带时反算子分解多目标声成像方法 |
CN106330794A (zh) * | 2015-06-29 | 2017-01-11 | 深圳市中兴微电子技术有限公司 | 一种均衡器系数的调整方法及装置 |
CN110456342A (zh) * | 2019-07-11 | 2019-11-15 | 西安电子科技大学 | 单发射天线雷达的远场多运动目标探测方法 |
CN112684437A (zh) * | 2021-01-12 | 2021-04-20 | 浙江大学 | 一种基于时域warping变换的被动测距方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040174770A1 (en) * | 2002-11-27 | 2004-09-09 | Rees Frank L. | Gauss-Rees parametric ultrawideband system |
CN102279394A (zh) * | 2011-05-17 | 2011-12-14 | 西安电子科技大学 | 低频窄带高分辨超声探测成像方法 |
CN102928844A (zh) * | 2012-11-08 | 2013-02-13 | 中北大学 | 一种水下亚波长分辨率三维成像方法 |
CN103513246A (zh) * | 2013-10-09 | 2014-01-15 | 中国科学院空间科学与应用研究中心 | 一种亚波长成像的***和方法 |
-
2014
- 2014-08-29 CN CN201410436820.5A patent/CN104280737B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040174770A1 (en) * | 2002-11-27 | 2004-09-09 | Rees Frank L. | Gauss-Rees parametric ultrawideband system |
CN102279394A (zh) * | 2011-05-17 | 2011-12-14 | 西安电子科技大学 | 低频窄带高分辨超声探测成像方法 |
CN102928844A (zh) * | 2012-11-08 | 2013-02-13 | 中北大学 | 一种水下亚波长分辨率三维成像方法 |
CN103513246A (zh) * | 2013-10-09 | 2014-01-15 | 中国科学院空间科学与应用研究中心 | 一种亚波长成像的***和方法 |
Non-Patent Citations (3)
Title |
---|
刘敏; 曹燕; 韦岗: "基于时反处理法的水下目标探测研究", 《科学技术与工程》, vol. 9, no. 21, 30 November 2009 (2009-11-30) * |
宫蕴瑞; 何迪; 何晨; 邓浩: "时间反转超宽带***的空间谱定位算法", 《上海交通大学学报》, vol. 46, no. 6, 30 June 2012 (2012-06-30) * |
李春晓; 宫先仪; 邹丽娜: "干扰环境下窄带信号时反算子分解聚焦性能分析", 《哈尔滨工程大学学报》, vol. 29, no. 8, 31 August 2008 (2008-08-31) * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106330794A (zh) * | 2015-06-29 | 2017-01-11 | 深圳市中兴微电子技术有限公司 | 一种均衡器系数的调整方法及装置 |
CN106330794B (zh) * | 2015-06-29 | 2019-10-22 | 深圳市中兴微电子技术有限公司 | 一种均衡器系数的调整方法及装置 |
US10461968B2 (en) | 2015-06-29 | 2019-10-29 | Sanechips Technology Co., Ltd. | Adjustment method and device for equalizer coefficient, and computer storage medium |
CN106019290A (zh) * | 2016-05-26 | 2016-10-12 | 浙江工业大学 | 加权宽带时反算子分解多目标声成像方法 |
CN106019290B (zh) * | 2016-05-26 | 2018-11-13 | 浙江工业大学 | 加权宽带时反算子分解多目标声成像方法 |
CN110456342A (zh) * | 2019-07-11 | 2019-11-15 | 西安电子科技大学 | 单发射天线雷达的远场多运动目标探测方法 |
CN110456342B (zh) * | 2019-07-11 | 2023-03-21 | 西安电子科技大学 | 单发射天线雷达的远场多运动目标探测方法 |
CN112684437A (zh) * | 2021-01-12 | 2021-04-20 | 浙江大学 | 一种基于时域warping变换的被动测距方法 |
CN112684437B (zh) * | 2021-01-12 | 2023-08-11 | 浙江大学 | 一种基于时域warping变换的被动测距方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104280737B (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gervaise et al. | Passive geoacoustic inversion with a single hydrophone using broadband ship noise | |
CN102183435B (zh) | 一种基于多路径反射理论的海底密度和声速测量方法 | |
CN104749568B (zh) | 一种基于水听器阵列的浅海目标深度的分类方法 | |
Battle et al. | Geoacoustic inversion of tow-ship noise via near-field-matched-field processing | |
CN101387701A (zh) | 一种基于前向预测的被动时间反转混响抑制方法 | |
CN105301114A (zh) | 一种基于多通道空时逆滤波技术的声学覆盖层***损失测量方法 | |
CN104280737A (zh) | 加权宽带时反算子分解声成像方法 | |
CN103076594A (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
Jiang et al. | Estimation of geoacoustic properties of marine sediment using a hybrid differential evolution inversion method | |
CN106019290B (zh) | 加权宽带时反算子分解多目标声成像方法 | |
CN101777954A (zh) | 基于本征声线视在搜索的水声信道有效声速估计方法 | |
CN115166817B (zh) | 一种基于冰层模态群慢度差特征的冰声定位方法 | |
Kargl et al. | Synthetic aperture sonar imaging of simple finite targets | |
CN104990994A (zh) | 基于脉冲压缩技术的沉积物中声速和衰减的宽带测量方法 | |
Kuperman et al. | Ocean acoustics, matched-field processing and phase conjugation | |
Gervaise et al. | Passive acoustic tomography: new concepts and applications using marine mammals: a review | |
CN115236592B (zh) | 一种基于单阵元时频曲线匹配的冰声定位方法 | |
Ballard et al. | A comparison between directly measured and inferred wave speeds from an acoustic propagation experiment in Currituck Sound | |
CN113126029A (zh) | 适用于深海可靠声路径环境的多传感器脉冲声源定位方法 | |
Culver et al. | Sonar signal processing using probabilistic signal and ocean environmental models | |
Reeder et al. | Walking on snow-covered Arctic sea ice to infer ice thickness | |
Waters et al. | Isolating scattering resonances of an air-filled spherical shell using iterative, single-channel time reversal | |
Chotiros et al. | Refraction and scattering into a sandy ocean sediment in the 30-40-kHz band | |
Li et al. | Feature extraction for acoustic scattering from a buried target | |
Shi et al. | Data‐based differential matched‐mode source depth estimation for a moving source |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |