WO2011071020A1 - 伝搬経路推定方法、プログラム及び装置 - Google Patents
伝搬経路推定方法、プログラム及び装置 Download PDFInfo
- Publication number
- WO2011071020A1 WO2011071020A1 PCT/JP2010/071845 JP2010071845W WO2011071020A1 WO 2011071020 A1 WO2011071020 A1 WO 2011071020A1 JP 2010071845 W JP2010071845 W JP 2010071845W WO 2011071020 A1 WO2011071020 A1 WO 2011071020A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- point
- reflection
- propagation path
- calculating
- rotation center
- Prior art date
Links
Images
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
- 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/10—Monitoring; Testing of transmitters
- H04B17/101—Monitoring; Testing of transmitters for measurement of specific parameters of the transmitter or components thereof
- H04B17/104—Monitoring; Testing of transmitters for measurement of specific parameters of the transmitter or components thereof of other parameters, e.g. DC offset, delay or propagation times
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
- H04W16/18—Network planning tools
Definitions
- the present invention relates to a propagation path estimation method, program, and apparatus.
- the present invention relates to a means for analyzing a scatterer that does not perform specular reflection in an analysis using GO or GTD using a geometric optical model, a ray tracing analysis, or an analysis using an imaging method.
- Non-Patent Document 1 an analysis method using a geometric optical model is widely known as a method for analyzing the propagation characteristics of radio waves (see, for example, Non-Patent Document 1).
- Fig. 1 shows a method for obtaining a radio wave propagation path by a conventional ray tracing method.
- an image point T ′ with respect to the wall surface 1 of the transmission point T is obtained, and then an image point T ′′ with respect to the wall surface 2 of the image point T ′ is obtained. .
- T ′′ and the reception point R are connected by a straight line, and the intersection of the straight line and the wall surface 2 is obtained. This intersection is the reflection point X2 on the wall surface 2.
- the reflection point X1 on the wall surface 1 can be obtained by connecting the reflection points X2 and T 'on the wall surface 2 with a straight line.
- the propagation distance has the same value as the length of the straight line connecting T ′′ and the reception point R.
- the electric field strength at the reception point R can be obtained using this propagation distance.
- FIG. 2 shows a conventional technique for obtaining a propagation path for a transmitted wave and a reflected wave (see, for example, Non-Patent Document 2).
- Patent Document 1 discloses a method for estimating propagation characteristics for a case where the transmission direction is not limited to the straight traveling direction and the reflection direction is not limited to the specular reflection direction.
- the traveling direction of the reflected wave is limited to the specular reflection direction (regular reflection direction), and the traveling direction of the transmitted wave is the straight traveling direction. limited.
- the above-described method cannot analyze the propagation characteristics in the case of scattering in a direction other than the specular reflection direction and the propagation characteristics in consideration of the refraction of the medium.
- FIG. 3 shows an example of a propagation path when scattering in a direction other than the specular reflection direction.
- a wall surface 3 and a wall surface 4 are scatterers that do not reflect in the specular direction (for example, direction control scatterers such as a reflect array), and the reflected wave is scattered in a direction different from the specular reflection direction. Then, the radio wave incident from the direction of the incident direction A1 ° is scattered in the direction of the reflection direction A2 ° at the reflection point x3. On the wall surface 4, the radio wave incident from the direction of the incident direction B1 ° is reflected at the reflection point x4. Scattered in the direction of B2 °.
- the reflection point on the wall surface 4 is x2, which is a place different from the reflection point x4. That is, the reflection point on the wall surface 4 cannot be obtained by the conventional method shown in FIG. Therefore, there is a problem that the propagation path using a mirror image cannot be calculated by the conventional imaging method.
- Non-Patent Documents 3 and 4 examples of applying a reflectarray or metamaterial as a direction control scatterer for improving the propagation environment have been reported in recent years.
- Analyzing the propagation characteristics of a radio wave when such a direction-controlled scatterer is in the middle of a radio wave propagation path between the transmission point T and the reception point R is an analysis of the effect of improving the propagation environment. It is important to do this, but it has been difficult with the conventional ray tracing technique.
- An object of the present invention is to provide a propagation path estimation method, program, and apparatus that can be used.
- a first feature of the present invention is a propagation path estimation method using an imaging method, wherein a reflection direction and a scattering direction are a specular reflection direction ⁇ ° on a propagation path from a radio wave transmission point TX to a reception point RX.
- the virtual reception point VRX is calculated by rotating the reception point RX by ⁇ ° around the rotation center point O, and the virtual reception point VRX
- the gist is to have a step of estimating a propagation path using
- a second feature of the present invention is a propagation path estimation method using an imaging method, wherein a reflection direction and a scattering direction are a specular reflection direction ⁇ ° on a propagation path from a radio wave transmission point TX to a reception point RX.
- the virtual transmission point VTX is calculated by rotating the transmission point TX by ⁇ ° around the rotation center point O, and the virtual transmission point VTX
- the gist is to have a step of estimating a propagation path using
- a third feature of the present invention is a propagation path estimation method, wherein the radio wave has a direction ( ⁇ ) ° in which the reflection direction and the scattering direction are different from the specular reflection direction ⁇ ° on the propagation path. After reflection or scattering by the object, when the reflection direction and the scattering direction are reflected or scattered by the second structure having the specular reflection direction ⁇ °, the first transmission point with respect to the first structure using the imaging method is used.
- the present invention includes a step and a step of estimating a propagation path using the second image point and the virtual reception point.
- a fourth feature of the present invention is a propagation path estimation method, in which radio waves are reflected or scattered by the first structure whose reflection direction and scattering direction are the specular reflection direction ⁇ ° on the propagation path, and then reflected.
- radio waves are reflected or scattered by the first structure whose reflection direction and scattering direction are the specular reflection direction ⁇ ° on the propagation path, and then reflected.
- FIG. 1 is a diagram for explaining a method for obtaining a radio wave propagation path (path) by a conventional ray tracing method.
- FIG. 2 is a diagram showing a conventional method for obtaining propagation paths (paths) for transmitted waves and reflected waves.
- FIG. 3 is a diagram for explaining the problems of the conventional method.
- FIG. 4 is a diagram for explaining the propagation path estimation method according to the first embodiment of the present invention.
- FIG. 5 is a flowchart for explaining the propagation path estimation method according to the first embodiment of the present invention.
- FIG. 6 is a diagram for explaining a propagation path estimation method according to the second embodiment of the present invention.
- FIG. 7 is a diagram for explaining a propagation path estimation method according to the third embodiment of the present invention.
- FIG. 8 is a diagram for explaining a propagation path estimation method according to the third embodiment of the present invention.
- FIG. 9 is a flowchart for explaining a propagation path estimation method according to the third embodiment of the present invention.
- FIG. 10 is a diagram for explaining a propagation path estimation method according to the fourth embodiment of the present invention.
- FIG. 11 is a diagram for explaining a propagation path estimation method according to the fifth embodiment of the present invention.
- FIG. 12 is a flowchart for explaining a propagation path estimation method according to the fifth embodiment of the present invention.
- FIG. 13 is a diagram for explaining a propagation path estimation method according to the sixth embodiment of the present invention.
- FIG. 14 is a diagram for explaining a propagation path estimation method according to the sixth embodiment of the present invention.
- FIG. 14 is a diagram for explaining a propagation path estimation method according to the sixth embodiment of the present invention.
- FIG. 15 is a flowchart for explaining a propagation path estimation method according to the sixth embodiment of the present invention.
- FIG. 16 is a diagram for explaining a propagation path estimation method according to the seventh embodiment of the present invention.
- FIG. 17 is a diagram for explaining a propagation path estimation method according to the eighth embodiment of the present invention.
- FIG. 18 is a diagram for explaining a propagation path estimation method according to the ninth embodiment of the present invention.
- FIG. 19 is a diagram for explaining a propagation path estimation method according to the tenth embodiment of the present invention.
- FIG. 20 is a diagram (before the reception point is rotated) for explaining the propagation path estimation method according to the eleventh embodiment of the present invention.
- FIG. 21 is a diagram (first time / 45 ° rotation) for explaining the propagation path estimation method according to the eleventh embodiment of the present invention.
- FIG. 22 is a diagram for explaining the propagation path estimation method according to the eleventh embodiment of the present invention (second time / 45 ° rotation).
- FIG. 23 is a diagram (second time / 45 ° rotation) for explaining the propagation path estimation method according to the eleventh embodiment of the present invention.
- FIG. 24 is a diagram (third time / 45 ° rotation) for explaining the propagation path estimation method according to the eleventh embodiment of the present invention.
- FIG. 25 is a diagram for explaining the propagation path estimation method according to the eleventh embodiment of the present invention (third time / 45 ° rotation).
- FIG. 22 is a diagram for explaining the propagation path estimation method according to the eleventh embodiment of the present invention (second time / 45 ° rotation).
- FIG. 23 is a diagram (second time / 45 ° rotation) for explaining the propagation path estimation method according to the eleven
- FIG. 26 is a diagram showing a state of convergence of the distance to the reflection point by the propagation path estimation method according to the eleventh embodiment of the present invention.
- FIG. 27 is a diagram (before the reception point is rotated) for explaining the propagation path estimation method according to the twelfth embodiment of the present invention.
- FIG. 28 is a diagram for explaining a propagation path estimation method according to the twelfth embodiment of the present invention (first time / 70 ° rotation).
- FIG. 29 is a diagram for explaining the propagation path estimation method according to the twelfth embodiment of the present invention (first time / 70 ° rotation).
- FIG. 30 is a diagram (second time / 70 ° rotation) for explaining the propagation path estimation method according to the twelfth embodiment of the present invention.
- FIG. 27 is a diagram (before the reception point is rotated) for explaining the propagation path estimation method according to the twelfth embodiment of the present invention.
- FIG. 28 is a diagram for explaining a propag
- FIG. 31 is a diagram for explaining the propagation path estimation method according to the twelfth embodiment of the present invention (fifth rotation of 70 °).
- FIG. 32 is a diagram showing a state of convergence of the distance to the reflection point by the propagation path estimation method according to the twelfth embodiment of the present invention.
- FIG. 33 is a diagram (before the reception point is rotated) for explaining the propagation path estimation method according to the thirteenth embodiment of the present invention.
- FIG. 34 is a diagram (first time) illustrating a propagation path estimation method according to the thirteenth embodiment of the present invention.
- FIG. 35 is a diagram (first time) illustrating a propagation path estimation method according to the thirteenth embodiment of the present invention.
- FIG. 36 is a diagram (second time) for explaining the propagation path estimation method according to the thirteenth embodiment of the present invention.
- FIG. 37 is a diagram for explaining the propagation path estimation method according to the thirteenth embodiment of the present invention (after convergence, sixth time).
- FIG. 38 is a diagram showing a state of convergence of the distance to the reflection point by the propagation path estimation method according to the thirteenth embodiment of the present invention.
- FIG. 39 is a diagram (first time / 45 ° rotation) for explaining the propagation path estimation method according to the fourteenth embodiment of the present invention.
- FIG. 40 is a diagram for explaining the propagation path estimation method according to the fourteenth embodiment of the present invention (second time / 45 ° rotation).
- FIG. 41 is a diagram (third time / 45 ° rotation) for explaining the propagation path estimation method according to the fourteenth embodiment of the present invention.
- FIG. 42 is a diagram (fourth rotation / 45 ° rotation) for explaining the propagation path estimation method according to the fourteenth embodiment of the present invention.
- FIG. 43 is a diagram (fifth rotation / 45 ° rotation) for explaining the propagation path estimation method according to the fourteenth embodiment of the present invention.
- FIG. 44 is a diagram (8th rotation / 45 ° rotation) for explaining the propagation path estimation method according to the fourteenth embodiment of the present invention.
- FIG. 45 is a diagram showing a convergence state of reflection points by the propagation path estimation method according to the fourteenth embodiment of the present invention.
- FIG. 46 is a diagram showing a convergence state of reflection point differences by the propagation path estimation method according to the fourteenth embodiment of the present invention.
- FIG. 47 is a diagram (first time / 70 ° rotation) for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 48 is a diagram (second time / 70 ° rotation) for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 49 is a diagram (third time / 70 ° rotation) for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 50 is a diagram (fourth rotation / 70 ° rotation) for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 47 is a diagram (first time / 70 ° rotation) for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 48 is a diagram (second time / 70 ° rotation) for explaining the propagation path estimation method according to
- FIG. 51 is a diagram (fifth rotation / 70 ° rotation) for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 52 is a diagram (sixth rotation / 70 ° rotation) for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 53 is a diagram (seventh rotation / 70 ° rotation) for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 54 is a view for explaining the propagation path estimation method according to the fifteenth embodiment of the present invention (after convergence, rotation by 70 degrees).
- FIG. 55 is a diagram showing a convergence state of reflection points by the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 56 is a diagram showing a convergence state of reflection point differences by the propagation path estimation method according to the fifteenth embodiment of the present invention.
- FIG. 4 shows that the radio wave incident on the reflect array K from the incident direction ⁇ i ° rotates by ⁇ ° in the normal direction of the reflecting surface 1A of the reflect array K passing through the reflection point RP from the specular reflection direction ⁇ r ° with respect to the reflect array K.
- An example in the case of scattering in the above direction is shown.
- the incident direction is defined by the incident angle
- the reflection direction is defined by the reflection angle.
- ⁇ i ° ⁇ r °
- an arbitrary point on the reflection surface 1A of the reflect array K is set as the rotation center point O in step S101.
- step S102 the virtual reception point VRX is calculated by rotating the reception point RX around the rotation center point O by ⁇ ° in the normal direction of the reflection surface 1A of the reflect array K passing through the reflection point RP.
- step S103 an image point TX 'with respect to the reflective surface 1A (structure) of the reflect array K at the transmission point TX is calculated using an imaging method.
- step S104 the image point TX 'and the virtual reception point VRX are connected by a straight line, and the intersection of the straight line and the reflection surface 1A of the reflect array K is defined as a reflection point RP.
- step S105 a transmission path and a propagation distance are calculated by connecting the transmission point TX, the reflection point RP, and the reception point RX with a straight line.
- the reception point RX uses the rotation center point O. Since the center is deviated from the direction of ⁇ ° from the specular reflection direction (regular reflection direction), a calculation error occurs.
- the rotation center point O is moved a certain distance (for example, O ⁇
- ), and steps S101 to S104 (see FIG. 5) of the first propagation path estimation method are repeated.
- a propagation path estimation method according to the third embodiment of the present invention will be described with reference to FIGS.
- the propagation path estimation method according to the present embodiment will be described focusing on differences from the propagation path estimation method according to the first or second embodiment described above.
- step S201 an image point TX ′ with respect to the reflecting surface 1A (structure) of the reflect array K of the transmission point TX is calculated using an imaging method.
- step S202 the nth rotation center point On is calculated.
- the initial value of “n” is “1”.
- the first rotation center point O1 is an intersection of the normal of the reflecting surface 1A passing through the reception point RX and the reflecting surface 1A
- the nth rotation center point On is the (n ⁇ 1) th rotation center point On. This is a point where ⁇ 1 is changed by ⁇ s.
- step S203 the nth virtual reception point VRXn is calculated by rotating the reception point RX around the rotation center point On by ⁇ ° in the normal direction of the reflecting surface 1A of the reflect array K passing through the rotation center point On. To do.
- step S204 the image point TX 'and the nth virtual reception point VRXn are connected by a straight line, and the intersection of the straight line and the reflection surface 1A of the reflect array K is defined as an nth reflection point RPn.
- step S205 it is determined whether or not the convergence condition is satisfied.
- the convergence condition may be determined to be satisfied when the distance between the nth rotation center point On and the nth reflection point RPn is smaller than the threshold value ⁇ s.
- step S207 If it is determined that the convergence condition is satisfied, this operation proceeds to step S207. If it is determined that the convergence condition is not satisfied, this operation sets “n” to “1” in step S206. Is incremented by one and the process returns to step S202.
- step S207 a transmission path and a propagation distance are calculated by connecting the transmission point TX, the nth reflection point RPn, and the reception point RX with a straight line.
- a propagation path estimation method according to the fourth embodiment of the present invention will be described with reference to FIG.
- the propagation path estimation method according to the present embodiment will be described focusing on differences from the propagation path estimation methods according to the first to third embodiments described above.
- an orthogonal coordinate system is used in which the direction of the reflective surface 1A of the reflect array K is the X-axis direction and the direction orthogonal to the reflective surface 1A is the Y-axis direction.
- an arbitrary point on the reflecting surface 1A of the reflect array K is set as the rotation center point O in step S301.
- step S302 the virtual transmission point VTX is calculated by rotating the transmission point TX around the rotation center point O by ⁇ ° in the normal direction of the reflecting surface 1A of the reflect array K passing through the rotation center O.
- step S303 an image point VTX ′ with respect to the reflecting surface 1A of the reflect array K of the virtual transmission point VTX is calculated using an imaging method.
- step S304 the image point VTX 'and the reception point RX are connected by a straight line, and the intersection of the straight line and the reflection surface 1A of the reflect array K is defined as a reflection point RP.
- step S305 a transmission path and a propagation distance are calculated by connecting the transmission point TX, the reflection point RP, and the reception point RX with a straight line.
- the rotation center point O is set at a certain distance (for example, O) until the distance between the reflection point RP and the rotation center point O is within a predetermined distance (for example, ⁇ or ⁇ s). ⁇
- step S401 the nth rotation center point On is calculated.
- the initial value of “n” is “1”.
- the first rotation center point O1 is an intersection of the normal of the reflecting surface 1A passing through the reception point RX and the reflecting surface 1A
- the nth rotation center point On is the (n ⁇ 1) th rotation center point On. This is a point where ⁇ 1 is changed by ⁇ s.
- step S402 the transmission point TX is rotated about the rotation center point On by ⁇ ° in the normal direction of the reflecting surface 1A of the reflect array K passing through the rotation center point On, thereby calculating the nth virtual transmission point VTXn. To do.
- step S403 an image point VTX′n with respect to the reflecting surface 1A of the reflect array K of the nth virtual transmission point VTXn is calculated using an imaging method.
- step S404 the image point VTX'n and the reception point RX are connected by a straight line, and the intersection of the straight line and the reflection surface 1A of the reflect array K is defined as an nth reflection point RPn.
- step S405 it is determined whether or not the convergence condition is satisfied.
- the convergence condition may be determined to be satisfied when the distance between the nth rotation center point On and the nth reflection point RPn is smaller than the threshold value ⁇ or ⁇ s.
- step S407 If it is determined that the convergence condition is satisfied, this operation proceeds to step S407. If it is determined that the convergence condition is not satisfied, this operation sets “n” to “1” in step S406. Is incremented by one and the process returns to step S401.
- step S407 the transmission path and the propagation distance are calculated by connecting the transmission point TX, the nth reflection point RPn, and the reception point RX with a straight line.
- a propagation path estimation method according to the seventh embodiment of the present invention will be described with reference to FIG.
- the propagation path estimation method according to the present embodiment will be described focusing on differences from the propagation path estimation methods according to the first to sixth embodiments described above.
- an orthogonal coordinate system is used in which the direction of the reflective surface 1A of the reflect array K is the X-axis direction and the direction orthogonal to the reflective surface 1A is the Y-axis direction.
- a propagation path estimation method according to the eighth embodiment of the present invention will be described with reference to FIG.
- the propagation path estimation method according to the present embodiment will be described focusing on differences from the propagation path estimation methods according to the first to seventh embodiments described above.
- the propagation path estimation method has a structure in which radio waves have a direction ( ⁇ ) ° in which the reflection direction and the scattering direction are different from the specular reflection direction ⁇ ° on the propagation path. This is applied to a case where the light is reflected or scattered by K1 (first structure) and then reflected or scattered by the structure K2 (second structure) whose reflection direction and scattering direction are the specular reflection direction ⁇ °.
- the propagation path estimation method according to the present embodiment is performed in the following steps.
- the first image point TX ′ for the structure K1 at the transmission point TX is calculated using the imaging method, and the second image point TX ′′ for the structure K2 at the first image point TX ′ is calculated.
- the virtual reception point VRX is calculated by rotating the reception point RX around the rotation center point O in the direction opposite to the normal line of the reflection surface 1A of the reflect array K passing through the rotation center point O by ⁇ °.
- the second image point TX ′′ and the virtual reception point VRX are connected by a straight line, and the intersection of the straight line and the structure K2 is a reflection point RP2, and the first image point TX ′ and the reflection point RP2 are straight lines.
- the intersection between the straight line and the structure K1 is defined as a reflection point RP1.
- the propagation path and propagation distance are calculated by connecting the transmission point TX, the reflection point RP1, the reflection point RP2, and the reception point RX with a straight line.
- the rotation center point O is shifted until the convergence condition is satisfied.
- the above steps may be repeated.
- radio waves are reflected on the propagation path by a structure K1 (first structure) whose reflection direction and scattering direction are the specular reflection direction ⁇ °.
- a structure K1 first structure
- the reflection direction and the scattering direction are reflected or scattered by the structure K2 (second structure) whose direction is different from the specular reflection direction ⁇ ° ( ⁇ ) ° after scattering.
- the propagation path estimation method according to the present embodiment is performed in the following steps.
- the virtual transmission point VTX is calculated by rotating the transmission point TX around the rotation center point O by ⁇ ° in the normal direction of the reflecting surface 1A of the reflect array K passing through the rotation center point O.
- a first image point VTX ′ for the structure K1 of the virtual transmission point VTX is calculated using an imaging method, and a second image point VTX ′ ′ for the structure K2 of the first image point VTX ′ is calculated. .
- the second image point VTX ′′ and the reception point RX are connected by a straight line
- the intersection of the straight line and the structure K2 is set as a reflection point RP2
- the first image point VTX ′ and the reflection point RP2 are connected by a straight line.
- the intersection between the straight line and the structure K1 is defined as a reflection point RP1.
- the propagation path and propagation distance are calculated by connecting the transmission point TX, the reflection point RP1, the reflection point RP2, and the reception point RX with a straight line.
- the rotation center point O is shifted until the convergence condition is satisfied.
- the above steps may be repeated.
- the rotation angle ⁇ ° described above may be a negative value.
- the reception point RX is rotated by ⁇ ° in the direction opposite to the normal line of the reflection surface 1A of the reflect array K passing through the reflection point RP with the rotation center point O as the center.
- the virtual reception point VRX can be calculated.
- the rotation angle ⁇ ° described above may be a negative value. Further, ⁇ ° may be a function of the incident angle ⁇ °. That is, in the propagation path estimation method according to the fifth embodiment described above, the transmission point TX with the rotation center point O as the center is set to ⁇ in the direction opposite to the normal line of the reflection surface 1A of the reflect array K passing through the reflection point RP. By rotating, the virtual transmission point VTX can be calculated.
- the virtual reception point VRX, the virtual transmission point VTX, the rotation center point O, and the reflection point RP can be obtained for each propagation path.
- ⁇ ° is not limited to a fixed value, and may be a value that varies depending on the incident angle ⁇ i , for example.
- FIG. 20 is a diagram for explaining how to obtain the reflection point related to the reflective surface 1A of the reflect array, and represents two-dimensional coordinates with one end of the reflective surface 1A of the reflect array as the origin (0, 0).
- the coordinates of the transmission point Tx are (0, ⁇ 30)
- the coordinates of the reception point Rx are (20, ⁇ 30)
- an image point Tx ′ for the reflection surface 1A of the reflect array of the transmission point Tx is created.
- the coordinates of Tx 'are (0, 30)
- the coordinates of the intersection PR1 with the reflecting surface 1A of the reflectarray are (10, 0).
- the subsequent steps will be described with reference to FIG.
- the rotation angle ⁇ is 45 degrees.
- the coordinates of the first virtual reception point VRx1 are (32.28, -14.14). It can be confirmed that the distance from the rotation center O 1 to the reception point Rx is equal to the distance from the rotation center O 1 to VRx 1 and is 31.62.
- the subsequent steps will be described with reference to FIG.
- the intersection point between the straight line connecting the first virtual reflection point Vx1 and the image point Tx ′ of the transmission point and the reflection surface 1A of the reflect array is defined as a first reflection point PR1.
- the coordinates of PR1 are (26.02, 0).
- the subsequent steps will be described with reference to FIG.
- the coordinates of VRx2 are (42.98, -25.47).
- it obtains a second reflection point PR2, the second reflection point PR2 and third rotational center O 3 of the.
- PR2 and O 3 is (23.24,0).
- Request VRx3 by rotating 45 degrees Rx around the rotational center O 3.
- the convergence condition of the distance is ⁇ ⁇ 0.03.
- FIG. 26 shows the convergence state of the distance from the rotation point to the reflection point.
- the vertical axis in FIG. 26 is the distance between the nth reflection point PRn and the (n + 1) th reflection point PR (n + 1), and the horizontal axis is the number of trials. It can be seen that the coordinates of the reflection point rapidly converge by increasing the number of trials.
- FIG. 27 is a diagram for explaining how to obtain the reflection point related to the reflective surface 1A of the reflect array.
- the reflective surface 1A of the reflect array is set on the xOz plane and is represented by two-dimensional coordinates.
- the coordinates of the transmission point Tx are set to (0, 3)
- the coordinates of the reception point Rx are set to (5, 5)
- an image point Tx ′ for the reflection surface 1A of the reflect array of the transmission point Tx is created.
- the coordinates of Tx ′ are (0, ⁇ 3).
- FIG. 28 is a diagram illustrating a method of determining the first rotation point O 1 .
- FIG. 29 is a diagram illustrating steps for calculating the first virtual reception point VRx1 and the first reflection point PR1.
- the coordinates of VRx1 are (11.01, 3.59).
- the first reflection point PR1 is calculated by connecting VRx1 and Tx ′.
- the coordinates of PR1 are (5.01,0).
- FIG. 30 is a diagram illustrating steps for calculating the second rotation point O 2 , the virtual reception point VRx2, and the reflection point PR2. If the midpoint between O 1 and PR 1 is the rotation point O 2 , the coordinates of O 2 are (6.01, 0).
- the distance between PR2 and O 2 is 0.51, because they do not meet the greater convergence condition than epsilon, followed by the third rotation point O 3, virtual reception point VRx3, to calculate the reflection point PR3. Repeat these steps until the convergence condition is met. Similarly, the steps of calculating the coordinates of the third and fourth rotation points, virtual reception points, and reflection points are repeated. At this time, the respective coordinates are O 3 (5.75,0), VRx3 (10.19, 2.42), PR3 (5.65,0), O 4 (5.7,0), VRx4 (10 .16, 2.37) and PR4 (5.68,0). Further, the distance from O 3 to PR 3 and the distance from O 4 to PR 4 are 0.11 and 0.02, respectively (the diagram on the way is omitted).
- FIG. 31 is a diagram illustrating the coordinates of the rotation point, the virtual reception point, and the reflection point in a converged state.
- the convergence state corresponds to the completion of the fifth reflection point calculation, and is O 5 (5.69,0), VRx5 (10.15, 2.36), and PR5 (5.68,0).
- the distance between the PR5 and O 5 is 0.004, or less epsilon.
- the incident angle ⁇ i and the reflection angle ⁇ r are 62.18 ° and ⁇ 7.8 °, respectively.
- FIG. 32 is a diagram illustrating a state of convergence of the distance from the rotation point to the reflection point.
- the vertical axis is the distance between the reflection point PRn a rotation point O n the n-th, the horizontal axis represents the number of trials.
- the incident angle ⁇ i , the reflection angle ⁇ r , and the rotation angle ⁇ satisfy the following conditions.
- sin ( ⁇ r ) sin ( ⁇ i ) ⁇ sin ( ⁇ p ) (1)
- ⁇ r ⁇ i ⁇ n (2)
- the rotation angle ⁇ is a function of the incident angle ⁇ i .
- FIG. 33 is a diagram for explaining how to obtain the reflection point related to the reflective surface 1A of the reflect array.
- the reflective surface 1A of the reflect array is installed on the xOz plane and is represented by two-dimensional coordinates.
- the coordinates of the transmission point Tx are set to (0, 3)
- the coordinates of the reception point Rx are set to (5, 5)
- an image point Tx ′ for the reflection surface 1A of the reflect array of the transmission point Tx is created.
- the coordinates of Tx ′ are (0, ⁇ 3).
- FIG. 34 is a diagram illustrating a method for determining the first rotation point O 1 .
- the incident angle ⁇ i is 18.44 °.
- the reflection angle ⁇ r and the rotation angle ⁇ 1 can be obtained as ⁇ 33.57 ° and 57.00 °, respectively.
- the angle ⁇ is 51.34 °. Since the angle condition ( ⁇ ⁇ ⁇ n ) of the rotation point is not satisfied, it is necessary to determine the rotation point O 1 again.
- FIG. 35 is a diagram illustrating steps for calculating the first virtual reception point VRx1 and the first reflection point PR1.
- FIG. 36 is a diagram illustrating steps for calculating the second rotation point O 2 , virtual reception point VRx2, and reflection point PR2. If the middle point between the rotation point O 1 and the reflection point PR1 and rotation point O 2, the coordinates of O 2 is (5.93,0).
- VRx2 is a point obtained by rotating Rx clockwise by a rotation angle ⁇ 2 around O 2 . At this time, the coordinates of VRx2 are (10.11, 2.39). Further, when the reflection point PR2 is calculated by connecting VRx2 and Tx ′, (5.15, 0) is obtained as the coordinate of PR2.
- O 3 (5.54, 0), VRx3 (9.84, 2.6), PR3 (5.27, 0) O 4 (5.41, 0), VRx4 (9.75, 2.5), PR4 (5.32, 0) O 5 (5.36, 0), VRx5 (9.72, 2.47), PR5 (5.33, 0)
- the distances of O 3 PR3, O 4 PR4, and O 5 PR5 are 0.265, 0.088, and 0.029, respectively (the middle figure is omitted).
- FIG. 37 is a diagram showing the coordinates of the rotation point, the virtual reception point, and the reflection point in the convergence state.
- the convergence state corresponds to the completion of the sixth reflection calculation, and is O 6 (5.35, 0), VRx6 (9.71, 2.46), PR6 (5.33, 0). In this case, the distance between PR6 and O 6 0.01, or less epsilon.
- the incident angle ⁇ i and the reflection angle ⁇ r are 60.66 ° and ⁇ 3.86 °, respectively.
- FIG. 38 is a diagram illustrating a state of convergence of the distance from the rotation point to the reflection point.
- the vertical axis is the distance between the reflection point PRn a rotation point O n the n-th, the horizontal axis represents the number of trials.
- FIG. 39 is a diagram for explaining how to obtain a propagation path when there are two reflecting surfaces including the reflecting surface of the reflect array and the number of times of reflection is two, and one end of the reflecting surface 1A of the reflect array is defined as the origin. This represents a two-dimensional coordinate (0, 0). Specifically, in FIG. 39, the coordinates of the transmission point Tx are (0, ⁇ 100), the coordinates of the reception point Rx are (40, ⁇ 40), and one reflection is performed on each of the two reflection surfaces from the transmission point Tx. That is, the calculation result for the case of performing reflection twice in total is shown.
- the first reflecting surface is the reflecting surface 1A of the reflect array
- the second reflecting surface is a wall surface that is normally regularly reflected.
- An image point Tx ′ for the reflection surface 1B of the transmission point Tx and an image point Tx ′′ for the reflection surface 1A of Tx ′ are created.
- the coordinates of Tx ′′ are ( ⁇ 60, ⁇ 100).
- the first rotation center O 1 is arbitrarily determined on the reflecting surface 1A.
- the first rotation center O 1 is (20, 0).
- a virtual reception point VRx1 is calculated by rotating the reception point Rx around the first rotation center O1.
- the rotation angle ⁇ is 45 degrees.
- the coordinates of the first virtual reception point VRx1 are (62.43, -14.14). It can be confirmed that the distance from the rotation center O 1 to the reception point Rx is equal to the distance from the rotation center O 1 to VRx 1 and is 44.72.
- An intersection point between a straight line connecting the first virtual reflection point Vx1 and the two-reflection image point Tx ′′ of the transmission point and the reflection surface 1A of the reflect array is defined as a first reflection point PR1.
- the coordinates of PR1 are (47.26, 0).
- the intersection point with the second reflection surface 1B is obtained as the intersection point between the reflection surface 1B and a straight line connecting the first reflection point PR1 and the image point Tx ′ of the transmission point Tx with respect to the reflection surface 1B.
- the coordinates of the reflection point are ( ⁇ 30, ⁇ 72.03).
- the coordinates of VRx2 are (70.41, ⁇ 33.42).
- it obtains a second reflection point PR2
- the second reflection point PR2 and third rotational center O 3 of the. PR2 and O 3 is (37.75,0).
- the virtual reception point VRx3 by rotating 45 degrees Rx around the rotational center O 3.
- FIGS. 42 and 43 by repeating the same steps, the coordinates of the reflection point PRn converge.
- FIG. 44 shows the situation after the eighth rotation.
- the incident angle of the reflect array on the reflection surface 1A is 45.00 degrees
- the reflection angle from the reflection point PR8 to the reflection point RX is 0 degrees with respect to the normal direction of the reflection surface. This coincides with the angle when 45 degree direction control is performed in the counterclockwise direction with respect to the original normal reflection direction of 45 degrees and converges to the result of calculating the propagation path to be finally obtained. It shows that you are doing.
- the vertical axis in FIG. 45 is the X coordinate of the reflective surface on the reflectarray reflective surface 1A
- the vertical axis in FIG. 46 is the distance between the nth reflection point PRn and the (n + 1) th reflection point PR (n + 1).
- the horizontal axis is the number of trials n. It can be seen that by increasing the number of trials n, the coordinates of the reflection points converge rapidly.
- FIG. 47 is a diagram for explaining a method of obtaining a propagation path when there are two reflecting surfaces including a reflecting surface of the reflect array and the number of reflections is two, and one end of the reflecting surface 1A of the reflect array is defined as an origin.
- This represents a two-dimensional coordinate (0, 0).
- the coordinates of the transmission point Tx are (0, -100), the coordinates of the reception point Rx are (30, -40), and one reflection is performed on the two reflecting surfaces from the transmission point Tx, that is, a total of two reflections are performed.
- the reflection rotation angle ⁇ is set to 70 degrees.
- the first reflecting surface is a reflecting surface 1A of the reflect array
- the second reflecting surface is a wall surface that is normally regularly reflected.
- An image point Tx ′ for the reflection surface 1B of the transmission point Tx and an image point Tx ′′ for the reflection surface 1A of Tx ′ are created.
- the coordinates of Tx ′′ are ( ⁇ 60, ⁇ 100).
- the first rotation center O 1 is arbitrarily determined on the reflecting surface 1A.
- the first rotation center O 1 is (20, 0).
- the virtual reception point VRx1 is calculated by rotating the reception point Rx around the first rotation center O1.
- the coordinates of the first virtual reception point VRx1 are (61.01, -4.28). It can be confirmed that the distance from the rotation center O 1 to the reception point Rx is equal to the distance from the rotation center O 1 to VRx 1 and is 41.23.
- the intersection point between the first virtual reflection point Vx1 and the two-time reflection image point Tx ′′ of the transmission point and the reflection surface 1A of the reflect array is defined as a first reflection point PR1.
- the coordinates of PR1 are (44.76, 0).
- the intersection point with the second reflection surface 1B is obtained as the intersection point between the reflection surface 1B and a straight line connecting the first reflection point PR1 and the image point Tx ′ of the transmission point Tx with respect to the reflection surface 1B.
- the subsequent steps will be described with reference to FIG.
- FIG. 54 shows the situation after the eighth rotation.
- the rotation angle from the virtual reception point VRX to the original reception point RX coincides with 70.00 degrees clockwise, and converges to the result of calculating the propagation path to be finally obtained. It shows that you are doing.
- the vertical axis in FIG. 55 is the x coordinate of the reflection point on the reflect array reflecting surface 1A in the n-th trial
- the vertical axis in FIG. 56 is the n-th reflection point PRn and the (n + 1) -th reflection point PR (n + 1).
- the horizontal axis is the number of trials n. It can be seen that the coordinates of the reflection point rapidly converge by increasing the number of trials.
- a first feature of the present embodiment is a propagation path estimation method using an imaging method, and is on a propagation path from a radio wave transmission point TX to a reception point RX.
- the reception point RX is centered on the rotation center point O.
- the gist is to calculate a virtual reception point VRX by rotating by ⁇ ° and to estimate a propagation path using the virtual reception point VRX.
- the above-described process calculates an image point TX ′ with respect to the reflection surface 1A (structure) of the transmission point TX using the imaging method.
- the above-described steps include a step A of calculating an image point TX ′ with respect to the reflection surface 1A of the transmission point TX using an imaging method, Step B for calculating the rotation center point O, Step C for calculating the virtual reception point VRX by rotating the reception point RX around the rotation center point O by ⁇ °, and the image point TX ′ and the virtual reception point VRX
- the step D for calculating the reflection point RP on the reflection surface 1A of the radio wave is calculated using “O ⁇
- the above-described steps include a step A of calculating an image point TX ′ with respect to the reflecting surface 1A of the transmission point TX using an imaging method, Step B for calculating the rotation center point O, Step C for calculating the virtual reception point VRX by rotating the reception point RX around the rotation center point O by ⁇ °, and the image point TX ′ and the virtual reception point VRX
- the step D for calculating the reflection point RP on the radio wave reflection surface 1A is defined as “O ⁇ s”.
- the process E to be updated, and the processes B to E may be repeated until the distance between the reflection point RP and the rotation center point O satisfies the convergence condition.
- the above-described steps are orthogonal in which the direction of the reflective surface 1A is the X-axis direction and the direction orthogonal to the reflective surface 1A is the Y-axis direction.
- the coordinates of the transmission point TX are (a, b)
- the coordinates of the image point TX ′ with respect to the reflection surface of the transmission point TX are (a, ⁇ b)
- the coordinates of the reception point RX are (c, d).
- a second feature (features according to the fifth to seventh embodiments) of the present embodiment is a propagation path estimation method using an imaging method, and is on a propagation path from a radio wave transmission point TX to a reception point RX.
- the transmission point TX is rotated by ⁇ ° around the rotation center point O.
- the gist is to calculate the virtual transmission point VTX by using the virtual transmission point VTX and estimate the propagation path using the virtual transmission point VTX.
- the above-described steps include a step A for calculating the rotation center point O, and a transmission point TX around the rotation center point O by ⁇ °.
- a process B for calculating a virtual transmission point VTX and a process AC for calculating an image point VTX ′ for the reflection surface 1A of the virtual transmission point VTX using an imaging method, an image point VTX ′ and a reception point RX
- the above-described steps include a step A for calculating the rotation center point O, and a transmission point TX around the rotation center point O by ⁇ °.
- a process B for calculating a virtual transmission point VTX and a process AC for calculating an image point VTX ′ for the reflection surface 1A of the virtual transmission point VTX using an imaging method, an image point VTX ′ and a reception point RX
- the rotation center point O is set to “O ⁇
- the above-described steps include a step A for calculating the rotation center point O, and a transmission point TX around the rotation center point O by ⁇ °.
- a process B for calculating a virtual transmission point VTX by rotation a process C for calculating an image point VTX ′ for the reflection surface 1A of the virtual transmission point VTX using an imaging method, an image point VTX ′ and a reception point RX
- the step D for calculating the reflection point RP on the radio wave reflection surface 1A is defined as “O ⁇ s”.
- the process E to be updated, and the processes A to E may be repeated until the distance between the reflection point RP and the rotation center point O satisfies the convergence condition.
- the above-described steps are orthogonal in which the direction of the reflective surface 1A is the X-axis direction and the direction orthogonal to the reflective surface 1A is the Y-axis direction.
- a third feature of the present embodiment is a propagation path estimation method in which a radio wave reflects and reflects in a direction different from the specular reflection direction ⁇ ° on the propagation path.
- a radio wave After being reflected or scattered by the structure K1 (first structure) having ( ⁇ ) °, reflected or scattered by the structure K2 (second structure) whose reflection direction and scattering direction are the specular reflection direction ⁇ °.
- a first image point TX ′ for the structure K1 at the transmission point TX is calculated using an imaging method
- a second image point TX ′′ for the structure K2 at the first image point TX ′ is calculated.
- calculating the virtual reception point VRX by rotating the reception point RX around the rotation center point O by ⁇ °, and estimating the propagation path using the second image point TX ′′ and the virtual reception point VRX. And having a process.
- the fourth feature of the present embodiment is a propagation path estimation method, in which a radio wave reflects and reflects in a specular reflection direction ⁇ ° on the propagation path.
- a radio wave reflects and reflects in a specular reflection direction ⁇ ° on the propagation path.
- the rotation center point O is the center.
- the first image point TX ′ for the structure K1 of the virtual transmission point VTX is calculated, and the first image point The gist is to include a step of calculating a second image point TX ′′ for the structure K2 of TX ′ and a step of estimating a propagation path using the second image point TX ′′ and the reception point RX.
- the gist of the fifth feature of the present embodiment is that it is a program or apparatus for causing a computer to realize the above-described propagation path estimation method.
- the above-described propagation path estimation method may be implemented by hardware, may be implemented by a software module executed by a processor, or may be implemented by a combination of both.
- the software modules include RAM (Random Access Memory), flash memory, ROM (Read Only Memory), ERPOM (Erasable Rotable ROM), EERPOM (Electronically Erasable, Random Disk, RD ROM Alternatively, it may be provided in a storage medium of an arbitrary format such as a CD-ROM.
- Such a storage medium is connected to the processor so that the processor can read and write information from and to the storage medium. Further, such a storage medium may be integrated in the processor. Such a storage medium and processor may be provided in the ASIC.
- a propagation path estimation method, program, and apparatus capable of estimating a propagation path in a propagation analysis model including a direction control scatterer and performing an analysis applying a ray trace analysis or a geometric optical model. This is useful in wireless communication and the like.
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Optical Communication System (AREA)
Abstract
Description
図4及び図5を参照して、本発明の第1の実施形態に係る伝搬経路推定方法について説明する。
次に、図6を参照して、本発明の第2の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
図7乃至図9を参照して、本発明の第3の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1又は2の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
図10を参照して、本発明の第4の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1乃至3の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
- 送信点TXの反射面1Aに対するイメージ点TX’の座標:(a,-b)
- 受信点RXの座標:(c,d)
- 回転中心点Oの座標:(e,0)
- 仮想受信点VRXの座標:(Xvrx,Yvrx)=(cos(-η°)×(c-e)-sin(-η°)×d+e,sin(-η°)×(c-e)+cos(-η°)×d)
(x1,y1)及び(x2,y2)を通る直線の方程式は、
で表されるから、イメージ点TX’と仮想受信点VRXとを通る直線の方程式は、「y-(-b)=((sin(-η°)×(c-e)+cos(-η°)×d-a)/(cos(-η°)×(c-e)-sin(-η°)×d-(-b))+e)×(x-a)」で表される。
図11及び図12を参照して、本発明の第5の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1乃至4の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
図13乃至図15を参照して、本発明の第6の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1乃至5の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
図16を参照して、本発明の第7の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1乃至6の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
- 受信点RXの座標:(c,d)
- 回転中心点Oの座標:(e,0)
- 仮想送信点VTX:(va,vb)=(cos(η°)×(a-e)-sin(η°)×b+e,sin(η°)×(a-e)+cos(η°)×b)
- 仮想送信点VTXの反射面1Aに対するイメージ点VTX’の座標:(va,-vb)
したがって、イメージ点VTX’と受信点VRXとを通る直線の方程式は、「y-(-b)=((sin(η°)×(a-e)+cos(η°)×b-va)/(cos(η°)×(a-e)-sin(-η°)×(d-(-vb))+e)×(x-va)」で表される。
図17を参照して、本発明の第8の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1乃至7の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
第3に、第2イメージ点TX’’と仮想受信点VRXとを直線で結び、かかる直線と構造物K2との交点を反射点RP2とし、第1イメージ点TX’及び反射点RP2を直線で結び、かかる直線と構造物K1との交点を反射点RP1とする。
図18を参照して、本発明の第9の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1乃至8の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
図19を参照して、本発明の第9の実施形態に係る伝搬経路推定方法について説明する。以下、本実施形態に係る伝搬経路推定方法について、上述の第1の実施形態に係る伝搬経路推定方法との相違点に着目して説明する。
図20乃至図26を参照して、本発明の第11の実施形態に係る伝搬路推定方法について説明する。本実施形態では、伝搬路推定方法の具体的な数値計算例について説明する。
図23を用いて、引き続くステップについて説明する。1回目の反射点PR1を2回目の回転中心O2とする。2回目の回転中心O2を中心に受信点Rxをη(=45度)回転させてVRx2を求める。このときVRx2の座標は(42.98,-25.47)である。同様のステップにより、2回目の反射点PR2を求め、2回目の反射点PR2を3回目の回転中心O3とする。それぞれの座標は、図24及び図25に示す通り、PR2及びO3は(23.24,0)となる。この回転中心O3を中心としてRxを45度回転することによりVRx3を求める。このように同様のステップを繰り返すことにより、反射点PRnの座標は収束する。なお、本実施形態では、当該距離の収束条件は、ε<0.03とした。
図27乃至図32を参照して、本発明の第12の実施形態に係る伝搬路推定方法について説明する。本実施形態では、伝搬路推定方法の具体的な数値計算例について説明する。
図32は、回転点から反射点までの距離の収束状況を示す図である。縦軸はn回目の回転点Onと反射点PRnとの距離であり、横軸は試行回数である。試行回数を増やすことにより,急速に回転点Onと反射点PRnとの距離が小さくなり、収束する様子が分かる。
図33乃至図38を参照して、本発明の第13の実施形態に係る伝搬路推定方法について説明する。本実施形態では、伝搬路推定方法の具体的な数値計算例について説明する。
sin(θr)=sin(θi)-sin(θp) …(1)
θr=θi-ηn …(2)
ここで、θpを固定の角度70°とすると、回転角ηは入射角θiの関数となる。収束条件としては、n回目の回転点Onと反射点PRnの距離が、ε(=0.01)より小さい場合に収束したものとみなして計算を完了するものとする。
O3(5.54、0),VRx3(9.84,2.6),PR3(5.27,0)
O4(5.41,0),VRx4(9.75,2.5),PR4(5.32,0)
O5(5.36,0),VRx5(9.72,2.47),PR5(5.33,0)
また、O3PR3、O4PR4、O5PR5の距離はそれぞれ、0.265、0.088、0.029である(途中の図は省略する)。
図39乃至図46を参照して、本発明の第14の実施形態に係る伝搬路推定方法について説明する。本実施形態では、伝搬路推定方法の具体的な数値計算例について説明する。
図47乃至図56を参照して、本発明の第14の実施形態に係る伝搬路推定方法について説明する。本実施形態では、伝搬路推定方法の具体的な数値計算例について説明する。
1A…反射面
TX…送信点
VTX…仮想送信点
RX…受信点
VRX…仮想受信点
Claims (15)
- イメージング法を用いた伝搬経路推定方法であって、
電波の送信点TXから受信点RXまでの伝搬経路上に、反射方向及び散乱方向が鏡面反射方向θ°と異なる方向(θ-η)°となる構造物が存在する場合に、回転中心点Oを中心にして該受信点RXをη°回転させることによって仮想受信点VRXを算出し、該仮想受信点VRXを用いて伝搬経路を推定する工程を有することを特徴とする伝搬経路推定方法。 - 前記工程は、
前記イメージング法を用いて、前記送信点TXの前記構造物に対するイメージ点TX’を算出する工程Aと、
前記回転中心点Oを算出する工程Bと、
前記回転中心点Oを中心にして前記受信点RXをη°回転することによって、前記仮想受信点VRXを算出する工程Cと、
前記イメージ点TX’及び前記仮想受信点VRXを用いて、前記電波の前記構造物における反射点RPを算出する工程Dとを有することを特徴とする請求項1に記載の伝搬経路推定方法。 - 前記工程は、
前記イメージング法を用いて、前記送信点TXの前記構造物に対するイメージ点TX’を算出する工程Aと、
前記回転中心点Oを算出する工程Bと、
前記回転中心点Oを中心にして前記受信点RXをη°回転することによって、前記仮想受信点VRXを算出する工程Cと、
前記イメージ点TX’及び前記仮想受信点VRXを用いて、前記電波の前記構造物における反射点RPを算出する工程Dと、
前記反射点RPと前記回転中心点Oとの間の距離が収束条件を満たさないとき、前記回転中心点Oを「O-|RP-O|」と更新する工程Eとを有し、
前記反射点RPと前記回転中心点Oとの間の距離が収束条件を満たすまで、前記工程B乃至Eを繰り返すことを特徴とする請求項1に記載の伝搬経路推定方法。 - 前記工程は、
前記イメージング法を用いて、前記送信点TXの前記構造物に対するイメージ点TX’を算出する工程Aと、
前記回転中心点Oを算出する工程Bと、
前記回転中心点Oを中心にして前記受信点RXをη°回転することによって、前記仮想受信点VRXを算出する工程Cと、
前記イメージ点TX’及び前記仮想受信点VRXを用いて、前記電波の前記構造物における反射点RPを算出する工程Dと、
前記反射点RPと前記回転中心点Oとの間の距離が収束条件を満たさないとき、前記回転中心点Oを「O-Δs」と更新する工程Eとを有し、
前記反射点RPと前記回転中心点Oとの間の距離が収束条件を満たすまで、前記工程B乃至Eを繰り返すことを特徴とする請求項1に記載の伝搬経路推定方法。 - 前記工程は、
前記構造物の反射面の方向をX軸方向とし、該反射面に直交する方向をY軸方向とする直交座標系において、前記送信点TXの座標を(a,b)とし、前記送信点TXの前記構造物に対するイメージ点TX’の座標を(a,-b)とし、前記受信点RXの座標を(c,d)とし、前記回転中心点O及び前記電波の前記構造物における反射点RPの座標を(x,0)とし、式「x=b/{((sin(-η°)×(c-x)+cos(-η°)×d-a)/(cos(-η°)×(c-x)-sin(-η°)×d-(-b))+x)}+a」から、xの値を算出する工程を有することを特徴とする請求項1に記載の伝搬経路推定方法。 - イメージング法を用いた伝搬経路推定方法であって、
電波の送信点TXから受信点RXまでの伝搬経路上に、反射方向及び散乱方向が鏡面反射方向θ°と異なる方向(θ-η)°となる構造物が存在する場合に、回転中心点Oを中心にして該送信点TXをη°回転させることによって仮想送信点VTXを算出し、該仮想送信点VTXを用いて伝搬経路を推定する工程を有することを特徴とする伝搬経路推定方法。 - 前記工程は、
前記回転中心点Oを算出する工程Aと、
前記回転中心点Oを中心にして前記送信点TXをη°回転することによって、前記仮想送信点VTXを算出する工程Bと、
前記イメージング法を用いて、前記仮想送信点VTXの前記構造物に対するイメージ点VTX’を算出する工程AC、
前記イメージ点VTX’及び前記受信点RXを用いて、前記電波の前記構造物における反射点RPを算出する工程Dとを有することを特徴とする請求項6に記載の伝搬経路推定方法。 - 前記工程は、
前記回転中心点Oを算出する工程Aと、
前記回転中心点Oを中心にして前記送信点TXをη°回転することによって、前記仮想送信点VTXを算出する工程Bと、
前記イメージング法を用いて、前記仮想送信点VTXの前記構造物に対するイメージ点VTX’を算出する工程AC、
前記イメージ点VTX’及び前記受信点RXを用いて、前記電波の前記構造物における反射点RPを算出する工程Dと、
前記反射点RPと前記回転中心点Oとの間の距離が収束条件を満たさないとき、前記回転中心点Oを「O-|RP-O|」と更新する工程Eとを有し、
前記反射点RPと前記回転中心点Oとの間の距離が収束条件を満たすまで、前記工程A乃至Eを繰り返すことを特徴とする請求項6に記載の伝搬経路推定方法。 - 前記工程は、
前記回転中心点Oを算出する工程Aと、
前記回転中心点Oを中心にして前記送信点TXをη°回転することによって、前記仮想送信点VTXを算出する工程Bと、
前記イメージング法を用いて、前記仮想送信点VTXの前記構造物に対するイメージ点VTX’を算出する工程Cと、
前記イメージ点VTX’及び前記受信点RXを用いて、前記電波の前記構造物における反射点RPを算出する工程Dと、
前記反射点RPと前記回転中心点Oとの間の距離が収束条件を満たさないとき、前記回転中心点Oを「O-Δs」と更新する工程Eとを有し、
前記反射点RPと前記回転中心点Oとの間の距離が収束条件を満たすまで、前記工程A乃至Eを繰り返すことを特徴とする請求項6に記載の伝搬経路推定方法。 - 前記工程は、
前記構造物の反射面の方向をX軸方向とし、該反射面に直交する方向をY軸方向とする直交座標系において、前記送信点TXの座標を(a,b)とし、前記受信点RXの座標を(c,d)とし、前記回転中心点O及び前記電波の前記構造物における反射点RPの座標を(x,0)とし、前記仮想送信点VTXの座標を(va(=cos(η°)×(a-x)-sin(η°)×b+x),vb(=sin(η°)×(a-x)+cos(η°)×b))とし、式「vb=((sin(-η°)×c+cos(-η°)×d-va/(cos(-η°)×c-sin(-η°)×d-(-vb))+x)×(x-va))から、xの値を算出する工程を有することを特徴とする請求項6に記載の伝搬経路推定方法。 - 伝搬経路上で、電波が、反射方向及び散乱方向が鏡面反射方向θ°と異なる方向(θ-η)°となる第1構造物で反射又は散乱した後、反射方向及び散乱方向が鏡面反射方向θ°となる第2構造物で反射又は散乱する場合に、
イメージング法を用いて、送信点の前記第1構造物に対する第1イメージ点を算出し、該第1イメージ点の前記第2構造物に対する第2イメージ点を算出する工程と、
回転中心点を中心にして受信点をη°回転させることによって仮想受信点を算出する工程と、
前記第2イメージ点及び前記仮想受信点を用いて伝搬経路を推定する工程とを有することを特徴とする伝搬経路推定方法。 - 伝搬経路上で、電波が、反射方向及び散乱方向が鏡面反射方向θ°となる第1構造物で反射又は散乱した後、反射方向及び散乱方向が鏡面反射方向θ°と異なる方向(θ-η)°となる第2構造物で反射又は散乱する場合に、
回転中心点を中心にして送信点をη°回転することによって、仮想送信点を算出する工程と、
イメージング法を用いて、前記仮想送信点の前記第1構造物に対する第1イメージ点を算出し、該第1イメージ点の前記第2構造物に対する第2イメージ点を算出する工程と、
前記第2イメージ点及び受信点を用いて伝搬経路を推定する工程とを有することを特徴とする伝搬経路推定方法。 - 前記η°は、前記θ°の関数であることを特徴とする請求項1乃至12のいずれか一項に記載の伝搬経路推定方法。
- コンピュータに、請求項1乃至13のいずれか一項に記載の伝搬経路推定方法を実現させるためのプログラム。
- コンピュータに、請求項1乃至13のいずれか一項に記載の伝搬経路推定方法を実現させるための装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/514,426 US9337942B2 (en) | 2009-12-07 | 2010-12-06 | Propagation path estimation method and program and apparatus using the same |
EP10835940.7A EP2511718B1 (en) | 2009-12-07 | 2010-12-06 | Propagation path estimation method, program, and device |
CN201080055588.0A CN102648419B (zh) | 2009-12-07 | 2010-12-06 | 传输路径推定方法 |
JP2011545204A JP5687206B2 (ja) | 2009-12-07 | 2010-12-06 | 伝搬経路推定方法、プログラム及び装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009278045 | 2009-12-07 | ||
JP2009-278045 | 2009-12-07 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2011071020A1 true WO2011071020A1 (ja) | 2011-06-16 |
Family
ID=44145565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2010/071845 WO2011071020A1 (ja) | 2009-12-07 | 2010-12-06 | 伝搬経路推定方法、プログラム及び装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9337942B2 (ja) |
EP (1) | EP2511718B1 (ja) |
JP (1) | JP5687206B2 (ja) |
CN (1) | CN102648419B (ja) |
WO (1) | WO2011071020A1 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10638479B2 (en) | 2015-11-17 | 2020-04-28 | Futurewei Technologies, Inc. | System and method for multi-source channel estimation |
WO2021026409A1 (en) | 2019-08-08 | 2021-02-11 | Regeneron Pharmaceuticals, Inc. | Novel antigen binding molecule formats |
WO2021091953A1 (en) | 2019-11-05 | 2021-05-14 | Regeneron Pharmaceuticals, Inc. | N-terminal scfv multispecific binding molecules |
WO2021178896A1 (en) | 2020-03-06 | 2021-09-10 | Go Therapeutics, Inc. | Anti-glyco-cd44 antibodies and their uses |
WO2023014863A1 (en) | 2021-08-05 | 2023-02-09 | Go Therapeutics, Inc. | Anti-glyco-muc4 antibodies and their uses |
EP4233909A2 (en) | 2016-05-17 | 2023-08-30 | AbbVie Biotherapeutics Inc. | Anti-cmet antibody drug conjugates and methods for their use |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102640009B (zh) | 2009-12-04 | 2016-03-09 | 株式会社Ntt都科摩 | 传输路径推定方法 |
US9654306B1 (en) * | 2015-11-17 | 2017-05-16 | Futurewei Technologies, Inc. | System and method for multi-source channel estimation |
WO2023184255A1 (en) * | 2022-03-30 | 2023-10-05 | Huawei Technologies Co., Ltd. | Methods and systems for sensing-based channel reconstruction and tracking |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH088846A (ja) * | 1994-06-21 | 1996-01-12 | Hitachi Ltd | 電波受信強度シミュレーション方法 |
JP2003318842A (ja) * | 2002-04-26 | 2003-11-07 | Nippon Telegr & Teleph Corp <Ntt> | 受信電界強度推定計算装置及び方法並びにプログラム及び記録媒体 |
JP2009168534A (ja) | 2008-01-11 | 2009-07-30 | Ntt Docomo Inc | 伝搬特性推定方法及びプログラム |
JP2009278045A (ja) | 2008-05-19 | 2009-11-26 | Sony Corp | 加工体およびその製造方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4201474A (en) | 1978-08-07 | 1980-05-06 | The United States Of America As Represented By The Secretary Of The Army | Variable angle of incidence reflectometer with a continuous read out |
US8160594B2 (en) * | 2006-12-28 | 2012-04-17 | Hitachi, Ltd. | Radio propagation estimating method and radio propagation estimating apparatus |
JP5074333B2 (ja) * | 2008-09-22 | 2012-11-14 | Kddi株式会社 | 電波伝搬特性推定装置及びコンピュータプログラム |
CN102640009B (zh) * | 2009-12-04 | 2016-03-09 | 株式会社Ntt都科摩 | 传输路径推定方法 |
-
2010
- 2010-12-06 EP EP10835940.7A patent/EP2511718B1/en active Active
- 2010-12-06 US US13/514,426 patent/US9337942B2/en active Active
- 2010-12-06 CN CN201080055588.0A patent/CN102648419B/zh active Active
- 2010-12-06 WO PCT/JP2010/071845 patent/WO2011071020A1/ja active Application Filing
- 2010-12-06 JP JP2011545204A patent/JP5687206B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH088846A (ja) * | 1994-06-21 | 1996-01-12 | Hitachi Ltd | 電波受信強度シミュレーション方法 |
JP2003318842A (ja) * | 2002-04-26 | 2003-11-07 | Nippon Telegr & Teleph Corp <Ntt> | 受信電界強度推定計算装置及び方法並びにプログラム及び記録媒体 |
JP2009168534A (ja) | 2008-01-11 | 2009-07-30 | Ntt Docomo Inc | 伝搬特性推定方法及びプログラム |
JP2009278045A (ja) | 2008-05-19 | 2009-11-26 | Sony Corp | 加工体およびその製造方法 |
Non-Patent Citations (6)
Title |
---|
JOUHOUSUURI KENKYUSHO, EEM-RTM RIRON SETSUMEISHO, Retrieved from the Internet <URL:http://www.imslab.co.jp/R Product/eem/doc/rtm_theory.pdf> |
L. LI; Q. CHEN; Q. YUAN; K. SAWAYA; T. MARUYAMA; T. FURUNO; S. UEBAYASHI: "Miscrostrip Reflectarray Using Crossed-Dipole with Frequency Selective Surface of Loops", ISAP2008, 2008 |
T. MARUYAMA; T. FURUNO; S. UEBAYASHI: "Experiment and Analysis of Reflect Beam Direction Control using a Reflector having Periodic Tapered Mushroom-like Structure", 1SAP2008, 2008, pages 9 |
TATSUO FURUNO ET AL.: "A Propagation Estimation Method for Using Beam Controlled Reflector, Tsushin 1, B-1-18", PROCEEDINGS OF THE IEICE CONFERENCE, vol. 2009, - 4 March 2009 (2009-03-04), pages 18, XP008167749 * |
TATSUO FURUNO ET AL.: "A Study of Area Coverage Improvement by using Beam Controlled Reflector, Tsushin Society 1, B-1-19", PROCEEDINGS OF THE IEICE CONFERENCE, vol. 2009, - 1 September 2009 (2009-09-01), pages 19, XP008167757 * |
YOSHIO HOSOYA: "DENPA DENPAN HANDBOOK", REALIZE INC., pages: 234 - 245 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10638479B2 (en) | 2015-11-17 | 2020-04-28 | Futurewei Technologies, Inc. | System and method for multi-source channel estimation |
EP4233909A2 (en) | 2016-05-17 | 2023-08-30 | AbbVie Biotherapeutics Inc. | Anti-cmet antibody drug conjugates and methods for their use |
WO2021026409A1 (en) | 2019-08-08 | 2021-02-11 | Regeneron Pharmaceuticals, Inc. | Novel antigen binding molecule formats |
WO2021091953A1 (en) | 2019-11-05 | 2021-05-14 | Regeneron Pharmaceuticals, Inc. | N-terminal scfv multispecific binding molecules |
WO2021178896A1 (en) | 2020-03-06 | 2021-09-10 | Go Therapeutics, Inc. | Anti-glyco-cd44 antibodies and their uses |
WO2023014863A1 (en) | 2021-08-05 | 2023-02-09 | Go Therapeutics, Inc. | Anti-glyco-muc4 antibodies and their uses |
Also Published As
Publication number | Publication date |
---|---|
US20120302177A1 (en) | 2012-11-29 |
US9337942B2 (en) | 2016-05-10 |
CN102648419B (zh) | 2016-03-16 |
JP5687206B2 (ja) | 2015-03-18 |
EP2511718A1 (en) | 2012-10-17 |
JPWO2011071020A1 (ja) | 2013-04-22 |
EP2511718B1 (en) | 2017-07-05 |
EP2511718A4 (en) | 2015-04-29 |
CN102648419A (zh) | 2012-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5687206B2 (ja) | 伝搬経路推定方法、プログラム及び装置 | |
JP6645346B2 (ja) | 電磁波伝搬経路の追跡方法、装置及びシステム | |
JP4457222B2 (ja) | チャネル特性解析装置及び方法 | |
US8244511B2 (en) | Radio wave propagation characteristic estimation apparatus and computer program | |
CN112154330A (zh) | 用于角度测量的系统和方法 | |
WO2017088814A1 (zh) | 室内定位方法和相关设备 | |
JP6478784B2 (ja) | 電波伝搬シミュレーション装置 | |
JP5355713B2 (ja) | 伝搬経路推定方法、プログラム及び装置 | |
WO2019105194A1 (zh) | 一种传播路径搜索方法以及装置 | |
JP5250062B2 (ja) | 伝搬経路推定方法、プログラム及び装置 | |
CN111801707A (zh) | 一种点云重采样的方法、装置和系统 | |
CN111406222A (zh) | 确定射频信号的到达角 | |
CN113078957A (zh) | 一种基于空间坐标变换技术的uwoc系统信道特征mcns方法 | |
JP2009005250A (ja) | 電波伝搬解析装置 | |
JP4972471B2 (ja) | 電波伝搬解析装置 | |
JP4372287B2 (ja) | 電波伝搬解析方法および装置ならびに記憶媒体 | |
Yang et al. | Robust sequential adaptive Kalman filter algorithm for ultrashort baseline underwater acoustic positioning | |
CN109885940A (zh) | 一种多径条件下波达分布的数学建模方法 | |
JP7511954B2 (ja) | レーダー測位方法、測位レーダー及び測位システム | |
Xing et al. | An identification method of LBL underwater positioning systematic error with optimal selection criterion | |
WO2022155433A1 (en) | Lidar-driven mmwave beam management | |
Honglin et al. | Numerical study of aeroheating predictions for hypersonic reentry bodies | |
JP4157856B2 (ja) | 音響反射経路判別方法、コンピュータプログラム、音響反射経路判別装置及び音響シミュレーション装置 | |
CN118642175A (zh) | 一种海洋声场计算方法、系统、设备及存储介质 | |
Asotov et al. | The algorithms for the ray tracing based on the analytical solutions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201080055588.0 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 10835940 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011545204 Country of ref document: JP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
REEP | Request for entry into the european phase |
Ref document number: 2010835940 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2010835940 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13514426 Country of ref document: US |