第 1 页 实验八： 差分方程模型
实验目的：
1. 掌握一阶和二阶常系数线性差分方程的求解方法；
2. 利用差分方程建模。
实验练习：
1. 求下列一阶常系数线性差分方程的通解：
(1) 134nnyy++=；
(2) 12cos()nnnyynp+-=.
注：若12()cossinfnbnbnww=+，则特解可设为
*cossinnynnawbw=+其中,ab为待定系数，求解关于,ab的方程组，若系数矩阵为0，则特解可设为
*(cossin)nynnnawbw=+
2. 教材P173，习题4.
3. 求下列二阶常系数线性差分方程的通解：
()()2222cossin.nnnnyypp++=-安徽师范大学
数计 学院实验报告
第 2 页
专业名称数学与应用数学
实 验 室2号实验楼#201
实验课程Matlab
实验名称差分方程模型
姓 名 张顺强
学 号 100701185
同组人员 无
实验日期 2013.
第 3 页 注：实验报告应包含（实验目的，实验原理，主要仪器设备和材料，实验过程和步骤，实验原始数据记录和处理，实验结果和分析，成绩评定）等七项内容。具体内容可根据专业特点和实验性质略作调整，页面不够可附页。
一.实验目的
1. 掌握一阶和二阶常系数线性差分方程的求解方法；
2. 利用差分方程建模。
二.实验原理
数学建模相关理论
三主要仪器设备和材料
计算机，数学建模相关书籍
四.实验过程和步骤
1. 求下列一阶常系数线性差分方程的通解：
(1) 134nnyy++=；
(2) 12cos()nnnyynp+-=.
注：若12()cossinfnbnbnww=+，则特解可设为
*cossinnynnawbw=+其中,ab为待定系数，求解关于,ab的方程组，若系数矩阵为0，则特解可设为
*(cossin)nynnnawbw=+
数学建模中的差分方程模型
数学建模是一种将实际问题转化为数学模型并寻求与之相连的数学方法的学科，不仅仅在理论研究上有很大的应用，也在实际生活中有着广泛的应用。在各种数学模型中，差分方程模型也是一种很重要的模型。本文将结合实例，介绍差分方程模型的定义、建立、求解以及应用。
差分方程模型定义
差分方程模型是一种通过离散化的方法，将连续时间问题转化为离散时间问题，来描述变量随时间的变化规律的数学模型。这种数学模型以时间为自变量，以某个状态量为因变量，由一定的关系式组成。例如：y(n+1)=ay(n)+b，式子中y(n)代表第n时刻系统状态，y(n+1)代表第n+1时刻系统状态，a和b为常数。
差分方程模型建立
建立差分方程模型的关键是将实际问题中的连续变化离散化。一般情况下，对于所建立的模型，首先要确定它的思路和范围，然后根据实际情况，确定差分方程的形式。此外，还需要进行参数的估计和参数变化的分析，以及对模型精确性的验证。
以物理学中的简谐振动为例，建立一个差分方程模型描述其运动，即一个质点在回复力作用下以简谐运动形式振动。设t为时间，y为质点的位移，v为质点的速度，a为质点的加速度，则有：
$$
y=n\Delta y \\
v=\dfrac{y(n+1)-y(n-1)}{2\Delta t} \\
a=\dfrac{y(n+1)-2y(n)+y(n-1)}{(\Delta t)^2}
$$
其中n为时间步长，$\Delta t$为时间间隔。我们利用受力平衡的原理，即简谐振动中的$F=-ky$得到：
$$
\dfrac{y(n+1)-2y(n)+y(n-1)}{(\Delta t)^2} = -\dfrac{k}{m}y(n)
$$
将$\alpha=\dfrac{k}{m}$带入上式得到：
$$
y(n+1)-2(1+\alpha)y(n)+y(n-1) = 0
$$
此时，我们便成功地建立了描述简谐振动的差分方程模型。
差分方法建模
经过一段时间的学习，我对《数学建模》从最初的认知，到最后系统性的了解，并且在老师的帮助下，我渐渐地喜欢上这么课程，由此把自己的从这门课程中学到的东西和认识给写下来，希望老师给予知道和讲评，我也会在以后的日子里加以改进。
首先说说自己对《数学建模》的认识。最初印象是从图书馆里开始认识的，那个时候自己还是大一刚来的小青年，对一切都充满着好奇和新鲜，看到《数学建模》这本书的时候，我随手拿过来看，发现里面对事物的分析方法和模型的建立都很详细，我就试着往下看，感觉这本书非常好，幸好今年有这门课程的建立，更加是我产生了浓厚的兴趣，在以后的日子里，每节课堂我都认真的听老师讲课，经过这半年的学习，自己多少对《数学建模》有了一些认识，利用课间时间我查阅了一些资料，是自己更加充分的了解它，下面是我对其的认识，数学建模是在20世纪60和70年代进入一些西方国家大学的，我国的几所大学也在80年代初将数学建模引入课堂。经过20多年的发展现在绝大多数本科院校和许多专科学校都开设了各种形式的数学建模课程和讲座，为培养学生利用数学方法分析、解决实际问题的能力开辟了一条有效的途径。
对于现在的大学上来说，大学生数学建模竞赛最早是1985年在美国出现的，1989年在几位从事数学建模教育的教师的组织和推动下，我国几所大学的学生开始参加美国的竞赛，而且积极性越来越高，近几年参赛校数、队数占到相当大的比例。可以说，数学建模竞赛是在美国诞生、在中国开花、结果的。 1992年由中国工业与应用数学学会组织举办了我国10城市的大学生数学模型联赛，74所院校的314队参加。教育部领导及时发现、并扶植、培育了这一新生事物，决定从1994年起由教育部高教司和中国工业与应用数学学会共同主办全国大学生数学建模竞赛，每年一届。十几年来这项竞赛的规模以平均年增长25%以上的速度发展。
全国大学生数学建模竞赛是全国高校规模最大的课外科技活动之一。本竞赛每年9月(一般在中旬某个周末的星期五至下周星期一共3天，72小时）举行，竞赛面向全国大专院校的学生，不分专业（但竞赛分本科、专科两组，本科组竞赛所有大学生均可参加，专科组竞赛只有专科生（包括高职、高专生）可以参加）。 2008 年全国有31个省/市/自治区(包括香港)1023所院校、12846个队（其中甲组10384队、乙组2462队）、3万8千多名来自各个专业的大学生参加竞赛，是历年来参赛人数最多 二、数学建模的定义 简单地说：数学模型就是对实际问题的一种数学表述。 具体一点说：数学模型是关于部分现实世界为某种目的的一个抽象的简化的数学结构。 更确切地说：数学模型就是对于一个特定的对象为了一个特定目标，根据特有的内在规律，做出一些必要的简化假设，运用适当的数学工具，得到的一个数学结构。数学结构可以是数学公式，算法、表格、图示等。 数学建模就是建立数学模型，建立数学模型的过程就是数学建模的过程（见数学建模过程流程图）。 数学建模是一种数学的思考方法，是运用数学的语言和方法，通过抽象、简化建立能近似刻画并"解决"实际问题的一种强有力的数学手段。
第八章 差分方程模型
差分方程是解决离散时间问题的常用的数学方法，本章介绍几个用差分方程建立的实际问题的数学模型。
8.1个人住房抵押贷款
随着经济的发展,金融问题正越来越多地进入普通市民的生活，贷款、保险、养老金和信用卡等都涉及金融问题，个人住房抵押贷款是其中最重要的一项。1998年12月,中国人民银行公布了新的存、贷款利率水平,其中贷款利率如下表所列:
表8.1 中国人民银行贷款利率表
贷款期限 半年 一年 三年 五年 五年以上
利率﹪ 6.12 6.39 6.66 7.20 7.56
当贷款期处于表中所列相邻年限之间时利率为对应相邻两数中较大者。其后,上海商业银行对个人住房商业性贷款利率作出相应调整。表8.2和表8.3分别列出了上海市个人住房商业抵押贷款年利率和商业抵押贷款(万元)还款额的部分数据(仅列出了五年)。
表8.2 上海市商业银行住房抵押贷款利率表
贷款期限 一年 二年 三年 四年 五年
利率﹪ 6.12 6.225 6.390 6.525 6.660
表8.3 上海市商业银行住房抵押贷款分期付款表
贷款期限 一年 二年 三年 四年 五年
月还款（元） 到期一次还清
本息总和（元） 10612.00 444.36
10664.54 305.99
11015.63 237.26
11388.71 196.41
11784.71
一个自然的问题是,表8.2和表8.3是如何依据中央人民银行公布的存、贷款利率水平制定的？
我们以商业贷款10000元为例,一年期贷款的年利率为6.12﹪,到期一次还本付息总计10612.00元,这很容易理解。然而二年期贷款的年利率为6.225﹪,月还款数444.36元为本息和的二十四分之一,这后两个数字究竟是怎样产生的？是根据本息总额算出月还款额，还是恰好相反？让我们稍微仔细一些来进行分析。由于贷款是逐月归还的,就有必要考察每个月欠款余额的情况。
