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

留言板

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

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

光电镊—一种用途广泛的微纳操作工具

张帅龙 李恭 李凤刚 徐冰睿 李航 符荣鑫

耿友林, 解成博, 尹川, 郭兰图, 王先义. 基于卡尔曼滤波的接收信号强度指示差值定位算法[J]. 电子与信息学报, 2019, 41(2): 455-461. doi: 10.11999/JEIT180268
引用本文: 张帅龙, 李恭, 李凤刚, 徐冰睿, 李航, 符荣鑫. 光电镊—一种用途广泛的微纳操作工具[J]. 电子与信息学报, 2023, 45(12): 4566-4575. doi: 10.11999/JEIT221315
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: ZHANG Shuailong, LI Gong, LI Fenggang, XU Bingrui, LI Hang, FU Rongxin. Optoelectronic Tweezers — A Versatile Micro/Nano Operation Technique[J]. Journal of Electronics & Information Technology, 2023, 45(12): 4566-4575. doi: 10.11999/JEIT221315

光电镊—一种用途广泛的微纳操作工具

doi: 10.11999/JEIT221315
基金项目: 国家自然科学基金(62103050, 62105177, 21904008),国家海外高层次人才计划青年项目
详细信息
    作者简介:

    张帅龙:男,博士,教授,博士生导师,主要研究方向为生物微纳操作系统、光电镊技术、数字微流控技术、生物检测

    李恭:男,硕士生,研究方向为基于光电镊技术的微操控

    李凤刚:男,博士生,主要研究方向为基于光电镊技术的微操控及数字微流控

    徐冰睿:女,博士生,主要研究方向为光操控技术与微生物行为学

    李航:女,博士,副研究员,硕士生导师,主要研究方向为质谱分析与蛋白质组学

    符荣鑫:男,博士,副研究员,硕士生导师,主要研究方向为光谱计算与微流控生物传感器

    通讯作者:

    符荣鑫 furongxin@bit.edu.cn

  • 中图分类号: TP212; O436

Optoelectronic Tweezers — A Versatile Micro/Nano Operation Technique

