闽江学院数学与数据科学学院(软件学院)，福建 福州

收稿日期：2023年6月18日；录用日期：2023年7月13日；发布日期：2023年7月20日

摘要

掌握市场趋势和规划收支费用对于财政部门而言具有极为重要的意义。本文选取1994~2021年福州市年财政收入相关数据，采用随机森林(RF)模型识别出财政收入的关键影响特征，随后建立GM(1,1)-BP组合模型，对2022~2025年福州市的年财政收入进行预测。预测结果表明RF-GM(1,1)-BP组合模型十分适合用于预测福州市财政收入；同时，福州市财政收入将在2021年之后稳步增长，并在2025年到达7,999,256万元。该结论能为相关部门实施的决策提供一定的理论参考。

关键词

财政收入，随机森林，GM(1,1)，BP，组合模型

Predicting Fuzhou City’s Fiscal Revenue Based on RF-GM(1,1)-BP Model

Wei Liu, Qiao Zhang, Wenbo Wang, Ke Xu, Huiyan Zhang, Xiuling Jin^*

School of Mathematics and Data Science (School of Software), Minjiang University, Fuzhou Fujian

Received: Jun. 18^th, 2023; accepted: Jul. 13^th, 2023; published: Jul. 20^th, 2023

ABSTRACT

Mastering market trends and planning revenue and expenditure expenses is of great significance for the financial department. This paper selects the relevant data of Fuzhou’s annual financial revenue from 1994 to 2021, uses the random forest (RF) model to identify the key impact characteristics of financial revenue, and then establishes a GM(1,1)-BP combination model to predict the annual financial revenue of Fuzhou from 2022 to 2025. The prediction results indicate that the RF-GM(1,1)-BP combination model is very suitable for predicting the fiscal revenue of Fuzhou City; Meanwhile, the fiscal revenue of Fuzhou City will steadily increase after 2021 and reach 79992.56 million yuan by 2025. This conclusion can provide a certain theoretical reference for decision-making implemented by relevant departments.

Keywords:Fiscal Revenue, Random Forest, GM(1,1), BP, Combining Model

Copyright © 2023 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

地方的财政收入作为国家财政收入的重要组成部分，它能够综合反映出一个地方的民众的经济生活水平。当政府对各项财政收支进行宏观调控时，地方政府的财政收入是其主要参考指标。因此，对于财政部门而言，对财政收入进行合理的预测能够充分有效地了解到地方财政、经济的发展状况，同时对掌握市场趋势和规划收支费用具有极为重要的意义。

国外对财政收入预测模型的研究，主要集中在对单一模型的选择和改进方面。Sexton (1987) [1] 使用经济计量模型、ARIMA模型、趋势模型三种方法对地方政府财政收入进行预测，认为前两种方法在预测效果上优于第三种。Mario等(2019) [2] 运用四个主要的欧元区经济体编制财政政策变量的季度数据集，建立四个主要欧元区国家的向量自回归模型，以此来预测各个国家的财政收入。

对于国内的相关研究，李伟(2011) [3] 运用神经网络与多元线性回归组合模型对国家财政收入进行了预测，得到了较为满意的结果。方博等(2015) [4] 采用ARMA-BP神经网络组合模型对国家财政收入进行预测，得到组合模型相较于单一模型有着更好预测效果的结论。赵海华(2016) [5] 则基于回归分析思想和相关算法，建立了一个多因素财政收入预测模型，该模型使用灰色系统理论结合RBF神经网络来实现数据预测。我国对单一省和地方性的财政收入预测研究相对较少。因此，本文针对福州市的财政收入进行预测研究。

综上所述，国内外学者较为常用的财政收入预测方法有神经网络、灰色理论和时间序列分析等。同时他们采用的单一预测模型得到的结果具有精度不高、泛化能力较差等局限性，组合模型未考虑到特征变量对财政收入的影响。由于随机森林(RF)在选择变量时具有较好的稳定性；GM(1,1)模型能够精准模拟指数序列且预测结果比较稳定；BP神经网络模型有着较好的容错性、很好的逼近特性和不俗的泛化能力；因此最终采用RF-GM(1,1)-BP组合模型对福州市财政收入进行科学的定量分析以及准确预测。

2. 算法模型

(一) 算法原理

