采用MVFSA算法的射电望远镜曲面优化研究 沈世云.pdf

该【采用MVFSA算法的射电望远镜曲面优化研究 沈世云 】是由【抱琴】上传分享，文档一共【9】页，该文档可以免费在线阅读，需要了解更多关于【采用MVFSA算法的射电望远镜曲面优化研究 沈世云 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。重庆理工大学学报(自然科学)
JournalofChongqingUniversityofTechnology(NaturalScience)
ISSN1674-8425,CN50-1205/T
《重庆理工大学学报(自然科学)》网络首发论文
题目:采用MVFSA算法的射电望远镜曲面优化研究
作者:沈世云,王扬帆,罗吉润,姜诲名
收稿日期:2022-04-01
网络首发日期:2022-08-01
引用格式:沈世云,王扬帆,罗吉润,
化研究[J/OL].重庆理工大学学报(自然科学).
.
网络首发:在编辑部工作流程中,稿件从录用到出版要经历录用定稿、排版定稿、整期汇编定稿等阶
段。录用定稿指内容已经确定,且通过同行评议、主编终审同意刊用的稿件。排版定稿指录用定稿按照期
刊特定版式(包括网络呈现版式)排版后的稿件,可暂不确定出版年、卷、期和页码。整期汇编定稿指出
版年、卷、期、页码均已确定的印刷或数字出版的整期汇编稿件。录用定稿网络首发稿件内容必须符合《出
版管理条例》和《期刊出版管理规定》的有关规定;学术研究成果具有创新性、科学性和先进性,符合编
辑部对刊文的录用要求,不存在学术不端行为及其他侵权行为;稿件内容应基本符合国家有关书刊编辑、
出版的技术标准,正确使用和统一规范语言文字、符号、数字、外文字母、法定计量单位及地图标注等。
为确保录用定稿网络首发的严肃性,录用定稿一经发布,不得修改论文题目、作者、机构名称和学术内容,
只可基于编辑规范进行少量文字的修改。
出版确认:纸质期刊编辑部通过与《中国学术期刊(光盘版)》电子杂志社有限公司签约,在《中国
学术期刊(网络版)》出版传播平台上创办与纸质期刊内容一致的网络版,以单篇或整期出版形式,在印刷
出版之前刊发论文的录用定稿、排版定稿、整期汇编定稿。因为《中国学术期刊(网络版)》是国家新闻出
版广电总局批准的网络连续型出版物(ISSN2096-4188,CN11-6037/Z),所以签约期刊的网络版上网络首
发论文视为正式出版。
网络首发时间:2022-08-0109:55:39
网络首发地址:.
收稿日期:2022-04-01
基金项目:国家自然科学基金项目(61906025);中国民航飞行技术与飞行安全重点实验室开放基金(FZ2021KF08)
作者简介:沈世云,男,副教授,主要从事数学建模研究,E-mail:******@;通讯作者王扬帆,男,主
要从事通信研究,E-mail:******@。
采用MVFSA算法的射电望远镜曲面优化研究
沈世云1,王扬帆2,罗吉润2,姜诲名2
(,重庆400065;,重庆400065)
摘要:依据FAST射电望远镜主动反射面相关结构和工作原理,针对射电望远镜面对不同方向天体观察时主动反射面需
要调整的问题,采用改进的非常快速模拟退火算法对射电望远镜曲面进行优化。以主动反射面单元面板的结点及其下
拉索促动器为研究对象,采用面形精度来反映反射面瞬时拟合程度。利用机理建模方式确定工作面范围和理想抛物面
方程;在满足约束条件的情况下,建立以工作面与理想抛物面拟合程度最高为目标的优化模型。为提高算法效率和拟
合质量,求解采用改进的非常快速模拟退火算法,并与其他算法求解的结果进行比较分析。Matlab仿真与求解结果表
明:改进的非常快速模拟退火算法较普通算法收敛速度更快,拟合结果更优。
关键词:改进的非常快速模拟退火算法;面形精度;FAST径射电望远镜;瞬时拟合程度
中图分类号:;:A文章编号:1671-0924(2007)xx-xxxx-x
Researchonsurfaceoptimizationofradiotelescopebasedonmodified
veryfastsimulatedannealing
SHENShiyun1,WANGYangfan2,LUOJirun2,JIANGHuiming2
(,ChongqingUniversityofPostsandTelecommunications,Chongqing400065,China;
,ChongqingUniversityofPostsandTelecommunications,
Chongqing400065,China)
Abstract:AccordingtotherelatedstructureandworkingprincipleoftheactivereflectorofFASTradiotelescope,aimingat
theproblemthattheactivereflectorneedstobeadjustedwhentheradiotelescopelooksatcelestialbodiesindifferentdirections,

