差分模型说课.ppt

一、 差分的概念与性质 在社会经济活动与自然科学研究中，我们经常遇到与时 间t 有关的变量，而人们往往又只能观察或记录到这些变量 在离散的t 时的值。对于这类变量，如何去研究它们的相互 关系，就离不开差分与差分方程的工具。 一般地，在连续变化的时间范围内，变量y关于时间t的 变化率是用dy/dx来刻画的；对离散型的变量y，我们常取在 规定的时间区间上的差商 来刻画变量的变化率。 如果选 择 ，则 可以近似表示变量y的变化率。 由此我们给出差分的定义： 以t 表示时间，规 定t只取非负整数。t=0 表示第一周期初，t=1表示第二周期初。 记 为变量y在时刻t 时的取值，则称 为 的一阶差分，称 为的二阶差分。 一般地，函数 的n-1阶差分的差分称为n阶差分， 记为 ，即  差分的性质： （1） （2） （3） （4） 例2求 的差分 二、差分方程的概念 含有未知函数 的差分的方程为差分方程。 差分方程的一般形式: 其中含的最高阶差分的阶数称为该差分方程 的阶。例如，二阶差分方程 也可改写成 例3 试改变差分方程 的形式. 例4试确定下列差分方程的阶. 例5指出下列等式哪一个是差分方程, 若是, 进一步指出是否为线性方程. rsolve(方程，解函数,选项) ，rsolve({方程组，初始条件}，{解函数}，选项) 选项为‘genfunc’(x)解以x为自变量；选项为‘makeproc’解为过程函数。 rsolve({F(n) = F(n-1) + F(n-2), F(1..2)=1}, F, genfunc(x)); rsolve({s(n) = 2*s(n-1), s(0)=1}, s,makeproc) Matlab求解 濒危物种的自然演变和人工孵化 问题 Florida沙丘鹤属于濒危物种，它在较好自然环境下，年均增长率仅为1.94%，而在中等和较差环境下年均增长率分别为 -3.24% 和 -3.82%，如果在某自然保护区内开始有100只鹤，建立描述其数量变化规律的模型，并作 数值计算。 模型建立 xk+1=(1+r)xk k=0,1,2······ 记第k年沙丘鹤的数量为xk,年均增长率为r，则第k+1年鹤的数量为 已知x0=100, 在较好，中等和较差的自然环境下 r=0.0194, -0.0324,和-0.0382 我们利用Matlab编程，递推20年后观察沙丘鹤的数量变化情况 Matlab实现 首先建立一个关于变量n ，r的m函数 function x=sqh(n,r) a=1+r; x=100; for k=1:n x(k+1)=a*x(k); end 利用plot 绘图观察数量变化趋势 可以用不同线型和颜色绘图 r g b c m y k w 分别表示 红绿兰兰绿洋红黄黑白色 ： + o * . X s d 表示不同的线型 plot(k,y1,k,y2,k,y3) 在同一坐标系下画图 筹措教育经费模型 某家庭从现在着手, 从每月工资中拿出一部分资 金存入银行, 用于投资子女的教育, 并计算20年后 开始从投资账户中每月支取1000元, 直到10年后 子女大学毕业并用完全部资金. 要实现这个投 资目标, 20年内要总共筹措多少资金? 每月要在 银行存入多少钱? 假设投资的月利率为0.5%, 为 此, 设第t个月, 投资账户资金为每月存资金为b元, 于是20年后, 关于的差分方程模型为 作业 参 考 文 献 1.数学模型 姜启源（第二版） P271-303 高等教育出版社 2.数学建模方法 刘承平 P95-105 高等教育出版社 3.数学模型引论 唐焕文 、贺明峰（第二版） P239-247 高等教育出版社 1，1，2，3，5，8，13，…… 裴波那契数列递推公式：F(n+2) = F(n+1) + F(n) ，F(1)=F(2)=1。 它的通项求解如下： F(n+2) = F(n+1) + F(n) = F(n+2) - F(n+1) - F(n) = 0 令 F(n+2) - aF(n+1) = b(F(n+1) - aF(n)) 展开 F(n+2) - (a+b)F(n+1) + abF