高分辨率微型星载相机光学系统的设计及应用

陈玉芳, 吴振聪, 王敏

陈玉芳, 吴振聪, 王敏. 高分辨率微型星载相机光学系统的设计及应用[J]. 应用光学, 2020, 41(2): 235-241. DOI: 10.5768/JAO202041.0201001
引用本文: 陈玉芳, 吴振聪, 王敏. 高分辨率微型星载相机光学系统的设计及应用[J]. 应用光学, 2020, 41(2): 235-241. DOI: 10.5768/JAO202041.0201001
CHEN Yufang, WU Zhencong, WANG Min. Design and application of optical system for high-resolution micro space-borne camera[J]. Journal of Applied Optics, 2020, 41(2): 235-241. DOI: 10.5768/JAO202041.0201001
Citation: CHEN Yufang, WU Zhencong, WANG Min. Design and application of optical system for high-resolution micro space-borne camera[J]. Journal of Applied Optics, 2020, 41(2): 235-241. DOI: 10.5768/JAO202041.0201001

高分辨率微型星载相机光学系统的设计及应用

基金项目: 中央引导地方科技发展专项(2017L3009);福建省科技厅项目(2018G32)
详细信息
    作者简介:

    陈玉芳(1997−),女,硕士,主要从事光学设计和光学检测方面的研究。E-mail:1364673876@qq.com

    通讯作者:

    王敏(1963−),女,教授级高级工程师,主要从事精密光学仪器方面的研究。E-mail:mwang@fjnu.edu.cn

  • 中图分类号: TN202; TH703

