专利名称：一种优化综合孔径阵构型的方法
技术领域：
本发明涉及一种对综合孔径阵构型进行优化的方法，具体地，本发明用来确定天文观测的天线阵中各单元天线的位置。
背景技术：
在天文观测中，成像分辨率是衡量望远镜性能的一个重要参数。对单孔径望远镜来说，其成像分辨率为λ/D，其中λ为观测波长，D为望远镜的口径。由于工程技术上的困难，单架望远镜的口径不能做得太大。例如目前最大的光学望远镜Keck的口径为10米，对应的角分辨率约为0. 1角秒。综合孔径成像技术可以克服以上限制。对综合孔径阵来说，其成像分辨率公式中的D为阵中望远镜之间的最大基线，进而可以通过延长阵列的基线来提高成像分辨率。综合孔径技术最早被应用在射电天文领域(Rylel952)。经过50多年的不断发展与完善，当今的射电综合孔径阵已经延伸到太空，其观测角分辨率可以超过0. 1毫角秒。该技术的创始人Martin Ryle也因此获得1974的诺贝尔物理奖。近年来该技术被逐步推广到光学波段(Baldwin 1996)，一些最大型的光学望远镜如两架口径为10米的Keck和四架口径为8米的VLT，都致力于在几年内使得光学综合孔径技术达到能够进行实际天文观测的水平。在可以预见的未来，综合孔径成像技术也将应用到X波段(Cash 2000)。NASA 预研的MAXIM卫星，正是要做到这一点。由于X射线的波长更短，MAXIM的角分辨率可以达到微角秒，足够用来探测一些临近星系中央的大质量黑洞。综合孔径阵的构型，也就是阵中望远镜之间的相对位置，直接影响到成像质量。综合孔径阵测量的是天体图像对应的空间谱，通过i^ourier变换把谱转化为天体图像。一个糟糕的构型会重复测量一些谱成份，同时又把另外一部分谱成份给弄丢了，相应的图像也会丢失部份信息。另外，如果构型是各向异性的，那么在望远镜数目比较少方向的成像分辨率会降低，还会带来较大的旁瓣，形成虚假的图像结构。因此，一个设计良好的综合孔径阵构型往往需要最大限度地获取图像的谱信息。这个要求又和我们所研究天体的类型紧密联系在一起如果我们所探测的天体是致密源，它们的空间谱成份主要集中在高频部分，相应的阵构型应该在长基线处(对应高频空间谱)布置更多的望远镜；反之，对于延展源，就需要在短基线处把望远镜布置得密一些。一个综合孔径阵的构型是否优化，还受到一些外部条件的限制。例如，一个望远镜分布在东西走向的利用地球自转来完成孔径综合的阵，对于高纬度的源会有比较好的观测效果，可以算是一个构型优化的阵；然而对于低纬度源来说， 该阵的观测效果会很差，不能算作构型优化的阵。由于综合孔径阵构型优化设计的复杂性，目前尚未有完美的解决方案。为了简化问题，许多优化设计方案一般只考虑其中几项重要指标，如u-v分布的均勻性、旁瓣的水平等等。澳大利亚的Bimton为澳大利亚SKA概念提出了由分布紧凑的多个单元和周边稀疏分布的非对称和旋转结构(Bunton 1999,2000)；为350个单元的Allen Telescope Array (ΑΤΑ), Borck初步考虑了面包圈和圆内均勻分布的干涉阵结构(Borck 2000,2001)； Cornwell考虑了小数目的干涉阵，得到对称花纹的u_v分布，如水晶一般非常漂亮，称为水晶阵列(Cornwell 1986)；为 ALMA(Atacama Large Millimeter/submillimeter Array)计划，Conway分析了和Bimton近似的有旋臂的旋转干涉阵结构，认为旋转干涉阵结构的u_v 分布更趋向于高斯型分布(Conway 2000a, 2000b, 2000c) ；Woody分析了位于圆周上对称的干涉阵(Woody 1999)。这些方案主要是通过个人的经验获得，缺乏严格的数学论证。除了上面的一些经验方案外，人们还就简化的综合孔径阵构型设计要求(如不考虑地球自转和观测源的位置等等)提出了一些优化方案。第一种方案是Keto在1997年为 SMA(Sub-Millimeter Array)项目提出的模拟退火法(Keto 1997)。该方法从一个随机的初始阵列构型出发，运用模拟退火法来搜索更加好的构型。当一个稳定构型出现时，用一个随机的，与之接近的新构型开始新的搜索。这样做可以跳出一些局部优化的解，使得解达到全局最优化。由于综合孔径阵构型的复杂性，在实际的应用时，该方法很难找到全局最优化的构型。而且该算法的计算量大，只适合望远镜数目比较少的阵构型的优化。第二种方法是Kogan为ALMA计划提出的最小旁瓣法(Kogan 1997)。该方法基于方向瓣与天线位置之间的解析关系，通过移动天线的位置来最小化综合孔径阵方向图旁瓣。由于这个方法只是从旁瓣角度出发，并且只适应于快照模式，其应用范围受到一定限制。另外，通过最小旁瓣法也不能直接得到特定的均勻或者呈高斯型的u-v分布。第三种方法是Boone为ALMA和 ATA提出的梯度压力法(Boone 2001)。该方法定义了一种虚拟的“压力”，该压力正比与一块u-v区域内实际u-v点数与理想u-v点数的差值。如果u-v点数过多，就对应一个正压力，该压力会驱动这些u-v点向邻近空旷的区域扩散。最终，实际的u-v点分布会很接近我们所想要的分布，而对应的阵构型就是优化构型。这种方法的灵活性要大一些，应用范围也相对广泛。它的不足之处在于该方案中的力和u-v点的数密度联系在一起，其大小和方向的精度都比较差，而且还随着计数区域的不同划分方法而变化，所以最终优化结果的可靠性和真实性较低。而且并不适合于望远镜数目少的情况，而对于望远镜数比较多的阵，其构型会相当复杂，我们很难确认是否已经获得最优化的构型。申请者2003年提出一种系统的综合孔径阵构型的优化设计算法，即筛选算法(Su et al. 2004)。它的基本思想是候选望远镜位置对应的u_v分布比最终综合孔径阵对应的 u-v分布密集的多，根据每个u-v点对u-v分布均勻性的贡献进行逐步筛选，使得最后剩下的u-v点符合我们所期望的分布，最终获得优化构型。筛选算法不但可以获得均勻u-v分布的优化阵构型，通过增加额外的权重因子，还可以获得Gauss型u-v分布的优化阵构型。 与其它方法不同的是，筛选算法可对以地球自转观测模式为主的综合孔径阵构型进行优化设计，还可针对不同纬度的观测源，进行综合孔径阵构型的折衷优化。因此，筛选算法是一种应用广泛的综合孔径阵构型的优化设计方法。目前，在下一代大型综合孔径阵如SKA (Square Kilometer Array)的预研究中，阵构型的优化设计是一个极为重要的部分，为此都成立了专门的工作组。该发明提出一种对综合孔径阵的构型进行优化的方法，基于严格的物理意义，经过完整的数学公式推导，可以得到优化的天线阵结构。
发明内容
针对目前国内外现有方法不普遍适用的缺点，本发明基于严格的物理意义，经过完整的数学公式推导，给出获得优化的天线阵结构的解决方案，本解决方案在不限制天线阵的数目的前提下，能够得到多天线组阵优化的准确结果。为了实现上述目的，本发明提出的一种优化综合孔径阵构型的方法，其特征在于， 所述综合孔径阵具有N个望远镜单元，所述N个望远镜单元具有N* (N-I)个U-V点，其中， 每个望远镜单元对应2* (N-I)个u-v点，该方法包括以下步骤为所述具有N个望远镜单元的综合孔径阵随机设置一初始构型；使用一个带电粒子来表示所述N* (N-I)个u-v点中的每一个；计算得到每个u-v点所受到的库伦力Fm，其中Fm表示第i个望远镜单元与第j 个望远镜单元对应的u-v点所受的库伦力；根据所述每个u-v点所受到的库伦力Fm计算得到每一个望远镜单元所受到的库伦合力Fi ；为每一个望远镜单元添加一耗散阻力& ；计算优化评价指标，所述优化评价指标达到某一条件时的构型即为优化的均勻分
布构型。本发明所提出的对综合孔径阵构型进行优化的方法的优点在于·首次将库伦场引入到综合孔径阵构型的优化设计中。库伦力的作用促使阵构型向优化的方向演化。通过设置额外的耗散阻力，使得阵构型停留在最优化的位置上。库伦场的引入使得阵构型的优化设计更为精确，并大大降低优化的计算量。·首次将势能当作评价综合孔径阵构型优劣的指标，从而将判断阵构型是否优化建立在合理的物理意义基础之上，能够将现有的一些优化构型进行比较，并进一步研究优化构型的规律性。·将最新的综合孔径阵构型的优化研究与标准软件包开发相结合。开发基于库伦场的构型优化设计软件包。这些软件包将在LOFAR(LOw Frequency Array for Radio astronomy)、SKA (Square Kilometer Array)等大型射电综合孔径阵构型的设计中发挥重要作用，对未来光学、X-ray波段等的综合孔径阵构型的优化设计具有参考和使用价值。
图1为本发明优化综合孔径阵构型的方法的流程图。图2为5单元综合孔径阵的望远镜位置分布及其u-v点分布的示意图。图3为基于库伦场的综合孔径阵构型优化设计的模拟结果示意图。
具体实施例方式为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。参考图1，本发明所提出的优化综合孔径阵构型的方法的具体步骤如下步骤1，在对综合孔径阵构型进行优化之前，先随机给定N个望远镜单元一个初始的阵构型。步骤2，使用一个带电粒子来表示每一个u-v点。在本发明的技术方案中，我们首次将电磁学中库伦场的概念引入综合孔径阵构型的优化设计中。假设一个综合孔径阵有N个望远镜单元，相应的u-v点有N* (N-I)个，每个望远镜单元对应2*(N-I)个u-v点。以具有5个望远镜单元的综合孔径阵为例，图2显示了该综合孔径阵望远镜的位置信息(左图)和相对应的u-v点分布(右图)。其中一个望远镜用粗线方块表示，其余的望远镜用细线方块表示；右图中共有20个u-v点，每个u-v点对应两个望远镜，u-v点的位置为这两个望远镜的位置差。与左图粗线方块望远镜相关的 u-v点(共8个)用粗线圆圈表示，其余的u-v点用细线圆圈表示。在本技术方案中，每个 u-v点均用一个带单位正电荷的粒子表示，这些粒子由于库伦力的作用，互相排斥，但又被限制在一个固定区域(如单位圆)内，所述库伦力可表示为 Γ ρ _ Ku
权利要求
1.一种优化综合孔径阵构型的方法，其特征在于，所述综合孔径阵具有N个望远镜单元，所述N个望远镜单元具有N* (N-I)个u-v点，其中，每个望远镜单元对应2* (N-I)个u_v 点，该方法包括以下步骤为所述具有N个望远镜单元的综合孔径阵随机设置一初始构型；使用一个带电粒子来表示所述N*(N-1)个u-v点中的每一个；计算得到每个u-v点所受到的库伦力Fm，其中Fm表示第i个望远镜单元与第j个望远镜单元对应的u-v点所受的库伦力；根据所述每个u-v点所受到的库伦力Fm计算得到每一个望远镜单元所受到的库伦合 ；为每一个望远镜单元添加一耗散阻力&；计算优化评价指标，所述优化评价指标达到某一条件时的构型即为优化的均勻分布构型。
2.如权利要求1所述的优化综合孔径阵构型的方法，其特征在于，所述带电粒子为带单位正电荷的粒子。
3.如权利要求1所述的优化综合孔径阵构型的方法，其特征在于，所述带电粒子由于库伦力的作用互相排斥，但又被限制在一个固定区域内。
4.如权利要求1所述的优化综合孔径阵构型的方法，其特征在于，所述每一个望远镜单元所受到的合力作用Fi表示为
5.如权利要求1所述的优化综合孔径阵构型的方法，其特征在于，所述耗散阻力的大小正比于望远镜单元的移动速度，方向与望远镜单元的移动速度相反，即所述耗散阻力表示为其中，Vi为第i个望远镜单元的移动速度，λ为耗散系数，设为常数。
6.如权利要求1所述的优化综合孔径阵构型的方法，其特征在于，所述望远镜单元在库伦合力与耗散阻力的作用下发生移动。
7.如权利要求1所述的优化综合孔径阵构型的方法，其特征在于，所述优化评价指标为整个系统的库伦势能。
8.如权利要求7所述的优化综合孔径阵构型的方法，其特征在于，所述整个系统的库伦势能达到最小时的构型即为优化构型，所述整个系统的库伦势能为所有u-v点之间距离的倒数和，即
9.如权利要求8所述的优化综合孔径阵构型的方法，其特征在于，所述方法还能够获得Gauss型的u-v分布。
10.如权利要求9所述的优化综合孔径阵构型的方法，其特征在于，所述获得Gauss型的u-v分布的步骤为将Gauss型负能量势井Mg与库伦势能结合作为优化评价指标以获得Gauss型的u_v分布，其中，Mg表示为M = -k e—O2g g ，其中，1^和σ为调节势井形状的参数，Rj为u-v点的位置坐标。
全文摘要
本发明公开了一种优化新型综合孔径阵构型的方法。综合孔径阵成像技术是实现极高分辨率天文观测的有效手段，对研究宇宙中一些天体如活动星系核、X射线双星等起着至关重要的作用；综合孔径阵的构型直接影响着成像的质量，因此构型的优化设计对建设大型的综合孔径阵十分关键。本优化方法首次将库伦场的概念引入到优化设计中，库伦力的作用促使阵构型向优化的方向演化，通过设置额外的耗散力，最终得到最优化的阵构型。库伦场的引入使得阵构型的优化设计更为精确，并大大降低优化的计算量。另外，本优化方法首次将势能作为评价综合孔径阵构型优劣的指标，使得将阵构型的优化建立在合理的物理基础之上。
文档编号G06F17/50GK102567573SQ20111040537
公开日2012年7月11日 申请日期2011年12月8日 优先权日2011年12月8日
发明者周建峰, 苏彦 申请人:中国科学院国家天文台, 清华大学