留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

芯片式光谱仪光束耦合与对准监测系统设计

刘智颖 姜馨 李明宇 贾文波

刘智颖, 姜馨, 李明宇, 贾文波. 芯片式光谱仪光束耦合与对准监测系统设计[J]. 红外与激光工程, 2020, 49(8): 20190532. doi: 10.3788/IRLA20190532
引用本文: 刘智颖, 姜馨, 李明宇, 贾文波. 芯片式光谱仪光束耦合与对准监测系统设计[J]. 红外与激光工程, 2020, 49(8): 20190532. doi: 10.3788/IRLA20190532
Liu Zhiying, Jiang Xin, Li Mingyu, Jia Wenbo. System design for beam coupling and alignment monitoring of chip spectrometer[J]. Infrared and Laser Engineering, 2020, 49(8): 20190532. doi: 10.3788/IRLA20190532
Citation: Liu Zhiying, Jiang Xin, Li Mingyu, Jia Wenbo. System design for beam coupling and alignment monitoring of chip spectrometer[J]. Infrared and Laser Engineering, 2020, 49(8): 20190532. doi: 10.3788/IRLA20190532

芯片式光谱仪光束耦合与对准监测系统设计

doi: 10.3788/IRLA20190532
基金项目: 国家自然科学基金(61805025)
详细信息
    作者简介:

    刘智颖(1981-),女,教授,博士,主要从事先进光学设计与成像方面的研究。Email:lzycccccc@126.com

  • 中图分类号: TH744