Design and application of optical system for high-resolution micro space-borne camera

  • 摘要: 为了克服太空环境的复杂性,满足航天工程的空间使用要求,研制一款2 500万像素宽光谱共焦成像的微型星载相机光学系统。该系统适应卫星发射和在轨道运行的恶劣环境,具有抗冲击震动、耐太空高温差强辐射,体积小,质量轻等优点。设计的系统可在450 nm~800 nm的谱段内清晰成像,焦距181 mm,入瞳口径45 mm,视场角10.4°,边缘相对照度0.81,轴上点MTF:0.57@55 lp/mm,0.33@110 lp/mm,畸变1.2%,镜头质量622 g,外形尺寸Φ58.3 mm×117 mm,抗辐照性能≥5 krad。通过温度适应性的模拟和优化,用户进行−30 ℃~+70 ℃光学镜头热真空试验,可正常工作。该系统已成功应用于天宫二号伴飞卫星相机中,获得的图像清晰稳定,为空间遥感实验观测发挥了重要的作用。
    Abstract: In order to overcome the complexity of the space environment and meet the space operating requirements of the space engineering, a micro space-borne camera optical system of 25-megapixel wide-spectrum confocal imaging was developed. This system adapted to the harsh environment of satellite launch and orbit operation, which had the advantages of impact shock resistance, high temperature and differential radiation resistance in space, small size and light weight. The designed system can be clearly imaged in the spectrum of 450 nm−800 nm, of which the focal length is 181 mm, the entrance pupil diameter is 45 mm, the field angle is 10.4º, the edge relative illumination is 0.81. The modulation transfer function (MTF) is 0.57 at 55 lp/mm and 0.33 at 110 lp/mm, the distortion is 1.2%, the quality of lens is 622 g, the overall dimension is Φ58.3 mm×117 mm and the anti-radiation performance is more than 5 krad. Through the simulation and optimization of the temperature adaptability, the user can perform the thermal vacuum test of −30℃~+70℃ optical lens normally. This system is successfully applied to the accompanying satellite camera of Tiangong-2 space lab, the obtained images are clear and stable, which plays an important role in the space remote sensing experiment observation.
  • 超连续谱光源在光纤传感、光谱测量、大气探测和生物医疗等领域得到了广泛的应用[1-6],基于光纤的超连续谱光源具有超宽带、高光谱功率和高光束质量等特点而被广泛研究。通常来讲,超连续谱光源主要包括2个部分,具有高功率的脉冲激光器和高非线性光纤。微结构光纤和硅基单模光纤均可以用来产生超连续谱[7-9],然而微结构光纤与泵浦激光器尾纤的耦合仍然是一项挑战。相比于微结构光纤,固态单模光纤与泵浦激光器的输出端能进行低损耗熔接,具有结构紧凑、系统高效的特点。此外,泵浦激光器也是产生超连续谱的关键部分。常用的泵浦激光器为脉冲激光器,锁模光纤激光器已经被证明能产生各种类型的脉冲,如常规孤子脉冲、耗散孤子脉冲和类噪声脉冲等。基于类噪声脉冲的锁模光纤激光器已经被证明在泵浦非线性光纤中产生超连续谱具有独特的优势。类噪声脉冲在频谱上表现出超宽带的特性,相比于传统脉冲,类噪声脉冲能更有效地实现非线性频率转换[10- 11]。此外,类噪声脉冲还具有高的峰值功率,有利于产生超连续谱。

    目前,基于类噪声脉冲的抽运非线性光纤产生超连续谱的实验已经被大量报道[12-16]。2012年J. C. HernandezGarciaa等[12]将获得的类噪声脉冲在腔外放大后泵浦0.75 km长的单模光纤得到的超连续谱范围是1.5 μm~1.75 μm,并且在1 640 nm~1 750 nm范围内光谱的不平坦度≤±1 dB。2013年Alexey Zaytsev等[13]利用ps级的类噪声脉冲泵浦长100 m的单模光纤,得到了1 050 nm~1 250 nm的平坦超连续谱,但其结构较为复杂。2014年Shih-Shian Lin等[14]报道了基于掺铒光纤放大器和高非线性光纤的超连续谱光源,相应的光谱范围为1.2 μm~2.1 μm,虽然光谱的范围较宽,但光谱不平坦。2015年Chen等[15]利用类噪声脉冲泵浦非线性光纤获得了500 nm~2 300 nm范围的超连续谱,在700 nm~1 500 nm范围内光谱的平坦度优于5 dB。2017年E. Hernández-Escobar等[16]报道了由功率放大的类噪声脉冲泵浦非线性光纤获得了超连续谱,超连续谱的覆盖范围为1 261 nm~2 261 nm,其光谱的平坦度为3 dB。在以上报道中,由类噪声脉冲泵浦非线性光纤获得的超连续谱的平坦度在3 dB以上,但获得平坦度高的超连续谱时结构又较复杂。高平坦度的超连续谱光源可以更好地满足光纤传感、光纤通信等领域的应用。因此,结构简单、紧凑,并且平坦度高的超连续谱光源值得研究。

    本文提出了一种由类噪声脉冲抽运的平坦超连续谱光源。在泵浦功率为450 mW时,实现了中心波长为1 600 nm的类噪声锁模脉冲输出,尖峰脉宽宽度为303 fs,3 dB光谱宽度为63.34 nm。利用掺铒光纤放大器将其谐振腔输出功率放大至338 mW,并且类噪声脉冲的光谱没有明显的变化。将功率放大后的类噪声脉冲耦合进一段长57 m的高非线性光纤,获得的超连续谱覆盖了1 530 nm~2 300 nm,超连续谱的最大输出功率为49.83 mW。其中,在1 736 nm~2 134 nm光谱范围内,光谱的平坦度优于0.5 dB。

    基于类噪声脉冲抽运全光纤超连续谱产生的实验结构如图1所示。由类噪声脉冲种子源、掺铒光纤放大器和一段高非线性光纤(HNLF)组成,其中类噪声脉冲种子源是基于非线性偏振旋转原理实现锁模的。由2个偏振控制器(PC1和PC2)和偏振相关隔离器(PD-ISO)组成重要的锁模器件,其中,PD-ISO还能保证光脉冲在环形腔内单向传输。工作波长为980 nm的半导体激光器(LD)通过波分复用器(WDM)泵浦一段长0.8 m的高掺杂掺铒光纤(EDF),耦合比为10:90的光纤耦合器的90%端用于腔内反馈,10%端用来输出脉冲。环形腔内加入长为1 m的色散补偿光纤(DCF)来控制腔内色散,使腔内净色散接近零色散;同时加入6 m色散位移光纤(DSF),用来增大腔长和增强非线性。腔内其他无源器件的尾纤均为SMF-28,尾纤长度估计为6.24 m,总腔长为14.04 m。在1.5 μm波段,SMF-28、DSF、DCF和EDF的色散值分别为−0.022 9 ps2/m、−0.005 1 ps2/m、+0.184 8 ps2/m和−0.02 ps2/m,因此,经过计算的环形腔的净色散为−0.0047 ps2/m,处于近零反常色散区,有利于孤子自频移等非线性效应的发生。

    图  1  基于类噪声脉冲抽运的全光纤超连续谱实验结构
    Figure  1.  Experimental structure diagram of all-fiber super-continuum spectrum based on noise-like pulse pumping

    由于超连续谱的产生存在阈值,因此,为了避免种子脉冲的输出功率达不到阈值,在激光器输出端搭建了掺铒光纤放大器。掺铒光纤放大器包括LD(Pump2)、2 m掺铒光纤(EDF2)、WDM2以及隔离器(ISO),ISO的隔离损耗为33 dB,采用背向泵浦的方式来提高泵浦效率,在放大器的输出端熔接一段HNLF来产生超连续谱。实验过程中采用最高分辨率为0.05 nm的光谱分析仪(YOKOGAWA AQ6375)测量光谱,时域信号由10 Gb/s光电探测器和2.5 GS/s示波器(OSC Agilent DSO9254A)来监测,采用自相关仪(femtochrome,FR-103XL)观测锁模脉冲的自相关轨迹。

    在实验中,通过调节抽运泵浦功率和偏振控制器,在泵浦功率为450 mW时出现了稳定的类噪声脉冲锁模,此时的输出功率约为3.13 mW,类噪声脉冲的光谱如图2(a)所示。从图2(a)可看出,光谱较平滑,和报道的类噪声脉冲光谱相似。类噪声脉冲的中心波长为1 600 nm,3 dB光谱宽度为63.34 nm。示波器观测的脉冲序列如图2(b)所示,脉冲的幅度基本均匀一致,脉冲间隔为71.22 ns,对应的重复频率为14.24 MHz,与14.04 m的腔长相对应。为了进一步确定是类噪声脉冲锁模,测量的自相关轨迹如图2(c)所示。图2(c)中插图为大范围的自相关轨迹,一个窄的尖峰位于一个大的基底上,这是典型的类噪声脉冲自相关轨迹,由一系列随机相位、强度和脉冲宽度的超短脉冲组成的脉冲包络。对尖峰采用高斯曲线拟合,相应的尖峰脉宽宽度为303 fs。类噪声脉冲的频谱图如图2(d)所示,信噪比为46 dB,其中插图为400 MHz范围的频谱图,频谱没有调制,表明类噪声脉冲锁模稳定。

    图  2  类噪声脉冲锁模
    Figure  2.  Noise-like pulse mode locking

    在获得类噪声脉冲锁模后,将泵浦功率增加至1 000 mW,类噪声脉冲仍能保持稳定的单脉冲运行。图3给出了不同泵浦功率下类噪声脉冲的输出光谱。从图3(a)可看出泵浦功率从450 mW增加至1 000 mW过程中,光谱的形状基本保持不变,而光谱的强度随着泵浦功率的增加有轻微的增加。此外,不同泵浦功率下,类噪声脉冲的自相关轨迹如图3(c)所示,与图2(c)中插图的形状非常相似,在基底上有一个尖峰,尖峰的宽度保持在303 fs没有改变。图3(b)显示了类噪声脉冲光谱的3 dB宽度随泵浦功率的变化。从图3(b)中可以看出,光谱宽度与泵浦功率之间呈线性关系。在泵浦功率达到最大1 000 mW时,类噪声脉冲的3 dB光谱宽度能达到69.66 nm。

    图  3  类噪声脉冲输出特性
    Figure  3.  Noise-like pulse output characteristics

    由于谐振腔的最大输出功率为8.6 mW,无法直接泵浦高非线性光纤,需要进行功率放大,为此设计了掺铒光纤放大器结构。为了确定类噪声脉冲经过掺铒放大器时光谱没有较大的变化,给出了种子脉冲经放大器放大后的光谱对比,如图4所示。图4中黑色虚线为没有经过放大的类噪声脉冲光谱,红色实线为经过放大后测得的类噪声脉冲光谱。从图4中可以看出,光谱强度增加,但是光谱没有较大的改变,说明光脉冲在放大过程中没有明显的非线性效应。图5给出了放大器的输出功率随泵浦功率的变化。从图5中可以看出,随着泵浦功率的增加,掺铒放大器的斜率有缓慢下降的趋势。当放大器的泵浦功率为1 000 mW时,放大器的最大输出功率为338 mW。由于泵浦2的最大功率为1 000 mW,因此,不能获得更高的输出功率。

    图  4  类噪声脉冲放大前后的光谱对比
    Figure  4.  Comparison of spectrum before and after noise-like pulse amplification
    图  5  放大器输出功率随泵浦功率变化
    Figure  5.  Amplifier output power varies with pump power

    将长为57 m的HNLF与放大器的输出端熔接可产生超连续谱。使用的HNLF在1 550 nm处非线性系数估计为10 W−1 km−1,截止波长为1 480 nm,此外HNLF的零色散波长也是1 550 nm。由于类噪声脉冲包络中飞秒脉冲的峰值功率较高,在高非线性光纤中受非线性效应影响使光谱展宽,如图6(a)所示。从图6(a)可看出,在泵浦2的功率为300 mW时,受到四波混频(FWM)和光孤子效应的作用,在HNLF的零色散波长附近光谱会急剧展宽;扩展到反常色散区的光谱在自相位调制和反常色散的作用下形成了高阶孤子。进一步提高泵浦功率,脉冲的峰值功率会超过受激拉曼散射的阈值,在受激拉曼效应的作用下会发生孤子自频移,导致光谱向更长的波长范围扩展。随着泵浦功率的进一步提高,光谱宽度会更宽[17]。从图6(a)可知,超连续谱覆盖范围为1 530 nm~2 300 nm。随着泵浦功率的增加,超连续谱的光谱强度有所增加,但是整体范围变化不明显,光谱截止在2 300 nm,这是由于光纤的非线性效应和石英玻璃光纤对2 300 nm以上的中红外波段具有很强的吸收损耗。光谱的10 dB宽度范围从1 547 nm~2 182 nm,宽度为634 nm。值得注意的是,在1 736 nm~2 084 nm范围内,光谱的不平坦度≤±0.5 dB,在1 562 nm处尖峰是由于泵浦光没有被充分吸收所致。此外,在1 400 nm和1 900 nm处光谱存在缺陷,这是HNLF中水吸收引起的。图6(b)给出了超连续谱的输出功率和10 dB带宽随放大器功率的变化,超连续谱的输出功率呈线性增加,其最大输出功率为49.83 mW,超连续谱的10 dB也随泵浦功率的增加而增加,最大带宽为639 nm。

    图  6  超连续谱和输出功率以及带宽
    Figure  6.  Output power and bandwidth of super-continuum spectrum

    设计了一个由类噪声脉冲抽运的全光纤结构的超连续谱光源。利用DCF和DSF控制腔内色散,提高非线性效应,在近零负色散区通过调节PCs和泵浦功率实现了类噪声脉冲锁模。类噪声脉冲的中心波长为1 600 nm,3 dB光谱宽度为69.66 nm,重复频率为14.04 MHz,脉冲尖峰宽度为303 fs。然后将此类噪声脉冲的直接输出功率放大,光谱没有明显的变化。将功率放大后的类噪声脉冲注入到一段长57 m的HNLF光纤中,实现了宽带的超连续谱,20 dB光谱范围覆盖了1 530 nm~2 300 nm,超连续谱的10 dB宽度为639 nm。此外,在1 736 nm~2 134 nm范围内光谱的平坦度优于0.5 dB。设计的结构简单、平坦的超连续谱光源在光传感、光通信和光谱学等领域具有较大的潜在应用价值。

  • 图  1   星载相机光学系统原理

    Figure  1.   Schematic diagram of space-borne camera optical system

    图  2   传递函数

    Figure  2.   MTF curves

    图  3   星点图

    Figure  3.   Spot diagrams

    图  4   球差和位置色差

    Figure  4.   Spherical aberration and position chromatism

    图  5   场曲和畸变

    Figure  5.   Field curvature and optical distortion

    图  6   倍率色差

    Figure  6.   Lateral color

    图  7   相对亮度

    Figure  7.   Relative illumination

    图  8   环围能量比

    Figure  8.   Diffraction encircle energy ratio

    图  9   离焦MTF曲线

    Figure  9.   Defocus MTF curve

    图  10   1 000 M近摄处MTF曲线

    Figure  10.   MTF curve at 1 000 m

    图  11   温度变化产生的像面漂移

    Figure  11.   Image plane drift caused by temperature changes

    图  12   真空下温度变化产生的像面漂移

    Figure  12.   Image plane drift caused by temperature changes in vacuum

    图  13   p=0,T=3 ℃处MTF

    Figure  13.   p =0, T =3 ℃ MTF

    图  14   p=0,T=18 ℃处MTF

    Figure  14.   p=0,T=18 ℃ MTF

    图  15   p=0,T=28 ℃处MTF

    Figure  15.   p=0,T=28 ℃ MTF

    图  16   实测MTF

    Figure  16.   Measured MTF

    图  17   2 500万像素星载相机光学系统

    Figure  17.   Space-borne camera optical system of 25-megapixel

    图  18   伴星拍摄的组合体可见光图像

    Figure  18.   Visible image by space-borne camera

    表  1   系统性能技术指标要求

    Table  1   System performance specification requirements

    ParametersValue
    Focal length/mm 179
    Pupil diameter/mm ≥44
    Distortion/% <3
    Lens weight/g ≤800
    Every lens weight/g ≤300
    Imaging band/nm 450~800
    Radiation resistance/krad ≥5
    MTF Point MTF on axis:≥0.5@55 lp/mm
    ≥0.3@110 lp/mm
    0.6 W field of MTF:≥0.4@55 lp/mm
    ≥0.25@110 lp/mm
    Boundary dimension/mm Φ70 ×150
    Field angle/(º) 10.4(8.35º×6.25º)
    Relative illumination ≥0.7
    下载: 导出CSV

    表  2   光学系统结构参数

    Table  2   Structural parameters of optical system

    SurfRadius/mmThickness/mmGlassSurfRadius/mmThickness/mmGlass
    1 5 SLLICA 11 48.24 6.26 H-ZF62
    2 0.1 12 −242.53 13.89 H-ZLAF76
    3 49.38 10.5 H-FK61 13 28.59 9.5
    4 −322.26 0.18 STO 12.95
    5 54.39 5.82 H-FK61 15 −98.86 4.2 F5
    6 185.58 1.09 16 −23.13 5.85 H-ZLAF68
    7 0.38 17 −56.19 63.94
    8 −637 2.6 H-ZLAF76 18 3 K9
    9 35.57 5.56 H-LAK2 19 0
    10 121.41 0.15 IMA 0
    下载: 导出CSV

    表  3   光学系统玻璃材料参数

    Table  3   Glass material parameters of optical system

    Materialndvd
    SLLICA 1.458 464 67.821 4
    H-FK61 1.497 81.61
    H-ZLAF76 1.850 13 30.06
    H-LAK2 1.692 11 54.54
    H-ZF62 1.922 86 20.88
    F5 1.624 35 35.92
    H-ZLAF68 1.883 39.22
    K9 1.516 8 64.2
    下载: 导出CSV
  • [1]

    HAN Changyuan. Recent earth imaging commercial satellites with high resolutions[J]. Chin. J. Opt. Appl. Opt.,2010,3(3):202-208.

    [2] 焦斌亮. 星载小型CCD相机发展现状及趋势[J]. 空间电子技术,1995(4):34-40.

    JIAO Binliang. Development status and trend of small spaceborne CCD cameras[J]. Spatial Electronics Technology,1995(4):34-40.

    [3] 兰丽艳, 黄颖. 星载大视场多光谱高分辨率CCD相机光学系统设计[J]. 航天返回与遥感,2002,23(4):34-37. doi: 10.3969/j.issn.1009-8518.2002.04.007

    LAN Liyan, HUANG Ying. Optical system design of wide-field, multi - spectrum and high-resolution CCD camera[J]. Journal of Astronautic Metrology,2002,23(4):34-37. doi: 10.3969/j.issn.1009-8518.2002.04.007

    [4] 张丹枫. 星载多光谱相机光学系统设计[D]. 长春: 长春理工大学, 2012: 6-9.

    ZHANG Danfeng, Design of optical system for spaceborne multi-spectral camera[D]. Changchun: Changchun University of Science and Technology, 2012: 6-9.

    [5]

    CHEN Y H, HSIEH Y F, HSU C W, et al. An image based optical lens eccentric error inspection system[C]// Instrumentation and Measurement Technology Conference (I2MTC), 2012 IEEE International. USA: IEEE, 2012.

    [6] 李晓彤, 岑兆丰. 几何光学•像差•光学设计[M]. 杭州: 浙江大学出版社, 2003: 159-166.

    LI Xiaotong, CEN Zhaofeng. Geometrical optics, aberrations and optical design[M]. Hangzhou: Zhejiang University Press, 2003: 159-166.

    [7] 刘乾, 杨维川, 袁道成, 等. 光谱共焦显微镜中色散物镜材料的优化选择[J]. 光电工程,2012,39(8):111-116. doi: 10.3969/j.issn.1003-501X.2012.08.017

    LIU Gan, YANG Weichuan, YUAN Daocheng, et al. Optimization and selection of materials for dispersive objective of chromatic confocal microscope[J]. Optoelectronic Engineering,2012,39(8):111-116. doi: 10.3969/j.issn.1003-501X.2012.08.017

    [8] 李博. 靶场紫外望远系统光学设计[J]. 中国激光,2014,41(10):267-272.

    LI Bo. Optical design of ultraviolet telescope system in shooting range[J]. Chinese Journal of Lasers,2014,41(10):267-272.

    [9] 巩盾, 王红, 田铁印. 温度对星载相机的影响和温控指标的制定[J]. 光学学报,2010,30(7):2017-2021. doi: 10.3788/AOS20103007.2017

    GONG Dun, WANG Hong, TIAN Tieyin. Temperature effect on imaging performance of three-line CCD satellite camera and temperature indicators affirmed[J]. Acta Optica Sinica,2010,30(7):2017-2021. doi: 10.3788/AOS20103007.2017

    [10] 李利, 张凯迪. 高分辨率超低畸变航天光学成像系统设计[J]. 应用光学,2019,40(3):363-368. doi: 10.5768/JAO201940.0301001

    LI Li, ZHANG Kaidi. Design of optical imaging system with high resolution, long-focus, wide angle and low distortion[J]. Journal of Applied Optics,2019,40(3):363-368. doi: 10.5768/JAO201940.0301001

    [11] 张建隆, 贺磊, 杨振. 高分辨率长焦广角低畸变光学成像系统设计[J]. 应用光学,2017,38(5):725-731.

    ZHANG Jianlong, HE Lei, YANG Zhen. Design of optical imaging system with high resolution, long-focus, wide angle and low distortion[J]. Journal of Applied Optics,2017,38(5):725-731.

    [12] 李博, 王淑荣, 王俊博. 天基紫外跟踪详查光学系统设计[J]. 中国激光,2013,40(9):252-256.

    LI Bo, WANG Shurong, WANG Junbo. Study on optical system design of space-based ultraviolet tracking[J]. Chinese Journal of Lasers,2013,40(9):252-256.

    [13] 温广武, 雷廷权, 周玉. 石英玻璃基复合材料的研究进展[J]. 材料工程,2002(1):40-43. doi: 10.3969/j.issn.1001-4381.2002.01.012

    WEN Guangwu, LEI Tingquan, ZHOU Yu. Progress in research on fused silica matrix composites[J]. Material Engineering,2002(1):40-43. doi: 10.3969/j.issn.1001-4381.2002.01.012

    [14] 张树青, 李韬. 含扫描反射镜的星载相机异速像移分析[J]. 光学学报,2014,34(2):292-297.

    ZHANG Shuqing, LI Tao. Analysis of all-speed image shift of on-board camera with scanning mirror[J]. Acta Optica Sinica,2014,34(2):292-297.

    [15] 焦斌亮, 王朝晖, 林可祥, 等. 星载多光谱CCD相机研究[J]. 仪器仪表学报,2004,25(2):146-148. doi: 10.3321/j.issn:0254-3087.2004.02.002

    JIAO Binliang, WANG Zhaohui, LIN Kexiang, et al. Study of spaceborne multispectral CCD camera[J]. Chinese Journal of Scientific Instrument,2004,25(2):146-148. doi: 10.3321/j.issn:0254-3087.2004.02.002

  • 期刊类型引用(0)

    其他类型引用(2)

图(18)  /  表(3)
计量
  • 文章访问数:  1186
  • HTML全文浏览量:  352
  • PDF下载量:  161
  • 被引次数: 2
出版历程
  • 收稿日期:  2019-05-19
  • 修回日期:  2019-06-23
  • 网络出版日期:  2020-03-31
  • 刊出日期:  2020-02-29

目录

/

返回文章
返回