美赛常用模型一副本PPT学习教案

1、会计学1 美赛常用模型一美赛常用模型一 副本副本 第1页/共43页 一个雨天，你有件急事需要从家中到学校去，学校离 家不远，仅一公里，况且事情紧急，你来不及花时间 去翻找雨具，决定碰一下运气，顶着雨去学校。假设 刚刚出发雨就大了，但你不打算再回去了，一路上， 你将被大雨淋湿。一个似乎很简单的事情是你应该在 雨中尽可能地快走，以减少雨淋的时间。这是最好的 策略吗？试建立数学模型来探讨如何在雨中行走才能 减少淋雨的程度。 例1 雨中行走 第2页/共43页 1 1 建模准备建模准备 建模目标：在给定的降雨条件下，设计一个雨中行走的策建模目标：在给定的降雨条件下，设计一个雨中行走的策 略，使得你被雨水

2、淋湿的程度最小。略，使得你被雨水淋湿的程度最小。 主要影响因素：主要影响因素： 淋雨量，淋雨量， 降雨的大小，降雨的方向（风），路程的远近，降雨的大小，降雨的方向（风），路程的远近， 行走的速度行走的速度 2）降雨大小用降雨强度 厘米/时来描述，降雨强度指单位 时间平面上的降下水的厚度。在这里可视其为一常量。 3）风速保持不变。 4）你一定常的速度 米/秒跑完全程 米。 h 2 模型假设及符号说明 1）把人体视为长方体，身高 米，宽度 米，厚度 米。 淋雨总量用 升来记。 w d C I v D 第3页/共43页 3 3 模型建立与计算模型建立与计算 1 1）不考虑雨的方向，此时，你的前后左右