System design for beam coupling and alignment monitoring of chip spectrometer

  • 摘要: 针对芯片式光谱仪的光束耦合与对准监测难的问题,提出了一种集成光学系统,避免光纤与芯片式光谱仪接收端接触产生磨损;有效解决光纤遮挡导致无法监测耦合效果的问题。光学系统由耦合前部分系统、监测后部分系统和复合共用系统三部分组成,复合共用系统需同时配合耦合前部分系统与监测后部分系统分别完成光束耦合与对准监测的功能。采用多重组合方式对整体系统进行设计,针对6 μm入射光纤与20 μm×20 μm的芯片式光谱仪接收端,在(1 550±50) nm的工作谱段对光纤光束与芯片式光谱仪接收端耦合系统与监测系统进行了设计,并通过LightTools软件对耦合系统进行能量分析,计算耦合效率为0.733。整体系统结构简单且无需手动调节,可同时进行光束耦合与对准监测,为芯片式光谱仪的耦合及监测提供了一种新方法。
  • 图  1  芯片式光谱仪入射光束耦合与对准监测系统结构示意图

    Figure  1.  Sketch of system for incident beam coupling and alignment monitoring of chip spectrometer

    图  2  根据Scheimpflug条件设置像面倾斜

    Figure  2.  Set the image plane tilt according to the Scheimpflug condition

    图  3  耦合系统结构图

    Figure  3.  Coupling system structure

    图  4  增大光纤与光谱仪距离示意图

    Figure  4.  Schematic diagram of increase the distance between the fiber and the spectrometer

    图  5  耦合系统能量图

    Figure  5.  Coupling system energy diagram

    图  6  监测系统结构图

    Figure  6.  Monitoring system structure

    图  7  不同波长和视场下的点列图

    Figure  7.  Spots diagram of different wavelengths and fields of views

    图  8  监测系统MTF曲线

    Figure  8.  MTF of monitoring system

    图  9  芯片式光谱仪光束耦合与对准监测集成光学系统结构图

    Figure  9.  Structure of integrated optical system for incident beam coupling and alignment monitoring of chip spectrometer

    表  1  光学设计参数

    Table  1.   Optical design parameters

    ParameterData source
    Spectral range/nm1 550±50
    Fiber NA0.12
    Fiber diameter/μm6
    Chip receiving area/μm220×20
    Angle between the optical axis and the chip normal/(°)7
    Coupling system vertical axis magnification−1×
    Monitoring system vertical axis magnification−3×
    下载: 导出CSV

    表  2  耦合系统透镜数据

    Table  2.   Coupling system lens data

    LensMaterialRadiusRadiusThicknessDiameter
    1H−K9L−15.214−5.5460.8002.536
    2H−K9L64.029−11.3880.8002.706
    Beam splitterH−K9LInfinityInfinity0.5003.748
    3H−ZK313.3856.1841.0002.326
    4H−K9L22.000−10.3321.0002.364
    5H−ZK35.40030.5451.0002.278
    下载: 导出CSV

    表  3  监测系统透镜数据

    Table  3.   Monitoring system lens data

    LensMaterialRadiusRadiusThicknessDiameter
    1H-ZK330.5455.4001.0002.278
    2H-K9L−10.33222.0001.0002.364
    3H-ZK36.18413.3851.0002.326
    4H-ZK3−7.211−20.4371.0002.336
    5H-K9L20.708−30.4531.0002.178
    6H-K9L65.43015.9771.0002.090
    下载: 导出CSV
  • [1] Wen Zhiyu, Chen Gang, Wen Zhongquan, et al. A hybrid integrated micro fiber spectrometer design, simulation and primary experiments [J]. Acta Optica Sinica, 2003(6): 740-744. (in Chinese)
    [2] Zhang Heli, Jiang Huan, Zhang Bing, et al. Design and properties approach of an intergrated interference chip made by optical waveguides [J]. Acta Photonica Sinica, 2019, 48(1): 0123002. (in Chinese) doi:  10.3788/gzxb20194801.0123002
    [3] Kern P, Coarer E L, Benech P. On-chip spectro-detection for fully integrated coherent beam combiners [J]. Optics Express, 2009, 17(3): 1976-1987. doi:  10.1364/OE.17.001976
    [4] Li Ming, Gao Chao, Chen Fang, et al. Structural design of miniature optical fiber spectrometer [J]. Optics & Optoelectronic Technology, 2013, 11(2): 100-103. (in Chinese)
    [5] Yang Minyue, Li Mingyu, He Jianjun. Polarization insensitive arrayed-input spectrometer chip based on silicon-on-insulator echelle grating [J]. Chinese Optics Letters, 2017, 25(8): 66-69.
    [6] Benedikovic D, Chenben P, Schmid J H, et al. Subwavelength index enguneered surface grating coupler with sub-decibel efficiency for 220-nm silicon-on-insulator waveguides [J]. Optics Express, 2015, 23(17): 22628-22635. doi:  10.1364/OE.23.022628
    [7] Zhang Shiqiang, Zhang Zheng, Cai Lei, et al. Laser injecting method from free space beam to single-mode fiber using single lens [J]. High Power Laser and Particle Beams, 2014, 26(3): 41-45. (in Chinese)
    [8] Deng Ke, Wang Bingzhong, Wang Xu, et al. Technology of coupling beam into single-mode fiber in free space optical communication system [J]. Electronic Science and Technology, 2007(5): 889-891. (in Chinese)
    [9] Wang Yuhua, Li Zhihua, Duan Fajie, et al. Research on coupling factor and receiving power of optical fiber [J]. Opto-Electronic Engineering, 2005(7): 42-44, 58. (in Chinese)
    [10] Christopher Richard Doerr; Katsunari Okamoto. Advances in silica planar lightwave circuits [J]. Journal of Lightwave Technology, 2006, 24(12): 4763-4789. doi:  10.1109/JLT.2006.885255
    [11] Richard Charles Alexander Pitwon, Lars Brusberg, Henning Schröder, et al. Pluggable electro-optical circuit board interconnect based on embedded graded-index planar glass waveguides [J]. Journal of Lightwave Technology, 2015, 33(4): 741-754. doi:  10.1109/JLT.2014.2361438
    [12] You Yang, Zhao Ming, Yang Zhenyu. Progress of study on waveguide and fiber coupling design [J]. Laser & Optoelectronics Progress, 2013, 50(2): 62-71. (in Chinese)
    [13] Yu Daoyin, Tan Hengying. Engineering Optics[M].Beijing:Mechanical Industry Press, 2011. (in Chinese)
    [14] Warren J S. Modern Optical Engineering[M].Beijing:Chemical Industry Press, 2017. (in Chinese)
  • [1] 赵海博, 刘彦丽, 杨雯铄, 苏云, 高大化, 孙权森, 赵慧洁.  双通道衍射计算成像光谱仪系统 . 红外与激光工程, 2022, 51(5): 20220077-1-20220077-8. doi: 10.3788/IRLA20220077
    [2] 于磊.  成像光谱仪的发展与应用(特邀) . 红外与激光工程, 2022, 51(1): 20210940-1-20210940-11. doi: 10.3788/IRLA20210940
    [3] 张玥婷, 谭毅, 王继红, 彭起, 杨智焜.  激光窗口杂散光固气耦合热效应对光束质量的影响 . 红外与激光工程, 2022, 51(9): 20210966-1-20210966-9. doi: 10.3788/IRLA20210966
    [4] 赵宇宸, 胡长虹, 吕恒毅, 孙铭, 李兆南.  紧凑型偏视场多光路耦合同轴四反光学系统设计 . 红外与激光工程, 2021, 50(3): 20200197-1-20200197-9. doi: 10.3788/IRLA20200197
    [5] 蔺淑琴, 张继前, 朱新蕾, 王飞, 蔡阳健, 余佳益.  非均匀关联光束在湍流中的单模光纤耦合效率(特邀) . 红外与激光工程, 2020, 49(12): 20201049-1-20201049-7. doi: 10.3788/IRLA20201049
    [6] 袁立银, 谢佳楠, 侯佳, 吕刚, 何志平.  紧凑型红外成像光谱仪光学设计 . 红外与激光工程, 2018, 47(4): 418001-0418001(6). doi: 10.3788/IRLA201847.0418001
    [7] 李金金, 孙晓兵, 康晴, 李树, 殷玉龙.  偏振光谱仪偏振探测精度分析 . 红外与激光工程, 2018, 47(1): 123002-0123002(7). doi: 10.3788/IRLA201847.0123002
    [8] 柯熙政, 罗静, 雷思琛.  空间光耦合自动对准方法与实现 . 红外与激光工程, 2018, 47(1): 103009-0103009(7). doi: 10.3788/IRLA201847.0103009
    [9] 柯熙政, 张雅.  高斯阵列光束自耦合特性的实验研究 . 红外与激光工程, 2017, 46(8): 822003-0822003(9). doi: 10.3788/IRLA201746.0822003
    [10] 周树道, 马忠良, 王敏, 彭舒龄.  采用光束扫描的透射仪测量光路自动准直系统 . 红外与激光工程, 2017, 46(10): 1017001-1017001(7). doi: 10.3788/IRLA201769.1017001
    [11] 谭亚运, 张合, 张祥金, 查冰婷.  水下单光束扫描探测系统光路参数优化设计 . 红外与激光工程, 2016, 45(6): 612001-0612001(6). doi: 10.3788/IRLA201645.0612001
    [12] 张宁, 刘宇龙, 吴嘉辉, 徐熙平.  微型光谱仪的CCD 数据采集系统设计 . 红外与激光工程, 2015, 44(1): 141-147.
    [13] 杨清波, 曹开锐, 赵生, 谭立英, 马晶, 于思源.  星地链路空间光至单模光纤耦合效率补偿特性研究 . 红外与激光工程, 2015, 44(9): 2789-2794.
    [14] 郝爱花, 胡炳樑.  非平行光束光栅照明C-T成像光谱仪消像散研究与设计 . 红外与激光工程, 2015, 44(7): 2132-2136.
    [15] 殷世民, 梁永波, 朱健铭, 梁晋涛, 陈真诚.  傅里叶变换红外成像光谱仪实时光谱复原FPGA芯片系统研究 . 红外与激光工程, 2015, 44(12): 3580-3586.
    [16] 刘子寒, 季轶群, 石荣宝, 陈宇恒, 沈为民.  机载红外推扫成像光谱仪光学设计 . 红外与激光工程, 2014, 43(9): 2941-2946.
    [17] 龚华平, 杨效, 屠于梦, 宋海峰, 董新永.  FBG传感器和电阻应变仪的振动监测特性 . 红外与激光工程, 2013, 42(3): 810-813.
    [18] 刘齐民, 阮萍, 李福, 潘海俊.  空间光谱仪光栅柔性支撑设计与分析 . 红外与激光工程, 2013, 42(9): 2457-2461.
    [19] 薛向尧, 高云国, 张文豹, 张磊.  库德光路光束指向稳定性分析 . 红外与激光工程, 2013, 42(6): 1514-1518.
    [20] 薛向尧, 高云国.  利用光束后向散射提高激光跟踪发射精度方法 . 红外与激光工程, 2013, 42(8): 2003-2007.
  • 加载中
