程定性理论与系统动力学模型分析方法.pdfVIP

微分方程定性理论与系统动力学 模型分析方法 刘宇熹1，杨波2 (1-上海理工大学管理学院上海200093； 2．上海交通大学安泰管理学院上海200052) 摘要应用微分方程定性理论结合系统动力学结构功能模拟方法对具体例子 进行了初步的实证研究，提出了几点思考，并给出了具体算法。作为应用的前提， 文中还解决了系统动力学模型对应微分方程组中离散表函数的连续化问题。 关键词 系统动力学；微分方程；表函数 0引言 率流位系分析研究复杂信息反馈系统的科学，以计算机为工具，通过结构功能模拟的方法 揭示复杂系统的行为特征，在诸多领域得到较好和广泛的应用。大量应用实践表明，系统动 力学建模、仿真和调试过程不是一个容易的过程。系统模型的行为特征结果与模型的参数 密切相关。如何有效地对大量参数进行控制和调节是困扰建模与分析人员的重要问题．这 方面一直以来缺乏有力的手段。由于系统动力学模型与微分方程组的对应关系．从微分方 程理论的角度进行系统动力学模型的分析也是一个受关注的议题，微分方程作为数学的·}t 心学科。经过约300年的发展，其理论和方法都El臻完善。微分方程理论在系统动力学中的 应用过去主要集中于对SD模型对应微分方程组的解的研究，但前提是要求出方程组的一 般解。事实上，随着系统动力学在众多复杂实际系统中的应用，其系统刻画模型对应的微分 方程组大多为非线性微分方程，由于非线性方程原则上已求不出解来，使得微分方程的传统 理论在SD模型分析中似乎失去价值，需要借鉴该领域中新的理论与方法。这就是本文应 用的微分方程定性理论。微分方程定性理论是1885年由法国著名数学力学家庞加莱创立 的，它是现代微分方程的主导分支，在揭示非线性科学中大量深刻而有趣的新现象的数学机 制上作用重大。借助这一理论对复杂的系统动力学模型进行深人分析，可以为参数的有效 调节和控制提供新的指导。 1离散值表函数的连续化问题 文献[1]中已给出了SD模型的微分方程组形式，见如下定理。 ：作者简介]刘宁矗(】977、．男．硕士，助教一研究方向为决簸支持、管理信皂系统。联系电话：02I一荫68185 Enmil：lv^yb@rolll．m，扬渡(1q7‘卜)．女．博士研究牛．研究方向为管理系统上程、复杂系统。 954 (i--1，2，…一)，则在此确定的流位流率系下，系统的微分方程组模型为： 堕掣一R。(}) dt ，，1(f)It-t．一Ll(如) ： (1) 掣一Rm(r) L。。(￡)I一一L。：(f) 以式(1)中第1个方程为例，右边各流位变量、外生变量和常量的组合式有解析表达式， 还可能是表函数形式，表函数提供的是离散点值，因此^此时不一定连续，这阻碍了微分方 程理论在SD中的应用，事实上，以前的相关文献忽略了此方面的问题。本文认为，严谨性 是进行科学研究的必备态度，因此有必要在此方面作一点说明。由于SD表函数提供了未 知函数关系，(z)在”+1个点上的值，(∞)，自然想到利用数值逼近中的捕值方法寻找某 数插值方法中的拉格朗日插值公式或者Newton插值公式，即逼近函数可取为： 一(，)一∑lⅡ(三≥)lY， (2) 120。i20、“J¨。 ，≠J 或者