activereflectorelementpanelanditslowercableactuatorastheresearchobject,thesurfaceshapeaccuracyisusedtoreflectthe

mechanismmodeling;Undertheconditionofmeetingtheconstraints,theoptimizationmodelaimingatthehighestfitting

thealgorithm,themodifiedveryfastsimulatedannealingalgorithmisused,andtheresultsarecomparedwithotheralgorithms.
ThesimulationandsolutionresultsofMATLABshowthatthemodifiedveryfastsimulatedannealingalgorithmhasfaster
convergencespeedandbetterfittingresultsthantheordinaryalgorithm.
Keywords:modifiedveryfastsimulatedannealing;surfaceshapeaccuracy;FASTradiotelescope;instantaneousfitting
degree
大、灵敏度最高的射电望远镜,其落成启用对中
0引言国加快建设科技强国、实现高水平科技自立自强有
重大意义。中国天眼FAST射电望远镜的控制系统
中国天眼FAST射电望远镜是国家“十一·五”立主要包括望远镜总控制、馈源支撑控制和反射面
项的国家大科学工程。作为目前世界上单口径最系统三大系统。主反射面设计的目的是调节反射
面的形状和方向,
波可以通过反射面集中在馈送舱的有效区域。将主主动反射面包括基准态与工作态。如图2的
动反射面调整为工作抛物面是主反射面技术的重要FAST剖面示意图所示,基准球面与焦面是以点C
过程,该过程通过将主索结点上的促动装置与下拉为圆点的部分同心球面,其中两面半径差为
索结合来实现。F=(R=),基准面可视作半径R=300
目前,已有众多学者对大口径射电望远镜的主m、口径D=500m的球面。当FAST进入观测状态
动反射面调节优化或智能算法优化进行了研究。朱并需要接收天体S发出的平行电磁波信号时,主动
丽春[1]通过主反射面控制系统控制策略,建立自适反射面处于工作态。为了反射面可以将来自目标
应、自学****机制节点位移控制模型,采用开闭环结观测天体的平行电磁波汇聚到馈源舱上直径为1m
合控制策略实现反射面高精度实时成型;薛建兴等的有效圆盘区域,馈源舱接受平面的中心将在焦
[2]基于反射面单元动态面形精度分析研究瞬
FAST面上移动,具体位置为直线lSC与焦面的交点P。同
[3]
时抛物面的拟合精度;王从思等根据促动器支撑时,主动反射面上有1个口径为300m的区域需要
节点位置确定该点促动器所在点的法线方向及促动调整为以点P为焦点的反射抛物面。
器所在直线方程,计算出面向抛物面的大型天线主
动面促动器最佳调整量;陈聪聪等[4]研究了一种确
定网状反射面天线原理误差计算的方法;张建春等
[5]基于传统模拟退火算法提出了更高求解效率的非
常快速模拟退火算法;王星童等[6]利用快速模拟退
火算法加快了路径规划问题的迭代收敛速度;陈华
根等[7]提出了一种改进的非常快速模拟退火算法,
相比传统算法有了更高的稳健性和求解效率。
由于传统方法全优化的时间冗长,本文基于现
有研究,依据FAST射电望远镜的主动反射面相关
结构和其主动反射面工作原理,针对反射面对不同
方向天体观察时工作面需要调整的问题,建立以工
作面与理想抛物面拟合程度最高为目标的优化模
型,提出了采用改进的非常快速模拟退火算法来优
化射电望远镜工作面,并对比不同算法的求解效
图2FAST剖面示意图
率。
2理想抛物面的确定
1天眼FAST的组成与反射面状态


