CN109412724A - A kind of high-accuracy water sound communication channel modeling method - Google Patents

A kind of high-accuracy water sound communication channel modeling method Download PDF

Info

Publication number
CN109412724A
CN109412724A CN201811002662.7A CN201811002662A CN109412724A CN 109412724 A CN109412724 A CN 109412724A CN 201811002662 A CN201811002662 A CN 201811002662A CN 109412724 A CN109412724 A CN 109412724A
Authority
CN
China
Prior art keywords
path
model
shifting
delay
communication channel
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.)
Withdrawn
Application number
CN201811002662.7A
Other languages
Chinese (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.)
715th Research Institute of CSIC
Original Assignee
715th Research Institute of CSIC
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 715th Research Institute of CSIC filed Critical 715th Research Institute of CSIC
Priority to CN201811002662.7A priority Critical patent/CN109412724A/en
Publication of CN109412724A publication Critical patent/CN109412724A/en
Withdrawn legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B11/00—Transmission systems employing sonic, ultrasonic or infrasonic waves
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Electromagnetism (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a kind of high-accuracy water sound communication channel modeling methods, analyze influence of the shifting to multi-path delay and Doppler, and this influence is considered in the time varying channel model of PARASTOO, to optimize its time-varying channel model, high-accuracy water sound communication Channel Modeling is realized in the case where not increasing calculation amount.And this method feasibility is demonstrated with WAFO modeling shifting.The first step of the present invention constructs channel large scale variation model, with the random motion of single order AR process simulation sound source and receiver；Second step requires to calculate multi-path delay and dampening information with bellhop model or mirror image model according to simulated environment；Third step analyzes influence of the shifting to multi-path delay and Doppler with a kind of approximate model according to shifting data；4th step analyzes small scale scattering model, and each path is scattered as a plurality of subpath, the time delay and Doppler contribution of subpath are analyzed；5th step analyzes the Doppler contribution of transmit-receive platform movement；6th step couples multipath, and considers the influence of shifting, constructs time varying channel shock response.

Description

A kind of high-accuracy water sound communication channel modeling method
Technical field
The present invention relates to field of acoustics, and in particular to underwater sound communication Channel Modeling technology, mainly a kind of high-precision underwater sound Communication channel modeling method.
Background technique
When constructing underwater sound communication system, such as the analysis design of the waveform of signal of communication, bandwidth and power factors Need to consider the influence of underwater acoustic channel.So high-precision underwater acoustic channel model is most important to communication system building. BELLHOP is a kind of model for analyzing waveguide sonic pressure field.The model is using very extensive, in the case where set marine environment, applies Ray theory analyzes channel function, propagation loss etc., but this model does not consider the time-varying characteristics of practical ocean environment parameter.Closely The attention rate of Nian Lai, time varying channel simulation model are higher and higher.PORTER etc. proposes a kind of time-varying letter based on BELLHOP Channel estimation method, this method analyze influence of the time varying channel to transmitting signal, but the algorithm by tracking multiple associated beams Calculation amount is huge.PARASTOO proposes a kind of underwater sound communication channel statistical modeling method, by the not true of actual environment parameter Reality is classified as two kinds, is large scale effect and Small-scale fading respectively, and with AR process model building.This method constructs time-varying Underwater acoustic channel model, and reduce calculation amount.Shifting is to cause a key factor of channel time-varying, but due to sea The fluctuations period is very short, if BELLHOP model is called to carry out analog sea variation repeatedly, will increase considerably calculation amount, so Shifting model is not considered in the method analytic process of PARASTOO, only one fluctuating velocity of setting is used to analyze Doppler, Influence of the relief height to time delay etc. is not considered.
Summary of the invention
It is an object of the invention to overcome the shortcomings of the prior art, and a kind of high-accuracy water sound communication channel is provided and is built Mould method can simulate ocean time varying channel, and analysis considers shifting to multichannel in the time varying channel model of PARASTOO The influence of diameter time delay and Doppler optimizes its time-varying channel model, and the high-precision underwater sound is realized in the case where not increasing calculation amount Communication channel modeling.
The object of the present invention is achieved by the following technical solutions.This high-accuracy water sound communication Channel Modeling side Method comprises the following steps:
(1) channel large scale variation model is constructed, with the random motion of single order AR process simulation sound source and receiver.
Transmit-receive platform is one of the main reason for causing large scale to change apart from mismatch.The influence reaction of large scale variation exists Each propagated distance and time delay.The variation of path length isΔlpRandom variation.The knots modification of time delay by ΔlpIt determines.Path gain hpUse lpIt is expressed as
Wherein a0For the absorption coefficient at signal center frequency.Use during actual analysis single order AR process simulation sound source and The large scale of receiver changes.
(2) it is required to calculate multi-path delay and dampening information with bellhop model or mirror image model according to simulated environment.
Etc. velocities of sound environment with mirror image model come analyzing multiple diameter time delay and dampening information, used when the velocity of sound is with change in depth Bellhop model is analyzed.
(3) according to shifting data, shifting is analyzed with approximate model to the shadow of multi-path delay and Doppler It rings.
Long dotted line indicates certain sea surface reflection under the smooth sea environment that BELLHOP model or mirror image model calculate in Fig. 1 Path, with time change, shifting, this path length changes, and knots modification is Δ, can be approximate with expression,
Δ(u,t)≈h(u,t)cos(θ) (3)
Wherein path is in the horizontal distance of sea surface reflection thus by u, and the angle on sea is incident in path to θ thus.Path length Knots modification determines time delay knots modification, so the propagation delay in this path is after shifting variation
The Doppler contribution of shifting determines that expression is as follows by fluctuating rate
W (t, u) is in t moment, and distance is the relief height on sea at u.
(4) analyze small scale scattering model, each path scatter analyze for a plurality of subpath the time delay of subpath with it is more Pu Le influences.
Small dimensional variation is mainly caused by the scattering of each path.Each path is scattering into a plurality of subpath when analysis, every The amplitude and time delay in path have fluctuating.Scattering subpath amplitude indicated with normal distribution, Shi Yanyong AR process description.Pth Scattering coefficient of the paths after small dimensional variation be
It is the random delay of each path scattering subpath.
(5) Doppler contribution of analysis transmit-receive platform movement.
Transmit-receive platform relative velocity expression formula is
vdp=vtdcos(θp-θtd)-vrdcos(θp+θrd) (7)
Wherein vtd,θtd,vrd,θrdFor the movement velocity and angle of target and receiver.θpFor the destination path direction of motion.
Doppler shift factor a is analyzed according to this relative velocityp=vdp/c。
(6) multipath is coupled, and considers the influence of shifting, construct time varying channel shock response.
After coupling,
Wherein hpFor the path gain of pth paths.For reference path.τ'pFor multi-path delay, wherein including Influence of the shifting to multi-path delay,For consider transmitting-receiving Doppler shift (platform movement and shifting) with Small dimensional variation coefficient afterwards.
The invention has the benefit that
(1) present invention analyzes shifting and reaches to multipath on the basis of the time varying channel model of PARASTOO The influence of signal time delay and Doppler, and received in construction and consider that this is influenced, it can be achieved that high-accuracy water sound communication in signal process Channel Modeling, construction high fidelity receive signal.
(2) the time varying channel model based on PARASTOO considers influence of the shifting to channel impulse response, proposes A kind of high-precision underwater acoustic channel modeling method.This method analyzes influence of the shifting to multi-path delay and Doppler, and This influence is considered in the time varying channel model of PARASTOO, to optimize its time-varying channel model, is not increasing calculation amount In the case of realize high-accuracy water sound communication Channel Modeling.And this method feasibility is demonstrated with WAFO modeling shifting.
Detailed description of the invention
Fig. 1: shifting approximation schematic diagram.
Fig. 2: b1s1 path proximity result and bellhop operation result propagation time are analyzed.
Fig. 3: the time varying channel shock response of Bellhop is run repeatedly.
Fig. 4: the shock response that approximation obtains.
The approximate shifting of Fig. 5: WAFO model.
Fig. 6: when constant channel impulse response.
Fig. 7: PARASTOO model channel impulse response.
Fig. 8: WAFO shifting model channel impulse response.
Fig. 9: Channel Modeling system schematic.
Specific embodiment
Below in conjunction with specific example and attached drawing, the present invention will be described in detail:
As shown, this high-accuracy water sound communication channel modeling method, comprises the following steps:
(1) channel large scale variation model is constructed, with the random motion of single order AR process simulation sound source and receiver.Greatly Scale effect single order AR process model building, AR coefficient are 0.9, each random offset Normal Distribution.
(2) it is required to calculate multi-path delay and dampening information with bellhop model according to simulated environment.
Bellhop model analysis, simulated conditions are as follows: waveguide depth 100m, depth of receiver are respectively are used in simulation process 50m and 20m, transmitting-receiving spacing are 2000m.The bottom parameters velocity of sound is 1800m/s, density 1.6g/cm3。
(3) according to shifting data, shifting is analyzed with approximate model to the shadow of multi-path delay and Doppler It rings.
The present invention simulates extra large surface random fluctuation movement with WAFO model.Its computation model is a stable Gaussian mistake Journey is calculated in certain location point, in its shifting height of set time.Expression formula is as follows
Wherein ΘjFor random phase fluctuating.RjFor L-S distribution amplitude, SjFor weighting coefficient.One obtained using this model A shifting variation diagram is as shown in Figure 5.The corresponding shifting situation single order AR process description of result shown in Fig. 2.According to Formula (4) and (5) calculate influence of the shifting to time delay and Doppler.
(4) analyze small scale scattering model, each path scatter analyze for a plurality of subpath the time delay of subpath with it is more Pu Le influences.
Each path is set in analytic process and scatters subpath 20, scatters the amplitude Gaussian distributed of subpath, Value is hp/ 20, variance 1e-6.Identical scattering subpath has correlation in the propagation delay of different time, is built with AR process Mould.Wherein random offset is simulated with Gaussian process.
(5) the Doppler contribution analysis of transmit-receive platform movement.
As long as analyzing the opposing stationary situation of transmit-receive platform in simulation process.
(6) multipath is coupled, and considers the influence of shifting, construct time varying channel shock response.
Fig. 2 provides the propagation delay by a sea surface reflection and a sub-bottom reflection path (b1s1) in different time, BELLHOP model is run with the approximation method in this achievement and repeatedly respectively to estimate multi-path delay.As can be seen from the figure Result with approximate method and theoretical calculation is very close, error very little.It is compared to each other for Different Results, Fig. 3-Fig. 6 does not have Consider that Small-scale fading influences.Fig. 7-Fig. 8 is influenced comprising Small-scale fading.Fig. 3 is that operation BELLHOP model obtains not repeatedly With time channel impulse response, Fig. 4 is the shock response obtained with approximate analysis as a result, p in figure0Represent direct sound wave, psIt represents Path through a sea surface reflection, pbRepresent the path through a sub-bottom reflection.Comparing two figures, to can be seen that two methods each Time point shock response coincide substantially.Fig. 6-8 constant channel impulse response, PARASTOO model Channel Impulse sound when being respectively It should be with WAFO shifting model channel impulse response.Constant channel when Fig. 6 is, so its channel impulse response is not at any time Variation.The corresponding shock response of Fig. 7 is by large scale variation and the processing of small dimensional variation, result fluctuations at any time, phase Than more tallying with the actual situation in Fig. 6, and comparison diagram 7 and Fig. 8 can be seen that under WAFO shifting model, through sea surface reflection Path fluctuations it is bigger, and p0And pbWithout sea surface reflection, this two paths is not influenced by shifting.So Fig. 8 Shock response more meet time varying channel characteristic.Fig. 9 is whole system operation schematic diagram.
The present invention is not limited to the above-described embodiments, all to use institute of the present invention no matter embodiments thereof make any variation The implementation structure design of offer, is all a kind of deformation of the invention, is considered as within the protection scope of invention.

Claims (7)

1. a kind of high-accuracy water sound communication channel modeling method, it is characterised in that: specific step is as follows for this method:
(1) channel large scale variation model is constructed, with the random motion of single order AR process simulation sound source and receiver；
(2) it is required to calculate multi-path delay and dampening information with bellhop model or mirror image model according to simulated environment；
(3) according to shifting data, shifting is analyzed with a kind of approximate model to the shadow of multi-path delay and Doppler It rings；
(4) small scale scattering model is analyzed, each path is scattered as a plurality of subpath, time delay and the Doppler of subpath are analyzed It influences；
(5) Doppler contribution of analysis transmit-receive platform movement；
(6) multipath is coupled, and considers the influence of shifting, construct time varying channel shock response.
2. high-accuracy water sound communication channel modeling method according to claim 1, it is characterised in that: in step (1), greatly The influence of dimensional variation is reacted in each propagated distance and time delay, and the variation of path length isΔlpAt random Variation, the knots modification of time delay is by Δ lpIt determines, path gain hpUse lpIt is expressed as
Wherein a0For the absorption coefficient at signal center frequency, with single order AR process simulation sound source and receiver in analytic process Large scale variation.
3. high-accuracy water sound communication channel modeling method according to claim 1, it is characterised in that: in step (2), etc. Velocity of sound environment with mirror image model come analyzing multiple diameter time delay and dampening information, the velocity of sound with when change in depth with bellhop model come Analysis.
4. high-accuracy water sound communication channel modeling method according to claim 1, it is characterised in that: in step (3), Certain sea surface reflection path under the smooth sea environment that BELLHOP model or mirror image model calculate, with time change, sea is risen Volt, this path length change, and knots modification is Δ, can be approximate with expression,
Δ(u,t)≈h(u,t)cos(θ) (3)
Wherein path is in the horizontal distance of sea surface reflection thus by u, and the angle on sea is incident in path to θ thus, and path length changes Amount determines time delay knots modification, so after shifting variation, the propagation delay in this path are as follows:
The Doppler contribution of shifting determines that expression is as follows by fluctuating rate
W (t, u) is in t moment, and distance is the relief height on sea at u.
5. high-accuracy water sound communication channel modeling method according to claim 1, it is characterised in that:, will in step (4) Each path is scattering into a plurality of subpath, and the amplitude and time delay of each path have fluctuating, scatters the amplitude normal state of subpath Distribution expression, Shi Yanyong AR process description, scattering coefficient of the pth paths after small dimensional variation are
It is the random delay of each path scattering subpath.
6. high-accuracy water sound communication channel modeling method according to claim 1, it is characterised in that: in step (5),
Transmit-receive platform relative velocity expression formula are as follows:
vdp=vtdcos(θp-θtd)-vrdcos(θp+θrd) (7)
Wherein vtd,θtd,vrd,θrdFor the movement velocity and angle of target and receiver, θpFor the destination path direction of motion, according to This relative velocity analyzes Doppler shift factor ap=vdp/c。
7. high-accuracy water sound communication channel modeling method according to claim 1, it is characterised in that: in step (6), coupling After conjunction,
Wherein hpFor the path gain of pth paths,For reference path, τ 'pFor multi-path delay, wherein including sea The influence to rise and fall to multi-path delay,To consider to receive and dispatch Doppler shift, later small dimensional variation coefficient.
CN201811002662.7A 2018-08-30 2018-08-30 A kind of high-accuracy water sound communication channel modeling method Withdrawn CN109412724A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811002662.7A CN109412724A (en) 2018-08-30 2018-08-30 A kind of high-accuracy water sound communication channel modeling method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811002662.7A CN109412724A (en) 2018-08-30 2018-08-30 A kind of high-accuracy water sound communication channel modeling method

Publications (1)

Publication Number Publication Date
CN109412724A true CN109412724A (en) 2019-03-01

Family

ID=65463615

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811002662.7A Withdrawn CN109412724A (en) 2018-08-30 2018-08-30 A kind of high-accuracy water sound communication channel modeling method

Country Status (1)

Country Link
CN (1) CN109412724A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117220808A (en) * 2023-10-08 2023-12-12 哈尔滨工程大学 Air-sea cross-medium channel modeling method based on environment driving and task driving

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103634262A (en) * 2012-08-29 2014-03-12 中国科学院声学研究所 Underwater acoustic coherent communication self-adaptive phase compensation method
KR101405247B1 (en) * 2014-02-17 2014-06-17 연세대학교 산학협력단 Apparatus and method for predicting transmission loss of underwater acoustic channel
CN104734795A (en) * 2015-02-11 2015-06-24 哈尔滨工业大学深圳研究生院 Time-variant underwater acoustic channel capacity simulation model
CN106411438A (en) * 2016-11-02 2017-02-15 东北农业大学 Shallow water time-varying multi-path underwater acoustic channel modeling method
CN108650043A (en) * 2018-06-29 2018-10-12 中国船舶重工集团公司第七〇五研究所 A kind of high-accuracy water sound communication channel modeling method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103634262A (en) * 2012-08-29 2014-03-12 中国科学院声学研究所 Underwater acoustic coherent communication self-adaptive phase compensation method
KR101405247B1 (en) * 2014-02-17 2014-06-17 연세대학교 산학협력단 Apparatus and method for predicting transmission loss of underwater acoustic channel
CN104734795A (en) * 2015-02-11 2015-06-24 哈尔滨工业大学深圳研究生院 Time-variant underwater acoustic channel capacity simulation model
CN106411438A (en) * 2016-11-02 2017-02-15 东北农业大学 Shallow water time-varying multi-path underwater acoustic channel modeling method
CN108650043A (en) * 2018-06-29 2018-10-12 中国船舶重工集团公司第七〇五研究所 A kind of high-accuracy water sound communication channel modeling method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIU SHENGXING: "Impact of Sea Waves on Performance of Shallow", 《IEEE/OES CHINA OCEAN ACOUSTICS SYMPOSIUM》 *
俞丹丽: "海上数字通信信道编解码技术研究", 《中国优秀硕士学位论文全文库 信息科技辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117220808A (en) * 2023-10-08 2023-12-12 哈尔滨工程大学 Air-sea cross-medium channel modeling method based on environment driving and task driving
CN117220808B (en) * 2023-10-08 2024-04-02 哈尔滨工程大学 Air-sea cross-medium channel modeling method based on environment driving and task driving

Similar Documents

Publication Publication Date Title
Hodges Underwater acoustics: Analysis, design and performance of sonar
CN104734795B (en) A kind of time- variant channel Capacity Simulation model
CN102590804B (en) Overland testing system of Doppler sonar and testing method thereof
CN110287535B (en) Simulation calculation method for motion response of underwater towing system under high sea condition
CN101482614A (en) Sound propagation velocity modeling method, apparatus and system
CN111983619B (en) Underwater acoustic target forward scattering acoustic disturbance positioning method based on transfer learning
CN110133627B (en) Method for optimizing array element position calibration measurement point spacing of underwater acoustic positioning navigation system
Nosal et al. Sperm whale three-dimensional track, swim orientation, beam pattern, and click levels observed on bottom-mounted hydrophones
CN106411438A (en) Shallow water time-varying multi-path underwater acoustic channel modeling method
CN105277932B (en) A kind of Doppler correction method in the Wave beam forming based on down coversion
CN108650043A (en) A kind of high-accuracy water sound communication channel modeling method
CN108572349B (en) Sound source depth setting method based on model calculation under deep sea environment
CN109412724A (en) A kind of high-accuracy water sound communication channel modeling method
Durofchalk et al. Data driven source localization using a library of nearby shipping sources of opportunity
CN111650574B (en) Underwater space-time self-adaptive processing method and system based on sparse recovery
CN103487793A (en) Broadband reverberation waveform simulation method based on normal mode theory
CN112305502A (en) Water surface and underwater sound source binary discrimination method based on array invariants
CN1995932A (en) Sound velocity profile inversion method based on buoyage
CN113126029B (en) Multi-sensor pulse sound source positioning method suitable for deep sea reliable acoustic path environment
Guo et al. Target depth estimation by deep neural network based on acoustic interference structure in deep water
CN111965601A (en) Underwater sound source passive positioning method based on nuclear extreme learning machine
Wang et al. Identification of NLOS Acoustic Signal Using CNN and Bi-LSTM.
CN111351561B (en) DSP-based multi-channel multi-path underwater acoustic channel real-time simulation method
CN115047408B (en) Underwater multi-sound-source positioning method based on single-layer large convolution kernel neural network
Hovem et al. Modeling Low-Frequency Anthropogenic Noise in the Oceans: A Comparison of Propagation Models

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
WW01 Invention patent application withdrawn after publication

Application publication date: 20190301

WW01 Invention patent application withdrawn after publication