面向天文多普勒差分测速的太阳_行星光谱对生成方法.pdf

第 45 卷第 2 期 2024 年 2 月宇 航 学 报Journal of AstronauticsNo.22024FebruaryVol.45面向天文多普勒差分测速的太阳/行星光谱对生成方法刘劲1，2，徐玉豪1，尤伟3，陈晓3，张子军1，马辛2（1.武汉科技大学信息科学与工程学院，武汉 430081；2.北京航空航天大学仪器科学与光电工程学院，北京 100191；3.上海卫星工程研究所，上海 201109）摘要：为了提供天文多普勒差分测速所需的同步太阳/行星光谱对，提出了一种变分自编码器（VAE）和对偶生成对抗网络（Dual GAN）相融合的VAE-Dual GAN。首先，实测太阳光谱经过VAE编码到隐空间，实现了光谱到光谱域的扩充；然后，由Dual GAN将隐空间映射到伪行星光谱；最后，利用伪行星光谱生成重构太阳光谱。此外，利用编码和生成重建损失加强对网络的约束。VAE-Dual GAN利用Dual GAN的转换学习能力完成了两个光谱域的转换，生成同步太阳/行星光谱对。实验结果表明，VAE-Dual GAN可生成高质量的太阳/行星光谱对，将天文多普勒差分测速精度提高60%以上。关键词：天文导航；测速导航；太阳/行星光谱对；生成对抗网络；变分自编码器中图分类号：P128.4；VS29.1 文献标识码：A 文章编号：1000-1328（2024）02-0273-10 DOI：10.3873/j.issn.1000-1328.2024.02.012Generative Method of Solar/Planetary Spectral Pair for Celestial Doppler Differential Velocity MeasurementLIU Jin1，2，XU Yuhao1，YOU Wei3，CHEN Xiao3，ZHANG Zijun1，MA Xin2（1.College of Information Science and Engineering，Wuhan University of Science and Technology，Wuhan 430081，China；2.School of Instrumentation and Optoelectronic Engineering，Beihang University，Beijing 100191，China；3.Shanghai Institute of Satellite Engineering，Shanghai 201109，China）Abstract:In order to provide synchronous solar/planetary spectral pairs required for the celestial Doppler differential velocimetry，a VAE-Dual GAN，which combines variational auto encoder（VAE）with dual generative adversarial networks（Dual GAN）is proposed.Firstly，the measured solar spectrum is encoded into the latent space by VAE，which realizes the expansion from spectrum to spectral domain.Then，the latent space is mapped to the pseudo-planetary spectrum by Dual GAN.Finally，the reconstructed solar spectrum is generated by using the pseudo-planetary spectrum.Furthermore，the encoder and generator reconstruction losses are used to strengthen the constraints on the network.VAE-Dual GAN utilizes the conversion learning ability of Dual GAN to complete the transformation of two spectral domains and generate synchronous solar/planetary spectral pairs.The experimental results show that the VAE-Dual GAN generates high-quality solar/planetary spectral pairs，the accuracy of celestial Doppler differential velocity measurement is improved by more than 60%.Key words:Celestial navigation；Velocimetry navigation；Solar/planetary spectral pairs；Generative adversarial networks；Variational auto encoder收稿日期：2023-05-11；修回日期：2023-08-27基金项目：国家自然科学基金（61873196，62373030，61772187）；科技创新2030-“量子通信与量子计算机”重大项目（2021ZD0303400）宇航学报第 45 卷0引言进入21世纪，航天大国纷纷开启各自的深空探测计划。中国航天筑梦苍穹：“嫦娥”奔月1、“祝融”探火2、“羲和”追日3，这些航天任务的成功实施标志着中国的深空探测技术日趋成熟。导航作为未来深空探测计划必不可少的关键技术，需满足的要求日趋增高。尤其是在行星探测的捕获阶段，航天器与地面站之间的通信延迟长，而捕获期短，地面站无法为航天器提供实时高精度的导航信息。天文导航不依赖地面站，成为解决这一问题的关键技术。天文导航有测角、测距、测速3种方式。测角导航和测距导航分别以航天器与近天体夹角4、脉冲星到达时间5作为测量量，通过卡尔曼滤波，间接获得航天器速度。为了直接获得速度信息，测速导航应运而生，它根据多普勒效应，测量天文光谱频移，获得速度信息，弥补了测角与测距导航的不足。通常情况下，测速会与测角或测距相结合6-8，以提高组合导航的精度。上海卫星工程研究所和中国科学院研制的太阳原子鉴频器9已完成在轨验证，证实了天文组合导航的可行性与精确性10。天文多普勒差分导航11通过差分技术消除了太阳活动的影响，成为研究热点。其测量量是太阳光谱频移和行星光谱频移之差。目前，针对天文光谱频移估计方法的研究进行得如火如荼，且已取得了一系列研究成果。欧洲南方天文台（ESO）的高精度径向速度行星搜索器（High accuracy radial velocity planet searcher，HARPS）利用互相关函数算法估计径向速度12。模板增强径向速度再分析应用算法比互相关函数算法具有更高的精度13。基于最小二乘匹配的光谱径向速度分析仪的径向测速精度14可达 1 ms-1，用其对HARPS 径向速度数据库的系统误差进行校正15。Silva等16提出了一种基于模板匹配的半贝叶斯方法。Rajpaul等17-18提出了无需高信噪比模板，仅用高斯过程对光谱建模测速的方法。文献 19 通过建立偏移和补偿模型降低太阳自转引起的多普勒偏差给测速造成的影响。文献 20 提出针对边界效应的天文光谱测速方法，解决天文光谱测速中的非线性傅立叶相移问题，使天文测速精度进一步提升。以上方法均是观测同一个天体，其光谱并不会受行星反射，自转等影响。天文差分测速则是利用太阳和行星光谱对比，行星光谱必然受到多种因素干扰。因此，不能照搬传统天文测速方法，而是要根据天文差分测速的特点，研究新的信号处理方法。太阳是一个宏观稳定，但微观极不稳定的星体。由太阳耀斑等引起的太阳光谱扰动21-22、光谱线深度的差异以及行星反射率都会对多普勒差分测速造成影响23。通过建立速度误差模型24，对太阳多普勒差分导航中直接和反射太阳光源的面积重叠率、时间色散以及太阳自转多普勒误差3个几何误差源进行分析25。多普勒差分测速原理上要求测速光谱产生于同一时刻离开太阳的太阳光子26。但目前缺乏时间同步的太阳直射和行星反射实测光谱数据。在图像生成方面，学者们先后提出变分自编码器（VAE）和生成对抗网络（GAN）。GAN自2014年被提出至今，在图像处理、机器视觉等领域被广泛引用，并改进出多个变体，如 Wasserstein 生成对抗网络30-31、条件对抗生成网络32以及双判别器对抗生成网络33。按照配对原则可将GAN分为非配对网络和配对网络，如循环生成对抗网络34（Cycle GAN）和对偶生成对抗网络35-36（Dual GAN）。由于行星和太阳光谱的轮廓大致相同，仅存在由宇宙中的各种影响因子引起的细节上的不同，符合 Dual GAN的主体不变、细节转换的特点，所以此次实验选用Dual GAN。同时Dual GAN继承了Wasserstein生成对抗网络损失函数不会梯度消失的优点，这也是选择Dual GAN的主要原因。为了解决同步光谱对缺失的问题，本文利用Dual GAN 将非同步的太阳/行星光谱对转换为同步的太阳/行星光谱对，用于多普勒差分测速。针对训练集有限，进而导致泛化性差这一问题，本文把Dual GAN与VAE结合起来，利用VAE的数据增强能力增广训练集。解码器与生成器的作用相似，同为生成过程。因此，用Dual GAN的生成器取代 VAE 中的解码器，使 VAE 和 Dual GAN 融合为VAE-Dual GAN。同时，利用Dual GAN网络的重建损失约束作用同时训练生成器、判别器和 VAE 的编码器。1VAE-Dual GAN为了生成同步的太阳/行星光谱对，考虑到VAE的数据增广能力，本文将VAE和Dual GAN结合起274第 2 期刘劲等：面向天文多普勒差分测速的太阳/行星光谱对生成方法来，提出了 VAE-Dual GAN，介绍了其网络结构，并从理论上分析了子网络的泛化性。1.1网络模型结构VAE-Dual GAN 结构如图 1 所示，模型由编码器、生成器和判别器各2个组成。生成器重构损失和编码器重构损失分别用虚线和点划线表示。黑色虚线框内为VAE-生成器模块。黄色和蓝色的实线箭头分别表示伪太阳光谱和伪行星光谱生成过程。在伪行星光谱生成过程中，首先，编码器sun将真实太阳光谱Ssun映射到太阳光谱隐空间Zsun，并在Zsun内采样得到太阳光谱样本zsun。然后，生成器sunplanet将zsun转化为与之同步的伪行星光谱Splanet。最后，伪行星光谱Splanet经过对偶编码器planet和对偶生成器planetsun生成重构太阳光谱Ssun。伪太阳光谱生成过程与之类似。Splanet表示真实行星光谱；Ssun和Splanet分别表示伪太阳光谱和重构行星光谱；Zplanet和Zsun分别代表真实行星光谱和伪太阳光谱经过编码器后的隐空间。图1中黑色虚线框中的VAE-生成器模块结构如图2所示。该模块利用生成器与VAE网络中的解码器作用相似的原理，用 Dual GAN 的生成器代替VAE解码器，构成VAE-生成器模块。Ssun域中的太阳光谱样本为ssun，隐藏层Hj上方的数字代表神经元数，隐藏层之间的FC代表全连接，编码器通过采样输出均值()i和方差2()i，编码后的隐空间样本zsun服从正态分布N()i，2()i)，然后通过生成器模块生成Splanet域的伪行星光谱。伪光谱真实光谱生成器重构损失生成器重构损失真实光谱伪光谱ZplanetZplanetSsunSsunZsunZsunSplanetSplanetSplanet编码器重构损失编码器重构损失VAE-生成器VAE-生成器SsunsunplanetsunplanetplanetDplanetDsunsun图1VAE-Dual GAN结构图Fig.1VAE-Dual GAN structureVAE生成器FCFCH1H2H3H4SsunzsunN(i),2(i)Splanet200100200100808080(i)2(i)图2VAE-生成器结构图Fig.2VAE-generator structure275宇航学报第 45 卷VAE 损失函数由 KL 损失LKL和编码器重构损失Ler构成，其表达式为：LVAE=-LKL+Ler=-2(1+ln21-21-21)+ssun Ssun sunplanet()sun()ssun-splanet2 -2(1+ln22-22-22)+splanet Splanet planetsun()planet()splanet-ssun2（1）式中：1和2分别为太阳和行星经过编码器输出的均值；21和22分别为二者的方差；(0，1)为控制KL散度的约束能力的常数。判别器损失Ld表达式为：Ld=ssunSsunDplanet(sunplanet(ssun)-splanetSplanetDplanet(splanet)+splanetSplanetDsun(planetsun(splanet)-ssunSsunDsun(ssun)（2）生成器损失Lg表达式为：Lg=-ssun SsunDplanet()sunplanet()ssun-splanet SplanetDsun()planetsun()splanet（3）将重构光谱S与真实光谱S的L1范数作为生成器重构损失Lgr，表达式为：Lgr=Lgr 1+Lgr 2=ssun Ssunssun-planetsun()sunplanet()ssun+splanet Splanet splanet-sunplanet()planetsun()splanet（4）式（1）（4）构成VAE-Dual GAN网络的总损失Ltotal，表达式为：Ltotle=LVAE+Ld+Lg+Lgr（5）式中：0，1为作用于生成器重构损失的常数，用来均衡生成器重构损失Lgr对网络的作用力，防止其造成的过拟合现象。1.2泛化性分析VAE-Dual GAN由2个互为对偶的子网络模型构成，两个VAE的作用是增广训练集。因此只需证明Dual GAN其中一个网络的泛化性即可。要证明泛化性，首先需要Lipschitz连续定义。定义1.对于任意2个样本x和y，如果损失函L(x)相对于距离度量满足：|L(x)-L(y)|L(x，y)（6）则称L(x)是 Lipschitz 连续。其中，L是 Lipschitz常数。在Dual GAN的子网络中，分布之间的Wasserstein距离表达式为：d(X，Y)=supf 1x Xf(x)-y Yf(y)（7）式中：f 1为连续函数，f满足1-Lipschitz连续；x和y分别为服从X和Y分布的样本。根据定义1，引入判别器函数f和生成器函数g后，由式（7）可以得到目标函数为：minU(f，g)=y Yf(g(y)-x Xf(x)（8）式（8）用来衡量判别器对真实样本和生成样本的区分质量。如果泛化效果好，随着样本数量的不断增加，无限分布与有限分布之间的差异会更小。利用给定的h个样本Xh=x1，x2，xh和Yh=y1，y2，yh，分别服从分布X和Y（分别是分布X和Y的有限分布），在固定生成器g的情况下，式（8）可以重写为：minUh(f，g)=1hi=1h()f()g()y-f()x（9）本文定义无限分布下和有限分布下的目标函数U和Uh为：U=minU(f，g)（10）Uh=minUh(f，g)（11）如果网络泛化成功，在输入样本增加到一定的数量后，两式的差值|Uh-U|会依概率收敛为0；如果泛化失败，它将不能依概率收敛为0，这也意味着生成的分布过拟合真实分布。文献 37 受GAN及其变体用有限的样本最小化真实分布X和生成分布Y的散度或距离的启发，对GAN的泛化性定义如下：定义2.对于生成分布Y，如果分布之间的距离或散度d(，)用含有h个训练样本的有限分布X和Y进行泛化，大概率式（6）成立，即：|d()X，Y-d()X，Y （12）则称GAN具有泛化性。根据定义1和定义2可以证明以下定理。定理 1.在与有界 Lipschitz 常数L有关的假设条276第 2 期刘劲等：面向天文多普勒差分测速的太阳/行星光谱对生成方法件下38：I.f(x)在参数上是Lipschitz连续，即对任意都有|f(x)-f(x)|L-；II.f(x)在样本x上是1-Lipschitz连续，即对任意x都有|f(x)-f(x)|x-x；III.两个样本的距离存在边界B，即|()x，xB。存在一个普适常数 b，有h baB22ln()La 2，使得：P(|Uh-U|)1-（13）式中：a为目标函数的参数数量。证.将式（10）（11）代入式（13）得式（14）：P(|Uh(f，g)-U(f，g)|)1-（14）因为需要用到 McDiarmid s不等式，所以必须界定样本发生变化前后目标函数的变化边界。本文用Uih(f，g)和Ujh(f，g)表示样本改变前后的目标函数。在假设条件下，结合式（9）可得式（15）：|Uih(f，g)-Ujh(f，g)|1h|f(g(yi)-f(g(yj)|+1h|f(xi)-f(xj)|1h()g()yi，g()yj+)()xi，xj2hB (15)需要说明：U(f，g)=s S m M i=1，hUh()f，g（16）由McDiarmid s不等式可以得到：P(|Uh(f，g)-U(f，g)|2)2exp()-2h82B2（17）接下来，通过构造一个()8L网络N37，将所有判别器函数f划分为2种：网络N内的f和网络N外的f。（1）在网络N内将式（15）表述的变化边界从适用于单个判别器函数f推广为联合边界，即对于网络N中的判别器函数f，总存在一个 N有-()8L。这个网络中所包含的损失函数满足|N|Q(aln(La)，其中a表示损失函数中的参数数量。因 此，对 于 所 有 的 N，当h baB22ln()La 2，有：P(|Uh(f，g)-U(f，g)|2)1-（18）（2）分析在网络N之外的判别器函数f的情况。|U(f，g)-U(f，g)|s S|f(g(y)-f(g(y)|+m M|f(x)-f(x)|2L-(19)同理可得：|Uh(f，g)-Uh(f，g)|2L-（20）因此，可以得出包括所有损失函数的联合边界：|Uh(f，g)-U(f，g)|2L-+2+2L-4+2+4=（21）所以有：P(|Uh(f，g)-U(f，g)|)1-（22）设最优判别器f使U(f，g)最小化，但不能保证可以最小化Uh(f，g)，所以有Uh Uh(f，g)。由式（22）的结论得出，有h baB22ln()La 2，以大于1-的概率使得：Uh-U Uh(f，g)-U(f，g)（23）定理1得证。这表明了VAE-Dual GAN有一定的泛化性。2天体光谱测速多普勒效应是由于波源与观测点的相对运动，使观测光谱波长与原波长存在频移。在只考虑多普勒效应下，航天器观测到的行星光谱波长p与太阳光谱波长s的关系如下：1p=1s(1+vc)（24）式中：c为光速；v为径向速度。对式（24）取对数，可以得到对数波长偏移为：277宇航学报第 45 卷=lns-lnp=ln(1+vc)（25）设太阳光谱ssun(lnn)和行星光谱splanet(lnn)的关系表达式为：splanet(lnn)=(lnn)ssun(lnn-)（26）式中：(lnn)为未知的反射率函数；n为离散波长。由于光谱的波长间隔不均匀，而原离散傅立叶变换无法处理非均匀间隔波长的光谱，所以采用非均匀离散傅立叶变换（NDFT）方法，表达式为：w(k)=n=0N-1w()lnnexp(-j2Lklnn)lnn（27）式中：L=lnmax-lnmin为对数波长的长度；k=1，2，K。式（27）经过NDFT后，splanet(lnn)和ssun(lnn)的幅值Splanet(k)与Ssun(k)的累积误差为：W2()=k=1N 21k|Splanet()k-()k exp()j()kSsun(k)exp j(k)-(k)+2kL)2（28）式中：k为傅立叶频率下噪声振幅的平方根；(k)和(k)分别为splanet(lnn)和ssun(lnn)的相位；(k)和(k)分别为(lnn)经过NDFT的幅值和相位；为循环卷积符号。当太阳/行星光谱对同步时，行星光谱仅受反射率影响；当太阳/行星光谱对不同步时，除反射率外，行星光谱还受太阳活动、太阳自转等因素影响，(lnn)的表达式复杂。因此，同步的太阳/行星光谱对可有效减小干扰因素的影响，提高测速精度。3实验结果与分析3.1实验参数初始学习率设为0.000 1，并且每100次迭代衰减1.5%。隐藏层使用Leaky ReLUs作为激活函数，泄露值设为 0.2。网络优化器为 RMSProp（Root mean square prop），并将每次更新的判别器梯度的绝对值修正到不超过截断常数，其值为0.1，优化器衰减速率为0.9。实验平台为宏碁笔记本，其处理器主频为3.10 GHz，内存为16 GB。3.2数据集本文数据集为ESO的HARPS观测的天文光谱，该仪器覆盖的光谱波长范围为 380690 nm，仪器的探测器由 2 个电荷耦合器件（CCD）组成，每个CCD分别可以读取到36个波长区间的光谱，共72段光谱。相较于第2个CCD，第1个CCD读取到的光谱的信噪比较低。在第2个CCD的35个波长区间的光谱中选择光谱进行测速。实验太阳和月球光谱的 4 个波长范围，包括：5 964.872-6 032.103；6 337.527-6 408.976；6 541.882-6 615.640；6 685.603-6 760.982。它们各有特点：5 964.872-6 032.103区间信噪比高而RV标准差低；6 337.527-6 408.976区间信噪比高且 RV 标准差高；6 541.882-6 615.640 区间是CHASE 工程的模板光谱；6 685.603-6 760.982区间是随机选择的，其信噪比高而RV标准差介于5 964.872-6 032.103 与 6 337.527-6 408.976之间。3.3光谱测速与相似度为了验证VAE-Dual GAN的可行性，将3种天文光谱的测速结果进行对比，如表1所示。3种天文光谱包括：生成的伪月球光谱，加高斯噪声的月球光谱，加泊松噪声的月球光谱。本文用Evr表示真实月球光谱径向速度标准差；Evp和Evg分别表示引入泊松噪声和高斯噪声的光谱测速标准差。表1显示，VAE-Dual GAN生成的伪月球光谱径向速度标准差Evgan低于Evr，高于Evp和Evg，表明噪声不是影响光谱测速精度的唯一因素。生成的同步光谱测速精度高于非同步光谱测速精度。为了说明真实太阳光谱和生成的伪月球光谱的时间同步性，本文引入了结构相似指标度量39（Structure similarity index measure，SSIM）。结构相似指标用来衡量两幅图像的相似度，计算公式为：VSSIM(Ssun，Splanet)=2sunplanet2sun_planet()2sun+2planet()2sun+2planet（29）式中：sun和planet分别为太阳光谱和行星光谱的均值；2sun和2planet分别为二者的方差；sun_planet为太阳光表1径向速度标准差Table 1Standard deviation of radial velocity波长/（）5 964.872-6 032.1036 337.527-6 408.9766 541.882-6 615.6406 685.603-6 760.982Evr11.473 6100.282 146.693 229.516 2Evp4.054 84.600 33.906 46.041 0Evg3.654 74.287 54.266 36.220 1Evgan4.146 76.536 218.308 68.363 7278第 2 期刘劲等：面向天文多普勒差分测速的太阳/行星光谱对生成方法谱与行星光谱的协方差。VSSIM的取值范围为0，1。VSSIM越大，表示二者相似度越高。表2给出了4种天文光谱对的VSSIM值，包括：真实太阳/真实月球光谱对（Real sun/Real moon，Rs/Rm）、真实太阳/伪月球光谱对（Real sun/Fake moon，Rs/Fm）、真实太阳/泊松噪声月球光谱对（Real sun/Poisson noise moon，Rs/Pm）、真实太阳/高斯噪声月球光谱对（Real sun/Gaussian noise moon，Rs/Gm）。表 2 可以看出，对于每个波长区间，Rs/Rm 的SSIM 低于 Rs/Fm，高于 Rs/Pm 和 Rs/Gm。这说明伪月球光谱与真实太阳光谱的SSIM高于真实月球光谱与真实太阳光谱的SSIM。太阳表面活动引起了光谱的变化。时间不同步造成真实月球光谱与真实太阳光谱的 SSIM 低。Rs/Fm 的 SSIM 高，这说明二者的时间是同步的。此外，引入噪声会降低光谱对的SSIM。这表明加噪方法无法得到高SSIM的光谱对，不是生成同步光谱对的有效手段。3.4网络层数为了考察网络层数对伪光谱质量的影响，本文在2层全连接层网络的基础上，将VAE编码器和生成器分别增加一层全连接层，并增加二者的神经元数量。VAE编码器的3层神经元数量分别为500，300，200；生成器的3层神经元数量分别为200，300，500。本节比较了 2 层 VAE-Dual GAN 和 3 层 VAE-Dual GAN的Rs/Fm的VSSIM值以及训练时间，结果如表 3所示，对于每个波长区间，3层 VAE-Dual GAN的Rs/Fm高于Rs/Rm，但低于2层VAE-Dual GAN的Rs/Fm。这说明 2 层 VAE-Dual GAN 模型生成的伪月球光谱质量高于3层的VAE-Dual GAN模型。究其原因，过多的全连接层和神经元数将造成过拟合，进而导致伪光谱质量的下降。此外，表3显示，2层VAE-Dual GAN 模型的训练时间比 3 层 VAE-Dual GAN模型短约10 000 s。这说明增加全连接层数和神经元数均会使训练时间延长。3.5卷积层与全连接层的对比为了体现全连接层的优越性，本文将其与卷积层进行对比。基于卷积层的VAE-Dual GAN的网络结构具体如下：VAE编码器为2层卷积层，生成器为5层转置卷积层和1层卷积层；卷积核大小为3 3，通道数为32，步幅为2。表4给出了2种VAE-Dual GAN的Rs/Fm的VSSIM值以及伪光谱的径向速度标准差。卷积层 VAE-Dual GAN的Rs/Fm略高于Rs/Rm和全连接层VAE-Dual GAN的Rs/Fm。卷积层VAE-Dual GAN的径向表2结构相似性（SSIM）Table 2Structural similarity index measure（SSIM）波长/（）5 964.872-6 032.1036 337.527-6 408.9766 541.882-6 615.6406 685.603-6 760.982VSSIM（Rs/Rm）0.903 0760.934 2730.919 6800.946 332VSSIM（Rs/Fm）0.906 9240.940 1260.921 6730.948 154VSSIM（Rs/Pm）0.902 6930.933 7980.919 3100.945 757VSSIM（Rs/Gm）0.902 7100.933 7990.919 3120.945 751表 3网络层数Table 3The number of network layers波长/（）5 964.872-6 032.1036 337.527-6 408.9766 541.882-6 615.6406 685.603-6 760.982VSSIMRs/Rm0.903 0760.934 2730.919 6800.946 332Rs/Fm（2层）0.906 9240.940 1260.921 6730.948 154Rs/Fm（3层）0.906 8220.937 8720.919 8190.947 599训练时间/s2层24 4443层34 194表4卷积层与全连接层Table 4Convolution layer and fully connection layer波长/（）5 964.872-6 032.1036 337.527-6 408.9766 541.882-6 615.6406 685.603-6 760.982VSSIMRs/Rm0.9030.9340.9200.946Rs/Fm（全连接层）0.9070.9400.9220.948Rs/Fm（卷积层）0.9080.9400.9220.948径向速度标准差Evr11.474100.28246.69329.516Evgan（全连接层）4.1476.53618.3098.364Evgan（卷积层）10.40343.38624.34393.301279宇航学报第 45 卷速度标准差高于全连接VAE-Dual GAN的径向速度标准差。究其原因，卷积层VAE-Dual GAN仅关注光谱特征，而非相位。这导致卷积层 VAE-Dual GAN Rs/Fm的结构相似度高于全连接层VAE-Dual GAN，测速精度却低于全连接层VAE-Dual GAN。3.6网络指标对比为了验证VAE-Dual GAN优于Dual GAN和VAE-Cycle GAN，本文采用SSIM和标准化相对鉴别分数40（Normalized relative discriminative score，NRDS）作为对比指标。光谱波长为 5 964.872-6 032.103。表5给出了生成的伪行星光谱与真实行星光谱的结构相似度。VAE-Dual GAN生成的伪月球光谱与真实太阳光谱的相似度最高，Rs/Rm最低，这说明3个模型的伪光谱对在相似度上都有不同程度的提升，并且VAE-Dual GAN 的伪光谱优于 VAE-Cycle GAN 和Dual GAN。NRDS是评判GAN模型优劣的指标。将相同个数的真实光谱和多个GAN模型生成的伪光谱混合投入一个分类器。经过足够多的训练迭代次数，使分类器可以完全区分出2类光谱样本。当生成伪光谱的概率分布与真实分布高度接近时，分类器很难完全区分。NRDS计算公式为：TNRDS(Ci)=A()CiA()Ci（30）式中：A（Ci）为第i个GAN模型输出曲线Ci围成的面积。TNRDS越高，对应的GAN模型越好。下面，考察 VAE-Dual GAN 与 Dual GAN 和VAE-Cycle GAN 的TNRDS，光谱波长为 5 964.872-6 032.103。图3给出了各种GAN模型生成的错误率曲线，表6给出了其对应的TNRDS。从以上结果可以看出，VAE-Dual GAN 的错误率曲线高于另 2 个GAN 模型，且 VAE-Dual GAN 的TNRDS最大。VAE-Dual GAN的伪光谱的区分难度最高，即该模型生成的光谱对质量优于Dual GAN和VAE-Cycle GAN。4结论天文多普勒差分导航需要同步的太阳光谱和行星光谱，但目前无法满足要求。为此，本文提出一种基于深度学习的太阳/行星光谱对生成方法。为充分利用VAE的数据增强能力和Dual GAN的转换学习能力，将二者融合，以生成同步的太阳/伪行星光谱对。实验结果表明，VAE 具有数据增强能力，解决了小样本情况下约束能力较弱的问题。VAE-Dual GAN的太阳/伪行星光谱对的测速精度和结构相似度均优于真实光谱对、含噪光谱对以及其他GAN生成的光谱对。这说明了VAE-Dual GAN的太阳/伪行星光谱对的质量高。参 考 文 献1 汪中生，孟占峰，高珊，等.嫦娥五号月球轨道交会对接远程导引轨道设计与飞行实践 J.宇航学报，2021，42（8）：939-952.WANG Zhongsheng，MENG Zhanfeng，GAO Shan，et al.Phasing orbit design and flight test of lunar orbit RVD in Chang E-5 mission J.Journal of Astronautics，2021，42（8）：939-952.2 饶炜，孙泽洲，董捷，等.天问一号火星进入、下降与着陆系统设计与实现 J.中国科学：技术科学，2022，52（8）：1162-1174.RAO Wei，SUN Zezhou，DONG Jie，et al.Design and implementation of the Mars entry，descent，and landing system for the Tianwen-1 mission J.Scientia Sinica（Technologica），表5GAN模型的VSSIMTable 5VSSIM of GAN models模型真实光谱（Rs/Rm）Dual GAN 生成光谱VAE-Cycle GAN 生成光谱VAE-Dual GAN 生成光谱VSSIM0.903 075 9250.903 966 5790.904 957 1980.906 923 56900.51.01.52.02.53.0迭代次数00.10.20.30.40.50.6错误率Dual GANVAE-Dual GANVAE-Cycle GAN104图3GAN模型的错误率曲线Fig.3Error rate curves of GAN models表6GAN模型的TNRDSTable 6TNRDS of GAN modelsGAN 模型VAE-Dual GANDual GANVAE-Cycle GANTNRDS0.392 70.315 70.291 6280第 2 期刘劲等：面向天文多普勒差分测速的太阳/行星光谱对生成方法2022，52（8）：1162-1174.3 ZHANG W，CHENG W Q，YOU W，et al.Levitated-body ultra-high pointing accuracy and stability satellite platform of the CHASE missionJ.Science China Physics，Mechanics&Astronomy，2022，65（8）：289604.4 NING X L，GUI M Z，FANG J C，et al.Differential X-ray pulsar aided celestial navigation for Mars explorationJ.Aerospace Science and Technology，2017，62：36-45.5 WANG Y D，ZHENG W，SUN S M，et al.X-ray pulsar-based navigation using time-differenced measurementJ.Aerospace Science and Technology，2014，36：27-35.6 陈晓，许俊，张晓，等.深空天文测角测速组合导航滤波算法研究 J.中国科学：物理学 力学 天文学，2022，52（1）：76-82.CHEN Xiao，XU Jun，ZHANG Xiao，et al.Research on the filter algorithm for integrated autonomous navigation based