随机森林(RF)算法是Breiman于2001年提出，是一种基于多棵树进行分类和预测的组合算法 [6] 。该算法通过自助抽样方法产生多个训练集，并使用这些训练集来构建多个分类回归树(CART)，作为基本的组件单元。在此过程中，Svetnik等 [7] 学者利用5折交叉验证计算训练集中变量的重要性得分，并对变量进行排序，重复20次后选取排名靠前的一半变量来构建RF模型，记录相应测试集上的错误率，并以最小平均误差来确定最佳变量的数目。随机森林模型选择变量的指标分为置换重要性和节点纯度增加的重要性，置换重要性 [8] 就是随机打乱变量 $\left(v=1,2,\cdots ,p\right)$ 在OOB中的值，如果变量重要，则打乱后增加，说明变量越重要。节点纯度增加的重要性则根据随机森林树模型中节点分裂前后的不纯度下降作为子节点纯度的增加，它们的公式如下：

$\begin{array}{l}VI{M}_v^{mse}=\frac{1}{ntree}{\displaystyle \underset{i=1}{\overset{ntree}{\sum }}\left\{MSE\left(OO{B^{\prime }}_{iv}\right)-MSE\left(OO{B}_{iv}\right)\right\}}\\ VI{M}_v^{purity}=\frac{1}{\left|T\right|}{\displaystyle \underset{X_V\in T}{\sum }\left\{i_v\left(t\right)-i_v\left(t_L\right)-i_v\left(t_R\right)\right\}}\end{array}$ (1)

其中 $OO{B}_{iv}$ 表示将第i棵树的第v个的变量OOB数据代入第i棵树，MSE表示平均平方误， $OO{B^{\prime }}_{iv}$ 的上标表示将OOB数据打乱， $i_v\left(t\right)$ 代表父节点的不纯度函数， $t_L$ 和 $t_R$ 分别代表父节点的左右子节点， $\left|T\right|$ 表示 $X_t$ 分裂变量的那些树的个数。建立随机森林回归模型选择变量，使用IncNodePurity作为重要性的判断指标，含义与节点纯度增加的重要性一致，IncNodePurity值越大，说明该变量越重要。

灰色系统理论 [9] 是我国学者邓聚龙教授于1982年提出的，适用于数据量有限、信息不完整时对不确定性问题进行研究。而灰色预测GM(1,1)模型是灰色系统理论的核心内容，其建模思想主要是通过已知的数据对系统信息的抽象和量化，推导出适合的模型，对未知数据进行预测。在该模型中，仅需少量样本数据即可建立数据模型，实现高精度的预测效果。灰色预测GM(1,1)模型是通过对原始数据进行一次累加后建立一阶线性微分方程，将得到的方程进行求解就可以得出预测模型：

${\hat{X}}^{\left(1\right)}\left(k+1\right)=\left[{\hat{X}}^{\left(1\right)}\left(0\right)-\frac{\hat{u}}{\hat{a}}\right]{\text{e}}^{-\hat{a}k}+\frac{\hat{u}}{\hat{a}}$ (2)

由于在建立微分方程时，采用的是一次累加后的数据，所以得到预测模型后，还需要进行累减还原得到原始序列的灰色预测GM(1,1)模型：

${\hat{X}}^{\left(0\right)}\left(k+1\right)=\left({\text{e}}^{-\hat{a}}-1\right)\left[{\hat{X}}^{\left(0\right)}\left(n\right)-\frac{\hat{u}}{\hat{a}}\right]{\text{e}}^{-\hat{a}k}$ (3)

BP神经网络是一种具有连续传递函数的前馈神经网络。BP神经网络的训练过程中使用误差反向传播算法来实现，最为常见的为梯度下降法 [10] 。通过最小化均方误差，不断地调整网络的权重和阈值以逐渐提高拟合数据的精度。这一过程可以帮助BP神经网络学习并优化输入数据之间的关系，从而更好地预测未知数据的结果。BP神经网络模型主要通过各层的权重、激活函数以及偏置来得到最终预测值，它将输入层的输出数据乘上一定的权重放入隐藏层中，作为隐藏层的起始值；在隐藏层中利用模型在该层中的偏置项和激活函数作用，得到新的值，然后将得到的新值乘上一定的权重放入输出层中，作为输出层的起始值；在输出层中同样利用模型在该层中的偏置项和激活函数作用，得到第一次迭代值。

