-
低噪声单频激光器是激光干涉仪空间引力波探测系统中的核心器件, 其噪声性能直接影响引力波探测器的灵敏度[1]. 由于空间引力波探测系统的探测目标是0.1 mHz—1 Hz频段内的引力波, 所以这对激光器在毫赫兹频段的噪声性能提出严苛的要求[2]. 例如, 在发展最久、成熟度最高的激光干涉仪空间天线(laser interferometer space antenna, LISA)项目中, 要求激光器在2 W输出时的相对强度噪声 (relative intensity noise, RIN) 不大于2×10–4 /Hz1/2@0.1 mHz—1 Hz [3]. 主控振荡器-功率放大器(master oscillator power amplifier, MOPA)结合低噪声种子源和高性能放大器的优点, 具有光学性能优越、易于实现等特点[4,5], 因此可用于空间引力波探测系统. 一般而言, 高性能种子源可以在大于弛豫振荡峰的频率范围获得低强度噪声, 但是由于泵浦噪声和环境因素等影响[6,7], 在低于弛豫振荡峰的频段, 尤其是在mHz—Hz频段内, 其噪声性能迅速劣化, 难以满足引力波探测需求. 如何有效地抑制毫赫兹频段的强度噪声是研制面向空间引力波探测激光器的关键问题之一.
以LISA项目为牵引, 国外研究人员对低噪声单频激光器开展了诸多研究并已取得了重要的进展. 2005年, 德国汉诺威激光中心研制了Nd:YAG非平面环形振荡器(non-planar ring oscillator, NPRO)注入到光纤放大器的MOPA激光器, 输出功率为1 W, RIN为3×10–3 /Hz1/2@1 mHz—1 Hz[8]. 2014年, 欧洲航天局(European space agency, ESA)为下一代引力任务(next generation gravity mission, NGGM)研制了高稳定性激光器, 种子激光器采用经过空间验证的NPRO, 注入功率为25 mW, 放大器输出功率>500 mW. 在该激光器中, 采用光电反馈控制技术稳定激光功率, 对泵浦激光二极管(laser diode, LD)和光纤器件进行温度控制, 获得了满足NGGM要求的低噪声激光[9]. 2018年, Dahl等[10]为LISA研制了平面波导外腔激光器(planar-waveguide external cavity laser, PW-ECL)注入的MOPA激光器, 种子激光经半导体光放大器(semiconductor optical amplifier, SOA)预放大后, 注入光纤放大器的激光功率约300 mW, 放大后得到2.1 W的输出功率, 但是低于3 mHz频段内的RIN不满足LISA的要求. 2019年, Numata等[1]研究了微型NPRO(m-NPRO)的性能, 结果表明, m-NPRO比NPRO更具优势, 采用m-NPRO注入光纤放大器的MOPA激光器是LISA中低噪声激光器的最优方案. 近年来, 国内也部署了空间引力波探测任务. 2015年, 中山大学启动了“天琴计划”[11,12]; 2016年, 中国科学院启动了“太极计划”[13,14], 其探测频段覆盖LISA的1 mHz—1 Hz. 空间引力波探测的相关研究迅速成为国内热点内容, 但是, 国内关于毫赫兹频段低噪声单频激光器的研究和报道较少.
本文采用低噪声、窄线宽NPRO激光器作为主控振荡器, 结合低噪声全保偏光纤放大技术, 利用光电负反馈回路控制放大器的泵浦电流, 从而对输出激光的强度噪声进行抑制, 获得了毫赫兹频段的低噪声激光输出, 在输出激光功率为2.13 W时测得1 mHz—1 Hz频段内的RIN均低于–60 dBc/Hz, 在1 mHz和1 Hz处的RIN分别为–63.4 dBc/Hz和–105.8 dBc/Hz.
-
图1为低噪声单频激光器的示意图, 激光器采用MOPA结构, 包括主控振荡器、全保偏光纤放大器和反馈控制回路三部分. 主控振荡器采用窄线宽单频NPRO激光器, 输出功率最大达500 mW, 线宽约为1 kHz; 全保偏光纤放大器由泵浦LD, (2+1)×1 保偏光纤合束器、2.1 m长的双包层掺镱增益光纤(PLMA-YDF-10/125-VIII, Nufern)、包层功率剥除器(cladding power stripper, CPS)以及光纤分束器等组成; 反馈控制回路由光纤衰减器、光电探测器以及比例-积分-微分(proportional- integral- derivative, PID)控制器组成. 主控振荡器输出的信号光经保偏光纤隔离器后与泵浦光合束, 共同注入到保偏增益光纤中进行功率放大, 输出激光经过CPS剥除包层中的残余泵浦光功率, 保证只有1064 nm的信号光输出. 保偏光纤分束器(分光比为1/99) 1%端的光束进入反馈控制回路, 经光电探测器转换为电信号, 与PID控制器的内部基准电压信号进行比较, 得到误差信号, 利用该误差信号实时调整泵浦LD的驱动电流, 通过调整泵浦功率实时改变放大器的输出激光功率, 实现输出激光强度噪声的抑制.
对于MOPA结构的激光器, 光纤放大器中泵浦LD的功率波动是造成输出激光RIN在低频段劣化的主要因素之一, 因而降低泵浦LD的功率波动是抑制激光RIN的重要途径[15,16]. 由于掺镱增益光纤在976 nm附近的吸收系数对泵浦波长敏感, 因此, 激光器的强度噪声与泵浦LD的波长稳定性有关. 通过实验比较了不同温度下的两种泵浦LD的光谱, 结果如图2所示, 图2(a)为不锁波长泵浦LD的光谱, 随着工作温度的升高, 其波长向长波方向漂移; 图2(b)为锁波长泵浦LD的光谱, 随温度变化漂移较小, 波长稳定性高. 因此, 为了降低泵浦源的噪声, 选用波长锁定的泵浦LD.
-
激光器输出功率的波动情况在时域上用功率稳定性描述, 在频域上用强度噪声描述. 在分析激光器的强度噪声性能时, 通常采用RIN来衡量, 定义为
$ {{\rm{RIN}}}= {{S}_{{\rm{p}}}(f)}/{{\overline P}^{2}} , $ 式中,
${S}_{{\rm{p}}}(f)$ 为输出激光功率波动的功率谱密度;$\overline{P}$ 为激光器输出的平均功率[17].在评估RIN时, 通常采用低噪声光电探测器和频谱分析仪直接测量[18]. 把输出激光衰减到合适的功率水平, 经光电探测器转化为电信号, 通过频谱分析仪测量电信号的功率谱密度. 受限于现有测量仪器的频率分辨率、测量范围和计算能力, 在测量低频段的RIN时, 测量时间长, 数据量大, 且需要分频段测量后进行拼接, 操作繁琐. 因此, 对于毫赫兹及更低频段的RIN评估, 一般采用数字万用表(digital multimeter, DMM)或高精度数字采集卡长时间记录电信号, 通过计算获得所测电信号的功率谱密度[19]. 为此, 设计和研制了高精度的4通道RIN测量系统, 采用24 bit的数据采集卡, 采样率最高为204.8 kS/s, 测量频率范围为2 μHz—102.4 kHz, 满足毫赫兹频段RIN的测量需求.
图3所示为激光RIN测量示意图, 分为内环(in-loop)和外环(out-of-loop)两路, 其中, 内环电信号也作为激光器强度噪声抑制的输入信号. 在激光器的内环控制回路中, 光信号和电信号的转化较为复杂, 经过“信号光→电信号→泵浦光→信号光”后, 各环节中的光学器件和电学器件均会引入噪声, 相当于在内环链路中引入多个技术噪声源, 其噪声经过光纤放大器后叠加到外环的光信号上, 因此, 采用光电负反馈控制的激光器的内、外环强度噪声存在较大的差别, 实测的外环噪声性能要比内环差[20].
-
测试采用铟镓砷光电探测器(DET01 CFC/M, Thorlabs)和动态信号采集卡(PXI4462, NI)进行长时间的记录, 并利用LPSD(logarithmic frequency axis power spectral density)算法得到对数坐标上平滑的功率谱密度[21]. 在室温条件((23±0.5) ℃)下, 当激光器输出功率为2.13 W时, 分别测得种子激光、放大器自由运转和反馈抑制时的信号, 结果如图4所示.
图4中曲线(a)为RIN测量系统的噪底, 在1 mHz—1 Hz内, 噪底低于–90 dBc/Hz; 曲线(b)和曲线(c)分别为种子激光和放大器自由运转时的RIN. 在低于50 mHz的频段内, 输出激光的RIN开始劣化, 10 mHz处劣化达10 dB, 在低于1 mHz频段劣化更为严重, 最大可达30 dB以上. 当反馈回路工作时, 测得内环和外环的RIN如曲线(d)和曲线(e)所示, 在低于0.1 Hz频段, 内、外环的RIN并不相同, 且内环明显优于外环, 随着频率的降低, 内外环的差别增大. 与自由运转时的RIN相比, 外环的RIN得到了有效的抑制, 在10 mHz处RIN抑制了约10 dB, 在1 mHz处RIN抑制了约20 dB. 但是, 在低于4 mHz频段, 外环RIN依然大于–60 dBc/Hz.
-
利用红外热成像仪(T540, FLIR)对光纤放大器内部进行拍照测温, 在放大器工作2 h后观察内部热分布情况, 结果如图5所示. 其中, 测温点①为泵浦LD, 测得LD的表面温度(26.2 ℃)高于设定的TEC控制温度(25 ℃). 测温点②为输出光纤分束器, 其温度随着输出激光功率的增大而升高, 在输出功率为2.13 W时, 测得温度为27.5 ℃. 测温点③为光纤放大器内部壳体, 温度为23 ℃, 与水冷板的水温设置值相同.
由于泵浦LD和输出光纤分束器的温度与放大器内部的温度差别较大, 存在热交换, 影响泵浦LD的泵浦功率和光纤分束器分光比的稳定性, 从而增大输出激光的RIN. 为降低温度对低频段RIN的影响, 对泵浦LD和光纤分束器设计了包裹更为严密、散热性能更好的紫铜夹具, 并对泵浦LD和光纤分束器进行主动精确控温. 在不同的工作温度下, 分别测量外环的RIN, 得到图6所示的RIN曲线. 从图6可知, 外环RIN受泵浦LD、光纤分束器和放大器内部温度的影响较大, 尤其是低频段的RIN更为严重. 因此, 精确控制关键器件温度和环境温度对进一步提高毫赫兹频段的强度噪声性能很关键. 通过实验测试和分析, 设置泵浦LD和光纤分束器的温度为27 ℃, 与长时间工作后放大器的内部温度十分接近, 在这种情况下, 测得外环RIN最低, 1 Hz处达–105.8 dBc/Hz(5×10–6/Hz1/2), 1 mHz处达–63.4 dBc/Hz(6.8×10–4/Hz1/2).
-
图7(a)是光纤放大器在输出功率为2.13 W时的功率稳定性. 当泵浦LD的注入电流为4.2 A, 泵浦光功率为6 W时, 光纤放大器的输出功率为2.13 W, 输出激光在4 h内的功率波动< 0.13%. 图7(b)是种子激光和保偏光纤放大器输出激光的光谱图, 可以看出, 种子激光和放大器输出激光的中心波长都是1064.35 nm, 信噪比均大于70 dB. 图7(c)是利用马赫-曾德尔干涉仪(Mach-Zehnder Interferometer, MZI)测量的2.13 W 输出功率下光纤放大器的附加线宽[22,23]. 由图7(c)可知, 光纤放大器产生的附加线宽在20 dB 处约为1.9 Hz, 因此光纤放大器对种子激光线宽的展宽非常小, 可忽略不计.
-
本文报道了一种面向空间引力波探测用的低强度噪声单频激光器. 激光器采用MOPA结构, 其中主控振荡器采用低噪声、窄线宽的NPRO激光器, 放大器采用抗外界扰动能力强的全保偏光纤放大器. 通过光电负反馈和对泵浦LD与输出光纤分束器的温度控制, 实现了低强度噪声的稳定运转. 通过自主搭建的基于高分辨率数据采集卡的强度噪声测量系统, 测得激光器RIN在1 mHz—1 Hz频段内低于–60 dBc/Hz, 在1 mHz和1 Hz处分别为–63.4 dBc/Hz和–105.8 dBc/Hz. 研究结果表明, 通过光电负反馈控制光纤放大器的泵浦电流可以有效地抑制输出激光的RIN. 另外, 在低于毫赫兹频段, 温度漂移是RIN劣化的主要原因. 因此, 对泵浦LD、光纤分束器等关键器件和环境温度进行更为精确的温度控制, 可以进一步提高激光器低频段的RIN性能.
-
低噪声单频激光器是空间引力波探测系统中的核心器件, 其噪声性能直接影响空间引力波探测器的灵敏度. 本文报道了一种面向空间引力波探测的低噪声单频激光器, 利用全保偏光纤结构的功率放大器对低功率、窄线宽、低噪声的非平面环形振荡器输出激光进行放大. 为降低激光的强度噪声, 比较了不同泵浦源的输出特性, 为光纤放大器选用波长锁定的泵浦源, 降低泵浦光波长随温度漂移对输出功率的影响, 利用光电负反馈控制技术抑制输出激光的强度噪声, 结合主动精确控温技术抑制关键器件的热噪声, 实现了毫赫兹频段强度噪声的抑制. 利用自主搭建的4通道相对强度噪声测量系统, 测得反馈控制后的激光器相对强度噪声在1 mHz—1 Hz频段内低于–60 dBc/Hz, 在1 mHz和1 Hz处分别为–63.4 dBc/Hz和–105.8 dBc/Hz. 研究结果表明, 通过放大器泵浦电流的反馈控制和关键器件的温度控制可以有效地抑制激光器在毫赫兹频段的强度噪声, 为进一步提高低频段强度噪声性能奠定基础.A low-noise single-frequency laser is a key component of the space-based gravitational wave detector, and the intensity noise of the laser directly affects the sensitivity of the space-based GW detector. In this work, we report a low-noise single-frequency laser designed for space-based gravitational wave detector. The laser is based on a master oscillator power amplifier (MOPA), which is designed to possess a low-power, narrow-linewidth seed laser acting as master oscillator (MO) and an all polarization-maintaining fiber amplifier acting as power amplifier (PA). The amplifier that uses a robust mechanical design consists of an Yb-doped double-clad fiber forward pumped by wavelength-locked 976 nm pump laser diode (LD) to achieve 2.13 W of output power and 70 dB of signal-noise ratio (SNR). To suppress the relative intensity noise (RIN) in a millhertz regime (1 mHz–1 Hz), we characterize the power stabilization of a pump diode laser based on a proportional-integral-derivative (PID) feedback control loop where an in-loop photodetector is used. The power fluctuation can be converted into the fluctuation of the current signal by the photodiode, the current signal is converted into the voltage signal and amplified by a transimpedance circuit. Then, the voltage signal is compared with the voltage reference signal, and the error signal is achieved to adjust real-timely the drive current of the pump laser diode. This is a good way to significantly suppress the RIN of a laser at low frequencies, but the measured RIN below 4 mHz is still higher than –60 dBc/Hz. In order to further suppress the RIN to lower than 4 mHz, an active precise temperature control technology is used to suppress the thermal noise from pump LD and fiber coupler. To assess the RIN milliertz regime, we design an RIN measurement system consisting of a high-precision signal acquisition card (24 bit) and a computer program based on LabVIEW. The measurement range of the system is 2 μHz–102.4 kHz and the frequency resolution up to 2 µHz, much better than the counterparts of commercial instruments. By stabilizing the fiber amplifier pump LD current and the temperature of pump LD and the temperature of fiber coupler, the out-of-loop RINs are measured to be –63.4 dBc/Hz@1 mHz and –105.8 dBc/Hz@1 Hz , and in a milliertz regime of 1 mHz–1 Hz, the RIN is below –60 dBc/Hz. The results show that the feedback control of the fiber amplifier pump LD current and the temperature control of the key devices can effectively suppress the RIN in the millihertz frequency band, which lays a foundation for further improving the intensity noise performance in the low frequency band.
-
Keywords:
- single-frequency laser /
- millihertz band /
- relative intensity noise /
- optoelectronic negative feedback
[1] Numata K, Yu A, Jiao H, Merritt S, Micalizzi F, Fahey M, Camp J, Krainak M 2019 Proceedings of SPIE: Solid State Lasers XXIII: Technology and Devices 10896 108961H
Google Scholar
[2] 罗子人, 白姗, 边星, 陈葛瑞, 董鹏, 董玉辉, 高伟, 龚雪飞, 贺建武, 李洪银, 李向前, 李玉琼, 刘河山, 邵明学, 宋同消, 孙保三, 唐文林, 徐鹏, 徐生年, 杨然, 靳刚 2013 力学进展 43 415
Google Scholar
Luo Z R, Bai S, Bian X, Chen G R, Dong P, Dong Y H, Gao W, Gong X F, He J W, Li H Y, Li X Q, Li Y Q, Liu H S, Shao M X, Song T X, Sun B S, Tang W L, Xu P, Xu S N, Yang R, Jin G 2013 Adv. Mech. 43 415
Google Scholar
[3] Peterseim M, Brozek O S, Danzmann K, Freitag I, Rottengatter P, Tünnermann A, Welling H 1998 AIP Conf. Proc. 456 148
Google Scholar
[4] O’Brien S, Welch D, Parke R, Mehuys D, Dzurko K, Lang R, Waarts R, Scifres D 1993 IEEE J. Quantum Electron. 29 2052
Google Scholar
[5] Yang C S, Xu S H, Chen D, Zhang Y F, Zhao Q L, Li C, Zhou K J, Feng Z M, Gan J L, Yang Z M 2016 J. Opt. 18 055801
Google Scholar
[6] Kane T 1990 IEEE Photonics Technol. Lett. 2 244
Google Scholar
[7] 李灿 2015 博士学位论文 (广州: 华南理工大学)
Li C 2015 Ph. D. Dissertation (Guangzhou: South China University of Technology) (in Chinese)
[8] Tröbs M 2005 Ph. D. Dissertation (Hannover: Leibniz University Hannover)
[9] Nicklaus K, Herding M, Wang X, Beller N, Fitzau O, Giesberts M, Herper M, Barwood G P, Williams R A, Gill P, Koegel H, Webster A, Gohlke M 2014 Proceedings of the International Conference on Space Optics 10563 105632T
Google Scholar
[10] Dahl K, Cebeci P, Fitzau O, Giesberts M, Greve C, Krutzik M, Peters A, Pyka S, Sanjuan J, Schiemangk M, Schuldt T, Voss K, Wicht A 2018 Proc. SPIE 11180 111800C
Google Scholar
[11] Milyukov V 2020 Astron. Rep. 64 1067
Google Scholar
[12] Gong Y G, Jun Luo, Wang B 2021 Nat. Astron. 5 881
Google Scholar
[13] Luo Z R, Guo Z K, Jin G, Wu Y L, Hu W R 2020 Results Phys. 16 102918
Google Scholar
[14] Luo Z R, Wang Y, Wu Y L, Hu W R, Jin G 2021 Prog. Theor. Exp. Phys. 2021 05A108
Google Scholar
[15] Foster S B, Tikhomirov A E. 2010 IEEE J. Quantum Electron. 46 734
Google Scholar
[16] Guiraud G, Traynor N, Santarelli G 2016 Opt. Lett. 41 4040
Google Scholar
[17] 杨中民, 徐善辉 2017 单频光纤激光器(北京: 科学出版社) 第59页
Yang Z M, Xu S H 2017 Single-frequency Fiber Laser (Beijing: Science Press) p59 (in Chinese)
[18] Hui R Q, O’Sullivan M 2009 Fiber Optic Test & Measurement (California: Elsevier Academic Press) pp263–264
[19] 李番, 王嘉伟, 高子超, 李健博, 安炳南, 李瑞鑫, 白禹, 尹王保, 田龙, 郑耀辉 2022 物理学报 17 209501
Google Scholar
Li F, Wang J W, Gao Z C, Li J B, An B N, Li R X, Bai Y, Yin W B, Tian L, Zheng Y H 2022 Acta Phys. Sin. 17 209501
Google Scholar
[20] Willke B 2010 Laser Photonics Rev. 4 780
Google Scholar
[21] Tröbs M, Heinzel G 2006 Measurement 39 120
Google Scholar
[22] Möller L 1998 IEEE J. Quantum Electron. 34 1554
Google Scholar
[23] Xue M Y, Gao C X, Niu L Q, Zhu A L, Sun C D 2020 Appl. Opt. 59 2610
Google Scholar
-
-
[1] Numata K, Yu A, Jiao H, Merritt S, Micalizzi F, Fahey M, Camp J, Krainak M 2019 Proceedings of SPIE: Solid State Lasers XXIII: Technology and Devices 10896 108961H
Google Scholar
[2] 罗子人, 白姗, 边星, 陈葛瑞, 董鹏, 董玉辉, 高伟, 龚雪飞, 贺建武, 李洪银, 李向前, 李玉琼, 刘河山, 邵明学, 宋同消, 孙保三, 唐文林, 徐鹏, 徐生年, 杨然, 靳刚 2013 力学进展 43 415
Google Scholar
Luo Z R, Bai S, Bian X, Chen G R, Dong P, Dong Y H, Gao W, Gong X F, He J W, Li H Y, Li X Q, Li Y Q, Liu H S, Shao M X, Song T X, Sun B S, Tang W L, Xu P, Xu S N, Yang R, Jin G 2013 Adv. Mech. 43 415
Google Scholar
[3] Peterseim M, Brozek O S, Danzmann K, Freitag I, Rottengatter P, Tünnermann A, Welling H 1998 AIP Conf. Proc. 456 148
Google Scholar
[4] O’Brien S, Welch D, Parke R, Mehuys D, Dzurko K, Lang R, Waarts R, Scifres D 1993 IEEE J. Quantum Electron. 29 2052
Google Scholar
[5] Yang C S, Xu S H, Chen D, Zhang Y F, Zhao Q L, Li C, Zhou K J, Feng Z M, Gan J L, Yang Z M 2016 J. Opt. 18 055801
Google Scholar
[6] Kane T 1990 IEEE Photonics Technol. Lett. 2 244
Google Scholar
[7] 李灿 2015 博士学位论文 (广州: 华南理工大学)
Li C 2015 Ph. D. Dissertation (Guangzhou: South China University of Technology) (in Chinese)
[8] Tröbs M 2005 Ph. D. Dissertation (Hannover: Leibniz University Hannover)
[9] Nicklaus K, Herding M, Wang X, Beller N, Fitzau O, Giesberts M, Herper M, Barwood G P, Williams R A, Gill P, Koegel H, Webster A, Gohlke M 2014 Proceedings of the International Conference on Space Optics 10563 105632T
Google Scholar
[10] Dahl K, Cebeci P, Fitzau O, Giesberts M, Greve C, Krutzik M, Peters A, Pyka S, Sanjuan J, Schiemangk M, Schuldt T, Voss K, Wicht A 2018 Proc. SPIE 11180 111800C
Google Scholar
[11] Milyukov V 2020 Astron. Rep. 64 1067
Google Scholar
[12] Gong Y G, Jun Luo, Wang B 2021 Nat. Astron. 5 881
Google Scholar
[13] Luo Z R, Guo Z K, Jin G, Wu Y L, Hu W R 2020 Results Phys. 16 102918
Google Scholar
[14] Luo Z R, Wang Y, Wu Y L, Hu W R, Jin G 2021 Prog. Theor. Exp. Phys. 2021 05A108
Google Scholar
[15] Foster S B, Tikhomirov A E. 2010 IEEE J. Quantum Electron. 46 734
Google Scholar
[16] Guiraud G, Traynor N, Santarelli G 2016 Opt. Lett. 41 4040
Google Scholar
[17] 杨中民, 徐善辉 2017 单频光纤激光器(北京: 科学出版社) 第59页
Yang Z M, Xu S H 2017 Single-frequency Fiber Laser (Beijing: Science Press) p59 (in Chinese)
[18] Hui R Q, O’Sullivan M 2009 Fiber Optic Test & Measurement (California: Elsevier Academic Press) pp263–264
[19] 李番, 王嘉伟, 高子超, 李健博, 安炳南, 李瑞鑫, 白禹, 尹王保, 田龙, 郑耀辉 2022 物理学报 17 209501
Google Scholar
Li F, Wang J W, Gao Z C, Li J B, An B N, Li R X, Bai Y, Yin W B, Tian L, Zheng Y H 2022 Acta Phys. Sin. 17 209501
Google Scholar
[20] Willke B 2010 Laser Photonics Rev. 4 780
Google Scholar
[21] Tröbs M, Heinzel G 2006 Measurement 39 120
Google Scholar
[22] Möller L 1998 IEEE J. Quantum Electron. 34 1554
Google Scholar
[23] Xue M Y, Gao C X, Niu L Q, Zhu A L, Sun C D 2020 Appl. Opt. 59 2610
Google Scholar
计量
- 文章访问数: 697
- PDF下载量: 22
- 被引次数: 0