图(9) / 表(3)
计量
  • 文章访问数:  589
  • HTML全文浏览量:  236
  • PDF下载量:  44
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-05
  • 修回日期:  2020-02-15
  • 网络出版日期:  2020-09-22
  • 刊出日期:  2020-08-28

芯片式光谱仪光束耦合与对准监测系统设计

doi: 10.3788/IRLA20190532
    作者简介:

    刘智颖(1981-),女,教授,博士,主要从事先进光学设计与成像方面的研究。Email:lzycccccc@126.com

基金项目:  国家自然科学基金(61805025)
  • 中图分类号: TH744

摘要: 针对芯片式光谱仪的光束耦合与对准监测难的问题,提出了一种集成光学系统,避免光纤与芯片式光谱仪接收端接触产生磨损;有效解决光纤遮挡导致无法监测耦合效果的问题。光学系统由耦合前部分系统、监测后部分系统和复合共用系统三部分组成,复合共用系统需同时配合耦合前部分系统与监测后部分系统分别完成光束耦合与对准监测的功能。采用多重组合方式对整体系统进行设计,针对6 μm入射光纤与20 μm×20 μm的芯片式光谱仪接收端,在(1 550±50) nm的工作谱段对光纤光束与芯片式光谱仪接收端耦合系统与监测系统进行了设计,并通过LightTools软件对耦合系统进行能量分析,计算耦合效率为0.733。整体系统结构简单且无需手动调节,可同时进行光束耦合与对准监测,为芯片式光谱仪的耦合及监测提供了一种新方法。