隐藏层公式：

$\begin{array}{l}{q}_j={\displaystyle \underset{i=1}{\overset{5}{\sum }}{x}_i}\ast w_{ij}-a_j\\ h_j=f\left(q_j\right)\end{array}$ (4)

输出层公式：

$\begin{array}{l}u={\displaystyle \underset{j=1}{\overset{7}{\sum }}{\displaystyle \underset{i=1}{\overset{5}{\sum }}{h}_i}\ast v_j-\beta }\\ \hat{y}=f\left(u\right)\end{array}$ (5)

其中 $x_i$ 为输入数据， $w_{ij}$ 为输入层权重， $a_j$ 为隐藏层偏置， $f\left(\right)$ 为激活函数， $v_j$ 为隐藏层权重，为输出层偏置。接着，通过最小化损失函数，不断更新权重和偏置，得到最优的预测值，损失函数公式如下：

$L\left(e\right)=\frac{1}{n}{\displaystyle \underset{k=1}{\overset{p}{\sum }}{e}_k^2}=\frac{1}{n}{\displaystyle \underset{k=1}{\overset{p}{\sum }}{\left({\hat{y}}_k-{\stackrel{\dot{}}{y}}_k\right)}^2}$ (6)

其中 ${\hat{y}}_k$ 为经过神经网络后的输出值， ${\stackrel{\dot{}}{y}}_k$ 为期望输出结果即真实值。

(二) 评价标准

均方根误差(RMSE)是指预测值和真实值之间差距的平方和的平均值的开根号，可以用于评估BP模型预测结果的精度，能够很好地反映预测值与实际值之间的误差且受异常值的影响程度较小，因此使用RMSE来验证BP模型准确性，所得到的结果更加直观。它的计算公式如下：

$X_{RMSE}=\sqrt{\frac{{\displaystyle {\sum }_{i=1}^N{\left(X_{obs,i}-X_{model,i}\right)}^2}}{N}}$ (7)

R方(R-squared)是指模型能够解释因变量的变异程度，可以消除观察值尺度的影响，对BP模型性能的好坏有着准确的评价，R方越大说明模型拟合的越好，预测结果越精确。它的计算公式如下：

$R^2=\frac{SSR}{SST}=\frac{{\displaystyle \underset{i=1}{\overset{n}{\sum }}{\left({\hat{y}}_i-\bar{y}\right)}^2}}{{\displaystyle \underset{i=1}{\overset{n}{\sum }}{\left(y_i-\bar{y}\right)}^2}}=1-\frac{{\displaystyle \underset{i=1}{\overset{n}{\sum }}{\left(y_i-\hat{y}\right)}^2}}{{\displaystyle \underset{i=1}{\overset{n}{\sum }}{\left({\hat{y}}_i-\bar{y}\right)}^2}}$ (8)

(三) 研究框架

图1. 福州市财政收入预测研究流程图

福州市财政收入预测研究流程如图1所示，具体为：一、从《福州市统计年鉴》收集整理得到1994~2021年福州市财政收入的相关年数据；二、对数据进行描述性统计分析、相关性分析，再根据随机森林(RF)重要性选择变量；三、建立所选取变量的GM(1,1)模型，再对数据进行对数标准化，最后分别建立SVR和BP神经网络模型。四、根据均方根误差和R方两种评价标准的结果，得到最优模型；五、利用GM(1,1)-BP组合模型，得到福州市在2022~2025年财政收入的预测值。

3. 实证分析

(一) 数据介绍

预测福州市未来4年财政收入的相关数据均来自于《福州市统计年鉴》(1994~2022年)，数据真实可靠。

财政收入(y)与许多因素相关，选取了12个对财政收入影响较大的变量来预测未来4年的财政收入。由于1994年我国的财政体制进行了重大的改革，实施了分税制的财政政策。因此，采用1994年之后的数据进行分析。

社会总从业人数(x₁)：单位为万人，社会总从业人数对市财政收入有显著的影响，从业的人数越多，税收越多，居民消费水平越高，间接影响市财政收入。

在岗职工工资总额(x₂)：单位为万元，在岗职工工资总额越高，政府所收的税收越多，进而影响市财政收入。

