Loading [MathJax]/jax/output/HTML-CSS/jax.js
高级搜索

留言板

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

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

基于分数阶傅里叶变换的低信噪比线性调频信号参数快速估计算法

刘利民 李豪欣 李琦 韩壮志 高振斌

耿友林, 解成博, 尹川, 郭兰图, 王先义. 基于卡尔曼滤波的接收信号强度指示差值定位算法[J]. 电子与信息学报, 2019, 41(2): 455-461. doi: 10.11999/JEIT180268
引用本文: 刘利民, 李豪欣, 李琦, 韩壮志, 高振斌. 基于分数阶傅里叶变换的低信噪比线性调频信号参数快速估计算法[J]. 电子与信息学报, 2021, 43(10): 2798-2804. doi: 10.11999/JEIT200973
Youlin GENG, Chengbo XIE, Chuan YIN, Lantu GUO, Xianyi WANG. Received Signal Strength Indication Difference Location Algorithm Based on Kalman Filter[J]. Journal of Electronics & Information Technology, 2019, 41(2): 455-461. doi: 10.11999/JEIT180268
Citation: Limin LIU, Haoxin LI, Qi LI, Zhuangzhi HAN, Zhenbin GAO. A Fast Signal Parameter Estimation Algorithm for Linear Frequency Modulation Signal under Low Signal-to-Noise Ratio Based on Fractional Fourier Transform[J]. Journal of Electronics & Information Technology, 2021, 43(10): 2798-2804. doi: 10.11999/JEIT200973

基于分数阶傅里叶变换的低信噪比线性调频信号参数快速估计算法

doi: 10.11999/JEIT200973
基金项目: 国家自然科学基金(61601496),河北省自然科学基金(F2019506037)
详细信息
    作者简介:

    刘利民:男,1971年生,教授,博士生导师,研究方向为雷达对抗、雷达信号处理

    李豪欣:女,1997年生,硕士,研究方向为雷达信号处理

    李琦:男,1974年生,教授,硕士生导师,研究方向为雷达信号处理、信息与编码理论

    韩壮志:男,1972年生,副教授,硕士生导师,研究方向为雷达信号处理、武器系统性能评估

    高振斌:男,1973年生,教授,硕士生导师,研究方向为雷达信号处理

    通讯作者:

    李琦 liqi@hebut.edu.cn

  • 中图分类号: TN957.51

A Fast Signal Parameter Estimation Algorithm for Linear Frequency Modulation Signal under Low Signal-to-Noise Ratio Based on Fractional Fourier Transform