通过文献[8]得知,FAST望远镜从基准态转变
主要是由反射面、望远镜总控、馈源支
FAST为工作态的过程是不断反馈、优化和调整的过
撑三大控制系统、信号接收系统以及支承系统组
程,通过馈源舱对平行电磁波的接收结果,控制
成。其中,主动反射面系统结构如图所示,是由
1系统对没有将光束射入馈源舱的反射面单元进行
背架、反射面板、主索网、促动器及下拉索等主要
调整和优化。整个系统完全优化完毕共耗时约
部件构成的一个可调节球面。每个主索节点连接一
11000s,调整优化的过程十分庞杂。故在求取理
根与固定在地表的促动器连接的下拉索,通过下拉
想抛物面时,选取FAST即将从基准态转变为工
索和促动器的调节实现对主索网的形态控制。
作态时的瞬间进行分析,尽可能忽略变化过程的
影响。

如图3所示,L1、L2、L3绕对称轴旋转360°
得到同样的抛物面,其上任一抛物线的准线共平
面称这个平面为抛物面的准面。对任意一条抛物
线而言,抛物线上的任意一点到焦点的距离等于
图1主动反射面系统结构示意图
该点到准线的距离。可以推得,抛物面上的任意一结合上文可知,由于点P在直线lSC上,所以
点到焦点的距离等于该点到准面的距离。考虑使用直线参数方程的形式来求解点P坐标。
根据图2可得CP=(1-)R=,所以
tsp=。将其代入式(1)可得点P的坐标为
(-,-,-)。

图5是基准面和理想抛物面的剖面图。从剖面
看,其抛物面的任意一条母线上的点到焦点的距
离等于其到准线的距离。,抛物面上
任意一点到焦点的距离等于其到准面的距离,所
以可根据此关系建立方程。
图3抛物面与抛物线的关系


图5理想抛物面与基准面示意图(剖面)
以图2中点C作为坐标空间系的原点。如图4
所示,设点M为直线lSC上的任意一点,直线方向设抛物面上任意一点B坐标为(x,y,z),已知
α=°,β=°。设m为直线SC的方向向点P的坐标为(-,-,-)。根据空间中
量,可以求得直线lSC的方向向量为:两点距离公式求出点B到点P距离dBP。
mcoscos,cossin,sin
求取lSC与准面的交点。设直线lSC与准面交于
lSC直线参数方程为:点T。由抛物面性质可知,直线lSC垂直于准面,
xt。联立直线
l的方向向量m与点T坐标,得到准面方程为:
ytt,为参数SC

zt
求得直线lSC方程的参数形式可以减少变量数
点B到准面的距离k以及点B到焦点P的距离
量,在仅知道t的情况下就可以得出直线上点的坐
d表示为:
标,简化了求取结点坐标以及理想抛物面的过程。BP
222
dxyzBP(1)
x
k(2)
12
令式(1)与式(2)相等,得到一个抛物隐式
方程。
理想抛物面顶点坐标说明如表1所示。
表1理想抛物面顶点坐标说明
x坐标y坐标z坐标
---
图4直线SC在空间中的示意图