社会消费品销售总额(x₃)：单位为亿元，社会消费品零售总额创造税收，而财政收入主要来自税收，所以社会消费品零售总额影响着市财政收入。

城镇居民人均可支配收入(x₄)：单位为元，城镇居民人均可支配收入提高，意味着居民的工作积极性提高，促进市财政收入快速增长。

城镇居民消费性支出(x₅)：单位为元，城镇居民消费性支出增加，意味着居民消费水平提升，政府可以从居民消费中获得更多税收，市财政收入也得到增加。

年末总人口(x₆)：单位为万人，在经济发展水平不变的情况下，一个地区的人口总数越多，则人均市财政收入越少。

地区生产总值(x₇)：单位为亿元，地区生产总值代表一个地区的经济发展水平，一个地区经济发展水平好，就会带动这个地区的财政收入，所以地区生产总值对市财政收入有着重大的影响。

第一产业产值(x₈)：单位为亿元，第一产业主要是农业方面，农业发展对市财政收入有着一定的影响。

居民消费价格指数(x₉)：将1990年的值视为100，得到的数据与之比较，单位为%，居民消费价格指数也叫做CPI，它反映了居民家庭所购买的消费品以及服务项目的价格变化的指标，影响着市财政收入。

第三产业与第二产业产值比(x₁₀)：单位为亿元，第三产业与第二产业产值比代表了产业结构的变动，产业结构得到优化，则市财政收入也会得到提升。

工业增加值(x₁₁)：单位为亿元，工业增加值反映了企业生产过程中新创造的价值，工业值变化，可以判断短期地区经济情况，与市财政收入息息相关。

建筑业增加值(x₁₂)：单位为万元，建筑业增加值是GDP的重要组成部分，反映了建筑业企业生产经营活动的最终成果，与市财政收入息息相关。

(二) 福州市财政收入相关变量描述性统计分析

对福州市财政收入相关变量进行描述性统计分析，得到福州市财政收入的最小值为209931.00万元，最大值为7498470.00万元，两者相差较大，说明福州市财政收入的变化比较明显，尤其是2010年之后，福州市财政收入大幅度增加；而通过福州市财政收入的标准差和均值，可以知道福州市各年财政收入数据较为离散。因此，通过建立模型来研究福州市财政未来的收入是有必要的。具体数据见表1。

表1. 福州市财政收入相关变量概括性描述

(三) 福州市财政收入各变量相关性分析

通过相关性分析得到各自变量与因变量间的相关关系，观察各自变量与财政收入的相关程度，对筛选变量做出初步判断。福州市财政收入各变量相关系数见表2，从中可以得到在岗职工工资总额(x₂)、社会消费品销售总额(x₃)、城镇居民人均可支配收入(x₄)、城镇居民消费性支出(x₅)、第一产业产值(x₈)、工业增加值(x₁₁)变量与财政收入(y)间均具有较高的相关关系且相关性最高的达到了0.99，第三产业与第二产业产值比(x₁₀)变量与财政收入(y)间具有最低的相关关系达到了0.70。另外，其他自变量间的相关程度也比较高，说明各变量间存在多重共线性。

表2. 福州市财政收入各变量相关系数表

4. RF-GM(1,1)-BP模型的建立

(一) RF选择变量

由于各变量间具有多重共线性，因此采用随机森林(RF)模型，筛选出对财政收入来说最重要的几个变量，确保在保留数据特征的前提下，减少数据的维度，使模型预测结果更加精确。通过5次重复十折交叉验证得到：当自变量保留4~7个时，所建立的模型最优，通过分析5次重复十折交叉验证法的结果并结合所学知识，决定保留对预测财政收入最重要的5个变量，将自变量按照下标顺序排序输出，具体结果见表3。

表3. 重要性原则选择变量表

