GPS/BDS Precise Point Positioning Model with Receiver DCB Parameters for Raw Observations
-
摘要: 在传统多系统非差非组合精密单点定位（precise point positioning，PPP）模型中，电离层延迟会吸收部分接收机码硬件延迟，其估计值可能为负数。提出了一种估计接收机差分码偏差（differential code bias，DCB）参数的GPS（Global Positioning System）/BDS（BeiDou Navigation Satellite System）非组合PPP模型，将每个系统第1个频率上的接收机码硬件延迟约束为零，对接收机DCB进行参数估计，达到了分离电离层延迟和接收机码硬件延迟的目的，降低了接收机钟差和电离层延迟的相关程度。利用4个多星座实验（multi-GNSS experiment，MGEX）跟踪站的GPS/BDS数据进行了静态和动态PPP试验，结果表明，与不估计DCB参数的PPP模型相比，采用估计DCB参数PPP模型后，静态模式下定位精度和收敛速度平均提高了29.3%和29.8%，动态模式下定位精度和收敛速度平均提高了15.7%和21.6%。Abstract: ue to ionospheric delays caused by the receiver code biases in the traditional multi-GNSS precise point positioning (PPP) using raw observations, the estimates can be negative values. An improved model of Global Positioning System/BeiDou Navigoction Satellite System (GPS/BDS) PPP with receiver differential code bias (DCB) parameters for raw observations is proposed in which receiver code biases on the first frequency of each system are constrained to zero and receiver DCB parameters are estimated. Ionospheric delays and receiver code biases are separated by the presented model. Additionally, the singularity between receiver clock offsets and ionospheric delays is reduced. GPS/BDS data from 4 stations of the Multi-GNSS experiment (MGEX) network are processed in static and kinematic modes. The results show that with the proposed PPP model, the average positioning accuracy and convergence time in static/kinematic mode are improved by 29.3% and 15.7%, 29.8% and 21.6%, respectively, in comparison with the traditional PPP model.
-
目前，北斗导航卫星系统(BDS)已实现局域覆盖，随着系统建设的不断完善和应用的不断拓展，与之相关的各类数据处理软件的开发成为重要的研究内容。因此，自主开发北斗高精度数据处理软件，成为发展高精度位置服务的迫切任务[1-8]。因北斗导航卫星系统与GPS在星座构造、坐标框架、时间系统、信号频率等方面具有明显差异[9-15]，现有的高精度GPS数据处理软件无法直接处理北斗数据。本文针对北斗高精度数据处理的系统设计、数据流、功能模块及高精度算法实现等进行了研究，研制开发了一套高精度北斗基线解算软件BGO(BeiDou Navigation Satellite System/Global Positioning System Office)，并将其用于高速铁路高精度控制测量建网。通过与商业软件TGO(Trimble Geomatics Office)和TBC(Trimble Business Center)，及高精度科研软件Bernese进行对比测试、性能分析，验证了该软件的正确性和有效性。
1 系统的设计与模块算法的实现
1.1 系统设计与数据流分析
北斗和GPS基线解算软件主要包含北斗基线处理、GPS基线处理及联合基线处理3大模块。各模块间相互独立，但使用相同的数据结构，且数据流基本一致。数据处理流程如图 1所示。
基线解算之前，需选择有效双频观测数据，具体包含低高度角卫星剔除、观测值粗差剔除、星历未获取观测数据剔除等。剔除质量较差的观测数据可通过可视化的方式实现。通过双频数据组合有效消除电离层延迟影响，伪距消电离组合能算出测站精确至10 m内的概略位置，从而形成网络拓扑图，便于用户查看站点的平面分布。基线解算时，北斗与GPS独立系统数据处理算法相同；联合处理需选择统一的坐标和时间框架，随着多余观测数的增加，还需设置合理的模糊度固定限值。基线解算后，进行网平差，应剔除不合格基线，直至平差结果满足要求。
1.2 高精度基线解算算法实现
高精度基线解算利用双差观测量建立误差方程，北斗双差观测量构造如式(1)：
$$ \mathit{\Delta} \nabla L^{{C_m}{C_n}}_{{S_i}{S_j}} = \left( {L^{{C_n}}_{{S_j}} - L^{{C_n}}_{{S_i}}} \right) - \left( {L^{{C_m}}_{{S_j}} - L^{{C_m}}_{{S_i}}} \right) $$ (1) 式中，Δ▽L表示双差观测量；Si和Sj表示任意站点；Cm和Cn表示任意北斗卫星。
依据式(1)构建的双差观测量，建立误差方程，如式(2)：
$$ \left[ \begin{array}{l} \mathit{\Delta} \nabla \boldsymbol{\varPhi} \\ \mathit{\Delta} \nabla \boldsymbol{P} \end{array} \right] = \boldsymbol{BX} + \boldsymbol{A}\mathit{\Delta} \nabla \boldsymbol{N} + \boldsymbol{V} $$ (2) 式中，Δ▽Φ和Δ▽P分别表示卫星载波相位和伪距双差观测量；X表示基线向量；Δ▽N表示双差整周模糊度；B和A为系数阵；V为残差向量。
利用式(2)构建的误差方程，解算基线向量和双差整周模糊度浮点解。利用LAMBAD方法[16, 17]固定双差整周模糊度后去除。再利用载波相位观测值获取高精度基线向量结果。基线解算过程中，主要利用抗差估计的切比雪夫多项式拟合法[18]及MW-GF组合法[19]探测与修复周跳。
对北斗和GPS双系统基线解算，只需将各系统的双差观测量误差方程叠加后平差计算，即可实现双系统联合基线解算。但需注意，星间差分需选择同一系统卫星，否则会引入系统间信号硬件延迟[20]，影响双差整周模糊度的固定。另外，北斗和GPS在时间框架、坐标框架等存在一定差异，双系统联合解算需保证框架的统一。
北斗和GPS时间转换公式如式(3)：
$$ {t_C} = {t_G}-14\;{\rm{s}} $$ (3) 式中，tC和tG分别表示北斗时和GPS时，两者均为原子时，起算原点不同[13]。
北斗和GPS坐标转换公式如式(4)：
$$ \begin{array}{c} \left[ {\begin{array}{*{20}{c}} {{X_C}}\\ {{Y_C}}\\ {{Z_C}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{X_G}}\\ {{Y_G}}\\ {{Z_G}} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{T_X}}\\ {{T_Y}}\\ {{T_Z}} \end{array}} \right] + \\ \left[ {\begin{array}{*{20}{c}} D&{ - {R_Z}}&{{R_Y}}\\ {{R_Z}}&D&{ - {R_X}}\\ { - {R_Y}}&{{R_X}}&D \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{X_G}}\\ {{Y_G}}\\ {{Z_G}} \end{array}} \right] \end{array} $$ (4) 式中，北斗坐标(XC，YC，ZC)与GPS坐标(XG，YG，ZG)可通过七参数TX、TY、TZ、D、RX、RY、RZ进行转换。北斗CGCS2000坐标系采用ITRF97框架2000历元的坐标和速度场，当前GPS WGS84坐标和ITRF08基本一致。因此，可利用ITRF97框架2000历元与ITRF08间转换的七参数(ITRF网站公布)实现北斗与GPS坐标框架的统一[11, 12]。
2 BGO数据处理实例与性能测试
2.1 高速铁路CPI控制网基线解算
处理高速铁路CPI控制网时，通过读取观测文件和星历文件，单点定位生成控制网的基线网络拓扑图，如图 2所示。基线解算前，设置相关参数包括卫星截止高度角、误差限差参数、框架、对流层模型、电离层模型、模糊度Ratio值、同步最小观测历元数等。设置完成后，可选择北斗、GPS、联合3种模式进行基线解算。基线解算完成后，软件界面中将显示解算的基线分量及其精度，并可显示残差向量检核基线解算效果。
2.2 BGO、TGO、Bernese软件处理GPS基线结果比较
为了测试BGO解算GPS基线的正确性，将其与TGO和Bernese软件处理结果进行了比较，得到57条GPS基线(基线最长6 667 m，最短446 m)的比较结果，如图 3所示。
![]()
图 3 BGO、TGO、Bernese软件处理GPS基线分量比较Figure 3. Comparing GPS Baseline Components from BGO, TGO and Bernese Software图 3(a)、3(b)分别表示BGO软件与TGO、Bernese软件处理GPS基线分量的差值ΔX、ΔY、ΔZ。图 3(a)中，BGO和TGO有52条基线在X、Y、Z方向的分量差值均在2 cm内，有48条基线各分量差值在mm级。TGO解算少量基线验后方差分量超限，与BGO基线分量差值较大。图 3(b)中，BGO和Bernese有55条基线在X、Y、Z方向的分量差值均在2 cm内，有49条基线各分量差值在mm级。
图 4(a)~4(c)分别表示BGO、TGO、Bernese软件处理GPS基线的内符合精度σX、σY、σZ(BGO、TGO、Bernese软件基线解算精度分别精确至0.1 mm、1 mm和0.1 mm)。整体上，约90%的基线3个软件的解算精度相当。
![]()
图 4 BGO、TGO、Bernese的GPS基线内符合精度比较Figure 4. Comparing GPS Baseline Precision from BGO, TGO and Bernese Software2.3 BGO、TBC软件处理北斗与GPS联合基线结果
为了测试BGO解算北斗与GPS联合基线的性能，本文选用美国Trimble的商业软件TBC与之进行比较。同上57条基线，每条基线观测数据均包含北斗与GPS观测数据。图 5展示了BGO和TBC处理北斗与GPS联合基线分量的差值ΔX、ΔY、ΔZ。图 5可见，98%的基线分量差值分布在mm级，表明BGO软件处理联合基线能达到与TBC软件相当的水平。另外，两者内符合精度绝大部分均在mm级，故图 5中未加以比较。
![]()
图 5 BGO与TBC软件处理北斗与GPS联合基线分量比较Figure 5. Comparing BDS and GPS Combined Baseline Components from BGO and TBC Software由此可知，BGO软件处理GPS基线、北斗与GPS联合基线的内外符合精度能达到TGO、Bernese、TBC相当的水平。因此，以BGO软件处理GPS、北斗与GPS联合基线结果为参考值，分析该软件处理北斗基线结果的正确性和可靠性，如图 6和图 7所示。图 6比较了北斗与GPS、联合基线分量的差值，图 7比较了北斗、GPS、联合基线解算的内符合精度。
![]()
图 6 BGO软件处理北斗与GPS、联合基线分量比较Figure 6. Comparing BDS, GPS and BDS/GPS Combined Baseline Components from BGO Software![]()
图 7 北斗、GPS、联合基线解的内符合精度统计Figure 7. The Statistics of Precision of BDS, GPS and BDS/GPS Combined Baseline Solutions图 6(a)表示BGO软件处理北斗与GPS基线分量的差值ΔX、ΔY、ΔZ，其中有43条基线在X、Y、Z方向上的分量差值Δx、Δy、Δz在2 cm内，有31条基线在X、Y、Z方向上的分量差值在mm级。图 6(b)表示BGO软件处理北斗与联合基线分量的差值，其中有54条基线在X、Y、Z方向上的分量差值在2 cm内，有38条基线在X、Y、Z方向上的分量差值在mm级(图 6中第6条基线北斗为浮点解，各分量差值结果较大，图中置为0)。
图 7中，93%的联合基线在X、Y、Z方向上的分量精度分别优于0.5 mm、1 mm、0.5 mm；约90%的北斗基线和95%的GPS基线在X、Y、Z方向上的分量精度分别优于1 mm、2 mm、1 mm。由北斗、GPS、联合基线3者精度比较可知，在北斗试运行阶段，GPS基线内符合精度略优于北斗，北斗与GPS联合系统基线内符合精度明显高于独立系统。
2.4 BGO基线网平差及其精度分析
BGO具备网平差功能，根据网平差后的基线分量改正数、相对中误差、点位精度等判断基线解算结果的可靠性。对上述解算的北斗、GPS、联合基线分别进行无约束网平差。
北斗、GPS、联合基线无约束网平差的平差改正数δX、δY、δZ绝大部分在±1 cm内，如图 8(a)~8(c)所示。最弱边相对中误差优于5.5 ppm(规范限值)，具体见表 1。据图 8、表 1及《高速铁路工程测量规范》[21]可知，BGO能合理稳定地解算北斗、GPS及联合基线，解算结果中的基线向量改正数、最弱边相对中误差、最弱点点位精度均满足CPI控制测量要求，各系统解算均能精确获得24个CPI控制点坐标。
![]()
图 8 GPS、北斗、联合无约束网平差基线向量改正数Figure 8. Baseline Vector Corrections from GPS, BDS and BDS/GPS Combined Unconstrained Adjustment表 1 GPS、北斗、联合无约束平差结果统计Table 1. The Statistics of GPS, BDS and BDS/GPS Combined Unconstrained Adjustment Results解算模式 独立基线 多余观测数 控制点个数 最弱边相对中误差/ppm 最弱点点位精度/mm GPS 55 66 24 3.6 23.6 北斗 51 57 24 3.1 26.9 联合 57 72 24 3.7 17.9 3 结语
本文系统地研究了北斗与GPS联合基线解算的算法，自主开发了北斗高精度基线解算软件BGO。通过实测高铁CPI控制网的数据处理测试表明：软件能进行高精度地处理北斗与GPS数据, 以及北斗与GPS联合数据处理；GPS基线解算性能与天宝TGO软件相当，能达到与Bernese软件一致的精度；北斗与GPS基线处理能达到与TBC相当的水平。BGO最大的优势在于能对北斗和GPS进行联合解算，从而提高北斗或GPS单系统的基线解算合格率和精度。经高速铁路CPI控制网实例测试，证明该软件处理基线结果可用于高精度北斗和GPS测量控制网的数据处理。
-
![]()
图 3 2016年021天MRO1站静态定位误差序列
Figure 3. Positioning Errors of PPP in Static Mode at Station MRO1 on Day 021 in 2016
![]()
图 4 2016年021天MRO1站动态定位误差序列
Figure 4. Positioning Errors of PPP in Kinematic Mode at Station MRO1 on Day 021 in 2016
![]()
图 5 无DCB参数模型2016年022天MRO1站的垂直电离层延迟
Figure 5. Vertical Ionospheric Delays for PPP Model Without DCB Parameters at Station MRO1 on Day 022 in 2016
![]()
图 6 2016年021天MRO1站静态PPP接收机DCB序列
Figure 6. Receiver DCBs of PPP in Static Mode at Station MRO1 on Day 021 in 2016
![]()
图 7 2016年021天MRO1站动态PPP接收机DCB序列
Figure 7. Receiver DCBs of PPP in Kinematic Mode at Station MRO1 on Day 021 in 2016
表 1 数据处理策略
Table 1 The Strategy of Data Processing
类别 处理策略 观测量 伪距和相位观测量；GPS：L1，L2；BDS：B1，B2 估计方法 卡尔曼滤波 截止高度角 7° 观测值加权 采用高度角定权模型 卫星天线
相位中心GPS的天线相位中心偏差(phase center offset，PCO)和相位中心变化(phase center variation, PCV)采用igs08.atx提供的改正值；BDS的PCO和PCV采用欧洲航天局提供的改正值 接收机天线
相位中心GPS的PCO和PCV采用igs08.atx提供的改正值；BDS采用与GPS相同的PCO和PCV改正值 相位缠绕 改正 潮汐效应 改正固体潮、海潮和极潮 卫星轨道
和钟差采用德国地学中心提供的精密轨道和精密钟差产品 卫星差分
码偏差采用MGEX提供的产品 对流层延迟 Saastamoinen模型和随机游走过程 电离层延迟 作为白噪声估计[15] 接收机坐标 在静态/动态模型下作为常数/白噪声估计 接收机钟差 作为白噪声估计 接收机系统间
偏差(inter
system bias，ISB)在一天中作为常数估计[15] 接收机DCB 在一天中作为常数估计 模糊度 在连续弧段作为常数估计 
下载: 导出CSV
表 2 静态平均定位误差RMS和收敛时间
Table 2 The Average of RMS of Positioning Error and Convergence Time in Static Mode
测站 模型 E/cm N/cm U/cm 收敛时间/min CUT0 无DCB 2.2 1.3 5.7 54.9 估计DCB 1.6 1.1 2.7 25.8 IF 1.2 0.6 1.8 24.8 JFNG 无DCB 1.9 1.0 5.5 60.6 估计DCB 2.6 0.7 1.7 56.4 IF 0.9 0.5 2.7 22.4 MAL2 无DCB 1.5 1.1 2.9 56.6 估计DCB 1.5 1.1 2.9 32. 7 IF 1.0 0.9 1.9 19.4 MRO1 无DCB 1.8 1.4 5.4 31.1 估计DCB 1.9 1.5 3.4 25.8 IF 1.1 0.5 1.3 19.3
下载: 导出CSV
表 3 动态平均定位误差RMS和收敛时间
Table 3 The Average of RMS of Positioning Error and Convergence Time in Kinematic Mode
测站 模型 E/cm N/cm U/cm 收敛时间/min CUT0 无DCB 3.9 3.5 8.3 18.6 估计DCB 3.9 2.2 7.4 17.8 IF 2.8 1.5 4.4 16.8 JFNG 无DCB 3.8 4.0 9.0 43.6 估计DCB 4.6 2.7 5.4 29.3 IF 2.9 1.8 3.6 23.9 MAL2 无DCB 3.2 2.1 8.8 22.1 估计DCB 2.9 2.0 8.5 12.5 IF 2.1 1.5 4.5 10.1 MRO1 无DCB 4.4 2.6 8.0 18.0 估计DCB 3.6 1.7 6.7 16.9 IF 3.2 1.6 5.1 15.9
下载: 导出CSV
-
[1] Zumberge J, Heflin M, Jefferson D, et al. Precise Point Positioning for the Efficient and Robust Ana-lysis of GPS Data from Large Networks[J]. Journal of Geophysical Research:Solid Earth, 1997, 102(B3):5005-5017 doi: 10.1029/96JB03860
[2] 张小红, 郭斐, 李星星, 等. GPS/GLONASS组合精密单点定位研究[J].武汉大学学报·信息科学版, 2010, 35(1):9-12 http://ch.whu.edu.cn/CN/abstract/abstract823.shtml Zhang Xiaohong, Guo Fei, Li Xingxing, et al. Study on Precise Point Positioning Based on Combined GPS and GLONASS[J]. Geomatics and Information Science of Wuhan University, 2010, 35(1):9-12 http://ch.whu.edu.cn/CN/abstract/abstract823.shtml
[3] 任晓东, 张柯柯, 李星星, 等. BeiDou、Galileo、GLONASS、GPS多系统融合精密单点[J].测绘学报, 2015, 44(12):1307-1313 http://d.old.wanfangdata.com.cn/Periodical/chxb201512002 Ren Xiaodong, Zhang Keke, Li Xingxing, et al. Precise Point Positioning with Multi-constellation Satellite Systems:BeiDou, Galileo, GLONASS, GPS[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(12):1307-1313 http://d.old.wanfangdata.com.cn/Periodical/chxb201512002
[4] Chen J, Zhang Y, Wang J, et al. A Simplified and Unified Model of Multi-GNSS Precise Point Positioning[J]. Advances in Space Research, 2015, 55(1):125-134 doi: 10.1016/j.asr.2014.10.002
[5] Tu R, Ge M, Zhang H, et al. The Realization and Convergence Analysis of Combined PPP Based on Raw Observation[J]. Advances in Space Research, 2013, 52(1):211-221 doi: 10.1016/j.asr.2013.03.005
[6] Guo F, Li X, Zhang X, et al. Assessment of Precise Orbit and Clock Products for Galileo, BeiDou, and QZSS from IGS Multi-GNSS Experiment (MGEX)[J]. GPS Solutions, 2017, 21(1):279-290 doi: 10.1007/s10291-016-0523-3
[7] Li X, Ge M, Dai X, et al. Accuracy and Reliability of Multi-GNSS Real-Time Precise Positioning:GPS, GLONASS, BeiDou, and Galileo[J]. Journal of Geodesy, 2015, 89(6):607-635 doi: 10.1007/s00190-015-0802-8
[8] 章红平, 高周正, 牛小骥, 等. GPS非差非组合精密单点定位算法研究[J].武汉大学学报·信息科学版, 2013, 38(12):1396-1399 http://ch.whu.edu.cn/CN/abstract/abstract2833.shtml Zhang Hongping, Gao Zhouzheng, Niu Xiaoji, et al. Research on GPS Precise Point Positioning with Un-differential and Un-combined Observations[J]. Geomatics and Information Science of Wuhan University, 2013, 38(12):1396-1399 http://ch.whu.edu.cn/CN/abstract/abstract2833.shtml
[9] 张小红, 左翔, 李盼.非组合与组合PPP模型比较及定位性能分析[J].武汉大学学报·信息科学版, 2013, 38(5):561-565 http://ch.whu.edu.cn/CN/abstract/abstract2636.shtml Zhang Xiaohong, Zuo Xiang, Li Pan. Mathematic Model and Performance Comparison Between Ionosphere-Free Combined and Uncombined Precise Point Positioning[J]. Geomatics and Information Science of Wuhan University, 2013, 38(5):561-565 http://ch.whu.edu.cn/CN/abstract/abstract2636.shtml
[10] Zhang H, Gao Z, Ge M, et al. On the Convergence of Ionospheric Constrained Precise Point Positioning (IC-PPP) Based on Undifferential Uncombined Raw GNSS Observations[J]. Sensors, 2013, 13(11):15708-15725 doi: 10.3390/s131115708
[11] 辜声峰.多频GNSS非差非组合精密数据处理理论及其应用[D].武汉: 武汉大学, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10486-1014135106.htm Gu Shengfeng. Research on the Zero-Difference Un-combined Data Processing Model for Multi-frequency GNSS and Its Applications[D]. Wuhan: Wuhan University, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10486-1014135106.htm
[12] Li X, Zhang X, Ren X, et al. Precise Positioning with Current Multi-constellation Global Navigation Satellite Systems:GPS, GLONASS, Galileo and BeiDou[J]. Scientific Reports, 2015, 5:8328-8341 doi: 10.1038/srep08328
[13] Lou Y, Zheng F, Gu S, et al. Multi-GNSS Precise Point Positioning with Raw Single-Frequency and Dual-frequency Measurement Models[J]. GPS Solutions, 2016, 20(4):849-862 doi: 10.1007/s10291-015-0495-8
[14] 陈华.基于原始观测值的GNSS统一快速精密数据处理方法[D].武汉: 武汉大学, 2015 http://cdmd.cnki.com.cn/Article/CDMD-10486-1015305222.htm Chen Hua. An Efficient and Unified GNSS Raw Data Processing Strategy[D]. Wuhan: Wuhan University, 2015 http://cdmd.cnki.com.cn/Article/CDMD-10486-1015305222.htm
[15] Liu T, Yuan Y, Zhang B, et al. Multi-GNSS Precise Point Positioning (MGPPP) Using Raw Observations[J]. Journal of Geodesy, 2017, 91(3):253-268 doi: 10.1007/s00190-016-0960-3
[16] Hakansson M, Jensen A B O, Horemuz M, et al. Review of Code and Phase Biases in Multi-GNSS Positioning[J]. GPS Solutions, 2017, 21(3):849-860 doi: 10.1007/s10291-016-0572-7
-
期刊类型引用(6)
1. 柯文清 ，陈业滨 ，赵志刚 ，韩德志 ，郭仁忠 . 基于文献计量的新世纪地图可视化研究演变和热点分析. 地理与地理信息科学. 2025(01): 15-23 . 
百度学术
2. 于峰一泽，汤国安，陆鼎阳，林晓芬，胡光辉，沈婕，吴明光. 语言学视角下的地图演化. 地理学报. 2024(01): 171-186 . 百度学术
3. 韩德志，郭仁忠，陈业滨，赵志刚，柯文清. 基于可视化维度理论的泛地图知识推荐方法. 地球信息科学学报. 2024(01): 110-120 . 百度学术
4. 刘强，刘金花，王磊斌，陈鑫，赵志斌，赵晓艳，李英奎. 使用虚拟现实技术提高冰川地貌表达维度的研究. 冰川冻土. 2024(03): 1087-1098 . 百度学术
5. 邓志钢，郭仁忠，陈业滨，马丁，赵志刚，朱维. 面向轨迹可视化的泛地图表达维度关联方法及应用. 测绘通报. 2024(11): 56-60+96 . 百度学术
6. 柯婷，杨品福，任福，李连营，杨晨. 内河航行参考图地图符号形式化表达. 地理空间信息. 2024(12): 102-105 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 1440
- HTML全文浏览量: 187
- PDF下载量: 313
- 被引次数: 10
