Loading [MathJax]/jax/output/SVG/jax.js

留言板

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

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

平行陶瓷棒介质阻挡放电分段模型等效参数预测

刘星亮 邱祁 王若宇 邓焰 何湘宁

高磊, 曾勇虎, 汪连栋, 等. 宽带成像雷达间歇采样转发干扰应用策略[J]. 强激光与粒子束, 2018, 30: 053203. doi: 10.11884/HPLPB201830.170430
引用本文: 刘星亮, 邱祁, 王若宇, 等. 平行陶瓷棒介质阻挡放电分段模型等效参数预测[J]. 强激光与粒子束, 2019, 31: 040007. doi: 10.11884/HPLPB201931.180385
Gao Lei, Zeng Yonghu, Wang Liandong, et al. Application strategy for intermittent sampling repeater jamming to wideband imaging radar[J]. High Power Laser and Particle Beams, 2018, 30: 053203. doi: 10.11884/HPLPB201830.170430
Citation: Liu Xingliang, Qiu Qi, Wang Ruoyu, et al. Equivalent parameters prediction of dielectric barrier discharge piecewise model with parallel ceramic rods[J]. High Power Laser and Particle Beams, 2019, 31: 040007. doi: 10.11884/HPLPB201931.180385

平行陶瓷棒介质阻挡放电分段模型等效参数预测

doi: 10.11884/HPLPB201931.180385
详细信息
    作者简介:

    刘星亮(1990-), 男,博士,主要从事高频高压特种电源研究;Liuxingliang@zju.edu.cn

  • 中图分类号: TM46

Equivalent parameters prediction of dielectric barrier discharge piecewise model with parallel ceramic rods

  • 摘要: 为实现对介质阻挡放电负载模型等效参数的有效预测,引入放电区域面积作为中间变量。以平行陶瓷棒为负载,通过Maxwell有限元仿真得到恒压静电场下的负载等效电容与放电区域的对应关系。结合恒压静电场下的电场分布,提出了放电区域逐步扩张下的气隙首次击穿电压、负载外加电压峰值及气隙放电维持电压的预估方法。进而,根据李萨如图形法求得各个工作点的功率。至此,建立起了放电区域同各个等效参数间的量化关系,并实现了对参数的预测。最后,对气隙间距为1,3,4 mm的工况进行了实验验证。实验结果表明:气隙放电维持电压预测值在局部变化趋势上与实测值存在一定差异;放电功率、负载外加电压峰值预测结果与实测值较为吻合。
  • 间歇采样转发干扰是一种近十年来受到较多关注的新型干扰方法,该方法通过对雷达信号进行低速率的间歇采样处理,巧妙利用对雷达信号的间歇性“欠采样”处理技术,可以产生一串相干假目标[1]。这种干扰方式对信号形式为线性调频信号的逆合成孔径雷达(ISAR),具有较好的干扰效果,无论是匹配滤波接收处理还是“去斜”接收处理[2],都能在高分辨距离像处理结果中形成一串相干假目标[3],达到对目标真实散射结构压制和欺骗的效果。从干扰效果上看,间歇采样转发干扰在真实目标的高分辨距离像中叠加了一串相干假目标,这一串假目标会导致在二维图像上形成沿距离维均匀分布,幅度由主假目标中心向两侧逐渐下降,沿方位维幅度快速下降的一系列条状假目标[4-5]。但在结合后续成像处理过程分析后,可以发现,在成像结果中,虽然真实目标和假目标叠加在一起,但这一串相干假目标在事实上保证了真实目标部分的成像质量。出现上述现象的主要原因是,在实际成像观测场景下,即使目标仅进行小角度转动,目标上的散射结构在高分辨距离像中的分布情况也会发生改变,从而使得包络对齐效果并不理想[2, 6]。在随队干扰方式下,在较短的观测时间内,通过间歇采样转发干扰产生的一串虚假散射点在高分辨距离像中分布情况保持基本不变,从而作为相邻回波高分辨距离像中强相关结构,保证了包络对齐的效果。针对常规间歇采样转发干扰在事实上提升了ISAR成像质量这一有违于干扰初衷的现象,本文考虑从应用策略层面上改变间歇采样转发干扰的使用方式,通过破坏成像处理中的包络对齐步骤,试图达到降低ISAR成像质量的目标。

    设理想间歇采样脉冲信号p1(t)为单位幅度的矩形包络脉冲串信号,表达式为

    p1(t)=rect(tτ)+n=δ(tnTs1) (1)

    式中:rect(t/τ)为矩形窗函数;τ为采样脉冲持续时间;Ts1为间歇采样信号周期;“*”表示卷积处理;δ(·)为单位冲激函数。矩形窗函数满足rect(t/τ)=1, |t/τ|≤0.5;rect(t/τ)=0, |t/τ|>0.5。

    设宽带雷达发射的线性调频信号中心频率为f0,脉冲持续时间为Tp,线性调频率为γ,信号带宽B=γTp, 脉冲重复间隔为TPRI,则第m(m=0, 1, …)个发射信号脉冲可以写为

    s0(ˆt,tm)=rect(ˆt/Tp)exp(j2π(f0t+γˆt2/2)) (2)

    式中:ˆt是计量电波传播的快时间;t=ˆt+tm为全时间;tm=mTPRI,表示计量脉冲发射时刻的慢时间。

    考虑间歇采样转发干扰场景模型如图 1所示,图中RO表示从雷达位置到目标坐标系原点的矢量,满足‖RO‖=R0ROO表示从目标坐标系原点到干扰机位置的矢量。图中假设干扰机位于目标上随目标运动,坐标位置为(xJ, yJ),目标相对于雷达运动等效的转动角速度为ω(rad/s)。在目标小角度转动情况下,信号从雷达天线经干扰机天线返回的信号延时可近似为

    τm=2(R0+yJ+xJωtm)/c (3)
    图  1  存在干扰情况下的ISAR成像观测模型
    Figure  1.  Observation model for inverse synthetic aperture radar(ISAR) imaging when jamming exists

    式中:c为光速。干扰机截获到雷达信号后,开始对信号进行采样并转发。设干扰机固有延迟为τr,则干扰信号由于采样转发和系统采样过程引起的时间延迟为τs=τr+τ,从而间歇采样转发形成的单位幅度转发干扰信号可以表示为

    sJ(ˆt,tm)=p1(t)s0(ˆtτmτs,tm)=[rect(tτ)+n=δ(tnTs1)]×exp[j2π(f0(tτΔJ)+γ2(ˆtτΔJ)2)] (4)

    式中:τΔJ=τs+τm为总的时间延迟量,由干扰机到雷达的距离、间歇采样脉冲持续时间及系统延迟共同决定。

    在ISAR中,可以采取匹配滤波处理或解线调频(Dechirping)处理两种方式完成回波信号的接收处理。这里考虑Dechirping处理方式,间歇采样转发干扰信号与参考信号按照Dechirping方式完成混频处理后,忽略残余视频相位,可以写为

    sJI(ˆt,tm)=[rect(tτ)+n=δ(tnTs1)]exp[j2πγ(trefτΔJ)ˆt]exp[j2πf0(trefτΔJ)] (5)

    式中:tref表示参考时间。相应地,干扰信号Dechirping处理结果经距离向傅里叶变换处理得到高分辨距离像,可以表示为

    sJF(fr,tm)=Fi(sJI(ˆt,tm))=+n=τfs1Tpsinc(nfs1τ)sinc(Tp(frnfs1γ(τΔJtref)))exp(j2πf0(trefτΔJ)) (6)

    式中:sinc(x)=sin(πx)/(πx)。注意到式(5),两个指数项实际上与没有干扰情况下目标回波的Dechirping结果具有相同的形式,差异主要体现在前面的[]中对应的两部分,矩形函数表示加窗,而求和项对应了周期性冲激函数。相应地在进行傅里叶变换之后,矩形函数对应了sinc函数,而周期性冲激函数则对应了频域的周期延拓,亦即式(6)中第二个sinc函数项中减去nfs1项。

    根据式(6)的特点可知,在高分辨距离像中,间歇采样转发干扰会形成一串等间隔排列的虚假散射点,其间隔由fs1决定,而幅度则由式(6)中的两个sinc函数项共同决定。

    在间歇采样转发干扰释放的场景下,干扰信号与真实目标回波信号叠加在一起,在散射点模型下,目标回波信号经Dechirping处理后,以目标上某个延时为τΔ的散射点为例,其对应的中频信号可以写为

    sTI(ˆt,tm)=rect(tTp)exp[2πγ(trefτΔ)ˆt]exp[j2πf0(trefτΔ)] (7)

    相应地,傅里叶变换得到高分辨距离像如下

    sTF(fr,tm)=Fi(sTI(ˆt,tm))=Tpsinc(Tp(frγ(τΔtref)))exp(j2πf0(trefτΔ)) (8)

    根据式(6),干扰信号在高分辨距离像中对应了一串均匀分布的虚假散射点,这些虚假散射点的位置为

    frJ=nfs1+γ(τΔJtref) (9)

    式中:τΔJ=τs+τm=τr+τ+τm,这里τm是由干扰机位置决定的延时常量,在较短的观测时间内,可以认为是常量;而fs1对应间歇采样周期,按照目前的常规应用策略,为常数。

    因此,在较短的观测时间内,不同脉冲回波的高分辨距离像中,虚假散射点序列均具有相同的排列方式,且虚假散射点序列内部的间隔相同,也就是具有很高的相似度。

    根据经典ISAR成像算法,在处理多帧回波数据得到高分辨距离像之后,需要进行平动补偿,从而使目标运动满足转台成像模型。平动补偿包括包络对齐和初相校正两部分,其中包络对齐是基础,其目标是将全体高分辨距离像中对应相同散射点的散射分量移至相同距离单元中。目前获得广泛应用的包络对齐方法是互相关法及其改进方法[7-8],如积累相关方法和全局最优对齐法等[9-10],这类方法的基础是相邻脉冲高分辨距离像中散射结构分布具有相似性,也就是高分辨距离像包络具有相似性。

    在实际成像观测中,目标真实散射点对应的高分辨距离像包络随着观测角度变化,通常会发生小幅度改变,尤其是在带宽较大或目标方位维尺寸过大的情况下[2]。而由前面对间歇采样转发干扰处理结果分析可知,其所产生的虚假散射点序列在前后脉冲高分辨距离像中,具有相同的位置间隔。

    根据上述分析,在真实目标回波中叠加间歇采样转发干扰信号后,高分辨距离像处理结果中,真实目标散射结构与虚假散射点序列叠加在一起。由于虚假散射点序列数量众多,且具有高度相似性,从而使得相邻回波高分辨距离像之间的相似性得到提高,实际上保证了包络对齐的效果,从而在一定程度上提高了目标真实结构部分的成像质量,这与干扰的目标是相悖的。也就是说,常规的间歇采样转发干扰信号的存在,虽然叠加了一些虚假散射点到成像结果中,但在一定程度上也提高了真实目标部分的成像质量。

    针对间歇采样转发干扰在一定程度上提高目标部分成像质量的问题,其根本原因是干扰信号形成的虚假散射点序列在成像观测过程中保持相同的分布结构和间隔,且这个特点在包络对齐过程中发挥了作用,提高了包络对齐的质量。根据上述对问题原因的分析,一个较为直接的解决思路是:针对包络对齐方法特点,改变虚假散射点序列在相邻脉冲回波对应的高分辨距离像中的分布情况,使得相邻脉冲回波的相似性变弱。由前述分析,虚假散射点序列位置可由fr=nfs1+γ(τΔJ-tref)描述,而由这个式子可知,虚假散射点之间的间隔由fs1决定,也就是由间歇采样信号周期Ts1决定。如果干扰机工作时,不断改变这个间歇采样信号周期,则能够使得相邻脉冲回波中干扰信号对应的虚假散射点序列具有不同的分布间隔,从而使得其相似性变弱,从而达到降低包络对齐质量的目的。实际应用时,可以考虑预设多个间歇采样周期参数,将预设的多个间歇采样周期参数编入伪随机序列中,按照伪随机序列设定的顺序,选择使用不同的间歇采样周期,从而破坏相邻回波高分辨距离像的相关性。

    目标模型如图 2(a)所示,是由19个散射点组成的目标模型,图中给出了全部散射点在目标坐标系下的坐标信息,每个散射点间隔为2 m。干扰机设置在目标模型的水平对称轴上左起第五个散射点上。在雷达坐标系下,雷达位于坐标系原点,目标坐标为(0 m,8000 m),目标在雷达坐标系下沿x轴方向水平运动,相应的运动速度为vx=230 m/s,vy=0 m/s。雷达信号与干扰信号参数如表 1表 2所示。图 2(b)中给出了没有干扰情况下的ISAR成像结果。按照表 1的参数,有1281帧回波数据,在实际处理时,进行了5倍抽取,实际用于成像处理的脉冲数为257帧。

    图  2  目标散射点模型与无干扰情况下成像结果
    Figure  2.  Scatter model of target and its imaging result when no jamming exists
    表  1  雷达工作信号参数
    Table  1.  Parameters for radar signal
    carrier frequency/GHz bandwidth/GHz pulse duration/μs pulse repetition frequency/Hz accumulation time/s
    10 1 10 500 2.56
    下载: 导出CSV 
    | 显示表格
    表  2  干扰信号参数
    Table  2.  Parameters for jamming signal
    duration time of sampling pulse/μs period of intermittent sampling signal/μs intermittent sampling duty ratio amplitude of intermittent sampling signal
    0.03 0.1 0.3 7
    下载: 导出CSV 
    | 显示表格

    在给出具体结果之前,说明一下对比成像处理结果所使用的量化参数。前面提及,变采样周期策略的主要目的是降低相邻回波脉冲高分辨距离像包络间的相似性,降低包络对齐的质量,进而降低成像效果。因此,拟定义包络对齐之后的高分辨距离像包络质量评估指标和ISAR图像指标,量化对比两种策略的干扰效果。

    受文献[11-12]中给出的评估指标启发,本文拟应用的高分辨距离像(HRRP)包络质量评估指标如下:设单次成像观测所应用的高分辨距离像包络矩阵为H(p, q),其中p=1, 2, …, P; q=1, 2, …, QP为距离单元数,Q为回波脉冲数。

    MH(q,l)=pp=1|H(p,q)||H(p,l)|Pp=1|H(p,q)|2Pp=1|H(p,l)|2,q,l=1,2,,Q (10)

    注意到MH(q, l)的物理含义是第q次回波与第l次回波高分辨距离像包络间的相关系数,遍历q, l则可以得到Q×Q维矩阵MH

    定义1. 相邻相关系数序列:将l=q+1,q=1, 2, …, Q-1的特殊序列定义为相邻相关系数序列,并记为MHN(q),它实际是矩阵MH中仅挨着主对角线的一条对角线元素,可写成

    MHN(q)=MH(q,q+1),q=1,2,,Q1 (11)

    定义2. 全局相关系数:将MH主对角线右上方全体元素的均值(不含主对角线元素)定义为高分辨距离像包络的全局相关系数,并记为GM,即可以写成

    GM=2Q2QQ1q=1Ql=q+1MH(q,l) (12)

    按照式(10)的特点,两帧回波高分辨距离像包络相似性越强,其相关系数越大,极端情况两帧回波高分辨距离像完全一致,则其相关系数为1。

    对于ISAR像,应用图像熵指标评估图像质量。对于图像I(p, q),1≤pP, 1≤qQ,图像熵定义为

    H(I)=Qq=1Pp=1|I(p,q)|2s(I)ln[|I(p,q)|2s(I)] (13)

    其中s(I)=Qq=1Pp=1|I(p,q)|2为图像能量。按照图像分析理论,熵值越小,成像质量越好。

    首先给出常规策略下的成像结果,如图 3所示,图 3(a)给出了ISAR像结果,图 3(b)给出了包络对齐之后的高分辨距离像包络的俯视图。这里所应用的成像处理方法是经典的距离-多普勒方法,包络对齐方法为全局相关法[13],初相校正方法为多散射点中心法[2]

    图  3  常规干扰策略下成像处理结果
    Figure  3.  Imaging results when normal jamming strategy was applied

    再给出变采样周期策略下结果,这里选择了两个间歇采样周期,第一个周期按照表 2给出的参数,第二个周期为第一个周期的一半,间歇采样占空比保持不变,幅度保持不变。相应的处理结果如图 4所示,图 4(a)为ISAR像结果,图 4(b)为高分辨距离像俯视图。对比图 3图 4的ISAR成像结果可知,在交替改变间歇采样周期参数后,ISAR成像结果出现显著下降。在图 3(a)中,虽然干扰信号产生的虚假散射点较强,但与图 2(b)对比可知,真实目标结构部分仍然实现了良好聚焦,也就是包络对齐部分是顺利完成的;而图 4(a)图 2(b)对比可以发现,此时真实目标结构部分未能实现良好聚焦。相应的图 3(a)的熵值为7.018 9,小于图 4(a)的熵值8.253 1,与直观判断相同。另外,图 2(b)对应的无干扰情况下成像结果的熵值为6.572 3,比两种干扰策略下的值都小。

    图  4  变采样周期策略下成像处理结果
    Figure  4.  Imaging results when changing sampling period strategy was applied

    图 3(b)图 4(b)的对比不够直观,下面先给出两种策略下相邻相关系数序列,如图 5所示,这项指标表明变采样周期策略下,相邻高分辨距离像之间的相似性受到严重破坏。相应地,全局相关系数结果分别为:常规策略0.781 3,变采样周期策略0.565 3。

    图  5  两种干扰策略下高分辨距离像包络相邻相关系数序列对比
    Figure  5.  Comparison of correlation coefficient for adjacent HRRP under two jamming strategies

    通过定量和定性对比两种干扰策略下成像处理结果,可以发现,在采取了变采样周期策略后,相对于常规干扰策略,成像质量显著下降,从而提高了干扰的效果。在上述仿真中,变采样周期策略选择了最简单的选择两个采样周期,交替使用的方式完成了仿真验证,显然,在采取更加复杂的采样周期变化方式后,可以更加有效地降低成像质量。

    本文首先从常规策略下间歇采样转发干扰原理出发,分析指出了该策略存在的问题:由于前后脉冲中干扰信号产生的虚假散射点序列的相似性,提高了包络对齐的质量,从而在一定程度上提高了目标真实结构部分的成像效果。针对此问题产生的原因,提出了在应用间歇采样转发干扰中,按照预设随机序列改变间歇采样周期的应用新策略。通过仿真实验,定量对比了常规策略和变采样周期策略下间歇采样转发干扰对ISAR成像处理过程的影响情况,结果显示,本文所提出的变采样周期策略干扰方式能够更加有效地降低ISAR成像质量,从而验证了变采样周期策略的有效性。

  • 图  1  负载等效模型及李萨如图

    Figure  1.  Load equivalent circuit and its Lissajous figures

    图  2  带负载的驱动源电路及DCM与CCM临界工作模式下的关键波形

    Figure  2.  Resonant converter with load and its key waveforms in the critical mode between DCM and CCM

    图  3  负载静电场分布及基于几何结构的参数等效方法示意图

    Figure  3.  Schematic diagram of electrostatic field distribution and equivalent parameter extraction methods based on geometric structure

    图  4  不同气隙间距下,放电阶段模型等效电容同放电区域的关系

    Figure  4.  Equivalent capacitances versus discharge area with different air gap distances

    图  5  放电延伸曲线上不同位置对应的均匀系数

    Figure  5.  Uniformity coefficient versus different position along the discharge expansion curve

    图  6  放电延伸曲线上不同位置对应的气隙电压占比

    Figure  6.  Air gap voltage versus applied voltage at different position

    图  7  放电区域扩张所需外加电压峰值的变化曲线

    Figure  7.  Relationship between applied voltage peak and the discharge area

    图  8  放电扩张到不同位置时的放电维持电压

    Figure  8.  Discharge maintaining voltage corresponding to different discharge area

    图  9  不同气隙间距下负载外加电压峰值、放电维持电压及放电功率的实验值同预测值比较

    Figure  9.  Experiment and prediction results of applied voltage peak value, gap discharge maintaining voltage, and discharge power with different air gap distances

    表  1  不放电阶段不同气隙间距下的的等效电容参数

    Table  1.   Equivalent capacitances during non-discharge stage with different air gap distances

    d/mm Cd-ch/pF Cg-ch/pF Cch/pF Cch-exp/pF
    1 120.06 53.79 37.15 36.27
    3 142.14 24.47 20.87 23.02
    4 148.43 20.30 17.86 21.62
    下载: 导出CSV
  • [1] 宋颖. 大气压非平衡等离子体在杀菌中的应用研究[D]. 大连: 大连理工大学, 2014.

    Song Ying. Inactivation applications of atmospheric pressure nonequilibrium plasmas. Dalian: Dalian University of Technology, 2014
    [2] Kogelschatz U. Dielectric-barrier discharges: Their history, discharge physics, and industrial applications[J]. Plasma Chemistry and Plasma Processing, 2003, 23(1): 1-46. doi: 10.1023/A:1022470901385
    [3] 胡小吐. AC/DC流光放电等离子体烟气脱硫实验研究[D]. 北京: 北京交通大学, 2007.

    Hu Xiaotu. Experimental research of AC/DC streamer plasmas in flue gas desulfurization. Beijing: Beijing Jiaotong University, 2007
    [4] Lieberman M A, Lichtenberg A J. Principles of plasma discharges and materials processing[M]. Hoboken: John Wiley & Sons, 2005.
    [5] 郝世强. 介质阻挡放电的能量压缩机理、实现及优化[D]. 杭州: 浙江大学, 2016.

    Hao Shiqiang. Mechanism, implementation and optimization of dielectric barrier discharge energy compression. Hangzhou: Zhejiang University, 2016
    [6] Kinnares V, Hothongkham P. Circuit analysis and modeling of a phase-shifted pulsewidth modulation full-bridge-inverter-fed ozone generator with constant applied electrode voltage[J]. IEEE Transactions on Power Electronics, 2010, 25(7): 1739-1752. doi: 10.1109/TPEL.2010.2042075
    [7] Eid A, Takashima K, Mizuno A. Experimental and simulation investigations of DBD plasma reactor at normal environmental conditions[J]. IEEE Transactions on Industry Applications, 2014, 50(6): 4221-4227. doi: 10.1109/TIA.2014.2315496
    [8] Guo Tangtang, Liu Xingliang, Hao Shiqiang, et al. Prediction of equivalent electrical parameters of dielectric barrier discharge load using a neural network[J]. Plasma Science and Technology, 2015, 17(3): 196-201. doi: 10.1088/1009-0630/17/3/05
    [9] Alonso J M, Valdés M, Calleja A J, et al. High frequency testing and modeling of silent discharge ozone generators[J]. Ozone Science & Engineering, 2003, 25(5): 363-376.
    [10] 王静, 蔡忆昔, 王军, 等. 介质阻挡放电等效电容的测量与分析[J]. 高电压技术, 2008, 34(2): 264-266, 308. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYJ200802013.htm

    Wang Jing, Cai Yixi, Wang Jun, et al. Measurement and analysis of equivalent capacitance in dielectric barrier discharge. High Voltage Engineering, 2008, 34(2): 264-266, 308 https://www.cnki.com.cn/Article/CJFDTOTAL-GDYJ200802013.htm
    [11] 王辉, 方志, 邱毓昌, 等. 介质阻挡放电等效电容变化规律的研究[J]. 绝缘材料, 2005, 38(1): 37-40. https://www.cnki.com.cn/Article/CJFDTOTAL-JYCT200501012.htm

    Wang Hui, Fang Zhi, Qiu Yuchang, et al. On the changing of equivalent capacitance in dielectric barrier discharge. Insulating Materials, 2005, 38(1): 37-40 https://www.cnki.com.cn/Article/CJFDTOTAL-JYCT200501012.htm
    [12] Xiao Dengming. Gas discharge and gas insulation[M]. Heidelberg: Springer-Verlag, 2016.
    [13] 张红. 高电压技术[M]. 北京: 中国电力出版社, 2009.

    Zhang Hong. High voltage technique. Beijing: China Electric Power Press, 2009
  • 期刊类型引用(2)

    1. 袁雅婷,曾冰倩,胡桐宁,李小飞,樊宽军. 可调电子束注入器调配调谐结构的设计及分析(英文). 强激光与粒子束. 2022(04): 129-135 . 本站查看
    2. 邱宇帆,李胜波,郑新建,符升平,白凤民. 多参数耦合下电子枪静电聚焦特性分析. 真空科学与技术学报. 2021(11): 1094-1100 . 百度学术

    其他类型引用(0)

  • 加载中
图(9) / 表(1)
计量
  • 文章访问数:  840
  • HTML全文浏览量:  297
  • PDF下载量:  107
  • 被引次数: 2
出版历程
  • 收稿日期:  2018-12-29
  • 修回日期:  2019-02-22
  • 刊出日期:  2019-04-15

目录

/

返回文章
返回