CN111162845A - Sea area sight distance channel generation method - Google Patents
Sea area sight distance channel generation method Download PDFInfo
- Publication number
- CN111162845A CN111162845A CN202010039118.0A CN202010039118A CN111162845A CN 111162845 A CN111162845 A CN 111162845A CN 202010039118 A CN202010039118 A CN 202010039118A CN 111162845 A CN111162845 A CN 111162845A
- Authority
- CN
- China
- Prior art keywords
- ray
- channel
- path
- parameters
- propagation
- 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.)
- Pending
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
- H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
-
- 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
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a sea area sight distance channel generation method, which realizes sea area sight distance channel generation through the steps of parameter setting, ray propagation type discrimination, formula rewriting, scattering area generation, height change of a receiving end, parameter calculation, channel generation and the like. The random ray tracing method is used for assuming actual complex factors such as ocean wave change and the like to improve the accuracy of generating the channel, reduce the complexity of generating the channel and enhance the applicability of the ray tracing method, so that the sea area sight distance channel which is more in line with the actual situation is generated.
Description
Technical Field
The invention relates to the field of communication, in particular to a sea area line-of-sight channel generation method.
Background
With the ever-increasing marine economy and the increasing activity of marine activities, sea area mobile communication has received a great deal of attention. In order to meet the requirements of future mobile communication applications, high-rate, wide-coverage and low-delay sea communication services need to be developed. As part of basic research, modeling of wireless channels at sea has been studied from different perspectives, and while most of the existing work has focused on studying large-scale propagation path loss in marine-based environments, fast fading channel models have not been well studied and developed.
In the prior art, the channel research is mainly modeled by Ray Tracing (Ray-Tracing). Although ray tracing methods can greatly reduce the complexity of modeling, it is often not possible to accurately describe the propagation path of each ray in a rich scattering environment. But in a marine propagation environment, scattering is typically sparse, which makes it possible to model the channel with low complexity using ray tracing techniques. The ray tracing method can calculate specific ray tracks in the multipath channel modeling, and has obvious advantages for small-scale parameters. Due to the lack of consideration of complex ocean waves on ray tracing method modeling in the prior art, the ray tracing method is still in a blank stage in the sea area sight distance channel modeling field.
Disclosure of Invention
The invention provides a sea area apparent distance channel generation method aiming at solving the technical problem of how to generate a sea area apparent distance channel with low complexity and strong applicability.
In order to solve the technical problems, the technical scheme of the invention is as follows: a method for generating a sea area sight distance channel comprises the following steps:
(1) setting parameters, namely setting natural environment parameters and working system parameters;
(2) judging the ray propagation type, dividing the ray propagation type, determining a ray tracing method calculation formula of each divided type of ray, and performing the next step;
(3) rewriting a formula, namely rewriting a ray tracing method calculation formula of each type of divided ray propagation into a formula taking the moment t as an unknown quantity according to the moving speed of the receiving antenna ship in the ocean;
(4) generating a scattering area, determining the number and the positions of effective scattering points in the sea level, and performing the next step after determining the number and the positions;
(5) the height change of the receiving end is evaluated to the height of the receiving end antenna according to the height change caused by the fluctuation of the sea waves, and the next step is carried out;
(6) calculating parameters, namely calculating to obtain the length of a direct path, the sea level reflection times and the length of a reflection ray propagation path, recording the parameters as the parameters of the Kth reflection arrival path, and performing the next step;
(7) and (4) generating a channel, namely bringing the channel parameters of the arrival path subjected to statistics into a channel generation module to form a sea area line-of-sight channel.
Further, the natural environment parameter includes a wave height.
Further, the working system parameters comprise the height of the transmitting antenna, the height of the receiving antenna, the distance between the transmitting antenna and the receiving antenna and the moving speed of the ship.
Further, the ray tracing method for each type of ray has a calculation formula as follows:
wherein d isLOS，dRef,0And dRef,kRespectively representing the lengths of a direct path, a specular reflection path and a scattering path; k is the total number of scattering points;
is the phase shift caused by the propagation delay corresponding to the length of each propagation path, f is the carrier frequency, c is the electricalThe magnetic wave propagation speed; gamma is the reflection loss factor, amplitude and the sum factor of gamma
In direct proportion.
Further, the effective scattering point position in the step (4) is calculated by the following formula: c. CRef,k(t)＝αrcosθu⊥+rsinθu+cRef,0(t) where R is a value of u (0, R), θ is a value of u (0,2 π), α is a compression factor that controls the scattering region.
By adopting the technical scheme, the random ray tracing method is utilized to assume actual complex factors such as ocean wave change and the like so as to improve the accuracy of generating the channel and reduce the complexity of generating the channel; the applicability of the ray tracing method is enhanced, and therefore the sea area sight distance channel which is more in line with the actual situation is generated.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention;
FIG. 2 is a diagram illustrating a multipath propagation map of a channel model and a determination of different ray propagation paths according to an embodiment of the present invention;
fig. 3 is a schematic diagram of determining a scattering region according to an embodiment of the present invention.
Detailed Description
The following further describes embodiments of the present invention with reference to the drawings. It should be noted that the description of the embodiments is provided to help understanding of the present invention, but the present invention is not limited thereto. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
Referring to fig. 1, the present embodiment discloses a method for generating a sea-area line-of-sight channel, which includes the following steps:
(1) setting parameters, namely setting natural environment parameters and working system parameters;
(2) judging the ray propagation type, dividing the ray propagation type, determining a ray tracing method calculation formula of each divided type of ray, and performing the next step;
(3) rewriting a formula, namely rewriting a ray tracing method calculation formula of each type of divided ray propagation into a formula taking the moment t as an unknown quantity according to the moving speed of the receiving antenna ship in the ocean;
(4) generating a scattering area, determining the number and the positions of effective scattering points in the sea level, and performing the next step after determining the number and the positions;
(5) the height change of the receiving end is evaluated to the height of the receiving end antenna according to the height change caused by the fluctuation of the sea waves, and the next step is carried out;
(6) calculating parameters, namely calculating to obtain the length of a direct path, the sea level reflection times and the length of a reflection ray propagation path, recording the parameters as the parameters of the Kth reflection arrival path, and performing the next step;
(7) and (4) generating a channel, namely bringing the channel parameters of the arrival path subjected to statistics into a channel generation module to form a sea area line-of-sight channel.
Further, the natural environment parameter includes a wave height.
Further, the working system parameters comprise the height of the transmitting antenna, the height of the receiving antenna, the distance between the transmitting antenna and the receiving antenna and the moving speed of the ship.
Further, the ray tracing method for each type of ray has a calculation formula as follows:
wherein d isLOS，dRef,0And dRef,kRespectively representing the lengths of a direct path, a specular reflection path and a scattering path; k is the total number of scattering points;
is a phase shift caused by propagation delay corresponding to the length of each propagation path, f is the carrier frequency, c is the electromagnetic wave propagation speed; gamma is the reflection loss factor, amplitude and the sum factor of gamma
In direct proportion.
Further, the effective scattering point position in the step (4) is calculated by the following formula: c. CRef,k(t)＝αrcosθu⊥+rsinθu+cRef,0(t) where R is a value of u (0, R), θ is a value of u (0,2 π), α is a compression factor that controls the scattering region.
In this embodiment, specifically taking the spatial multipath channel (SCM) model as an example, the simulation needs to generate a channel described by the following formula:
wherein:
dLOS，dRef,0and dRef,kThe lengths of the direct path, specular reflection path, and scattering path are represented, respectively.
K is the total number of scattering points and is influenced by the height of the transmitting and receiving ends, the distance between the transmitting and receiving ends and the wave factor.
Gamma is the reflection loss factor, amplitude and the sum factor of gamma
In direct proportion.
The multipath propagation diagram of the channel model and the schematic diagram for distinguishing different ray propagation paths shown in fig. 2; under the condition that a single antenna is adopted at both the transmitting end and the receiving end, the method for calculating the positions of the receiving end, the specular reflection point and the scattering point and the lengths of all paths in the sea area stadium channel comprises the following steps:
the method comprises the following steps: coordinates of the transmitting and receiving end positions, defining a transmitting end position C under the three-dimensional coordinatesTX＝(xTX,yTX,zTX) Receiving end vessel initial position CRX(0)＝(xRX(0),yRX(0),hRX)，hRXIs the height of the receiving antenna.
Step two: when the ship moves, the real-time position of the receiving end changes according to the moving speed of the ship as followsRX(t)＝CRX(0) + vt, where v ═ v (v)x,vy,0)。
Step three: adding wave factors, wherein the height of the receiving end randomly changes along with the fluctuation of waves:
Step four: position coordinates of specular reflection points, its position CRef,0(t)＝(xRef,0(t),yRef,0(t),zRef,0(t)) the coordinates of the transmitting/receiving end positions are calculated by a geometric method:
zRef,0(t)＝0
step five: determining scattering area and scattering point positionThe scattering point area is shown in FIG. 3, the large area is divided into many small areas dA, each dA is accompanied by a random β defined as the angle between the normal vector and the vertical vector on dA, and is Gaussian-distributed, so that when β value is large, the probability that the corresponding dA is effective is small, the scattering point area is roughly an ellipse, and when the transceiving distance is larger, the scattering point area is correspondingly compressed into a narrow band, the K diffuse reflection point coordinates are generated as follows, firstly, the K diffuse reflection point coordinates are expressed in (x) xRef,0(t),yRef,0(t)) as a center, randomly generating a point in a circular area with R as a radius, wherein the value of R determines the size and the boundary of the scattering area; the coordinates of this diffuse reflection point are then at u⊥The "compression" in the direction, which is a unit vector perpendicular to the TX-RX direction, on the basis of which the coordinates c of the K diffuse reflection points are obtainedRef,k(t)＝αrcosθu⊥+rsinθu+cRef,0(t) where R to u (0, R) and θ to u (0,2 π), α are the compression factors that control the shape of this scattering region.
Step six: the length of each propagation path is calculated and can be obtained by the following formula:
dLOS(t)＝||cTX-cRX(t)||
dref,k(t)＝||cTX-cref,k(t)||+||cRX(t)-cref,k(t)||,k＝0,1,...,K.
the random ray tracing method is used for assuming actual complex factors such as ocean wave change and the like to improve the accuracy of generating channels and reduce the complexity of generating the channels. The applicability of the ray tracing method is enhanced, and therefore the sea area sight distance channel which is more in line with the actual situation is generated.
The embodiments of the present invention have been described in detail with reference to the accompanying drawings, but the present invention is not limited to the described embodiments. It will be apparent to those skilled in the art that various changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, and the scope of protection is still within the scope of the invention.
Claims (5)
1. A method for generating a sea area line-of-sight channel is characterized by comprising the following steps:
(1) setting parameters, namely setting natural environment parameters and working system parameters;
(2) judging the ray propagation type, dividing the ray propagation type, determining a ray tracing method calculation formula of each divided type of ray, and performing the next step;
(3) rewriting a formula, namely rewriting a ray tracing method calculation formula of divided three types of ray propagation into a formula taking the moment t as an unknown quantity according to the moving speed of the receiving antenna ship in the ocean;
(4) generating a scattering area, determining the number and the positions of effective scattering points in the sea level, and performing the next step after determining the number and the positions;
(5) assigning the height change caused by the fluctuation of the sea waves into the height of the receiving end antenna according to the height change of the receiving end, and carrying out the next step;
(6) calculating parameters, namely calculating to obtain the length of a direct path, the sea level reflection times and the length of a reflection ray propagation path, recording the parameters as the parameters of the Kth reflection arrival path, and performing the next step;
(7) and (4) generating a channel, namely bringing the channel parameters of the arrival path subjected to statistics into a channel generation module to form a sea area line-of-sight channel.
2. The method according to claim 1, wherein the natural environment parameter comprises sea wave height.
3. The method as claimed in claim 2, wherein the operating system parameters include transmitting antenna height, receiving antenna height, transmitting and receiving antenna spacing and vessel moving speed.
4. The method as claimed in claim 3, wherein the ray tracing method for each type of ray has a formula:
wherein d isLOS，dRef,0And dRef,kRespectively representing the lengths of a direct path, a specular reflection path and a scattering path; k is the total number of scattering points;
is a phase shift caused by propagation delay corresponding to the length of each propagation path, f is the carrier frequency, c is the electromagnetic wave propagation speed; gamma is the reflection loss factor, amplitude and the sum factor of gamma
In direct proportion.
5. The method according to claim 4, wherein the effective scattering point position in step (4) is calculated by the following formula: c. CRef,k(t)＝αrcosθu⊥+rsinθu+cRef,0(t) where R is a value of u (0, R), θ is a value of u (0,2 π), α is a compression factor that controls the scattering region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010039118.0A CN111162845A (en) | 2020-01-14 | 2020-01-14 | Sea area sight distance channel generation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010039118.0A CN111162845A (en) | 2020-01-14 | 2020-01-14 | Sea area sight distance channel generation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111162845A true CN111162845A (en) | 2020-05-15 |
Family
ID=70563398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010039118.0A Pending CN111162845A (en) | 2020-01-14 | 2020-01-14 | Sea area sight distance channel generation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111162845A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114244459A (en) * | 2021-11-29 | 2022-03-25 | 西安电子科技大学 | Method, system, medium, equipment and terminal for simulating wireless channel in marine environment |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101504289B1 (en) * | 2014-02-26 | 2015-03-30 | 연세대학교 산학협력단 | Apparatus and method for generating underwater channel response |
| CN104656078A (en) * | 2013-11-25 | 2015-05-27 | 北京环境特性研究所 | Sea surface ship method based on ray tracing |
| CN105388465A (en) * | 2015-12-17 | 2016-03-09 | 西安电子科技大学 | Sea clutter simulation method based on sea wave spectrum model |
| CN108540248A (en) * | 2018-03-16 | 2018-09-14 | 西安电子科技大学 | Maritime wireless communication dynamic multi-path channel model and method, wireless communication system |
| CN109408919A (en) * | 2018-10-11 | 2019-03-01 | 南通先进通信技术研究院有限公司 | A kind of sea area evaporation waveguide channel generation method |
| CN110208764A (en) * | 2019-05-05 | 2019-09-06 | 南京航空航天大学 | Dynamic sea large scene echo simulation method based on electromagnetism Scattering Calculation |
-
2020
- 2020-01-14 CN CN202010039118.0A patent/CN111162845A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104656078A (en) * | 2013-11-25 | 2015-05-27 | 北京环境特性研究所 | Sea surface ship method based on ray tracing |
| KR101504289B1 (en) * | 2014-02-26 | 2015-03-30 | 연세대학교 산학협력단 | Apparatus and method for generating underwater channel response |
| CN105388465A (en) * | 2015-12-17 | 2016-03-09 | 西安电子科技大学 | Sea clutter simulation method based on sea wave spectrum model |
| CN108540248A (en) * | 2018-03-16 | 2018-09-14 | 西安电子科技大学 | Maritime wireless communication dynamic multi-path channel model and method, wireless communication system |
| CN109408919A (en) * | 2018-10-11 | 2019-03-01 | 南通先进通信技术研究院有限公司 | A kind of sea area evaporation waveguide channel generation method |
| CN110208764A (en) * | 2019-05-05 | 2019-09-06 | 南京航空航天大学 | Dynamic sea large scene echo simulation method based on electromagnetism Scattering Calculation |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114244459A (en) * | 2021-11-29 | 2022-03-25 | 西安电子科技大学 | Method, system, medium, equipment and terminal for simulating wireless channel in marine environment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109039506B (en) | A kind of underwater mobile channel emulation mode | |
| US20190250244A1 (en) | Underwater source node positioning method | |
| CN112070894B (en) | Real environment navigation multipath real-time simulation method, device, medium and electronic equipment | |
| CN105430664B (en) | It is a kind of to be fitted the method and apparatus that path loss is propagated in prediction based on classification | |
| Fan et al. | A structured bidirectional LSTM deep learning method for 3D terahertz indoor localization | |
| Cheng et al. | Node selection algorithm for underwater acoustic sensor network based on particle swarm optimization | |
| CN110399680B (en) | Shallow sea elastic structure radiation sound field calculation method | |
| Qiao et al. | Distributed localization based on signal propagation loss for underwater sensor networks | |
| CN113532437B (en) | Ultra-wideband indoor target positioning method based on multipath utilization | |
| CN108540248B (en) | Dynamic multi-path channel model for offshore wireless communication and construction method | |
| WO2023213052A1 (en) | Deep-sea sound source localization method, computer device and storage medium | |
| CN111162845A (en) | Sea area sight distance channel generation method | |
| Li et al. | Cramer-rao lower bound analysis of data fusion for fingerprinting localization in non-line-of-sight environments | |
| CN109408919B (en) | Sea area evaporation waveguide channel generation method | |
| CN115629389A (en) | Target positioning and error analysis method based on three-dimensional space multi-base sonar | |
| CN113534219B (en) | Beidou positioning outdoor target method based on multipath utilization | |
| CN112511199B (en) | Low-complexity 3D-MIMO related channel fast generation method | |
| Nuckols | Implementation of geometrically based single-bounce models for simulation of angle-of-arrival of multipath delay components in the wireless channel simulation tools, SMRCIM and SIRCIM | |
| CN115397012A (en) | Realization method of UWB positioning tracking system based on TWR-TDOA estimation and MPGA layout optimization | |
| CN106918809A (en) | Quick interference synthetic aperture sonar original echo time-domain simulation method | |
| CN110865359A (en) | Underwater acoustic ranging method based on received signal strength | |
| Lv et al. | Signal Design and Processing for Underwater Acoustic Positioning and Communication Integrated System | |
| Wu et al. | A Weighted Random Forest Based Positioning Algorithm for 6G Indoor Communications | |
| Masmitjà Rusiñol et al. | Range-only underwater target localization: error characterization | |
| CN102035608B (en) | Method for analyzing spatial difference of group path near receiving point |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Jue Inventor after: Ding Rui Inventor after: Ji Xuefei Inventor after: Li Ye Inventor after: Sun Qiang Inventor after: Zhang Yanqiu Inventor after: Li Kai Inventor before: Wang Jue Inventor before: Ding Rui Inventor before: Ji Xuefei Inventor before: Li Ye Inventor before: Sun Qiang Inventor before: Zhang Yanqiu Inventor before: Li Kai |
|
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200515 |