Loading [MathJax]/jax/output/SVG/jax.js
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

卫星通信中高鲁棒性自适应波束赋形算法

张颉 徐厚东 王海 贺洪星 李里 傅宁 涂斌

付善强, 吴冬华, 韩伟涛, 周颖. 基于非线性材料的高速磁浮电磁铁建模与分析[J]. 西南交通大学学报, 2023, 58(4): 879-885. doi: 10.3969/j.issn.0258-2724.20220741
引用本文: 张颉, 徐厚东, 王海, 贺洪星, 李里, 傅宁, 涂斌. 卫星通信中高鲁棒性自适应波束赋形算法[J]. 西南交通大学学报, 2024, 59(3): 556-563. doi: 10.3969/j.issn.0258-2724.20210744
FU Shanqiang, WU Donghua, HAN Weitao, ZHOU Ying. Modeling and Analysis of High-Speed Maglev Electromagnets Based on Nonlinear Materials[J]. Journal of Southwest Jiaotong University, 2023, 58(4): 879-885. doi: 10.3969/j.issn.0258-2724.20220741
Citation: ZHANG Jie, XU Houdong, WANG Hai, HE Hongxing, LI li, FU Ning, TU Bin. Highly Robust Adaptive Beamforming Algorithm in Satellite Communications[J]. Journal of Southwest Jiaotong University, 2024, 59(3): 556-563. doi: 10.3969/j.issn.0258-2724.20210744

卫星通信中高鲁棒性自适应波束赋形算法

doi: 10.3969/j.issn.0258-2724.20210744
基金项目: 国网四川省电力公司科技研发计划(52199720002G)
详细信息
    作者简介:

    张颉(1983—),男,高级工程师,博士,研究方向为北斗电力通信增强技术,E-mail:zhangjie2093@163.com

    通讯作者:

    李里(1983—),男,副教授,博士,研究方向为多天线赋形技术,E-mail:ll5e08@swjtu.edu.cn

  • 中图分类号: TN911.7