3优化模型的建立
确定一个理想抛物面后,可以通过调节促动器
的径向伸缩量将反射面调节为工作抛物面,
工作抛物面尽量拟合理想抛物面,以此使得平行电如图8所示,由于促动器只能沿径向运动,所
磁波尽可能多地反射汇聚到馈源舱上的有效圆盘区
以主索结点的运动轨迹是沿直线lAC运动的。P1为
域。由于工作区口径只有,所以需要先确定
300m主索结点的最理想位置,位于理想抛物面上;P2为
索网中促动器需要伸缩改变的主索结点。
通过促动器调节后的主索结点的新位置;点P3为基
准球面上的主索结点,即主索结点的初始状态;点
P4为促动器;Δx1为主索节点与促动器相连的距离;
Δx2为初始位置与新位置的径向差值;Δx3为新位置
与最理想位置的径向差值。
图6工作态口径范围剖面图示意图
图8主索结点与理想抛物面的位置关系图(剖面)
如图6所示,以点C为原点,平行于地平线的
当Δx3越小时,证明经过促动器调节后的主索
直线作为x轴建立平面直角坐标系,口径范围的剖
结点新位置更靠近理想抛物面上的位置,此时调节
面图即为扇形CWK区域。同样,选取FAST即将的效果最好。为此,引入面形精度RMS来衡量反
从基准态转变为工作态时的瞬间进行分析,得到射面形精度。反射面面形精度是衡量射电望远镜综
CW=300m,WK=300m,WV=WK/2=150m,所以合性能的重要指标之一,直接影响望远镜的工作效
。对索网中全部个主索结点进行
2226率。本文中用反射面主索结点实际坐标与理想抛物
遍历,分别计算第i个主索结点与C点的直线li的面拟合误差来衡量反射面的RMS:
方向向量与CG方向向量的夹角,利用如下公式:
695
abXxYyZz222
cosiiiiii
abRMSi1
695
计算得到各个夹角θ。若θ∈[0°,ɛ+30°],则主索节
其中:Xi、Yi、Zi分别表示第i个节点调节成理想抛
点处在口径范围内;若θ在此范围内,则主索节点
物面时的理想坐标,位于理想抛物面上;xi、yi、zi
不属于口径范围。为了保证工作面完整性,引入1表示为第i个节点调整后的实际坐标。计算695个
个任意小的数ɛ。待调整的主索结点如图7所示。实际坐标与理想坐标的距离平方和的平均值,然后
通过计算得到,共有695个主索结点需要。再开平方得到RMS。

:伸缩长度限制
各结点促动器伸缩调整长度为ti,其中
i=1,2,3,……,695。根据实际意义,认为趋向球心时
ti取正,远离球心时ti取负。长度的调整需要在一
定范围内,因此可以得到第1个关于促动器伸缩调
整长度的约束条件:
图7待调整的主索结点示意图
ti
:伸缩方向限制成功使用SA解决了组合优化问题。SA作为一种
有效的非线性组合寻优算法,在理论上已经被证
促动器仅可以在径向方向上变化,所以采用直
明是严格的,且有大量的实际应用可以证明该算
线方程的参数形式进行径向变化量表示:
法有效,可靠,能够处理规模较大的背包问题,且能
uiiiiu'tm[9-10]
够鲁棒地得到满意解。SA算法的关键思想是:
viiiiv'tn算法会有一定概率接受一个劣解并跳出某个局部
最优解空间并尝试寻找更优解。从全局搜索的角度
wiiiiw'tq
来看,算法会尽力朝最优解进行搜索寻优。

算法原理和金属固体退火的原理相似,就
设调整后节点之后,第i个结点原坐标调整为SA
是将热力学的理论运用到统计学中。从机理角度
ui、vi、wi,在调整坐标时,通过调整各促动器的伸
缩量来实现各结点向理想坐标位置逼近,从而实现解释算法思想就是先将金属固体加温到一个足够
工作面向理想抛物面逼近。综合考虑RMS精度值高的温度,再逐渐降温冷却;在降温的过程中粒
越小越好的目标,以及需要满足调整长度的范围和子会自发地进行热运动以使得固体状态朝着自由
调整方向的约束条件,建立如下优化模型:能减小,以此朝着能达到自由能最小的平衡态进
行单点寻优。
Xutm2
iiii具体步骤如下:
695
2步骤1给定一个较高且合适的初始温度T,
Yvtniiii0
i1在符合模型给定变量区间内随机选择一个初始模
2
Zwtqiiii型,并计算初始模型对应的目标函数值;
minRMSm0E(m0)
695步骤2采用合适的扰动方法,在符合模型给
定变量变化区间内,扰动当前模型产生一个新模

