差分方程模型-

1、差分方程模型数学建模讲座一、关于差分方程模型简单的例子1. 血流中地高辛的衰减地高辛用于心脏病。考虑地高辛在血流中的衰减问题以开出能使地高辛保持在可接受(安全而有效的水平上的剂量处方。假定开了每日0.1毫克的剂量处方,且知道在每个剂量周期(每日末还剩留一半地高辛,则可建立模型如下:设某病人第n 天后血流中地高辛剩余量为n a , 则2. 养老金问题对现有存款付给利息且允许每月有固定数额的提款, 直到提尽为止。月利息为1,月提款额为1000元,则可建模型如下:设第n 月的存款额为n a ,则100001.11=+n n a a (一阶非齐次线性差分方程3. 兔子问题(Fibonacci 数设第一

2、月初有雌雄各一的一对小兔,假定两月后长成成兔,同时(即第三个月开始,每月初产雌雄各一的一对小兔, 新增小兔也按此规律繁殖,设第n 月末共有n F 对兔子,则建模如下:=+=12121F F F F F n n n (二阶线性差分方程初值问题 34232143212211F F F F F F F F F F +=+ 注意上月新生的小兔不产兔(因第n 月末的兔子包括两部分, 一部分上月留下的为1n F , 另一部分为当月新生的,而新生的小兔数=前月末的兔数4.车出租问题 A ,B 两地均为旅游城市,游客可在一个城市租车而在另一个城市还车。 A , B 两汽车公司需考虑置放足够的车辆满足用车需要,

3、以便估算成本。分析历史记录数据得出:n x : 第n 天营业结束时A 公司的车辆数 n y :第n 天营业结束时B 公司的车辆数则 +=+=+n n n n1 (一阶线性差分方程组(问题模型可进一步推广二 差分方程的解法类比: 差分方程是数列间关系; 微分方程是函数间关系 定义1. 形如02211=+k n k n n n a b a b a b a L 的差分方程,称为n a 的k 阶常系数线性齐次差分方程,其中i b 为常数,0k b 且k n . 02211=+k k k k b x b x b x L 称为差分方程的特征方程,其根称为特征根。定理1(单根情形若特征方程恰有k 个相异的特

4、征根k x x x ,21L ,则差分方程的通解为 nk k n n n x c x c x c a +=L 2211.例1 求解兔子问题 =+=12121F F F F F n n n解:差分方程的特征方程为 012=x x特征根 251,25121=+=x x 通解为 nnn c c F+ +=25125121, 由初始条件 121=F F 得: = +=+2(12512511(1251251222121c c c c 解得 =515121c c 故 +=nn n F 25121551定理 2 (重根情形 若特征方程的相异特征根为12,t x x x L , 重数依次为t m m m ,2

5、1L , 其中k m m m t =+L 21,则差分方程的通解为:112(t t m n t t tm t c c n c n x +L L(定理1包含在定理2之中定理 3 若差分方程的特征方程的特征根出现一对共轭虚根,=+=ivu x ivu x 21 和相异的2k 个根k x x ,3L , 则差分方程的通解为:nk k nnnn x c x c n c n c a +=L 3321sin cos .定义2 形如(2211n f a b a b a b a k n k n n n =+L(k b b b ,21L 为常数,0k b ,0(n f ,k n 的差分方程为k 阶常系数线性非齐

6、次差分方程。称011=+k n k n n a b a b a L为其对应的齐次方程。定理4 非齐次差分方程的通解等于对应齐次差分方程的通解加上非齐次方程的特解,即n nn a a a +=*,其中*n a 为通解, n a 为特解.齐次特征方程 05.0=, 齐次方程通解 n nnc a 例3 (养老金问题解法1 100001.11=+n n a a齐次特征方程001.1=, 齐次方程通解nnc a 01.1(*=. 设特解为 D a n =, 代入代入原方程得 100000=D . 通解为 10000001.1(+=n n c a .解法2(化齐 100001.11=+n n a a ,1

7、00001.11=n n a a .=+.01(01.1(=,通解为1(1.01n n a c c =+,代入原方程得故 10000001.1(+=nnc a . 例4 求非齐次差分方程nn n n a a a 24421=+的通解. 解法1 齐次特征方程0442=+, 二重根2=,对应齐次方程的通解为n n n n c c a 2221*+=. 因nn f 2(=中, 2 是2 重根,故设特解为nn n A a 22=代入得12A =, 故通解为1221222+=n n n n n n c c a 方法2(化齐 nn n n a a a 24421=+,13212244(2=+n n n n a a a相减得12361280n n n