搜索

x

留言板

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

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

基于高斯混合模型的无向网络重构

何瑞辉 张海峰 王欢 马闯

引用本文:
Citation:

基于高斯混合模型的无向网络重构

何瑞辉, 张海峰, 王欢, 马闯

Gaussian mixture model based reconstruction of undirected networks

He Rui-Hui, Zhang Hai-Feng, Wang Huan, Ma Chuang
PDF
HTML
导出引用
  • 从数据中推断网络的结构作为复杂网络中一个重要科学问题已得到广泛关注. 现有的网络重构方法大多将网络重构问题转化为一系列线性方程组的求解问题, 然后通过某种截断方法对每个方程组的解进行截断, 从而确定每个节点的局部结构. 然而现有的截断方法大多存在着精度不足的问题, 且少有方法衡量每个方程组解的可截断性, 即节点的可重构性. 为了解决这些问题, 本文提出了一种基于高斯混合模型的无向网络重构方法. 该方法首先将节点间连接关系的推断问题转化为一个聚类问题, 然后利用高斯混合模型进行求解, 得到每个节点与其他节点的连接概率, 并根据概率定义一个基于信息熵的可重构指标, 从而在真实网络结构未知的情况下衡量每个节点的可重构性. 将该方法用于无向网络中, 可以利用无向网络的对称特征, 将可重构性高的节点作为训练集指导可重构性低的节点进行结构推断, 从而更好地重构出无向网络. 最后, 通过在合成数据和真实数据上与现有的截断方法进行比较, 证明了该方法可以更有效地重构出网络结构.
    The reconstruction of network structure from data represents a significant scientific challenge in the field of complex networks, which has attracted considerable attention from the research community. The most of existing network reconstruction methods transform the problem into a series of linear equation systems, to solve the equations. Subsequently, truncation methods are used to determine the local structure of each node by truncating the solution of each equation system. However, truncation methods frequently exhibit inadequate accuracy, and lack methods of evaluating the truncatability of solutions to each system of equations, that is to say, the reconstructability of nodes. In order to address these issues, in this work an undirected network reconstruction method is proposed based on a Gaussian mixture model. In this method, a Gaussian mixture model is first used to cluster the solution results obtainedby solving a series of linear equations, and then the probabilities of the clustering results are utilized to depict the likelihood of connections between nodes. Subsequently, an index of reconstructibility is defined based on information entropy, thus the probability of connections between each node and other nodes can be used to measure the reconstructibility of each node. The proposed method is ultimately applied to undirected networks. Nodes identified with high reconstructibility are used as a training set to guide the structural inference of nodes with lower reconstrucibility, thus enhancing the reconstruction of the undirected network. The symmetrical properties of the undirected network are then employed to infer the connection probabilities of the remaining nodes with other nodes. The experiments on both synthetic and real data are conducted and a variety of methods are used for constructing linear equations and diverse dynamical models. Compared with the results from a previous truncated reconstruction method, the reconstruction outcomes are evaluated. The experimental results show that the method proposed in this work outperforms existing truncation reconstruction methods in terms of reconstruction performance, thus confirming the universality and effectiveness of the proposed method.
      通信作者: 马闯, chuang_m@126.com
    • 基金项目: 国家自然科学基金(批准号: 12005001, 61973001)资助的课题.
      Corresponding author: Ma Chuang, chuang_m@126.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 12005001, 61973001).
    [1]

    Li X, Sun L, Ling M J, Peng Y 2023 Neurocomputing 549 126441Google Scholar

    [2]

    张彦超, 刘云, 张海峰, 程辉, 熊菲 2011 60 050501Google Scholar

    Zhang Y C, Liu Y, Zhang H F, Cheng H, Xiong F 2011 Acta Phys. Sin. 60 050501Google Scholar

    [3]

    Gardner T S, Di Bernardo D, Lorenz D, Collins J J 2003 Science 301 102Google Scholar

    [4]

    Geier F, Timmer J, Fleck C 2007 BMC Syst. Biol. 1 1Google Scholar

    [5]

    Gao C, Fan Y, Jiang S H, Deng Y, Liu J M, Li X H 2021 IEEE Trans. Intell. Transp. Syst. 23 6509Google Scholar

    [6]

    Zhou Y M, Li S P, Kundu T, Bai X W, Qin W 2021 IEEE Trans. Network Sci. Eng. 8 2249Google Scholar

    [7]

    张海峰, 王文旭 2020 69 088906Google Scholar

    Zhang H F, Wang W X 2020 Acta Phys. Sin. 69 088906Google Scholar

    [8]

    Wang J Y, Zhang Y J, Xu C, Li J Z, Sun J C, Xie J R, Feng L, Zhou T S, Hu Y Q 2024 Nat. Commun. 15 2849Google Scholar

    [9]

    康玲, 项冰冰, 翟素兰, 鲍中奎, 张海峰 2018 67 198901Google Scholar

    Kang L, Xiang B B, Zhai S L, Bao Z K, Zhang H F 2018 Acta Phys. Sin. 67 198901Google Scholar

    [10]

    Xiang B B, Bao Z K, Ma C, Zhang X Y, Chen H S, Zhang H F 2018 Chaos: An Interdisciplinary Journal of Nonlinear Science 28 013122Google Scholar

    [11]

    Zhao J, Cheong K H 2024 IEEE Trans. Syst. Man Cybern. Part A Syst. Humans 54 6Google Scholar

    [12]

    Guo Q T, Jiang X, Lei Y J, Li M, Ma Y F, Zheng Z M 2015 Phys. Rev. E 91 012822Google Scholar

    [13]

    Li D D, Qian W Q, Sun X X, Han D, Sun M 2023 Appl. Math. Comput. 458 128233

    [14]

    Lv X J, Fan D M, Li Q, Wang J L, Zhou L 2023 Physica A 627 129131Google Scholar

    [15]

    徐翔, 朱承, 朱先强 2021 70 088901Google Scholar

    Xu X, Zhu C, Zhu X Q 2021 Acta Phys. Sin. 70 088901Google Scholar

    [16]

    Wang H, Ma C, Chen H S, Lai Y C, Zhang H F 2022 Nat. Commun. 13 3043Google Scholar

    [17]

    Ma C, Wang H, Zhang H F 2023 Europhys. Lett. 144 21002Google Scholar

    [18]

    杨浦, 郑志刚 2012 61 120508Google Scholar

    Yang P, Zheng Z G 2012 Acta Phys. Sin. 61 120508Google Scholar

    [19]

    Ma C, Chen H S, Li X, Lai Y C, Zhang H F 2020 SIAM J. Appl. Dyn. Syst. 19 124Google Scholar

    [20]

    Shen Z S, Wang W X, Fan Y, Di Z R, Lai Y C 2014 Nat. Commun. 5 4323Google Scholar

    [21]

    Liu Q M, Ma C, Xiang B B, Chen H S, Zhang H F 2019 IEEE Trans. Syst. Man Cybern. Part A Syst. Humans 51 4639Google Scholar

    [22]

    Zhang A B, Fan Y, Di Z R, Zeng A 2023 Chaos, Solitons Fractals 173 113712Google Scholar

    [23]

    Wang W X, Lai Y C, Grebogi C, Ye J P 2011 Phys. Rev. X 1 021021

    [24]

    Li G J, Li N, Liu S H, Wu X Q 2019 Chaos: An Interdisciplinary Journal of Nonlinear Science 29 053117Google Scholar

    [25]

    Mei G F, Wu X Q, Wang Y F, Hu M, Lu J A, Chen G R 2017 IEEE Trans. Cybern. 48 754Google Scholar

    [26]

    Pandey P K, Adhikari B 2017 IEEE Trans. Knowl. Data Eng. 29 2072Google Scholar

    [27]

    Pandey P K, Adhikari B, Mazumdar M, Ganguly N 2020 IEEE Trans. Knowl. Data Eng. 34 3377Google Scholar

    [28]

    Ma C, Chen H S, Lai Y C, Zhang H F 2018 Phys. Rev. E 97 022301

    [29]

    Zhang Z, Zhao Y, Liu J, Wang S, Tao R, Xin R, Zhang J 2019 Appl. Network Sci. 4 1Google Scholar

    [30]

    Xu X, Zhu X Q, Zhu C 2023 Complex Intell. Syst. 9 3131Google Scholar

    [31]

    Mignone P, Pio G, D’ Elia D, Ceci M 2020 Bioinformatics 36 1553Google Scholar

    [32]

    Reynolds D A 2009 Encyclopedia of Biometrics 741 659

    [33]

    Wang Y, Chakrabarti D, Wang C X, Faloutsos C 2003 In 22nd International Symposium on Reliable Distributed Systems Florence, Italy, October 6–8, 2003 pp25–34

    [34]

    Perotti J I, Tessone C J, Clauset A, Caldarelli G 2018 arXiv: 1806.07005 (Physics and Society

    [35]

    Erds P, Rényi A 1960 Publ. Math. Inst. Hungar. Acad. Sci. 5 17

    [36]

    Watts D J, Strogatz S H 1998 Nature 393 440Google Scholar

    [37]

    Barabási A L, Albert R 1999 Science 286 509Google Scholar

    [38]

    Li J W, Shen Z S, Wang W X, Grebogi C, Lai Y C 2017 Phys. Rev. E 95 032303

  • 图 1  网络重构结果 (a)和(b)分别为求解两组一系列线性方程组的求解结果, 红色点表示存在连接关系的求解结果, 蓝色点表示不存在连接关系的求解结果; (c)为(a)中红色方框内节点求解结果的直方图分布; (d)为(b)中红色方框内节点求解结果的直方图分布. (a)和(b)中的横坐标表示网络中节点编号, 纵坐标表示线性方程组的求解结果, (c)和(d)中的横坐标表示线性方程组的求解结果, 纵坐标表示分布的数量

    Fig. 1.  Network reconstruction results: (a) and (b) represent two different solution results for solving a series of linear equation systems, respectively, with red dots indicating solutions with connectivity and blue dots indicating solutions without connectivity; (c) represents histogram distribution of the solution result for the node within the red box in panel (a); (d) represents histogram distribution of the solution result for the node within the red box in panel (b). The horizontal axes in panels (a) and (b) represent the node number in the network, the vertical axes represent the solution results of the linear equation system, the horizontal axes in panels (c) and (d) represent the solution results of the linear equation system, and the vertical axes represent the number of distributions.

    图 2  $ {H(i)} $与节点重构效果的关系 横坐标表示节点编号, 左纵坐标表示节点的可重构性指标负值$ {-H(i)} $, 右纵坐标表示重构效果F1

    Fig. 2.  The relationship between $ {H(i)} $ and node reconstruction effect: The horizontal axis represents the node number, the left vertical axis represents the negative value of node’s reconfigurability index $ {-H(i)} $, and the right vertical axis represents the reconstruction effect F1

    图 3  UNRGMM与TTM在合成网络中的重构效果比较 (a)和(d)为ER网络上的重构效果; (b)和(e)为WS网络上的重构效果; (c)和(f)为BA网络上的重构效果. 误差棒表示10次独立实验的标准差

    Fig. 3.  Comparison of reconstruction effects between UNRGMM and TTM in synthetic networks: (a) and (d) represent the reconstruction effect on the ER network; (b) and (e) represent the reconstruction effect on the WS network; (c) and (f) represent the reconstruction effect on the BA network. The error bar represents standard deviation over ten independent trials.

    图 4  UNRGMM与TTM在噪声干扰下的重构效果 (a)和(d)为ER网络上的重构效果; (b)和(e)为WS网络上的重构效果; (c)和(f)为BA网络上的重构效果. 误差棒表示10次独立实验的标准差

    Fig. 4.  The reconstruction effect of UNRGMM and TTM under noise interference: (a) and (d) represent the reconstruction effect on the ER network; (b) and (e) represent the reconstruction effect on the WS network; (c) and (f) represent the reconstruction effect on the BA network. The error bar represents standard deviation over ten independent trials.

    图 5  UNRGMM与TTM在不同平均度的合成网络中的重构效果 (a)和(d)为ER网络上的重构效果; (b)和(e)为WS网络上的重构效果; (c)和(f)为BA网络上的重构效果. 误差棒表示10次独立实验的标准差

    Fig. 5.  The reconstruction effect of UNRGMM and TTM in synthetic networks with different average degrees: (a) and (d) represent the reconstruction effect on the ER network; (b) and (e) represent the reconstruction effect on the WS network; (c) and (f) represent the reconstruction effect on the BA network. The error bar represents standard deviation over ten independent trials.

    图 6  UNRGMM与TTM在不同动力学下的重构效果 (a)和(d)为Ising动力学的重构效果; (b)和(e)为Game动力学的重构效果; (c)和(f)为Majority动力学的重构效果. 误差棒表示10次独立实验的标准差

    Fig. 6.  The reconstruction effect of UNRGMM and TTM under different dynamics: (a) and (d) represent the reconstruction effect for Ising dynamics; (b) and (e) represent the reconstruction effect for Game dynamics; (c) and (f) represent the reconstruction effect for Majority dynamics. The error bar represents standard deviation over ten independent trials.

    表 1  真实网络结构特征NE分别是节点和连边数量; ${\langle k\rangle}$表示平均度; Cr分别是聚类系数和分类系数; H是异质性程度, 定义为${H= {\langle k^2\rangle}/{{\langle k\rangle}^{2}}}$

    Table 1.  Real networks structure characteristics. N and E are the number of nodes and edges, respectively; ${\langle k\rangle}$ indicates the average degree; C and r are clustering coefficients and classification coefficients, respectively; H is the degree heterogeneity, defined as ${H= {\langle k^2\rangle}/{{\langle k\rangle}^{2}}}$.

    真实网络NE$ {\langle k\rangle} $CrH
    Karate34784.5880.59–0.4761.693
    Dolphins621595.1290.29–0.0711.326
    Football11561310.6610.40.1621.007
    Polbooks1054418.400.49–0.1281.421
    下载: 导出CSV

    表 2  UNRGMM与TTM在真实网络中的重构结果比较

    Table 2.  Comparison of reconstruction results between UNRGMM and TTM on real networks.

    真实网络 TTM UNRGMM
    F1 Accuracy F1 Accuracy
    Karate 0.723 0.889 0.989 0.997
    Dolphins 0.776 0.951 0.970 0.995
    Football 0.484 0.874 0.774 0.963
    Polbooks 0.571 0.896 0.838 0.978
    下载: 导出CSV
    Baidu
  • [1]

    Li X, Sun L, Ling M J, Peng Y 2023 Neurocomputing 549 126441Google Scholar

    [2]

    张彦超, 刘云, 张海峰, 程辉, 熊菲 2011 60 050501Google Scholar

    Zhang Y C, Liu Y, Zhang H F, Cheng H, Xiong F 2011 Acta Phys. Sin. 60 050501Google Scholar

    [3]

    Gardner T S, Di Bernardo D, Lorenz D, Collins J J 2003 Science 301 102Google Scholar

    [4]

    Geier F, Timmer J, Fleck C 2007 BMC Syst. Biol. 1 1Google Scholar

    [5]

    Gao C, Fan Y, Jiang S H, Deng Y, Liu J M, Li X H 2021 IEEE Trans. Intell. Transp. Syst. 23 6509Google Scholar

    [6]

    Zhou Y M, Li S P, Kundu T, Bai X W, Qin W 2021 IEEE Trans. Network Sci. Eng. 8 2249Google Scholar

    [7]

    张海峰, 王文旭 2020 69 088906Google Scholar

    Zhang H F, Wang W X 2020 Acta Phys. Sin. 69 088906Google Scholar

    [8]

    Wang J Y, Zhang Y J, Xu C, Li J Z, Sun J C, Xie J R, Feng L, Zhou T S, Hu Y Q 2024 Nat. Commun. 15 2849Google Scholar

    [9]

    康玲, 项冰冰, 翟素兰, 鲍中奎, 张海峰 2018 67 198901Google Scholar

    Kang L, Xiang B B, Zhai S L, Bao Z K, Zhang H F 2018 Acta Phys. Sin. 67 198901Google Scholar

    [10]

    Xiang B B, Bao Z K, Ma C, Zhang X Y, Chen H S, Zhang H F 2018 Chaos: An Interdisciplinary Journal of Nonlinear Science 28 013122Google Scholar

    [11]

    Zhao J, Cheong K H 2024 IEEE Trans. Syst. Man Cybern. Part A Syst. Humans 54 6Google Scholar

    [12]

    Guo Q T, Jiang X, Lei Y J, Li M, Ma Y F, Zheng Z M 2015 Phys. Rev. E 91 012822Google Scholar

    [13]

    Li D D, Qian W Q, Sun X X, Han D, Sun M 2023 Appl. Math. Comput. 458 128233

    [14]

    Lv X J, Fan D M, Li Q, Wang J L, Zhou L 2023 Physica A 627 129131Google Scholar

    [15]

    徐翔, 朱承, 朱先强 2021 70 088901Google Scholar

    Xu X, Zhu C, Zhu X Q 2021 Acta Phys. Sin. 70 088901Google Scholar

    [16]

    Wang H, Ma C, Chen H S, Lai Y C, Zhang H F 2022 Nat. Commun. 13 3043Google Scholar

    [17]

    Ma C, Wang H, Zhang H F 2023 Europhys. Lett. 144 21002Google Scholar

    [18]

    杨浦, 郑志刚 2012 61 120508Google Scholar

    Yang P, Zheng Z G 2012 Acta Phys. Sin. 61 120508Google Scholar

    [19]

    Ma C, Chen H S, Li X, Lai Y C, Zhang H F 2020 SIAM J. Appl. Dyn. Syst. 19 124Google Scholar

    [20]

    Shen Z S, Wang W X, Fan Y, Di Z R, Lai Y C 2014 Nat. Commun. 5 4323Google Scholar

    [21]

    Liu Q M, Ma C, Xiang B B, Chen H S, Zhang H F 2019 IEEE Trans. Syst. Man Cybern. Part A Syst. Humans 51 4639Google Scholar

    [22]

    Zhang A B, Fan Y, Di Z R, Zeng A 2023 Chaos, Solitons Fractals 173 113712Google Scholar

    [23]

    Wang W X, Lai Y C, Grebogi C, Ye J P 2011 Phys. Rev. X 1 021021

    [24]

    Li G J, Li N, Liu S H, Wu X Q 2019 Chaos: An Interdisciplinary Journal of Nonlinear Science 29 053117Google Scholar

    [25]

    Mei G F, Wu X Q, Wang Y F, Hu M, Lu J A, Chen G R 2017 IEEE Trans. Cybern. 48 754Google Scholar

    [26]

    Pandey P K, Adhikari B 2017 IEEE Trans. Knowl. Data Eng. 29 2072Google Scholar

    [27]

    Pandey P K, Adhikari B, Mazumdar M, Ganguly N 2020 IEEE Trans. Knowl. Data Eng. 34 3377Google Scholar

    [28]

    Ma C, Chen H S, Lai Y C, Zhang H F 2018 Phys. Rev. E 97 022301

    [29]

    Zhang Z, Zhao Y, Liu J, Wang S, Tao R, Xin R, Zhang J 2019 Appl. Network Sci. 4 1Google Scholar

    [30]

    Xu X, Zhu X Q, Zhu C 2023 Complex Intell. Syst. 9 3131Google Scholar

    [31]

    Mignone P, Pio G, D’ Elia D, Ceci M 2020 Bioinformatics 36 1553Google Scholar

    [32]

    Reynolds D A 2009 Encyclopedia of Biometrics 741 659

    [33]

    Wang Y, Chakrabarti D, Wang C X, Faloutsos C 2003 In 22nd International Symposium on Reliable Distributed Systems Florence, Italy, October 6–8, 2003 pp25–34

    [34]

    Perotti J I, Tessone C J, Clauset A, Caldarelli G 2018 arXiv: 1806.07005 (Physics and Society

    [35]

    Erds P, Rényi A 1960 Publ. Math. Inst. Hungar. Acad. Sci. 5 17

    [36]

    Watts D J, Strogatz S H 1998 Nature 393 440Google Scholar

    [37]

    Barabási A L, Albert R 1999 Science 286 509Google Scholar

    [38]

    Li J W, Shen Z S, Wang W X, Grebogi C, Lai Y C 2017 Phys. Rev. E 95 032303

  • [1] 张赞, 黄北举, 陈弘达. 基于可重构硅光滤波器的计算重建片上光谱仪.  , 2024, 73(14): 140701. doi: 10.7498/aps.73.20240224
    [2] 罗恺明, 管曙光, 邹勇. 基于相位同步动力学重构网络单纯复形的相互作用.  , 2024, 73(12): 120501. doi: 10.7498/aps.73.20240334
    [3] 刘宇航, 林曈, 李少波, 于文琦, 马向, 梁晓东, 恽斌峰. 可调反射器辅助的可重构微环光滤波器.  , 2023, 72(8): 084208. doi: 10.7498/aps.72.20222384
    [4] 谷同凯, 王兰兰, 国阳, 蒋维涛, 史永胜, 杨硕, 陈金菊, 刘红忠. 光盘上集成的液体微透镜阵列与可重构超分辨成像.  , 2023, 72(9): 099501. doi: 10.7498/aps.72.20222251
    [5] 黄晓俊, 高焕焕, 何嘉豪, 栾苏珍, 杨河林. 动态可调谐的频域多功能可重构极化转换超表面.  , 2022, 71(22): 224102. doi: 10.7498/aps.71.20221256
    [6] 王明照, 王少杰, 许河秀. 基于剪纸方法的一种可重构线极化转换空间序构超表面.  , 2021, 70(15): 154101. doi: 10.7498/aps.70.20210188
    [7] 张娜, 赵健民, 陈克, 赵俊明, 姜田, 冯一军. 编码超构表面实现双波束独立可重构.  , 2021, 70(17): 178102. doi: 10.7498/aps.70.20210344
    [8] 徐翔, 朱承, 朱先强. 一种基于离散数据从局部到全局的网络重构算法.  , 2021, 70(8): 088901. doi: 10.7498/aps.70.20201756
    [9] 王晨阳, 段倩倩, 周凯, 姚静, 苏敏, 傅意超, 纪俊羊, 洪鑫, 刘雪芹, 汪志勇. 基于遗传算法优化卷积长短记忆混合神经网络模型的光伏发电功率预测.  , 2020, 69(10): 100701. doi: 10.7498/aps.69.20191935
    [10] 谭索怡, 祁明泽, 吴俊, 吕欣. 复杂网络链路可预测性: 基于特征谱视角.  , 2020, 69(8): 088901. doi: 10.7498/aps.69.20191817
    [11] 张海峰, 王文旭. 复杂系统重构.  , 2020, 69(8): 088906. doi: 10.7498/aps.69.20200001
    [12] 于惠存, 曹祥玉, 高军, 杨欢欢, 韩江枫, 朱学文, 李桐. 一种宽带可重构反射型极化旋转表面.  , 2018, 67(22): 224101. doi: 10.7498/aps.67.20181041
    [13] 孔江涛, 黄健, 龚建兴, 李尔玉. 基于复杂网络动力学模型的无向加权网络节点重要性评估.  , 2018, 67(9): 098901. doi: 10.7498/aps.67.20172295
    [14] 徐明, 许传云, 曹克非. 度相关性对无向网络可控性的影响.  , 2017, 66(2): 028901. doi: 10.7498/aps.66.028901
    [15] 梁潇, 钱志鸿, 田洪亮, 王雪. 基于马尔可夫决策模型的异构无线网络切换选择算法.  , 2016, 65(23): 236402. doi: 10.7498/aps.65.236402
    [16] 谢文科, 高穹, 马浩统, 魏文俭, 江文杰. 基于模型选择的模式波前重构算法研究.  , 2015, 64(14): 144201. doi: 10.7498/aps.64.144201
    [17] 黄锦旺, 李广明, 冯久超, 晋建秀. 一种无线传感器网络中的混沌信号重构算法.  , 2014, 63(14): 140502. doi: 10.7498/aps.63.140502
    [18] 李黎, 郑庆华, 管晓宏. 基于有限资源提升网络可生存性的拓扑重构方法.  , 2014, 63(17): 170201. doi: 10.7498/aps.63.170201
    [19] 陈帝伊, 柳烨, 马孝义. 基于径向基函数神经网络的混沌时间序列相空间重构双参数联合估计.  , 2012, 61(10): 100501. doi: 10.7498/aps.61.100501
    [20] 梁玉, 郭立新, 王蕊. 粗糙面重构问题的混合算法研究.  , 2011, 60(3): 034102. doi: 10.7498/aps.60.034102
计量
  • 文章访问数:  556
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-04-22
  • 修回日期:  2024-07-22
  • 上网日期:  2024-08-05
  • 刊出日期:  2024-09-05

/

返回文章
返回
Baidu
map