CN106411439B - Based on the ray emission method for weighting received simple match formula - Google Patents
Based on the ray emission method for weighting received simple match formula Download PDFInfo
- Publication number
- CN106411439B CN106411439B CN201611054816.8A CN201611054816A CN106411439B CN 106411439 B CN106411439 B CN 106411439B CN 201611054816 A CN201611054816 A CN 201611054816A CN 106411439 B CN106411439 B CN 106411439B
- Authority
- CN
- China
- Prior art keywords
- ray
- tubular body
- formula
- power
- antenna
- 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
-
- 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
- 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
-
- 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/3913—Predictive models, e.g. based on neural network models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The invention belongs to field intensity prediction fields to make improvements received field calculating section, field strength is avoided to compute repeatedly problem for the ray emission new method for proposing simple match formula, the final field intensity prediction precision for improving prediction model.Thus, the technical solution adopted by the present invention is that, based on the ray emission method for weighting received simple match formula, a multipath component is represented with a ray tubular body, power is uniformly distributed in tubular body, direct projection and the search of reflection path are completed with the path search algorithm of Ray Launching, it is wherein handled using received method is weighted near receiving antenna using the process for receiving ball screening path, reception specifically is weighted to power entrained by every effective rays, weight is to represent the tubular body of ray and receive ball great circle intersecting area to account for the ratio between tubular body cross-sectional area.Present invention is mainly applied to field intensity predictions.
Description
Technical field
The invention belongs to field intensity prediction fields, specifically, are related to based on the ray emission for weighting received simple match formula
Method.
Background technique
The ray shooting method (Ray Launching) of simple match formula is a kind of than more typical positive ray trace side
Method.RayLaunching method indicates path with ray, (sufficiently small) is equiangularly spaced from source point penetrates to the total space
Then line tracks the geometric path for calculating every ray and its corresponding field intensity information, finally with each effective at receiving antenna
The field strength vector sum in path receives field strength as total, and the method generally can in advance screen path using reception ball, if penetrating
To intersect then decision paths effective with ball is received for line.
Ray Launching has algorithm simple, scene strong applicability, receiving point number unrelated with route searching number etc. excellent
Point.But the field strength as caused by reception ball computes repeatedly problem and is worth noting in its receive process.Field strength compute repeatedly refer to by
The ray caused by the geometrical property for receiving ball is possible to the problem of being received twice.
Summary of the invention
In order to overcome the deficiencies of the prior art, the present invention is directed to propose the ray emission new method of simple match formula, to reception
Field computation part makes improvements, and field strength is avoided to compute repeatedly problem, the final field intensity prediction precision for improving prediction model.For this purpose,
The technical solution adopted by the present invention is that based on the ray emission method for weighting received simple match formula, with a ray tubulose
Body represents a multipath component, and power is uniformly distributed in tubular body, completes direct projection with the path search algorithm of Ray Launching
And the search of reflection path, wherein being handled using received method is weighted near receiving antenna using reception ball screening path
Process, reception specifically is weighted to power entrained by every effective rays, weight is the tubular body for representing ray
The ratio between tubular body cross-sectional area is accounted for ball great circle intersecting area is received.
Specific steps refine in an example are as follows:
A piece ray tubular body represents a multipath component, the horizontal and vertical angular spacing of every ray tubular bodyΔ
θ is identical, and ray tracing models are according to antenna transmission power PtAnd gain GtThe initial field strength of every root canal is directly set, l indicates current
Distance of the position to transmitting antenna, initial fields E0The maximum field intensity being defined as at range transmission antenna Tx central point l=1m
The vertical direction field strength E of amplitude, the neighboring area Tx meets formula (1),It is directions of rays and horizontal plane i.e. antenna method plane holder
Angle, E are unrelated with horizontal azimuth θ:
Formula (2) indicates that the isotropically radiated power near Tx, the π of η=120 Ω are free space impedance,For Poynting vector
Amount, initial fields E0With the relationship of the general power of the spherical surface by the inclusion of Tx:Individual ray field is obtained by formula (1)
It is strong:Enabling the ds in formula (2) is the cross section of tubular body,Then single is penetrated
Shown in the power of line such as formula (3):
It indicates that i-th ray directly reaches i.e. order of reflection n=0, or is reached near receiving antenna after n times reflect
Field strength, form meets formula (4), and j indicates the jth rank reflection in the reflection of this n times, n > 0, j=1,2 ... n, wherein Rj
For jth rank reflection coefficient, d is total path length,For i-th RX path direction of an electric field vector, ray is partially received
Situation, using intersecting area accounting method distribution power, as shown in formula (5), wherein Si,crossExpression represents the tubulose of ray
Body cross section and the intersecting area for receiving spherical surface, SiExpression represents the tubular body cross-sectional area of ray, Si,crossWith SiThe ratio between be
Weight has S if the tubular body for representing ray is received ball and completely includes i.e. complete receptioni,cross=Si, weight 1, Pi
Indicate i-th ray carry whole power, be withCorresponding power, PriIndicate i-th ray being received by the antenna
Effective power, in formula (6),For PriEffective reception field strength of corresponding i-th ray tube, formula (7) indicate final
The reception field strength acquired
The effective reception field strength for the i ray that will all intersect with reception ballIt adds up.
The features of the present invention and beneficial effect are:
Compared with existing Ray Launching method, the present invention has the following advantages and beneficial effects: improved
RayLaunching method is compared than traditional Ray Launching method, and field intensity prediction precision is higher, can more reflect local field strength
Variation tendency.
Detailed description of the invention:
Fig. 1 is the transmitting antenna signal of ray geometry distribution nearby；
Fig. 2 is common polygon trigonometric ratio flow chart；
Fig. 3 is indoor scene model and realistic picture；
Fig. 4 is the simulation result and measured result of conventional method, improved method.
Specific embodiment
For the above problem in existing prediction model, the improved model of the present invention is represented with a ray tubular body
One multipath component (power is uniformly distributed in tubular body) completes direct projection with the path search algorithm of traditional Ray Launching
And the search of reflection path, it thes improvement is that and introduces the received method of weighting to handle near receiving antenna using reception
The process in ball screening path.The method is to be weighted reception to power entrained by every effective rays, and weight is penetrated for representative
The tubular body and reception ball great circle intersecting area of line account for the ratio between tubular body cross-sectional area.
A multipath component is represented with a ray tubular body herein.Transmitting antenna Tx ray geometry distribution nearby signal is such as
Shown in Fig. 1, the horizontal and vertical angular spacing of every lineΔ θ is identical and smaller.Ray tracing models are generally sent out according to antenna
Penetrate power PtAnd gain GtThe initial field strength of every root canal is directly set, l indicates distance of the current location to transmitting antenna, initial fields
E0The maximum field intensity amplitude being defined as at range transmission antenna Tx central point l=1m.The vertical direction field of the neighboring area Tx
Strong E meets formula (1),It is directions of rays and horizontal plane (antenna normal plane) angle, E is unrelated with horizontal azimuth θ.
Formula (2) indicates that the isotropically radiated power near Tx, the π of η=120 Ω are free space impedance,For Poynting vector
(Poynting Vector).Initial fields E0With the relationship of the general power of the spherical surface by the inclusion of Tx:By formula
(1) single ray field strength is obtainedEnabling the ds in formula (2) is the cross section of tubular body,The then power of single rayAs shown in formula (3).
The present invention handles that may be present penetrate using the intersecting area for receiving ball and tubular body when calculating reception power
Line is by the received situation in part, it is contemplated that antenna body is to the influence factor of field, and ray is practical receives for more accurate simulation
Process.
Indicate that i-th ray directly reaches (order of reflection n=0) or reach near receiving antenna after n times reflect
Field strength, form meets formula (4), j indicate in the reflection of this n times the reflection of jth rank (n > 0, j=1,2 ... n).Wherein, R
For reflection coefficient, d is total path length,For i-th RX path direction of an electric field vector.Fig. 2 is the part for receiving ball and ray
Intersect situation and all intersection situation schematic diagram.For ray by the received situation in part, intersecting area accounting method is used herein
Distribution power, as shown in formula (5).Wherein, Si,crossIndicate the tubular body cross section for representing ray and the intersection for receiving spherical surface
Area, SiExpression represents the tubular body cross-sectional area of ray, Si,crossWith SiThe ratio between be weight, if representing the tubulose of ray
Body is received ball and completely includes and (receive completely), then has Si,cross=Si, weight 1.Pi(Corresponding power) indicate i-th
Whole power that ray carries, PriIndicate the effective power for i-th ray being received by the antenna.In formula (6),For
PriEffective reception field strength of corresponding i-th ray tube.Formula (7) indicates the reception field strength finally acquiredIt will be whole
With the effective reception field strength for receiving the i ray that ball intersectsIt adds up.
Below in conjunction with attached drawing and specific embodiment, the present invention is described further:
Below in conjunction with the specific example of two indoor scenes.The present invention in two and certain corridor inside certain room one
Item emulate and has been surveyed to the sample point being located on path under totally three RX paths, and to emulate and measured result into
Comparison gone to verify the accuracy of this method, has finally again compared this method and tradition Ray Launching method.
Improved Ray Launching method emulation is to utilize the software realization independently write.The software be
The visualization window application program write under 2010 platform of VisualStudio using C# language and OpenGL tool.The program
The graphic file that pre-rendered DXF format can be loaded, setting, dual-mode antenna, dielectric constant, reception radius of a ball etc. are each
After item parameter, route searching is carried out by test for intersection, field computation is reflected in direct projection field, and calculated result is stored, shows, lead
Result etc. out.Fig. 3 is the indoor scene illustraton of model and its realistic picture for carrying out emulation experiment.Wherein, the model of place in room shows
Meaning and realistic picture such as Fig. 3 a, shown in Fig. 3 b.Model of place signal and realistic picture such as Fig. 3 c in corridor, shown in Fig. 3 d.
The present invention respectively indoors with emulation experiment has been carried out to the point on three different paths under two scenes of corridor.Its
In, the launch point coordinate of scene is (14.5,1.85,2.35) in room, in the room path 1---y=1.85m, x=5.15
~9.95m, z=1.4m；Path 2---y=1.17m in room, x=6.05~11.45m uniformly take on z=1.4m a little and will
These points are predicted as receiving point by emulation；The launch point coordinate of corridor scene is (11.1,0.2,2), on corridor road
It is uniformly taken on diameter 3---y=0.9m, x=2~10m, z=1.4m a little and these points are pre- by emulation progress as receiving point
It surveys.
Under same two indoor scenes, the present invention three paths identical to launch point and routing information carry out real
It surveys.The launch point coordinate of scene is (14.5,1.85,2.35) in room, in the room path 1---y=1.85m, x=5.15
It is uniformly taken on path 2---y=1.17m in~9.95m, z=1.4m and room, x=6.05~11.45m, z=1.4m a little and will
These points are measured as receiving point；The launch point coordinate of corridor scene is (11.1,0.2,2), in corridor path 3---y=
It uniformly takes on 0.9m, x=2~10m, z=1.4m a little and is measured by these points as receiving point.
During actual measurement, transmitting terminal and receiving end are all made of vertical polarization single antenna, and transmitting 2.5GHz narrow band signal carries out
Actual measurement.Hardware use software radio USRP equipment (daughter board model RFX2400, maximum radio frequency power 50mW), this equipment it is soft
Part running environment is Linux system.Measured data is stored in database in real time.
In order to verify the superiority of improved method proposed by the present invention, the transmitting and reception of ray is respectively adopted in we
Traditional Ray Launching method and improved Ray Launching method are realized.The selection of remaining electrical parameter, field
The part such as scape Triangulation Algorithm, test for intersection algorithm, mirror field calculation method is consistent.Fig. 4 a, Fig. 4 b, Fig. 4 c are shown respectively
Field intensity prediction result and actual measured results that two methods obtain on above-mentioned three different RX paths are shown.
Upper table lists the root-mean-square error numerical value of traditional Ray Launching method and measured value and improves Ray
The root-mean-square error numerical value of Launching method and measured value.Improved method and the root-mean-square error value of measured value and biography
System method, which is compared, averagely reduces 3.97dB.It can be seen that improved Ray Launching method is than traditional Ray
Launching method precision of prediction is higher, can more reflect local field strength variation tendency.
Claims (2)
1. a kind of based on the ray emission method for weighting received simple match formula, characterized in that with a ray tubular body generation
One multipath component of table, power is uniformly distributed in tubular body, completes direct projection and anti-with the path search algorithm of Ray Launching
The search of rays diameter, wherein being handled using received method is weighted near receiving antenna using the mistake for receiving ball screening path
Journey is specifically weighted reception to power entrained by every effective rays, and weight is to represent the tubular body of ray and connect
It receives ball great circle intersecting area and accounts for the ratio between tubular body cross-sectional area.
2. as described in claim 1 based on the ray emission method for weighting received simple match formula, characterized in that at one
Specific steps refine in example are as follows:
A piece ray tubular body represents a multipath component, the horizontal and vertical angular spacing of every ray tubular bodyΔ θ phase
Together, ray tracing models are according to antenna transmission power PtAnd gain GtThe initial field strength of every root canal is directly set, l indicates present bit
Set the distance of transmitting antenna, initial fields E0The maximum field intensity width being defined as at range transmission antenna Tx central point l=1m
Value, the vertical direction field strength E of the neighboring area Tx meet formula (1),It is directions of rays and horizontal plane i.e. antenna method plane included angle,
E is unrelated with horizontal azimuth θ:
Formula (2) indicates that the isotropically radiated power near Tx, the π of η=120 Ω are free space impedance,For Poynting vector, just
Beginning field E0With the relationship of the general power of the spherical surface by the inclusion of Tx:Individual ray field strength is obtained by formula (1):Enabling the ds in formula (2) is the cross section of tubular body,Then individual ray
Shown in power such as formula (3):
It indicates that i-th ray directly reaches i.e. order of reflection n=0, or reaches the field near receiving antenna after n times reflect
By force, form meets formula (4), and j indicates that the jth rank in the reflection of this n times is reflected, n > 0, j=1,2 ... n, wherein
R is reflection coefficient, RjFor jth rank reflection coefficient, d is total path length,It is right for i-th RX path direction of an electric field vector
In ray by the received situation in part, using intersecting area accounting method distribution power, as shown in formula (5), wherein Si,crossTable
Show the tubular body cross section for representing ray and receives the intersecting area of spherical surface, SiExpression represents the tubular body cross-sectional area of ray,
Si,crossWith SiThe ratio between be weight have if the tubular body for representing ray is received ball and completely includes i.e. complete reception
Si,cross=Si, weight 1, PiIndicate i-th ray carry whole power, be withCorresponding power, PriIndicate quilt
The effective power for i-th ray that antenna receives, in formula (6),For effectively connecing for corresponding i-th ray tube of Pri
End is strong, and formula (7) indicates the reception field strength finally acquired
The effective reception field strength for the i ray that will all intersect with reception ballIt adds up.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611054816.8A CN106411439B (en) | 2016-11-25 | 2016-11-25 | Based on the ray emission method for weighting received simple match formula |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611054816.8A CN106411439B (en) | 2016-11-25 | 2016-11-25 | Based on the ray emission method for weighting received simple match formula |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106411439A CN106411439A (en) | 2017-02-15 |
| CN106411439B true CN106411439B (en) | 2019-04-05 |
Family
ID=58081535
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611054816.8A Active CN106411439B (en) | 2016-11-25 | 2016-11-25 | Based on the ray emission method for weighting received simple match formula |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106411439B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107276705B (en) * | 2017-07-28 | 2020-07-31 | 中铁上海设计院集团有限公司 | High-speed railway tunnel wireless communication channel modeling method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102111874A (en) * | 2009-12-24 | 2011-06-29 | 杨槐 | Cellular mobile communication positioning technology |
| CN105160698A (en) * | 2015-08-21 | 2015-12-16 | 天津大学 | Triangulation ray tracing path searching method |
-
2016
- 2016-11-25 CN CN201611054816.8A patent/CN106411439B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102111874A (en) * | 2009-12-24 | 2011-06-29 | 杨槐 | Cellular mobile communication positioning technology |
| CN105160698A (en) * | 2015-08-21 | 2015-12-16 | 天津大学 | Triangulation ray tracing path searching method |
Non-Patent Citations (3)
| Title |
|---|
| Improved 3D Ray Launching Method for Wireless Propagation Prediction;Durgin G,Patwari N等;《IEEE》;19970618;第1412-1413页 |
| 基于射线追踪法的3D无线信道研究与仿真;余鑫;《中国优秀硕士学位论文全文数据库科技信息辑》;20150430;全文 |
| 用SBR法研究室内微蜂窝的电波传播;张保华等;《南京邮电学院学报》;19980831;第18卷(第4期);第31-35页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106411439A (en) | 2017-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109239684B (en) | Radar target echo simulation system based on electromagnetic scattering model and simulation method thereof | |
| CN107238825B (en) | The test method of RCS when a kind of realization antenna transmitting using vector network instrument | |
| CN104619020A (en) | RSSI and TOA distance measurement based WIFI indoor positioning method | |
| CN105510876B (en) | A kind of indoor distance-measuring and positioning method based on propagation properties | |
| CN109188355B (en) | Method for optimizing and optimally distributing receiving antennas of multi-point positioning system | |
| CN106912105A (en) | 3-D positioning method based on PSO_BP neutral nets | |
| CN108802674A (en) | It is a kind of for the combined method for searching and device that directly position | |
| CN105467369B (en) | A kind of target echo simulation method and apparatus | |
| CN105005039A (en) | Satellite signal positioning method and system based on 3D modeling scene dynamic fingerprints | |
| CN110456655B (en) | Microwave landing simulation system | |
| CN109917333A (en) | Merge the passive location method of AOA observed quantity and TDOA observed quantity | |
| CN109151714A (en) | A kind of three-dimensional Robust Estimation localization method | |
| CN110068803A (en) | A kind of aerial bracketing device and method of radar equipment | |
| CN109901106A (en) | A kind of TDOA/AOA hybrid locating method | |
| CN110058204A (en) | One kind being based on the matched satellite antenna beam center calibrating method of directional diagram | |
| CN106842241A (en) | A kind of high dynamic satellite navigation receiver performance estimating method | |
| CN108732535A (en) | A kind of localization method, device and equipment | |
| CN109239686A (en) | A kind of transmitter and receiver layout method for the positioning of distributed MIMO radar target | |
| KR102319704B1 (en) | Measurement system and a method | |
| CN106411439B (en) | Based on the ray emission method for weighting received simple match formula | |
| CN105372676A (en) | Multi-path prediction method for three-dimensional scene navigation signal | |
| CN105160698B (en) | A kind of trigonometric ratio ray trace method for searching path | |
| CN113156418B (en) | Monte Carlo simulation-based radar target tracking precision prediction method | |
| CN110133588A (en) | A kind of antenna positioning method, device and equipment | |
| CN109579793A (en) | Mapping method, apparatus, flying platform, computer equipment and storage medium |
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 |