型m,并计算新模型的目标函数值E(m),同时计
uutmiiii'算目标函数差值ΔE=E(m)-E(m);
0
.'iiii步骤3在当前温度T下,若ΔE<0,则接受新
模型m;否则,按照概率P接受当前较劣模型;当
wwtqiiii'
模型被接受时,用新模型m替换初始模型m;
i1,2,3,,6950
步骤4在温度T下,重复进行多次扰动和接
695个主索结点对应关系见表2。其中,单位方受判别,即重复步骤2与步骤3共Lk次;
向向量为原坐标指向球心的向量单位化后所得向步骤5通过一定的降温退火计划下降温度T;
量。步骤6重复步骤2—5,直到温度下降到终止
表2695个主索结点对应关系温度或者搜索达到收敛条件,认为本次已经搜索
到一个较优模型解。
序号原坐标方向向量理想坐标

1uvw111,,mnq111,,XYZ111,,
Ingher在1989年提出了一种模拟退火算法的
2uvw,,m,,nqXYZ,,
222222222改进算法,非常快速模拟退火算法(VeryFast
3
uvw333,,mnq333,,X3,,Y3Z3SimulatedAnnealing,VFSA)。VFSA算法较SA算
…………………法具有更高的效率,使得模拟退火理论真正走上
u,,vwm,,nqXYZ,,
695695695695695695695695965695了实际应用并可以解决更多的实际问题。VFSA算
法的流程与SA算法流程一致,只在模型扰动、接
通过求解上述优化模型可以得到最为拟合理想受概率或者退火计划上做出了改进。相比传统算
抛物面的工作面和各个主索结点的变化量ti。法,该算法在高温时有着更高的全局搜索效率。
但由于在低温下仍然进行全局搜索,搜索的扰动
4优化模型的求解幅度过大,难以继续搜寻更优模型,因此在低温
状态下搜索方式可以作出改进以增加低温搜索效

率。本文基于改进的非常快速模拟退火算法
最早是在年,提出模拟退火算
1953Metropolis(modifiedveryfastsimulatedannealing,
法(SimulatedAnnealing,SA);Kirkpatrick在1983年
[7]
MVFSA),利用Matlab对模型进行求解。
(8)
TTjj01exp/

式(7)和(8)由式(5)和式(6)改进而来,其中
模拟退火算法是在当前温度下扰动当前模型产L(j)为算法搜索空间限制因子,与算法迭代次数有
生一个新模型,算法扰动由一个随机分布函数控关。当迭代次数j增大时,L(j)也增大,这意味着搜
制。VFSA算法使用依赖温度的似Cauchy分布扰动索范围将变小,通过搜索空间限制因子来达到算
产生新模型。模型扰动规则如下:法在一个更小的搜索区间内寻找更优模型的目
的。j1为上一阶段迭代次数;ρ为降温形状因子,
21
rTT1(3)与当前温度成反比关系,意味着当越小时越
iTTρ
大。算法的改进目的是使得搜索初期即在高温时
有较大的搜索范围和劣解接受概率,在低温时可
mrmm'()(4)
iiiimaxmin以在已找到的最优解区间继续以较小的搜索范围
’
式(3)和式(4)中:mi表示扰动后的模型,mimax和寻找更优解。
mimin分别表示mi取值范围的上、下限;
子,其中γ表示(0,1)区间上的随机数;sgn(x)
函数。模型进行全局搜索和局部搜索。设置MVFSA算法
的温度衰减系数,初始温度,终止
==200
-7
温度为Tend=10。
模拟退火算法中新模型根据目标函数的差值通过随机函数的生成,选取1个在模型取值区
,按照准则来接受新模型,目的是跳
ΔEMetropolis间内的合法初始模型m0(t1,t2,……,t695),并根据面
出局部最优解。该准则是传统模拟退火的重要环
形精度计算该初始模型的函数值RMS0。每次达到
节。平衡后都会进行降温,,
若(新模型较优),则新模型无条件被
1)ΔE<0直到T<Tend则认为阶段1的全局搜索结束,已经
接受;获得全局最优解。此时,对模型进行适当回火使得
2)若ΔE>0(新模型较劣),计算接受概率P:模型处于低温,
PETexp/(5)索,当该阶段收敛或者达到终止温度时认为已找
在ΔE>0的情况下,产生一个分布在[0,1]内的到最优解。
随机数ξ,若P>ξ则接受新模型,否则舍弃新模型。在相同初始解下,进行50次的MVFSA算法
搜索,取平均值得到RMS=。695个主索

