李林红 王娟 马亚梅

摘要：云南高原山地农业的边缘性和分散性使得云南山地农业系统在很大程度上陷入了一种孤立封闭的自给自足状态。传统小农户经营生产方式与社会化大市场存在巨大鸿沟，信息不对称导致市场失灵，单靠高原山地农户自身很难实现可持续发展，需要公众、专家、政府共同参与规划，三者达成共识是制定可持续发展规划的前提。以云南高原山地农业可持续发展为研究对象，结合网络上公众、专家和政府观点的响应关系建立群体认知系统动力学模型，模拟三者对云南高原山地农业可持续发展的认知共识过程。借鉴种群生态学Logistic模型对模型的稳定性进行分析，并用Matlab对稳定性分析结果进行模拟仿真，进一步验证模型的稳定性，最后给出实现云南高原山地可持续发展的政策建议。

关键词：可持续发展;群体认知;系统动力学;稳定性分析

中图分类号： F327  文献标志码： A  文章编号：1002-1302（2020）03-0311-06

高原山地农业耕地呈块多、面小、分布零散，距离遥远、交通设施落后、环境封闭、流动性有限，导致高原山地农业依然处于自给自足的传统农业阶段，农村经济发展滞后、农业资源短缺、抗灾能力低、农业环境污染等问题更加凸显[1]。山地农业的脆弱性，使得各种以山地资源为开发对象的物质生产活动极易导致水土流失等环境灾害[2]，探寻如何在保护生态环境前提下发展山地农业，一直是难题[3]。可持续发展的核心是满足当代人的需求，又不损害子孙后代的利益满足其需求能力，实现永远稳定生存与发展[4]。然而，由于存在市场失灵，单靠高原山地农户自身很难实现可持续发展，需要公众（包括农户和消费者）、专家、政府共同参与规划，三者达成共识是制定可持续发展规划的前提。

系统的稳定性是系统分析和设计的基本问题[5]，是实现可持续发展的理论依据，只有当所建立的非线性动力学模型处于稳定点或稳定区间时，云南高原山地农业才可实现长期可持续发展。高原山地农业是一个极其复杂的生态生产系统，因其所处偏远的地理位置和近年来日益严峻的农业环境污染问题，使得高原山地农业系统不断受到外界的干扰和挑战，只有充分认识系统的承受能力和稳定区间，才可以更好地在发展经济的同时保证长期可持续发展。

国内外已有许多专家致力于研究非线性动力系统的稳定性，现有关于解决非线性动力学系统稳定性问题的方法的研究可分为两大类：（1）基于数学模型和知识推理的客观分析方法[6]，其中最常用的是Logistic模型和Lyapunov-Krasovskii模型，Logistic模型源于生物种群共生理论，主要解决多要素在竞争和协同共同作用下实现系统长期稳定性的问题[7];Lyapunov-Krasovskii泛函模型可以研究复杂时滞切换系统，但包含时间平移或者状态导数，容易导致稳定性分析失去因果性[8]。为了将稳定性分析与复杂的农业生产系统紧密联系起来，并避免政策结论推广中的局限性，选择Logistic模型作為稳定性分析的方法。（2）基于仿真模拟来分析非线性动力学系统的稳定性。传染病动力系统模型对于求解非线性系统稳态解的条件具有明显优势，但是模型中会涉及与时间有关的复杂速率函数和过程函数[9]，不便于本研究进行仿真模拟。Vensim PLE仿真模拟软件对于证明系统存在稳定性和寻找系统稳定点十分便捷，但是却不便于准确求解复杂系统的稳定区间[10]，本研究运用Matlab软件对非线性代数方程组进行求解和稳定区间的模拟，系统地分析高原农业可持续发展的稳定区间，为促进高原农业长期可持续发展提供见解和指导。

1 构建可持续发展网络群体认知系统动力学模型

可持续发展网络群体认知是基于互联网的公众、专家和政府人员对可持续发展的认知，从系统动力学的角度而言，公众、专家和政府之间的认知相互影响和相互关联，可用系统动力学（SD）方法予以分析研究。因此，选取网络上的公众、专家和政府人员为认识主体，高原农业可持续发展问题为认知对象，运用SD方法对公众、专家和政府群体的共识过程进行模拟。

1.1 构建模型的网络图

根据群体认知和系统动力学理论构建“广义专家”群体研讨发言的基本结构单元见图1。

