第二讲简介_差分方程模型与案例.docVIP

第二讲 简介_差分方程模型与案例 1.1 差分方程的基本定义 差分方程是在离散时段上描述现实世界中变化过程的数学模型。现实中的问题通常是连续变化的，但我们常常只能在离散的时间点上对其进行观测和描述。为了表述这一类的数学模型，我们引入了差分方程的方法。 1.2 一阶线性常系数差分方程 一阶线性常系数差分方程的一般形式 差分方程的平衡点 差分方程的解 平衡点稳定的条件 1.3高阶线性常系数差分方程 n阶线性常系数差分方程的一般形式 称方程 (*) 为对应的齐次方程 设方程(*)有形如的解,则有 特征方程 特征根(假定有k个不同的实根,其它情形参见资料《差分方法建模理论与案例》) 平衡点 差分方程的解 ???? 平衡点稳定的条件 所有特征值的模均小于1 1.4非线性差分方程 建模案例： 题目1? 濒危物种的自然演变和人工孵化 问题：Florida 沙丘鹤属于濒危物种，生态学家估计它在自然环境下，年平均增长率为 -3.24%。如果在某自然保护区内开始有100只鹤，每年人工孵化5只鹤放入该保护区，建立描述其数量变化规律的模型，并作数值计算。 模型及其求解 记第k年沙丘鹤的数量为xk，自然环境下的年平均增长率为r，记a=1+r，每年孵化的数量为b，则第k+1年鹤的数量为 模型分析 讨论时间充分长以后沙丘鹤数量的变化趋势, 即k→∞时xk的极限状态。 当a1即r0时xk→x=b/(1-a)。 　 MATLAB演示计算 ? 计算并作图，程序如下： 子程序： function x=exf11(x0,n,r,b)?????? % 建立名为exf11的函数M文件， x0,n,r,b可调节 a=1+r; x=x0;?????????????????????? % 赋初值 for k=1:n ??? x(k+1)=a*x(k)+b;?????? ????% 按照（3）迭代计算 end ? 主程序: clc; clf; clear all % x0:初始值；r=-0.0324;? b:人工孵化数 x0=100;n=20;b=5; r=-0.0324; % 给定x0,n,r，b,调用exf11计算 k=(0:n); y=exf0201(x0,n,r,b); plot(k,y,r*-),?? title(Florida 沙丘鹤数量变化趋势);% 在图上做标记 ??? （运行结果显示） 　 ?题目2? 汽车租赁公司的运营 问题:一家汽车租赁公司在3个相邻的城市运营，为方便顾客起见公司承诺，在一个城市租赁的汽车可以在任意一个城市归还。根据经验估计和市场调查，一个租赁期内在A市租赁的汽车在A, B, C市归还的比例分别为0.6, 0.3, 0.1；在B市租赁的汽车在A, B, C市归还的比例分别为0.2, 0.7, 0.1；在C市租赁的汽车在A, B, C市归还的比例分别为0.1, 0.3, 0.6。若公司开业时将N辆汽车按一定方式分配到3个城市，建立运营过程中汽车数量在3个城市间转移的模型，并讨论时间充分长以后的变化趋势。 模型及其求解 记第k个租赁期末公司在A, B, C市的汽车数量分别为x1(k), x2(k), x3(k) ,容易写出第k+1（k=0,1,2, …）个租赁期末公司在A, B, C市的汽车数量为? ????????? ????????????? 记向量x(k)=[x1(k), x2(k), x3(k)]T，矩阵 ???????????? ???????????????????????????????????? 则????????? ?????????????????????????? 给定初始值x(0)，可以计算各个租赁期3个城市汽车数量的变化。 模型分析 猜想：时间充分长以后3个城市的汽车数量趋向稳定，并且稳定值与汽车的初始分配无关。 为了证实该猜想，记稳定值为x，由x应满足?表明矩阵A的一个特征根=1，且x是对应的特征向量。 MATLAB演示计算 初始分配：城市1：200，城市2：200 ，城市3：200 ? 计算并作图，程序如下： ? A=[0.6,0.2,0.1;0.3,0.7,0.3;0.1,0.1,0.6]; x(:,1)=[200,200,200];???????????????? % 赋初值 n=10; for k=1:n ??? x(:,k+1)=A*x(:,k);???????????????? % 迭代计算 end round(x), k=0:10; plot(k,x),grid, ? （运行结果显示） 表:计算结果 k 0 1 2 3 4 5 6 7 8 9 10 x1(k) 200 180 176 176 178 179 179 180 180 180 180 x2(k) 200 260 284