结点对应关系见表3。
诸多实践证明,采取温度以指数变化方式降
T表3695个主索结点对应关系
温的计划较符合固体退火的实质:
序号理想坐标点优化后坐标伸缩量

TTjexp1(6)
01(0,0,-)(0,0,-)
其中:ν为模拟退火算法中的温度衰减系数;j为算2(,,-)(,,-)-
法迭代次数。3(,-,-)(,-,-)
…………………
局部搜索阶段
(-,-,-)(-,-,-)

全局搜索以找到最优区间解。当完成全局搜索算法
达到收敛时,进行回火升温操作,此时从一较低温
度开始退火,继续在已找到的最优解区间内进行局5结果分析与算法比较
部搜索寻优。本阶段中改进的模型扰动与退火计划
表示为:
如图9所示,在MVFSA算法下,算法迭代13
m'mr(mm)/Lj(7)
iiiimaximin次找到1个全局最优解区间,在温度很低时再经过
4次迭代在该区间内找到更优解,提高了结果的精
度和准确性。但是解的优化效率并没有提高很多,
原因是模型变量多达695个,且每个变量取值区间
相对搜索范围较狭窄,很难大幅度提高寻优结果。
图10不同算法各代最佳函数值曲线图
上述结果表明:求解本文中建立的优化模
型,采用MVFSA算法有更高的搜索效率和更优的
搜索结果。
图9MVFSA算法各代最佳函数值曲线

为了验证算法的高效性,将MVFSA算法、
针对FAST射电望远镜的主动反射面调节优化
VFSA算法和传统的SA算法进行比较。结果如表4
问题进行研究。通过机理建模的方式从下拉索的调
与图10所示。比较在同一初解、相同初始温度T0
-7节量出发,基于改进的非常快速模拟退火的智能
=200、终止温度Tend=10和温度衰减系数ν=
算法对衡量拟合效果的面形精度优化模型进行求
时3种算法的搜索效率。在迭代开始时,VFSA算
解。以695个待调整主索结点的下拉索促动器伸缩
法和MVFSA算法迭代曲线斜率的绝对值远大于传
量为研究变量。该算法在非常快速模拟退火算法全
统SA算法,且在迭代约8次时就快速找到了1个
局快速搜索的基础上,继续在局部进行快速搜
较优解区间。而在相同迭代次数下,传统SA算法
索,在几乎相同且短暂的总搜索时间下得到更优
仍在做低效或者无效的搜索。
的搜索结果,使得主动反射面的工作面更加拟合
MVFSA算法在适当回火至一个较低温度时进
理想抛物面。
行局部搜索,相比VFSA算法多进行了5次迭代搜
实际中,采用该算法或智能算法来解决射电
索,收敛也十分迅速。二者总搜索时间几乎一致,
望远镜的主动反射面调节案例较少。本文中提出的
但采用MVFSA算法得到的RMS值较原来下降了
优化模型可以推广到其余口径天线或者射电望远
%,意味着拟合理想抛物面程度更高。
镜的主动反射面调节优化中,提供了一种较好的
表4退火计划表与寻优结果优化思路,有利于馈源舱更好地接受天体发出的
迭代次数/次,时
算法初始温度温度衰减系数RMS电磁波信号。后续研究将考虑更多影响因素,进
间/s
一步提高算法的效率和准确性,提高主动反射面
,
的工作效率。
,
,
参考文献:
[1](FAST)主动反射面
整网变形控制[J].科研信息化技术与应用,2012,3(4):
67-75.
[2]薛建兴,
时抛物面拟合精度的预估与改善[J].光学精密工程,
201