低频射电频谱仪观测目标源流量密度反演处理方法技术

本发明专利技术公开了一种低频射电频谱仪观测目标源流量密度反演处理方法，包括：对低频射电频谱仪接收到的科学数据进行预处理；将预处理后的数据分为目标源爆发前与目标源爆发时两组，并对目标源爆发前的数据进行平均处理；用目标源爆发时数据的信号频谱减去平均处理后的目标源爆发前数据的信号频谱，计算得到降噪后的观测数据；利用低频射电频谱仪的定标参数将观测数据反演成空间电场强度；利用空间电场强度计算出目标源的流量密度。本发明专利技术提供的低频射电频谱仪观测目标源流量密度反演处理方法，结合低频射电频谱仪收到的科学数据，并结合改进的谱减法减去天线接收到的信号中的着陆器噪声，反演出目标源的空间电场强度以及目标源的流量密度。

Inversion Processing Method of Flow Density of Observed Target Source by Low Frequency Radio Spectrometer

【技术实现步骤摘要】
低频射电频谱仪观测目标源流量密度反演处理方法
本专利技术涉及数据处理
，具体涉及一种将低频射电频谱仪的观测数据反演目标源流量密度的方法，用于将低频射电频谱仪的观测数据反演成目标源流量密度的方法。
技术介绍
地球大气的存在使得很多射电波段的信号特别是频率低于10MHz的甚低频波段难以穿透到地面，因此在地面进行低频射电观测非常困难。为了弥补甚低频波段的观测空白，国际上许多国家或组织发射过卫星到空间进行观测，比如WIND和Ulysses。由于地球以及众多在轨卫星的电磁干扰，空间低频射电观测也面临挑战。因为月球自身可以遮挡来自地球方向的射电干扰，所以月球背面的射电宁静区是十分理想的观测地点。欧航局(ESA)于2015年1月提出了在2025年对月球背面发射一颗载有低频射电频谱仪的着陆器的FARSIDE计划。我国提出嫦娥四号探月计划。嫦娥四号是世界首颗在月球背面软着陆和巡视探测的航天器。2018年5月21日，嫦娥四号中继星“鹊桥”成功发射，着陆器和巡视器已经在2018年12月8号发射并于2019年1月3号成功着陆，该着陆器上搭载了一台科学载荷低频射电频谱仪，是人类在月球背面首次实现对宇宙和太阳的低频段(0.1MHz～40MHz)进行观测。根据着陆器整器在微波暗室的电磁兼容性测试(EMC)试验结果，虽然来自外界的信号干扰可以被很好的屏蔽掉，但是着陆器本身搭载的电子设备以及其它的载荷对低频观测存在电磁干扰，由于无法对平台进行工程上的调整，所以本专利技术利用改进的谱减法对低频射电频谱仪观测数据进行了降噪处理并反演成可以直接用于科学研究的空间电场强度。
技术实现思路
(一)要解决的技术问题有鉴于此，本专利技术的目的在于提供一种低频射电频谱仪观测目标源流量密度反演处理方法，有效降低嫦娥四号着陆器的平台噪声，并将低频射电频谱仪的观测数据反演成可以用于科学研究的目标源的流量密度。(二)技术方案本专利技术提供了一种低频射电频谱仪观测目标源流量密度反演处理方法，包括：对低频射电频谱仪接收到的科学数据进行预处理；将预处理后的数据分为目标源爆发前与目标源爆发时两组，并对目标源爆发前的数据进行平均处理；用目标源爆发时数据的信号频谱减去平均处理后的目标源爆发前数据的信号频谱，计算得到降噪后的观测数据；利用低频射电频谱仪的定标参数将观测数据反演成空间电场强度；利用空间电场强度计算出目标源的流量密度。其中，所述对低频射电频谱仪接收到的科学数据进行预处理的步骤中预处理包括：科学数据提取、数据拼接、物理量转换和傅里叶变换。具体地，所述科学数据以二进制存储，以包为单位，每包数据包括8组，每组数据包含若干道，每道数据包含了道头和通道科学数据；所述数据拼接是将每道数据中道头去掉，并将同一组的科学数据拼接在一起；所述物理量转换采用如下公式实现，式中，UAD(t)为物理量转换前的科学数据，Uout(t)为物理量转换后的科学数据；所述傅里叶转换是对每一组物理量转换后的科学数据进行傅里叶变换，变换公式为式中，Uout(t)为物理量转换前的科学数据；Uout(f)为物理量转换后的科学数据；N为采样点数。其中，所述将预处理后的数据分为目标源爆发前与目标源爆发时两组的步骤包括：将预处理后的数据分为目标源爆发前的数据以及目标源爆发时的数据两组，所述目标源为太阳；对目标源爆发前的数据在频域上对对应频点的幅值进行平均处理；其中，所述用目标源爆发时数据的信号频谱减去目标源爆发前数据的信号频谱，计算得到降噪后的观测数据的步骤，是采用以下公式：式中，a为过减因子；V′(f)为目标源爆发前信号；V(f)目标源正在爆发时的信号；U(f)为降噪后的观测数据；b(b＞0)是补偿因子；其中，所述利用低频射电频谱仪的定标参数将观测数据反演成空间电场强度的步骤中，将观测数据反演成空间电场强度所采用的反演公式为：式中，是空间中低频射电信号电场；he是天线的有效长度；是天线的归一化接收方向图；是天线A电场极化方向的单位矢量；Za为天线阻抗；ZPreamp为前置放大器的输入阻抗；GPreamp为前置放大器的放大倍数；ILCoax为高频电缆损耗；GReceiver为接收机增益；U为降噪后的观测数据。其中，所述利用空间电场强度计算出目标源的流量密度的步骤中，利用电磁波电场强度计算目标源流量密度的公式为：式中，ε和μ分别是电介质常数和磁介质常数；为利用低频射电频谱仪的定标参数将观测数据反演得到的空间电场强度；为目标源的流量密度；为电磁波传播方向上的单位矢量。(三)有益效果从上述技术方案可以看出，本专利技术提供的低频射电频谱仪观测目标源流量密度反演处理方法具有以下的有益效果：(1)本专利技术提供的低频射电频谱仪观测目标源流量密度反演处理方法，在谱减的过程中使用信号的幅值或者功率谱，并采用平均谱值进行计算噪声的谱减，优化了谱减方法；(2)本专利技术提供的低频射电频谱仪观测目标源流量密度反演处理方法，结合低频射电频谱仪搭载的三根长天线接收到的科学数据，并结合改进的谱减法减去天线接收到的信号中的着陆器噪声，反演出目标源的空间电场强度以及目标源的流量密度。附图说明图1是本专利技术低频射电频谱仪观测目标源流量密度反演处理方法的流程图；图2是本专利技术实施例中5MHz点频信号的试验数据降噪处理结果图；图3是本专利技术实施例中22.5MHz点频信号的试验数据降噪处理结果图；图4是本专利技术实施例中1.75MHz宽带信号的试验数据降噪处理结果图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本专利技术进一步详细说明。本专利技术提供了一种低频射电频谱仪观测目标源流量密度反演的处理方法。图1是该处理方法的流程图，参阅图1，该方法具体包括：本专利技术提供了一种低频射电频谱仪观测目标源流量密度反演处理方法，包括：对低频射电频谱仪接收到的科学数据进行预处理；将预处理后的数据分为目标源爆发前与目标源爆发时两组，并对目标源爆发前的数据进行平均处理；用目标源爆发时数据的信号频谱减去平均处理后的目标源爆发前数据的信号频谱，计算得到降噪后的观测数据；利用低频射电频谱仪的定标参数将观测数据反演成空间电场强度；利用空间电场强度计算出目标源的流量密度。其中，对低频射电频谱仪接收到的科学数据进行预处理的步骤中预处理包括：科学数据提取、数据拼接、物理量转换和傅里叶变换，具体地：科学数据以二进制存储，以包为单位，每包数据包括8组，每组数据包含若干道，每道数据包含了道头和通道科学数据；数据拼接是将每道数据中道头去掉，并将同一组的科学数据拼接在一起；物理量转换采用如下公式实现，式中，UAD(t)为物理量转换前的科学数据，Uout(t)为物理量转换后的科学数据；傅里叶转换是对每一组物理量转换后的科学数据进行傅里叶变换，变换公式为式中，Uout(t)为物理量转换前的科学数据；Uout(f)为物理量转换后的科学数据；N为采样点数。其中，将预处理后的数据分为目标源爆发前与目标源爆发时两组的步骤包括：将预处理后的数据分为目标源爆发前的数据以及目标源爆发时的数据两组，所述目标源为太阳；对目标源爆发前的数据在频域上对对应频点的幅值进行平均处理；其中，用目标源爆发时数据的信号频谱减去目标源爆发前数据的信号频谱，计算得到降噪后的观测数据的步骤，是采用以下公式实现：式中本文档来自技高网...