“广义专家”指公众、专家和相关政府人员三大群体，通过网络参与主题研讨活动。“参与者”是指参与群体研讨发言的相关人员，专家采用真实姓名，政府人员采用政府部门名称，公众采用网络的国际互联协议（interhet protocol，IP）。“信息名”指每条发言信的编号，“响应对象”指诱发本条发言信息出现的另一条已存在的发言信息，指向其“信息名”。“响应属性”反映发言信息与其响应对象之间的关系属性，是支持、反对、中立、质疑、补充等属性之一。“关键词”指发言信息中出现的高频、关键性词语[11]。收集公众、专家和政府人员对如何实现高原山地农业可持续发展的看法，依据图1整理每条发言的“参与人”“信息名”“响应对象”“响应属性”“关键词”等信息，建立“广义专家”群体研讨发言信息资料（表1）。其中Ax、Bx和Cx分别代表公众、专家和政府人员发言的编号。

以每次发言为1个节点，“序号”和“信息名”标识该节点，以“响应属性”作为连接不同节点的有向弧，指向其“响应对象”。建立“广义专家”群体研讨互动过程的有向链接属性图，即模型的网络见图2。

1.2 构建模型的因果关系

从表1可以看出，图2模型的网络图中节点与响应对象节点之间的响应属性有质疑、反对、中立、支持、补充5种。当响应属性为中立时，二者不构成因果关系，其响应关系忽略不计;当响应属性为质疑或反对时，二者构成负向因果关系，用负性因果键标识;当响应属性为支持或补充时，二者构成正向因果关系，用正性因果键标识。根据以上原则构建模型的因果关系，结果见图3。

1.3 构建模型的反馈结构

根据图3构建模型的反馈结构见图4，图3中因果键指向的节点在反馈图中转化为存量，无因果键指向的节点转化为辅助变量。图4中矩形表示存量，其值表示发言信息的观点强度，每个存量都有输入流和输出流，输入量表示存量所代表观点强度的增加量，输出量表示存量所代表观点的减少量，输入流和输出流由因果关系键的指向节点所转化的存量和自身存量共同决定，用流线连接。图中每个箭头连接的2个变量间的数学关系经人为设定就能进行仿真模拟[12]。该模型共包括25个变量，其中存量8个，流量16个，辅助变量1个。模型的反馈结构图中各参数的含义见表2。

1.4 模型仿真结果

假设0～1表示观点的反对强度，越靠近0，观点的反对强度越强，0表示完全反对;1～2表示观点的赞同强度，越靠近2，观点的赞同强度越强，2表示完全赞同，1表示观点持中立状态。

图4中各存量的初始值设为1，即在初始状态下存量所代表观点持中立态度，取值范围为[0，2]。图4中观点所代表的存量的流入量和流出量由响应该观点的存量决定，如B5的输入流IN 1和输出流OUT 1由A9决定。A9和B5是正向因果关系，当A9>1时，IN 1等于（A9-1）×B5;当A9     运用系统动力学软件Stella进行仿真模拟，得出各个观点随时间的变化观点强度的变化情况见表3。

2 模型稳定性分析

从表3可以看出，公众、专家和政府三大群体对高原山地农业可持续发展问题的看法所持观点强度变化类似于种群进化Logistic规律变化[13-14]，本研究借助种群生态学的Logistic模型对可持续发展网络群体认知系统动力学模型进行抽象建立微分方程组，找出模型的均衡点进行稳定性分析。

2.1 建立微分方程

将模型中的存量看作各个独立的种群，它们相互之间存在于网络研讨环境中，存量所代表观点的强度变化服从Logistic模型规律。根据Logistic模型作如下假设[13-14]：假设1，辅助变量和存量所代表观点强度是关于时间t的函数，用xi（t）表示。假设2，在独立存在的状态下，辅助变量和存量所代表观点强度存在1个最大值，用Mi表示。假设3，在独立存在的状态下，辅助变量和存量所代表观点强度的平均增长率，用ri表示。假设4，辅助变量和存量所代表观点强度对响应对象所代表观点强度的影响，用φi表示。

2.3 稳定性结果仿真

根据上述稳定性分析结果，在Matlab中将Mi值和求得的φi稳定区间带入（1）～（9）微分方程，得到系统動力学模型的稳定性结果仿真图（图5）。观察各个辅助变量和存量xi（t）随时间t的变化趋势，发现在现有的φi稳定区间内，A1、A2、A3、B1、B3、B4、B5、C1都在第12个时间节点之后逐渐趋于稳定，仿真结果较好地验证了上述稳定性分析过程。

