微分方程模型及差分模型.ppt

经济不稳定时政府的干预办法 1. 使 ? 尽量小，如 ?=0 以行政手段控制价格不变 2. 使 ? 尽量小，如 ? =0 靠经济实力控制数量不变 x y 0 y0 g f x y 0 x0 g f 结果解释 需求曲线变为水平 供应曲线变为竖直 2012-3-17 * Anna 模型的推广 生产者根据当前时段和前一时段的价格决定下一时段的产量。 生产者管理水平提高 设供应函数为 需求函数不变 二阶线性常系数差分方程 x0为平衡点 研究平衡点稳定，即k??, xk?x0的条件 2012-3-17 * Anna 方程通解 (c1, c2由初始条件确定) ?1, 2~特征根，即方程 的根 平衡点稳定，即k??, xk?x0的条件: 平衡点稳定条件 比原来的条件 放宽了 模型的推广 2012-3-17 * Anna 1.5 按年龄分组的种群增长 不同年龄组的繁殖率和死亡率不同 建立差分方程模型，讨论稳定状况下种群的增长规律 假设与建模 种群按年龄大小等分为n个年龄组，记i=1,2,… , n 时间离散为时段，长度与年龄组区间相等，记k=1,2,… 以雌性个体数量为对象 第i 年龄组1雌性个体在1时段内的繁殖率为bi 第i 年龄组在1时段内的死亡率为di, 存活率为si=1- di 2012-3-17 * Anna 假设 与 建模 xi(k)~时段k第i 年龄组的种群数量 ~按年龄组的分布向量 预测任意时段种群按年龄组的分布 ~Leslie矩阵(L矩阵) (设至少1个bi0) 2012-3-17 * Anna 稳定状态分析的数学知识 L矩阵存在正单特征根?1， 若L矩阵存在bi, bi+10, 则 P的第1列是x* 特征向量 , c是由bi, si, x(0)决定的常数 且 解释 L对角化 2012-3-17 * Anna 稳态分析——k充分大种群按年龄组的分布 ~ 种群按年龄组的分布趋向稳定，x*称稳定分布, 与初始分布无关。 ~ 各年龄组种群数量按同一倍数增减， ?称固有增长率 与基本模型 比较 3）?=1时 ~ 各年龄组种群数量不变 2012-3-17 * Anna ~ 1个个体在整个存活期内的繁殖数量为1 稳态分析 ~存活率 si是同一时段的 xi+1与 xi之比 （与si 的定义 比较） 3）?=1时 2012-3-17 * Anna 相关Matlab知识 命令：常微分方程的符号解 函数：dsolve 格式：r=dsolve(‘eq1,eq2,…’,’cond1,cond2,…’,’v’) 说明：对给定的常微分方程（组）eq1,eq2,…中指定的符号自变量v，与给定的边界条件和初始条件cond1,cond2,…，求符号解（即解析解）r；若没有指定变量v，则默认变量为t；在微分方程（组）的表达式eq中，大写字母D表示对自变量（设为x）的微分算子：D=s/sx,D2=d2/x2,….微分算子D后面的字母则表示为因变量，即待求解的未知函数. 2012-3-17 * Anna 相关Matlab知识 初始和边界条件由字符串表示：y(a)=b,Dy(c)=d,D2y(e)=f等，分别表示 若边界条件少于方程（组）的阶数，则返回的结果r中会出现任意常数C1,C2,…； dsolve命令最多可以接受12个输入变量（包括方程组与定解条件个数，当然可以做到输入的方程个数多于12个，只要将多个方程置于一个字符串内即可）. 若没有给定输出参量，则在命令窗口显示解列表. 2012-3-17 * Anna 相关Matlab知识 若该命令找不到解析解，则返回一警告信息，同时返回一空的sym对象. 这时，用户可以用命令ode23或ode45求解方程组的数值解. 例1：解常微分方程： （注：【程序】『输出结果』） （1）求 的通解； 【 s=dsolve(Dy=a*y+b)】 『s = -b/a+exp(a*t)*C1』 2012-3-17 * Anna 相关Matlab知识 （2）初值问题 【dsolve(Dy=y-2*t/y,y