差分方程模型的稳定性分析及其应用信息与计算科学毕业设计论文(1).doc 32页

'差分方程模型的稳定性分析及其应用The Stability Analysis and Application of the Differential Equation Model专业：2011信息与计算科学姓名：郭甜甜指导教师：申请学位级别：学士论文提交日期：2015年5月25日学位授予单位：天津科技大学 摘要本文首先对差分方程这一门课程进行全面深入的研究，了解差分方程的背景，学习差分方程的理论知识，在此基础上对差分方程的稳定性进行学习.并研究相应的数学模型，不仅使这一类常见问题更容易得到解决，更增加了人们的实践经验.差分方程模型作为一种重要的数学模型可以使复杂的生活问题准确、形象地反映出来，并通过对结果的分析对问题进行评估与改善.我主要研究了五个差分模型，分别为金融问题：其一贷款问题研究了欠款，利率，还款额等的关系，其二养老保险问题研究了交保费，保险收益，利率等的关系；减肥计划模型：此模型研究了节食与运动对维持体重的影响关系，制定了减肥方案并给出了维持体重的办法；市场经济中的蛛网模型:通过产品的销售价格和生产产量建立数学模型，在得出市场经济趋于稳定的条件，并且对结果进行分析,最后探讨了在市场经济不稳定时政府能够采用的干预措施；人口控制与预测模型：研究了人口总数的变化状况；军备力量模型：研究了军备力量参与预测战争时的影响.在整个过程中，用MATLAB软件进行计算和画图.关键词：差分方程；稳定性；数学建模；MATLAB ABSTRACTThis paper firstly gives a profound and systematic overview of the course of differential equation, including its background and the related theories, and uses these knowledge as a basis to study the stability of the differential equations. Then the study of the related mathematical model is given, which will not only make solving this kind of common problems easier, but also provide more practical experience for future studies. As an important mathematical model, the differential equation can reflect the complicated problems in people"s life both accurately and vividly. And analyzing the outcome of the differential equation will lead to the evaluation and improvement of the problem. In this paper, I mainly analyze five types of models by using differential equations: the first one is the financial model which can be further divided into two parts --- the loan model which studies the relationship of the debt, the interest rate, the repayment and other related elements and the endowment insurance model which studies the relationship of the premium, the insurance proceeds, the interest rate and other related elements; the second one is the weight-loss plan model which studies the influence of diet and exercise on keeping fit and has created a plan to lose weight and keep fit; the third one is the cobweb model in market economy which is created by taking into account the price and production volume of the product and whose outcome is studied after the condition in which the market economy is heading to stability is achieved, and then discusses about the measures the government can take to enhance its intervention when the market is unstable; the forth one is the population control and prediction model which studies the pattern of the variation in population; the fifth one is the arms race model which studies the impact of arms race on predicting wars. In the whole process, I have used the software --- MATLAB to do calculation and drawings.Key words: differential calculation; stability; creating mathematic models; MATLAB 目录1基础知识11.1差分方程11.2MATLAB介绍31.3数学建模32金融问题模型52.1贷款问题52.2养老保险模型63减肥计划模型93.1问题重述93.2问题分析93.3模型假设93.4符号说明93.5建立模型94蛛网模型114.1问题重述124.2问题分析124.3符号说明124.4蛛网模型124.5差分方程模型14 4.6干预办法154.7模型的推广165人口的预测与控制模型185.1问题重述185.2问题分析185.3建立模型185.4模型的扩展206军备力量模型216.1问题重述216.2问题分析216.3建立模型21结论24参考文献25致谢27 天津科技大学2015届本科生毕业论文1基础知识差分方程表达的为有关离散变量的取值与变换的规律.它是根据所要解决的问题，引进过程中或系统的离散变量，依照实际问题中背景的本质、规律、相关联系，写出离散变量符合的关系等式，进而建立差分方程.得到方程的解后利用分析方程的解，或分析方程的解的某些特性（如稳定性、周期性等），进而明确这些离散变量的变换进程的规律，然后再连同其他分析，从而得到原问题的解.1.1差分方程差分方程的使用范围十分普遍，因为能够使离散变量的逼近与近似来表示连续变量，所以许多模型就可以类似于差分方程模型来解决.所以差分方法既可以在建立离散的数学模型进程中使用，也可以在连续模型化为离散模型的数值计算中广泛的使用.一般来说，但凡涉及到有关变量的规律、本质，便能够使用差分方程模型去表达与分析求解.1.1.1差分方程的概念差分：对于数列，称在处的前向差分为差分算子：.并且称在处的后向差分为差分算子：.本文皆是只前向差分.可知是关于的函数.进而可以定义为处的二阶差分为的差分：，它反映的为量的增量.同理可以定义为处的阶差分.差分方程：由某个函数多个不同时期值的符号或者某个函数的差分组成的方程称为差分方程，其中大多形式为或或通过差分方程的性质和定义能够知道，各种表达形式的差分方程能够互相变换，各自相通.差分方程的解：若将某函数代入差分方程，让方程两边相等，那么就称此函数为差分方程的解，要是差分方程的阶数与此方程的所有解中拥有互相独立的任意常数的个数相同，那么就称此解释差分方程的通解,以便体现在变化过程中某一事物的客观规律性，通常依据此事物在初始时刻所处情况，在差分方程上添加一定的条件，称此为初始条件，24 天津科技大学2015届本科生毕业论文若初始条件确定了通解中任意常数后，此解称之为差分方程的特解[1].1.1.1差分方程常用解法常系数线性差分方程的解方程（1-1）其中是常数，则称方程（1-1）为常系数线性方程.并且称方程（1-2）是方程（1-1）相对应的齐次方程.若（1-2）的解形式为，代近方程中可以得到（1-3）则称方程（1-3）是方程（1-1）和（1-2）的特征方程.可见，只要能够得到方程（1-3）的根，就能够求出方程（1-2）的解.一般结果为：如果方程（1-3）存在个不相同的实根，那么方程（1-2）有通解：如果方程（1-3）存在重根,那么方程（1-2）通解可以表示为：如果方程（1-3）存在两个单复根，记，，,那么方程（1-2）通解可以表示为:如果方程（1-3）存在重复根，记，那么方程（1-2）通解可以表示为:由上可知，由于方程（1-3）恰好有个根，所以方程（1-2）的通解定有个相互独立的任意常数.记方程（1-2）的通解为：,若可以得到方程（1-1）的一个特解：，那么方程（1-1）定有通解：差分方程的变换解法24 天津科技大学2015届本科生毕业论文在差分方程的左右取有关的变换，然后写出的变换，利用的变换，最后利用求解代数方程的方法求出，同时将展开成洛朗级数在解析圆环域里，此系数即为所要求的.1.1.1差分方程稳定性阶常系数线性差分方程（1-1）稳定的充分必要条件为它所相应的特征方程（1-3）所有的特征根满足[1].一阶非线性差分方程（1-4）的平衡点由方程所决定，展开为泰勒形式将在点处，因此：当时，方程（1-4）的解是稳定的.当时，方程（1-4）的平衡点是不稳定的.1.1MATLAB介绍MATLAB为一个面向科学与工程计算的高级语言，一个具有超强能力的数值计算和可视化特点的软件.相对于别的计算机软件，MATLAB的运行方式与人们计算公式时的思考方法非常类似，它编写程序的过程就如同人们在演算纸上罗列出公式进行求解，这避免了较多的重复、繁琐的机械性的编写程序细节，把重点放在有创造性问题上，在最短的时间内得到更具价值的结果.MATLAB具有许多特点，比如功能性强、容易学懂、效率高、应用面广泛、操作简单、节约时间等.MATLAB不仅简单好用，而且数据和图像处理能力很是强大并且可以完成数值分析、管理与调度优化计算、通讯系统设计与仿真、工程与科学绘图等众多功能.现在MATLAB已经演变成为一种大型软件应用在多科学、多工作平台，被各个国家所接收和认可并在一定程度上体现了国际上计算机软件的总体水平，也成为了众多大学生应该熟练掌握的一项基本技能.本文在研究进程中将会使用到制作图像和求解功能.1.2数学建模数学建模是通过数学的知识和思想来简单清晰的表示现实问题中的重点方面，以此完成现实中的问题，即通过使用各种数学办法来完成现实问题.数学建模是一个模拟过程，它是用程序、图像、数学的公式等对实际问题进行抽象、假设、简化后用数学方式描述出来，它可以预料将来的进行情况，可以说明一些客观存在的现象，也可以为有些现象的未来走向提供在特定环境中最合适的方案或相对好的方法.24 天津科技大学2015届本科生毕业论文数学模型建立不但要求对现实问题谨小慎微的分析和观察，而且要求熟练应用各个方面的数学知识.建立数学模型多数应有如下几个阶段：最开始应该明确研究的角色、目标以及问题的类型是确定型还是随机型；把问题简单化后列出将要研究的因素，并把这些因素用参量和变量的方式表现出来；应用数学知识和方法表达出问题中变量之间的联系，一般是列成数学表达式，进而建立了数学模型；通过各种数学知识、数学软件等解出模型的解；把模型的结果转换为与实际问题相适应的清晰易懂的语言；最后进行模型的检测与评估.24 天津科技大学2015届本科生毕业论文1金融问题模型2.1贷款问题2.1.1问题重述由于社会经济的飞速增长，人们生活水平的持续攀升，人们的经济需要更加增多.越来越多的人尤其是工薪阶层需要通过贷款来实现一些经济活动，比如买房、买车、向银行贷款等等.作为贷款人必须清楚贷款的整个运行过程，了解每一个细节尤其是要知道还款方式是等额本息还款法、等额本金还款法或是等本等息等额还款法亦或是其他方法.还应清楚贷款总额，各种还款方式每期应还款额，贷款时间以及贷款利率等.2.1.2问题分析在日常生活中较为常用的就是等额本息还款法，因此在本文中只讨论此贷款方法.等额本息还款法即为每期的所要还的钱数是一定的，而每期所还的本金在逐渐增多，利息越来越少,将通过贷款总值、贷款时间、贷款利率、每期还款额这些因素相互的联系建立模型，再通过数学的递推思想得出第期的欠款额，令欠款额为零时即可得到每期还款额的表达式[6].2.1.3模型假设l贷款期间贷款利率一直不变；l贷款人能如期偿还每期的还款额；l贷款期间不考虑其他的经济问题影响.2.1.4符号说明:贷款总额；:贷款期限（以月计算）；:第个月的欠款额()；:贷款月利率；：每月还款额.2.1.5建立模型等额本息还款法中每月所还的钱应等于每月所还的本金加上每月利息，即为则有第个月还款后欠款额：第1个月还款后欠款额：第2个月还款后欠款额：第3个月还款后欠款额：24 天津科技大学2015届本科生毕业论文...第个月还款后欠款额：应用数学归纳法和等比级数求和公式可得（2-1）当到达最后期限即时，有，带入（2-1）式可得（2-2）(2-2)式即为等额本息还款法中每月还款额.2.1.6举例买一辆11万元的汽车，首付.分12个月还完，年利率为，分别用等额本息和等额本金还款法计算，并进行分析.贷款金额为元；月利率为；等额本息还款法：贷款期限为12个月；每月还款额为元所以还款总额元,其中总利息为3361元.等额本金还款法：还款总额：元其中总利息为2740.24元.由上可知，等额本金还款法所付的利息相对等额本息还款法要少些，并且还款时间越长，利息差值越大.2.2养老保险模型2.2.1问题重述随着社会的不断发展，人们生活水平的持续攀升，人们平均寿命也有所增加，以至于我国慢慢进入老龄化阶段.24 天津科技大学2015届本科生毕业论文为了确保人们老年后的生活有所保障，使得养老保险问题备受关注.现有一保险公司提出了一个养老保险策略，为投保人每月缴费200元一直到59岁末，从60岁开始领取养老金.如果投保人从25岁开始投保，那么60岁以后每月可得2282元养老金，如果投保人从35岁开始投保，那么60岁以后每月可得1056元养老金.2.2.1问题分析本文要研究此保险公司每月至少要有多少投资收益率才能确保保险责任.即保险公司为确保保险人的保险收益必需利用保险人所交的保费最少收获多少利润.通过缴纳的保费和收益的总值，每月收益率，60岁前每月缴费额，60岁后每月领取额，终止缴纳保险费与终止领取养老金的月份之间的关系建立数学模型.2.1.3模型假设l投保人能按期缴纳保险费2.1.4符号说明：截止到第个月所交保费和收益的总额;：每月收益率;：60岁前每月缴费额;：60岁后每月领取额;:停缴保险费的月份;:停领养老金的月份.2.1.5建立模型在全部过程中，可知：(2-3)其中代表的是从投保人开始交保费月后算起的.所要研究的是在第个月时，的数值为多少.若为正数，那么代表保险公司最终获；若为负数，那么代表保险公司最终亏损；若为零，那么代表保险公司最终一无所有，投保人最终获益.2.1.6举例某男子从25岁开始投保，假设男子活到75岁，所以,由（2-3）式可得：24 天津科技大学2015届本科生毕业论文(2-4)在（2-4）式中，分别取,可得设利用MATLAB软件编写代码如下：symsxF=x^600-12.14*x^180+11.41;x=solve(F)由于一定大于1，对众多的根进行分析可得，，即求出每月收益率为：用同样的方法也可求出，35岁开始投保的每月收益率为：24 天津科技大学2015届本科生毕业论文1减肥计划模型3.1问题重述在现代社会中，越来越多的人们尤其女性认为瘦是衡量美的一种重要标准，因此许多自感肥胖的人开始尝试用各种方法减肥，但是减肥药和节食等方法都是存在安全隐患的.专家表明：想要在健康的条件下达到减肥的效果并且维持下去，只有利用控制饮食和进行适当的运动.通常用体重指标（简记）来衡量体重，为体重（千克）除以身高（米）的平方.当时，体重为正常；当时，体重为超重；当时，体重为肥胖.3.2问题分析一般，但凡人体内的能量守恒被破坏就将会导致体重的变化.人们在饮食过程中吸收热量，以至体重增加；人们又通过运动以及代谢消耗热量，以至体重减少.当然减肥的前提是不伤害身体，所以要求每天吸收的热量不能过多，体重减少的也不能过快了.由此就可以通过体重，吸收热量，消耗热量的关系建立数学模型.3.3模型假设（Ⅰ）增加的体重与吸收的热量成正比，每吸收8000千卡热量体重增加1千克,由代谢导致的体重减少与体重成正比，一般一公斤体重每周消耗200千卡到320千卡的热量（每人不同），即为一个70千克的人每天消耗2000千卡到3200千卡的热量[7].（Ⅱ）运动导致的体重减少与体重成正比，并且与运动的时间和形式相关.（Ⅲ）为保证身体的健康，一周内吸收的热量不能小于10000千卡，一周内体重减少不能超过1.5千卡.3.4符号说明：第周末体重；：第周吸收的热量；：热量转换系数；：每小时每千克体重运动消耗的热量（千卡）：每周运动的时间（小时）：代谢消耗系数（因人而异）；：运动消耗系数3.5建立模型可知体重变化的方程为24 天津科技大学2015届本科生毕业论文（3-1）3.5.1减肥计划的提出现为一个具体的人制定减肥计划来探讨此模型的应用.某人高为1.7米，体重为100千克，，现在每周平均吸收20000千卡热量，并保证体重不发生改变.现若让此人体重减到75千克并且保持下去，请依照以下三点制定减肥计划：l当不进行任何运动时计划划分为两个阶段，第一阶段：每周控制饮食慢慢减少吸收的热量，使每周体重减1千克，一直到所吸收热量的最低点（10000千卡）；第二阶段：每周吸收的热量维持在下限，直至达到减肥的目标.l在第二阶段添加运动以加速减肥速度，重新制定第二阶段方案.l制定一个达到目标体重后保持体重的策略.3.5.2减肥计划的制定（Ⅰ）在不进行运动时，可知，已知千卡，千克，（千克/千卡），由（5-1）式可得也就是每周每千克体重消耗千卡的热量.第一阶段：需要每周体重减1千克，一直到所吸收的热量成为最低点（10000千卡），可得带入（5-1)式可得再将带入上式，又因吸收热量的下限为10000千卡，可得说明第一阶段为10周，热量的吸收是依照使得每周体重减少1千克，到第10周末体重减为90千克.第二阶段：每周吸收的热量维持在下限，要将体重减到75千克，由（3-1）可得24 天津科技大学2015届本科生毕业论文（3-2）对（5-2）式进行递推并用等比数列求和可得（3-3）将代入（3-3）可得（3-4）说明第二阶段为19周，在吸收的热量每周维持在10000千卡时，依照减到目标体重75千克.（Ⅱ）依据查询资料可知每小时每千克体重各项运动消耗的热量如下：表5-1各项运动消耗的热量运动跑步跳舞乒乓自行车（中速）游泳热量消耗（千卡）7.03.04.42.57.9在第二阶段添加运动以加速减肥进程，其中，在此取，故，那么（3-4）式中的应改为,则（3-4）式为说明如果在第二阶段增加的运动（如一周骑10小时自行车或跳8小时的舞蹈），那么第二阶段将会减为14周.（Ⅲ）若想达到目标后保持体重，那么要使每一周吸收的热量都维持某常数，并让体重维持不变，由（3-1）式可得可得出：如果不运动，千卡；如果运动，千卡.24 天津科技大学2015届本科生毕业论文1蛛网模型4.1问题重述在处于完全自由的经济市场里，许多商品的销售和生产明显表现出周期性.主要体现在：在一定时期里商品的生产产量、销售价格和销售量是稳定的，所以这些经济数据在某个时期里是离散变量的形式.商品的销售价格和生产产量是最为关注的两个因素，若要做好经营，获得较好的经济效益，必须掌握好这两个经营过程中的最重要的因素.4.2问题分析由于本期产品的销售价格决定于消费者的需求关系，产品数量越少就会导致价格越高.然而下一期产品的数量决定于供应关系，产品的价格越高生产的数量就越多.市场经济中的产品数量与价格产生的振荡决定于这种供求关系.事实上，存在各种形式的振荡，既有可能振幅越来越小直至趋于平稳，也有可能振幅越来越大，此时若没没有外界的干预（如政府）极有可能致使经济崩溃.通过产品的销售价格和生产产量建立数学模型，在得出市场经济趋于稳定的条件，并且对结果进行分析，再探讨政府能够采用的干预措施在市场经济不稳定时.4.3符号说明：第时段产品的数量：第时段产品的数量:平衡点在函数的斜率的绝对值：平衡点在函数的斜率的绝对值4.4蛛网模型将时间离散化划分为若干段，产品的一个生产周期即为一个时段，由于在一个时间段中产品的销售价格由产品产量决定，因此可设：（4-1）它是需求函数，体现的是此商品与消费者的需求关系.由于产品的销售产量与价格成反比，故是单调递减的函数.由于上一个时段的销售价格决定了下一个时段产品的产量，因此可设：或（4-2）24 天津科技大学2015届本科生毕业论文为的反函数，它们都是供应函数，体现的是生产者的供应关系.由于本时段价格与下时段生产产量成正比，故是单调递增的函数.通过函数和反映和的变化过程，把点列和利用对应的几何关系画出来，即将点列连接起来（见图4-1），则将连成折线形似蛛网，因此这种用图形来研究市场经济的稳定性称为蛛网模型.图4-1图4-2可见，若点列最终收敛于点，即而且点是函数和的交点，则代表市场经济在未来的一段时间里将会趋向稳定.若没有收敛于一点（见图4-2），则代表市场经济将会趋向不稳定.通常，是由消费者的消费能力和需求程度决定的，是由生产者的经营能力和生产能力等因素决定的[8].通过分析图形可知：当时，点是稳定的；当时，点是不稳定的.举例说明蛛网模型：设：产品的本期产品数量由上期的销售价格决定，那么供给函数是，产品本期的需求量由本期产品销售价格，那么需求函数是，那么结合动态供需均衡模型，蛛网模型可以表达为：其中都是正值.由上述三式可得：24 天津科技大学2015届本科生毕业论文（4-3）因此能够知道第期的产品价格是：由于市场是均衡的，故有均衡价格，带入（4-3）式得，将其带入上式有（4-4）对（4-4）式进行分析可得：l当时，那么，称为收敛型蛛网；l当时，那么，称为发散型蛛网；l当时，那么是常数，称为封闭型蛛网.4.5差分方程模型分别取函数和在点附近的近似曲线，可得：（4-3）（4-4）24 天津科技大学2015届本科生毕业论文将（4-3）和（4-4）合并后能得：（4-5）对（4-5）进行递推可得：（4-6）由（4-6）可得，当时，则当或时点稳定；当时，则当或时点稳定；由于是点在上的切线斜率，是点在上的切线斜率，则有，可见差分方程模型与蛛网模型结果是相同的.从（4-3）可得，的意义是产品的数量下降一单位时销售价格的上升幅度，因此代表的是购买者对产品需要的灵敏度，若是生活必需的产品，并且消费者的状态是持币待购，一旦产品的数量缺少，人们就会抢购，则相对较大.的意义为这期销售价格上升一单位是产品数量的增加量，因此代表生产者对产品价格的灵敏度，若生产者贪图当下的高利润，一旦价格上升就增多生产，则相对较大.4.6干预办法综上可知，当一定时，越小，代表购买者对产品需要的灵敏度就越小，越对经济稳定有利;当一定时，越小，代表生产者对产品价格的灵敏度就越小，越对经济稳定有利.相反的，当，越大时，越对经济稳定不利.图4-3图4-4存在两种干预办法在市场经济倾向不稳定时，第一种是让尽可能小，为了更加明显研究的情况，也就是的图像为水平直线（见图4-3），此刻市场经济永远是稳定的无论如何变化（也就是无论多大）.24 天津科技大学2015届本科生毕业论文现实中就相当于控制价格不能变化，不管产品数量为多少，即政府控制物价.第二种是让尽可能小，为了更加明显研究的情况，也就是的图像为竖直直线（见图4-4），此刻市场经济永远是稳定的无论如何变化（也就是无多大）.现实中就相当于不管产品的价格为多少，产品数量不能变化，当供不应求时将从其他地方购买或调货过来，当供应多于需要时，收购多于部分.4.7模型的推广为了更加谨慎生产者在计算下一期的产品数量时，不但考虑这期的销售价格也考虑前一期的销售价格，则（4-2）式将变为：（4-7）（4-2）式的近似直线（4-4）相应的改为：（4-8）由于（4-1）式和（4-3）式没有变化，所以合并（4-3）式和（4-8）式可得：（4-9）只要方程的特征根都在单位圆里，那么当时，即点稳定.（4-9）式的特征方程为：（4-10）并得出（4-10）的特征根为（4-11）当时，有因而，故不在单位圆内，所以舍去.当时，可由（4-11）式得：如果让所有特征根在单位圆里，也就是，所有24 天津科技大学2015届本科生毕业论文（4-12）（4-12）式即为点稳定的条件.与之前点稳定的条件相比，这个模型的的使用范围都放宽了，即稳定性条件变宽了.若要更深一步的研究这个模型，在计算下一期的产品数量时，可考虑最近三年来的价格，即.24 天津科技大学2015届本科生毕业论文1人口的预测与控制模型5.1问题重述随着生活水平的不断提高，人口问题备受人们的关注.人口的均衡发展是与每一个国家息息相关的，有些国家的出生率过高，持续发展下去会对人们的正常生活产生很大的影响；又有些国家的自然增长率接近于零更有甚者成为负数，对国家的发展也很是不利的，如缺少劳动力.我国在人口问题上就有很大问题，不但人口总数增长率太高，并且人口老龄化问题日益明显.所以现在最重要的任务就是在控制人口增长率的情况下有效的控制人口老龄化问题，尽量将年龄结构调节成适当的水平.5.2问题分析每个地方人口数量的变化都是由多种因素所决定的，例如国家政策，社会发展水平，迁移，自然灾害等.本文主要研究性别比例，女性生育率与死亡率对人口数量的影响.通过总人数，死亡率，性别比例，生育率的关系建立差分模型.5.3建立模型为了便捷定义时间单位为年，即将时间离散化.设第年为岁的人数是(其中岁为最大年龄)；第年为岁的人口死亡率是，也就是第年为岁的人中死亡总数与总人数的比：由此可得到：第年为岁女性的生育率是，生育区间是；第年为岁的人口女性比是，也就是第年为岁的人中女性总数与第年总人数的比.可见，在第年出生的总人数是：其中，第年婴儿存活下来的数量为24 天津科技大学2015届本科生毕业论文有代表岁女性的总生育率，那么有，若假定女性生育率一直维持不变，那么：可知：代表平均生育孩子的数量每个女性一生中，体现了人口数量变化的基本因素，定义为总和生育率.综上可得：记可得：并且有以便深入的表达人口数量的状况，用向量说明此问题.用代表人口的分布向量，即记存活率矩阵为：24 天津科技大学2015届本科生毕业论文记生育模式矩阵为：那么可知即这里是状态变量，是控制变量.在实际应用中，可以通过查阅资料了解人口最原始的分布与存活率矩阵,同时知道了生育模式矩阵,那么再知道总和生育率后就可以预测未来的人口数，若想控制未来人口数就通过改变总和生育率.5.4模型的扩展以便于更加全面的表达人口的相关因素，还可以加入以下因素：人口总数：平均年龄：平均寿命：利用这三个因素与关系，在原有的模型上进行扩展，对此模型更加具体科学的分析.24 天津科技大学2015届本科生毕业论文1军备力量模型6.1问题重述历史证明武器会触发一个国家宣布战争.大型军械库的存在会增加暴力冲突的可能性.如果没有破坏性武器，有时国与国之间也许会以其他方式解决争端.因此对军备竞赛进行研究是至关重要的，军备竞赛即为在不战争时期与对立国家或者可能成为对立国家互相当成假想敌，进行质量与数量在军备部分的竞争.每个国家考虑到将来有可能触发的战争，扩大军备，提高军事水平.6.2问题分析一个国家的军备力量是由多方面所决定的，本文中将通过理查森军备竞赛模型来建立数学模型对军备力量进行研究.军备竞赛往往是战争的前兆.如果两个国家都增加国防支出，那么即使是一件小事也会触发战争.而如果两国减少国防支出，小事不会触发战争.因此将用一个具体的例子用此模型来检验这场军备竞赛是“稳定”还是“不稳定”.6.3建立模型现假设国在时间的军备力量为，国从时到时的军备力量的变化表示为：这个模型同样适用于国,国在时间的军备力量为,国从时到时的军备力量的变化表示为：考虑到一国军备力量对另一国的影响，需引入对方军备力量对自己影响程度的相关系数，即有：再考虑到自己以往军备力量对现在的影响，需引入自身制约程度的系数，即有：最后考虑到其他因素，如管理者的野心，政治格局等，需引入一个常量，即有：（6-1）24 天津科技大学2015届本科生毕业论文（6-2）仔细思考表6-1中的数据，该表显示的是伊拉克和伊朗在1975年战争之前的军备力量.表中的数据是两国从1954年到1974年的军备费用/国防支出.表6-1军备费用/国防支出年费伊朗伊拉克195478751955107671956126941957151102195824311019592711296196029214519613201851962345206196338727119644253591965435402196646045019674734801968498513196953454919706127231971732781197284092119739801292197413081632利用表6-1中的数据计算出（6-1）式和（6-2）式中的系数，可以通过MATLAB软件编写代码如下：data=[7875;10767;12694;151102;243110;271129;292145;320185;345206;387271;425359;435402;460450;473480;498513;534549;612723;732781;840921;24 天津科技大学2015届本科生毕业论文9801292;13081632];>>x=data(2:end,1);>>y=data(1:end-1,:)y=787510767126941511022431102711292921453201853452063872714253594354024604504734804985135345496127237327818409219801292>>LinearModel.fit(y,x)ans=Linearregressionmodel:y~1+x1+x2EstimatedCoefficients:EstimateSEtStatpValue(Intercept)37.06326.3531.40640.17761x10.650770.16513.94160.0010523x20.431690.120413.58530.0022805因此可得：24 天津科技大学2015届本科生毕业论文(6-3)同理可得：(6-4)使用线性代数，可以解出系统稳定性.现将（6-3）式和（6-4）式化为矩阵式，模型如下：（6-5）现用向量表示，向量表示，该模型可以简化为：找出并使用特征值和特征向量，解出（6-5）式齐次部分的通解，得出以下的解：用公式解出（6-5）式非齐次的部分特解,如下：综上，最终的通解是由此可知当时，(1.266798)这项的增长是趋于无限的.因此这个系统不稳定，将会导致战争.一场稳定的军备竞赛中至少存在一个或两个国家都满意的平衡点.双方都没必要加强军备力量，超出这个平衡点.军备竞赛中的平衡点代表了军备力量的水平，所以这时军备竞赛停止.由于军备力量的价值，所以两国都没必要改变军备力量.而一场不稳定的军备竞赛中不存在平衡点.军备费用持续增加，破坏性武器的数量也持续升高.军备竞赛的力度加大，一国军备力量加大后，另一国也会随之加大.这种不稳定的状态下，极有可能迅速爆发战争.24 天津科技大学2015届本科生毕业论文结论差分方程是数学中重要的一部分，因为能够用离散变量的逼近和近似来表示连续变量，所以许多模型就可以类似于差分方程模型来解决.所以差分方法既可以在建立离散的数学模型过程应用，也可以在连续模型化为离散模型的数值计算中广泛的应用.一般来说，只要涉及到有关变量的规律、本质，就可以适当地用差分方程模型来表达与分析求解.本文中的五个模型在实际问题中应用十分广泛，对于问题的解决更加方便快捷.金融问题是与生活是密切相关的，通过贷款模型可以算出实际贷款时每月所应支付的金额，清楚总够支付了多少利息；养老保险模型能够知道方案是否划算.减肥是现代女性最为关注的问题，减肥模型给出了健康的减肥计划和维持的方案.在市场中价格与产量是有着密不可分的关系，蛛网模型清晰的反映了二者之间的关系，给出了经济趋于稳定的条件，并讨论了在不稳定时的调整对策，在市场经济中十分受用.我国作为人口大国，对人口的研究是至关重要的，本文中探究了一些因素对人口数量的影响.每个国家的军备竞赛是必不可少的，本文给出了与军备力量相关的因素以及计算方式.2424 参考文献[1]吴传生.经济数学微积分第二版[M].北京：高等教育与出版社，2014:406-410.[2]杨德平.赵维加.管殿柱等.MATLAB基础教程[M].北京：机械工业出版社，2013：193-196.[3]姜启源.数学建模第四版[M].北京：高等教育与出版社，2011:5-26.[4]颜士新.高等数学在住房贷款中的应用[N].钦州学院学报，2007,6:21-25.{5]贺晓东.差分方程在数学建模中的应用[J].中国科技博览，2012,37:100—101.[6]张瑞.张新洋.商业养老保险的精算模型[N].郑州大学学报，2006,4:117-120.[7]项海飞.分阶段减肥的数学模型[J].内江科技，2007,5:54-55.[8]白彦峰.叶菲.基于蛛网模型的我国小宗农产品价格调控研究[J].创新，2011,3:41-46.[9]张婧.数学建模的几种常用方法[N].太原大学教育学院报，2012，3:38-40.[10]姜辉.我国人口预测分析[J].科技管理研究，2005，11:142—145[11]徐新荣.差分方程数学建模分析[J].中国科技信息，2012,14:46-46.[12}张功盛.康光清.差分方程在数学建模中的几个应用实例[N].江西电力职业技术学院学报，2009,1:75-77.[13]D.D.Mooney,R.J.Swift.AcourseinMathenamaticalModeling.TheMathematicalAcssociationofAmerica,2006.[14]S.Chatterjee,A.S.Hadi,B.Price.RegressionAnalysisByExample(thirdedition).JohnWiley&Sons,2000.[15]Chiarell,C.TheCobwebModel:ItsInstabilityandTheOnsetofChaos.EconomicModelling,1988(8:377-384). 致谢今年六月份的到来代表着我的科大校园时光就要结束了.在这四年里我哭过笑过疯狂过，挥洒着青春，洋溢着朝气，学习了知识，收获了友谊，遇到了良师.在这即将分别的时刻我不想去感伤只想去感恩，感恩大学校园的一切美好的人与物，教师孜孜不倦的教诲、同窗、朋友的关怀与关心.这篇论文是在刘老师孜孜不倦的教诲中写完的，使我在知识的海洋中汲取了更多的营养，从而让我更加深刻了解了数学.从选题、搜集资料、制定纲领、中期预备、确定初稿到定稿的每一个细节，刘老师都给了我很多意见与建议，常常帮我疏通思路，解答不懂的问题，主动关心我们的论文进展并对论文提出宝贵意见，在刘老师的帮助下我最终完成了本论文.借此机会对刘寅立老师表达最诚心的谢意和无限的敬佩.在此，我也要感谢在本次书写论文的过程中过去向我提供过帮助的教师、家人和同窗，你们使我论文的完成的更加顺利优秀.我会用未来的成就回报每一个爱我和向我提供过帮助的人.'