CN104573368B - 基于面元投影的三角形截面射线管电磁射线追踪方法 - Google Patents
基于面元投影的三角形截面射线管电磁射线追踪方法 Download PDFInfo
- Publication number
- CN104573368B CN104573368B CN201510017244.5A CN201510017244A CN104573368B CN 104573368 B CN104573368 B CN 104573368B CN 201510017244 A CN201510017244 A CN 201510017244A CN 104573368 B CN104573368 B CN 104573368B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- bin
- ray tube
- 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
Landscapes
- Image Generation (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
一种基于面元投影的三角形截面射线管电磁射线追踪方法，步骤如下：1、将目标表面划分成三角形面元来近似；2、判断该面元的朝向，面元法向朝向称亮面元，背向为暗面元；3、亮面元之间进行遮挡判断，剔除被遮挡的亮面元；4、设定经过原点、与入射方向垂直的平面，称为投影面；5、将亮面元向投影面投影，以投影三角形为横截面，形成沿电磁波入射方向的射线管；6、判断射线管与目标表面面元的求交情况，计算反射射线，进而计算每次反射所形成的RCS，第一次反射面就是形成射线管的亮面元；7、将射线管的反射形成的RCS进行叠加，即得总的RCS。该方法将射线管与目标面元加以关联，消除了传统矩形射线管与三角形面元相交形成的误差。
Description
一、技术领域
本发明一种基于面元投影的三角形截面射线管电磁射线追踪方法，属于电磁散射仿真分析领域。
二、背景技术
射线追踪法是计算目标高频区电磁散射的一种电磁仿真方法，主要方法是将入射波划分为许多射线管，追踪每条射线管在目标表面的反射情况，射线管在目标表面的反射方向、幅度、相位基于几何光学法进行计算。射线管每次照射到目标表面，都会形成电磁散射，将所有射线管在目标表面形成的电磁散射进行叠加，就得到总的电磁散射，进而可以求出目标的雷达散射截面(radar cross section，RCS)。射线追踪法在计算大尺寸雷达散射目标时非常有效，能够计算目标各个部件的耦合散射。
传统的射线追踪采用的方法是首先根据目标的形状、尺寸设定一个垂直于入射方向的矩形口面，再将此矩形口面划分成许多很小的正方形，每个正方形沿电磁波入射方向形成一条射线管。射线追踪的过程就是计算此正方形射线管在目标表面的多次反射路线，进而计算出目标的雷达散射截面。在射线追踪的过程中，计算量最大的部分就是射线与面元的求交判断，所有射线管与所有面元都要进行，最后还要判断是否有遮挡关系。一般来说，要计算的目标的形状是不规则的，因而许多射线管没有与目标相交，但在计算过程中仍然要与面元进行求交判断，耗费了大量的计算时间。
本发明提出一种基于三角形面元的射线管方法，射线管的设置不再通过入射口面划分的方式，而是依据目标表面的三角形面元进行设定。将目标表面用许多小的三角形面元近似，再将朝向雷达方向的面元沿入射方向投影在平面上，即得到三角形射线管，这些射线管会覆盖目标全部。如果目标各面元之间没有互相遮挡关系，则三角形射线管也不会重叠。如果面元之间互相遮挡，则剔除被遮挡的面元之后，形成的三角形射线管也不会重叠。这些射线管恰好覆盖了散射目标，不再有多余的射线管进行求交判断。对每个三角形射线管的反射射线进行追踪，再经进一步物理光学积分计算，即可得到目标的雷达散射截面。
三、发明内容
本发明的目的是针对传统射线追踪法射线管较多、射线追踪判断计算量大的特点，提出了一种基于面元投影的三角形截面射线管电磁射线追踪方法，它能够降低计算量、降低内存耗费，同时提高精度。
本发明一种基于面元投影的三角形截面射线管电磁射线追踪方法，具体包括以下步骤：
步骤1：将要计算的目标表面划分成许多三角形面元来近似，用三角形面元的优点是网格划分较为灵活，而且任何三个顶点都在一个平面上，方便后续物理光学积分计算；
步骤2：判断三角形面元的朝向，面元法向朝向入射电磁波方向的面元称为亮面元，背向入射电磁波朝向的面元称为暗面元；
步骤3：所有的亮面元之间进行遮挡判断，剔除被遮挡的亮面元；
步骤4：根据电磁波入射方向设定一个经过原点、与入射方向垂直的平面，称为投影面；
步骤5：将经过步骤3剔除后的亮面元向投影面投影，投影也是三角形形状，这些投影三角形恰好覆盖了目标的投影轮廓，以这些投影三角形为横截面，形成沿电磁波入射方向的射线管；
步骤6：判断每条射线管与目标表面各个面元的求交情况，计算反射射线，进而计算每次反射所形成的RCS，其中第一次反射面就是形成射线管的亮面元，无需进行求交判断；
步骤7：将所有射线管的所有反射形成的RCS进行叠加，即得到目标总的RCS。
其中，在步骤1中所述的“将要计算的目标表面划分成许多三角形面元来近似”，所述的“划分成许多三角形面元”，是指目标形状用封闭的多面体表示，每个面都是三角形；所述的“近似”，是指三角形面元边长小于目标表面曲率半径的1/8，同时也要小于入射电磁波波长的1/5，这时能够较好的体现目标形状。
其中，在步骤2中所述的“判断三角形面元的朝向”，其说明如下：看单独一个面元是否能被入射波照射，如可以照射，说明其法向与入射波方向同向，点积小于0，即其中为面元法向，为入射波方向，此面元归类为亮面元；如不能照射，说明其法向与入射波方向反向，点积大于0，即此面元归类为暗面元；对的临界情形，按照暗面元处理。
其中，在步骤3中所述的“进行遮挡判断”，其说明如下：每两个亮面元之间都要通过计算看是否相对于入射电磁波有遮挡关系，如一个面元的重心被另一个面元遮挡，则电磁波不能照射到该重心；用图3表示了遮挡示意图，如图3所示，r0是待判断遮挡的点。设L1＝AB，L2＝AC，三角形ABC三个顶点的位置矢量分别为rA、rB、rC，则三角形所在平面上任一点r可以表示为r＝rA+α1L1+α2L2。如果同时满足条件α1＞0、α2＞0、α1+α2＜1，则r在三角形内，否则r在三角形外。
由于r是是经过r0的入射射线上的一点，r也可以用r0、入射方向表示，即其中s表示从r到r0的位移。s＞0时，r沿入射方向在r0的后面，此时三角形面元可能遮挡r0；s＜0时，r沿入射方向在r0的前面，此时三角形面元不可能遮挡r0；s＝0时，为临界情况，r0在三角形面元ABC上。在遮挡判断过程中，临界情况认为面元遮挡无效。
根据r的两种表达式，可以得出如下关系：进一步可以写为线性方程组的形式：
其中各矢量都是列向量，因此系数矩阵是一个方阵，在矩阵非奇异情况下可以求得三个未知数α1、α2、s。当同时满足条件α1＞0、α2＞0、α1+α2＜1、s＞0时，认为r0被面元遮挡，否则r0不被面元遮挡。
再考虑系数矩阵奇异的情形。当奇异，即行列式为零时，表示三个矢量在同一平面上，此时入射方向与三角形面元ABC平行，这种情况下r0也不会被面元遮挡；所述的“剔除被遮挡的亮面元”，是指如果面元的重心被遮挡，则其不参与物理光学积分计算，对电磁散射没有贡献。
其中，在步骤5中所述的“向投影面投影”，其说明如下：投影面是过坐标原点并与入射方向垂直的平面，各面元沿入射方向延伸，与投影面的交点即是投影；设需要投影的三角形面元三个顶点位置矢量是rA,rB,rC，重心是r0，将三角形各顶点和重心在投影面上的投影位置矢量分别是rA′,rB′,rC′，r0′。以r0的投影计算为例进行说明，r0′满足的关系为由此可以计算得三角形的顶点位置矢量在投影面上的投影同样可以计算得到
其中，在步骤6中所述的“判断”，其说明如下：三角形射线管沿射线方向延伸，与面元的交点计算；所述的“计算反射射线”，其说明如下：通过面元法线入射方向求出反射方向结果为所述的“进而计算”，其说明如下：计算射线管照射到面元上的RCS贡献，第一次反射采用公式
式中Lm是第m条边的矢量，rmc是第m条边中点的位置矢量，A是面元的面积，r0是三角形面元的重心。当时，意味着散射方向是入射方向对于反射面的镜面反射方向；第二次反射采用公式
式中rmc1是第m条边的中点rm沿与二次反射面的交点。是三角形射线管在二次反射面上投影的中心。
其中，在步骤7中所述的“进行叠加”，其说明如下：所有照射到的面元RCS按照相位求和，得到总的RCS为其中N为所有面元总的反射次数，表示每次反射的RCS贡献。
本发明的方法相对于原有射线追踪法，有以下创新及优点：1)避免了传统方式的射线管与面元的求交判断，大大降低了计算量；2)消除了传统方法中与目标不相交的大量射线管，节省了内存耗费；3)射线管的形状根据目标表面的面元进行划分，恰好覆盖了目标全部，而传统方法的正方形射线管在目标轮廓处求交时一般总会有误差存在。
四、附图说明
图1是本发明基于面元投影的三角形截面射线管电磁射线追踪方法整体步骤流程图。
图2是二面角划分成三角面元示意图。
图3是面元遮挡判断示意图。
图4是面元遮挡判断流程图。
图5是三角形射线管多次反射示意图。
图6(a)是HH极化计算结果。
图6(b)是VV极化计算结果。
附图中符号说明如下：
A,B,C——三角形面元三个顶点
L1——三角形面元AB边的矢量
L2——三角形面元AC边的矢量
r——照射到r0的电磁波射线与三角形面元ABC平面的交点
α1——三角形面元ABC平面上从顶点A到r沿CA在AB边投影的长度
α2——三角形面元ABC平面上从顶点A到r沿BA在AC边投影的长度
——电磁波入射方向
r0——待判断遮挡的三角形面元重心
s——r0到r的长度
五、具体实施方式
下面将结合附图和实施例对本发明做进一步的详细说明。
如图1所示，本发明一种基于面元投影的三角形截面射线管电磁射线追踪方法，包括如下步骤：
步骤1：目标表面划分三角形面元：将要计算的目标表面划分成许多三角形面元进行近似，用三角形面元的优点是网格划分较为灵活，而且任何三个顶点都在一个平面上，方便后续物理光学积分计算。
步骤2：判断三角面元朝向，分成亮面元和暗面元：根据面元法向判断其朝向，法向朝向入射电磁波方向的面元称为亮面元，背向入射电磁波方向的面元称为暗面元。通过这种朝向判断，暗面元就不需要参与计算，也不需要进行遮挡判断，计算量降低一半。
步骤3：所有亮面元之间进行遮挡判断，剔除被遮挡的亮面元：目标表面的面元的遮挡判断只需要考虑两面元之间的遮挡关系就行，因为目标是一个封闭体，被亮面元遮挡就意味着被某个部件所遮挡。
步骤4：设定垂直于入射方向的投影面：根据电磁波入射方向设定一个经过原点、与入射方向垂直的平面，称为投影面。传统方法中该投影面要离开目标一定距离，不能与目标相交，但本方法中投影面可以经过原点，并不影响计算结果。
步骤5：没有被遮挡的亮面元向投影面投影，形成三角形射线管：将经过步骤3剔除后的亮面元向投影面投影，形状也是三角形，这些三角形恰好覆盖了目标的投影轮廓，以这些投影三角形为横截面，形成沿电磁波入射方向的射线管。这些射线管形状、大小不一，根据几何光学和物理光学假设，只要射线管能够覆盖目标整体即可，形状、尺寸不影响最终计算结果。
步骤6：每条射线管与每个面元进行求交判断，计算每次反射形成的RCS：判断每条射线管与目标表面各个面元的求交情况，计算反射射线，进而计算每次反射所形成的RCS，其中第一次反射面就是形成射线管的亮面元，无需进行求交判断。用传统正方形射线管的方法在第一次反射时也需要求交判断，而本方法则省去了第一次反射的求交判断，这是本方法能够减少计算量的关键所在。
步骤7：所有RCS贡献叠加，得到目标总的RCS：将所有射线管的所有反射波形成的RCS进行叠加，即得到目标总的RCS。在计算过程中，需要保留相位信息，以便在最后叠加的时候能够反映出各个面元散射波的相干效果。
针对上述步骤的进一步说明如下：
步骤1：目标表面划分三角形面元
目标表面划分三角形面元的优点是简单、灵活，而且任意的三角形面元都是平面形状，便于后续的反射射线、物理光学积分计算。图2是二面角剖分成三角面元的示意图。
目标表面划分成三角面元后，网格信息用顶点数组和面元数组表示。顶点数组为N行3列的浮点数，其中N是顶点数量，每一行代表一个顶点，3列分别表示每个顶点的x、y、z坐标。面元数组为M行3列的整数，其中M是面元数量，每一行代表一个面元，3列分别表示面元3个顶点的序号。
步骤2：判断三角面元朝向，分成亮面元和暗面元
判断面元的亮暗可以用法向与入射方向的点积进行。设是面元的法向单位矢量，是入射方向单位矢量，则如果表示面元是暗面元，入射方向与法线方向朝向同一方向，如果表示面元是暗面元，入射方向与发现方向朝向相反方向。如果表示入射方向和面元平行，这是一种临界状况，本发明在处理过程中，将这种情况归类为暗面元，不参与电磁散射计算。在物理光学的RCS计算公式中，也可看出如果则RCS计算结果为0，说明这种处理方法合理。
步骤3：所有亮面元之间进行遮挡判断，剔除被遮挡的亮面元
在此步骤中，需要进行面元之间的遮挡判断，这也是该方法计算量较大的部分之一。由于三角形面元尺寸较小，在处理过程中可近似认为只要三角形面元的中心被遮挡，则整个面元被遮挡。这种近似处理在面元较小时，可以得到较高的精度。因此，问题就归结为一个点被一个三角形面元遮挡的问题。图3是点被面元遮挡判断的示意图。
如图3所示，r0是待判断遮挡的点。设L1＝AB，L2＝AC，三角形ABC三个顶点的位置矢量分别为rA、rB、rC，则三角形所在平面上任一点r可以表示为r＝rA+α1L1+α2L2。如果同时满足条件α1＞0、α2＞0、α1+α2＜1，则r在三角形内，否则r在三角形外。
由于r是是经过r0的入射射线上的一点，r也可以用r0、入射方向表示，即其中s表示从r到r0的位移。s＞0时，r沿入射方向在r0的后面，此时三角形面元可能遮挡r0；s＜0时，r沿入射方向在r0的前面，此时三角形面元不可能遮挡r0；s＝0时，为临界情况，r0在三角形面元ABC上。在遮挡判断过程中，临界情况认为面元遮挡无效。
根据r的两种表达式，可以得出如下关系：进一步可以写为线性方程组的形式：
其中各矢量都是列向量，因此系数矩阵是一个方阵，在矩阵非奇异情况下可以求得三个未知数α1、α2、s。当同时满足条件α1＞0、α2＞0、α1+α2＜1、s＞0时，认为r0被面元遮挡，否则r0不被面元遮挡。
再考虑系数矩阵奇异的情形。当奇异，即行列式为零时，表示三个矢量在同一平面上，此时入射方向与三角形面元ABC平行，这种情况下r0也不会被面元遮挡。
图4为面元遮挡判断方法流程图。经过遮挡判断之后，剔除了被遮挡的亮面元，这些面元不参与形成射线管。
步骤4：设定垂直于入射方向的投影面
在这个步骤中，设定垂直于入射方向的投影面，此投影面经过坐标原点。经过步骤3余下的亮面元都投影到此投影面上，形成诸多三角形射线管，参与射线追踪计算。
设需要投影的三角形面元三个顶点位置矢量是rA,rB,rC，重心是r0，将三角形各顶点和重心在投影面上的投影位置矢量分别是rA′,rB′,rC′，r0′。以r0的投影计算为例进行说明，r0′满足的关系为 由此可以计算得三角形的顶点位置矢量在投影面上的投影同样可以计算得到至此，三角形面元在投影面上的投影计算完毕。
步骤5：计算第一次反射的RCS贡献
计算第一反射的RCS贡献时，因为参与计算的亮面元已经全部得以确认，因此不必再通过射线管与面元求交再确定反射面，直接将亮面元进行物理光学积分即可。物理光学积分计算理想导体的RCS公式为：
其中σ为RCS，k＝2π/λ为波数，λ为波长，为法向，为接收天线电场方向，为入射波磁场方向， 为散射波传播方向，r为被积面元上位置矢量，S为积分三角形面元。对于三角形面元，由于法向恒定，则RCS计算公式成为
上式面积积分可以化成沿三角形面元边的积分，结果如下：
式中Lm是第m条边的矢量，rmc是第m条边中点的位置矢量，A是面元的面积，r0是三角形面元的重心。当时，意味着散射方向是入射方向对于反射面的镜面反射方向。
步骤6：计算反射射线管及第二次反射的RCS贡献
经第一次反射后，三角形射线管会反射到其他方向。由第一次反射面的法线方向可计算得到反射方向的表达式相应的电场方向、磁场方向也可以根据电磁场边界条件进行计算，结果为其中是入射波电场方向，是反射电磁波的电场、磁场方向。将电磁波反射后新的作为入射波，就可以计算二次反射。反射后的射线管横截面仍然保持为三角形，要与各面元进行求交判断，方法同前。图5表示三角形射线管经二面角一个面反射后，照射到另一个面的示意图。
设第一次反射面重心r0沿反射线在二次反射面上交点的位置矢量为r1，三角形面元顶点经射线管在二次反射面上的交点分别为rA1,rB1,rC1。以rA1的计算为例进行说明，设反射面的法向为rA1在二次反射面上，因此满足关系rA1又是由rA沿延伸而来，rA1-rA与同一方向，因此满足关系由这两个关系，即可计算得出rA1的表达式，即其他两个顶点以及r1的计算也按照这种方法进行。
得到三角形射线管在二次反射面上的投影之后，就可以计算RCS，但是需要对物理光学计算RCS的公式(6)进行相应修改，将反射后相位考虑进去，方法是将公式(6)乘以因子再将式(6)中相应的入射波方向改为最终二次反射的RCS计算公式如下：
式中rmc1是第m条边的中点rm沿与二次反射面的交点。是三角形射线管在二次反射面上投影的中心。
步骤7：将所有RCS贡献叠加，得到目标总和的RCS
将所有照射到的面元RCS进行叠加，得到总的RCS为其中N为所有面元总的反射次数，表示每次反射的RCS贡献。
以二面角的RCS计算为例说明本方法。二面角边长1m，三角形面元划分和角度规定如图5所示，入射电磁波频率为3GHz，方位角从0～90°，RCS计算结果如图6(a)、(b)所示。图6(a)表示HH(水平-水平)极化电磁波的RCS，入射波电场方向与地面平行，图6(b)表示VV(垂直-垂直)极化电磁波的RCS，入射波电场方向与地面垂直。可见，采用三角形射线管计算二次反射精度能够保证。
综上所述，本发明提出一种基于面元投影的三角形截面射线管电磁射线追踪方法，应用于电磁场仿真分析领域。基于目标表面划分的面元形成三角形截面射线管，计算的效率、精度更高，可用于分析大尺寸目标的雷达散射截面(RCS)。首先将目标表面划分成许多三角形面元，近似目标的表面形状，然后将三角形面元沿入射方向投影，形成三角形截面的射线管。根据几何光学假设，电磁波在三角形射线管中传播，三角形的射线管经反射之后，截面仍旧是三角形，反射后的射线管再照射到其他面元上，形成多次反射。以三角形面元的中心作为相位变化的参考点，用于计算相位的影响。该方法将射线管与目标面元加以关联，射线管恰好覆盖了目标的横截面，消除了传统矩形射线管与三角形面元相交形成的误差。
Claims (3)
1.一种基于面元投影的三角形截面射线管电磁射线追踪方法，其特征在于：具体包括以下步骤：
步骤1：将要计算的目标表面划分成许多三角形面元来近似；
所述的“划分成许多三角形面元”，是指目标形状用封闭的多面体表示，每个面都是三角形；所述的“近似”，是指三角形面元边长小于目标表面曲率半径的1/8，同时也要小于入射电磁波波长的1/5，这时能够较好的体现目标形状；
步骤2：判断三角形面元的朝向，面元法向朝向入射电磁波方向的面元称为亮面元，背向入射电磁波朝向的面元称为暗面元；
步骤3：所有的亮面元之间进行遮挡判断，剔除被遮挡的亮面元；
每两个亮面元之间都要通过计算看是否相对于入射电磁波有遮挡关系，如一个面元的重心被另一个面元遮挡，则电磁波不能照射到该重心；r0是待判断遮挡的点；设L1＝AB，L2＝AC，三角形ABC三个顶点的位置矢量分别为rA、rB、rC，则三角形所在平面上任一点r能表示为r＝rA+α1L1+α2L2；如果同时满足条件α1＞0、α2＞0、α1+α2＜1，则r在三角形内，否则r在三角形外；
由于r是是经过r0的入射射线上的一点，r也能用r0、入射方向表示，即其中s表示从r到r0的位移；s＞0时，r沿入射方向在r0的后面，此时三角形面元可能遮挡r0；s＜0时，r沿入射方向在r0的前面，此时三角形面元不可能遮挡r0；s＝0时，为临界情况，r0在三角形面元ABC上；在遮挡判断过程中，临界情况认为面元遮挡无效；
根据r的两种表达式，能得出如下关系：进一步写为线性方程组的形式：
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mover>
<mi>i</mi>
<mo>^</mo>
</mover>
<mo>&rsqb;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>s</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>A</mi>
</msub>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中各矢量都是列向量，因此系数矩阵是一个方阵，在矩阵非奇异情况下能求得三个未知数α1、α2、s；当同时满足条件α1＞0、α2＞0、α1+α2＜1、s＞0时，认为r0被面元遮挡，否则r0不被面元遮挡；
再考虑系数矩阵奇异的情形，当奇异，即行列式为零时，表示三个矢量L1,L2,在同一平面上，此时入射方向与三角形面元ABC平行，这种情况下r0也不会被面元遮挡；步骤3中所述的“剔除被遮挡的亮面元”，是指如果面元的重心被遮挡，则其不参与物理光学积分计算，对电磁散射没有贡献；
步骤4：根据电磁波入射方向设定一个经过原点、与入射方向垂直的平面，称为投影面；
步骤5：将经过步骤3剔除后的亮面元向投影面投影，投影也是三角形形状，这些投影三角形恰好覆盖了目标的投影轮廓，以这些投影三角形为横截面，形成沿电磁波入射方向的射线管；
投影面是过坐标原点并与入射方向垂直的平面，各面元沿入射方向延伸，与投影面的交点即是投影；设需要投影的三角形面元三个顶点位置矢量是rA,rB,rC，重心是r0，将三角形各顶点和重心在投影面上的投影位置矢量分别是r′A,r′B,r′C，r′0；以r0的投影计算为例进行说明，r′0满足的关系为由此计算得三角形的顶点位置矢量在投影面上的投影同样能计算得到
步骤6：判断每条射线管与目标表面各个面元的求交情况，计算反射射线，进而计算每次反射所形成的RCS，其中第一次反射面就是形成射线管的亮面元，无需进行求交判断；
计算第一反射的RCS贡献时，因为参与计算的亮面元已经全部得以确认，因此不必再通过射线管与面元求交再确定反射面，直接将亮面元进行物理光学积分；物理光学积分计算理想导体的RCS公式为：
<mrow>
<msqrt>
<mi>&sigma;</mi>
</msqrt>
<mo>=</mo>
<mfrac>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<munder>
<mo>&Integral;</mo>
<mi>S</mi>
</munder>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>e</mi>
<mo>^</mo>
</mover>
<mi>s</mi>
</msub>
<mo>&times;</mo>
<msub>
<mover>
<mi>h</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mi>w</mi>
<mo>&CenterDot;</mo>
<mi>r</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>S</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中σ为RCS，k＝2π/λ为波数，λ为波长，为法向，为接收天线电场方向，为入射波磁场方向， 为散射波传播方向，r为被积面元上位置矢量，S为积分三角形面元；对于三角形面元，由于法向恒定，则RCS计算公式成为
<mrow>
<msqrt>
<mi>&sigma;</mi>
</msqrt>
<mo>=</mo>
<mfrac>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>e</mi>
<mo>^</mo>
</mover>
<mi>s</mi>
</msub>
<mo>&times;</mo>
<msub>
<mover>
<mi>h</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<munder>
<mo>&Integral;</mo>
<mi>S</mi>
</munder>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mi>w</mi>
<mo>&CenterDot;</mo>
<mi>r</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>S</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
上式面积积分可以化成沿三角形面元边的积分，结果如下：
<mrow>
<msqrt>
<mi>&sigma;</mi>
</msqrt>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>^</mo>
</mover>
<mi>s</mi>
</msub>
<mo>&times;</mo>
<msub>
<mover>
<mi>h</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&times;</mo>
<mi>w</mi>
</mrow>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mo>|</mo>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&times;</mo>
<mi>w</mi>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<msub>
<mi>L</mi>
<mi>m</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mi>w</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</msup>
<mi>sin</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>k</mi>
<mi>w</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>L</mi>
<mi>m</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&times;</mo>
<mi>w</mi>
<mo>|</mo>
<mo>&NotEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>^</mo>
</mover>
<mi>s</mi>
</msub>
<mo>&times;</mo>
<msub>
<mover>
<mi>h</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mi>w</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
</mrow>
</msup>
<mi>A</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&times;</mo>
<mi>w</mi>
<mo>|</mo>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中Lm是第m条边的矢量，rmc是第m条边中点的位置矢量，A是面元的面积，r0是三角形面元的重心；当时，意味着散射方向是入射方向对于反射面的镜面反射方向；
经第一次反射后，三角形射线管会反射到其他方向，由第一次反射面的法线方向可计算得到反射方向的表达式相应的电场方向、磁场方向也可以根据电磁场边界条件进行计算，结果为其中是入射波电场方向，是反射电磁波的电场、磁场方向；将电磁波反射后新的作为入射波，就可以计算二次反射；反射后的射线管横截面仍然保持为三角形，要与各面元进行求交判断，方法同前；
设第一次反射面重心r0沿反射线在二次反射面上交点的位置矢量为r1，三角形面元顶点经射线管在二次反射面上的交点分别为rA1,rB1,rC1；以rA1的计算为例进行说明，设反射面的法向为rA1在二次反射面上，因此满足关系rA1又是由rA沿延伸而来，rA1-rA与同一方向，因此满足关系由这两个关系，就能计算得出rA1的表达式，即剩下两个顶点以及r1的计算也按照这种方法进行；
得到三角形射线管在二次反射面上的投影之后，就可以计算RCS，但是需要对物理光学计算RCS的公式(6)进行相应修改，将反射后相位考虑进去，方法是将公式(6)乘以因子再将式(6)中相应的入射波方向改为最终二次反射的RCS计算公式如下：
<mrow>
<msqrt>
<mi>&sigma;</mi>
</msqrt>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<msub>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>^</mo>
</mover>
<mi>s</mi>
</msub>
<mo>&times;</mo>
<msub>
<mover>
<mi>h</mi>
<mo>^</mo>
</mover>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mi>w</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
</mrow>
</msup>
<mfrac>
<mrow>
<msub>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mo>|</mo>
<msub>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<msub>
<mi>L</mi>
<mi>m</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>jkw</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mi>c</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msup>
<mi>sin</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>kw</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>L</mi>
<mi>m</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<msub>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>&NotEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<msub>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>^</mo>
</mover>
<mi>s</mi>
</msub>
<mo>&times;</mo>
<msub>
<mover>
<mi>h</mi>
<mo>^</mo>
</mover>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<mover>
<mi>i</mi>
<mo>^</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
</mrow>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
<mi>w</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
<mi>A</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<msub>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式中rmc1是第m条边的中点rm沿与二次反射面的交点；
步骤7：将所有射线管的所有反射形成的RCS进行叠加，即得到目标总的RCS。
2.根据权利要求1所述的一种基于面元投影的三角形截面射线管电磁射线追踪方法，其特征在于：在步骤2中所述的“判断三角形面元的朝向”，其说明如下：看单独一个面元是否能被入射波照射，如能照射，说明其法向与入射波方向同向，点积小于0，即其中为面元法向，为入射波方向，此面元归类为亮面元；如不能照射，说明其法向与入射波方向反向，点积大于0，即此面元归类为暗面元；对的临界情形，按照暗面元处理。
3.根据权利要求1所述的一种基于面元投影的三角形截面射线管电磁射线追踪方法，其特征在于：在步骤7中所述的“进行叠加”，其说明如下：所有照射到的面元RCS按照相位求和，得到总的RCS为其中N为所有面元总的反射次数，表示每次反射的RCS贡献。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510017244.5A CN104573368B (zh) | 2015-01-13 | 2015-01-13 | 基于面元投影的三角形截面射线管电磁射线追踪方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510017244.5A CN104573368B (zh) | 2015-01-13 | 2015-01-13 | 基于面元投影的三角形截面射线管电磁射线追踪方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104573368A CN104573368A (zh) | 2015-04-29 |
| CN104573368B true CN104573368B (zh) | 2018-02-23 |
Family
ID=53089414
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510017244.5A Active CN104573368B (zh) | 2015-01-13 | 2015-01-13 | 基于面元投影的三角形截面射线管电磁射线追踪方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104573368B (zh) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105405144B (zh) * | 2015-11-16 | 2017-11-10 | 徐灿 | 基于OpenGL双重拾取的空间目标OCS快速计算方法 |
| CN105403881A (zh) * | 2015-12-14 | 2016-03-16 | 北京环境特性研究所 | 基于大面元模型的引信多普勒频率获取方法 |
| CN105953925A (zh) * | 2016-06-03 | 2016-09-21 | 华中科技大学 | 一种基于快速射线追踪的毫米波辐射亮温获取方法 |
| CN106355641B (zh) * | 2016-08-26 | 2019-07-23 | 西安电子科技大学 | 三维粗糙目标雷达散射截面的仿真方法 |
| CN106556833B (zh) * | 2016-11-24 | 2019-01-08 | 上海无线电设备研究所 | 基于时域弹跳射线法快速近场计算的isar成像仿真方法 |
| CN108286968B (zh) * | 2017-01-10 | 2021-11-09 | 中兴通讯股份有限公司 | 一种反向射线追踪方法和装置 |
| CN107300645B (zh) * | 2017-06-06 | 2019-07-23 | 华中科技大学 | 一种快速的射线追踪方法及系统 |
| CN108446430B (zh) * | 2018-02-05 | 2021-08-06 | 西安电子科技大学 | 基于投影法的高频电磁遮挡判断方法 |
| CN110580742A (zh) * | 2018-06-07 | 2019-12-17 | 南京理工大学 | 基于gpu并行的sbr实现目标电磁散射特性建模与分析的方法 |
| CN109663358A (zh) * | 2018-12-21 | 2019-04-23 | 成都四方伟业软件股份有限公司 | 动态场景剔除遮挡物体的方法、装置及电子设备 |
| CN110705058B (zh) * | 2019-09-19 | 2021-08-31 | 电子科技大学 | 一种针对超电大尺寸规模目标的近场电磁散射仿真方法 |
| GB2618213A (en) * | 2021-03-23 | 2023-11-01 | Imagination Tech Ltd | Intersection testing in a ray tracing system |
| GB2599184B (en) | 2021-03-23 | 2022-11-23 | Imagination Tech Ltd | Intersection testing in a ray tracing system |
| CN113030900B (zh) * | 2021-03-26 | 2022-08-09 | 中国人民解放军国防科技大学 | 基于面元分布的动态匹配反射系数缩比测量方法和装置 |
| CN114387594B (zh) * | 2022-01-17 | 2024-04-19 | 中国人民解放军63660部队 | 一种高频电磁散射遮挡识别方法 |
| CN118135228A (zh) * | 2024-04-30 | 2024-06-04 | 中国人民解放军战略支援部队航天工程大学 | 基于面元法向和射线三角形相交算法的电磁遮挡判断方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103439698A (zh) * | 2013-08-27 | 2013-12-11 | 西北工业大学 | 获取雷达散射面积的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7646332B2 (en) * | 2007-08-09 | 2010-01-12 | Raytheon Company | Method and apparatus for interleaved gridding in distributed multiple computing for real-time RCS prediction |
-
2015
- 2015-01-13 CN CN201510017244.5A patent/CN104573368B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103439698A (zh) * | 2013-08-27 | 2013-12-11 | 西北工业大学 | 获取雷达散射面积的方法 |
Non-Patent Citations (4)
| Title |
|---|
| 基于通用CAD几何模型的复杂目标RCS计算;李向军 等;《探测与控制学报》;20081231;第30卷(第6期);第46-50页 * |
| 基于面元边缘法的直升机RCS计算与分析;蒋相闻 等;《南京航空航天大学学报》;20110630;第43卷(第3期);第430-434页 * |
| 改进的物理光学技术计算目标的雷达散射截面;彭兰;《中国优秀硕士学位论文全文数据库 基础科学辑》;20131215(第S2期);第15-19页 * |
| 雷达散射截面计算方法分析与研究;刘英;《中国优秀硕士学位论文全文数据库 信息科技辑》;20020615(第1期);第34页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104573368A (zh) | 2015-04-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104573368B (zh) | 基于面元投影的三角形截面射线管电磁射线追踪方法 | |
| CN101216556B (zh) | 电大复杂体目标与粗糙面背景复合电磁散射数值仿真方法 | |
| EP2887092A1 (en) | Computing radar cross section | |
| CN106772301B (zh) | 一种基于介质射线追踪的多层非平行界面介质电磁散射仿真算法 | |
| CN104200074B (zh) | 快速获取目标电磁散射特性的多层复波束方法 | |
| CN104865562B (zh) | 基于混合模型的雷达非合作目标的识别方法 | |
| CN103487789B (zh) | 应用于复杂目标多次散射的雷达散射截面快速获取方法 | |
| CN110907901A (zh) | 一种基于tdsbr的海战无源干扰计算及评估方法 | |
| CN107315881A (zh) | 用于电磁散射仿真模型的半空间格林函数与射线追踪方法 | |
| CN105953925A (zh) | 一种基于快速射线追踪的毫米波辐射亮温获取方法 | |
| CN104573257B (zh) | 基于面元空间分集的电磁遮挡判断方法 | |
| Fang | The simulation of quantum radar scattering for 3D cylindrical targets | |
| CN103217675A (zh) | 多个不共轴旋转对称体电磁散射特性的仿真方法 | |
| CN109696672B (zh) | 一种基于空间结构关联性的高分辨率穿墙雷达成像方法 | |
| He et al. | A vector meshless parabolic equation method for three-dimensional electromagnetic scatterings | |
| CN103295188B (zh) | X射线蒙特卡洛模拟的路径积分方法 | |
| CN109125952A (zh) | 基于核模型的卷积叠加剂量计算方法 | |
| CN103439698A (zh) | 获取雷达散射面积的方法 | |
| CN107300645A (zh) | 一种快速的射线追踪方法及系统 | |
| Gao et al. | Fast RCS prediction using multiresolution shooting and bouncing ray method on the GPU | |
| Ye et al. | A hybrid KA-MoM algorithm for computation of scattering from a 3-D PEC target above a dielectric rough surface | |
| CN103853914B (zh) | 一种快速预估水下圆形角反射体散射声场的方法 | |
| CN110083904A (zh) | 基于gpu加速的量子雷达散射截面计算方法 | |
| CN111339606A (zh) | 一种基于绕射原理的机翼遮挡效应的计算方法 | |
| Gao et al. | Efficient GPU implementation of SBR for fast computation of composite scattering from electrically large target over a randomly rough surface |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180716 Address after: 100041 room 2961, 3 building, 3 West well road, Badachu high tech park, Shijingshan District, Beijing. Patentee after: Beijing boo Rui Technology Co., Ltd. Address before: 100191 Xueyuan Road, Haidian District, Beijing, No. 37 Patentee before: Beihang University |