城市系统非线性动力学分析的基本理论与方法.docVIP

城市系统非线性动力学分析的基本理论与方法 信阳师范学院(自然科学版) 第11卷第4期1996年lO月 JomalofrTeach~c口l (N曲1ralscjerHEdi0n)Ⅵ.儿№.4Oct.1998 城市系统非线性动力学分析的基本理论与方法 单纬东陈彦光 (信阳师范学院地理系.河南,信阳46舢) 摘耍奉文提出城市系统非裁性动力学分析的基本思路:1,一般稳定性分析;2,分数堆分 析,即支起相空闻,计算关联堆,确定系统控制变量的数目;3,模拟分析,即找出主导变量,重建城 市曲力乐蟪,基此借助计算札模拟分析城市系境的演化机嘲. 关■词城市蛄构;城市体系;滩沌;分彤 分粪号I;一06 城市系统(本文就广义而言,包括citys血厦 systemdd曲)具有.先验的分形(a咖1)几何特征 (阵彦光,199S)…,而分形结构与混沌(d蚺)过程内 在相美:混沌是决定论系统中的内在随机性,分形在 某种意义上乃是随机现象中的确定型规剜,二者似乎 是一对.共轭性问题.因此可以判断,城市系统的演 化具有棍沌性质,必须借助非线性理论与方法才能揭 示其动力学机制(李后强,艾南山,1996)Ⅲ.本文试从 一 般城市系统的稳定性分析出发,提出通过分散维 (岛.nd1.mio.)分析,重建城市动力系统,基此模 拟其演化机制的基本思路. 1基本理论 1.1城市系统的数学描述 . 设城市是由种元素组成的非线性系统,剐其自 治方程可写成如下初值同题 r也..,,. J…,… L(b)=0),(i=1…2..,H) 很定系统演化处于静态状况,即弓}=0,此时=五= …==o.由此得H个变量为n的方程,其定态解为 屯=知(i=1,2,…,),设{自}附近的小扰动为, 郾有=.并且ll《1,则系统(1)在定态 I}附近经Tl0r级数展开后可线性化为如下方程 组 警=喜c叶 可以证明:如果方程组(2)的零解(毗=0)是渐近稳定 的,剜式(1)的定态解{轴】也是新近稳定的;如果式 (2)的零解不稳定,剜式(1)的定态解{1也不稳 定.这个结论是应用线性分析方法分析城市演化 非线性美系的重要依据之一. 1.2城市系统的稳定性分析 为了简化同题,只考虑两个状态变量(如群体中 的两个城市.或个体城市中的两种元素)的非线性动 力系统.式(1)化为 f警(3) 【鲁=(t, 如果存在一点(知,知)使得(,)=o且(, 知)=0,划此点非正则,人们称之为奇点,它是相空同 中的轨道变汇点,实划系统的定态解,故可借助奇点 附近的轨道分布图式分析系统的稳定性,为此令, *国家自然科学基4#厦河南省自然科学基金(984町10o0)资助项目 收稿日期199s-_o6—09 ● 第ll卷第4期单纬东等:城市系统非线性动力学分析的基奉理论与方法379 为(m,和)附近的小扰动.即设施=+,而= 和叶,且有l等l《1,ll《1,则根据式(2)井考虑 到奇点的性质有 『警=十叶1{) 叶 式中=(差)0B-=(罄,(甏:(罄 方程(4)的解可表示为=,=,代人上 式,化为(,)的一阶齐欢代敷方程组,表为矩阵 形式如下. 【==】=【::㈣ 可以看出.式(5)具有非零解的条件为 三一】=.㈤【一】 展开得特征方程(一^)(一)一.唧=0,即 ^.一a+D=0(7) 式中=,D=口j一吩:特征方翟可有两 个不同的特征根^,,即 ^=告±√()一D={(c±/)(8)●■● 显然E=一4D.研究表明:当c《0,Dgt;0.E)0 时,特征方程有两个负解,由此得到一个结点,当e一 =0,第二项以后各项依欢变小,系统趋向于静止状态 (图1一口); 当c《0,Dgt;0,E(0时,特征方程的两个解为带 有负实部的复敷,由此得到一个环,点(t)沿着曩 旋线趋向于定志(,如)(图1—6); 当=0,口gt;0,《0耐,两个解都是虚数.因此 方程组包含周期项,围绕静止点摆动或麓转,点(而. ≈)围绕(,却)画出一条封闭曲线(图1一c); 在cgt;0,D)0,Egt;0时.两个解都是正敷,不存 在静态(L_youBem曲,1968).[43此时系统可瞻会进 人混沌状态.在前三种情况下,系统的相图中出现的 都是较为简单的吸引子,如不动点或授限环等(图1): 系统进人混沌志以后则会出现奇异吸引子(蛐学砒. 血ac.而城市,在某种意义上,其自身便是空间的 混沌吸引子(D.s.I)endfin~,1996)Is3研究城市系统 的混沌过程.仅作稳定性分析是不够的,必须重建动 力系统,揭示其时空结构的分敷雏特征. 王 国1城市系统非线性相关作用的时空围式(相固 2分维分析 2,l建立城市系统的相空间 如果城市系统的滇化过程是混j电的,则可计算其 关联维数,然后进行分析.首先须将状态变量转换成 时问序列,以一种