Funds: The National Natural Science Foundation of China (61601496), The Nature Science Foundation of Hebei Province (F2019506037)
  • 摘要: 针对低信噪比线性调频信号参数估计精度低且运算量大的问题,该文提出一种基于高效分数阶傅里叶变换(FRFT)和分数阶频谱4阶原点矩的快速估计算法。该算法通过判断调频斜率的正负,以确定旋转阶次所在初始区间;进而应用高效FRFT获得初始旋转阶次;最终利用分数阶频谱4阶原点矩,进一步确定搜索区间和步长,实现精准搜索,从而满足参数精度的要求。实验结果表明,该算法尤其适合用于低信噪比情况下的线性调频(LFM)信号检测与参数的准确估计,而且运算量较低。
  • 电磁频谱监测[1]是实现无线电管理的技术手段之一。随着通信技术发展和应用普及,无线电频谱管理工作越来越复杂,为更加高效管理无线电资源,避免未知信号源对现有电磁环境干扰,根据频谱监测数据及时对信号进行定位和跟踪一直是电磁频谱研究热点之一。

    对信号源定位技术主要有以下几种实现方法,例如信号到达时间差(TDOA)[26]、信号到达角度(AOA)[5,6]、信号到达频率差(FDOA)[7,8]以及接收信号强度指示(RSSI)[911]等。其中TDOA, AOA, FDOA及其联合定位的定位精度较高,但需要监测设备拥有额外的硬件[10]。综合考虑,在频谱监测系统中频谱数据所包含的信号相关参数较少,且地面监测站接收机并没有都配备昂贵的测向设备等特点,因而RSSI定位技术是目前频谱监测中一种比较有效的方法。

    由于被监测信号无法控制且没有任何先验知识,只能通过对信号的被动监测,即接收与处理信号来估计信号源位置的情况,无法直接利用RSSI测距定位。本文提出了一种新的定位算法,即将监测站间的接收信号强度指示差值(RSSID)转换成信号源到两监测站的距离之比进行定位。此外,在已知信号源到监测站距离的情况下,距离之比定位(Localization with Ratio-Distance, LRD)比三边测量法(Trilateration)具有更加优越的定位性能[12]。结合LRD思想,本文提出了一种适用于频谱监测系统对信号源进行有效定位的方法,即RSSID定位算法。

    在研究信号源位置估算方法之前,有必要选择一种适合当前环境的无线电波传播的数学模型。本文采用适合解决市区无线电信号源位置估计的Egli电波传播模型[13],其经验公式如式(1)所示。

    L=78+20lgf+40lgd20ht20hr (1)

    式中,L为路径衰减,单位dB;f为信号频率,单位MHz;d为信号源到监测站的距离,单位km;ht为信号源发射天线高度,单位m;hr监测站接收天线的高度,单位m。

    发射功率与接收功率之间的关系为

    RSSI=PtL (2)

    式中,RSSI为监测站接收功率,单位dB;Pt为信号源发射功率,单位dB。

    为了消除信号源的相关参数,通过式(3)消去信号源发射天线高度ht以及信号源发射功率Pt

    RSSI0RSSIi=LiL0=40lgdid0+20lghr0hri (3)

    式中,RSSI0为主监测站接收功率,RSSIi为第i个次监测站接收功率,因此信号源到各次监测站距离与到主监测站的距离之比ki0

    ki0=did0=10RSSI0RSSIi20lghr0hri40 (4)

    由解析几何可知,信号源到两个监测站的距离之比(比值不等于1)为常数时,其轨迹是一个圆,此圆的圆心位于两监测站位置的连线上并且靠近接收信号大的监测站一侧。假设有3个监测站同时接收到同一信号电平,就可以得到2个距离比圆,得到2个可能的信号源位置解,再根据信号源辐射范围即可确定信号源位置。

    设主监测站位置为B0(x0,y0),接收天线的高度为h0,次监测站位置为B1(x1,y1), B2(x2,y2),接收天线高度分别为h1, h2,信号源的实际位置为S(xs,ys),两个距离比圆的圆心为C1(xc1,yc1), C2(xc2,yc2),所以根据距离之比和已知监测站位置,可得到

    ki0=10RSSI0RSSIi20lghr0hri40=did0=(xxi)2+(yyi)2(xx0)2+(yyi)2=λi,i=1,2 (5)

    整理式(5)得

    (xxiλix01λi)2+(yyiλiy01λi)2=λi(1λi)2[(x0xi)2+(y0yi)2]=λi(1λi)2d2i0 (6)

    式(6)中,di0为次监测站到主监测站的距离,可以得到距离比圆的信息。圆心坐标为

    xci=xiλix01λi,yci=yiλiy01λi (7)

    距离比圆半径:

    Ri=|ki01λi|di0 (8)

    因此,至少使用3个监测站才可得到2个距离比圆,两圆在监测站共同覆盖范围的交点即为目标位置。实际情况中,多个监测站进行监测时,可以得到多个距离比圆,由于误差的存在,多个距离比圆相交将产生多个可能的目标解,式(9)对应的目标序号j即为所要求的目标解:

    min[ni=0(xixj)2+(yiyj)2],j[0,n1] (9)

    本文算法的定位原理如图1所示。

    图 1  RSSID定位原理图

    从式(4)可以看出,信号源到次监测站与主监测站的距离之比与RSSI信号密切相关。但是,在信号的实际传播过程中RSSI信号容易受到电波传播环境影响,例如障碍物尺寸、天气的好坏以及接收节点所在位置的高度等,因此,需要对RSSI信号进行预处理得到平稳变化的RSSI信号。

    卡尔曼滤波是基于线性最小均方误差预测和线性递归更新的优化算法[14],根据前一时刻的预测值和当前时刻的测量值更新当前时刻的预测值,可以有效解决干扰噪声服从正态分布的滤波问题,适用于RSSI信号的滤波处理[1518]。在这项工作中,RSSI信号噪声被假定为高斯噪声,虽然这个假设并不完全准确,但实际测量中RSSI值的误差波动比较集中,此时卡尔曼滤波器可以在一定程度上解决信号漂移和冲击问题[18]。其次,卡尔曼滤波在求解时不需要储存大量的观测数据,并且当得到新的观测数据时,可随时算得新的参数滤波值,便于实时地处理观测数据。

    对RSSI信号进行“测量-预测-更新”,消除系统随机噪声,平滑RSSI信号。主要包括两个阶段:预测和更新。预测阶段为[16]

    k时刻状态预测:

    x(k|k1)=x(k1|k1) (10)

    协方差预测:

    P(k|k1)=P(k1|k1)+Q (11)

    式中,x(k1|k1)k1时刻RSSI信号的状态值,x(k|k)k时刻RSSI信号的状态值, x(k|k1)k时刻RSSI信号的预测值,P(k1|k1)是对应x(k1|k1)的系统误差协方差,P(k|k1)对应x(k|k1)的系统误差协方差,Q状态噪声协方差。递归估计更新阶段[16]

    卡尔曼增益:

    Kg(k)=P(k|k1)P(k|k1)+R (12)

    k时刻状态更新:

    x(k|k)=x(k|k1)+Kg(k)[y(k)x(k|k1)] (13)

    协方差更新:

    P(k|k1)=(1Kg(k)P(k|k1)) (14)

    式中,R是测量噪声协方差,Kg(k)k时刻的卡尔曼增益,y(k)k时刻RSSI信号的测量值。当对象模型足够准确且系统状态和参数不发生突变时,卡尔曼滤波性能较好。

    实际使用中发现,监测站距离信号源较近时,由于受到传播路径的影响较小,RSSI测量值相对较准确[19],因此选取RSSI值最大的监测站作为主监测站。根据式(5)可以求得信号源到其余各监测站与主监测站的距离之比。

    另外由信号源S(xs,ys)到达监测站Bi(xi,yi), i=1,2,···,N+1的欧式距离公式可得

    k2i0=d2id20=SBi2SB02=(xsxi)2+(ysyi)2(xsx0)2+(ysy0)2 (15)

    将式(15)整理展开可以得到

    (1k2i0)(xs2+ys2)+2(k2i0x0xi)xs+2(k2i0y0yi)ys =ki02(x02+y02)xi2yi2 (16)

    S=[xs,ys,xs2+ys2]为未知矢量建立方程,N个上述方程构成定位矩阵:

    GS=b (17)

    式中,

    G=[2(k210x0x1)2(k210y0y1)1k2102(k220x0x2)2(k220y0y2)1k2202(k2N0x0xN)2(k2N0y0yN)1k210],b=[k210(x20+y20)x21y21k220(x20+y20)x22y22k2N0(x20+y20)x2Ny2N] (18)

    于是得到了S的最小二乘估计值:

    S=(GTG)1GTb (19)

    则信号源位置S=S[1:3]

    本文使用Java编程实现RSSI和RSSID定位算法的仿真实验,并利用Python语言对实验结果进行了分析和处理。本文算法已使用Java编程成功应用于国内某频谱大数据监测平台。

    为分析本文算法的定位性能,在相同的仿真环境下进行以下仿真分析实验。仿真环境设置为4个监测站,分别位于20 km×20 km的正方形区域的4个顶点上,其坐标为(0, 0), (0, 20), (20, 0)和(20, 20)。另外,为对比本文算法与RSSI测距定位算法的定位性能,以下实验中已知信号源的相关参数,如频率、发射天线高度及发射功率等。

    分别验证接收信噪比SNR以及监测站间距对定位精度的影响,采用定位结果与真实值之间的均方根误差衡量定位精度,

    RMSE(u)=Ni=0uobj,iumea,iN (20)

    其中,N是信号源个数,uobj,i表示第i个信号源的真实位置,umea,i表示i个信号源的定位结果。参与比较的算法包括:传统RSSI定位算法,本文所提的RSSID定位算法。

    4.1.1   接收信噪比SNR对定位精度的影响

    为验证接收信噪比SNR对定位精度的影响,在相同的监测站位置配置下比较不同定位算法的定位性能,设置区域内随机点数N=400,接收信噪比SNR∈[–16, 20] dB。

    定位精度随信噪比不同的变化情况如图2所示,对比传统RSSI定位算法,当SNR∈[0, 20] dB时本文算法性能明显提升,这是因为即使估计的距离与真实距离相差较大,距离的比例也没有太大的影响。当SNR≤–4 dB时,噪声在接收信号中起决定性作用,本文算法与RSSI定位算法误差均较大,无法正确估计信号源位置;当SNR≥–4 dB时,本文算法定位性能优于RSSI定位算法。

    图 2  随信噪比变化的定位性能
    4.1.2   监测站间距对定位精度的影响

    为验证监测站间距对定位精度的影响,设置信噪比SNR=16 dB,监测区域内信号源个数N=1000,监测站间距d[1,20] km。

    定位精度随监测站间距不同的变化情况如图3所示,对比传统RSSI定位算法,随着监测站间距的增大,本文算法定位误差明显较小,性能明显提升。

    图 3  随监测站间距变化的定位性能

    在此分析监测站接收信号信噪比为SNR=16 dB时,RSSID定位与RSSI定位的定位性能,以及卡尔曼滤波算法对RSSID定位的优化能力。

    4.2.1   数据预处理

    RSSID定位与RSSI定位的定位精度均受到信号接收强度的影响,因此有必要对RSSI信号进行预处理。本文采用卡尔曼滤波处理RSSI信号,如图4(a)图4(b)所示,卡尔曼滤波可以实现对RSSI信号的平滑处理,减少RSSI信号的测量误差。

    图 4  卡尔曼滤波效果图
    4.2.2   单目标定位

    为对比RSSID定位算法与RSSI定位算法对单个目标的定位性能,对区域内随机一个目标进行200次定位。两种不同算法的定位误差如图5(a)图5(b)所示,可以看出滤波前后,RSSID定位算法与RSSI定位算法都可以对信号源进行有效定位,RSSID定位算法均比RSSI定位算法定位误差较小。

    图 5  单目标定位误差对比

    为了避免仿真实验的随机性,对1个随机生成的信号源进行10000次定位,分析其误差,如表1所示。

    表 1  单目标定位10000次误差统计分析(km)
    是否预处理定位方法最大误差最小误差平均误差
    RSSI定位5.08020.02391.3993
    RSSID定位4.62240.00760.8527
    RSSI定位1.45370.22730.6249
    RSSID定位0.88010.00680.2683
    下载: 导出CSV 
    | 显示表格
    4.2.3   多目标定位

    为对比RSSID定位算法与RSSI定位算法对多个目标同时定位的定位性能,对区域内随机生成的25个信号源进行200次定位。两种不同算法的定位误差如图6(a)图6(b)所示,可以看出滤波前后,本文RSSID算法定位性能均明显优于RSSI定位算法。

    图 6  多目标定位平均误差对比

    为避免因为数据量较少出现偶然误差,对200个信号源进行10000次定位,分析每次对200个信号源定位的平均误差,统计结果如表2所示。

    表 2  多目标定位10000次平均误差统计分析(km)
    是否预处理定位方法最大误差最小误差平均误差
    RSSI定位1.86021.35991.5911
    RSSID定位1.10150.76200.9170
    RSSI定位1.53120.94701.1530
    RSSID定位0.82840.19480.2930
    下载: 导出CSV 
    | 显示表格

    本文提出了一种运用于频谱监测系统对监测信号有效定位的方法。在频谱监测系统中,被监测信号无法控制并且没有任何先验知识,只能通过对信号的被动监测、接收与处理来估计信号源的位置,此时不适用传统RSSI定位法,而本文RSSID算法是根据监测站地理位置以及各监测站接收信号强度对信号源定位的,不需要任何先验知识。因此,本文算法适用性更广,可用于频谱监测系统对监测到的未知信号进行实时定位。

    仿真分析表明,在信号相关参数已知的情况下,本文RSSID算法比传统RSSI定位法受到信噪比的影响较小,定位性能更加优越。另外引入卡尔曼滤波对RSSI值进行预处理可以极大的提高本文算法的定位性能。

  • 图  1  FRFT频谱幅度与阶次关系图

    图  2  归一化幅度对比图

    图  3  FRFT与W-V分布关系图

    图  4  LFM信号时频分布

    图  5  改进算法流程图

    图  6  算法检测性能曲线

    图  7  两种算法对比

    表  1  3种算法对比仿真结果

    估计方法允许阶次误差ˆk(Hz/t)ˆf(Hz)kerror(%)ferror(%)FRFT运算次数
    高效FRFT算法1537.20001580.353.72445.35653
    FRFT 2维搜索<0.01001104.00001511.110.40100.7417201
    FRFT 2维搜索<0.00101024.60001501.82.45800.12142001
    FRFT 2维搜索<0.00011003.90001501.50.39500.100720001
    改进搜索算法<0.0100940.20611510.25.97940.68007
    改进搜索算法<0.0010993.16241501.10.68380.075420
    改进搜索算法<0.00011002.30001500.80.22810.053143
    下载: 导出CSV
  • [1] DUAN Yu, WANG Jinzhen, SU Shaoying, et al. Detection of LFM signals in low SNR based on STFT and wavelet denoising[C]. 2014 International Conference on Audio, Language and Image Processing, Shanghai, China, 2014: 921–925.
    [2] YIN Qingbo, SHEN Liran, LU Mingyu, et al. Selection of optimal window length using STFT for quantitative SNR analysis of LFM signal[J]. Journal of Systems Engineering and Electronics, 2013, 24(1): 26–35. doi: 10.1109/JSEE.2013.00004
    [3] XU Fenfei, BAO Qinglong, CHEN Zengping, et al. Parameter estimation of multi-component LFM signals based on STFT+Hough transform and fractional fourier transform[C]. The 2nd IEEE Advanced Information Management, Communicates, Electronic and Automation Control Conference (IMCEC), Xi’an, China, 2018: 839–842.
    [4] ZHANG Zhichao. Linear canonical Wigner distribution based noisy LFM signals detection through the output SNR improvement analysis[J]. IEEE Transactions on Signal Processing, 2019, 67(21): 5527–5542. doi: 10.1109/TSP.2019.2941071
    [5] WU Yushuang and LI Xiukun. Elimination of cross-terms in the Wigner–Ville distribution of multi-component LFM signals[J]. IET Signal Processing, 2017, 11(6): 657–662. doi: 10.1049/iet-spr.2016.0358
    [6] ZHANG Zhichao. The optimal linear canonical Wigner distribution of noisy linear frequency-modulated signals[J]. IEEE Signal Processing Letters, 2019, 26(8): 1127–1131. doi: 10.1109/LSP.2019.2922510
    [7] HUANG Xiang, ZHANG Linrang, ZHANG Juan, et al. Efficient angular chirp-Fourier transform and its application to high-speed target detection[J]. Signal Processing, 2019, 164: 234–248. doi: 10.1016/j.sigpro.2019.06.011
    [8] YANG Tiantian, SHAO Jie, CHEN Yongliang, et al. Parameter estimation of multi component LFM signals based on nonlinear mode decomposition and FRFT[C]. The 10th International Conference on Advanced Computational Intelligence (ICACI), Xiamen, China, 2018: 204–209.
    [9] MIAO Hongxia, ZHANG Feng, and TAO Ran. Fractional Fourier analysis using the Möbius inversion formula[J]. IEEE Transactions on Signal Processing, 2019, 67(12): 3181–3196. doi: 10.1109/TSP.2019.2912878
    [10] LIU Yifei, ZHAO Yuan, ZHU Jun, et al. Iterative high-accuracy parameter estimation of uncooperative OFDM-LFM radar signals based on FRFT and fractional autocorrelation interpolation[J]. Sensors, 2018, 18(10): 3550. doi: 10.3390/s18103550
    [11] 赵兴浩, 邓兵, 陶然. 分数阶傅里叶变换数值计算中的量纲归一化[J]. 北京理工大学学报, 2005, 25(4): 360–364. doi: 10.3969/j.issn.1001-0645.2005.04.019

    ZHAO Xinghao, DENG Bing, and TAO Ran. Dimensional normalization in the digital computation of the fractional fourier transform[J]. Transactions of Beijing Institute of Technology, 2005, 25(4): 360–364. doi: 10.3969/j.issn.1001-0645.2005.04.019
    [12] 仇兆炀, 陈蓉, 汪一鸣. 基于FRFT的线性调频信号欠采样快速检测方法[J]. 电子学报, 2012, 40(11): 2165–2170.

    QIU Zhaoyang, CHEN Rong, and WANG Yiming. Fast detection of LFM signal based on FRFT and sub-nyquist sampling[J]. Acta Electronica Sinica, 2012, 40(11): 2165–2170.
    [13] ALDIMASHKI O and SERBES A. Performance of chirp parameter estimation in the fractional fourier domains and an algorithm for fast chirp-rate estimation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2020, 56(5): 3685–3700. doi: 10.1109/TAES.2020.2981268
    [14] 黄响, 唐世阳, 张林让, 等. 一种基于高效FRFT的LFM信号检测与参数估计快速算法[J]. 电子与信息学报, 2017, 39(12): 2905–2911.

    HUANG Xiang, TANG Shiyang, ZHANG Linrang, et al. A fast algorithm of LFM signal detection and parameter estimation based on efficient FRFT[J]. Journal of Electronics &Information Technology, 2017, 39(12): 2905–2911.
    [15] 宋耀辉, 黄仰超, 张衡阳, 等. 基于FRFT的多分量LFM信号检测与参数估计方法[J]. 北京航空航天大学学报, 2020, 46(6): 1221–1228. doi: 10.13700/j.bh.1001-5965.2019.0430

    SONG Yaohui, HUANG Yangchao, ZHANG Hengyang, et al. Multicomponent LFM signal detection and parameter estimation method based on FRFT[J]. Journal of Beijing University of Aeronautics and Astronautics, 2020, 46(6): 1221–1228. doi: 10.13700/j.bh.1001-5965.2019.0430
    [16] YIN Zhiping, ZHANG Dongchen, CHEN Weidong, et al. LFM signal detection using the origin moment of fractional spectrum[C]. 2008 9th International Conference on Signal Processing, Beijing, China, 2008: 191–194.
  • 期刊类型引用(12)

    1. 余修武,彭威,余员琴,刘永. 基于SSO的铀尾矿库无线传感器网络定位算法. 中国安全科学学报. 2023(04): 84-90 . 百度学术
    2. 束仁义,朱家兵,沈晓波,蔡俊,夏泐. 基于KF-BA-LSSVR的无线传感器网络节点定位研究. 重庆科技学院学报(自然科学版). 2022(03): 57-61 . 百度学术
    3. 吴红生,黄良永. 激光传输过程中的信号捕获与自动跟踪研究. 激光杂志. 2022(07): 139-143 . 百度学术
    4. 路韬,黄友朋,党三磊,张捷. 基于改进ELM的计量终端定位算法仿真. 计算机仿真. 2022(09): 377-383 . 百度学术
    5. 林逢春,谭文武. 人工智能技术的高速光纤通信链路智能选择研究. 激光杂志. 2022(12): 159-163 . 百度学术
    6. 赵书红,董绍武,白杉杉,高喆. 一种优化的频率驾驭算法研究. 电子与信息学报. 2021(05): 1457-1464 . 本站查看
    7. 张旭,姚善化. 未调制可见光室内定位. 湖南文理学院学报(自然科学版). 2021(03): 58-61+80 . 百度学术
    8. 孙小君,周晗,沈海滨,闫广明. 加权融合鲁棒增量Kalman滤波器. 电子与信息学报. 2021(12): 3680-3686 . 本站查看
    9. 郭永宁,李燕龙. 无人机无线电监测混合定位算法. 桂林电子科技大学学报. 2021(05): 362-367 . 百度学术
    10. 梁拥军. 基于卡尔曼滤波的无线电台电磁干扰源定位方法研究. 电子设计工程. 2020(09): 112-116 . 百度学术
    11. 孙小君,周晗,闫广明. 基于新息的自适应增量Kalman滤波器. 电子与信息学报. 2020(09): 2223-2230 . 本站查看
    12. 康颜爽,苏胜君,施伟斌,乐燕芬,韩承毅. 一种基于核岭回归与卡尔曼滤波的定位算法. 软件导刊. 2020(09): 74-77 . 百度学术

    其他类型引用(16)

  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  1158
  • HTML全文浏览量:  661
  • PDF下载量:  149
  • 被引次数: 28
出版历程
  • 收稿日期:  2020-11-12
  • 修回日期:  2021-08-16
  • 网络出版日期:  2021-08-27
  • 刊出日期:  2021-10-18

目录

/

返回文章
返回