基于轨道不平度的大型天线任意位姿时波束指向偏差计算方法
【专利摘要】本发明公开了一种基于轨道不平度的大型天线任意位姿时波束指向偏差计算方法，包括确定天线轨道采样点间距，测量天线服役时轨道各点的高度值；确定天线座架坐标系，对轨道各点高度值进行坐标变换；计算轨道各点方位坐标的拟合系数矩阵和轨道各点高度的拟合系数向量；确定轨道不平度的最佳拟合函数；计算天线处于任意方位角时所有滚轮对应的轨道高度值；计算天线在任意方位角、俯仰角时的波束指向偏差；输出波束指向偏差至天线控制系统，调整补偿天线。本发明不仅真实反映天线轨道表面，且能准确计算大型天线在任意位姿下因轨道不平度而引起的波束指向偏差，可用于定量评价轨道不平度对大型天线波束指向的影响，具有重要的学术意义和工程应用价值。
【专利说明】
基于轨道不平度的大型天线任意位姿时波束指向偏差计算 方法
技术领域
[0001] 本发明属于天线技术领域，具体是一种基于轨道不平度的大型天线任意位姿时波 束指向偏差计算方法，可用于定量评价轨道不平度对大型天线波束指向的影响，从而指导 服役时天线性能的调整与补偿，具有重要的学术意义和工程应用价值。
【背景技术】
[0002] 射电望远镜是主要接收天体射电波辐射的精密仪器，广泛应用于射电天文、雷达、 通讯和太空探测等领域。为了满足射电天文的发展需求，天线口径越造越大、结构越来越复 杂，随之天线体积与重量都不断增大，为了能更好地支撑大型天线，并使其进行预定的旋转 运动，天线座架是不可或缺的重要装置。天线座架共有三种基本形式，分别是立轴式、转台 式和轮轨式，其中轮轨式主要适用于大型天线，低速、重载工况，轮轨式天线的座架一般是 桁架结构，天线方位转动部分由滚轮和轨道支撑，这样可以省略大齿轮、大转台等构件，具 有结构简单、造价低、安装维修方便，且有较高结构精度等优点。考虑到设计、制造、运输、成 本等方面，目前国内外的大型射电望远镜多采用轮轨式座架设计，由于轨道制造与安装精 度有限，处于恶劣环境中服役的大型天线还要进行方位俯仰运动，通过轮轨来支撑如此庞 然大物，避免不了会造成轨道表面的不平度，这将严重影响着天线波束指向。
[0003] 随着天线工作频段的增高，即使是微小的轨道不平度，天线波束指向精度也变得 敏感，通常轨道不平度造成的天线波束指向偏差高达2角秒，这对高精度指向要求的天线是 不可忽视的。在已有的一些基于轨道不平度的相关研究中，例如姜正阳等的文献《考虑轨道 不平度的射电望远镜指向修正方法》中，采用了简单的三角函数来表示轨道表面，这种方法 过于粗糙，拟合误差过大，拟合残差影响了高指向精度要求的大型天线的波束指向；李永江 等的文献《天线轨道变形精密测量与指向偏差模型研究》中，指出了轨道表面的不平度可能 导致天线倾斜及挠性变形，从而降低天线指向精度，然而其建立的指向偏差模型中，指向偏 差是通过方位轴四个固定方向的倾斜误差拟合而来，拟合误差较大，分析较为片面，计算结 果不能准确反映轨道不平度对天线波束指向的影响。
[0004] 因此有必要给出轨道不平度的最佳拟合函数，并基于轨道不平度建立天线指向模 型，准确计算大型天线在任意位姿下因轨道不平度而引起的波束指向偏差，用于定量评价 轨道不平度对大型天线波束指向的影响，从而指导服役时天线性能的调整与补偿。这一过 程即为基于轨道不平度的大型天线任意位姿时波束指向偏差计算方法。
【发明内容】
[0005] 针对以前计算方法存在的不足，本发明提供了一种基于轨道不平度的大型天线任 意位姿时波束指向偏差计算方法，该方法针对大型轮轨式天线，通过最佳拟合轨道不平度， 并定量给出轨道不平度对大型天线波束指向的影响，以对服役时的天线进行调整补偿。
[0006] 为了实现上述目的，本发明提供的计算方法包括如下步骤：
[0007] (1)根据大型天线的轨道尺寸，确定轨道采样点的间距，并利用高精度测量仪器测 量大型天线在服役一段时间后的轨道各点的高度值；
[0008] (2)确定天线座架坐标系，并对测量得到的天线轨道各点高度值进行坐标变换；
[0009] (3)根据天线轨道各点坐标和高度值，计算轨道各点方位坐标的拟合系数矩阵和 轨道各点高度的拟合系数向量；
[00?0] (4)基于最小二乘原理，确定轨道不平度的最佳拟合函数；
[0011] (5)根据大型天线的方位角和轨道不平度的最佳拟合函数，计算天线处于任意方 位角时所有滚轮对应的轨道高度值；
[0012] (6)根据天线滚轮对应的轨道高度值，计算大型天线在任意方位角、俯仰角时的波 束指向偏差；
[0013] (7)输出天线波束指向偏差至天线控制系统(ACU)，对天线进行调整补偿。
[0014] 所述步骤(2)中，天线座架坐标系的确定，是以大型天线的轨道中心为原点，Z轴与 方位轴重合，垂直于地面指向天空，Y轴指向正南方向。
[0015] 所述步骤(3)计算轨道各点方位坐标的拟合系数矩阵和轨道各点高度的拟合系数 向量，包括如下过程：
[0016] (3a)假设轨道不平度的最佳拟合函数为：
[0017] f (x) =ao+aix+a2X2+. . .+akxk
[0018] 其中x表示轨道位置，f(x)表示相应的轨道高度值，k为拟合函数的最高次数，ao， ai，a2. . .ak为拟合函数中的待定系数；
[0019] (3b)根据测量得到的天线服役时轨道各点的高度值(χ^?η)···^，^)···^，!^)， 计算轨道各点方位坐标的拟合系数矩阵X:
[0021]其中，Ν为采样点个数；
[0022 ] (3c)计算轨道各点高度的拟合系数向量Η:
[0024]所述步骤(4)，确定轨道不平度的最佳拟合函数，包括如下过程：
[0025] (4a)根据最小二乘原理，构造方程组ΧΑ = Η，即：
[0027] (4b)根据轨道各点方位坐标的拟合系数矩阵X和轨道各点高度的拟合系数向量H， 可以计算矩阵A为：
[0028] A=(XTX)-VH
[0029] 其中，Χτ为拟合系数矩阵的转置；
[0030] 从而确定轨道不平度最佳拟合函数中各参数的数值，即计算出这些参数（ao， ai, ··· ,ak) 〇
[0031] 所述步骤(5)中，结合天线方位座架与滚轮的分布，按照四个滚轮对应的轨道各点 方位角X分别为
，根据大型天线服役过程中所处的方 位角Az和轨道不平度的最佳拟合函数f (X) =ao+aix+a2X2+. . .+akxk，计算天线滚轮对应的轨 道高度值：
[0033]所述步骤（6)中，根据天线滚轮对应的轨道高度值
偏差： ，计算天线服役过程中所处的方位角Az、俯仰角E1时的波束指向
[0035]其中波束指向单位为弧度，r为轨道半径，s为俯仰轴在天线座架坐标系中的高度。 [0036]本发明具有以下特点：
[0037] (1)本发明基于测量得到的轨道表面各点高度值，能快速给出轨道不平度的最佳 拟合函数，很好地描述了轨道不平度的实际轮廓。
[0038] (2)本发明建立了基于轨道不平度的天线指向模型，能准确计算出大型天线在任 意位姿下因轨道不平度而引起的波束指向偏差，可用于定量评价轨道不平度对大型天线波 束指向的影响，从而指导服役时天线性能的调整与补偿，具有较高的工程实用价值。
【附图说明】
[0039] 图1为本发明的流程图；
[0040] 图2为大型轮轨式天线整体示意图；
[0041] 图3为天线座架坐标系以及天线滚轮、方位架示意图；
[0042] 图4为轨道不平度对天线方位架的影响示意图；
[0043] 图5为轨道不平度测量值与最佳拟合曲线图。
【具体实施方式】
[0044] 下面结合附图及实施例对本发明作进一步说明。
[0045] 如图1所示，一种基于轨道不平度的大型天线任意位姿时波束指向偏差计算方法， 具体步骤如下：
[0046] 步骤1，确定天线轨道采样点间距，测量天线服役时轨道各点的高度值
[0047]根据大型天线的轨道尺寸，确定轨道采样点的间距，并利用徕卡数字水准仪或其 他高精度测量仪器，测量大型天线在服役一段时间后的轨道各点的高度值。
[0048] 步骤2,确定天线座架坐标系，对轨道各点高度值进行坐标变换
[0049] 以大型天线的轨道中心为原点，Z轴与方位轴重合，垂直于地面指向天空，Y轴指向 正南方向，确定天线座架坐标系，并基于步骤1中测量得到的天线服役时轨道各点的高度 值，对其进行坐标变换。
[0050] 步骤3,计算轨道各点方位坐标的拟合系数矩阵和轨道各点高度的拟合系数向量
[0051] 3.1假设轨道不平度的最佳拟合函数为：
[0052] f (x) =ao+aix+a2X2+. . .+akxk
[0053] 其中x表示轨道位置，f(x)表示相应的轨道高度值，k为拟合函数的最高次数，ao， ai，a2. . .ak为拟合函数中的待定系数；
[0054] 3.2基于步骤2中得到的经过坐标变换后的天线服役时轨道各点的高度值，分别是 (xi，hi)-_(Xi，hi)-_(XN，hN)，计算轨道各点方位坐标的拟合系数矩阵X:
[0056]其中，N为采样点个数；
[0057] 3.3计算轨道各点高度的拟合系数向量H:
[0059] 步骤4，确定轨道不平度的最佳拟合函数
[0060] 4.1根据最小二乘原理，构造方程组XA=H，即：
[0062] 4.2基于步骤3中得到的轨道各点方位坐标的拟合系数矩阵X和轨道各点高度的拟 合系数向量H，可以计算矩阵A为：
[0063] A=(XTX)-VH
[0064]其中，Χτ为拟合系数矩阵的转置；
[0065] 从而确定轨道不平度最佳拟合函数中各参数的数值，即计算出这些参数（ao, ai, ··· ,ak) 〇
[0066] 步骤5,计算天线处于任意方位角时所有滚轮对应的轨道高度值
[0067] 基于步骤3、4中确定的大型天线服役过程中所处的方位角Az和轨道不平度的最佳拟 合函数f (X) =ao+aix+a2X2+.. . +akxk，针对多数大型轮轨式天线的滚轮为四个，结合天线方位座
架与滚轮的分布，这里取四个滚轮对应的轨道各点方位角X分别为 ，计算天线滚轮对应的轨道高度值：
[0069 ]步骤6，计算天线在任意方位角、俯仰角时的波束指向偏差 [0070]基于步骤5中得到的天线滚轮对应的轨道高度值
，计算天线服役过程中所处的方位角Az、俯仰角E1时的波束指向 偏差：
[0072] 其中波束指向单位为弧度，r为轨道半径，s为俯仰轴在天线座架坐标系中的高度。
[0073] 步骤7,输出波束指向偏差至天线控制系统，调整补偿天线
[0074]基于步骤6中得到的天线服役过程中所处的方位角Az、俯仰角E1时的波束指向偏 差，将该波束指向偏差输出至天线控制系统(ACU)，对天线进行调整补偿。
[0075] 本发明的优点可通过一下仿真进一步说明：
[0076] 1.确定天线轨道采样点间距，测量天线服役时轨道各点的高度值
[0077] 本实施例中，以某50米口径的大型轮轨式天线为案例进行分析，见图2所示，其方 位轨道B直径为32.5m，四个天线滚轮A沿方位轨道设置，轨道不平度优于0.8mm。轨道方位架 C整体由21段拼接而成，轨道接触角度呈45°，接头处设置大型垫板以保证接缝位置两端轨 道等高，每段轨道中部设置小垫板，板厚均为30_。轨道方位架C上设有俯仰齿轮D，反射体E 架设在轨道方位架C上;轨道采用地脚螺栓和压板固定，呈倒T型结构。轨道上部采用42CrMo 锻造加工而成，经过表面调质处理后硬度能达到HRC22~HRC26,具有优异的接触强度及抗 弯强度。
[0078] 本实施例中，轨道不平度的测量选用了徕卡NI007型高精度水准仪，采样点越密 集，越能反映轨道形貌的真实性。根据测量要求，天线轨道采样点间距设为3°，全轨共设120 个测点，即采样点个数N=120,待起始位置的读数稳定之后，记录数值，然后将天线转动3° 到达下一个测点。
[0079] 2.确定天线座架坐标系，对轨道各点高度值进行坐标变换
[0080] 以大型天线的轨道中心为原点，Z轴与方位轴重合，垂直于地面指向天空，Y轴指向 正南方向，确定天线座架坐标系，如图3所示，为天线座架坐标系以及天线滚轮、方位架示意 图。当轨道存在误差时，会引起方位架的倾斜和扭曲，如图4所示为轨道不平度对天线方位 架的影响示意图。并根据上一步骤中测量得到的天线服役时轨道各点的高度值，对其进行 坐标变换。轨道不平度的测量结果和最佳拟合曲线见图5。
[0081] 3.计算轨道各点方位坐标的拟合系数矩阵和轨道各点高度的拟合系数向量
[0082] 本实施例中，多项式最高次幂k=20,基于得到的经过坐标变换后的天线服役时轨道各点 的高度值，可计算得到轨道各点方位坐标的拟合系数矩阵
，和轨 道各点高度的拟合系数向量
[0083] 4.确定轨道不平度的最佳拟合函数
[0084] 利用上一步骤中得到的轨道各点方位坐标的拟合系数矩阵X和轨道各点高度的拟 合系数向量H，通过软件编程计算，可得到轨道不平度最佳拟合函数中各参数的数值，即矩 阵A=[0…0 -0.0001 -0.0004 0.0098 -0.0678 4.5806]，轨道不平度的最佳拟合函数f (x) =-0.0001x16-0.0004x17+0.0098x18-0.0678x 19+4.5806x20〇
[0085] 5.计算天线处于任意方位角时所有滚轮对应的轨道高度值
[0086] 利用前一步骤中确定轨道不平度的最佳拟合函数，取天线所处的方位角Az = 95°，
，可计算得到天线滚轮对应的轨道高度值：
[0088] 6.计算天线处于任意方位角、俯仰角时的波束指向偏差
[0089] 本实施例中，轨道半径r=16.25m，俯仰轴在天线座架坐标系中的高度s = 27.3m。 利用上一步骤中得到的天线滚轮对应的轨道高度值，计算天线服役过程中所处的方位角Az = 95°、俯仰角El = 45°时的波束指向偏差为
（角秒）。下表给出基于轨 道不平度的大型天线不同位姿时波束指向偏差：
[0090] 表1不同位姿时的天线波束指向偏差
[0092] 7.输出波束指向偏差至天线控制系统，调整补偿天线
[0093] 将天线波束指向偏差输出至天线控制系统(ACU)，通过主控计算机判断并对天线 进行调整补偿，使天线达到服役性能要求。
[0094] 通过上述仿真可以看出，采用本发明的方法可以对轨道不平度进行快速、准确描 述，给出轨道不平度的最佳拟合函数，并基于轨道不平度，直接计算出大型天线在任意位姿 时的波束指向偏差，可用于天线服役过程中波束指向的调整与补偿，保证了天线服役过程 中的高指向精度要求。
【主权项】
1. 一种基于轨道不平度的大型天线任意位姿时波束指向偏差计算方法，其特征在于， 包括如下过程： (1) 根据大型天线的轨道尺寸，确定轨道采样点的间距，并利用高精度测量仪器测量大 型天线在服役一段时间后的轨道各点的高度值； (2) 确定天线座架坐标系，并对测量得到的天线轨道各点高度值进行坐标变换； (3) 根据天线轨道各点坐标和高度值，计算轨道各点方位坐标的拟合系数矩阵和轨道 各点高度的拟合系数向量； (4) 基于最小二乘原理，确定轨道不平度的最佳拟合函数； (5) 根据大型天线的方位角和轨道不平度的最佳拟合函数，计算天线处于任意方位角 时所有滚轮对应的轨道高度值； (6) 根据天线滚轮对应的轨道高度值，计算大型天线在任意方位角、俯仰角时的波束指 向偏差； (7) 输出天线波束指向偏差至天线控制系统，对天线进行调整补偿。2. 根据权利要求1所述的一种基于轨道不平度的大型天线任意位姿时波束指向偏差计 算方法，其特征在于，步骤(2)中，所述天线座架坐标系的确定，是W大型天线的轨道中屯、为 原点，Z轴与方位轴重合，垂直于地面指向天空，Y轴指向正南方向。3. 根据权利要求1所述的一种基于轨道不平度的大型天线任意位姿时波束指向偏差计 算方法，其特征在于，步骤(3)按如下过程进行： (3a)假设轨道不平度的最佳拟合函数为： f ( X )=曰0+曰 1X+曰2χ2+·.. +曰kxk 其中，X表示轨道位置，f(x)表示相应的轨道高度值，k为拟合函数的最高次数，ao，ai, 曰2…ak为拟合函数中的待定系数； (3b)根据测量得到的天线服役时轨道各点的高度值(XI山）…佔山）一(^加），计算 轨道各点方位坐标的拟合系数矩阵X:其中，N为采样点个数； (3c)计算轨道各点高度的拟合系数向量H:4. 根据权利要求1所述的一种基于轨道不平度的大型天线任意位姿时波束指向偏差计 算方法，其特征在于，步骤(4)按如下过程进行： (4a)根据最小二乘原理，构造方程组XA = H，即：(4b)根据轨道各点方位坐标的拟合系数矩阵X和轨道各点高度的拟合系数向量H，可W 计算矩阵A为： Λ=(χΤχ)-νΗ 其中，χτ为拟合系数矩阵的转置； 从而确定轨道不平度最佳拟合函数中各参数的数值，即计算出运些参数(a〇，ai，···， 祉）。5. 根据权利要求1所述的一种基于轨道不平度的大型天线任意位姿时波束指向偏差计 算方法，其特征在于，步骤(5)中，结合天线方位座架与滚轮的分布，按照四个滚轮对应的轨 道各点方位角X分别为，根据大型天线服役过程中所 处的方位角Az和轨道不平度的最佳拟合函数f(x)=a日+aix+a2x2+. . .+akxk，计算天线滚轮对 应的轨道高度值：6. 根据权利要求1所述的一种基于轨道不平度的大型天线任意位姿时波束 指向偏差计算方法，其特征在于，步骤（6)中，根据天线滚轮对应的轨道高度值，计算大型天线服役过程中所处的方位角 Az、俯仰角E1时的波束指向偏差：其中波束指向单位为弧度，r为轨道半径，s为俯仰轴在天线座架坐标系中的高度。
【文档编号】G06F19/00GK106092047SQ201610371031
【公开日】2016年11月9日
【申请日】2016年5月30日
【发明人】王从思, 肖岚, 项斌斌, 李素兰, 许谦, 保宏, 段宝岩, 吴江, 冯树飞, 王艳, 蒋力, 王娜
【申请人】西安电子科技大学