Funds: The National Natural Science Foundation of China (62103050, 62105177, 21904008), The Recruitment Program of Global Experts and National Special Support Plan for High-level Talents
  • 摘要: 光电镊(OET)是一种基于光致介电泳效应的微尺度光操控技术,可在流体场、光电场、生物力场耦合的复杂环境下对微小目标进行精准操控,在细胞操作、微机械系统等领域有着重要的应用。光电镊技术可以单独使用或与其他技术协同使用,应用十分广泛。迄今为止,基于光电镊的研究主要集中在:微纳米材料的操作、组装和合成;单个细胞/分子的操作、分离和分析;细胞固有特性的分析和获取;细胞的电穿孔、融合和裂解;细胞封装生物材料和生物结构的制备;流体传输的光流体器件的开发。这些工作展示了光电镊技术优越的性能和独特的通用性和灵活性。该文系统地介绍了光电镊技术的现有应用,总结了该技术的应用前景、局限性及发展趋势。
  • 电磁频谱监测[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  3种光电镊设备[10]

    图  2  使用微生物运输微小物体[11]

    图  3  光诱导单细胞电穿孔[12]

    图  4  细胞在接收点进行光诱导细胞融合[13]

    图  5  细胞基于光电镊的1阶响应[14]

    图  6  利用收缩光圈富集细胞外囊泡[15]

    图  7  光电镊结合数字微流控技术[16]

    图  8  Beacon光电镊平台[17]

    图  9  基于光电镊技术印制微电路[19]

    图  10  可进行并行操作的光电镊微机器人[20]

    图  11  光电镊驱动微齿轮组实现物体的3维运输[21]

  • [1] LI Jinxing, DE ÁVILA B E F, GAO Wei, et al. Micro/nanorobots for biomedicine: Delivery, surgery, sensing, and detoxification[J]. Science Robotics, 2017, 2(4): eaam6431. doi: 10.1126/scirobotics.aam6431
    [2] ZHANG Yong, CHEN B K, LIU Xinyu, et al. Autonomous robotic pick-and-place of microobjects[J]. IEEE Transactions on Robotics, 2009, 26(1): 200–207. doi: 10.1109/TRO.2009.2034831
    [3] ZHU Wei, LI Jinxing, LEONG Y J, et al. 3D-printed artificial microfish[J]. Advanced Materials, 2015, 27(30): 4411–4417. doi: 10.1002/adma.201501372
    [4] XU Xiaobin, LIU Chao, KIM K, et al. Electric-driven rotation of silicon nanowires and silicon nanowire motors[J]. Advanced Functional Materials, 2014, 24(30): 4843–4850. doi: 10.1002/adfm.201303505
    [5] XIE Hui, SUN Mengmeng, FAN Xinjian, et al. Reconfigurable magnetic microrobot swarm: Multimode transformation, locomotion, and manipulation[J]. Science Robotics, 2019, 4(28): eaav8006. doi: 10.1126/scirobotics.aav8006
    [6] WANG Wei, LI Sixing, MAIR L, et al. Acoustic propulsion of nanorod motors inside living cells[J]. Angewandte Chemie International Edition, 2014, 53(12): 3201–3204. doi: 10.1002/anie.201309629
    [7] PALIMA D and GLÜCKSTAD J. Gearing up for optical microrobotics: Micromanipulation and actuation of synthetic microstructures by optical forces[J]. Laser & Photonics Reviews, 2013, 7(4): 478–494. doi: 10.1002/lpor.201200030
    [8] ASHKIN A and DZIEDZIC J M. Optical trapping and manipulation of viruses and bacteria[J]. Science, 1987, 235(4795): 1517–1520. doi: 10.1126/science.3547653
    [9] WU M C. Optoelectronic tweezers[J]. Nature Photonics, 2011, 5(6): 322–324. doi: 10.1038/nphoton.2011.98
    [10] VALLEY J K, OHTA A T, HSAN-YIN H, et al. Optoelectronic tweezers as a tool for parallel single-cell manipulation and stimulation[J]. IEEE Transactions on Biomedical Circuits and Systems, 2009, 3(6): 424–431. doi: 10.1109/TBCAS.2009.2031329
    [11] LIANG Shuzhang, GAN Chunyuan, DAI Yuguo, et al. Interaction between positive and negative dielectric microparticles/microorganism in optoelectronic tweezers[J]. Lab on A Chip, 2021, 21(22): 4379–4389. doi: 10.1039/D1LC00610J
    [12] VALLEY J K, NEALE S, HSU H Y, et al. Parallel single-cell light-induced electroporation and dielectrophoretic manipulation[J]. Lab on A Chip, 2009, 9(12): 1714–1720. doi: 10.1039/b821678a
    [13] HSIAO Y C, WANG C H, LEE W B, et al. Automatic cell fusion via optically-induced dielectrophoresis and optically-induced locally-enhanced electric field on a microfluidic chip[J]. Biomicrofluidics, 2018, 12(3): 034108. doi: 10.1063/1.5028158
    [14] ZHAO Yuliang, LIANG Wenfeng, ZHANG Guanglie, et al. Distinguishing cells by their first-order transient motion response under an optically induced dielectrophoretic force field[J]. Applied Physics Letters, 2013, 103(18): 183702. doi: 10.1063/1.4827300
    [15] CHEN Y S, LAI C P K, CHEN C, et al. Isolation and recovery of extracellular vesicles using optically-induced dielectrophoresis on an integrated microfluidic platform[J]. Lab on A Chip, 2021, 21(8): 1475–1483. doi: 10.1039/D1LC00093D
    [16] PEI Shaoning, VALLEY J K, WANG Yilun, et al. Distributed circuit model for multi-color light-actuated opto-electrowetting microfluidic device[J]. Journal of Lightwave Technology, 2015, 33(16): 3486–3493. doi: 10.1109/JLT.2015.2405076
    [17] Berkeley Lights[EB/OL]. https://www.berkeleylights.com/systems/beacon/, 2022.
    [18] CHO H, GONZALES-WARTZ K K, HUANG Deli, et al. Bispecific antibodies targeting distinct regions of the spike protein potently neutralize SARS-CoV-2 variants of concern[J]. Science Translational Medicine, 2021, 13(616): eabj5413. doi: 10.1126/scitranslmed.abj5413
    [19] ZHANG Shuailong, LI Weizhen, ELSAYED M, et al. Integrated assembly and photopreservation of topographical micropatterns[J]. Small, 2021, 17(37): 2103702. doi: 10.1002/smll.202103702
    [20] ZHANG Shuailong, SCOTT E Y, SINGH J, et al. The optoelectronic microrobot: A versatile toolbox for micromanipulation[J]. Proceedings of the National Academy of Sciences of the United States of America, 2019, 116(30): 14823–14828. doi: 10.1073/pnas.1903406116
    [21] ZHANG Shuailong, ELSAYED M, PENG Ran, et al. Reconfigurable multi-component micromachines driven by optoelectronic tweezers[J]. Nature Communications, 2021, 12: 5349. doi: 10.1038/s41467-021-25582-8
  • 期刊类型引用(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)

  • 加载中
图(11)
计量
  • 文章访问数:  1012
  • HTML全文浏览量:  1120
  • PDF下载量:  190
  • 被引次数: 28
出版历程
  • 收稿日期:  2022-10-19
  • 修回日期:  2023-01-16
  • 网络出版日期:  2023-02-22
  • 刊出日期:  2023-12-26

目录

/

返回文章
返回