English Abstract

    • 光谱仪器是物质化学组成及含量的重要分析仪器,已在现代科学实验、生物医学、工农业生产、国防等领域得到广泛的应用。传统光谱仪器使用条件苛刻、体积庞大,因而大大地限制了其应用范围。芯片式光谱仪的出现解决了尺寸过大带来的系列问题,但光束的传输主要还是依靠光纤来进行[1-5]

      为解决芯片式光谱仪与光纤的耦合问题,国内外研究人员提出许多方案,目前应用最广的方案是通格衍射原理,将光纤输入的光束衍射到水平方向,耦合进波导中去。但上述耦合方式由于光栅耦合器与光纤接触会造成磨损;由于光纤本身的遮挡导致无法观测到其对准情况;且不同类型的芯片式光谱仪的光栅耦合器与光纤端面尺寸不匹配导致能量损耗[6-11]

      若增大光纤与芯片式光谱仪之间的距离,虽然能够解决芯片式光谱仪的磨损和遮挡问题,但是由于光纤传输的光发散角较大,而光栅耦合器尺寸在微米级别,当光传输到芯片式光谱仪时光斑远大于光栅耦合器尺寸,耦合效果不好,且对准精度要求高,每次进行光束耦合及监测时都需利用六维调节架反复调节光纤位置,成本极高[12]

      针对芯片式光谱仪的耦合及监测,文中提出了一种芯片式光谱仪光束耦合与对准监测集成光学系统。该系统有效解决了尺寸不匹配、芯片磨损、调节困难等问题,且能够对耦合效果进行实时监测。为芯片式光谱仪的耦合及监测提供了一种新方法。

    • 芯片式光谱仪入射光束耦合与对准监测系统结构如图1所示,系统由耦合前部分系统、监测后部分系统和复合共用系统三部分组成,复合共用系统需要同时配合耦合前部分系统与监测后部分系统分别完成光束耦合与耦合效果监测的功能,成功解决了监测系统的遮挡问题。即耦合前部分系统与复合共用系统构成了耦合系统,监测后部分系统与复合共用系统构成了监测系统。所设计的监测后部分系统与耦合系统为非共轴光学系统,不影响耦合系统工作,可在实现光斑与光谱仪对准后去除监测后部分系统,实现监测后部分系统共用。

      图  1  芯片式光谱仪入射光束耦合与对准监测系统结构示意图

      Figure 1.  Sketch of system for incident beam coupling and alignment monitoring of chip spectrometer

      由于芯片式光谱仪接收端为光栅结构,在接收光束的同时存在反射,且光栅耦合需满足布拉格条件:当入射光矢量与光栅矢量叠加后的矢量方向刚好和波导的某一模式的传输方向上的分量相等时,可以通过光栅耦合进入波导中,因此将耦合系统倾斜一定角度提高耦合效率,并避免光在光栅结构上反射回发射端。耦合系统将光纤发出的光束以最大效率导入到芯片上,为了降低由于对准产生的能量损耗,要求经耦合系统在芯片式光谱仪上成像尺寸至少小于芯片式光谱仪接收端的尺寸,这里取接收端的一半。监测系统将芯片接收端放大成像在探测器上,利用分光装置同时实现反射和透射。

    • 耦合系统满足拉赫不变量:

      $${{J}} = nyu = n'y'u'$$ (1)
      $$\beta = \frac{{y'}}{y} = \frac{{nl'}}{{n'l}} = \frac{{nu}}{{n'u'}}$$ (2)

      式中:n为物方介质折射率;y为物高;u为物方孔径角;n′为像方介质折射率;y′为像高;u′为像方孔径角;β为垂轴放大率。即光学系统在近轴区内成像时,在一对共轭平面内,物高y、孔径角u和介质折射率n的乘积为一常数。y′u′是两个成反比例关系的量,增大像高,必然减小像方孔径角。耦合系统要求在芯片接收范围内光斑较大以便进行对准,且像方孔径角较大以便后方监测系统提高放大倍率,综合考虑取β=−1×左右[13]

      考虑到要与监测系统进行拼接,耦合系统与监测系统共用光路。因此先将光纤发出的光束准直,便于设置分光板对平行光进行透射与反射。后接复合共用系统对光线汇聚,成像在芯片的接收端面上。如此复合共用系统可同时配合耦合前部分系统与监测后部分系统分别完成光束耦合与耦合效果监测的功能。

    • 由于复合共用系统需要同时配合耦合前部分系统与监测后部分系统分别完成光束耦合与耦合效果监测的功能,因此耦合系统的像面与监测系统的物面共面,耦合系统的像方孔径角即为监测系统的物方孔径角,但是成像在芯片接收端时要求光斑小于接收端面,监测时要求可以观察到整个接收端面,因此监测系统的物面要比耦合系统的像面大。监测系统需要将芯片接收端放大成像在探测器上以观察光束耦合情况,由于耦合系统需留出较大的后截距来更换芯片式光谱仪及调整接收端位置,即监测系统物距l较大,根据公式(2)可知,β值越大,像距l′越大,因此可以设置较低的放大倍率来降低系统的总长,现有的近红外探测器像元大小基本在25 μm×25 μm以内,因此设置β=−3×即可以使探测器分辨芯片式光谱仪的对准情况。通过复合共用系统将芯片式光谱仪的光线准直,再通过监测后部分系统成像在探测器上。

      由于芯片式光谱仪接收端为光栅耦合器,而光栅耦合需满足布拉格条件,因此需耦合系统整体倾斜一定角度,监测系统与耦合系统共用一部分光路,则监测系统物面与光轴不垂直。根据Scheimpflug条件:

      $$\tan \theta ' = \frac{{s'}}{s}\tan \theta = \beta \tan \theta $$ (3)

      式中:β为放大倍率,物面、透镜主平面和像面延长线交于一点,设置像面倾斜来校正,如图2所示[14]

      图  2  根据Scheimpflug条件设置像面倾斜

      Figure 2.  Set the image plane tilt according to the Scheimpflug condition

    • 本系统的指标参数如表1所示。

      表 1  光学设计参数

      Table 1.  Optical design parameters

      ParameterData source
      Spectral range/nm1 550±50
      Fiber NA0.12
      Fiber diameter/μm6
      Chip receiving area/μm220×20
      Angle between the optical axis and the chip normal/(°)7
      Coupling system vertical axis magnification−1×
      Monitoring system vertical axis magnification−3×
    • 表2为耦合系统透镜数据,图3为耦合系统结构图,光纤光束通过耦合系统高效的耦合到芯片式光谱仪上。左侧为输入光纤,右侧为芯片接收端,前两片透镜为耦合前部分系统,后三片透镜为复合共用系统。

      表 2  耦合系统透镜数据

      Table 2.  Coupling system lens data

      LensMaterialRadiusRadiusThicknessDiameter
      1H−K9L−15.214−5.5460.8002.536
      2H−K9L64.029−11.3880.8002.706
      Beam splitterH−K9LInfinityInfinity0.5003.748
      3H−ZK313.3856.1841.0002.326
      4H−K9L22.000−10.3321.0002.364
      5H−ZK35.40030.5451.0002.278

      图  3  耦合系统结构图

      Figure 3.  Coupling system structure

      图4(a)、(b)为耦合系统的点列图,表明耦合系统在不同波长和视场下点列斑分布均匀,单个弥散斑半径在3 μm内,全视场光斑直径大小在12 μm内,小于集成芯片式光谱仪接收端尺寸。因此集成芯片式光谱仪接收端可接收到光纤发出的全部光束,避免能量浪费。

      图  4  增大光纤与光谱仪距离示意图

      Figure 4.  Schematic diagram of increase the distance between the fiber and the spectrometer

    • 将耦合系统导入LightTools中进行能量分析,图5 为耦合系统能量图,可见能量集中在光栅耦合器中央位置且分布均匀。光源发出的总功率为0.005 665 7 W,通过LightTools模拟可得接收端总功率为0.004 153 3,则在光栅面处耦合效率为0.733。

      图  5  耦合系统能量图

      Figure 5.  Coupling system energy diagram

      通过LightTools模拟及计算,可知光束耦合系统达到满足测试需求的耦合效率。

    • 表3为耦合系统透镜数据,图6为监测系统结构图,底端为芯片接收端,左侧为探测器,下端为复合共用系统,上端为监测后部分系统。由于物面与光轴成7°倾斜角,通过公式计算得出像面倾斜角度,如图6所示。

      表 3  监测系统透镜数据

      Table 3.  Monitoring system lens data

      LensMaterialRadiusRadiusThicknessDiameter
      1H-ZK330.5455.4001.0002.278
      2H-K9L−10.33222.0001.0002.364
      3H-ZK36.18413.3851.0002.326
      4H-ZK3−7.211−20.4371.0002.336
      5H-K9L20.708−30.4531.0002.178
      6H-K9L65.43015.9771.0002.090

      图  6  监测系统结构图

      Figure 6.  Monitoring system structure

      图7图8为监测系统的点列图和MTF曲线,表明监测系统在不同波长和视场下点列斑均小于一个像元大小,MTF曲线均贴近衍射极限,成像效果良好,可通过监测系统清晰地观测到芯片式光谱仪接收端情况。

      图  7  不同波长和视场下的点列图

      Figure 7.  Spots diagram of different wavelengths and fields of views

      图  8  监测系统MTF曲线

      Figure 8.  MTF of monitoring system

    • 将设计优化的子系统进行拼接,完成最终的芯片式光谱仪入射光束耦合与对准监测系统的设计,如图9所示。系统用于芯片式光谱仪入射光的耦合时,可将光纤光束完整的成像在芯片接收端面上,提高能量的利用率;同时将光纤与芯片分离,解决了光纤与芯片的磨损问题;确定好芯片位置后添加一个卡槽,同尺寸同类型的芯片无需再次调整校准即可进行耦合;无需六维调节架来进行反复校准,提高效率的同时大大缩减了成本。且由于芯片式光谱仪入射光束耦合与对准监测系统设计了复合共用系统,因此可同时对芯片式光谱仪入射光的耦合效果进行实时监测,将耦合效果清晰的成像在探测器上。系统结构简单,易于固定,无需手动调节。

      图  9  芯片式光谱仪光束耦合与对准监测集成光学系统结构图

      Figure 9.  Structure of integrated optical system for incident beam coupling and alignment monitoring of chip spectrometer

    • 文中设计的芯片式光谱仪入射光束耦合与对准监测集成光学系统结构简单,像质良好。此芯片式光谱仪入射光束耦合与对准监测集成光学系统解决了芯片式光谱仪的耦合及监测难的问题,提高能量利用率,实现了光纤与芯片式光谱仪的分离,可在光束的耦合的同时对耦合效果进行实时监测,为同类型的芯片式光谱仪的耦合与监测提供了参考,具有一定现实意义。

参考文献 (14)

目录

    /

    返回文章
    返回