10~50 MHz地基低频射电天文数字接收机及上位机软件研制

许丙强, 白宇, 路光, 平劲松, 王明远, 杨文军, 王震, 严发宝. 2023. 10~50 MHz地基低频射电天文数字接收机及上位机软件研制. 地球物理学报, 66(3): 891-904, doi: 10.6038/cjg2022Q0015
引用本文: 许丙强, 白宇, 路光, 平劲松, 王明远, 杨文军, 王震, 严发宝. 2023. 10~50 MHz地基低频射电天文数字接收机及上位机软件研制. 地球物理学报, 66(3): 891-904, doi: 10.6038/cjg2022Q0015
XU BingQiang, BAI Yu, LU Guang, PING JingSong, WANG MingYuan, YANG WenJun, WANG Zhen, YAN FaBao. 2023. Development of 10~50 MHz ground based low frequency radio astronomical digital receiver and PC software. Chinese Journal of Geophysics (in Chinese), 66(3): 891-904, doi: 10.6038/cjg2022Q0015
Citation: XU BingQiang, BAI Yu, LU Guang, PING JingSong, WANG MingYuan, YANG WenJun, WANG Zhen, YAN FaBao. 2023. Development of 10~50 MHz ground based low frequency radio astronomical digital receiver and PC software. Chinese Journal of Geophysics (in Chinese), 66(3): 891-904, doi: 10.6038/cjg2022Q0015

10~50 MHz地基低频射电天文数字接收机及上位机软件研制

  • 基金项目:

    国家自然科学基金(41904158,11941003),中科院与荷兰科技组织射电天文国际合作项目(GJHZ1867),嫦娥四号月球着陆探测NCLE载荷专项(2016年)和山东大学(威海)青年学者未来计划(20820201005)联合资助

详细信息
    作者简介:

    许丙强, 男, 1996年生, 硕士研究生, 主要从事低频段射电天文接收机研究.E-mail: xubingqiangwy@163.com

    通讯作者: 严发宝, 男, 1985年生, 副教授, 博士生导师, 主要从事微波探测、射电观测与系统关键技术研究.E-mail: yanfabao2022@163.com
  • 中图分类号: P353

Development of 10~50 MHz ground based low frequency radio astronomical digital receiver and PC software

