Research on Flow Scheduling of Wired and Wireless Converged Satellite Time Sensitive Network
-
摘要: 随着空间通信任务日趋复杂化，尤其是对时间敏感的需求不断提升，一方面要求星内系统的高带宽、可靠性和实时性；另一方面星间无线链路也应具备低时延和高可靠性。但由于卫星内部有线链路与星间无线链路差异大，业务数据经过有线和无线链路联合传输时，容易引发节点拥塞，而无法保障时敏业务的时延有界需求。为了提升数据在空间网络传输的实时性，该文提出了一种有线无线融合的时间敏感网络(TSN)流调度方案，首先对有线和无线链路资源分配与终端时延关系进行分析建模，并通过TSN控制器收集终端时敏需求，构建以全网时敏业务端到端最小平均时延为优化目标，然后采用基于增强精英保留遗传算法进行方案的快速求解。通过Pycharm对比测试时隙分配算法的性能，同时设计实现基于EXata网络仿真平台的低轨卫星TSN系统，并搭建实验场景进行试验验证。测试结果表明，该文所提出的流联合调度方案能够为空间时敏任务提供有界、稳定的时延保障。Abstract: With the increasing complexity of space communication tasks, especially the increasing demand for time sensitivity, on the one hand, the performances of high bandwidth, high reliability and real time in intra satellite systems are required; on the other hand, low latency and high reliability are demanded for inter satellite wireless links. However, due to the large difference between the satellite internal wired link and the inter-satellite wireless link, when the data is transmitted from the wired link to wireless link, it is easy to cause congestion problem on node, and can not guarantee the bounded latency for time-sensitive service. To improve the real time performance for data transmission in space information network, a flow scheduling scheme in wired and wireless converged Time-Sensitive Network (TSN) is proposed. Firstly, the relationship between terminal delay requirement and wired/wireless link resources is analyzed. Then, the time-sensitive requirements of terminals are collected by TSN controller and the optimization goal is determined by the end-to-end minimum average delay of time-sensitive flow in the entire network. Moreover, the enhanced elite retention genetic algorithm is adopted to solve quickly the scheme. Finally, the performance of proposed time slot allocation algorithm is verified through Pycharm. Meanwhile, a low-orbit satellite network scenario is implemented under EXata network simulation platform to evaluate the proposed scheme in the further. The results demonstrate that the proposed wired and wireless joint stream scheduling scheme can provide bounded and stable delay for space time-sensitive tasks.
-
1. 引言
作为2维平面结构的超材料，超表面[1,2]通过精心设计亚波长尺寸的单元结构可以调控电磁波的幅度、相位和极化方式，具有低剖面、重量轻、易于制作等优点，被广泛应用于极化转化器[3,4]、涡旋发生器[5,6]、图像处理器[7]、天线[8,9]等领域。而随着现代集成系统的迅速发展，功能单一的超表面已无法满足集成系统的需求。相对于有源超表面(如可重构超表面)[10,11]为实现多功能进行的复杂控制系统设计，无源超表面通过改变入射波的极化、频率和角度等实现多种功能的集成，极大降低了系统的损耗和成本。2018年，文献[12]基于P-B相位在微波波段提出了一款3层金属结构堆叠的双功能超表面；同年，Zhuang等人[13]利用一种基于反射和透射结构的多层超表面实现3种不同的功能；2019年，文献[14]研制了一款集成共振和几何相位的级联多功能超表面，以独立操控反射和透射电磁波；2021年，Shang等人[15]利用一组双层交叉型反射单元设计了一款2 bit编码超表面，在不同极化波作用下实现两幅独立的全息图像。上述研制的多功能超表面，大多数采用多层结构堆叠的形式实现360°相位覆盖，并不是研制低成本和集成化功能性器件的最佳解决方案。
惠更斯超表面[16-19]通过合理构建电偶极子和磁偶极子，能够对电场和磁场产生响应，从而实现对入射电磁波的自由调控，相对于只调控极化电流的超表面，它既能保证高的透射效率，又能保证全相位覆盖，且制作工艺简单，能显著提高超表面效率，一直是人工电磁领域的研究热点。文献[20]提出了一款极化不敏感的惠更斯超表面，在远场产生高效率的全息图像。2018年，Wang等人[21]在微波波段研制了4款惠更斯超表面全息图，并在y极化电磁波的激励下呈现不同的图像。之后，该研究团队设计了一款宽频段的惠更斯超表面，y极化电磁波经过超表面的调制实现3维全息图像[22]。目前大多数研究的惠更斯超表面极化模式单一，并且需要在物理结构上同时构建电、磁偶极子，设计过程相对复杂，极大限制了在实际工程中的应用。因此，围绕惠更斯超表面开展结构简单的新型极化复用电磁集成器件的研究方兴未艾。
本文提出一种极化复用的聚焦惠更斯超表面，可通过改变相应的物理参数实现x极化和y极化独立调控，35 GHz频率下具有较小的传输损耗和较高的聚焦质量。超表面单元由一层厚度为0.17λ的介质基板和位于两侧的两个不对称电偶极子元件组成，物理结构上没有单独的磁偶极子元件，而是通过反向流动的表面电流诱导磁流，使单元结构更加紧凑。利用全息理论，设计了在x极化波和y极化波作用下具有独立聚焦特性的极化复用超表面，仿真和实测结果基本一致，验证了所提出的惠更斯超表面具有极化复用的可行性。
2. 超表面单元设计
惠更斯超表面来源于等效原理，该原理通过无穷小的薄型表面电流密度Js和磁流密度Ms来满足边界条件[16,17]。当入射波激励惠更斯表面时，会同时感应表面电流和磁流密度，产生依赖于电流和磁流密度的反射波和透射波，如图1所示。对于区域1中的入射场，区域2中的任意场可以通过在超表面上引入电流和磁流密度来确定。电磁场和表面电磁流密度满足边界条件为
Js=n×(H2−H1) (1) Ms=−n×(E2−E1) (2) 其中，n为朝向区域2的单位法向向量，E和H为电场和磁场。在实际应用中，通过取表面电磁流密度与表面平均切向场的比值，得到特定表面电磁流密度与表面电导纳
ˉˉYes 和磁阻抗ˉˉZms 之间的关系Js=ˉˉYes⋅Etan,av|S (3) Ms=ˉˉZms⋅Htan,av|S (4) 其中，Etan,av和Htan,av分别是表面上的平均切向电场和磁场。
ˉˉYes 和ˉˉZms 为张量，表示为ˉˉYes = (YxxesYxyesYyxesYyyes) (5) ˉˉZms = (ZxxmsZxymsZyxmsZyyms) (6) 由场分布得到表面阻抗后，将超表面离散为单元。根据反射系数(R)和透射系数(T)确定特定单元的表面阻抗
Yes=2η⋅1−T−R1+T+R (7) Zms=2η⋅1−T+R1+T−R (8) 其中，η为自由空间的波阻抗，对于无反射、全透射的惠更斯超表面，即
R=0 ,T=ejφt 满足Yes=−j2η⋅tan(φt2) (9) Zms=−j2η⋅tan(φt2) (10) 由式(9)和式(10)可以看出，表面阻抗与传输相位
φt 有关，对于每个特定的传输相位，都有相应的表面阻抗。此外，如果归一化的表面电导纳和磁阻抗相等且为纯虚数(Yesη = Zms/η)，则会激发电磁共振，获得高传输幅度。图2说明了本文利用的透射型惠更斯单元结构，分别由x方向和y方向的金属结构线性叠加构成，具有4种不同叠加方式。当y极化波入射时，调控y方向的金属结构可产生轴向表面电流，实现传输幅度接近1的宽带传输谱和360°相位覆盖，而不改变x极化波的传输相位和幅度。同理，x方向的金属结构可对x极化波进行波束调控，而对y极化波不产生响应，单元实现极化复用特性。
图3(a)和图3(b)为对y极化波响应的单元结构E_1和E_2，由印刷在1.5 mm厚的介质基板(介电常数2.22，tanδ = 0.001)两侧的不对称金属贴片(电偶极子)构成，其中P = 4.8 mm, h = 1.5 mm, w = 0.2 mm。通过建立归一化电导纳Yesη与磁阻抗Zms/η同单元传输幅度之间的关系分析其工作原理，如图4(a)所示，谐振点31.5 GHz和34.2 GHz处，Yesη与Zms/η的虚部相等，实部趋于0，实现电磁共振，形成两个传输幅度峰值。两电磁共振频率之间形成一个通频带，传输幅度优于0.9。谐振点34.2 GHz处的表面电流在矩形金属贴片上流动，假设上、下层金属贴片的电流分别为Iup和Ilow，他们具有相同的电流密度但具有不同的相位，即Ilow =Iupexp(jΔ
φ )，Δφ 代表延迟相位。图4(b)给出了随时间相位因子θ = ωt变化的表面电流分布，灰色部分的上、下金属贴片的表面电流方向相反，在单元侧面形成电流环，即构成磁偶极子，诱导出垂直于电流的磁流，与表面电流相互作用激发电磁共振。因此，可通过调控单元E_1的贴片长度L1和单元E_2的贴片弧长L2获得所需的电磁共振，实现360°相位覆盖和高传输幅度。y极化电磁波垂直入射到超表面单元，L1从1 mm变化到4.8 mm，共极化透射波Tyy的传输幅度均优于0.90，传输相位在445°～105°之间变化，实现340°相位积累，L1的变化对交叉极化透射波Txy的传输相位几乎没有影响，可独立实现对x极化波和y极化波的相位调控，如图5(a)所示。由图5(b)可知，L2从0.95 mm变化到1.40 mm，共极化透射波Tyy的传输幅度优于0.90，传输相位在105°～85°之间变化，相位积累达到20°，L2的变化使交叉极化透射波Txy的传输相位偏移7°左右，偏移量较小，可以在后续超表面的阵列排布中进行修正。结合图5(a)和图5(b)的分析，通过改变单元E_1和单元E_2的金属参数L1和L2，在中心频率35 GHz处实现360°的相位覆盖。
3. 超表面设计与测试
本文利用全息理论[18,23]，设计了在x极化和y极化波作用下具有独立聚焦特性的极化复用惠更斯超表面。为实现任意波前的再现，惠更斯超表面需要如下步骤得到相位补偿：首先，将所需要的波前看作物波，并记录相位补偿；其次，利用相同的参考波对全息图进行激励，再现物波。基于此思想，在超表面设计过程中，将期望的目标场分布和入射平面波分别视为物波和参考波。
图6是基于超表面的全息原理图，详细的相位补偿推导如下：超表面上的场分布来自目标场辐射的球面波的叠加。为了再现目标场，需要将目标场离散为多个虚拟点源。假设在Fn位置产生n个点源，则超表面的叠加E场表示为
E(xi,yi)=N∑n=1An(xi,yi)⋅exp(jk0Rn) (11) An(xi,yi) = AnRn (12) Rn=|Fn−ri| (13) (xi,yi)表示第i个超表面单元的中心坐标。常数An是第n个点源的振幅，Rn是第n个点源到第i个超表面单元的距离。k0为自由空间中的波数。Fn是第n个点源到原点的距离，ri是第i个单元到原点的距离，N为离散的点源总数目。参考波为正入射的平面波，到超表面上的振幅和相位恒定。假设相位是
φr ，超表面的相位补偿为φ(xi,yi)=φr−angleE(xi,yi) (14) 为了验证所设计的惠更斯超表面在35 GHz处具有极化复用特性，设置x极化波作用下，惠更斯超表面在焦距为40 mm和80 mm分别实现单焦点和4个焦点的3维聚焦。在y极化波作用下，焦距70 mm平面处再现字母“I”，图7为极化复用超表面的相位分布和尺寸分布。利用商用软件CST Microwave Studio进行模拟，得到在线极化波垂直入射下，经过惠更斯超表面调制的透射场分布。图8(a)为x极化波作用下，焦距为40 mm和80 mm的电场强度分布，实现了单个焦点和4个焦点的3维聚焦。而y极化波作用下，在焦距为70 mm的平面中心再现字母“I”，如图8(b)所示。焦点和图像“I”附近存在少量的衍射场，是电磁波通过惠更斯超表面产生了基于衍射光图案的高次谐波，可采取一些优化算法改善效果。
采用印刷电路板(Printed Circuit Board, PCB)工艺制作了超表面样品，以实验验证模拟的极化复用惠更斯超表面在中心频率处的性能，如图9(b)。设计的惠更斯超表面由729个单元组成，覆盖面积为129.6 mm×129.6 mm。利用喇叭天线发射线极化准平面波(通过调整端口方向以实现不同的极化方式)，并使其与超表面中心保持一致，用近场探针在指定焦距处测量透射电场，如图9(a)所示。图10分别为x极化波和y极化波入射下，频率为35 GHz的电场强度分布，从图中可以看出，x极化波经过惠更斯超表面调制后，在焦距为40 mm和80 mm平面处分别实现单焦点和4个焦点的3维聚焦。而y极化波经过超表面的调制后再现字母“I”。
为进一步验证实测与模拟结果的一致性，提取特定平面处x方向和y方向的电场强度，如图11所示。图11(a)显示利用极化复用惠更斯超表面调制x极化波后，提取焦距为40 mm，中心焦点x方向的模拟和实测电场强度分布，以及焦距为80 mm, 4个焦点x方向和y方向的电场强度分布，结果表明实测的单焦点聚焦和模拟结果一致，而4个焦点的实测结果相对于模拟结果有轻微的偏差。这样的偏差主要来自喇叭天线辐射的入射波并不是严格的平面波，在超表面上存在轻微离轴偏差。此外，相位补偿误差、加工缺陷和测量误差等也会造成焦点之间的偏差和微弱的强度变化。对于y极化，图11(b)展现了焦距为70 mm，x方向和y方向的模拟和实测结果，实测的峰值强度与模拟结果有一定的偏差，存在的主要原因是喇叭天线提供的非理想的平面波前，并且加工缺陷、测试误差也会造成一定影响。实测和模拟结果表明，所设计的惠更斯超表面能够实现对x极化波和y极化波的独立调控，具有极化复用特性。
超表面的聚焦效率η可以由聚焦平面处聚焦区域的能量总和P1与无超表面时该平面处的辐射功率P2的比值来确定，即η=P1/P2，P1和P2可通过相应区域的坡印廷矢量积分得到。对于x极化波，计算得中心频率35 GHz的聚焦效率为29.73%，y极化波的聚焦效率达到33.64%。
表1为本文设计的超表面与已有超表面的对比。通过比较发现，所设计的超表面具有双极化独立调控功能。与文献[25,26]相比，本文设计的惠更斯单元实现360°全相位覆盖和优于0.9的高传输幅度，且单层无金属过孔，具有低剖面和结构简单等特点，可印刷在价格低廉的PCB板两侧，降低加工难度。此外，所设计的惠更斯单元没有单独构建磁偶极子，在物理结构上消除了对磁性元件的需求，使超表面更加紧凑。
表 1 文献对比分析4. 结论
本文采用4种简单的单元编码实现一款极化复用的聚焦惠更斯超表面。通过调整电磁共振，实现了惠更斯超表面的高透射率、360°相位覆盖和极化控制。并在35 GHz处进行了实验验证，通过操纵入射电磁波的极化分量，超表面根据透射波的极化状态实现独立聚焦，实测和仿真结果吻合良好。所提出的超表面具有结构紧凑、极化调控、易加工等优点，可应用于信息复用、多通道图像显示、加密防伪、存储显示等领域。
-
表 1 TSN-GA算法
输入：网络设备数N，每个设备权值αi，无线传输速率v，设备
业务量大小Bsi，种群数量Np，最大迭代数maxCycle，决策变量
维度Dim，变异概率Pm，交叉概率Pc输出：网络最低时延、时隙最佳位置 (1) while Cycle≤maxCycle do (2) for m=1 to Np (3) vnm=Crossover(αnm) (4) end for (5) for m=1 to Np (6) μnm=Mutation(vnm) (7) end for (8) for m=1 to 2Np (9) ρnm=Select(αnm, μnm) (10) end for (11) Cycle=Cycle+1; (12) end while (13) for Cycle=1 to maxCycle (14) min\_delay(Cycle) = F(Cycle)； (15) best\_position(Cycle)= F(Cycle) (16) end for 
下载: 导出CSV
表 2 编队卫星的参数设置
参数 数值 TSN控制卫星轨道参数[a,e,i,Ω,w,f] 8×106 m, 0°, 30°, 15°, 15°, 10° 空间圆半径(R/km) 30 编队卫星个数 7 成员卫星1初始相位角θ1(°) 0 相位之差Δθ(°) 72
下载: 导出CSV
表 3 实验参数设置
业务类型 源端 目的端 周期(ms) 业务量(bit)
快速位置报告信息终端11 终端33 6 500 终端18 终端33 6 500 终端14 终端33 6 500
北斗短报文终端13 终端37 12 2000 终端8 终端15 12 2000 终端23 终端29 12 2000
遥感勘测终端12 终端15 24 4000 终端18 终端17 24 4000 终端29 终端12 24 4000
下载: 导出CSV
-
[1] XIE Renchao, TANG Qinqin, WANG Qiuning, et al. Satellite-terrestrial integrated edge computing networks: Architecture, challenges, and open issues[J]. IEEE Network, 2020, 34(3): 224–231. doi: 10.1109/MNET.011.1900369 [2] 蔚保国, 鲍亚川, 魏海涛. 面向时间同步业务的空间信息网络拓扑聚合图模型[J]. 电子与信息学报, 2017, 39(12): 2929–2936. doi: 10.11999/JEIT170252YU Baoguo, BAO Yachuan, and WEI Haitao. Time synchronization service oriented topology aggregation model of space information network[J]. Journal of Electronics &Information Technology, 2017, 39(12): 2929–2936. doi: 10.11999/JEIT170252 [3] 吕梦昭. 空间时延敏感通信关键技术研究[D]. [硕士论文], 电子科技大学, 2020.LV Mengzhao. Research on key technologies of space time delay sensitive communication[D]. [Master dissertation], University of Electronic Science and Technology of China, 2020. [4] NASRALLAH A, THYAGATURU A S, ALHARBI Z, et al. Ultra-low latency (ULL) networks: The IEEE TSN and IETF DetNet standards and related 5G ULL research[J]. IEEE Communications Surveys & Tutorials, 2019, 21(1): 88–145. doi: 10.1109/COMST.2018.2869350 [5] HUANG Li, LIANG Yun, ZHANG Yajian, et al. Time-sensitive network technology and its application in energy internet[C]. 2019 IEEE International Conference on Energy Internet (ICEI), Nanjing, China, 2019: 211–216. [6] ZHOU Zefeng and SHOU Guochu. An efficient configuration scheme of OPC UA TSN in industrial internet[C]. 2019 Chinese Automation Congress (CAC), Hangzhou, China, 2019: 1548–1551. [7] ZHOU Zifan, LEE J, BERGER M S, et al. Simulating TSN traffic scheduling and shaping for future automotive Ethernet[J]. Journal of Communications and Networks, 2021, 23(1): 53–62. doi: 10.23919/JCN.2021.000001 [8] SEIJO Ó, LÓPEZ-FERNÁNDEZ J A, and VAL I. w-SHARP: Implementation of a high-performance wireless time-sensitive network for low latency and ultra-low cycle time industrial applications[J]. IEEE Transactions on Industrial Informatics, 2021, 17(5): 3651–3662. doi: 10.1109/TII.2020.3007323 [9] BHATTACHARJEE S, SCHMIDT R, KATSALIS K, et al. Time-sensitive networking for 5G fronthaul networks[C]. The ICC 2020 - 2020 IEEE International Conference on Communications (ICC), Dublin, Ireland, 2020: 1–7. [10] WANG Xinheng, XU Chuan, ZHAO Guofeng, et al. Tuna: An efficient and practical scheme for wireless access point in 5G networks virtualization[J]. IEEE Communications Letters, 2018, 22(4): 748–751. doi: 10.1109/LCOMM.2017.2768511 [11] CHAINE P J, BOYER M, PAGETTI C, et al. TSN support for quality of service in space[C]. The 10th European Congress on Embedded Real Time Software and Systems (ERTS 2020), Toulouse, France, 2020: 1–7. [12] SEIJO O, FERNÁNDEZ Z, VAL I, et al. SHARP: Towards the integration of time-sensitive communications in legacy LAN/WLAN[C]. 2018 IEEE Globecom Workshops (GC Wkshps), Abu Dhabi, United Arab Emirates, 2018: 1–7. [13] LARRAÑAGA A, LUCAS-ESTAÑ M C, MARTINEZ I, et al. Analysis of 5G-TSN integration to support industry 4.0[C]. The 25th IEEE International Conference on Emerging Technologies and Factory Automation (ETFA), Vienna, Austria, 2020: 1111–1114. [14] GINTHÖR D, GUILLAUME R, VON HOYNINGEN-HUENE J, et al. End-to-end optimized joint scheduling of converged wireless and wired time-sensitive networks[C]. The 25th IEEE International Conference on Emerging Technologies and Factory Automation (ETFA), Vienna, Austria, 2020: 222–229. [15] 王鹏宇. 飞行器编队网络测距技术研究[J]. 遥测遥控, 2019, 40(2): 22–30. doi: 10.3969/j.issn.2095-1000.2019.02.004WANG Pengyu. Research on ranging technique for spacecraft formation flying network[J]. Journal of Telemetry,Tracking and Command, 2019, 40(2): 22–30. doi: 10.3969/j.issn.2095-1000.2019.02.004 [16] ZHENG Ge, YAO Yiping, WANG Dongdong, et al. Study of an access protocol for satellite network with open TDMA[C]. 2020 International Conference on Wireless Communications and Smart Grid (ICWCSG), Qingdao, China, 2020: 37–42. [17] SUHAIL N A, LUBEGA J T, and MAIGA G. Optimization based multimedia performance to enhance blended learning experience in constrained low bandwidth environment[C]. The 4th International Conference on Hybrid Learning, Hong Kong, China, 2011: 188–199. [18] 苏春, 李乐. 基于隐半马尔科夫退化模型的非等周期预防性维修优化[J]. 东南大学学报:自然科学版, 2021, 51(2): 342–349. doi: 10.3969/j.issn.1001-0505.2021.02.022SU Chun and LI Le. Optimization of non-equal periodic preventive maintenance based on hidden semi-Markov degradation model[J]. Journal of Southeast University:Natural Science Edition, 2021, 51(2): 342–349. doi: 10.3969/j.issn.1001-0505.2021.02.022 [19] WEI Y H, LENG Quan, HAN Song, et al. RT-WiFi: Real-time high-speed communication protocol for wireless cyber-physical control applications[C]. The 34th Real-Time Systems Symposium, Vancouver, Canada, 2013: 140–149. [20] SHOAEI A D, DERAKHSHANI M, PARSAEEFARD S, et al. Learning-based hybrid TDMA-CSMA MAC protocol for virtualized 802.11 WLANs[C]. The 26th Annual International Symposium on Personal, Indoor, and Mobile Radio Communications (PIMRC), Hong Kong, China, 2015: 1861–1866. [21] 严冰, 张进, 罗亚中. 面向编队卫星的空间系绳在轨服务[J]. 系统工程与电子技术, 2021, 43(3): 806–813. doi: 10.12305/j.issn.1001-506X.2021.03.26YAN Bing, ZHANG Jin, and LUO Yazhong. On-orbit service for formation satellites with space tether[J]. Systems Engineering and Electronics, 2021, 43(3): 806–813. doi: 10.12305/j.issn.1001-506X.2021.03.26 [22] 范本尧, 刘天雄, 徐峰, 等. 全球卫星导航系统数据传输业务发展研究[J]. 航天器工程, 2016, 25(3): 1–8. doi: 10.3969/j.issn.1673-8748.2016.03.001FAN Benyao, LIU Tianxiong, XU Feng, et al. Research on data transmission service development of global satellite navigation system[J]. Spacecraft Engineering, 2016, 25(3): 1–8. doi: 10.3969/j.issn.1673-8748.2016.03.001 期刊类型引用(5)
1. 孙健健，徐建华，成海峰，祝庆霖，韩煦. 基于金属销钉封装的Ka波段固态功率放大模块研究. 电子与信息学报. 2020(12): 3074-3080 . 
本站查看2. 张德磊，沃得良，陈少华. EBG结构在微波器件设计中的应用. 信息通信. 2019(02): 33-34 . 
百度学术3. 王鹏，芦浩，王子元，刘金枝，薛茜男. 新型折线形共面EBG电源层结构的噪声抑制特性验证与分析. 北京理工大学学报. 2019(05): 514-519 . 百度学术4. 郭慧丽，杨晓邦，王素玲，王志军. 螺旋地板式电磁带隙结构的传输特性分析. 电子元件与材料. 2016(10): 58-62 . 百度学术5. 刘峰，杨曙辉，陈迎潮. 一种抑制同步开关噪声的新型S桥电磁带隙结构. 科学技术与工程. 2016(16): 70-72+94 . 百度学术其他类型引用(3)
-
百度学术
计量
- 文章访问数: 1324
- HTML全文浏览量: 933
- PDF下载量: 214
- 被引次数: 8
下载:
