2021年数学建模国赛A题优秀论文（Word）（FAST”工作抛物面的优化设计）

“FAST”工作抛物面的优化设计

摘要

本文基于FAST的工作原理，通过机理分析、坐标変换、非线性最小二乘优化等方 法，建立了反射面板谟节优化模型・并利用BFGS算法、蒙特卡洛积分算法等算法，对 不同条件下反射光线吸收比率进行了研究。

问题一中，首先基于固定的仰角观测目标S、圆心C和焦点P・利用旋转抛物 面的中心对祢性，选取焦距作为自由度控制变量，构建在极坐标系下开口竖直向上的二 维抛物线方程.得到不同偏转角度下原点到抛物线的距离.进而导出三维下的旋转掀物 面方程。其次，以焦距为决策变量，将口径300米的拋物面作为积分域•将理想抛物面 到原点的距离与基准球面半径差值平方作为被积函数进行积分作为最小化目标函数.建 立了确定理想抛物面的优化模型。最后，使用二分法求得目标函数导函数在定义区间上 的零点.得到理想抛物面焦距的精确值为280.854,误差平方积分的最小偵为10.112o 此时对应理想抛物面的解析式为z= Q + #)2/561.708 300.841,

问题二中，首先利用球坐标下不同轴线方向抛物面的旋转不变性.在原坐标系和问 题一的坐标系之间建立了双向可逆的变换关系.得到了不同方位角下理想抛物面到原点 的距离。其次.以主索节点的工作坐标和促动器的伸缩长度为决策变量:.以积分域覆盖 的主索节点到原点的距离与理想抛物面到原点的距厲之差的平方和为最小化目标函数. 分别考虑下拉索长度固定、相邻节点的距离变化幅度不超过0.07%.促动器的伸缩范围 在±0.6m为约束条件.建立反射面板调节优化模型。最后，使用拉格朗日乗子法和BFGS 算法进行求解.得到误差平方在抛物面口径上的积分的最小值为5.1353X10 9.理想抛 物线的顶点坐标为(-49.392,-36.943,-294.450).调节后反射面300米口径内的主 索节点编号、位置坐标、各促动器的伸缩量等结果见文件result.xlsxe

问题三中，首先通过旋转变换.将反肘问题的倾斜入射光线转化为垂直入射光线。 其次.通过求解线性方程组.依次确定出入射光线与三角形面板的相交判定式、交点坐 标、三角形面板的法线向量，并利用光线垂直入射的性质，使用法线向最简化计算得到 出射光线的方向角。再次，通过联立射线方程与馈源舱所在的目标髙度.得出射线方程 的步长和出射光到达目标高度时的坐标.并与馈源鲍的有效区域进行比对.作为入射光 线是否被有效接收的判定式。最后，将300米口径内实际接收区域作为积分域，将入射 光线有效接收判别式作为被枳函数.使用蒙特卡洛算法进行积分，得到有效接收的光源 面积，并计算出调节昉接收比0.811%.调节后接收比1.103%.提升了 36%,调节工作 抛物面的过程使有效光源的光斑可以尽量完整地出现在每个三角形面板内(图】1)。

本文的特色在于将机理分析与非线性最小二乘优化相结合，并灵活釆用二分法、 BFGS算法和蒙特卡洛积分算法进行求解，在大规模、髙维且带有非线性约束的求解中 仍然以接近二阶的速度收敛至最优解.为FAST在不同情况下的调节与设计提供了参考 依据。

关键词：坐标旋转、非线性最小二乘、BFGS算法、蒙特卡洛积分

一、  问题重述

1.1问题背景

中国天眼(FAST)由主动反射面、馈源舱及其他系统组成,其中主动反射面是由主 索网、反射面板、下拉索、促动器及支承结构等构成的可调节球面。主索冋由柔性主索 按短程线三角冋格方式构成，每个主索节点连接一根下端与固定在地表的促动器连接的 下拉索。促动器沿基准球面径向安装，可沿径向伸缩。

主动反射面有两个状态：基准态时反射面为半径约300米、口径为500米的球面： 工作态时反射面为一个300米口径的近似旋转抛物面。馈源館接收中心在与基准球面同 心、半径差为的一个球面上移动。当观测某个方向的目标S时.馈源舱接收中心被移 动到直线SC与焦面的交点P处，调节基准球面上的反射板形成以直线SC为对称轴、 以P为焦点的近似旋转抛物面.从而将来自S的平行电磁波汇聚到有效区域。

1.2问题提出

在反射面板调节约束下，确定一个理想抛物面，调节促动器，使工作拋物面尽量贴 近理想掀物面，以获得反射后的最佳接收效果。建立模型解决以下问题：

•问题一：当待观测天体S位于基准球面正上方，即。=0°,。= 90°时，结合考 虑反射面板调节因素.确定理想拋物面。

•问题二：当待观测天体S位于。=36.765°, 0 = 78.169°时.确定理想抛物面。 建立反射面调节模型，调节相关促动器的伸缩量.使反射面尽量贴近该理想抛 物面O

•问题三：基于问题二方案，计算调节后馈源能接收比，与基准球面接收比作比 较。

二、  问题假设

1. 假设光线在反射面不存在二次反射。
2. 假设不考虑板间间隙给反射带来的影响。
3. 假设电磁波在大气中的传播为直线传播,
4. 假设反射面板不会发生弯曲，始终维持一个平面。
5. 假设电磁波的反射为全反射,不考虑反射时的损耗。
6. 假设主索节点调节后，相邻节点之间的距离会发生微小变化.变化幅度不超过 0.07%。

三、符号说明