CN102640009A - 传输路径推定方法、程序以及装置 - Google Patents

传输路径推定方法、程序以及装置 Download PDF

Info

Publication number
CN102640009A
CN102640009A CN2010800548976A CN201080054897A CN102640009A CN 102640009 A CN102640009 A CN 102640009A CN 2010800548976 A CN2010800548976 A CN 2010800548976A CN 201080054897 A CN201080054897 A CN 201080054897A CN 102640009 A CN102640009 A CN 102640009A
Authority
CN
China
Prior art keywords
mentioned
reflection
transmission path
middle point
rotary middle
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.)
Granted
Application number
CN2010800548976A
Other languages
English (en)
Other versions
CN102640009B (zh
Inventor
丸山珠美
古野辰男
大矢智之
沈纪恽
小田恭弘
陈玉豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NTT Docomo Inc
Original Assignee
NTT Docomo Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NTT Docomo Inc filed Critical NTT Docomo Inc
Publication of CN102640009A publication Critical patent/CN102640009A/zh
Application granted granted Critical
Publication of CN102640009B publication Critical patent/CN102640009B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • 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
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Quality & Reliability (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)
  • Optical Communication System (AREA)
  • Aerials With Secondary Devices (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明的利用成像法的传输路径推定方法具有如下工序：在传输路径上存在反射方向以及散射方向成为与镜面反射的方向θ°不同的预定的方向（θ－η）°的反射阵列1（构造物）时，通过以旋转中心点O为中心，使反射阵列1旋转η/2°来设定假想旋转反射阵列2（假想构造物），使用假想旋转反射阵列2推定传输路径。

Description

传输路径推定方法、程序以及装置
技术领域
本发明涉及传输路径推定方法、程序以及装置。
尤其，本发明涉及在使用几何光学模型的GO或GTD等的分析或使用光线跟踪分析或成像法的分析中，对不进行镜面反射的散射体进行分析的方法。
背景技术
以前，作为分析电波的传输特性的方法，已知使用几何光学模型的分析方法（例如，参照非专利文献1）。
图1中表示基于现有的光线跟踪法求出电波的传输路径（路线）的方法。
如图1所示，在现有的光线跟踪法中，第一，求出与发送点T的壁面相对应的成像点T’，然后求出与成像点T’的壁面2相对应的成像点T”。
第二，用直线连接T”和接收点R，求出这样的直线与壁面2的交点。这样的交点成为壁面2中的反射点X2。
第三，用直线连接壁面2中的反射点X2和T’，能够求出壁面1中的反射点X1。
此时，传输距离成为与连接T”和接收点R的直线的长度相同的值。能够使用该传输距离来求出接收点R的电场强度。
此外，在图2中表示求出关于现有的透过波以及反射波的传输路径（路线）的方法（例如，参照非专利文献2）。
此外，在专利文献1中，公开了推定透过方向并不限于直线传输方向、且反射方向并不限于镜面反射方向的情形时的传输特性的方法。
在先技术文献
专利文献
专利文献1：日本专利特开第2009－168534号公报
非专利文献1：細矢良雄監修、「電波伝搬ハンドブツク」、234～245頁、リアライズインク社
非专利文献2：「EEM-RTM理論説明書」、(株)情報数理研究所、http://www.imslab.co.jp/Rroduct/eem/doc/rtm_theory.pdf
非专利文献3：L.Li,Q.Chen,Q.Yuan,K.Sawaya,T.Maruyama,T.Furunoand S.Uebayashi:“Microstrip Reflectarray Using Crossed-Dipole with FrequencySelective Surface of Loops,”ISAP2008,TP-C05,1645278,2008
非专利文献4：T.Maruyama,T.Furuno and S.Uebayashi:“Experiment andAnalysis of Reflect Beam Direction Control using a Reflector having PeriodicTapered Mushroom-like Structure,”ISAP2008,MO-IS1,1644929,p.9,2008
发明内容
但是，在上述方法中，如图1以及图2所示，反射波的行进方向局限于镜面反射方向（正规反射方向），此外透过波的行进方向限于传输直线前进方向。
从而，在上述的方法中不能考虑到向镜面反射方向以外的方向散射时的传输特性以及介质的折射地进行传输特性的分析。
图3中表示向镜面反射方向以外的方向散射时的传输路径的例子。
在图3中，壁面3以及壁面4分别是方向控制散射体，反射波向与镜面反射方向不同的方向散射，在壁面3从入射方向A1°的方向入射的电波，在反射点×3向反射方向A2°的方向散射，在壁面4，从入射方向B1°的方向入射的电波在反射点×4向反射方向B2°的方向散射。
与此相对，当以现有的成像法求出路径时，在壁面4中的反射点成为×2，成为与反射点×4不同的地方。即、通过现有的图1的方法不能求出壁面4中的反射点。从而，对于方向控制反射板，散射波不沿着镜面反射方向，因此存在不能通过现有的成像法来计算出使用镜像的传输路径的问题点。
另一方面，报告了近年来用于改善传输环境的应用反射阵列或超颖材料来作为传输方向控制散射体的例子（非专利文献3、4），
对在发送点T和接收点R之间的电波的传输路径（路线）的途中存在这样的方向控制散射体时的电波的传输特性进行分析，在分析传输环境的改善效果的方面显得重要，但很难用以前的光线跟踪的方法进行这样的分析。
另一方面，根据专利文献1的方法，通过利用求出路线并删除路线的算法，能够推定传输路径，但无法只通过分析前的构造预先推定。因此，有计算次数增加的缺点。
因此，本发明是鉴于上述的课题而提出的，其目的在于：提供一种传输路径推定方法、程序以及方法，能够推定包含有方向控制散射体的传输分析模型中的传输路径，并进行光线跟踪分析或应用了几何光学模型的分析传输。
本发明的第一特征为，在利用成像法的传输路径推定方法中，具有如下工序：在传输路径上存在反射方向以及散射方向成为与镜面反射的方向θ°不同的方向（θ－η）°的构造物时（η可以是θ的函数），通过以旋转中心点为中心使该构造物旋转η/2°来设定假想构造物，使用该假想构造物推定传输路径。
附图说明
图1是用于说明求出基于现有的光线跟踪法的电波的传输路径（路线）的方法的图。
图2是表示关于现有的透过波以及反射波求出传输路径（路线）的方法的图。
图3是用于说明现有的方法的问题点的图。
图4是用于说明本发明的第一实施方式的传输路径推定方法的图。
图5是用于说明本发明的第二实施方式的传输路径推定方法（第一步骤）的图。
图6是用于说明本发明的第二实施方式的传输路径推定方法（第二步骤）的图。
图7是表示本发明的第二实施方式的传输路径推定方法的流程图。
图8是用于说明本发明的第三实施方式的传输路径推定方法（第一步骤）的图。
图9是用于说明本发明的第三实施方式的传输路径推定方法（第二步骤）的图。
图10是用于说明本发明的第四实施方式的传输路径推定方法（第一）的图。
图11是用于说明本发明的第四实施方式的传输路径推定方法（第二）的图。
图12是用于说明本发明的第四实施方式的传输路径推定方法（第三）的图。
图13是用于说明本发明的第四实施方式的传输路径推定方法（收敛后）的图。
图14是表示基于本发明的第四实施方式的传输路径推定方法的、从反射板旋转中心到反射点为止的距离的收敛状况的图。
图15是用于说明本发明的第五实施方式的传输路径推定方法（反射板旋转前）的图。
图16是用于说明本发明的第五实施方式的传输路径推定方法（第一、旋转45度）的图。
图17是用于说明本发明的第五实施方式的传输路径推定方法（第二、旋转45度）的图。
图18是用于说明本发明的第五实施方式的传输路径推定方法（第三、旋转45度）的图。
图19是用于说明本发明的第五实施方式的传输路径推定方法（收敛后、旋转45度）的图。
图20是表示本发明的第五实施方式的传输路径推定方法的反射点的收敛状况（其1、旋转45度）的图。
图21是表示本发明的第五实施方式的传输路径推定方法的反射点的收敛状况（其2、旋转45度）的图。
图22是用于说明本发明的第五实施方式的传输路径推定方法（第一次、旋转70度）的图。
图23是用于说明本发明的第五实施方式的传输路径推定方法（第二次、旋转70度）的图。
图24是用于说明本发明的第五实施方式的传输路径推定方法（第三次、旋转70度）的图。
图25是用于说明本发明的第五实施方式的传输路径推定方法（第四次、旋转70度）的图。
图26是用于说明本发明的第五实施方式的传输路径推定方法（第五次、旋转70度）的图。
图27是用于说明本发明的第五实施方式的传输路径推定方法（第五次、旋转70度）的图。
图28是表示本发明的第五实施方式的传输路径推定方法的反射点的收敛状况（其1、旋转70度）的图。
图29是表示本发明的第五实施方式的传输路径推定方法的反射点的收敛状况（其2、旋转70度）的图。
具体实施方式
（本发明的第一实施方式）
参照图4，说明本发明的第一实施方式的传输路径推定方法。
图4表示从入射方向θi°（在此，θi°指与一般的电磁场理论中的入射角相同地、反射阵列1的反射面1A的法线方向和入射方向形成的角度）入射到反射阵列1的电波从与反射阵列1相对应的镜面反射方向θr°（在此，θr°指与一般的电磁场理论中的反射角相同地、反射阵列1的反射面1A的法线方向和镜面反射方向形成的角度）旋转η°之后的方向散射的情形。
另外，在本说明书中，通过入射角来规定入射方向，通过反射角来规定反射方向。在图4中，“θi°＝θr°”。
根据现有的成像法，通过用直线连接发送点TX的与反射阵列1相对应的成像点TX’和接收点RX来求出镜面反射时的传输距离。在此，镜面反射时的反射点是连接TX’和接收点RX的直线与反射阵列1相交的点PI。
此外，镜面反射时的传输路径由从发送点TX到PI的直线、和从PI到接收点RX的直线构成。
但是，实际的传输路径由从发送点TX到PR的直线与从PR到接收点RX的直线构成，所以不能通过现有的成像法来求出这样的传输路径以及传输距离。在此，PR是反射阵列1时的反射点。
在此，通过使反射阵列1（构造物）以PR为中心旋转ф°＝η°/2来设定假想旋转反射阵列2（假想构造物）。
在这样的情况下，通过连接发送点TX的与假想旋转反射阵列2相对应的成像点T”和接收点RX，能够求出反射阵列1时的传输距离。在此，ζi°是与假想旋转反射阵列2相对应的入射角，ζr°是与假想旋转反射阵列2相对应的镜面反射角，“ζi°＝ζr°”。
在此，通过以下顺序求出假想旋转反射阵列的旋转角ф°。
将反射阵列1的反射面1A与连接接收点RX和反射阵列1的反射面1A上的反射点PR的线段形成的角设为Ang1时，成为以下公式（1）：
Ang1＝（90－θr+η）°…（1）
将假想旋转反射阵列2的反射面2A与连接接收点RX和反射点PR的线段形成的角设为Ang2时，成为以下公式（2）。
Ang2＝（90－ζr）°…（2）
反射阵列1的反射面1A和假想旋转反射阵列2的反射面2A形成的角是反射阵列的旋转角ф°，所以成立以下公式（3）。
ф°+Ang2＝Ang1…（3）
从而，通过（1）～（3），成为如下公式（4）。
ф°＝ζr°－θr°+η°…（4）
在此，在使ζr°以及θr°满足如下公式（5）、（6）时，公式（7）成立。
ζ°＝ζr°＝ζi°…（5）
θ°＝θr°＝θi°…（6）
ζ°＝（2θ°－η°）/2…（7）
从而，根据公式（4）和（7）的关系，求出“ф°＝η°/2”的关系。
（本发明的第二实施方式）
下面叙述方向控制反射面，例如反射阵列1（方向控制散射体）的反射面1A中的反射点的求出方法。
如上述的第一实施方式中所示，通过使反射阵列1的反射面1A以旋转中心点O为中心仅旋转（η°/2），设置作为假想旋转反射阵列2的反射面的假想旋转反射面2A，能够求出传输路径。
但是，在一般的成像法中求出传输路径，将这样的传输路径和反射面的交点作为反射点。因此，一般很难先求出旋转中心点O（即、反射点）。在本实施方式中，使用图5至图7说明决定这样的反射点的步骤。
<第一步骤>
如图5所示，第一，将反射阵列1的反射面1A的中央点O作为旋转中心OC，使反射面1A仅旋转（η/2）°，来设定作为假想旋转反射阵列2的反射面的假想旋转反射面2A。
第二，对于，假想旋转反射面2A，制作发送点TX的成像点（镜像）TX”。
第三，将连接TX”和接收点RX的直线与假想旋转反射面2A的焦点作为第一步骤的反射点PR1。
能够求出如图5所示那样电波从发送点TX以入射方向ζi1°入射到假想旋转反射面2A，面向接收点RX以反射方向（假想旋转反射面2A中的镜面反射方向）ζr1°反射时的传输路径。
在此，ζi1°是假想旋转反射面2A的法线方向与入射方向形成的角度，ζr1°是假想旋转反射面2A的法线方向与镜面反射方向形成的角度，“ζi1°＝ζr1°”。
第四，如果经过第一步骤的反射点PR1，设定与原来的反射阵列1平行的假想平行移动反射阵列10的反射面即假想平行移动反射面1B，则假想平行移动反射面1B与从发送点TX到第一步骤的反射点PR1的线段之间的角度成为（90－ζ－ф）°（在此，ζ°＝ζi1°＝ζr1°）。
此外，从第一步骤的反射点PR1到接收点RX的线段和假想平行移动反射面1B之间的角度成为（ф+90－ζ）°。这是由于“ф°＝η°/2”，所以与在成为镜面反射方向的角度（90－ζ－ф）°上加上作为反射方向的控制角的η°的值相等。即、通过第一步骤，能够对于假想平行移动反射面1B求出相对于对反射方向控制η°的反射阵列1的传输路径。
但是，在此，第一步骤的反射点PR1位于距离反射阵列1的地方，所以在与反射阵列1相对应的传输路径之间不产生误差。
因此，如以下所示，作为第二步骤，使用图6中所示的方法。
<第二步骤>
如图6所示，第一，通过以将步骤1的旋转中心点OC加上预先决定的△s后的点作为旋转中心点OC，使反射阵列1的反射面1A仅旋转（η/2）°，由此设定第二假想旋转反射面2B，通过与第一步骤相同的方法求出第二反射点PR2。
即、成为“（第二次旋转中心OC）＝（第一次旋转中心OC）+△s＝O+△s”。
在此，在这样的反射点PR2和反射阵列1的反射面1A之间的距离大于预定的收敛条件（例如△s）的情况下，第二，以点（O+2×△s）作为旋转中心点OC，使反射阵列1的反射面1A仅旋转（η/2）°，由此设定第三假想旋转反射面1C，到满足收敛条件为止重复上述工序。
即、“（第三次旋转中心OC）＝（第二次旋转中心OC）+△s＝O+2△s”。
以下，参照图7中所示的流程图，整理本实施方式的传输路径推定方法中的动作。
如图7所示，在步骤S101中，决定旋转中心点OC（例如反射阵列1的中央点O）。
在步骤S102中，以旋转中心点O为中心，使目前已设定的反射面（反射阵列1的反射面1A或第n假想平行移动反射面）旋转，由此设定假想旋转反射面。
在步骤S103中，根据使用这样的假想旋转反射面的成像法的光线跟踪，决定传送路径、传送距离或假想旋转反射面上的反射点。
在步骤S104中，计算出所决定的假想旋转反射面上的反射点和反射阵列1的反射面1A之间的距离d。
在步骤S105中，判定是否满足“d<判定值（例如△s）”。
在“YES”的情况下，结束本动作，在“NO”的情况下，本动作回到步骤S101。
（本发明的第三实施方式）
参照图8以及图9，说明本发明的第三实施方式的传输路径推定方法。
本发明的第三实施方式的传输路径推定方法的特征在于，将在上述第二实施方式的传输路径推定方法中求出的第n步骤中的反射点PRn和第n步骤中的旋转中心点（例如，反射阵列1的中央点）O之间的距离设为△Ln，将第（n+1）次旋转中心点设为“O+△Ln（O为0的情况下，为△Ln）”。
图8是用于说明第三实施方式的传输路径推定方法中的第一步骤的图，求出△L1的值来作为与第一步骤的旋转中心点（例如，反射阵列1的中央点）O之间的距离。
图9是用于说明第三实施方式的传输路径推定方法中的第二步骤的图，使反射面1A以点（O+△Ln）为旋转中心点仅旋转“ф°＝η°/2”。在此，如果原来的反射阵列1的反射面1A和第二步骤中的反射点PR2之间的距离d满足收敛条件，则结束处理。
另外，可以将上述的实施方式的传输路径推定方法应用到存在反射阵列1以外的构造物的情形。此外，上述的实施方式的传输路径推定方法可以在存在多个经由反射阵列1的传输路径（路线）时，对于每个路线设定假想旋转反射面。
（本发明的第四实施方式）
参照图10至图13，说明本发明的第四实施方式的传输路径推定方法。在本实施方式中说明传输路径推定方法的具体的数值计算例子。
图10表示相对于反射阵列1的反射面1A，以旋转中心O1为中心旋转η/2角度的情形。在本计算例子中考虑将最初的旋转中心O1作为坐标的原点的二维坐标。将发送点的坐标设为（－10，－30），将接收点的坐标设为（30，－30），将旋转角η设为－45度（向逆时针方向45度）。
对于反射阵列1的反射面1A，将旋转中心O1设为中心，将相对于反射阵列1旋转旋转角η/2、即旋转－22.5度而得的面设为假想旋转反射阵列的反射面2A。
接着，对于假想旋转反射阵列的反射面2A，制作发送点Tx的镜像Tx’，将连接Tx’和接收点Rx的直线与2A相交的点设为PR1。此时，PR1的坐标是（-6.21，－2.57）。计算第一次旋转中心O1和第一反射位置PR1的距离的话是6.73。
利用图11，说明下一个步骤。在下一个步骤，求出从旋转中心O1沿着旋转反射阵列的反射面1A仅移动-6.73的距离之后的点（-6.73，0），将该点作为第二次旋转中心O2。制作以O2为旋转中心旋转η/2的角度之后的假想旋转反射阵列。即、以假想旋转反射阵列的旋转中心与O2重叠的方式使假想旋转反射阵列在横向平行移动。对于以该O2为中心旋转后的假想反射阵列，制作发送点Tx的镜像。另外，假想旋转反射阵列的旋转中心不必一定与假想旋转反射阵列的反射面的长度方向中的中点一致，从旋转中心到假想旋转反射阵列的一端的距离可以与从旋转中心到假想旋转反射阵列的另一端的距离不同。
此外，可以将从假想反射板的中点仅移动－6.73之后的点作为假想反射板的旋转中心。
将与连接Tx的镜像和接收点Rx的直线相交的交点设为PR2。此时，PR2的坐标是（－6.89，0.066）。计算第二次旋转中心O2和第二反射位置PR2的距离的话，是0.17。
利用图12说明后续的步骤。将使旋转中心从O2沿着反射阵列的反射面1A仅移动－0.17的距离之后的点（－6.90，0）设为第三次旋转中心O3。制作以O3为旋转中心，旋转η/2的角度之后的假想旋转反射阵列。即、以假想旋转反射阵列的旋转中心与O3重叠的方式使假想旋转反射阵列在横向平行移动。对于以该O3为中心旋转后的假想反射阵列，制作发送点Tx的镜像。
将与连接Tx的镜像和接收点Rx的直线相交的交点设为PR3。此时，PR3的坐标是（－6.90，0.002）。计算第二次旋转中心O2和第二次反射位置PR3的距离的话0.005。另外，在本实施方式中，该距离的收敛条件是ε<0.03。
图13表示本实施方式的计算收敛的状况的发送点、接收点以及反射点的位置关系。可知能够求出如图13所示那样从发送点Tx发射的电波经过发射点PR3，向从正规反射方向仅旋转45度的角度之后的方向行进，到达Rx的传输路线。
图14表示从旋转中心到反射点的距离的收敛状况。可知如图14所示，伴随着试行次数的增加，从旋转中心On（n为1、2、3的整数）到反射点PRn的距离急速收敛。
如以上说明的那样，通过重复实施上述的步骤，能够推定传输路径。在上述的计算例子中使用二维坐标系，但即使是三维坐标系，也可以用同样的步骤进行计算。
（本发明的第五实施方式）
参照图15至图29，说明本发明的第五实施方式的传输推定方法。在本实施方式中说明在发送点和接收点之间电波两次反射时的传输路径推定方法的具体的数值计算例子。
图15是说明存在包含反射阵列的反射面的两个反射面、且反射次数为两次时的传输路径的求出方法。在图15中表示了将反射阵列的反射面1A的一端设为原点（0，0）的二维坐标。在本实施方式中说明以下的情形，即将发送点Tx的坐标设为（0，－100），将接收点Rx的坐标设为（40，－40），从发送点Tx开始在两个反射面中分别反射一次、即进行合计两次的反射。
将反射阵列的反射面1A设为第一反射面，将通常的正规反射的壁面设为第二反射面。制作与发送点Tx的反射面1B相对应的成像点Tx’以及与Tx’的反射面1A相对应的成像点Tx”。此时，Tx”的坐标为（－60，－100）。
图15表示使反射阵列的反射面旋转之前的状态。图16表示使反射面1A旋转时的计算例子。在本实施方式中说明将旋转角η设为45度，使假想反射面2A仅旋转η/2=22.5度时的情形。此时，将反射面1A的中点（35，0）设为第一次旋转中心。旋转后的反射点的坐标是（39.69，－1.94）。
使用图16说明后续的步骤。第一次旋转中心（35，0）和第一反射点的坐标（39.69，－1.94）的距离是5.07。第二次旋转中心是向与第一反射点之间的距离变小的方向仅移动5.07而得的点，具体地说将（40.07，0）设为第二次旋转中心。
如图17所示，计算以第二次旋转中心为中心仅旋转η/2=22.5度而得的假想反射面2A，求出此时的反射点坐标的话成为（40.00，0.03）。
接着，同样地对第三次的旋转进行计算。第二次旋转中心（40.07，0）和第二反射点（40.00，0.003）之间的距离是0.073。将从第二次旋转中心向反射点的距离变小的方向仅移动0.073而得的点设为第三次旋转中心。此时，第三次旋转中心的坐标是（40.00，0）。如图18所示，此时反射点的坐标是（40.00，0.00）。
图19表示计算收敛后的传输路线。与进行正规反射时的通常的传输路线相比，通过计算求出的反射阵列1A上的反射点到Rx的传输路线成为使反射面向顺时针方向仅旋转45度后的传输路线，这与设计值η相等。
图20以及图21表示传输路线的收敛状况。图20的纵轴是旋转中心的X坐标，横轴是试行次数。图21的纵轴是在第n次试行中的旋转中心的坐标与第n+1次旋转中心的坐标之间的距离，横轴是试行次数n。
如图20以及图21所示，可知通过两次的试行次数，已经大致收敛。另外，在本实施方式中，相应距离的收敛条件是ε<0.03。
图22～29表示对于旋转角η＝70度进行与上述同样的计算的情形。当把第一次旋转中心作为反射阵列1A的反射面的中心（35，0），进行η/2=35度的旋转，求出第一次反射点的坐标时，该坐标成为（50.15，－10.61）。
从第一次旋转中心到反射点的距离是18.50。通过将第一次旋转中心向与反射点的距离变短的方向仅移动18.50，能够求出第二次旋转中心。此时，第二次旋转中心是（53.50，0）。
在后续的步骤中将第二次旋转中心作为中心，使反射面旋转η/2=35°来求出反射点。通过重复这些步骤，即使在如图27所示，η＝70度的情况下，也能够同样地求出旋转角70的传输路线。
图28以及图29与图20以及图21同样地表示传输路线的收敛状况，可知通过增加试行次数而急速地收敛传输。
可以如下地表现以上所述的本实施方式的特征。
本实施方式的第一特征是使用成像法的传输路径推定方法，其主旨在于，具有如下工序：在传输路径上存在反射方向以及散射方向成为与镜面反射的方向θ°不同的方向（θ－η）°的构造物时，通过以旋转中心点O为中心使反射阵列1旋转η/2°来设定假想旋转反射阵列2（假想构造物），使用假想旋转反射阵列2推定传输路径。
在本实施方式的第一特征中，可以到反射阵列1上的反射点PRn与旋转中心点O一致为止，移动旋转中心点O，重复上述工序。
在本实施方式的第一特征中，上述工序可以具有以下工序：决定旋转中心点O的工序A；通过以旋转中心点O为中心，使反射阵列1旋转γ次（例如1次）来设定假想旋转反射阵列2的工序B；应用成像法，对于假想旋转阵列2求出假想旋转阵列2上的反射点PR的工序C；以及在反射点PR和旋转中心点O之间的距离不满足收敛条件时，将旋转中心点O更新为“O－∣PR-O∣”的工序D，到反射点PR和旋转中心点O之间的距离满足收敛条件为止重复工序B至D。
在本实施方式的第一特征中，上述工序可以具有以下工序：决定旋转中心点O的工序A；通过以旋转中心点O为中心，使反射阵列1旋转γ次（例如1次）来设定假想旋转反射阵列2的工序B；应用成像法，对于假想旋转阵列2求出假想旋转阵列2上的反射点PR的工序C；以及在反射点PR和旋转中心点O之间的距离不满足收敛条件时，将旋转中心点O更新为“O－△s×n”的工序D，到反射点PR和旋转中心点O之间的距离满足收敛条件为止重复工序B至D。
在本实施方式的第一特征中，η°可以是θ°的函数。
本实施方式的第二特征的要旨在于使计算机实现上述的传输路径推定方法的程序或装置。
即、可以通过硬件来实施上述的传输路径推定方法，也可以通过处理器执行软件模块来实施上述的传输路径推定方法，也可以通过两者的组合来实施。
软件模块可以设在RAM（Random Access Memory）、闪存器、ROM（ReadOnly Memory）、EPROM（Erasable Programmable ROM）、EEPROM（Electronically Erasable and Programmable ROM）、寄存器、硬盘、可移动磁盘、或CD-ROM中的任意形式的存储介质内。
这样的存储介质与相应处理器连接，以便处理器能够对于相应存储介质读写信息。此外，这样的存储介质可以集成在处理器内。此外，这样的存储介质以及处理器可以被设在ASIC内。
以上使用上述的实施方式，详细地说明了本发明，但本领域技术人员明了本发明并不限于在本说明书中说明的实施方式。本发明可以在不脱离通过权利要求请求记载而确定的本发明的主旨以及范围的前提下进行修正以及作为变更形态来实施。从而，本说明书的记载以例示说明作为目的，并不对本发明具有任何限制的意思。
另外，通过参照而将日本专利申请的2009－277070号（2009年12月4日申请）的全部内容组合到本说明书中。
产业上的利用可能性
根据本发明，能够提供一种传输路径推定方法、程序以及装置，能够推定在包含有方向控制散射体的传输分析模型中的传输路径，并进行光线跟踪分析或应用几何光学模型的分析传输，因此在无线通信等中有用。
符号说明
1…反射阵列
1A…反射面
1B…假想平行移动反射面
2…假想旋转反射阵列
2A、2B…假想旋转反射面

Claims (7)

1.一种传输路径推定方法，是利用成像法的传输路径推定方法，其特征在于，
具有如下工序：在传输路径上存在反射方向以及散射方向成为与镜面反射的方向θ°不同的方向（θ－η）°的构造物时，通过以旋转中心点为中心使该构造物旋转η/2°来设定假想构造物，使用该假想构造物推定传输路径。
2.根据权利要求1所述的传输路径推定方法，其特征在于，
到上述构造物上的反射点与上述旋转中心点一致为止，移动该旋转中心点，重复上述工序。
3.根据权利要求1或2所述的传输路径推定方法，其特征在于，
上述工序具有以下工序：
决定上述旋转中心点O的工序A；
通过以上述旋转中心点O为中心，使上述构造物旋转γ次来设定上述假想构造物的工序B；
应用上述成像法，对于上述假想构造物算出该假想构造物上的反射点PR的工序C；以及
在上述反射点PR和上述旋转中心点O之间的距离不满足收敛条件时，将上述旋转中心点O更新为“O－∣PR-O∣”的工序D，
到上述反射点PR和上述旋转中心点O之间的距离满足收敛条件为止重复上述工序B至D。
4.根据权利要求1或2所述的传输路径推定方法，其特征在于，
上述工序具有以下工序：
决定上述旋转中心点O的工序A；
通过以上述旋转中心点为中心，使上述构造物旋转γ次来设定上述假想构造物的工序B；
应用上述成像法，对于上述假想构造物求出该假想构造物上的反射点PR的工序C；以及
在上述反射点PR和上述旋转中心点O之间的距离不满足收敛条件时，将上述旋转中心点O更新为“O－△s×n”的工序D，
到上述反射点PR和上述旋转中心点O之间的距离满足收敛条件为止重复上述工序B至D。
5.根据权利要求1至4中的任意一项所述的传输路径推定方法，其特征在于，
上述η°是上述θ°的函数。
6.一种程序，其特征在于，
使计算机实现权利要求1至5中的任意一项所述的传输路径推定方法。
7.一种装置，其特征在于，
使计算机实现权利要求1至5中的任意一项所述的传输路径推定方法。
CN201080054897.6A 2009-12-04 2010-12-03 传输路径推定方法 Active CN102640009B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2009277070 2009-12-04
JP2009-277070 2009-12-04
PCT/JP2010/071742 WO2011068224A1 (ja) 2009-12-04 2010-12-03 伝搬経路推定方法、プログラム及び装置

Publications (2)

Publication Number Publication Date
CN102640009A true CN102640009A (zh) 2012-08-15
CN102640009B CN102640009B (zh) 2016-03-09

Family

ID=44115069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201080054897.6A Active CN102640009B (zh) 2009-12-04 2010-12-03 传输路径推定方法

Country Status (5)

Country Link
US (1) US9264157B2 (zh)
EP (1) EP2511717B1 (zh)
JP (1) JP5355713B2 (zh)
CN (1) CN102640009B (zh)
WO (1) WO2011068224A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5355713B2 (ja) 2009-12-04 2013-11-27 株式会社エヌ・ティ・ティ・ドコモ 伝搬経路推定方法、プログラム及び装置
WO2011071020A1 (ja) * 2009-12-07 2011-06-16 株式会社エヌ・ティ・ティ・ドコモ 伝搬経路推定方法、プログラム及び装置
JP2022108025A (ja) * 2021-01-12 2022-07-25 電気興業株式会社 リフレクトアレー、リフレクトアレーの設計方法、および、リフレクトアレーシステム

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5689812A (en) * 1994-06-21 1997-11-18 Hitachi, Ltd. Radio propagation simulation method
CN1235435A (zh) * 1998-04-23 1999-11-17 三星电子株式会社 使用城区峡谷模型的无线电波传播预测方法
JP2003318842A (ja) * 2002-04-26 2003-11-07 Nippon Telegr & Teleph Corp <Ntt> 受信電界強度推定計算装置及び方法並びにプログラム及び記録媒体
JP2009168534A (ja) * 2008-01-11 2009-07-30 Ntt Docomo Inc 伝搬特性推定方法及びプログラム
CN101592690A (zh) * 2009-05-05 2009-12-02 上海大学 基于射线跟踪法的电磁波传播预测方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
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
JP5344457B2 (ja) 2008-05-15 2013-11-20 キャタピラー エス エー アール エル 操作体装置用カバーおよび操作弁装置
JP5074333B2 (ja) 2008-09-22 2012-11-14 Ｋｄｄｉ株式会社 電波伝搬特性推定装置及びコンピュータプログラム
JP5355713B2 (ja) 2009-12-04 2013-11-27 株式会社エヌ・ティ・ティ・ドコモ 伝搬経路推定方法、プログラム及び装置
WO2011071020A1 (ja) 2009-12-07 2011-06-16 株式会社エヌ・ティ・ティ・ドコモ 伝搬経路推定方法、プログラム及び装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5689812A (en) * 1994-06-21 1997-11-18 Hitachi, Ltd. Radio propagation simulation method
CN1235435A (zh) * 1998-04-23 1999-11-17 三星电子株式会社 使用城区峡谷模型的无线电波传播预测方法
JP2003318842A (ja) * 2002-04-26 2003-11-07 Nippon Telegr & Teleph Corp <Ntt> 受信電界強度推定計算装置及び方法並びにプログラム及び記録媒体
JP2009168534A (ja) * 2008-01-11 2009-07-30 Ntt Docomo Inc 伝搬特性推定方法及びプログラム
CN101592690A (zh) * 2009-05-05 2009-12-02 上海大学 基于射线跟踪法的电磁波传播预测方法

Also Published As

Publication number Publication date
WO2011068224A1 (ja) 2011-06-09
EP2511717A4 (en) 2015-04-29
JPWO2011068224A1 (ja) 2013-04-18
US20120303328A1 (en) 2012-11-29
US9264157B2 (en) 2016-02-16
EP2511717B1 (en) 2017-07-05
CN102640009B (zh) 2016-03-09
JP5355713B2 (ja) 2013-11-27
EP2511717A1 (en) 2012-10-17

Similar Documents

Publication Publication Date Title
CN102648419A (zh) 传输路径推定方法、程序以及装置
CN102053245B (zh) 车载雷达装置和控制车载雷达装置的方法
Rizzo et al. Signal-based deployment planning for robot teams in tunnel-like fading environments
CN105430740A (zh) 基于WiFi信号强度仿真与位置指纹算法的室内无线定位方法
CN102640009A (zh) 传输路径推定方法、程序以及装置
CN110018440B (zh) 基于定向天线的l型直角阵列测向方法、装置和系统
Petrov et al. Mobile near-field terahertz communications for 6G and 7G networks: Research challenges
Gubelli et al. Ray-tracing simulator for radar signals propagation in radar networks
CN107241150B (zh) 一种两径无线信道相关时间的测量装置及测量方法
Lemic et al. Localization in power-constrained terahertz-operating software-defined metamaterials
CN115877328A (zh) 一种阵列雷达的信号收发方法及阵列雷达
JP5250062B2 (ja) 伝搬経路推定方法、プログラム及び装置
CN104977560B (zh) 移动设备、定位方法和定位系统
CN109117575A (zh) 表面等离激元波导系统的结构参数确定方法及设备
Jakobus et al. Recent advances of FEKO and WinProp
CN105550436B (zh) 一种融合波模和射线理论的弯曲巷道电波传播建模方法
KR20080000936A (ko) 삼각형 셀을 이용한 3차원 광선 추적 방법 및 그를 이용한전파 특성 예측 방법
US20160274714A1 (en) Touch Identification Device on the Basis of Doppler Effect, Touch Identification Method on the Basis of Doppler Effect and Touch Screen
Shu et al. Link failure rate and speed of nodes in wireless network
Bai et al. Channel modeling by ray tracing for city underground power pipe gallery
Tomasi Reliable communication in mine environments for autonomous vehicles
JP2004356680A (ja) マイクロ波帯移動電波伝搬推定方法
Serghiou Scattering Models for Fixed sub-THz Wireless Links
Costanzo et al. Through-the-wall short-range sensing by Bessel beams sources
Taghvaee et al. Fully Autonomous Reconfigurable Metasurfaces with Integrated Sensing and Communication

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant