差分方程模型.ppt

简单例子 例1: 汉诺塔问题：n个大小不同的圆盘依其半径大小依次套在桩A上，大的在下，小的在上，现要将n个盘移到空桩B或C上，但要求一次只能移动一个盘且移动过程中，始终保持大盘在下，小盘在上，移动过程中桩A也可利用，设移动n个盘的次数为 ，试建立关于的差分方程。 最新研究进展 （数学的实践与认识，2010,40（17）：40-46.） 利用 也就是 等与时间k无关， 化简得 下面给出MATLAB计算程序：%最优捕鱼策略ch431%文件名：ch431.mx=sym(x); b3=0.42*x;d=1.22*10^11;r=0.8;q=d*exp(-(3*r+2/3*b3))*(32529.55*exp(r)+ 65059.1*exp(-2/3*x)/(1-exp(-(r+2/3*x))));N10=d*q/(d+q);N40=d*q/(d+q)*exp(-(3*r+2/3*b3))/(1-exp(-(r+2/3*x)));a3=b3/(r+b3)*(1-exp(-2/3*(r+b3)))*d*q/(d+q)*exp(-2*r);a4=x/(r+x)*(1-exp(-2/3*(r+x)))*d*q/(d+q)*exp(-(3*r+2/3*b3))/ (1-exp(-(r+2/3*x))); M3=17.86;M4=22.99;M=M3*a3+M4*a4;M1=-M;M10=char(M);M11=char(M1);fplot(M10,[0,100])b4=fmin(M11,0,100);b3=0.42*b4;d=1.22*10^11;r=0.8;q=d*exp(-(3*r+2/3*b3))*(32529.55*exp(r)+ 65059.1*exp(-2/3*b4)/(1-exp(-(r+2/3*b4))));N10=d*q/(d+q);N40=d*q/(d+q)*exp(-(3*r+2/3*b3))/(1-exp(-(r+2/3*b4))); a3=b3/(r+b3)*(1-exp(-2/3*(r+b3)))*d*q/(d+q)*exp(-2*r); a4=b4/(r+b4)*(1-exp(-2/3*(r+b4)))*d*q/(d+q)*exp(-(3*r+2/3*b3))/ (1-exp(-(r+2/3*b4)));M3=17.86;M4=22.99;b3b4M=M3*a3+M4*a4;Max=MN1=N10 %各年龄组鱼的数量N2=N1*exp(-r)N3=N2*exp(-r)N4=N40 执行后输出 b3 = 7.3359b4 = 17.4664 Max =3.8886e+011各年龄组数为N1 = 1.1958e+011N2 = 5.3730e+010N3 = 2.4142e+010N4 = 8.1544e+007 2 市场经济中的蛛网模型 问 题 供大于求 现 象 商品数量与价格的振荡在什么条件下趋向稳定 当不稳定时政府能采取什么干预手段使之稳定 价格下降 减少产量 增加产量 价格上涨 供不应求 描述商品数量与价格的变化规律 数量与价格在振荡 蛛 网 模 型 g x0 y0 P0 f x y 0 xk~第k时段商品数量；yk~第k时段商品价格 消费者的需求关系 生产者的供应关系 减函数 增函数 供应函数 需求函数 f与g的交点P0(x0,y0) ~ 平衡点 一旦xk=x0，则yk=y0, xk+1,xk+2,…=x0, yk+1,yk+2, …=y0 x y 0 f g y0 x0 P0 设x1偏离x0 x1 x2 P2 y1 P1 y2 P3 P4 x3 y3 P0是稳定平衡点 P1 P2 P3 P4 P0是不稳定平衡点 x y 0 y0 x0 P0 f g ? ? 曲线斜率 蛛 网 模 型 ? 在P0点附近用直线近似曲线 P0稳定 P0不稳定 方 程 模 型 方程模型与蛛网模型的一致 ? ~ 商品数量减少1单位, 价格上涨幅度 ? ~ 价格上涨1单位, (下时段)供应的增量 考察? , ? 的含义 ? ~ 消费者对需求的敏感程度 ? ~ 生产者对价格的敏感程度 ?小, 有利于经济稳定 ? 小, 有利于经济稳定 结果解释 xk~第k时段商品数量；yk~第k时段商品价格 经济稳定 结果解释 经济不稳定时政府的干预办法 1. 使 ? 尽量小，如 ?=0 以行政手段控制价格不变 2. 使 ? 尽量小，如 ? =0 靠经济实力控制数量不变 x y 0 y0 g f x y 0 x0 g f 结果解释 需求曲线变为水平 供应曲线变为竖直 模型的推广