【技术保护点】
1.一种低频射电频谱仪观测目标源流量密度反演处理方法，包括：对低频射电频谱仪接收到的科学数据进行预处理；将预处理后的数据分为目标源爆发前与目标源爆发时两组，并对目标源爆发前的数据进行平均处理；用目标源爆发时数据的信号频谱减去平均处理后的目标源爆发前数据的信号频谱，计算得到降噪后的观测数据；利用低频射电频谱仪的定标参数将观测数据反演成空间电场强度；利用空间电场强度计算出目标源的流量密度。

【技术特征摘要】
1.一种低频射电频谱仪观测目标源流量密度反演处理方法，包括：对低频射电频谱仪接收到的科学数据进行预处理；将预处理后的数据分为目标源爆发前与目标源爆发时两组，并对目标源爆发前的数据进行平均处理；用目标源爆发时数据的信号频谱减去平均处理后的目标源爆发前数据的信号频谱，计算得到降噪后的观测数据；利用低频射电频谱仪的定标参数将观测数据反演成空间电场强度；利用空间电场强度计算出目标源的流量密度。2.根据权利要求1所述的处理方法，其特征在于，所述对低频射电频谱仪接收到的科学数据进行预处理的步骤中预处理包括：科学数据提取、数据拼接、物理量转换和傅里叶变换。3.根据权利要求1或2所述的方法，其特征在于，所述科学数据以二进制存储，以包为单位，每包数据包括8组，每组数据包含若干道，每道数据包含了道头和通道科学数据。4.根据权利要求3所述的方法，其特征在于，所述数据拼接是将每道数据中道头去掉，并将同一组的科学数据拼接在一起。5.根据权利要求2所述的方法，其特征在于，所述物理量转换采用的公式为其中，UAD(t)为物理量转换前的科学数据，Uout(t)为物理量转换后的科学数据。6.根据权利要求5所述的方法，其特征在于，所述傅里叶转换是对每一组物理量转换后的科学数据进行傅里叶变换，变换公式为其中，Uout(t)为物理量转换前的科学数据；Uout(f)为物理量转换后的科学数据；N为采样点数。7.根据权利要求1所述的处理方法，其特征在于，所述将预处理...

【专利技术属性】
技术研发人员：苏彦，朱新颖，邢树果，孔德庆，李俊铎，薛喜平，戴舜，李臣，李春来，张洪波，张韬，
申请(专利权)人：中国科学院国家天文台，
类型：发明
国别省市：北京,11