3 结论

根据模拟结果可以看出，所构建的可持续发展网络群体认知模型是合理的、有效的，可以较好地表达公众、专家和政府3个群体的相互作用。模拟结果与高原山地农业的发展现状是比较符合的。通过模型的构建和稳定性分析发现，在竞争和协同作用的双重作用下，复杂系统中群体通过相互涌现使得观点的平衡状态保持长期稳定。因此，为了促进云南高原山区农业的可持续发展，结合本研究结论提出以下政策建议。

（1）促进高原山区农业的可持续发展不是一个可以通过扩大资本投资渠道和增加政府支持来解决的问题。有必要找出更深入、更有效解决方案。近年来，国家不断加大对于高原农业发展的资金扶持力度， 但是农业发展资金的大量投入并没有收到

显著的回报，当前政府提出的为提高农民收入、促进农业可持续发展而呈现出投资大、回报小的状况得到了政府、专家、公众的一致赞同，这与当前高原农业发展的现状是相符的。

（2）公众、专家和政府人员意识到，实现高原山区农业的可持续发展不仅取决于政府投资，还取决于全民参与。因此，政府应该在促进高原山区农业可持续发展的过程中听取高原山区农民的意见，确保科学、正确的决策。促进高原农业可持续发展不能一味地进行资金投入，也不能只依靠政府干预，群众才是从事农业生产的主力军，只有广泛听取专家的意见和群众在农业生产中遇到的问题，从而切实地解决和落实，才能使得政府政策发挥更大的效用。

（3）公众、专家和政府人员希望将支持资金用于合理的地方，防止支持资金的使用不合理，降低支持高原山区农业可持续发展的有效性。因此，在支持可持续发展的过程中对高原山区农业，政府应加大对支持资金的监管力度，完善相应的扶持资金审查制度，在促进高原农业可持续发展的过程中接受群众的监督。

参考文献：

[1]杨庆华，杨世先，马文彬，等. 山地农业区耕地资源空间分布状况——以云南省玉溪市为例[J]. 山地学报，2005，23（6）：749-755.

[2]程厚思，卢亚鹏. 云南山地农业的环境特征及其可持续发展[J]. 云南社会科学，1999（6）：1-7.

[3]富兰克林·H·金. 四千年农夫[M]. 北京：东方出版社，2011.

[4]李林红. 滇池流域可持续发展投入产出系统动力学模型[J]. 系统工程理论与实践，2002（8）：89-94.

[5]闫红梅，张 鸣，李远征. 基于Matlab的系统稳定性分析实验设计[J]. 实验技术与管理，2018，35（4）：144-146，150.

[6]Tai L J，Li L H，Du J. Evaluation of scale effect of fragmented agricultural land transfer based on neural network[J]. NeuroQuantology，2018，16（5）：839-846.

[7]阮平南，张敬文. 战略网络的自组织机制及稳定性分析[J]. 科技进步与对策，2008，25（12）：29-31.

[8]丛 屾，费树岷，李 涛. 时滞切换系统指数稳定性分析：Lyapunov-Krasovskii泛函方法[J]. 系统工程理论与实践，2007（7）：111-115.

[9]李先利，胡顺菊. 一类非线性动力学系统的稳定性分析[J]. 系统工程理论与实践，1998（2）：66-71.

[10]邢海龙，高长元，张树臣. 基于系统动力学的大数据联盟稳定性模型构建与仿真研究[J]. 情报杂志，2017，36（10）：159-165.

[11]马亚梅，李林红. 基于AHP的网络公众对可持续发展认知的研究——以昆明信息港“彩龙论坛”为例[J]. 昆明理工大学学报（社会科学版），2010，10（4）：55-58，64.

[12]李金兵，韩玉启，罗建强. 基于系统动力学的企业复杂系统资源负熵管理模型仿真[J]. 系统管理学报，2011，20（5）：600-605.

[13]胡静静，严广乐. 从生态学分析中小企业与大型企业自愿性信息披露[J]. 商业时代，2009（5）：93-94.

[14]胡静静，严广乐. 上市公司自愿性信息披露动因机制的生态学角度研究[J]. 科技与管理，2009，11（1）：25-27，30.