表3中，%IncMSE就是increase in MSE，指使用随机值代替各个变量，变量替换后模型的误差，而IncNodePurity就是increase in node purity，指每个变量对回归树的每个节点观测值的影响，两者值越大说明该变量越重要，都可以作为重要性的判断指标，采用IncNodePurity来选择变量。在经济意义方面，社会总从业人数(x₁)、在岗职工工资总额(x₂)、社会消费品销售总额(x₃)、城镇居民人均可支配收入(x₄)、城镇居民消费性支出(x₅)这5个变量对财政收入的影响均反映在税收上，因此，只需保留一个变量；由于2006年国家取消了农业税，导致第一产业产值(x₈)对财政收入影响较小，因此剔除该变量；居民消费价格指数(x₉)影响了居民的消费额，对财政收入的影响也反映在税收上，所以剔除该变量；第三产业与第二产业产值比(x₁₀)代表了产业结构的变动，对财政收入有一定影响，但相较于其余变量权重较小，因此剔除该变量。

综上所述，结合随机森林重要程度识别出5个影响福州市财政收入的变量分别为：社会消费品销售总额(x₃)、年末总人口(x₆)、地区生产总值(x₇)、工业增加值(x₁₁)、建筑业增加值(x₁₂)。

(二) GM(1,1)-BP模型的建立

未来的财政收入情况关系到政府工作规划，需要十分精确的结果。本文所得到的数据单位和大小都存在一定的差距，直接进行建模预测，得到的结果误差极大，不具备可信度，而对于经济数据来说，最适合的标准化方法为对数标准化 [11] 。为了得到精确的预测结果，对数据进行对数标准化，消除数据间大小和单位的影响。再建立BP神经网络模型。

第一步，将数据划分为训练集和测试集，利用测试集检验模型优劣，防止出现过拟合现象，提高模型的泛化能力。由于数据为时间序列，具有连续性，数据的划分不能采用随机形式，选取1994~2012年福州市财政收入相关数据作为训练集，2013~2021年福州市财政收入相关数据作为测试集。

第二步，根据得到的训练集建立BP神经网络模型，神经网络具有多种不同的模型，采用三层神经网络，该模型共分为三层，分别为输入层、隐藏层和输出层。再将BP神经网络中最大迭代次数设置为10,000，误差精度设置为0.00001；同时，由于适合的隐藏层节点数会大大增加模型的精确度。因此，通过调参优化的方式，将隐藏层神经元个数设置为7，建立最优的BP神经网络模型。

第三步，将建立的支持向量机回归模型与BP神经网络模型相比较，再构建对比模型SVR。

(三) BP模型检验及优化

将划分好的测试集分别代入SVR和BP神经网络中，利用R中的函数计算两个模型的RMSE (均方根误差)和R-squared (R方)，根据得到的结果判断模型的优劣，选择出最适合的模型，检验结果见表4。

表4. 模型检验及对比表

从表4得到的结果来看：RF-BP模型的RMSE (均方根误差)为0.0354，R-squared (R方)为0.9497；RF-SVR模型的RMSE (均方根误差)为0.4903，R-squared (R方)为0.9022。均方根误差(RMSE)是指预测值和真实值之间差距的平方和的平均值的开根号，用来评价模型的准确性；R方(R-squared)是指模型能够解释因变量的变异的程度，用来评价模型的拟合程度，均方根误差越小、R方越大说明模型拟合的越好，预测结果越精确。因此，使用RF-BP模型对福州市财政收入进行预测，会得到精确的结果。

5. 财政收入预测

(一) GM(1,1)预测自变量的值

GM(1,1)模型有着几个优点：1) 具有时效性。对于变化较快的数据，可以更好地满足预测需求；2)适用性强。对于不规则时间序列的预测能力强，可以更准确地反映数据的发展规律；3) 对异常值鲁棒性强。即使数据中存在异常情况，也可以保证预测的准确性；4) 中短期预测准确。对于中短期的预测更加精确。而福州市财政收入各变量的数据没有明显的规律，预测未来4年的值也属于短期预测。因此，与其他模型相比，在预测未来4年福州市社会消费品零售总额、年末总人口、地区生产总值、工业增加值以及建筑业增加值的情况下，采用GM(1,1)模型具有更大的优势。

为预测未来4年福州市财政收入总值，利用GM(1,1)模型对通过随机森林法选出的5个关键变量进行预测，得到影响福州市财政收入各变量未来4年的值。

首先，将社会消费品销售总额(x₃)、年末总人口(x₆)、地区生产总值(x₇)、工业增加值(x₁₁)、建筑业增加值(x₁₂)这5个特征变量分别进行一次累加生成一次累加序列，建立出5个灰色微分方程进行求解，运用后减运算还原得模型输入序列预测序列，具体结果见表5。

