专利名称：一种粒子发射模拟的电磁求解方法
技术领域：
本发明涉及一种电磁求解方法，特别是一种通过模拟发射电磁粒子进行电磁分布态势计算的方法，用于大范围区域的电磁计算求解及电磁态势可视化。
背景技术：
不论是在军事领域，还是民用领域，电磁都发挥着重大的作用，对于大范围区域， 利用专门的设备进行接收或探测将是一个浩大的不可能实现的工程，因此通过计算的手段将电磁的分布态势计算出来具有重大意义。射线跟踪方法(Ray Tracing)是一种比较常见的电磁计算方法，射线跟踪方法分为正向跟踪与反向跟踪两种。正向跟踪是指射线跟踪由发射源出发，向周围空间等角度间隔发射出大量的射线，然后跟踪每一条射线的直射、反射和绕射，并在待测点处用接收球判定到达的射线，将这些射线进行叠加，从而得到该点处的电磁场强度。但是接收球的定义可能会带来较大的误差，定义小了会漏掉应该统计的射线，定义大了会造成过大估计，这样造成了正向算法预测的可靠性低。反向跟踪与正向跟踪相似，区别在于反向跟踪从待测点发出射线，然后测得到达发射源的射线，根据传播路径的可逆性得知发射源到达待测点的射线。反向跟踪算法比正向跟踪算法精确度高，可是正向算法根据一次发射的射线束就可计算出整个空间中的态势分布，反向算法一次只能算出一点的情况，因此反向跟踪算法计算量较大。在三维可视化中，常使用粒子方法来模拟一些云、雾、火焰等模糊现象，通过大量粒子的运动变化来实现逼真的视觉效果。本文结合粒子方法对正向射线跟踪方法进行改进，构建电磁粒子模型，计算电磁在空间的分布态势。
发明内容
本发明的技术解决问题是克服现有技术中由于定义接受球导致的预测误差大的不足，提供一种粒子发射模拟的电磁求解方法，该方法通过模拟发射电磁粒子进行电磁分布态势计算。不仅在进行电磁态势可视化时，可以根据粒子的位置属性直接对其进行绘制渲染得出电磁分布态势，而且粒子的运动与电磁空间衰减的扩散一致，粒子运动导致的疏密程度的变化本身就表征了扩散导致的电磁场强度的变化，因此省去了电磁扩散因子的计笪
o本发明的技术解决方案为一种粒子发射模拟的电磁求解方法，其特征在于通过模拟发射电磁粒子进行电磁分布态势计算，其具体步骤如下步骤a :从发射源处发射粒子，初始化每个粒子的初始状态，将发射的新粒子加入到粒子列表中；步骤b :更新粒子列表中粒子的下一时刻的属性；步骤c :判断粒子是否位于地面以下，若位于地面下，说明应该经过反射或绕射， 根据粒子是否位于地形尖劈边缘附近判断是绕射还是反射；
步骤d :若粒子应该进行反射，根据电磁波的反射定律更正此粒子的实际位置，并根据反射系数来求得粒子的反射消亡概率，通过此概率判定粒子是否消亡；若粒子应该进行绕射，同反射类似，根据电磁波的绕射定律更正此粒子的实际位置，并根据绕射系数来求得粒子的绕射消亡概率，通过此概率判定粒子是否消亡；步骤e :判断粒子是否位于边界区域外，剔除边界外粒子，剔除反射和绕射中消亡的粒子；步骤f:若列表中的粒子数目未到达稳定，即计算的空间区域内粒子还未达到平衡状态，重复步骤a至e ;步骤g:若列表中的粒子数目到达稳定，将空间区域进行网格划分，统计每个网格内的粒子数目，并输出到文件中。本发明的原理是通过模拟发射的电磁粒子进行电磁分布态势计算，用场源处间隔发射的粒子来模拟电磁传播，粒子的疏密程度表示电磁能量，粒子的反射、绕射遵循电磁波的反射与绕射定律，根据反射系数和绕射系数计算粒子的消亡概率，删除消亡的粒子，使粒子在空间的疏密分布符合电磁场强度的分布规律。本发明与现有技术相比的优点在于本发明使用粒子代替射线来跟踪算法中的射线，每个粒子代表的能量是一定的，粒子运动与空间衰减的扩散一致，因此粒子运动本身就表征了扩散导致的电磁场强度的变化。现有的正向射线跟踪算法采用接收球的方法，定义接收球半径，统计叠加进入接收球的射线，接收球的定义可能会带来较大的误差。假设无限大的自由空间中，只考虑直射的情况，接收球半径等于两条射线在该点附近的间距，一个接收球只能接收到一条射线，若接收球半径略大于此间距，可能出现一个接收球接收到两条射线的情况，若略小于此间距，又可能出现某接收球未接收到射线的情况，误差会很大。如果宏观观察电磁态势，可直接使用粒子系统，以粒子进行渲染模拟实现可视化；如果需要精确求得电磁能量，可根据需要的分辨率对整个空间进行划分，统计每个网格内的粒子数量， 除以该网格的体积，即可得出该点的电磁坡印廷矢量平均值。而且此方法计算可并行度很高，可以大大提高计算速度，并且特别适用于在计算同时绘制渲染，进行电磁态势可视化观察。
图I为本发明使用的粒子发射模拟的电磁求解方法流程图；图2为粒子运动时根据粒子间距离求解粒子速度的示意图。
具体实施例方式如图I所示，本发明的具体实施方法如下I、从发射源发射粒子，将发射的新粒子加入到粒子列表中。粒子在各个方向的发射密度主要取决于天线的方向图因子/⑷,的，6和0为球坐标系中点的角度坐标分量。以发射源为顶点，向各个方向进行空间等间隔切分，给定0的取样间隔A 0，炉的取样间隔Ap，一般来说AS = 。切分后的各角度为9 j = i A 9 i = 1,2, 3... int (Ji / A 9 )(I)
权利要求
1.一种粒子发射模拟的电磁求解方法，其特征在于包括以下步骤步骤a :从发射源处发射粒子，初始化每个粒子的初始状态，将发射的新粒子加入到粒子列表中；步骤b :更新粒子列表中粒子的下一时刻的属性；步骤c :判断粒子是否位于地面以下，若位于地面下，说明应该经过反射或绕射，根据粒子是否位于地形尖劈边缘附近判断是绕射还是反射；步骤d :若粒子应该进行反射，根据电磁波的反射定律更正此粒子的实际位置，并根据反射系数来求得粒子的反射消亡概率，通过此概率判定粒子是否消亡；若粒子应该进行绕射，同反射类似，根据电磁波的绕射定律更正此粒子的实际位置，并根据绕射系数来求得粒子的绕射消亡概率，通过此概率判定粒子是否消亡；步骤e :判断粒子是否位于边界区域外，剔除边界外粒子，剔除反射和绕射中消亡的粒子;步骤f :若列表中的粒子数目未到达稳定，即计算的空间区域内粒子还未达到平衡状态,重复步骤a至e ；步骤g:若列表中的粒子数目到达稳定，将空间区域进行网格划分，统计每个网格内的粒子数目，并输出到文件中。
2.根据权利要求I所述的粒子发射模拟的电磁求解方法，其特征在于所述的步骤a 中，粒子发射的初始位置定义为Am)，其中R0为定义的粒子发射的初始位置与场源的距离，为定值，假定场源位于原点+亨 P = W-TiiCDu年*⑵2 Tlj为以发射源为顶点向各个方向进行空间等间隔切分后的角度，给定e的取样间隔A 0，炉的取样间隔A炉，其中A沒=A炉9j = i A 9 i = l,2,3... int (Ji / A 9 )(3)(Pj = jAg) 7 = 1,2,3 int(2^- / Ag))(4)i与j均为正整数，int (Ji/A 0 )表示对Ji / A 0取整数；IVij 为定义在{矽，炉)|#+ A0/2],(p Gicpj -Acp/2,(Pj +A炉/2]}范围内一次发射的粒子数；叫与rij均为正整数，riin-表示Iii与rij的乘积；n「-歷 N'(5)DUj=m^N2(6)N1N2为给定的此范围内一次发射的粒子数最大值，N1N2为N1与N2的乘积，N1与N2均为给定常数，若M = #，则N1 = N2，并且Iii = n」，/队妁为球坐标系中随0和小角度变化的函数，被称为天线的方向图因子，/max@，妁为函数最大值。
3.根据权利要求I所述的粒子发射模拟的电磁求解方法，其特征在于所述的步骤b中，对粒子的速度!>更新具体模型如下离散的粒子运动速度为
4.根据权利要求I所述的粒子发射模拟的电磁求解方法，其特征在于所述的步骤d 中，反射或绕射后一部分粒子产生消亡，等价于发生反射或绕射存活下来的粒子的相邻粒子发射间隔角度发生改变，更改后的粒子发射间隔S'为5' = S R(9)或6r = 8D(10)6为更改前相邻粒子发射间隔角度，R为电磁波的反射系数，D为绕射系数。
全文摘要
本发明涉及一种粒子发射模拟的电磁求解方法，将电磁传播用粒子的运动来模拟，粒子的疏密程度表示电磁能量。每个粒子具有位置、速度、移动总距离、间隔角度、寿命、是否死亡等属性，粒子遵循电磁的反射与绕射定律。这样，对粒子绘制渲染就可实现可视化；也可根据需要的分辨率对空间进行划分，统计网格内的粒子数量，求解出该点的电磁能量；另外粒子运动表征了辐射的扩散，省去了扩散因子的计算。本发明具有精度高、方法简便直观、可并行度高的特点，可应用于各种电磁计算求解及可视化系统。
文档编号G01R29/08GK102608437SQ20121010501
公开日2012年7月25日 申请日期2012年4月11日 优先权日2012年4月11日
发明者任磊, 张霖, 穆兰, 陶飞 申请人:北京航空航天大学