Highly Robust Adaptive Beamforming Algorithm in Satellite Communications

  • 摘要:

    在卫星通信中,采用自适应波束赋形技术能极大提升导航接收机的抗干扰性能,然而,传统自适应波束赋形算法对模型失配极为敏感. 针对传统自适应波束赋形算法在期望信号导向矢量和协方差矩阵失配时,输出性能急剧下降的问题,本文提出了一种高鲁棒性自适应波束赋形算法. 该算法首先估计期望信号的输入信噪比(SNR),再根据SNR构造投影矩阵来估计期望信号导向矢量和干扰信号导向矢量;接着,基于不确定集优化方法再次校正期望信号导向矢量和干扰信号导向矢量;最后,估计信号功率并重构干扰加噪声协方差矩阵. 仿真结果表明,相比于协方差矩阵重构波束赋形算法,本文算法具有更优异的输出信干噪比(SINR)性能,当SNR为10.0 dB并存在期望信号波达方向(DOA)估计失配时,本文算法的输出SINR增益可达1.9 dB;当存在波前扰动失真时,输出SINR增益约1.5 dB;当存在局部相干散射时,输出SINR增益约1.6 dB.

     

  • 磁悬浮技术具有无摩擦、微振动、长寿命及高精度等优点,被广泛地应用于高速机械领域. 磁轴承作为核心部件,一般与高速旋转电机配合用于高速旋转机械,如磁悬浮飞轮、磁悬浮控制力矩陀螺及磁悬浮电机等[1-3]. 对于直线运动,通常采用能够实现悬浮及导向功能的电磁铁与直线电机的组合,如磁浮列车、磁浮电梯等. 目前为止,世界很多国家已经对磁浮列车技术研究了较长时间,而且部分实现了商业运行,如中国、德国、日本、韩国等[4-6]. 磁浮列车由电磁铁与轨道功能件之间产生的电磁力支撑,从而实现无接触运行. 电磁铁为悬浮系统的核心执行部件,其电磁力特性决定了列车的承载能力,并与控制器、传感器相互配合实现列车的稳定悬浮. 尤其对于高速运行的磁浮列车,运行载荷更加复杂、苛刻,对承载能力及稳定性要求更高. 因此,电磁力特性分析作为基础研究,对悬浮系统的设计及优化起至关重要作用[7-9].

    目前,电磁力建模分析方法主要包含等效磁路法(EMC)及有限元法(FEM),其中,FEM计算精度较高,但效率较低,很难与控制模型联合用于分析系统实时特性. EMC计算速度较高,能够与控制模型联合用于系统实时动态特性分析,但其计算精度较低. 因此,在传统EMC模型中通常会考虑加入补偿系数,通过调整系数校正电磁力结果,使其与FEM结果接近[10-11]. 然而,传统EMC模型仅考虑线性工作区,导磁材料采用恒定的相对磁导率,忽略磁饱和影响,甚至忽略导磁材料磁阻. 这样会导致EMC模型结果在小电流区间内较准确,而在大电流区间就会出现较大偏差[12-13]. 本文在搭建高速磁浮悬浮电磁铁EMC模型时,考虑了导磁材料的磁阻及其非线性. 通过导磁材料B-HB为磁感应强度;H为磁场强度)曲线的拟合及引入,求解电磁力的准确性大幅度提高,适用范围增加.

    高速磁浮列车的悬浮电磁铁共有12个磁极,极性为NS交替,相邻磁极之间通过磁轭连接,磁场经过长定子铁芯形成回路,磁极与长定子之间的磁场产生电磁吸力,实现悬浮功能. 12个磁极分为左、右两组,分别由两个悬浮控制器单独控制,从而形成两个控制回路,每个回路对应两个间隙传感器. 传感器实时监测电磁铁与长定子之间的间隙,并反馈给悬浮控制器,经过控制策略计算,悬浮控制器输出相应电压给悬浮电磁铁,实现动态稳定悬浮.

    悬浮电磁铁与长定子的物理模型如图1所示. 其中电磁铁分为左、右两个回路单独控制,为简化计算模型,仅对半个电磁铁进行建模,并忽略两个回路间磁场的影响. 搭建的等效磁路模型如图2所示,磁路中包含了气隙磁阻、漏磁磁阻及导磁材料磁阻. 图中:Raj为磁极与长定子间的气隙磁阻;Rsi为长定子铁芯磁阻;Rei为磁极铁芯及磁轭磁阻;Rli为相邻磁极之间的漏磁磁阻;ϕajϕsiϕeiϕli分别为磁阻RajRsiReiRli对应的磁通;ϕpj 为磁极磁通;θj为磁极磁动势nIn为磁极匝数,I为控制回路电流;i=1,2,…,5,j=1,2,…,6.

    图  1  悬浮电磁铁及长定子模型
    Figure  1.  Model of maglev electromagnet and long stator
    图  2  半悬浮电磁铁等效磁路
    Figure  2.  EMC of half-maglev magnet

    根据建模需求,定义磁通、磁动势向量为

    {ϕ=(ϕ1,ϕl1,ϕ2,ϕl2,,ϕ5,ϕl5)T,ϕs=(ϕ1,ϕ2,,ϕ5)T,ϕa=(ϕ1,ϕ1+ϕ2,,ϕ4+ϕ5,ϕ5)T,ϕe=(ϕ1+ϕl1,ϕ2+ϕl2,,ϕ5+ϕl5)T,ϕp=(ϕ1+ϕl1,ϕ1+ϕl1+ϕ2+ϕl2,,ϕ4+ϕl4+ϕ5+ϕl5,ϕ5+ϕl5)T,θ=(θ1,θ2,,θ6)T. (1)

    磁极磁通向量ϕp与磁通向量ϕ的转换关系为

    ϕp =Tϕ

    式中:

    T=[11111111111111111111].

    在等效磁路中,根据基尔霍夫电压定律,建立关于磁通的方程组,如式(3).

    {Ra(j1)ϕa(j1)+Rsiϕsi+Rajϕaj+Reiϕei=θj1+θj,Rliϕli+Reiϕei=θj1+θj, (3)

    式中:j=2,3,,6.

    将式(1)表达为矩阵及向量形式为

    Aϕ=TTθ=TTnI (4)

    式中:向量n=(n, n, n, n, n, nTA为磁阻矩阵,AR10×10

    A=[Ra1+Ra2+Rs1+Re1Re1Ra200Re1Rl1+Re1000Ra20Ra2+Ra3+Rs2+Re2Re2Ra300Re2Rl2+Re20].

    1) 磁极与长定子间气隙磁阻

    受磁极直线发电机(linear generator, LIG)槽与长定子齿槽结构的影响,磁极与长定子之间的气隙磁通分布较为复杂,如图3所示. 因此,将气隙磁通等效分为主磁通、槽磁通及LIG磁通,分别对应3种磁阻.

    图  3  气隙磁场分布
    Figure  3.  Magnetic field distribution of air gap

    气隙磁阻为3种磁阻并联,即

    1Raj=1Ra,mj+1Ra,nj+1Ra,Lj, (5)

    式中:Ra,mj为主磁通对应的磁阻;Ra,nj为槽磁通对应的磁阻;Ra,Lj为LIG磁通对应的磁阻.

    每种磁阻可由式(6)计算.

    Ra,oj=sj+hojμ0Aa,oj,o[m,n,L] (6)

    式中:sj为磁极与长定子齿之间的间隙;μ0为空气磁导率;Aa,oj为相应气隙面积;hoj为额外气隙长度;m、n、L分别对应主磁通、槽磁通和LIG磁通.

    2) 相邻磁极间漏磁磁阻

    磁极磁通大部分经过长定子回到相邻磁极,小部分未经过长定子而直接回到相邻磁极,该部分磁通为相邻磁极之间的漏磁,对应的磁阻称为漏磁磁阻,可由式(7)计算.

    Rli=hliμ0Ali (7)

    式中:Ali为等效气隙面积;hli为相邻磁极间的等效气隙长度.

    3) 导磁材料磁阻

    对于长定子铁芯、磁极铁芯及磁轭的磁阻,采用分段方式进行求解,且尽可能保证每段的截面积相同. 具体分段如图4所示,其中,磁极及磁轭共分为5段(1~5),长定子分为3段(6~8). 该部分磁阻计算时考虑导磁材料的非线性.

    图  4  铁芯的分段
    Figure  4.  Iron core sections

    磁极铁芯及磁轭磁阻为

    Rei=5k=1lei,kμ0μr(ϕei,k)Aei,k (8)

    长定子铁芯磁阻为

    Rsi=8k=6lsi,kμ0μr(ϕsi,k)Asi,k, (9)

    式(8)、(9)中:lei,k, lsi,k分别为悬浮电磁铁侧、长定子侧每段铁芯的长度;Aei,k, Asi,k分别为悬浮电磁铁侧、长定子侧每段铁芯的截面积;μr为每段铁芯的相对磁导率;ϕsi,k为第i个磁回路中第k段长定子内的磁通量;ϕei,k为第i个回路中第k段磁极铁芯或磁轭内的磁通量.

    为在EMC模型中引入导磁材料的非线性,不再将μr简单地设置为恒定值,而是根据每段铁芯的磁通进行计算. 磁极铁芯、磁轭及长定子铁芯均采用硅钢片,牌号为M530-50A,导磁材料的非线性可通过B-H曲线体现,如图5所示.

    图  5  铁芯B-H曲线- M530-50A
    Figure  5.  B-H curve of iron core-M530-50A

    根据文献[14],B-H曲线可采用式(10)函数进行拟合.

    H=α1B+α2Bα3, (10)

    式中:α1α2α3为拟合函数自变量B的系数,可通过对图5B-H曲线的拟合确定.

    长定子铁芯的相对磁导率可表示为

    μr1=Bμ0H=1μ0(α1+α2Bα31). (11)

    根据磁密、面积及磁通的关系,式(11)可表达为

    μr1=1μ0(α1+α2(ϕsiAsi,k)α31). (12)

    同理可求解磁极铁芯及磁轭的相对磁导率为

    μr2=1μ0(α1+α2(ϕeiAei,k)α31). (13)

    将式(12)、(13)代入式(8)、(9)求解磁极铁芯、磁轭及长定子铁芯磁阻. 通过式(5)~(9)可知,磁阻矩阵A与间隙s及磁通ϕ相关,因此,将其记为As,ϕ).

    悬浮电磁铁分为两个控制回路,分别由一个控制器进行供电. 控制器的输出为电压,根据控制回路的负载特性转变为相应的负载电流. 在进行电磁力模型与控制模型联合分析时,需要搭建电磁铁的电路模型,其功能是将控制模型的输入电压转变为负载电流,再结合EMC模型计算电磁力. 与EMC模型类似,仅搭建一个回路的模型,6个磁极串联实际可等效为6个电阻与6个电感的串联,如图6所示. 图中:RMj为单个磁极电阻;LMj为单个磁极电感;ULj为单个磁极电感电压;UM为单个控制回路的输入电压.

    图  6  悬浮电磁铁控制回路
    Figure  6.  Control loop of maglev electromagnet

    根据电路模型,单个控制回路的输入电压为

    UM=RMI+6j=1ULj=RMI+nTϕp, (14)

    式中:RM为单个控制回路的总电阻.

    整理式(14),控制回路电流可表达为

    I=(UMnTϕp)/RM. (15)

    式(4)给出了磁通与电流关系,式(15)给出了磁通、电流与电压关系. 将式(15)代入式(4),可得磁通与电压的关系为

    A(s,ϕ)ϕ=TTn(UMnTϕp)/RM. (16)

    根据磁通关系ϕp=Tϕ,式(16)整理为

    TTnnTTϕ = TTnUMRMA(s,ϕ)ϕ. (17)

    经过计算发现,矩阵TTnnTT为奇异矩阵,常微分式(17)很难进行求解. 因此,为便于求解方程,采用中间变量替换原变量.

    根据磁通关系ϕp=Tϕ,将式(16)整理成变量为ϕp的方程,如式(18)

    A(s,ϕ)T1ϕp=TTn(UMnTϕp)/RM. (18)

    进一步整理为

    nTϕp=nTTA1(s,ϕ)TTn(UMnTϕp)/RM. (19)

    将6个磁极的磁通之和β作为中间变量,则β可表示为

    β=6j=1ϕpj. (20)

    将式(20)代入式(19),得到关于β的常微分方程为

    nβ=nTTA1(s,ϕ)TTn(UMnβ)/RM, (21)

    式中:nTTA−1(s,ϕ)TTn不再是一个矩阵或向量,而是一个关于sϕ的变量,记为1/M (s,ϕ).

    对式(21)进行整理得

    ˙β=UM/nRMM(s,ϕ). (22)

    依据式(15)、(22),可得Iβ间的关系为

    I=nM(s,ϕ)β. (23)

    依据式(19)、(22)及磁通关系ϕp=Tϕ,得到ϕβ间的关系为

    ϕ=A1(s,ϕ)TTnM(s,ϕ)nβ. (24)

    通过式(22)~(24)计算出悬浮电磁铁电流以及电磁力求解所需的磁通.

    根据式(5),计算各部分气隙的磁通为

    ϕa,oj=ϕajRajRa,oj. (25)

    基于虚功原理,各部分气隙对应的电磁力为

    Fmag,oj=ϕ2a,oj2μ0Aa,oj. (26)

    单个磁极电磁力为式(26)中各部分电磁力之和,即

    Fmag,j=Fmag,mj+Fmag,nj+Fmag,Lj.

    最后,求解半个悬浮电磁铁的电磁力为

    Fmag=6j=1Fmag,j.

    根据电磁力解析过程,EMC模型可简化为图7所示的结构框图,输入为电压及间隙,输出为电流及电磁力. 首先进行磁阻计算,并组建磁阻矩阵,而导磁材料磁阻的计算需要将磁通作为输入. 采用磁阻矩阵A进行电流及磁通计算,磁通需通过常微分方程及代数方程求解,磁通求解结果一方面用于电磁力计算,一方面反馈给导磁材料磁阻计算.

    图  7  电磁力模型框图
    Figure  7.  Magnetic force model

    根据悬浮电磁铁及长定子的尺寸及参数(如表1所示),对本文EMC模型进行量化. 此外,将导磁材料的相对磁导率设为恒定值,搭建基于线性导磁材料的传统EMC模型. 对两个模型电磁力的计算结果进行对比分析,如图8所示.

    表  1  悬浮电磁铁及长定子参数
    Table  1.  Parameters of maglev electromagnet and long stator
    项点取值项点取值
    定子极距/mm258.0铁芯厚度/mm170.0
    电磁铁极距/mm266.5磁极匝数300
    定子齿宽度/mm43.0额定磁间隙/mm12.5
    定子槽宽度/mm43.0恒定相对磁导率7 000
    下载: 导出CSV 
    | 显示表格
    图  8  线性与非线性材料EMC电磁力
    Figure  8.  Electromagnetic forces of EMC models with linear and nonlinear materials

    磁间隙为12.5 mm,电流在0~35 A内,两个模型的电磁力结果非常接近;当电流超过35 A时,随着电流增加,计算结果偏差增大. 原因是实际工作状态下,随着电流增大,导磁部件的磁密增大;当达到材料饱和磁密时,磁密随电流的增加率大幅降低,电磁力也相应地出现饱和现象,而传统EMC模型并未考虑材料的磁饱和. 因此,传统EMC模型适用于小电流区间,一般应用在工作点附近的线性区间20~30 A. 而大电流区间与实际工作情况不符,例如在故障、起浮、降落等特殊工况时,模型精度大幅降低,无法用于电磁力计算及系统特性分析.

    为验证本文提出的EMC模型准确性,搭建了悬浮电磁铁与长定子FEM模型,如图9所示,两者电磁力的计算结果如图10所示. 磁间隙12.5 mm与16.0 mm,电流0~80 A内,两者电磁力的计算结果具备非常高的一致性,均存在饱和现象;磁间隙12.5 mm,电流50 A时,电磁力偏差最大,EMC计算结果为115 kN,FEM计算结果为110 kN,偏差仅为4.5%,这表明了基于非线性材料的EMC模型具有较高的准确性.

    图  9  悬浮电磁铁及长定子有限元模型
    Figure  9.  FEM model of maglev electromagnet and long stator
    图  10  EMC与FEM电磁力结果
    Figure  10.  Electromagnetic force results of EMC and FEM

    针对高速磁浮电磁铁特性研究,搭建了地面试验平台,对悬浮电磁铁的静态电磁力进行测试,如图11所示. 试验台通过液压系统调整长定子与悬浮电磁铁的间隙;采用两路电源供电,但受最大输出电流限制,仅对0~50 A电流进行测试,步长为5 A;通过力传感器检测电磁铁与长定子间的电磁力. 磁间隙12.5 mm下电磁力的测试结果与EMC及FEM的计算结果对比如图12所示.

    图  11  悬浮电磁铁静态电磁力测试
    Figure  11.  Static electromagnetic force test of maglev electromagnet
    图  12  电磁力计算及测试结果
    Figure  12.  Electromagnetic force calculation and test results

    额定工况下,悬浮电磁铁的工作点:磁间隙为12.5 mm,电流为25 A,电磁力约为46 kN. 在工作点处,EMC、FEM及试验测试的电磁力结果几乎相同,在其他电流值下,电磁力结果偏差也极小,从而进一步验证了本文EMC模型以及所搭建FEM模型的准确性.

    本文基于非线性材料搭建了高速磁浮悬浮电磁铁的磁路、电路及磁力模型,将计算结果与传统EMC模型进行对比分析,并通过有限元及试验验证,通过对模型研究分析,得到以下结论:

    1) 搭建悬浮电磁铁EMC模型时,采用了非线性导磁材料,通过引入B-H曲线的拟合函数,将导磁材料的非线性及饱和特性体现在模型中.

    2) 无论小电流区,还是大电流区,本文EMC模型求解的电磁力均与实际情况接近,相比传统EMC模型,结果更加准确,适用范围更广.

    3) 本文EMC模型能够快速、准确地求解电磁力,且通过电路模型实现与控制模型的良好匹配,因此,可通过联合仿真对悬浮系统动态特性进行深入分析,为悬浮系统设计及参数优化提供了依据.

  • 图 1  均匀直线阵自适应波束形成器结构

    Figure 1.  Adaptive beamformer structure of uniform linear array

    图 2  存在信号DOA估计误差时各算法的输出性能

    Figure 2.  Output performances of algorithms with signal DOA mismatch

    图 3  存在波前扰动失真时各算法输出性能

    Figure 3.  Output performances of algorithms with wavefront distortion

    图 4  存在局部相干散射时各算法输出性能

    Figure 4.  Output performances of algorithms with coherent local scattering

    图 5  各算法计算复杂度随天线阵元数变化曲线

    Figure 5.  Computation complexity curve varying with number of antenna elements

  • [1] 卢丹. 稳健的全球卫星导航系统抗干扰技术研究[D]. 西安: 西安电子科技大学,2013.
    [2] 朱亮. 北斗卫星导航系统干扰识别与测向技术的研究与实现[D]. 北京: 北京交通大学,2019.
    [3] 刘江波. 稳健的收/发波束形成方法研究[D]. 成都: 电子科技大学,2018.
    [4] 黄磊. 非理想条件下的自适应波束形成算法研究[D]. 合肥: 中国科学技术大学,2016.
    [5] COX H, ZESKIND R, OWEN M. Robust adaptive beamforming[J]. IEEE Transactions on Acoustics, Speech, and Signal Processing, 1987, 35(10): 1365-1376. doi: 10.1109/TASSP.1987.1165054
    [6] FELDMAN D D, GRIFFITHS L J. A projection approach for robust adaptive beamforming[J]. IEEE Transactions on Signal Processing, 1994, 42(4): 867-876. doi: 10.1109/78.285650
    [7] LI J, STOICA P, WANG Z. Doubly constrained robust Capon beamformer[C]//The Thrity-Seventh Asilomar Conference on Signals, Systems & Computers. Pacific Grove: IEEE, 2004: 1335-1339.
    [8] JIAN L, STOICA P, WANG Z S. On robust Capon beamforming and diagonal loading[J]. IEEE Transactions on Signal Processing, 2003, 51(7): 1702-1715. doi: 10.1109/TSP.2003.812831
    [9] GU Y J, LESHEM A. Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation[J]. IEEE Transactions on Signal Processing, 2012, 60(7): 3881-3885. doi: 10.1109/TSP.2012.2194289
    [10] HUANG L, ZHANG J, XU X, et al. Robust adaptive beamforming with a novel interference-plus-noise covariance matrix reconstruction method[J]. IEEE Transactions on Signal Processing, 2015, 63(7): 1643-1650. doi: 10.1109/TSP.2015.2396002
    [11] HUANG Y W, ZHOU M K, VOROBYOV S A. New designs on MVDR robust adaptive beamforming based on optimal steering vector estimation[J]. IEEE Transactions on Signal Processing, 2019, 67(14): 3624-3638. doi: 10.1109/TSP.2019.2918997
    [12] 毛晓军. 高性能阵列天线稳健自适应波束形成技术研究[D]. 哈尔滨: 哈尔滨工程大学,2017.
    [13] SHEN F, CHEN F F, SONG J Y. Robust adaptive beamforming based on steering vector estimation and covariance matrix reconstruction[J]. IEEE Communications Letters, 2015, 19(9): 1636-1639. doi: 10.1109/LCOMM.2015.2455503
    [14] 毛英. 北斗空时自适应抗干扰子空间投影方法研究[D]. 哈尔滨: 哈尔滨工程大学,2017.
    [15] ZHENG Z, ZHENG Y, WANG W Q, et al. Covariance matrix reconstruction with interference steering vector and power estimation for robust adaptive beamforming[J]. IEEE Transactions on Vehicular Technology, 2018, 67(9): 8495-8503. doi: 10.1109/TVT.2018.2849646
    [16] ZHOU C W, GU Y J, HE S B, et al. A robust and efficient algorithm for coprime array adaptive beamforming[J]. IEEE Transactions on Vehicular Technology, 2018, 67(2): 1099-1112. doi: 10.1109/TVT.2017.2704610
    [17] 鲁郁. 北斗/GPS双模软件接收机原理与实现技术[M]. 北京: 电子工业出版社,2016:82-85.
    [18] KHABBAZIBASMENJ A, VOROBYOV S A, HASSANIEN A. Robust adaptive beamforming based on steering vector estimation with as little as possible prior information[J]. IEEE Transactions on Signal Processing, 2012, 60(6): 2974-2987. doi: 10.1109/TSP.2012.2189389
  • 加载中
图(5)
计量
  • 文章访问数:  526
  • HTML全文浏览量:  144
  • PDF下载量:  32
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-09-16
  • 修回日期:  2022-01-04
  • 网络出版日期:  2024-04-17
  • 刊出日期:  2022-04-01

目录

/

返回文章
返回