More Information
  • 低频段(< 50 MHz)射电观测,是研究射电爆发、恒星形成、宇宙初期状态的重要手段,而我国对该频段相关仪器研究较少.本文从关键器件选型、仿真与设计、接口和通讯、功能模块等方面,对组成地基低频射电观测系统的核心部件——数字接收机与上位机软件进行论述.数字接收机由采集卡和工控机构成,采集卡配置包含AD9265芯片和K7系列现场可编程门阵列(Field Programmable Gate Array,FPGA).数字接收机探测频率覆盖1~62.5 MHz,接收功率范围-78~-6 dBm,拥有先进的频率分辨率(15.3 kHz)和时间分辨率(0.5~32 ms).为满足射电观测显示与存储需求,研发了具有针对性的上位机软件解决方案.软件支持高速串行计算机扩展总线标准(Peripheral Component Interconnect express,PCIe)数据传输协议,支持实时频谱显示、数据存储和交互操作,具有全频带内强度图显示与图像自动保存功能.自主研制的数字接收机与上位机软件已成功应用于新疆奇台观测站低频(10~50 MHz)射电观测设备,该设备自2021年6月投入运行,捕获多次太阳射电暴事件.将观测结果与澳大利亚Learmonth观测设备对比分析,结果表明,新疆奇台站低频射电观测系统探测精度更为精细、探测频率更低,性能指标位于国际前列.系统后期可开展星-地协同观测,推动我国低频段射电研究.

  • 加载中
  • 图 1 

    数字接收机与上位机软件位置示意图

    Figure 1. 

    Digital receiver and upper computer software position diagram

    图 2 

    采集卡主要结构

    Figure 2. 

    Main structure of acquisition card

    图 3 

    MATLAB仿真ADC采样数据图片

    Figure 3. 

    MATLAB simulation ADC sampling data picture

    图 4 

    ADC无杂散动态范围测试图

    Figure 4. 

    ADC spurious free dynamic range test diagram

    图 5 

    FPGA内部数据处理步骤框图

    Figure 5. 

    Block diagram of FPGA internal data processing steps

    图 6 

    MATLAB仿真数据加入Hanning窗之后的影响

    Figure 6. 

    Influence of adding Hanning window to MATLAB simulation data

    图 7 

    数字接收机相关接口连接示意图

    Figure 7. 

    Digital receiver correlation interface connection diagram

    图 8 

    软件设计

    Figure 8. 

    Software design

    图 9 

    DMA数据通讯示意图

    Figure 9. 

    DMA data communication diagram

    图 10 

    数据采集线程处理步骤

    Figure 10. 

    Procedure for processing data acquisition threads

    图 11 

    上位机软件主界面

    Figure 11. 

    Main interface of upper computer software

    图 12 

    上位机软件强度图界面

    Figure 12. 

    Upper computer software strength chart interface

    图 13 

    数据存储线程

    Figure 13. 

    Data storage thread

    图 14 

    数据标定示意图

    Figure 14. 

    Data calibration diagram

    图 15 

    澳大利亚观测设备提供的观测数据,右侧色标单位为(dBm/15.3 kHz),数据来自World Data Center,https://downloads.sws.bom.gov.au/wdc/wdc_spec/data/learmonth/raw/20/

    Figure 15. 

    Observation data provided by Australian observation equipment, right colorbar unit (dBm/15.3 kHz), data from World Data Center, https://downloads.sws.bom.gov.au/wdc/wdc_spec/data/learmonth/raw/20/

    图 16 

    本系统提供的观测数据,右侧色标单位为(dBm/15.3 kHz)

    Figure 16. 

    Observation data provided by the system, right colorbar unit (dBm/15.3 kHz)

    图 17 

    本系统提供的观测数据,右侧colorbar单位为(dBm/15.3 kHz)

    Figure 17. 

    Observation data provided by the system, right colorbar unit (dBm/15.3 kHz)

    表 1 

    国内外部分低频段射电观测设备及参数

    Table 1. 

    Some low-frequency radio observation equipment and parameters at home and abroad

    名称 观测频段 时间分辨率 频率分辨率
    LOFAR 10~80 MHz 1 s 0.76 kHz
    LWA1 24~87 MHz 0.84 ms 4.79 kHz
    Learmonth 25~180 MHz 3 s ~128 kHz
    NDA 10~120 MHz 3 ms 35 kHz
    明安图VLF 30~70 MHz ≥10 ms ~39 kHz
    山东大学CSO 150~500 MHz (扩展了25~110 MHz) 1~10 ms 3~320 kHz
    奇台低频射电观测设备 10~50 MHz 0.5~32 ms 15.3 kHz
    下载: 导出CSV

    表 2 

    AD9265主要参数

    Table 2. 

    Main parameters of AD9265

    指标 参数 指标 参数
    采样率 ≤125 MSPS 输入带宽 ≤300 MHz
    位数 16 bit SNR 79.0 dBFS
    通道数 2 模拟输入 < ±1 V
    下载: 导出CSV

    表 3 

    XC7K410T部分资源

    Table 3. 

    XC7K410T partial resources

    指标 参数 指标 参数
    Logic Cells 406000 DSP Slices 1540
    CLB Slices 63550 Block RAM (Max) 28620 kb
    Max Distributed RAM 5663 kb I/O Pin 500
    下载: 导出CSV
  •  

    Bentum M J, Verma M K, Rajan R T, et al. 2020. A roadmap towards a space-based radio telescope for ultra-low frequency radio astronomy. Advances in Space Research, 65(2): 856-867, doi: 10.1016/j.asr.2019.09.007.

     

    Chen L J, Yan Y H, Fan Q X, et al. 2021. An agile very low frequency radio spectrum explorer. Research in Astronomy and Astrophysics, 21(4): 085, doi: 10.1088/1674-4527/21/4/85.

     

    Du G H, Kong X L, Chen Y, et al. 2015. An observational revisit of band-split solar type-ii radio bursts. The Astrophysical Journal, 812(1): 52, doi: 10.1088/0004-637X/812/1/52.

     

    Du Q F, Chen L, Zhao Y C, et al. 2017. A solar radio dynamic spectrograph with flexible temporal-spectral resolution. Research in Astronomy and Astrophysics, 17(9): 098, doi: 10.1088/1674-4527/17/9/98.

     

    Feng S W, Chen Y, Li C Y, et al. 2018. Harmonics of solar radio spikes at metric wavelengths. Solar Physics, 293(3): 39, doi: 10.1007/s11207-018-1263-z.

     

    Jester S, Falcke H. 2009. Science with a lunar low-frequency array: From the dark ages of the Universe to nearby exoplanets. New Astronomy Reviews, 53(1-2): 1-26, doi: 10.1016/j.newar.2009.02.001.

     

    Kaiser M L. 2005. The STEREO mission: an overview. Advances in Space Research, 36(8): 1483-1488, doi: 10.1016/j.asr.2004.12.066.

     

    Kennewell J A, Cornelius D W. 1983. Learmonth solar observatory. Aust. Phys. , 20(11): 276-280.

     

    Lecacheux A. 2000. The Nançay Decameter Array: A useful step towards giant, new generation radio telescopes for long wavelength radio astronomy. //Stone R G, Weiler K W, Goldstein M L, et al eds. Radio Astronomy at Long Wavelengths. Washington, DC: American Geophysical Union, 119: 321-328, doi: 10.1029/GM119p0321.

     

    Liu W Z, Wang D H. 2020. Influence of front-end circuit design on SFDR of ADC. Shipboard Electronic Countermeasure (in Chinese), 43(5): 108-111, doi: 10.16426/j.cnki.jcdzdk.2020.05.025.

     

    Mei L, Su Y, Zhou J F. 2018. The history and development of the low-frequency radio observation. Astronomical Research & Technology (in Chinese), 15(2): 127-139, doi: 10.14005/j.cnki.issn1672-7673.20180103.012.

     

    Peng L. 2008. Calibration precision analysis for spaceborne microwave radiometer [master′s thesis] (in Chinese). Beijing: University of Chinese Academy of Sciences.

     

    Ping J S, Wang M Y, Zhang M, et al. 2021. Introduction of space exploration progress for planetary radio burst emission. Journal of Deep Space Exploration (in Chinese), 8(1): 80-91, doi: 10.15982/j.issn.2096-9287.2021.20200016.

     

    Shang Z Q, Xu K, Liu Y, et al. 2022. A broadband solar radio dynamic spectrometer working in the millimeter-wave band. The Astrophysical Journal Supplement Series, 258(2): 25, doi: 10.3847/1538-4365/AC4257.

     

    Shi F, Xie Z Y, Yin Y N, et al. 2015. Research on fast discrete fourier transform algorithm and FPGA implementation. Electrical Measurement & Instrumentation (in Chinese), 52(5): 15-19, doi: 10.3969/j.issn.1001-1390.2015.05.004.

     

    Taylor G B, Ellingson S W, Kassim N E, et al. 2012. First light for the first station of the long wavelength array. Journal of Astronomical Instrumentation, 1(1): 1250004, doi: 10.1142/S2251171712500043.

     

    Tian H, Qu Z Q, Chen Y J, et al. 2017. Observations of the solar corona during the total solar eclipse on 21 August 2017. Earth and Planetary Physics, 1(1): 68-71. doi: 10.26464/epp2017010.

     

    Van Haarlem M P, Wise M W, Gunst A W, et al. 2013. LOFAR: The low-frequency array. Astronomy & Astrophysics, 556: A2, doi: 10.1051/0004-6361/201220873.

     

    Wang Y, Dong Z, Lu G, et al. 2022a. Design of a high-resolution antenna array solar observing system for radio frequencies of 25-110 MHz. Publications of the Astronomical Society of Australia, Volume 39, e030. doi: https://doi.org/10.1017/pasa.2022.27.

     

    Wang Y M, Jia X Z, Wang C B, et al. 2020. Locating the source field lines of Jovian decametric radio emissions. Earth Planet. Phys. , 4(2), 95-104 doi: 10.26464/epp2020015.

     

    Wang Y M, Zheng R B, Jia X Z, et al. 2022b. Reply to Comment by Lamy et al. on "Locating the source field lines of Jovian decametric radio emissions". Earth and Planetary Physics, 6(1), 13-17, doi: 10.26464/epp2022019.

     

    Xu K, Shang Z Q, Yan F B, et al. 2021. Compensation method of signal flatness for a broadband solar millimeter radio observation system. Scientia Sinica Technologica (in Chinese), 51(4): 413-423, doi: 10.1360/SST-2020-0283.

     

    Yan F B, Liu Y, Xu K, et al. 2020. A broadband digital receiving system with large dynamic range for solar radio observation. Research in Astronomy and Astrophysics, 20(9): 156, doi: 10.1088/1674-4527/20/9/156.

     

    Yan F B, Liu Y, Xu K, et al. 2021. Study of the truncation strategy in the FPGA of a solar radio digital receiver. Publications of the Astronomical Society of Japan, 73(2): 439-449, doi: 10.1093/pasj/psab010.

     

    Yang L J, Zhu X F, Liu S M. 2013. Measurement techniques for spurious free dynamic range (SFDR) of analog to digital converter (ADC). Journal of Beijing Institute of Graphic Communication (in Chinese), 21(6): 77-79.

     

    Ye P J, Sun Z Z, Zhang H, et al. 2019. Mission design of Chang′e-4 probe system. Scientia Sinica Technologica (in Chinese), 49(2): 123-137, doi: 10.1360/N092018-00400.

     

    Yu X D, Yuan Z G, Huang S Y, et al. 2019. Excitation of extremely low-frequency chorus emissions: The role of background plasma density. Earth and Planetary Physics, 3(1): 1-7. doi: 10.26464/epp2019001.

     

    Zarka P. 1998. Auroral radio emissions at the outer planets: Observations and theories. Journal of Geophysical Research: Planets, 103(E9): 20159-20194, doi: 10.1029/98je01323.

     

    Zhang Y Y, Zhang L, Shang Z Q, et al. 2022. A new multichannel parallel real-time FFT algorithm for a solar radio observation system based on FPGA. Publications of the Astronomical Society of the Pacific, 134(1033): 034502, doi: 10.1088/1538-3873/AC5212.

     

    刘文政, 王德恒. 2020. 前端电路设计对ADC的SFDR影响. 舰船电子对抗, 43(5): 108-111, doi: 10.16426/j.cnki.jcdzdk.2020.05.025.

     

    梅丽, 苏彦, 周建锋. 2018. 极低频射电天文观测现状与未来发展. 天文研究与技术, 15(2): 127-139, doi: 10.14005/j.cnki.issn1672-7673.20180103.012.

     

    彭璐. 2008. 星载微波辐射计的定标精度分析[硕士论文]. 北京: 中国科学院研究生院.

     

    平劲松, 王明远, 张墨等. 2021. 行星低频射电爆发的空间探测进展. 深空探测学报, 8(1): 80-91, doi: 10.15982/j.issn.2096-9287.2021.20200016.

     

    时方, 谢志远, 尹亚南等. 2015. 快速离散傅里叶变换算法研究与FPGA实现. 电测与仪表, 52(5): 15-19, doi: 10.3969/j.issn.1001-1390.2015.05.004.

     

    徐珂, 尚自乾, 严发宝等. 2021. 毫米波宽带太阳射电观测系统的信号平坦度补偿方法. 中国科学: 科学技术, 51(4): 413-423, doi: 10.1360/SST-2020-0283.

     

    杨莉军, 朱晓峰, 刘书明. 2013. 模数转换系统无杂散动态范围的测量技巧. 北京印刷学院学报, 21(6): 77-79.

     

    叶培建, 孙泽洲, 张熇等. 2019. 嫦娥四号探测器系统任务设计. 中国科学: 技术科学, 49(2): 123-137, doi: 10.1360/N092018-00400.

  • 加载中

(17)

(3)

计量
  • 文章访问数: 
  • PDF下载数: 
  • 施引文献:  0
出版历程
收稿日期:  2022-01-07
修回日期:  2022-06-21
上线日期:  2023-03-10

目录