表5. 相关因素预测结果及检验

通过GM(1,1)模型得到2025年社会消费品销售总额(x₃)的预测值为10892.05亿元，年末总人口(x₆)的预测值为748.75万人，地区生产总值(x₇)的预测值为20984.12亿元，工业增加值(x₁₁)的预测值为5270.40亿元，建筑业增加值(x₁₂)的预测值为3601.48万元。接着，比较后验差比值(C)以及预测图来检验模型预测精度。从表5中可以得到社会消费品销售总额、年末总人口、地区生产总值、工业增加值、建筑业增加值的后验差比值(C)均小于0.35，说明模型精度等级为好，模型预测结果可靠，可用于预测福州市财政收入。

(二) RF-GM(1,1)-BP模型预测福州市未来4年财政收入

将自变量在福州市未来4年的预测值代入建立隐藏层节点为7、最大迭代次数为10,000、误差精度为0.00001的BP神经网络模型中，得到福州市未来4年财政收入的预测值，具体结果见表6。

表6. 福州市财政收入未来4年预测结果

从表6中，可以得到2022年福州市财政收入预测值为7,959,923万元，2023年福州市财政收入预测值为801,524万元，2024年福州市财政收入预测值为8,025,953万元，2025年福州市财政收入预测值为7,999,256万元，可以看到福州市财政收入将在2021年之后稳步增长，表明福州市经济发展状态良好。

6. 结论与建议

(一) 结论

1) 福州市财政收入将在2021年之后稳步增长，表明福州市针对经济发展的政策十分优秀。

2) 工业增加值对福州市财政收入而言最为重要，年末总人口数次之，然后是地区生产总额，社会消费品销售总额，建筑业增加值。政府的财政收入大多来自税收，工业增加值代表着企业创造的价值，说明企业创造的价值在税收收入中所占的比例较大；年末总人口数对财政收入起着积极作用，说明福州市人民的收入和人民生活水平较高。

3) 社会消费品零售总额在2022、2023、2024、2025年预测值为7666.86亿元、8618.83亿元、9689.00亿元、10892.05亿元；年末总人口数在2022、2023、2024、2025年预测值为727.40万人、734.45万人、741.57万人、748.75万人；地区生产总额在2022、2023、2024、2025年预测值为14654.55亿元、16517.53亿元、18617.35亿元、20984.12亿元；工业增加值在2022、2023、2024、2025年预测值为3962.96亿元、4358.08亿元、4792.58亿元、5270.40亿元；建筑业增加值在2022、2023、2024、2025年预测值为2268.40万元、2646.31万元、3087.17万元、3601.48万元。

4) 福州市财政收入总值在2022、2023、2024、2025年的预测值为7,959,923万元、8,015,241万元、8,025,953万元、7,999,256万元，表明福州市财政收入将在2021年后不断增长。在疫情结束后，各个地区都在加快恢复经济建设，福州市在未来的财政收入有极大的增长，福州市政府的治理能力和管理水平较高。

(二) 建议

1) 福州市针对经济发展的政策还有可以提升的地方，如改善财政支出结构等。

2) 企业创造的价值在税收收入中所占的比例较大，政府应该加强对地方产业政策的制定和实施，引导产业结构调整和转型升级，提高地方财政收入；提高政策影响，提高福州市人民生活水平，保持年末总人口数对财政收入增长的积极作用。

3) 政府可以加大社会消费品零售、工业和建筑业方面的投入，并通过政策鼓励，增加福州市人口总数，提高福州市地区生产总值，来增加地区的财政收入。

4) 地区财政收入增长较多，并不一定代表地区真正的繁荣，相关部门应该更加注意，避免出现资源的过度消耗以及税收的过度增加。

基金项目

项目来源：福建省科技厅。

项目名称：基于ATOT技术的智能养老系统设计与开发。

项目编号：2023350104000282。

文章引用

刘 威,张 巧,王文博,许 可,张慧妍,金秀玲. 基于RF-GM(1,1)-BP模型预测福州市财政收入
Predicting Fuzhou City’s Fiscal Revenue Based on RF-GM(1,1)-BP Model[J]. 应用数学进展, 2023, 12(07): 3240-3249. https://doi.org/10.12677/AAM.2023.127323

参考文献