数学建模例题1
2016年数学建模论文
第 1 套
论文题目：人口增长模型的确定组别：54
姓名：俞嘉艺吕游姜飞龙提交日期：2016.7.4
人口增长模型的确定
摘要
本文根据某地区的人口统计，建立模型该地区1980年后每隔十年预测五次人数量。首先，通过直接观察人口的变化规律后，我们假设该地区人口是时间的指数模型，建立一个指数模型，并用最小二乘法进行数据拟合，得到数据的具体参数，从而对人口数量进行预测。然后我们发现从1880年以后人口增长。。。。。。
关键字：人口预测指数函数模型
一、问题重述
1790-1980年间美国每隔10年的人口记录如下表所示。
1.试用以上数据建立马尔萨斯(Malthus)人口二次函数增长模型，并对接下来的每隔十年预测五次人口数量，并查阅实际数据进行比对分析。
2.如果数据不相符，再对以上模型进行改进，寻找更为合适的模型进行预测，并对两次预测结果进行对比分析。
3.查阅资料找出中国人口与表1同时期的人口数量，用以上建立的两个模型进行人口预测与分析。
二、变量说明
X(t) t时刻的人口数量
初始时刻的人口数量
r 人口增长率
环境所能容纳的最大人口数量
三、问题分析
首先，我们运行matlab软件编程（附件一），绘制出1790年到1980年的的人口数据图，如图所示：……
图1
结论
四、模型建立
模型以：二次函数模型，即：
我们假设该地区t时刻的人口数量X(t)是时间的儿媳函数，我们柑橘最小二乘法，利用已有数据拟合得到具体的参数，即，a/b/c，使得以下函数达到最小值：
表示该地区时间的人口数
五、模型求解
令, , ,可得到关于a/b c的一次方程，用matlab编程得到，即
a=….b=……c=……
二次函数模型为：。。。。。。
图。。。
六、结果分析
结论。。。。。。。
七、参考文献
[1]刘卫国，陈兆平。MATLAB程序设计与应用M,北京高等教育出版社，2002年。
[2] 姜美元，谢金星。数学建模M,,北京；高等教育出版社，2004年符箓（程序）
数学建模典型例题
一、人体重变化 某人的食量是10467焦/天，最基本新陈代谢要自动消耗其中的5038焦/天。每天的体育运动消耗热量大约是69焦/（千克?天）乘以他的体重（千克）。假设以脂肪形式贮存的热量100% 地有效，而1千克脂肪含热量41868焦。试研究此人体重随时间变化的规律。 一、问题分析 人体重W（t）随时间t变化是由于消耗量和吸收量的差值所引起的，假设人体重随时间的变化是连续变化过程，因此可以通过研究在△t时间内体重W的变化值列出微分方程。 二、模型假设 1、以脂肪形式贮存的热量100%有效 2、当补充能量多于消耗能量时，多余能量以脂肪形式贮存 3、假设体重的变化是一个连续函数 4、初始体重为W0 三、模型建立 假设在△t时间内： 体重的变化量为W（t+△t）-W(t)； 身体一天内的热量的剩余为（10467-5038-69*W（t）） 将其乘以△t即为一小段时间内剩下的热量； 转换成微分方程为：d[W(t+△t)-W(t)]=(10467-5038-69*W(t))dt； 四、模型求解 d(5429-69W)/(5429-69W)=-69dt/41686 W(0)=W0 解得： 5429-69W=（5429-69W0）e(-69t/41686) 即： W（t）=5429/69-（5429-69W0）/5429e(-69t/41686) 当t趋于无穷时，w=81; 二、投资策略模型 一、问题重述 一家公司要投资一个车队并尝试着决定保留汽车时间的最佳方案。5年后，它将卖出所有剩余汽车并让一家外围公司提供运输。在策划下一个5年计划时，这家公司评估在年i 的开始买进汽车并在年j的开始卖出汽车，将有净成本a ij（购入价减去折旧加上运营和维修成本）ij
数学建模1例题解析
1.贷款问题 小王夫妇计划贷款20万元购买一套房子，他们打算用20年的时间还清贷款。目前，银行的利率是%/月。他们采用等额还款的方式(即每月的还款额相同)偿还贷款。 (1)在上述条件下，小王夫妇每月的还款额是多少共计付了多少利息 (2)在贷款满5年后，他们认为他们有经济能力还完余下的款额，打算提前还贷，那么他们在第6年初，应一次付给银行多少钱，才能将余下全部的贷款还清 (3)如果在第6年初，银行的贷款利率由%/月调到%/月，他们仍然采用等额还款的方式，在余下的15年内将贷款还清，那么在第6年后，每月的还款额应是多少 (4)某借贷公司的广告称，对于贷款期在20年以上的客户，他们帮你提前三年还清贷款。但条件是： (i)每半个月付款一次，但付款额不增加，即一次付款额是原付给银行还款额的1/2； (ii)因为增加必要的档案、文书等管理工作，因此要预付给借贷公司贷款总额10%的佣金。 试分析，小王夫妇是否要请这家借贷公司帮助还款。 解答： (1)贷款总月数为N=20*12=240，第240个月的欠款额为0，即。 利用式子 (元)，即每个月还款元，共还款(元)，共计付利息元。 (2)贷款5年(即5*12=60个月)后的欠款额为， 利用公式：， 所以，
(元) (3)元，即第六年初，贷款利率，所以余下的15年，每个月还款额为：(元) (4)按照借贷公司的条件(i)每半个月付款一次，但付款额不增加，即一次付款额是原付给银行还款额的，付款的时间缩短，但是前17年的付款总额不变。帮忙提前三年还清需要资金数： 。 对于条件(ii)佣金数： 分析：因为预付佣金20000元，按照银行存款利率/月，17年的存款本息为 即在第17年需要给付借贷公司的钱少于给付银行的钱。所以建议请这家借贷公司帮助还款。 2.冷却定律与破案 按照Newton冷却定律，温度为T的物体在温度为的环境中冷却的速度与温差成正比。用此定律建立相应的微分方程模型。 凌晨某地发生一起凶杀案，警方于晨6时到达案发现场，测得尸温26℃，室温10℃，晨8时又测得尸温18℃。若近似认为室温不变，估计凶杀案的发生时间。 解答： 根据Newton冷却定律，可知温度T的微分方程为：
19191-数学建模-3.1
微分方程模型 浙江大学数学建模实践基地
§3.1 微分方程的几个简单实例 在许多实际问题中，当直接导出变量之间的函数关系较为困难，但导出包含未知函数的导数或微分的关系式较为容易时，可用建立微分方程模型的方法来研究该问题， 本节将通过一些最简单的实例来说明微分方程建模的一般方法。在连续变量问题的研究中，微分方程是十分常用的数学工具之一。
例1（理想单摆运动）建立理想单摆运动满足的微 分方程，并得出理想单摆运动的周期公式。 从图3-1 中不难看出，小球所受的合力为mgsin θ，根据牛顿第二定律可得：sin ml mg θ θ=-从而得出两阶微分方程：0sin 0(0)0,(0)g l θθθθθ+==?=????（3.1）这是理想单摆应满足的运动方程 （3.1）是一个两阶非线性方程，不 易求解。当θ很小时，sin θ≈θ，此时，可 考察（3.1）的近似线性方程： 0(0)0,(0)g l θθθθθ+==?=?? ??（3.2）由此即可得出2g T l π=（3.2）的解为: θ(t )=θ0cosωt g l ω=其中当时,θ(t )=04T t =42g T l π =故有M Q P mg θl 图3-1 （3.1）的 近似方程
例2我方巡逻艇发现敌方潜水艇。与此同时敌方潜水艇也发现了 我方巡逻艇，并迅速下潜逃逸。设两艇间距离为60哩，潜水艇最大航速为30节而巡逻艇最大航速为60节，问巡逻艇应如何追赶潜水艇。 这一问题属于对策问题，较为复杂。讨论以下简单情形：敌潜艇发现自己目标已暴露后，立即下潜，并沿着直线方向全速逃逸，逃逸方向我方不知。 设巡逻艇在A 处发现位于B 处的潜水艇，取极坐标，以B 为极点，BA 为极轴，设巡逻艇追赶路径在此极坐标下的方程为r =r (θ)，见图3-2。 B A A1 dr ds dθ θ图3-2 由题意，，故ds =2dr 2ds dr dt dt =图3-2可看出， 2 2 2 ()()()ds dr rd θ=+
数学建模典型例题(二)
6 小行星的轨道模型 问题 一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立以太阳为原点的直角坐标系，在两坐标轴上取天文测量单位（一天文单位为地球到太阳的平均距离：1.4959787×1011m ）.在5个不同的时间对小行星作了5次观察，测得轨道上5个点的坐标数据如表6.1. 表6.1 坐标数据 由Kepler （开普勒）第一定律知，小行星轨道为一椭圆.现需要建立椭圆的方程以供研究（注：椭圆的一般方程可表示为 012225423221=+++++y a x a y a xy a x a . 问题分析与建立模型 天文学家确定小行星运动的轨道时，他的依据是轨道上五个点的坐标数据： （x 1, y 1）, （x 2, y 2）, （x 3, y 3）, （x 4, y 4）, （x 5, y 5）. 由Kepler 第一定律知，小行星轨道为一椭圆.而椭圆属于二次曲线，二次曲线的一般方程为012225423221=+++++y a x a y a xy a x a .为了确定方程中的五个待定 系数，将五个点的坐标分别代入上面的方程，得 ???? ?????-=++++-=++++-=++++-=++++-=++++.122212221222122212225554253552251454424344224 135342 3333223125242 232222211514213112211y a x a y a y x a x a ， y a x a y a y x a x a ，y a x a y a y x a x a ，y a x a y a y x a x a ，y a x a y a y x a x a 这是一个包含五个未知数的线性方程组，写成矩阵
数学建模题目及答案
09级数模试题 1. 把四只脚的连线呈长方形的椅子往不平的地面上一放，通常只有三只脚着地，放不稳，然后稍微挪动几次，就可以使四只脚同时着地，放稳了。试作合理的假设并建立数学模型说明这个现象。（15分） 解：对于此题，如果不用任何假设很难证明，结果很 可能是否定的。 因此对这个问题我们假设： （1）地面为连续曲面 （2）长方形桌的四条腿长度相同 （3）相对于地面的弯曲程度而言，方桌的腿是足够长的 （4）方桌的腿只要有一点接触地面就算着地。 那么，总可以让桌子的三条腿是同时接触到地面。 现在，我们来证明：如果上述假设 条件成立，那么答案是肯定的。以长方 桌的中心为坐标原点作直角坐标系如图 所示，方桌的四条腿分别在A、B、C、D 处，A、、D的初始位置在与x轴平行，再 假设有一条在x轴上的线,则也与A、B，C、D平行。当方桌绕中心0旋转时，对角线与x轴的夹角记为θ。 容易看出，当四条腿尚未全部着地时，腿到地面的距离是不确定的。为消除这一不确定性，令() fθ为A、B离地距离之和，
()g θ为C 、D 离地距离之和，它们的值由θ唯一确定。由假设（1）， ()f θ,()g θ均为θ的连续函数。又由假设（3） ，三条腿总能同时着地， 故()f θ()g θ=0必成立（?θ）。不妨设(0)0f =(0)0g >（若(0)g 也为0，则初始时刻已四条腿着地，不必再旋转），于是问题归结为： 已知()f θ,()g θ均为θ的连续函数，(0)0f =,(0)0g >且对任意θ有00()()0f g θθ=，求证存在某一0θ，使00()()0f g θθ=。 证明：当θ=π时，与互换位置，故()0f π>，()0g π=。作()()()h f g θθθ=-，显然，()h θ也是θ的连续函数，(0)(0)(0)0h f g =-<而()()()0h f g πππ=->，由连续函数的取零值定理，存在0θ，00θπ<<，使得0()0h θ=，即00()()f g θθ=。又由于00()()0f g θθ=，故必有00()()0f g θθ==，证毕。 2.学校共1000名学生，235人住在A 宿舍，333人住在B 宿舍，432人住在C 宿舍。学生 们要组织一个10人的委员会，试用合理的方法分配各宿舍的委员数。（15分） 解：按各宿舍人数占总人数的比列分配各宿舍的委员数。设：A 宿舍的委员数为x 人，B 宿舍的委员数为y 人，C 宿舍的委员数为z 人。计算出人数小数点后面的小数部分最大的整数进1，其余取整数部分。 则 10； 10=235/1000；
2003全国大学生数学建模竞赛B题优秀论文(出题人亲作)
2003高教社杯全国大学生数学建模竞赛 B 题参考答案 注意：以下答案是命题人给出的，仅供参考。各评阅组应根据对题目的理解及学生的解答，自主地进行评阅。 问题分析： 本题目与典型的运输问题明显有以下不同： 1． 运输矿石与岩石两种物资； 2． 产量大于销量的不平衡运输； 3． 在品位约束下矿石要搭配运输； 4． 产地、销地均有单位时间的流量限制； 5． 运输车辆每次都是满载，154吨/车次； 6． 铲位数多于铲车数意味着最优的选择不多于7个产地； 7． 最后求出各条路线上的派出车辆数及安排。 运输问题对应着线性规划，以上第1、2、3、4条可通过变量设计、调整约束条件实现； 第5条使其变为整数线性规划；第6条用线性模型实现的一种办法，是从1207 10 C 个整数规划中取最优的即得到最佳物流；对第7条由最佳物流算出各条路线上的最少派出车辆数（整数），再给出具体安排即完成全部计算。 对于这个实际问题，要求快速算法，计算含50个变量的整数规划比较困难。另外，这是一个二层规划，第二层是组合优化，如果求最优解计算量较大，现成的各种算法都无能为力。于是问题变为找一个寻求近优解的近似解法，例如可用启发式方法求解。 调用120次整数规划可用三种方法避免：（1）先不考虑电铲数量约束运行整数线性规划，再对解中运量最少的几个铲位进行筛选；（2）在整数线性规划的铲车约束中调用sign 函数来实现；（3）增加10个0－1变量来标志各个铲位是否有产量。 这是一个多目标规划，第一问的目标有两层：第一层是总运量（吨公里）最小，第二层是出动卡车数最少，从而实现运输成本最小。第二问的目标有：岩石产量最大；矿石产量最大；运量最小，三者的重要性应按此序。 合理的假设主要有： 1. 卡车在一个班次中不应发生等待或熄火后再启动的情况； 2. 在铲位或卸点处因两条路线（及以上）造成的冲突时，只要平均时间能完成任务即 可，不进行排时讨论； 3. 空载与重载的速度都是28km/h ，耗油相差却很大，因此总运量只考虑重载运量； 4. 卡车可提前退出系统。 符号：x ij ~ 从i 号铲位到j 号卸点的石料运量 单位 吨； c ij ~ 从i 号铲位到j 号卸点的距离 公里； T ij ~ 从i 号铲位到j 号卸点路线上运行一个周期平均所需时间 分； A ij ~ 从i 号铲位到j 号卸点最多能同时运行的卡车数 辆； B ij ~ 从i 号铲位到j 号卸点路线上一辆车最多可以运行的次数 次； p i ~ i 号铲位的矿石铁含量。 % p =(30,28,29,32,31,33,32,31,33,31) q j ~ j 号卸点任务需求 吨 q =(1.2,1.3,1.3,1.9,1.3)*10000
数学建模例题及解析
。 例１差分方程—-资金的时间价值 问题1：抵押贷款买房——从一则广告谈起 每家人家都希望有一套(甚至一栋)属于自己的住房，但又没有足够的资金一次买下，这就产生了贷款买房的问题。先看一下下面的广告（这是1991年1月1日某大城市晚报上登的一则广告）,任何人看了这则广告都会产生许多疑问,且不谈广告中没有谈住房面积、设施等等,人们关心的是：如果一次付款买这栋房要多少钱呢?银行贷款的利息是多少呢？为什么每个月要付12０0元呢？是怎样算出来的?因为人们都知道,若知道了房价(一次付款买房的价格），如果自己只能支付一部分款,那就要把其余的款项通过借贷方式来解决,只要知道利息，就应该可以算出五年还清每月要付多少钱才能按时还清贷款了,从而也就可以对是否要去买该广告中所说的房子作出决策了。现在我们来进行数学建模。由于本问题比较简单无需太多的抽象和简化。 a。明确变量、参数，显然下面的量是要考虑的： 需要借多少钱，用记; 月利率(贷款通常按复利计）用R记； 每月还多少钱用x记; 借期记为N个月。 b.建立变量之间的明确的数学关系。若用记第k个月时尚欠的款数,则一个月后（加上利息后)欠款 , 不过我们又还了x元所以总的欠款为 ｋ=0，1，２，３， 而一开始的借款为.所以我们的数学模型可表述如下 (1) c. （1)的求解。由
（2）这就是之间的显式关系。 d.针对广告中的情形我们来看（１）和（2）中哪些量是已知的。Ｎ=5年=60个月，已知;每月还款x＝１20０元，已知Ａ.即一次性付款购买价减去７０000元后剩下的要另外去借的款，并没有告诉你，此外银行贷款利率Ｒ也没告诉你,这造成了我们决策的困难.然而,由（2）可知6０个月后还清,即,从而得 (3） A和x之间的关系式，如果我们已经知道银（３)表示Ｎ=6０,x=1２00给定时0 A。例如，若Ｒ=0．01,则由(3)可算得行的贷款利息R，就可以算出0 5３94６元。如果该房地产公司说一次性付款的房价大于70０0０十53946＝12３946元的话，你就应自己去银行借款。事实上，利用图形计算器或Maｔheｍａtiｃa这样的 数学软件可把（3）的图形画出来,从而可以进行估算决策。以下我们进一步考虑下面两个问题。 注1问题1标题中“抵押贷款”的意思无非是银行伯你借了钱不还，因而要你用某种不动产（包括房子的产权）作抵押，即万一你还不出钱了，就没收你的不动产。 例题１某高校一对年青夫妇为买房要用银行贷款６000０元，月利率0.01，贷款期25年=300月，这对夫妇希望知道每月要还多少钱，25年就可还清。假设这对
全国数学建模大赛题目
2010高教社杯全国大学生数学建模竞赛题目 A题储油罐的变位识别与罐容表标定 通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。 请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为α=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 （2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度α和横向偏转角度β）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。 附件1：小椭圆储油罐的实验数据 附件2：实际储油罐的检测数据 地平线油位探针
数学建模优化问题经典练习
1、高压容器公司制造小、中、大三种尺寸的金属容器，所用资源为金属板、劳 万元，可使用的金属板有500t，劳动力有300人/月，机器有100台/月，此外，不管每种容器制造的数量是多少，都要支付一笔固定的费用：小号为100万元，中号为150万元，大号为200万元，现在要制定一个生产计划，使获得的利润为最大， max=4*x1+5*x2+6*x3-100*y1-150*y2-200*y3; 2*x1+4*x2+8*x3<=500; 2*x1+3*x2+4*x3<=300; 1*x1+2*x2+3*x3<=100; @bin(y1); @bin(y2); @bin(y3); y1+y2+y3>=1; Global optimal solution found. Objective value: 300.0000 Extended solver steps: 0 Total solver iterations: 0 Variable Value Reduced Cost X1 100.0000 0.000000 X2 0.000000 3.000000 X3 0.000000 6.000000 Y1 1.000000 100.0000 Y2 0.000000 150.0000 Y3 0.000000 200.0000 Row Slack or Surplus Dual Price 1 300.0000 1.000000 2 300.0000 0.000000 3 100.0000 0.000000 4 0.000000 4.000000 5 0.000000 0.000000
数学建模试题
2012-2013第一学期 《数学建模》试题卷 班级：2010级统计 姓名：石光顺 学号：20101004025 成绩：
一、用Matlab 求解以下优化问题（10分） 用Matlab 求解下列线性规划问题： 解：首先化Matlab 标准型，即 123min 3w x x x =-++ 123121114123x x x ?? -??????≤??????---???? ???? ， [][]1 2 32011T x x x -?= 然后编写Matlab 程序如下： f=[-3,1,1]; a=[1,-2,1;4,-1,-2]; b=[11,-3]; aeq=[-2,0,3]; beq=1; [x,y]=linprog(f,a,b,aeq,beq,zeros(3,1)); x,y=-y 运行结果： x = 0.0000 2.3333 0.3333 y = -2.6667 即当1230, 2.3333,0.3333x x x ===时，max 2.6667z =-。
二、求解以下问题，列出模型并使用Matlab求解（20分） 某厂生产三种产品I，II，III。每种产品要经过A, B两道工序加工。设该厂有两种规格的设备能完成A工序，它们以A1, A2表示；有三种规格的设备能完成B工序，它们以B1, B2, B3表示。产品I可在A, B任何一种规格设备上加工。产品II可在任何规格的A设备上加工，但完成B工序时，只能在B1设备上加工；产品III 只能在A2与B2设备上加工。已知在各种机床设备的单件工时，原材料费，产品销售价格，各种设备有效台时以及满负荷操作时机床设备的费用如表1，求安排最优的生产计划，使该厂利润最大。 表1 解：（1）根据题意列出所有可能生产产品I、II、III的工序组合形式，并作如下假设： x ; 按(A1，B1)组合生产产品I，设其产量为 1 x; 按(A1，B2)组合生产产品I，设其产量为 2 x; 按(A1，B3)组合生产产品I，设其产量为 3 x; 按(A2，B1)组合生产产品I，设其产量为 4 x; 按(A2，B2)组合生产产品I，设其产量为 5
matlab数学建模实例
第四周 3. 中的三个根。 ，在求8] [0,041.76938.7911.1-)(2 3=-+=x x x x f function y=mj() for x0=0:0.01:8 x1=x0^3-11.1*x0^2+38.79*x0-41.769; if (abs(x1)<1.0e-8) x0 end end 4.分别用简单迭代法、埃特金法、牛顿法求解方程，并比较收敛性与收敛速度（ε分别取10-3、10-5、10-8）。 简单迭代法： function y=jddd(x0) x1=(20+10*x0-2*x0^2-x0^3)/20; k=1; while (abs(x1-x0)>=1.0e-3) x0=x1; x1=(20+10*x0-2*x0^2-x0^3)/20;k=k+1; end x1 k 埃特金法： function y=etj(x0) x1=(20-2*x0^2-x0^3)/10; x2=(20-2*x1^2-x1^3)/10; x3=x2-(x2-x1)^2/(x2-2*x1+x0); k=1; while (abs(x3-x0)>=1.0e-3) x0=x3; x1=(20-2*x0^2-x0^3)/10; x2=(20-2*x1^2-x1^3)/10; x3=x2-(x2-x1)^2/(x2-2*x1+x0);k=k+1; end 2 ,020102)(023==-++=x x x x x f
x3 k 牛顿法： function y=newton(x0) x1=x0-fc(x0)/df(x0); k=1; while (abs(x1-x0)>=1.0e-3) x0=x1; x1=x0-fc(x0)/df(x0);k=k+1; end x1 k function y=fc(x) y=x^3+2*x^2+10*x-20; function y=df(x) y=3*x^2+4*x+10; 第六周 1.解例6-4（p77)的方程组，分别采用消去法（矩阵分解）、Jacobi迭代法、Seidel迭代法、松弛法求解，并比较收敛速度。 消去法： x=a\d 或 [L,U]=lu(a); x=inv(U)inv(L)d Jacobi迭代法： function s=jacobi(a,d,x0) D=diag(diag(a)); U=-triu(a,1); L=-tril(a,-1); C=inv(D); B=C*(L+U); G=C*d; s=B*x0+G; n=1; while norm(s-x0)>=1.0e-8 x0=s; s=B*x0+G;
数学模型经典例题
一、把椅子往地面一放，通常只有三只脚着地，放不稳，然而只需稍挪动几次，就可以使四只脚同时着地放稳了，就四脚连线成长方形的情形建模并加以说明。（15分） 解：一、模型假设： 1. 椅子四只脚一样长，椅脚与地面的接触可以看作一个点，四脚连线呈长方形。 2. 地面高度是连续变化的，沿任何方向都不会出现间断，地面可以看成一张光滑曲面。 3. 地面是相对平坦的，使椅子在任何位置至少有三只脚同时着地。 （3分） 二、建立模型： 以初始位置的中位线为坐标轴建立直角坐标系，用θ表示椅子绕中心O 旋转的角度，椅子的位置可以用θ确定： ()f θ记为A 、B 两点与地面的距离之和 ()g θ记为C 、D 两点与地面的距离之和 由假设3可得，()f θ、()g θ中至少有一个为0。 由假设2知()f θ、()g θ是θ的连续函数。 （3分） 问题归结为： 已知()f θ和()g θ是θ的连续函数，对任意θ， ()()0f g θθ=，且设()()00,00g f =>。证明存在0θ， 使得()()000f g θθ== （3分） 三、模型求解： 令()()()h f θθθ=－g 若()()000f g =，结论成立 若()()000f g 、不同时为，不妨设()()00,00g f =>，椅子旋转()180π或后，AB 与CD 互换，即()()0,0g f ππ>=，则()(0)0,0h h π><。 （3分） 由f g 和的连续性知h 也是连续函数。根据连续函数的基本性质，必存在 ()000θθπ<<使000()0,()()h f g θθθ==即。 最后，因为00()()0f g θθ=，所以00()()0f g θθ==。 （3分） 图 5
数学建模习题及问题详解
第一部分课后习题 1.学校共1000名学生，235人住在A宿舍，333人住在B宿舍，432人住在C宿舍。学生 们要组织一个10人的委员会，试用下列办法分配各宿舍的委员数： （1）按比例分配取整数的名额后，剩下的名额按惯例分给小数部分较大者。 （2）2.1节中的Q值方法。 （3）d’Hondt方法：将A，B，C各宿舍的人数用正整数n=1，2，3，…相除，其商数 将所得商数从大到小取前10个（10为席位数），在数字下标以横线，表中A，B，C行有横线的数分别为2，3，5，这就是3个宿舍分配的席位。你能解释这种方法的道理吗。 如果委员会从10人增至15人，用以上3种方法再分配名额。将3种方法两次分配的结果列表比较。 （4）你能提出其他的方法吗。用你的方法分配上面的名额。 2.在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗。比如洁银牙膏50g 装的每支1.50元，120g装的3.00元，二者单位重量的价格比是1.2：1。试用比例方法构造模型解释这个现象。 （1）分析商品价格C与商品重量w的关系。价格由生产成本、包装成本和其他成本等决定，这些成本中有的与重量w成正比，有的与表面积成正比，还有与w无关的因素。 （2）给出单位重量价格c与w的关系，画出它的简图，说明w越大c越小，但是随着w 的增加c减少的程度变小。解释实际意义是什么。 3.一垂钓俱乐部鼓励垂钓者将调上的鱼放生，打算按照放生的鱼的重量给予奖励，俱乐部 只准备了一把软尺用于测量，请你设计按照测量的长度估计鱼的重量的方法。假定鱼池中只有一种鲈鱼，并且得到8条鱼的如下数据（胸围指鱼身的最大周长）： 先用机理分析建立模型，再用数据确定参数 4.用宽w的布条缠绕直径d的圆形管道，要求布条不重叠，问布条与管道轴线的夹角 应 多大（如图）。若知道管道长度，需用多长布条（可考虑两端的影响）。如果管道是其他形状呢。
数学建模题目及其答案(疾病诊断)
数学建模疾病的诊断 现要你给出疾病诊断的一种方法。 胃癌患者容易被误诊为萎缩性胃炎患者或非胃病者。从胃癌患者中抽 取5人（编号为1-5），从萎缩性胃炎患者中抽取5人（编号为6-10），以及非胃病者 中抽取5人（编号为11-15），每人化验4项生化指标：血清铜蓝蛋白（ X）、 1 蓝色反应（ X）、尿吲哚乙酸（3X）、中性硫化物（4X）、测得数据如表1 2 所示： 表1. 从人体中化验出的生化指标
* 根据数据，试给出鉴别胃病的方法。 论文题目：胃病的诊断 摘要 在临床医学中，诊断试验是一种诊断疾病的重要方法。好的诊断试验方法将对临床诊断的正确性和疾病的治疗效果起重要影响。因此，对于不同疾病不断发现新的诊断试验方法是医学进步的重要标志。传统的诊断试验方法有生化检测、DNA检测和影像检测等方法。而本文则通过利用多元统计分析中的判别分析及SPSS软件的辅助较好地解决了临床医学中胃病鉴别的问题。在临床医学上，既提高了临床诊断的正确性，又对疾病的治疗效果起了重要效果，同时也减轻了病人的负担。 判别分析是在分类确定的条件下，根据某一研究对象的各种特征值判别其类型归属问题的一种多变量统计分析方法。 其基本原理是按照一定的判别准则，建立一个或多个判别函数，用研究对象的大量资料确定判别函数中的待定系数，并计算判别指标。 , 首先，由判别分析定义可知，只有当多个总体的特征具有显著的差异时，进行判别分析才有意义，且总体间差异越大，才会使误判率越小。因此在进行判别分析时，有必要对总体多元变量的均值进行是否不等的显著性检验。 其次，利用判别分析中的费歇判别和贝叶斯判别进行判别函数的建立。 最后，利用所建立的判别函数进行回判并测得其误判率，以及对其修正。 本文利用SPSS软件实现了对总体间给类变量的均值是否不等的显著性检验并根据样本建立了相应的费歇判别函数和贝叶斯判别函数，最后进行了回判并测得了误判率，从而获得了在临床诊断中模型，给临床上的诊断试验提供了新方法和新建议。 关键词：判别分析；判别函数；Fisher判别；Bayes判别 一问题的提出 在传统的胃病诊断中，胃癌患者容易被误诊为萎缩性胃炎患者或非胃病患者，为了
经典的数学建模例子1
经典的数学建模例子 一、摘要 SARS SARS就是传染性非典型肺炎,全称严重急性呼吸综合症(Severe Acute Respiratory Syndromes),简称SARS,是一种因感染SARS相关冠状病毒而导致的以发热、干咳、胸闷为主要症状,严重者出现快速进展的呼吸系统衰竭,是一种新的呼吸道传染病,传染性极强、病情进展快速。 当一种传染病流行的时候,会给人们的工作学习带来很大的不变,能有效地进行隔离、预防,会大大减少人员的得病率,当一种传染病开始流行时,在一定的条件下其趋势就像真菌的繁殖曲线,如果能通过计算预测但大概推算出其发病率高峰时期,及时的隔离预防。那会给社会人力带来很大的方便,当年SARS的爆发给我们带来和大的不便和损失,因此本论文就以SARS为例,来研究传染病的传播规律、为预测和控制传染病蔓延创造条件和帮助。 1 二、正文 1、模型的背景问题描述 SARS(Severe Acute Respiratory Syndrome,严重急性呼吸道综合症, 俗称:非典型肺炎)是21世纪第一个在世界范围内传播的传染病。SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响,我们从中得到了许多重要的经验和教训,认识到定量地研究传染病的传播规律、为预测和控制传染病蔓延创造条件的重要性。 要求:(1)建立传染病传播的指数模型,评价其合理性和实用性。 (2)建立一个适合的模型,说明为什么优于问题1中的模型;特别要说明怎样才能 3 建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型,这样做的困难在哪里?对于卫生部门所采取的措施做出评论,如:提前或延后5天采取严格的隔离措施,对疫情传播所造成的影响做出估计。表中提供的数据供参考。 (3)说明建立传染病数学模型的重要性。 2、模型假设 (一)答;
数学建模考试题(开卷)及答案
2010年上学期2008级数学与应用数学，信息与计算科学专业 《数学建模》课程考试供选试题 第1题 4万亿投资与劳动力就业： 2008以来，世界性的金融危机席卷全球，给我国的经济发展带来很大的困难。沿海地区许多中小企业纷纷裁员，造成大量的人员失业。据有关资料估计，从2008年底，相继有2000万人被裁员，其中有1000万人是民工。部分民工返乡虽然能够从一定程度上缓解就业压力，但2009年的600多万毕业大学生给我国就业市场带来巨大压力。但可喜的是，我国有庞大的外汇储备，民间资本实力雄厚，居民储蓄充足。中国还是发展中国家，许多方面的建设还处于落后水平，建设投资的潜力巨大。为保持我国经济快速发展，特别是解决就业问题带来希望，实行政府投资理所当然。在2009年两代会上，我国正式通过了4万亿的投资计划，目的就是保GDP增长，保就业，促和谐。但是有几个问题一直困扰着我们，请你运用数学建模知识加以解决。问题如下： 1、GDP增长8%，到底能够安排多少人就业？如果要实现充分就业，2009年的GDP到底要增长多少？ 2、要实现GDP增长8%，4万亿的投资够不够？如果不够，还需要投资多少？ 3、不同的产业(或行业)吸纳的劳动力就业能力不同，因此投资的流向会有所不同。请你决策，要实现劳动力就业最大化，4万亿的投资应该如何分配到不同的产业(或行业)里？ 4、请你给出相关的政策与建议。 第2题 深洞的估算：假如你站在洞口且身上仅带着一只具有跑秒功能的计算器，你出于好奇心想用扔下一块石头听回声的方法来估计洞的深度，假定你捡到一块质量是1KG的石头，并准确的测定出听到回声的时间T=5S，就下面给定情况，分析这一问题，给出相应的数学模型，并估计洞深。 1、不计空气阻力； 2、受空气阻力，并假定空气阻力与石块下落速度成正比，比例系数k1=0.05； 3、受空气阻力，并假定空气阻力与石块下落速度的平方成正比，比例系数k2=0.0025； 4、在上述三种情况下，如果再考虑回声传回来所需要的时间。 第3题 优秀论文评选：在某数学建模比赛的评审过程中，组委会需要在一道题目的150 篇参赛论文中选择4 篇论文作为特等奖论文。评审小组由10 名评委组成，包括一名小组组长(出题人)，4 名专业评委(专门从事与题目相关问题研究的评委)，5 名普通评委(从事数学建模的教学和组织工作，参与过数学建模论文的评审)。组委会原先制定的评审步骤如下： step1:首先由普通评委阅读所有150 篇论文，筛选出20 篇作为候选论文。 Step2:然后由小组内的所有评委阅读这些候选论文，每人选择4 篇作为推荐的论文。 Step3:接着进入讨论阶段，在讨论阶段中每个评委对自己选择的 4 篇论文给出理由，大家进行讨论，每个评委对论文的认识都会受到其他评委观点的影响。 Step4:在充分讨论后，大家对这些推荐的论文进行投票，每个评委可以投出4票，获得至少6 票的论文可以直接入选，如果入选的论文不足，对剩余的论文(从20篇候选论文中除去已经入选的论文)重复step2至step4 步的评审工作。如果三轮讨论后入选的论文仍然不够，则由评选小组组长确定剩下名额的归属。 如果有超过4 篇的论文获得了至少6票，则由评选小组组长确定最终的名额归属。问题:
matlab数学建模实例
第四周3. 中的三个根。 ，在求8] [0,041.76938.7911.1-)(2 3=-+=x x x x f function y=mj()for x0=0:0.01:8 x1=x0^3-11.1*x0^2+38.79*x0-41.769;if (abs(x1)<1.0e-8)x0 end end 4.分别用简单迭代法、埃特金法、牛顿法求解方程，并比较收敛性与收敛速度（ε分别取10-3、10-5、10-8）。 简单迭代法： function y=jddd(x0) x1=(20+10*x0-2*x0^2-x0^3)/20;k=1; while (abs(x1-x0)>=1.0e-3) x0=x1; x1=(20+10*x0-2*x0^2-x0^3)/20;k=k+1;end x1k 埃特金法： function y=etj(x0) x1=(20-2*x0^2-x0^3)/10;x2=(20-2*x1^2-x1^3)/10; x3=x2-(x2-x1)^2/(x2-2*x1+x0);k=1; while (abs(x3-x0)>=1.0e-3) x0=x3; x1=(20-2*x0^2-x0^3)/10;x2=(20-2*x1^2-x1^3)/10; x3=x2-(x2-x1)^2/(x2-2*x1+x0);k=k+1;end 2 ,020102)(023==-++=x x x x x f
x3 k 牛顿法： function y=newton(x0) x1=x0-fc(x0)/df(x0); k=1; while(abs(x1-x0)>=1.0e-3) x0=x1; x1=x0-fc(x0)/df(x0);k=k+1; end x1 k function y=fc(x) y=x^3+2*x^2+10*x-20; function y=df(x) y=3*x^2+4*x+10; 第六周 1.解例6-4（p77)的方程组，分别采用消去法（矩阵分解）、Jacobi迭代法、Seidel迭代法、松弛法求解，并比较收敛速度。 消去法： x=a\d 或 [L,U]=lu(a); x=inv(U)inv(L)d Jacobi迭代法： function s=jacobi(a,d,x0) D=diag(diag(a)); U=-triu(a,1); L=-tril(a,-1); C=inv(D); B=C*(L+U); G=C*d; s=B*x0+G; n=1; while norm(s-x0)>=1.0e-8 x0=s; s=B*x0+G;
一些基本的数学建模示例
1.3 一些基本的数学建模示例 1.3.1椅子的摆放问题 1.3.2 双层玻璃的功效问题 1.3.3 搭积木问题 1.3.4 四足动物的身长和体重关系问题 1.3.5 圆杆堆垛问题 1.3.6 公平的席位分配问题 1.3.7 中国人重姓名问题 1.3.8实物交换问题 椅子能在不平的地面上放稳吗？下面用数学建模的方法解决此问题。 模型准备 仔细分析本问题的实质，发现本问题与椅子腿、地面及椅子腿和地面是否接触有关。如果把椅子腿看成平面上的点，并引入椅子腿和地面距离的函数关系就可以将问题1与平面几何和连续函数联系起来，从而可以用几何知识和连续函数知识来进行数学建模。为讨论问题方便，我们对问题进行简化，先做出如下3个假设： 模型假设 1、椅子的四条腿一样长，椅子脚与地面接触可以视为一个点，四脚连线是正方形（对椅子的假设） 2、地面高度是连续变化的，沿任何方向都不出现间断。（对地面的假设） 3、椅子放在地面上至少有三只脚同时着地，（对椅子和地面之间关系的假设） 根据上述假设做本问题的模型构成： 模型构成Array用变量表示椅子的位置，引入平面图形及坐 标系如图1-1。图中A、B、C、D为椅子的四只脚， 坐标系原点选为椅子中心，坐标轴选为椅子的四 只脚的对角线。于是由假设2，椅子的移动位置 可以由正方形沿坐标原点旋转的角度θ来唯一表 示，而且椅子脚与地面的垂直距离就成为θ的函 数。注意到正方形的中心对称性，可以用椅子的 相对两个脚与地面的距离之和来表示这对应两 个脚与地面的距离关系，这样，用一个函数就可 以描述椅子两个脚是否着地情况。本题引入两个 函数即可以描述椅子四图 1-1
数学建模案例分析线性代数建模案例例
线性代数建模案例汇编 目录
案例一. 交通网络流量分析问题 城市道路网中每条道路、每个交叉路口的车流量调查，是分析、评价及改善城市交通状况的基础。根据实际车流量信息可以设计流量控制方案，必要时设置单行线，以免大量车辆长时间拥堵。 【模型准备】 某城市单行线如下图所示, 其中的数字表示该路段每小时按箭头方向行驶的车流量(单位: 辆). 图3 某城市单行线车流量 (1) 建立确定每条道路流量的线性方程组. (2) 为了唯一确定未知流量, 还需要增添哪几条道路的流量统计? (3) 当x 4 = 350时, 确定x 1, x 2, x 3的值. (4) 若x 4 = 200, 则单行线应该如何改动才合理? 【模型假设】 (1) 每条道路都是单行线. (2) 每个交叉路口进入和离开的车辆数目相等. 【模型建立】 根据图3和上述假设, 在①, ②, ③, ④四个路口进出车辆数目分别满足 500 = x 1 + x 2 ① 400 + x 1 = x 4 + 300 ② x 2 + x 3 = 100 + 200 ③ x 4 = x 3 + 300 ④ 【模型求解】根据上述等式可得如下线性方程组 12142334500100300300x x x x x x x x +=??-=-??+=??-+=? 其增广矩阵 (A , b ) =1100500100110001103000011300?? ?-- ? ? ?-??????→初等行变换10011000101600001130000000--?? ? ?-- ? ?? ? 由此可得
142434 100600300x x x x x x -=-??+=??-=-? 即 14243 4100600300x x x x x x =-??=-+??=-?. 为了唯一确定未知流量, 只要增添x 4统计的值即可. 当x 4 = 350时, 确定x 1 = 250, x 2 = 250, x 3 = 50. 若x 4 = 200, 则x 1 = 100, x 2 = 400, x 3 = ?100 < 0. 这表明单行线“③?④”应该改为“③?④”才合理. 【模型分析】(1) 由(A , b )的行最简形可见, 上述方程组中的最后一个方程是多余的. 这意味着最后一个方程中的数据“300”可以不用统计. (2) 由142434100600300x x x x x x =-??=-+??=-?可得213141500200100x x x x x x =-+??=-??=+?, 123242500300600x x x x x x =-+??=-+??=-+?, 13234 3200300300x x x x x x =+??=-+??=+?, 这就是说x 1, x 2, x 3, x 4这四个未知量中, 任意一个未知量的值统计出来之后都可以确定出其他三个未知量的值. Matlab 实验题 某城市有下图所示的交通图, 每条道路都是单行线, 需要调查每条道路每小时的车流量. 图中的数字表示该条路段的车流数. 如果每个交叉路口进入和离开 图4 某城市单行线车流量 (1)建立确定每条道路流量的线性方程组. (2)分析哪些流量数据是多余的. (3)为了唯一确定未知流量, 需要增添哪几条道路的流量统计.
数学建模例题及解析
。 例1差分方程——资金的时间价值 问题1:抵押贷款买房——从一则广告谈起 每家人家都希望有一套(甚至一栋)属于自己的住房，但又没有足够的资金一次买下，这就产生了贷款买房的问题。先看一下下面的广告(这是1991年1月1日某大城市晚报上登的一则广告)，任何人看了这则广告都会产生许多疑问，且不谈广告中没有谈住房面积、设施等等，人们关心的是：如果一次付款买这栋房要多少钱呢?银行贷款的利息是多少呢?为什么每个月要付1200元呢?是怎样算出来的?因为人们都知道，若知道了房价(一次付款买房的价格)，如果自己只能支付一部分款，那就要把其余的款项通过借贷方式来解决，只要知道利息，就应该可以算出五年还清每月要付多少钱才能按时还清贷款了，从而也就可以对是否要去买该广告中所说的房子作出决策了。现在我们来进行数学建模。由于本问题比较简单无需太多的抽象和简化。 a.明确变量、参数，显然下面的量是要考虑的： 需要借多少钱，用记； 月利率(贷款通常按复利计)用R记； 每月还多少钱用x记； 借期记为N个月。 b．建立变量之间的明确的数学关系。若用记第k个月时尚欠的款数，则一个月后(加上利息后)欠款，不过我们又还了x元所以总的欠款为 k=0，1，2，3， 而一开始的借款为。所以我们的数学模型可表述如下 (1) c. (1)的求解。由
(2) 这就是之间的显式关系。 d．针对广告中的情形我们来看(1)和(2)中哪些量是已知的。N=5年＝60个月，已知；每月还款x＝1200元，已知A。即一次性付款购买价减去70000元后剩下的要另外去借的款，并没有告诉你，此外银行贷款利率R也没告诉你，这造成了我们决策的困难。然而，由(2)可知60个月后还清，即，从而得 (3) A和x之间的关系式，如果我们已经知道银行(3)表示N＝60，x＝1200给定时0 A。例如，若R ＝0．01，则由(3)可算得的贷款利息R，就可以算出0 53946元。如果该房地产公司说一次性付款的房价大于70000十53946＝123946元的话，你就应自己去银行借款。事实上，利用图形计算器或Mathematica这样的 数学软件可把(3)的图形画出来，从而可以进行估算决策。以下我们进一步考虑下面两个问题。 注1问题1标题中“抵押贷款”的意思无非是银行伯你借了钱不还，因而要你用某种不动产(包括房子的产权)作抵押，即万一你还不出钱了，就没收你的不动产。例题1某高校一对年青夫妇为买房要用银行贷款60000元，月利率0．01，贷款期25年＝300月，这对夫妇希望知道每月要还多少钱，25年就可还清。假设这对夫妇每月可有节余900元，是否可以去买房呢?