3、和上方都将淋雨。）不考虑雨的方向，此时，你的前后左右和上方都将淋雨。 淋雨的面积 )( 22 2 米wddhwhS 雨中行走的时间 )(秒 v D t 降雨强度 )/()3600/01. 0()/(01. 0)/(smIII时米时厘米 （升）米SIvDSItC3600/)/(10)(01. 0)3600/( 3 模型中 为变量。为参数，而vSID, 结论，淋雨量与速度成反比。这也验证了尽可能快跑能 减少淋雨量。 第4页/共43页 。米即米米米 小时厘米米若取参数 2 2 . 2,20. 0,50. 0,50. 1 ,/2,1000 Sdwh ID 秒。分秒，即你在雨中行走了 每秒，则计算得米度

4、你在雨中行走的最大速 472167 /6v 从而可以计算被淋的雨水的总量为2.041（升）。 经仔细分析，可知你在雨中只跑了2分47 秒，但被淋了 2 升的雨水，大约有4 酒瓶的水量。这是不可思议的。 表明：用此模型描述雨中行走的淋雨量不符合实际。 原因：不考虑降雨的方向的假设， 使问题过于简化。 第5页/共43页 2）考虑降雨方向。 v h w d 人前进的方向 若记雨滴下落速度为 （米/秒） r 雨滴的密度为 雨滴下落的反方向 1 ,pp 表示在一定的时刻 在单位体积的空间 内，由雨滴所占的 空间的比例数，也 称为降雨强度系数。 所以， rpI 因为考虑了降雨的方向，淋湿的部位只有顶部和

5、前面。分两部分计算淋雨量。 第6页/共43页 顶部的淋雨量 )sin()/( 1 prwdvDC 度。表示雨滴垂直下落的速 表示顶部面积，表示在雨中行走的时间 sin ,/ r wdvD 前表面淋雨量 )cos()/( 2 vrpwhvDC 总淋雨量（基本模型） )cos(sin( 21 vrhdr v pwD CCC 6 1039. 1,/23600,/4 pscmIsmr取参数 第7页/共43页 )5 . 1cos6sin8 . 0( 1095. 6 4 v v C 可以看出：淋雨量与降雨的方向和行走的速度有关。 问题转化为给定 ，如何选择 使得 最小。v C 情形1 90 )5 . 1

6、8 . 0 (1095. 6 4 v C 结果表明：淋雨量是速度的减函数，当速度尽可能大时 淋雨量达到最小。 假设你以6米/秒的速度在雨中猛跑，则计算得 第8页/共43页 升13. 1103 .11 34 mC 情形2 60 / )334 . 0(5 . 1 1095. 6 4 vC 结果表明：淋雨量是速度的减函数，当速度尽可能大时 淋雨量达到最小。 假设你以6米/秒的速度在雨中猛跑，则计算得 升47. 1107 .14 34 mC 情形3 18090 此时，雨滴将从后面向你身上落下。 5 . 1/ )cos6sin8 . 0(1095. 6 4 vC 第9页/共43页 5 . 1/)90co

7、s(6)90sin(8 . 0(1095. 6 4 vC 5 . 1/ )sin6cos8 . 0(1095. 6 4 vC 能的。可能取负值，这是不可时，当C 900 出现这个矛盾的原因：我们给出的基本模型是针对雨从出现这个矛盾的原因：我们给出的基本模型是针对雨从 你的前面落到身上情形。你的前面落到身上情形。 因此，对于这种情况要另行讨论。因此，对于这种情况要另行讨论。 当行走速度慢于雨滴的水平运动速度，即 sinrv 这时，雨滴将淋在背上，而淋在背上的雨水量是 vvrpwDh/ )sin( 淋雨总量为 vvrhdrpwDC/)sin(cos 。，则令 90090 第10页/共43页 取到最

8、小值。时，当Crvsin cos sin wdpr r D C 再次代如数据，得 )sin4/()cos8 . 0(1095. 6 4 C 结果表明：当行走速度等于雨滴下落的水平速度时，淋 雨量最小，仅仅被头顶上的雨水淋湿了。 若雨滴是以 的角度落下，即雨滴以 的角从 背后落下，你应该以 120 30 的速度行走smv/230sin4 此时，淋雨总量为 升24. 02/ )2/38 . 0(1095. 6 34 mC 这意味着你刚好跟着雨滴前进，前后都没淋雨。 第11页/共43页 当行走速度快于雨滴的水平运动速度，即 sinrv 你不断地追赶雨滴，雨水将淋湿你的前胸。被淋得雨量是 vrvpwD

9、h/ )sin( 淋雨总量为vrvhdrpwDC/)sin(cos / )sincos(rhvrdpwDrC 才可能小。尽可能大，当Cvrd, 0sincos 才可能小。尽可能小，当Cvrd, 0sincos ，而sinrv ，所以sinrv 才可能小。C 升。 时，取 77. 06/ )634 . 0(1095. 6 30,/6 34 mC smv 第12页/共43页 若雨是迎着你前进的方向向你若雨是迎着你前进的方向向你 落下落下,这时的策略很简单这时的策略很简单,应以应以 最大的速度向前跑；最大的速度向前跑； 若雨是从你的背后落下若雨是从你的背后落下,你应控你应控 制你在雨中的行走速度制你

10、在雨中的行走速度,让它刚让它刚 好等于落雨速度的水平分量。好等于落雨速度的水平分量。 第13页/共43页 例二例二 森林救火森林救火 森林失火后，要确定派出消防队员的数量森林失火后，要确定派出消防队员的数量. 队员多，森林损失小，救援费用大；队员多，森林损失小，救援费用大； 队员少，森林损失大，救援费用小队员少，森林损失大，救援费用小. 综合考虑损失费和救援费，确定队员数量综合考虑损失费和救援费，确定队员数量. 问问 题题 分分 析析 问问 题题 记队员人数记队员人数x, 失火时刻失火时刻t=0, 开始救火时刻开始救火时刻 t1, 灭火时刻灭火时刻t2, 时刻时刻t森林烧毁面积森林烧毁面积B(

11、t). 损失费损失费f1(x)是是x的减函数的减函数, 由烧毁面积由烧毁面积B(t2)决定决定. 救援费救援费f2(x)是是x的增函数的增函数, 由队员人数和救火时间决定由队员人数和救火时间决定. 存在恰当的存在恰当的x，使，使f1(x), f2(x)之和最之和最 小小. 第14页/共43页 关键是对关键是对B(t)作出合理的简化假设作出合理的简化假设. 问题问题 分析分析 失火时刻失火时刻t=0, 开始救火时刻开始救火时刻t1, 灭火时刻灭火时刻 t2, 画出时刻画出时刻t森林烧毁面积森林烧毁面积B(t)的大致图的大致图 形形. t1t20t B B(t2) 分析分析B(t)比较困难比较困难

12、, 转而讨论单位时间转而讨论单位时间 烧毁面积烧毁面积 dB/dt (森林烧毁的速度森林烧毁的速度). 第15页/共43页 模型假设模型假设 3）f1(x)与与B(t2)成正比，系数成正比，系数c1 (烧毁单位面积损失费烧毁单位面积损失费 ） 1）0 t t1, dB/dt 与与 t成正比，系数成正比，系数 (火势蔓延速度火势蔓延速度 ). 2）t1 t t2, 降为降为 - x ( 为队员的平均灭火为队员的平均灭火速度速度). 4）每个）每个队员的单位时间灭火费用队员的单位时间灭火费用c2, 一次性费用一次性费用c3 . 假设假设1）的解释）的解释 r B 火势以失火点为中心，均匀向四火势以

13、失火点为中心，均匀向四 周呈圆形蔓延，半径周呈圆形蔓延，半径 r与与 t 成正成正 比比. 面积面积 B与与 t2 成正比成正比 dB/dt与与 t 成正比成正比 第16页/共43页 x b tt 12 2 0 2) ( t dt dt dB tB 模型建立模型建立 dt dB b 0 t1t t2 x 假设假设1） , 1 tb xcttxcxftBcxf 31222211 )()(),()( 目标函数目标函数总费总费 用用 )()()( 21 xfxfxC 假设假设3）4 ） x t tt 1 12 假设假设2 ） )(222 2 1 22 12 x ttbt 第17页/共43页 0 dx

14、 dC xc x xtc x tctc xC 3 12 2 1 2 1 2 11 )(22 )( 模型建立模型建立目标函数目标函数总费总费 用用 模型求模型求 解解 求求 x使使 C(x)最小最小 2 3 12 2 11 2 2 c tctc x 结果解结果解 释释 / 是火势不继续蔓延的最少队员数是火势不继续蔓延的最少队员数 dt dB b 0 t1t2 t x 其中其中 c1,c2,c3, t1, , 为已知参数为已知参数 第18页/共43页 模模 型型 应应 用用 c1,c2,c3已知已知, t1可估计可估计, c2 x c1, t1, x c3 , x 结果结果 解释解释 2 3 12

15、 2 11 2 2 c tctc x c1烧毁单位面积损失费烧毁单位面积损失费, c2每个每个队员单位时间灭火队员单位时间灭火 费费, c3每个每个队员一次性费用队员一次性费用, t1开始救火时刻开始救火时刻, 火火势蔓延速度势蔓延速度, 每个每个队员平均灭火队员平均灭火速度速度. 为什么为什么? ? , 可可设置一系列数设置一系列数 值值 由模型决定队员数量由模型决定队员数量 x 第19页/共43页 50万元基金用于投资三种股票万元基金用于投资三种股票A、B、C： A每股年期望收益每股年期望收益5元元(标准差标准差2元元)，目前市价，目前市价20元；元； B每股年期望收益每股年期望收益8元元

16、(标准差标准差6元元)，目前市价，目前市价25元；元； C每股年期望收益每股年期望收益10元元(标准差标准差10元元)，目前市价，目前市价30元元 ； 股票股票A、B收益的相关系数为收益的相关系数为5/24; 股票股票A、C收益的相关系数为收益的相关系数为0.5; 股票股票B、C收益的相关系数为收益的相关系数为0.25。 假设：假设：1、基金不一定要用完（不用不计利息或贬值）、基金不一定要用完（不用不计利息或贬值） 2、风险通常用收益的方差或标准差衡量、风险通常用收益的方差或标准差衡量 第20页/共43页 A、B、C每手每手(百股百股)的收益分别记为的收益分别记为S1,S2和和S3(百元百元)

17、 ： ES1=5, ES2=8, ES3=10， DS1=4, DS2=36, DS3=100， r12=5/24, r13=-0.5，r23=-0.25 121212 131313 232323 cov( ,)25 cov( ,)10 cov(,)15 S SrDSDS S SrDSDS S SrDSDS 决策向量决策向量 x1 、x2和和 x3 分别表示投资分别表示投资A、B、C的数的数 量（国内股票通常以量（国内股票通常以“一手一手”（100股）为最小单位股）为最小单位 出售，这里以出售，这里以100股为单位，期望收益以百元为单位股为单位，期望收益以百元为单位 ） 总收益总收益 S=x1

18、S1+x2S2+x3S3 ：是一个随机变量：是一个随机变量 第21页/共43页 总期望收益为总期望收益为 Z1=ES= x1ES1+x2ES2+x3ES3=5x1+8x2+10 x3 投资风险（总收益的方差）为投资风险（总收益的方差）为 21 122331 12233 1 1221 1332233 222 1122331 212 1 3132 323 222 1231 21 32 3 ()()()() 2cov(,) 2cov(,) 2cov(,) 2cov( ,) 2cov( ,) 2cov(,) 43610052030 ZD xSx Sx SD xSD x SD x S xS x SxS

19、x Sx S x S x DSx DSx DSx xS S x xS Sx xS S xxxx xx xx x 第22页/共43页 222 21231 21 32 3 min43610052030Zxxxxxxxx x s.t. 5x1 +8x2+10 x3 1000 20 x1+25x2+30 x3 5000 x1，x2，x3 0 解得解得x = 1.0e+002 *（1.3111，0.1529，0.2221） 如果一定要整数解，可以四舍五入到（如果一定要整数解，可以四舍五入到（131，15，22 ） 如利用如利用LINGO软件软件,可得整数最优解可得整数最优解(132，15，22) 用去资

20、金为用去资金为132 20+15 25+22 30 = 3675（百元）（百元） 期望收益为期望收益为132 5+15 8+22 10 = 1000（百元）（百元） 风险风险(方差方差)为为 68116，标准差约为，标准差约为261（百元）（百元） 第23页/共43页 例四例四 男生追女生模型男生追女生模型 问问 题题 某男生某男生A对于某女生对于某女生B非常喜欢，但是刚开始的时候该非常喜欢，但是刚开始的时候该 女生对该男生并没有好感，该男生想采取一些行动来女生对该男生并没有好感，该男生想采取一些行动来 改变二者之间的关系，但是男女之间的过多接触势必改变二者之间的关系，但是男女之间的过多接触势

21、必 会对学习成绩造成影响，试问该男生能否在保持学习会对学习成绩造成影响，试问该男生能否在保持学习 成绩不下降的前提下追到该女生？成绩不下降的前提下追到该女生？ 要要 求求 建立适当的数学模型分析男生建立适当的数学模型分析男生A的学习成绩与女生的学习成绩与女生 B对该男生的好感之间的关系，并对模型作出解释。对该男生的好感之间的关系，并对模型作出解释。 第24页/共43页 模模 型型 假假 设设 1.A男生的学习成绩与男生的学习成绩与B女生对于女生对于A男生的疏远度均为男生的疏远度均为 时间时间 t的函数，分别设为的函数，分别设为Y(t)和和X(t)。 2. 初始时刻初始时刻X(t)是随着时间是随

22、着时间t 增长的（增长的（B女生发现了女生发现了A 男生的一些缺点），假设增长符合男生的一些缺点），假设增长符合Malthus模型，模型， 即：即：dX/dt=aX(t)dX/dt=aX(t) 其中其中a为增长率。为增长率。 3. 随着随着A A男生对男生对B B女生发动追求攻势后，女生发动追求攻势后，A A男生的学习男生的学习 成绩成绩Y(t)Y(t)呈现自然下降，假设也符合呈现自然下降，假设也符合Malthus模型，模型， 即：即：dY/dt=-eY(t) dY/dt=-eY(t) 其中其中e e为增长率。为增长率。 4. 当当Y(t)存在时，单位时间内存在时，单位时间内X(t)的减少值与

23、的减少值与X(t)成正成正 比，比例系数为常数比，比例系数为常数b。 5. 假定假定A男生对男生对B女生发动追求攻势后，立即转化成女生发动追求攻势后，立即转化成B女女 生对生对A男生的好感，对学习有帮助，设转化系数为男生的好感，对学习有帮助，设转化系数为 。 第25页/共43页 模模 型型 建建 立立 被食者被食者-食者食者 Volterra模型模型 dY(t) ( bX(t)e)Y(t) dt dX(t) (bY(t)X(t) dt a 这样就得到了一个在无外界干扰的条件下，学习成绩这样就得到了一个在无外界干扰的条件下，学习成绩 与疏远度相互作用的模型。这个模型在生物学中称为与疏远度相互作用

24、的模型。这个模型在生物学中称为 被食者和食者的被食者和食者的Volterra模型。模型。 初始条件：初始条件： 0 0 (0) (0) Xx Yy 按照前面的假设列出按照前面的假设列出Y(t)和和X(t)符合的关系式符合的关系式: 第26页/共43页 模型求解模型求解 这个方程组是一个非线性方程组，不易直接求解，将这个方程组是一个非线性方程组，不易直接求解，将 两个方程相除得微分方程两个方程相除得微分方程 ( )( ) ( ) ( )( )( ) dY tbX te Y t dX tabY tX t 分离变量积分后得到隐式解：分离变量积分后得到隐式解： lnlnlnbYbXaYeXCC为任意常

25、数为任意常数 以初始条件代入不难确定以初始条件代入不难确定C的值，从而可以得到一个特的值，从而可以得到一个特 解，它是解，它是X-Y平面上的一条闭曲线，只要初始值不为零平面上的一条闭曲线，只要初始值不为零, 这条闭曲线就永远不通过零点。这条闭曲线就永远不通过零点。 第27页/共43页 令：令：(, )lnlnF X YbYbXaYeX 模型分析模型分析 容易求出函数容易求出函数F有唯一的极小点有唯一的极小点(,) ea M cb 同时易见：当同时易见：当X （B女生对女生对A男生恨之入骨）或男生恨之入骨）或 Y （A男生是一块只会学习的男生是一块只会学习的“木头木头”）时均有）时均有 F ，而

26、：当，而：当0X 0Y （A男生不学无术）时男生不学无术）时 F （A男生属于天皇巨星，男生属于天皇巨星，B女生女生 对对A男生毫无防备）或男生毫无防备）或 也有也有，由此不难看出，由此不难看出F的图像是以的图像是以M为最小值为最小值 在第一卦限向上无限延伸的曲面，而在第一卦限向上无限延伸的曲面，而 (, ),0F X YZ Z是环绕点是环绕点M的闭曲线簇。的闭曲线簇。 第28页/共43页 模型应用模型应用 通过上面的分析可以知道通过上面的分析可以知道A男生的学习成绩与男生的学习成绩与B女生女生 对他的疏远度是呈周期性变化的，从生态意义上可以对他的疏远度是呈周期性变化的，从生态意义上可以 理解

27、为：当理解为：当A男生的学习成绩下降时，男生的学习成绩下降时，B女生会远离女生会远离A 男生，于是男生，于是A男生又开始发奋图强，学习成绩男生又开始发奋图强，学习成绩Y(t)又开又开 始上升，于是始上升，于是B女生又开始和女生又开始和A男生来往，疏远度降低；男生来往，疏远度降低； 交往多了，自然又分散了学习的时间，交往多了，自然又分散了学习的时间，A男生的学习成男生的学习成 绩绩Y(t)又开始下降。这样周而复始，形成了一个动态平又开始下降。这样周而复始，形成了一个动态平 衡。我们还可以证明，虽然对于不同的初始值可能出衡。我们还可以证明，虽然对于不同的初始值可能出 现不同的闭轨线，但在一个周期内

28、现不同的闭轨线，但在一个周期内X和和Y的平均数量都的平均数量都 分别是一个常数，而且恰为平衡点分别是一个常数，而且恰为平衡点M的两个坐标，这说的两个坐标，这说 明初始情况并不是决定明初始情况并不是决定A男生能否追到男生能否追到B女生的决定因女生的决定因 素。素。 第29页/共43页 模型的进一步讨论模型的进一步讨论 前面的结果都是在不考虑其他外界因素影响的前提下进前面的结果都是在不考虑其他外界因素影响的前提下进 行的，如果存在一些外界影响会对结果有些什么影响呢？行的，如果存在一些外界影响会对结果有些什么影响呢？ 考虑两种外界影响：考虑两种外界影响： 1.A男生的朋友对于男生的朋友对于A男生非常

29、支持，并且对于男生非常支持，并且对于A男生男生 追追B女生提供便利条件。女生提供便利条件。 2.出现一个出现一个C男生也在追男生也在追B女生，对于女生，对于A男生能否追上男生能否追上B 女生造成极大的威胁。女生造成极大的威胁。 根据根据Volterra原理，上面两种情况都会使得原理，上面两种情况都会使得A男生的学男生的学 习成绩习成绩Y(t)下降，同时下降，同时B女生对于女生对于A男生的疏远程度男生的疏远程度X(t) 增加。增加。 第30页/共43页 对于男生的一点儿忠告对于男生的一点儿忠告 通过上面的分析可以看出，初始情况对于结果的影响通过上面的分析可以看出，初始情况对于结果的影响 并不大，

30、一些成绩不好的同学也不要自卑，另外即使并不大，一些成绩不好的同学也不要自卑，另外即使 女同学对于你的某些缺点极为反感也不能决定最终的女同学对于你的某些缺点极为反感也不能决定最终的 结果，也许努力去追求就会得到接受。切忌强大的爱结果，也许努力去追求就会得到接受。切忌强大的爱 情攻势是不一定能达到满意的效果的，反而不利于学情攻势是不一定能达到满意的效果的，反而不利于学 业。有时通过慢慢的接触，慢慢的了解，再加上适当业。有时通过慢慢的接触，慢慢的了解，再加上适当 的追求行动，女生的疏远程度会慢慢降低，你的学习的追求行动，女生的疏远程度会慢慢降低，你的学习 成绩还不会下降！成绩还不会下降！ 注：以上观

31、点均属于个人看法，不具有指导意义！注：以上观点均属于个人看法，不具有指导意义！ 第31页/共43页 v1能源利用量能源利用量, v2能源价格能源价格, v3能源生产率能源生产率, v4环境质量环境质量, v5工业产值工业产值, v6就业机会就业机会, v7人口总数人口总数. 例五例五 社会经济系统的冲量过程社会经济系统的冲量过程 系统的元素系统的元素图的顶点图的顶点 元素间的直接影响元素间的直接影响有方向的弧有方向的弧 正面影响正面影响弧旁的弧旁的+号；负面影响号；负面影响弧旁的弧旁的号号 带符号的有向图带符号的有向图 符号、符号、 客观规律；客观规律；方针政策方针政策 例例 能源利用系统的预

32、测能源利用系统的预测 + - + - + + + + - - + v2 v1 v3 v4 v6 v7 v5 第32页/共43页 Evv vv vv a ji ji ji ij 若， 为若 为若， 0 , 1 1 0000001 1000000 0100001 1000000 0010010 0000001 0001110 A 带符号有向图带符号有向图G1=(V,E)的邻接矩阵的邻接矩阵A V顶点集顶点集 , E弧弧 集集 定性模型定性模型 - vivj + 某时段某时段vi 增加导致增加导致 下时段下时段vj 增加增加(减少减少 ) 带符号的有向图带符号的有向图G1 + - + - + + +

33、 + - - + v2 v1 v3 v4 v6 v7 v5 第33页/共43页 0000005 . 1 1000000 05 . 100002 . 1 3 . 0000000 0010020 0000007 . 0 0002 . 18 . 05 . 00 W 加权有向图加权有向图G2及其邻接矩阵及其邻接矩阵W 定量模型定量模型 某时段某时段vi 增加增加1单位导单位导 致下时段致下时段vj 增加增加wij单位单位 j w i vv ij 的特例视为 WA v7 0.3 1 1.5 1 1.5 1.2 0.8 -2 -2 -0.7 -0.5 v1 v2 v3 v4 v5v6 加权有向图加权有向图

34、G2 第34页/共43页 , 2 , 1 , 0, 2 , 1),1()() 1(tnitptvtv iii n i n i iijjiijj tpatptpwtp 11 )()1(),()1(或 ) 1()() 1(tptvtv 冲量过程冲量过程（Pulse Process) 研究由某元素研究由某元素vi变化引起的系统的演变过程变化引起的系统的演变过程 vi(t) vi在时段在时段t 的的值值； pi(t) vi在时段在时段t 的的改变量改变量(冲量冲量 ) )(),(),()(),(,),(),()( 2121 tptptptptvtvtvtv nn j w i vv ij 冲量过程模型冲

35、量过程模型 Wtptp)()1(Atptp)()1(或或 第35页/共43页 能源利用系统的预测能源利用系统的预测 简单冲量过程简单冲量过程初始冲量初始冲量p(0)中中 某个分量为某个分量为1，其余为，其余为0的冲量过的冲量过 程程. 若开始时能源利用量有突然增加，预测系统的演变若开始时能源利用量有突然增加，预测系统的演变. )0()0(pv ) 1()() 1(tptvtv Atptp)() 1( 设设 能源利用系统的能源利用系统的 p(t)和和v(t) -1 1 0-11-100011-10000 t 4 p 3 p 5 p 6 p 7 p 2 p 4 v 3 v 2 v 1 v 5 v

36、6 v 7 v 0100000 0 100000 0 0 1 p 2 3 1 -10010-12-21-110-1 1-11-10103-32-21 1 -1 第36页/共43页 简单冲量过程简单冲量过程S的稳定性的稳定性 任意时段任意时段S的各元素的值和冲量是否为有限的各元素的值和冲量是否为有限(稳定稳定 )? S不稳定时如何改变可以控制的关系使之变为稳定不稳定时如何改变可以控制的关系使之变为稳定? S冲量稳定冲量稳定对任意对任意 i,t, | pi(t) |有有 界界 S值稳定值稳定对任意对任意 i,t, | vi(t) |有界有界 值稳定值稳定 冲量稳定冲量稳定 ) 1()() 1(tp

37、tvtvWtptp)()1( t Wptp)0()(S的稳定性取决于的稳定性取决于W的特征根的特征根 记记W的的非零非零特征根为特征根为 第37页/共43页 S冲量稳定冲量稳定 | | 1 S冲量稳定冲量稳定 | | 1且均为单根且均为单根 S值稳定值稳定 S冲量稳定冲量稳定且且 不等于不等于1 0000001 1000000 0100001 1000000 0010010 0000001 0001110 A 对于能源利用系统的邻接矩阵对于能源利用系统的邻接矩阵A ) 1()( 2352 f 特征多项式特征多项式 76)2(, 2) 1 (ff) 2 , 1 ( 能源利用系统存在能源利用系统存在冲量冲量 不稳定不稳定的简单冲量过程的简单冲量过程 简单冲量过程简单冲量过程S的稳定性的稳定性 第38页/共43页 简单冲量过程的稳定性简单冲量过程的稳定性 改进的玫瑰形图改进的玫瑰形图S* 带符号的带符号的 有向图双向连通，且存在一有向图双向连通，且存在一 个位于所有回路上的中心顶个位于所有回路上的中心顶 点点.