CN104062645B - 一种测量参量阵差频波与同频小振幅波信号相位差的方法 - Google Patents
一种测量参量阵差频波与同频小振幅波信号相位差的方法 Download PDFInfo
- Publication number
- CN104062645B CN104062645B CN201410298876.9A CN201410298876A CN104062645B CN 104062645 B CN104062645 B CN 104062645B CN 201410298876 A CN201410298876 A CN 201410298876A CN 104062645 B CN104062645 B CN 104062645B
- Authority
- CN
- China
- Prior art keywords
- signal
- hydrophone
- small amplitude
- parametric array
- amplitude wave
- 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
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52004—Means for monitoring or calibrating
-
- 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)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明属于水声测量领域,具体涉及一种测量参量阵差频波与同频小振幅波信号相位差的方法。本发明包括:在水中吊放参量阵声源、小振幅波声源、第一水听器和第二水听器,触发信号发生器产生两路触发信号;参量阵声源和小振幅波声源在Δt时间间隔后分别发射信号,采集器采集并保存第一水听器和第二水听器接收到的信号;进行相关处理,找出处理结果中的相关峰;得到参量阵差频波与同频小振幅波信号相位差。本发明测量数据的处理方法简单,不需要复杂的信号处理算法。
Description
技术领域
本发明属于水声测量领域,具体涉及一种测量参量阵差频波与同频小振幅波信号相位差的方法。
背景技术
参量阵广泛应用于水声定位、水声导航、水中通讯、水下探测等许多方面。参量阵是利用了介质的非线性作用,由声源发射两个频率相近的高强声波(原波),在其互作用区内产生所需的差频声波。参量阵差频声波的传播规律依赖于原波频率下的介质参数,而海水介质通常含有气泡和微小沙粒等,其声学参数是频率的函数。因而,海水中参量阵差频声波和同频小振幅波在传播过程中相位变化会有所差别。在参量阵的应用中需要测量这种相位差的变化,
实际应用中参量阵多发射正弦波填充的脉冲信号。关于相位差的测量在信号分析、电力电子、通讯和导航定位等技术中有广泛应用,如直接闸门计数法、过零点法、固定相移法等,这些方法均是测量两个连续正弦信号的相位差。本发明针对测量参量阵差频波与同频小振幅波信号相位差的问题,发明了一种测量脉冲信号相位差的方法。
发明内容
本发明的目的在于提供一种计算方法简单的测量参量阵差频波与同频小振幅波信号相位差的方法。
为达到上述目的,本发明采用以下技术方案:
(1)在水中吊放参量阵声源、小振幅波声源、第一水听器和第二水听器,触发信号发生器产生两路触发信号,分别对参量阵声源和小振幅波声源的信号源进行触发,触发信号的间隔Δt在1~5倍的信号源发射脉冲信号脉宽范围内,触发信号对采集器进行触发,即在触发信号到达时刻开始接收信号的采集并保存;
(2)参量阵声源和小振幅波声源在Δt时间间隔后分别发射信号,采集器采集并保存第一水听器和第二水听器接收到的信号;
(3)第一水听器和第二水听器接收到的信号均为两个脉冲信号,截取第一水听器接收信号中的第一个脉冲信号,并把截取的第一个脉冲信号与第一水听器接收的全部信号进行相关处理,找出处理结果中第一个相关峰的最大值位置,记为N1;截取第二水听器接收信号中的第一个脉冲信号,把截取的第一个脉冲信号与第二水听器接收的全部信号进行相关处理,找出处理结果中第一个相关峰的最大值位置,记为N2;
(4)截取第一水听器接收信号中的第二个脉冲信号,并把截取的第二个脉冲信号与第一水听器接收的全部信号进行相关处理,找出处理结果中第二个相关峰的最大值位置,记为N3;截取第二水听器接收信号中的第二个脉冲信号,截取的第二个脉冲信号与第二水听器接收的全部信号进行相关处理,找出处理结果中第二个相关峰的最大值位置,记为N4;
(5)对N1、N2、N3、N4进行计算,得到参量阵差频波与同频小振幅波信号相位差ΔΦ
其中:fs为采样频率;T为暂态信号参量阵差频波或同频小振幅波信号的周期;ΔΦ为所求的相位差,单位为度。
本发明具有以下有益效果:
1、可以测量两个脉冲信号传播过程中的相位差;
2、水下吊放的声源和水听器在水中的位置变化对测量结果的影响较小,对实验条件的要求不高;
3、由于采用了短时间内分别发射参量阵信号和小振幅波信号,水介质随时间的变化对测量结果影响较小;
4、实验测量数据的处理方法简单,不需要复杂的信号处理算法。
附图说明
图1是实验设备连接示意图。
图2是两个水听器的接收信号。
图3是截取的第一个脉冲信号。
图4是截取的第一个脉冲信号与第一个水听器接收信号相关处理结果。
图5是截取的第一个脉冲信号与第二个水听器接收信号相关处理结果。
图6是截取的第二个脉冲信号。
图7是截取的第二个脉冲信号与第一个水听器接收信号相关处理结果。
图8是截取的第二个脉冲信号与第二个水听器接收信号相关处理结果。
具体实施方式
下面结合附图1~6和实例对本发明作进一步详细说明。
本发明包括:
1、如图1所示连接仪器设备,并在水中吊放参量阵声源1、小振幅波声源2、水听器3和水听器4。触发信号发生器的作用是产生两路触发信号,分别对参量阵声源和小振幅波声源对应的信号源进行触发,触发信号的间隔Δt在1~5倍的信号源发射脉冲信号脉宽范围内。其中一路触发信号对采集器进行触发,即在触发信号到达时刻开始接收信号的采集并保存。
2、系统工作时,参量阵声源和小振幅波声源在Δt时间间隔后分别发射信号,采集器采集并保存水听器3和水听器4接收到的信号。
3、水听器3和水听器4接收到的信号均为两个脉冲信号,截取水听器3接收信号中的第一个脉冲信号,并把截取的脉冲信号与整个水听器3接收信号进行相关处理,找出处理结果中第一个相关峰的最大值位置,记为N1;再把截取的脉冲信号与整个水听器4接收信号进行相关处理,找出处理结果中第一个相关峰的最大值位置,记为N2。
4、截取水听器3接收信号中的第二个脉冲信号,并把截取的脉冲信号与整个水听器3接收信号进行相关处理,找出处理结果中第二个相关峰的最大值位置,记为N3;再把截取的脉冲信号与整个水听器4接收信号进行相关处理,找出处理结果中第二个相关峰的最大值位置,记为N4。
5、对N1、N2、N3、N4进行以下计算,得到参量阵差频波与同频小振幅波信号相位差ΔΦ。
其中:fs为采样频率;T为暂态信号参量阵差频波或同频小振幅波信号的周期;ΔΦ为所求的相位差,单位为度。
实施例
第一步,结合图1,触发信号发生器产生两个间隔为Δt=1.05s的触发信号,参量阵声源1发射的差频波信号和小振幅波声源2发射信号均为10Hz,发射信号脉冲宽度为0.7s。
第二步,结合如图2,水听器3和水听器4接收到差频波脉冲信号分别为5、6,接收信号采样频率为1000Hz。这里设定参量阵差频波与同频小振幅波信号在水听器3和水听器4之间距离内传播的相位差为39度。
第三步,结合图3,截取水听器3接收信号5中的第一个脉冲信号;结合图4,把截取的脉冲信号7与整个水听器3接收信号5进行相关处理,找出相关处理结果8中第一个相关峰的最大值位置N1=8442;结合图5,把截取的脉冲信号7与整个水听器4接收信号6进行相关处理,找出相关处理结果中第一个相关峰的最大值位置N2=11452。
第四步,结合图6,截取水听器3接收信号5中的第二个脉冲信号;结合图7,把截取的脉冲信号10与整个水听器3接收信号5进行相关处理,找出相关处理结果11中第二个相关峰的最大值位置N3=10194;结合图8,把截取的脉冲信号10与整个水听器4接收信号6进行相关处理,找出相关处理结果12中第二个相关峰的最大值位置N4=13214。
第五步,利用公式 得到参量阵差频波与同频小振幅波信号相位差ΔΦ=39.6°。
测量结果存在误差的主要原因是由采样频率决定,此实例中采样频率为1000Hz,则对于频率为10Hz信号而言,理论上测量的误差范围为±3.6°。为了提高测量精度,可以增加信号采样频率,例如此实例中采样频率如果提高为2000Hz时,则测量的误差范围为±1.8°。
Claims (1)
1.一种测量参量阵差频波与同频小振幅波信号相位差的方法,其特征在于:
(1)在水中吊放参量阵声源(1)、小振幅波声源(2)、第一水听器(3)和第二水听器(4),触发信号发生器产生两路触发信号,分别对参量阵声源和小振幅波声源的信号源进行触发,两路触发信号的间隔Δt在1~5倍的信号源发射脉冲信号脉宽范围内,触发信号对采集器进行触发,即在触发信号到达时刻开始接收信号的采集并保存;
(2)参量阵声源和小振幅波声源以Δt的时间间隔分别发射信号,采集器采集并保存第一水听器和第二水听器接收到的信号;
(3)第一水听器和第二水听器接收到的信号均为两个脉冲信号,截取第一水听器接收信号中的第一个脉冲信号,并把截取的第一个脉冲信号与第一水听器接收的全部信号进行相关处理,找出处理结果中第一个相关峰的最大值位置,记为N1;截取第二水听器接收信号中的第一个脉冲信号,把截取的第一个脉冲信号与第二水听器接收的全部信号进行相关处理,找出处理结果中第一个相关峰的最大值位置,记为N2;
(4)截取第一水听器接收信号中的第二个脉冲信号,并把截取的第二个脉冲信号与第一水听器接收的全部信号进行相关处理,找出处理结果中第二个相关峰的最大值位置,记为N3;截取第二水听器接收信号中的第二个脉冲信号,截取的第二个脉冲信号与第二水听器接收的全部信号进行相关处理,找出处理结果中第二个相关峰的最大值位置,记为N4;
(5)对N1、N2、N3、N4进行计算,得到参量阵差频波与同频小振幅波信号相位差ΔΦ
其中:fs为采样频率;T为暂态信号参量阵差频波或同频小振幅波信号的周期;ΔΦ为所求的相位差,单位为度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298876.9A CN104062645B (zh) | 2014-06-26 | 2014-06-26 | 一种测量参量阵差频波与同频小振幅波信号相位差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298876.9A CN104062645B (zh) | 2014-06-26 | 2014-06-26 | 一种测量参量阵差频波与同频小振幅波信号相位差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104062645A CN104062645A (zh) | 2014-09-24 |
CN104062645B true CN104062645B (zh) | 2016-12-07 |
Family
ID=51550445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410298876.9A Active CN104062645B (zh) | 2014-06-26 | 2014-06-26 | 一种测量参量阵差频波与同频小振幅波信号相位差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104062645B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104539569A (zh) * | 2015-01-09 | 2015-04-22 | 哈尔滨工程大学 | 一种基于参量阵正交频分复用编码水声通信的方法 |
CN105973447A (zh) * | 2016-05-05 | 2016-09-28 | 哈尔滨工程大学 | 直线水听器阵列幅度和相位一致性测量方法 |
CN106324278A (zh) * | 2016-08-31 | 2017-01-11 | 电子科技大学 | 一种基于声参量阵的风速测量方法 |
CN106556827B (zh) * | 2016-11-18 | 2018-10-19 | 西北工业大学 | 基于前向声散射的双发双收组网式目标探测系统及方法 |
CN110703186A (zh) * | 2019-10-22 | 2020-01-17 | 哈尔滨工程大学 | 水池近场的双水听器相位差校准测量系统及其相位差校准测量方法 |
CN112198515B (zh) * | 2020-10-13 | 2021-06-29 | 湖南国天电子科技有限公司 | 一种参量阵浅剖差频转换性能优化方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1474026A (en) * | 1973-08-10 | 1977-05-18 | Marconi Co Ltd | Pulse compression |
US4071821A (en) * | 1976-04-28 | 1978-01-31 | Honeywell Inc. | Quadrature correlation phase determining apparatus |
JP2001051046A (ja) * | 1999-08-10 | 2001-02-23 | Oki Electric Ind Co Ltd | 車間距離測定方法 |
JP4844380B2 (ja) * | 2006-12-19 | 2011-12-28 | 株式会社日立プラントテクノロジー | 立軸ポンプの内部観察装置 |
CN103760418A (zh) * | 2014-01-20 | 2014-04-30 | 中国人民解放军后勤工程学院 | 一种基于希尔伯特变换的互相关相位差测量方法 |
-
2014
- 2014-06-26 CN CN201410298876.9A patent/CN104062645B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104062645A (zh) | 2014-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104062645B (zh) | 一种测量参量阵差频波与同频小振幅波信号相位差的方法 | |
Perelli et al. | Acoustic emission localization in plates with dispersion and reverberations using sparse PZT sensors in passive mode | |
RU2474793C1 (ru) | Способ параметрического приема волн различной физической природы в морской среде | |
CN104501938B (zh) | 一种在置于空气中的矩形混响水槽内测量水下声源低频辐射声功率的方法 | |
CN102997988B (zh) | 大型潜标矢量水听器低频声学指向性的水池测试方法 | |
MX2011010162A (es) | Procesamiento de datos sismicos. | |
Jing et al. | Channel characterization of acoustic waveguides consisting of straight gas and water pipelines | |
CN104360251B (zh) | 一种变压器局部放电的超声波信号时延估计方法 | |
CN106290580B (zh) | 一种真空高低频声学测量装置及方法 | |
RU2018140060A (ru) | Способы и системы для определения направлений быстрых и медленных поперечных волн в анизотропном пласте с применением прибора для каротажа в процессе бурения | |
US10458878B2 (en) | Position determination device, leak detection system, position determination method, and computer-readable recording medium | |
Besedina et al. | Estimation of the depth of a stationary sound source in shallow water | |
Kafle et al. | Active acoustic leak detection and localization in a plastic pipe using time delay estimation | |
Kushnir et al. | Robustness of statistical algorithms for location of microseismic sources based on surface array data | |
CN106770668A (zh) | 一种用于单孔的基桩质量声波透射法检测方法 | |
AlTheyab* et al. | Reflection full-waveform inversion for inaccurate starting models | |
CN106066472A (zh) | 一种二维振速梯度水听器的被动目标相关检测方法 | |
Tashmukhambetov et al. | Three-dimensional seismic array characterization study: Experiment and modeling | |
CN104393933B (zh) | 声表面波电力测温采集器的移相叠加抗干扰方法 | |
Bentoumi et al. | Welsh DSP estimate and EMD applied to leak detection in a water distribution pipeline | |
Digulescu et al. | Underwater object tracking using time frequency signatures of acoustic signals | |
Taylor et al. | Development of a matched filter detector for acoustic signals at local distances from small explosions | |
Kwak et al. | Temporal windowing and inverse transform of the wavefield in the Laplace-Fourier domain | |
RU2561997C1 (ru) | Стенд для поверки доплеровского измерителя скорости течений | |
Angela et al. | Underwater object tracking using time frequency signatures of acoustic 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 |