数学模型课程期末大作业题
要求:
1）选题方式：共53题，每个同学做一题，你要做的题目编号是你的学号mod52所得的值+1。(例如：你的学号为119084157，则你要做的题为mod(119084157,52)+1=50)。
2）该类题目基本为优划问题，要求提交一篇完整格式的建模论文，文字使用小四号宋体，公式用word的公式编辑器编写，正文中不得出现程序以及程序冗长的输出结果，程序以附录形式附在论文的后面，若为规划求解必须用lingo 集合形式编程，其它可用Matlab或Mathmatica编写。
3）论文以纸质文档提交，同时要交一份文章和程序电子文档，由班长统一收上来，我要验证程序。
1、生产安排问题
某厂拥有4台磨床，2台立式钻床，3台卧式钻床，一台镗床和一台刨床，用以生产7种产品，记作p1至p7。工厂收益规定作产品售价减去原材料费用之余。每种产品单件的收益及所需各机床的加工工时（以小时计）列于下表（表1）：
表
到6月底每种产品有存货50件。
工厂每周工作6天，每天2班，每班8小时。
不需要考虑排队等待加工的问题。
在工厂计划问题中，各台机床的停工维修不是规定了月份，而是选择最合
适的月份维修。除了磨床外，每月机床在这6个月中的一个月中必须停工维修；6个月中4台磨床只有2台需要维修。扩展工厂计划模型，以使可作上述灵活安排维修时间的决策。停工时间的这种灵活性价值若何？
注意，可假设每月仅有24个工作日。
5、生产计划
某厂有4台磨床，2台立钻，3台水平钻，1台镗床和1台刨床，用来生产7种产品，已知生产单位各种产品所需的有关设备台时以及它们的利润如表所示：
台镗床，4月—1台立钻，5月—1台磨床和1台立钻，6月—1台刨床和1台水平钻，被维修的设备在当月内不能安排生产。又知从1月到6月份市场对上述7种产品最大需求量如表所示：
量均不得超过100件。现在无库存，要求6月末各种产品各贮存50件。若该厂每月工作24天，每天两班，每班8小时，假定不考虑产品在各种设备上的加工顺序，要求：
（a）该厂如何安排计划，使总利润最大；
（b）在什么价格的条件下，该厂可考虑租用或购买有关的设备。
34、瓶颈机器上的任务排序
在工厂车间中，经常会出现整个车间的生产能力取决于一台机器的情况（例如，仅有一台的某型号机床，生产线上速度最慢的机器等）。这台机器就称为关键机器或瓶颈机器。此时很重要的一点就是尽可能地优化此机器将要处理的任务计划。
此问题的目的是为在单台机器上的任务调度提供一个简单的模型，此模型可以结合多种不同的目标函数进行使用。在这里我们将看到如何最小化总处理时间，平均处理时间，以及总超时时间。
在一台机器上将要处理一组任务。任务的执行不具有抢先性（即一旦一个任务开始执行，就不允许被打断）。对于每个任务i ，都给出了它的发布时间和持续时间。
对于最后一个优化目标（总超时时间），也需要使用截止时间（规定的最后完成时间）来对系统的超时长度进行度量，即度量任务完成时间超出规定时间的长度。下表中列出了我们的问题要使用的各种数据。
我们希望求出下面这些目标的最优值：计划总需时（makespan）的最小值，平均处理时间的最小值，或总超时时间的最小值。
35、油画制造
有一家油画公司有一些大型客户，它们一直有稳定的需求，此公司每周需要为这些客户制造 5 批油画，每批油画都完全相同。每批油画都在同一个制造过程中完成，所有批油画都要使用同一支调和画笔，在绘制两批油画之间必须清洗此画笔。第1到5 批油画的绘制时间分别为40，35，45，32，和50 分钟。清洗时间取决于所使用的颜色和颜料类型。例如，如果在使用水性颜料使用油性颜料，或者在使用深色后使用浅色，则需要较长的清洗时间。下表中给出了清洗时间数组CLEAN ，其中CLEANij表示在第i 批油画之后绘制第j 批油画所需的清洗时间。
的时间（绘制时间和清洗时间）。那么应采取什么顺序绘制这些批次的油画？所指定的顺序将每周重复执行，因此总清洗时间中也应计入一周的最后一批油画与下周的第一批油画之间所需的清洗时间。
36、生产线平衡
有一家电子工厂有一条生产线生产一种放大器，此生产线由四个工作台组成。生产放大器要经过12 道工序，这些工序之间存在先决关系约束。下表列出了每道工序需要花费的时间（分钟），并列出了所有的直接先决关系（表格中PCB 是印刷电路板的简称）。
制造管理人员希望在满足先决关系的条件下将这些工序分配到四个工作台
上，以使生产线得到平衡，从而使生产周期尽可能缩短，即缩短组装一台放大器所需的总时间。每道工序都需要分配到一个工作台上，并且在进行此工序时不许打断。每个工作台在一个时刻都只能进行一道工序。由于每个工作台上的每个工序都对每个放大器重复执行一次，因此我们称一台放大器组装所需的总时间为一个生产周期。当一台放大器完成组装之后，则工作台1 到3 上的放大器都将移动到下一个工作台上，并且在第一个工作台上开始组装新一台放大器。
37、自行车生产规划
有一家公司生产儿童自行车。在下表中给出了明年预期的销售量（以千辆为单位计）。此公司的生产能力为每个月30,000 辆自行车。通过工人加班，可以将产量提高50%，但是会将每辆自行车的生产成本从30 欧元提高到40 欧元。
底都需要支出 5 欧元的存储费用。我们假定此公司的库存能力是无限的（即虽然此公司的实际库存能力是有限的，但不会给我们这个例子带来限制）。现在是一月一日，在下面的十二个月里面每个月应生产和存储多少辆自行车才能够满足此销售预期，并最小化总成本？
38、考试日程安排
有一所工程大学每个学期三年级学生都需要根据其希望在第四年内学习的内容（可以从“生产规划”和“质量和安全管理”中选择）从11 门课程中选择8 门。在学生选择了下一个学年内的学习方向后，则在此学期内有些课程即变为必修课。这些必修课程是统计学（S），图模型与算法（GMA），生产管理（PM），离散系统与事件（DSE）。其他可选的课程为：数据分析（DA），数值分析（NA），数学规划（MP），C++，Java（J），逻辑规划（LP），以及软件工程（SE）
时间。安排两天用于考试，每天可用安排考试的时间段为：8:00-10:00，10:15 - 12:15， 14:00 -16:00，16:15- 18:15。在安排考试时间时需要考虑到有些学生同时选定了一些课程，因此这些课程的考试就不能安排在同时进行。表14.4.1 列出了这些有冲突的考试。
请帮助Edeetee 女士制定一份考试时间表，以使得每个学生每个时刻都只需参加一门考试。
39 玻璃杯生产规划
在法国北部有一家公司主要生产饮水用玻璃杯。当前此公司销售六种不同类型的玻璃杯（V1 到V6），这些杯子按照1000 个为一批进行生产，此公司希望为下面12个月的生产做出规划。可以生产小数数目批（不足1000 个）。下表列出了未来12 个星期预期每种类型的杯子的需求量。
子为单位）。已知每种杯子的每一批的生产成本和库存成本（单位为欧元），以及需要的工人和机器工作时间（小时）以及需要的储存空间（单位为货箱数）。
限制工人每周的总工时不能超过390 小时，机器每周的总工时不能超过850 小时。库存空间足够保存1000 个货箱。那么在每个时期内每种类型的杯子应各生产多少个才能够使生产和库存成本最小化？
40、物料需求规划
Minorette 公司生产两种大型儿童玩具车：蓝色集装箱卡车和红色油罐卡车。每种类型的玩具车都由13 个部件组装而成。图8.1 列出了这些部件的分解
图8.1：组件分解（Gozinto 图）
其中轮轴，底盘，蓝色或红色驾驶室可以由此公司自己组装，也可以转包给别的公司。下表列出了自己组装和转包这些组件的成本以及此公司的产能。在组装成本中未计入购买原料的费用。
量为0。那么Minorette 公司购买或转包这些组件的数量各为多少才能够在满足需求的同时又能够最小化生产成本？
41、电子元件生产规划
为增加利润，一家小公司希望提高其最畅销产品的产量。这家公司的主要业务之一是生产IC 卡和电子徽章。这家公司也生产这些卡和徽章的各个元件。因此对这些元件的生产进行良好的规划就成了此工厂取得成功的决定因素。在这个
例子中对这些元件的需求都来源于此公司内部，因此很容易对其进行预测。
在未来的六个月内，计划生产四种元件，型号分别为X43-M1，X43-M2，Y54-N1， Y54-N2。这些元件的产量受到产能变化的影响，并且每次产能改变后都需要重新进行控制和调整，因此会带来不可忽略的费用。因此公司希望最小化这些改变带来的费用，以及生产和库存的成本。
在下表中列出了每种产品每个时期内的需求量，生产和库存成本，初始库存量，以及最后希望保留的库存量。当产量发生变化时，需要对机器和控制系统进行重新调整。由此带来的费用与产量较前一月的改变量（提高或上升）成正比。产量每提高一个产品单位，则需要支出 1 欧元；产量每降低一个产品单位，只需要支出0.50 欧元。
为最小化由于产量改变引起的费用，以及生产和库存成本，应采取何种生产方案？
42、玻璃纤维生产规划
有一家公司生产玻璃纤维，产量以立方米为单位计算。这家公司希望对未来六个星期的生产进行规划。产能有一定上限，且在每个时期产能的上限都不同。规划所覆盖的整个期间的每周需求量都已知。不同时期的生产和存储的费用也不相同。下表中列出了这些数据。
43、货车载荷平衡
有三节铁路货车车厢，其最大允许载重均为100 公担（1 公担=100 千克），将使用这三节车厢运输16 个箱子。下表中列出了这些箱子的重量，单位为公担。应如何将箱子分配到各个货车上，才能使每节货车实际载重均不超过最大允许载重，且使装载量最大的车厢的装载量最小？
表格9.1：箱子属性
对问题进行求解：在所有箱子都装到车上之前，我们选择最重的尚未装车的箱子，并将它放到实际载重最轻的货车上。
44、储存罐注入
5 艘油船抵达了一个化工厂。这些油船中装有液体载荷，这些液体不允许混合： 1200 吨苯，700 立方米丁醇，1000 立方米丙醇，450 立方米苯乙烯，1200 立方米四氢呋喃（THF）。在此化工厂里有九个不同大小的存储罐。有一些存储罐已经装有一些液体。下表列出了每储存罐的属性（单位为吨）。应将这些船中装载的液体分别装载到哪个储存罐中才能使未使用的储存罐容量最大（问题1），或怎样才能使保留未用的储存罐数目最多？
45、汽车租赁
有一家小型汽车租赁公司，此公司有94 辆可供出租的汽车，分布于10 个代理点中。每个代理点的位置都将以地理坐标X 和Y 的形式给出，单位为千米。我们假定两个代理点之间的距离约为它们之间欧氏距离（即最短距离）的 1.3 倍。下表给出了各个代理点的位置坐标，以及第二天早晨汽车租赁的需求量和前一天晚上各个代理点拥有的汽车数。
之间调度分配汽车才能够满足各处的需求，并且使转运成本最低。
46、选择运输方式
在法国西南部有一家公司，这家公司需要将180 吨存放于仓库D1 到D4 中的化学产品运输到3 个回收中心C1，C2 和C3。仓库D1 到D4 分别储存有50，40，35，和65 吨化学产品，总计为190 吨。可以选用两种运输方式：公路运输和铁路运输。仓库D1 只能通过公路向回收中心C1 和C2 进行运输，运费分别为12 欧元/吨和14欧元/吨。仓库D2 只能向回收中心C2 运输，可以选择通过铁路或公路，运费分别为12 欧元/吨和14 欧元/吨。仓库D3 可以通过公路向回收中心C2 运输（9 欧元/吨），或通过铁路或公路向回收中心C3 运输，运费
分别为4 欧元/吨和5 欧元/吨。仓库D4可以通过铁路或公路向回收中心C2 运输，运费分别为11 欧元/吨和14 欧元/吨，或者通过铁路或公路向回收中心C3 运输，运费分别为10 欧元/吨和14 欧元/吨。
此公司与铁路公司签订的化学物品运输合同规定，每次运输量至少应为10 吨，最多为50 吨。除了标准的安全规章之外，对公路运输不存在其他特殊的限制。那么此公司应如何运输这180 吨化学物品才能够使总运费最低？
47、仓库位置选取
有一家大公司希望开设一些新的仓库，以向销售中心供货。每开设一个新仓库都有一些固定费用。货物将从仓库运输到附近的销售中心。每次运输的运费取决于运输的距离。这两种类型的费用非常不同：仓库开设费用属于投资支出，通常在若干年后将勾销，而运输费用属于运营成本。如何结合这两种费用不属于本书的讨论范围，我们假定这两种费用可比，为此可能需要以年为单位计算运营费用。
有12 个可以建造新仓库的位置，并且需要从这些仓库向12 个销售中心供货。
下表10.3 给出了每个仓库完全满足每个客户（销售中心）需求所需的总成本（千欧元，不是单位成本）。因此，例如从仓库1 向客户9（根据表10.5 可以看到此客户总需求量为30 吨）供货的单位成本为60000 欧元/30 吨，即2000 欧元/吨。如果无法进行送货，则对应的成本标记为无穷大∞。
数）和仓库的容量上限，这些信息都列于表10.4 中。
哪些位置开办仓库才能使总的建设成本以及运输成本最低，同时仍然能够满足所有客户需求？
48、燃油运输
有一个运输商需要将一些燃油从位于S处的炼油厂运输到他的一些客户那里。这些客户分别位于A1，A2，A3，A4，A5，和A6。下表列出了每个地方的需求量升数。
向所有客户运输的总里程数最少。
49、组合使用不同运输方式
有20 吨货物需要沿着一条路径运输到五个城市，可以选择三种不同的运输方式：铁路，公路，和航空运输。在三个位于路途中间的城市里可以更改运输方式，但是在相邻的两个城市之间只能采取一种运输方式。表10.9 列出了在每一对城市之间运输1 吨货物的成本。
下面的表（10.10）列出了在更换运输方式时每吨货物需要的额外支出。此支出与地点无关。表格10.10：更换运输工具费用
50、货车车队规划
有一个连锁商店从不同的汽车租赁商那里租赁货车，从而组成一支车队。此公司预测未来六个月内的货车的需求如下（表10.12）：
为满足需求，此连锁店可以选择三种类型的租赁合同，每个合同都将在每个月1日生效，这些合同为：三个月的租赁合同，每辆车租金总计1700 元；4 个月的租赁合同，每辆车租金总计2200 元；以及5 个月的租赁合同，每辆车租金总计2600 元。
此公司每个月应签订每种类型的合同各多少份才能够满足业务需要，并使总支出最小，且在六月结束时所有车辆租赁期都结束？
图10.4：第5 个月（五月）内仍然在租赁期的车辆
51机组编成
在第二次世界大展中，英国皇家空军（RAS）中有很多说不同语言的飞行员，飞行员所学习驾驶的飞机也不尽相同。皇家空军希望为每架飞机安排一对飞行员-副飞行员（一个机组），他们必须语言相通，并且都熟悉此机型。在我们这个例子中有8 名飞行员。下面的表中用0（最差）到20（最好）表示每名飞行员对各种语言（英语，法语，荷兰语，挪威语）的掌握程度，以及对不同类型的双座飞机（侦察机，运输机，轰炸机，战斗轰炸机，补给运输机）的驾驶经验。
表格11.3：飞行员各项评分
都超过10/20，这样才能组成一个机组。
问题1：是否有可能让所有飞行员都编入机组？
然后计算每个可行机组人员对其掌握得分均超过10/20 的每种类型飞机的掌握得分之和。这样我们可以定义每个机组的最高得分。例如，飞行员5 和6 对轰炸机的掌握分别为13 和10，对补给运输机的掌握分别为12 和18。因此它们组成的机组的最高得分为max(13+10, 12+18) = 30。
问题2：应如何编组才能让所有机组的得分之和最大？
52、航班着陆调度(具有时间窗口的调度问题)
在大型机场中，飞机的降落要受到很多安全约束条件的限制。本节中将研究如何对单条跑道上的飞机降落进行调度。已经有人研究了更一般的问题，但这些问题相当复杂（动态案例，例如航班晚点，同时有多条跑道，等等），因此我们只讨论一个简单的情形。
有十个航班需要降落。每个航班都有一个最早到达时间（飞机以最高速度到达降落区域的时间）和最晚到达时间（可能受其他因素如燃油量等的影响）。在这个时间窗口内，航空公司需要选择一个目标时间，并将它作为航班到达时间公布出去。如果比此目标时间迟到或早到，则可能会引起机场秩序混乱并带来额外的费用支出。为将这些费用计入考虑，并方便进行比对，每个航班都定义了早到每分钟的惩罚和晚到每分钟的惩罚。下表列出了每个航班的时间窗口（以从当天零时起分钟数计）和惩罚值。
一段安全时间。在表11.5 中第p 行第q 列即表示在航班p 和q 降落之间需要等待的最短时间（分钟），即便这两个航班实际上不是连续降落的。应采取何种降落调度方案才能够在使总惩罚最小，同时航班又都在指定的时间窗口中降落，并且满足两个航班降落之间的时间间隔？
HIMCM 2014美国中学生数学建模竞赛试题 Problem A: Unloading Commuter Trains Trains arrive often at a central Station, the nexus for many commuter trains from suburbs of larger cities on a “commuter” line. Most trains are long (perhaps 10 or more cars long). The distance a passenger has to walk to exit the train area is quite long. Each train car has only two exits, one near each end so that the cars can carry as many people as possible. Each train car has a center aisle and there are two seats on one side and three seats on the other for each row of seats.To exit a typical station of interest, passengers must exit the car, and then make their way to a stairway to get to the next level to exit the station. Usually these trains are crowded so there is a “fan” of passengers from the train trying to get up the stairway. The stairway could accommodate two columns of people exiting to the top of the stairs.Most commuter train platforms have two tracks adjacent to the platform. In the worst case, if two fully occupied trains arrived at the same time, it might take a long time for all the passengers to get up to the main level of the station.Build a mathematical model to estimate the amount of time for a passenger to reach the street level of the station to exit the complex. Assume there are n cars to a train, each car has length d. The length of the platform is p, and the number of stairs in each staircase is q. Use your model to specifically optimize (minimize) the time traveled to reach street level to exit a station for the following: 问题一:通勤列车的负载问题 在中央车站,经常有许多的联系从大城市到郊区的通勤列车“通勤”线到达。大多数火车很长(也许10个或更多的汽车长)。乘客走到出口的距离也很长，有整个火车区域。每个火车车厢只有两个出口,一个靠近终端, 因此可以携带尽可能多的人。每个火车车厢有一个中心过道和过道两边的座椅，一边每排有两个座椅，另一边每排有三个座椅。走出这样一个典型车站,乘客必须先出火车车厢,然后走入楼梯再到下一个级别的出站口。通常情况下这些列车都非常拥挤,有大量的火车上的乘客试图挤向楼梯，而楼梯可以容纳两列人退出。大多数通勤列车站台有两个相邻的轨道平台。在最坏的情况下,如果两个满载的列车同时到达,所有的乘客可能需要很长时间才能到达主站台。建立一个数学模型来估计旅客退出这种复杂的状况到达出站口路上的时间。假设一列火车有n个汽车那么长,每个汽车的长度为d。站台的长度是p，每个楼梯间的楼梯数量是q。使用您的模型具体来优化(减少)前往主站台的时间，有如下要求: Requirement 1. One fully occupied train's passengers to exit the train, and ascend the stairs to reach the street access level of the station. 要求1.一个满载乘客的火车，所有乘客都要出火车。所有乘客都要出楼梯抵达出主站台的路上。 Requirement 2. Two fully occupied trains' passengers (all passengers exit onto a common platform) to exit the trains, and ascend the stairs to reach the street access level
西南大学2014年春《数学建模》作业及答案(已整理) 第一次作业 1：[填空题] 名词解释: 1．原型 2．模型 3．数学模型 4．机理分析 5．测试分析 6．理想方法 7．计算机模拟 8．蛛网模型 9．群体决策 10．直觉 11．灵感 12．想象力 13．洞察力 14．类比法 15．思维模型 16．符号模型 17．直观模型 18．物理模型19．2倍周期收敛20．灵敏度分析21．TSP问题22．随机存储策略23．随机模型24．概率模型25．混合整数规划26．灰色预测 参考答案： 1．原型：原型指人们在现实世界里关心、研究或者从事生产、管理的实际对象。2．模型：指为某个特定目的将原形的某一部分信息简缩、提炼而构造的原型替代物。3．数学模型：是由数字、字母或其它数字符号组成的，描述现实对象数量规律的数学公式、图形或算法。4．机理分析：根据对客观事物特性的认识，找出反映内部机理的数量规律，建立的模型常有明显的物理意义或现实意义。5．测试分析：将研究对象看作一个"黑箱”系统，通过对系统输入、输出数据的测量和统计分析，按照一定的准则找出与数据拟合得最好的模型。6．理想方法：是从观察和经验中通过想象和逻辑思维，把对象简化、纯化，使其升华到理状态，以其更本质地揭示对象的固有规律。7．计算机模拟：根据实际系统或过程的特性，按照一定的数学规律用计算机程序语言模拟实际运行情况，并依据大量模拟结构对系统或过程进行定量分析。8．蛛网模型：用需求曲线和供应曲线分析市场经济稳定性的图示法在经济学中称为蛛网模型。9．群体决策：根据若干人对某些对象的决策结果，综合出这个群体的决策结果的过程称为群体决策。10．直觉：直觉是人们对新事物本质的极敏锐的领悟、理解或推断。11．灵感：灵感是指在人有意识或下意识思考过程中迸发出来的猜测、思路或判断。12．想象力：指人们在原有知识基础上，将新感知的形象与记忆中的形象相互比较、重新组合、加工、处理，创造出新形象，是一种形象思维活动。13．洞察力：指人们在充分占有资料的基础上，经过初步分析能迅速抓住主要矛盾，舍弃次要因素，简化问题的层次，对可以用那些方法解决面临的问题，以及不同方法的优劣作出判断。14．类比法：类比法注意到研究对象与以熟悉的另一对象具有某些共性，比较二者相似之处以获得对研究对象的新认识。15．思维模型：指人们对原形的反复认识，将获取的知识以经验的形式直接储存于人脑中，从而可以根据思维或直觉作出相应的决策。16．符号模型：是在一定约束条件或假设下借助于专门的符号、线条等，按一定形式组合起来描述原型。17．直观模型：指那些供展览用的实物模型以及玩具、照片等，通常是把原型的尺寸按比例缩小或放大，主要追求外观上的逼真。18．物理模型：主要指科技工作者为一定的目的根据相似原理构造的模型，它不仅可以显示原型的外形或某些特征，而且可以用来进行模拟实验，间接地研究原型的某些规律。19．2倍周期收敛：在离散模型中，如果一个数列存在两个收敛子列就称为2倍周期收敛。20．灵敏度分析：系数的每个变化都会改变线性规划问题，随之也会影响原来求得的最优解。为制定一个应付各种偶然情况的全能方法，必须研究以求得的最优解是怎样随输入系数的变化而变化的。这叫灵敏性分析。21．TSP问题：在加权图中寻求最佳推销员回路的问题可以转化为在一个完备加权图中寻求最佳哈密顿圈的问题，称为TSP问题。22．随机存储策略：商店在订购货物时采用的一种简单的策略，是制定一个下界s和一个上界S，当周末存货不小于s时就不定货；当存货少于s 时就订货，且定货量使得下周初的存量达到S，这种策略称为随机存储策略。23．随机模型：如果随机因素对研究对象的影响必须考虑，就应该建立随机性的数学模型，简称为随机模型。24．概
2011年数学建模集训小题目 1.求下列积分的数值解 ? +∞ +-?23 2 2 3x x x dx 2.已知)s i n ()()c o s (),(2h t h t h t e h t f h t ++++=+，dt h t f h g ?=10 ),()(，画出 ]10,10[-∈h 时，)(h g 的图形。 3.画出16)5(2 2=-+y x 绕x 轴一周所围成的图形，并求所产生的旋转体的体积。 4.画出下列曲面的图形 （1）旋转单叶双曲面 14 92 22=-+z y x ； （2）马鞍面xy z =； 5.画出隐函数1cos sin =+y x 的图形。 6.（1）求函数x x y -+=12 ln 的三阶导数； 法一：syms x y dy; >> y=log((x+2)/(1-x)); >> dy=diff(y,3) dy = (6/(1-x)^3+6*(x+2)/(1-x)^4)/(x+2)*(1-x)-2*(2/(1-x)^2+2*(x+2)/(1-x)^3)/(x+2)^2*(1-x)-2*(2/(1-x)^2+2*(x+2)/(1-x)^3)/(x+2)+2*(1/(1-x)+(x+2)/(1-x)^2)/(x+2)^3*(1-x)+2*(1/(1-x)+(x+2)/(1-x)^2)/(x+2)^2 （2）求向量]425.00[=a 的一阶向前差分。 7.求解非线性方程组 （1）?????=-+=-+060622x y y x （2）???=+=++5 ln 10tan 10cos sin y x y e y x 8.求函数186)(2 3-++=x x x x f 的极值点，并画出函数的图形。 9.某单位需要加工制作100套钢架，每套用长为2.9m ，2.1m 和1m 的圆钢各一根。已知原料长6.9m ，问应如何下料，使用的原材料最省。 10. 某部门在今后五年内考虑给下列项目投资，已知： 项目A ，从第一年到第四年每年年初需要投资，并于次年末回收本利115％； 项目B ，从第三年初需要投资，到第五年末能回收本利125％，但规定最大投资额不超过4万元；
2003高教社杯全国大学生数学建模竞赛题目 （请先阅读“对论文格式的统一要求”） A题 SARS的传播 SARS（Severe Acute Respiratory Syndrome，严重急性呼吸道综合症, 俗称：非典型肺炎）是21世纪第一个在世界范围内传播的传染病。SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响，我们从中得到了许多重要的经验和教训，认识到定量地研究传染病的传播规律、为预测和控制传染病蔓延创造条件的重要性。请你们对SARS 的传播建立数学模型，具体要求如下： （1）对附件1所提供的一个早期的模型，评价其合理性和实用性。 （2）建立你们自己的模型，说明为什么优于附件1中的模型；特别要说明怎样才能建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型，这样做的困难在哪里？对于卫生部门所采取的措施做出评论，如：提前或延后5天采取严格的隔离措施，对疫情传播所造成的影响做出估计。附件2提供的数据供参考。
（3）收集SARS对经济某个方面影响的数据，建立相应的数学模型并进行预测。附件3提供的数据供参考。 （4）给当地报刊写一篇通俗短文，说明建立传染病数学模型的重要性。 附件1： SARS疫情分析及对北京疫情走势的预测 2003年5月8日 在病例数比较多的地区，用数理模型作分析有一定意义。前几天，XXX老师用解析公式分析了北京SARS疫情前期的走势。在此基础上，我们加入了每个病人可以传染他人的期限（由于被严格隔离、治愈、死亡等），并考虑在不同阶段社会条件下传染概率的变化，然后先分析香港和广东的情况以获得比较合理的参数，最后初步预测北京的疫情走势。希望这种分析能对认识疫情，安排后续的工作生活有帮助。 1 模型与参数 假定初始时刻的病例数为N0，平均每病人每天可传染K个人（K
练习1基础练习 一、矩阵及数组操作: 1．利用基本矩阵产生3×3和１5×８的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵([－１，1］之间）、正态分布矩阵（均值为1,方差为4）。 A=eye(3) Ｂ＝ｅyｅ(15,８)C=ones(3)D=oｎes(15,8）E=ｚeros （3) Ｆ=zeroｓ(1５,８) Ｇ＝(-1+(１-(-1))*ｒand(３)) H=1+ｓqrt（4)*ｒaｎdn(5) ２．利用fｉx及ｒanｄ函数生成[0,10]上的均匀分布的10×10的整数随机矩阵a,然后统计a中大于等于5的元素个数 ａ=fix（０＋(１０-0）*rand(1０)）; K=finｄ(a>＝5); Num=lenｇｔｈ(K)或者ｎｕm=suｍ（ｓum(a>=５)) num = 53 3.在给定的矩阵中删除含有整行内容全为０的行,删除整列内容全为0的列。 如已给定矩阵A在给定的矩阵中删除含有整行内容全为０的行 在命令窗口中输入Ａ(find(ｓum(ａｂｓ(A'))==０）,:)=[]; 删除整列内容全为0的列。A(：，finｄ(suｍ(ａbｓ(Ａ'))＝=0))=［];
二、绘图： 4．在同一图形窗口画出下列两条曲线图像: ｙ1=２x+5； ｙ2=x ＾2-3x+1， 并且用ｌegen ｄ标注 ｘ=０:０.01：10; y1=2*x+5; y ２=x.^2-３*x+１； plot(ｘ,y １，x,ｙ2,'ｒ') ｌeg ｅn ｄ('y １', 'y2') 12345678910 -1001020304050607080 5.画出下列函数的曲面及等高线： ｚ=x^２+y^2+ｓｉn(xy). 在命令窗口输入:
大致符合f=0.0069*x.^5-0.2669*x.^4+3.7108*x.^3-20.8889*x.^2+36.2143*x+60.0136;的分布，呈现夏季短日照冬季长日照，春秋居中的规律。1-6月份递减，
X=[0 400 800 1200 1600 2000 2400 2800]; Y=[0 400 800 1200 1600 2000 2400]; [x,y]=meshgrid(X,Y); z=[1180 1320 1450 1420 1400 1300 700 900;... 1230 1390 1500 1500 1400 900 1100 1060;... 1270 1500 1200 1100 1350 1450 1200 1150;... 1370 1500 1200 1100 1550 1600 1550 1380;... 1460 1500 1550 1600 1550 1600 1600 1600;... 1450 1480 1500 1550 1510 1430 1300 1200;... 1430 1450 1470 1320 1280 1200 1080 940]; figure(1); meshz(x,y,z) title('源数据点') xi=[0:40:2800]; yi=[0:40:2400]; zc=interp2(x,y,z,xi,yi','spline'); figure(2); surfc(xi,yi,zc); title('地貌图'); figure(3); contour(xi,yi,zc); 【2】结果：
习 题 1 1. 请编写绘制以下图形的MA TLAB 命令，并展示绘得的图形. (1) 221x y +=、224x y +=分别是椭圆2241x y +=的内切圆和外切圆. (2) 指数函数x y e =和对数函数ln y x =的图像关于直线y=x 对称. (3) 黎曼函数 1, (0)(0,1) 0 , (0,1), 0,1 q x p q q x y x x x =>∈?=? ∈=?当为既约分数且当为无理数且或者 的图像（要求分母q 的最大值由键盘输入）. 3. 两个人玩双骰子游戏，一个人掷骰子，另一个人打赌掷骰子者不能掷出所需点数，输赢的规则如下：如果第一次掷出3或11点，打赌者赢；如果第一次掷出2、7或12点，打赌者输；如果第一次掷出4、5、6、8、9或10点，记住这个点数，继续掷骰子，如果不能在掷出7点之前再次掷出该点数，则打赌者赢. 请模拟双骰子游戏，要求写出算法和程序，估计打赌者赢的概率. 你能从理论上计算出打赌者赢的精确概率吗？请问随着试验次数的增加，这些概率收敛吗？
4. 根据表1.14的数据，完成下列数据拟合问题： (1) 如果用指数增长模型0()0()e r t t x t x -=模拟美国人口从1790年至2000年的变化过程，请用MATLAB 统计工具箱的函数nlinfit 计算指数增长模型的以下三个数据拟合问题： (i) 取定0x =3.9，0t =1790，拟合待定参数r ； (ii) 取定0t =1790，拟合待定参数0x 和r ； (iii) 拟合待定参数0t 、0x 和r . 要求写出程序，给出拟合参数和误差平方和的计算结果，并展示误差平方和最小的拟合效果图. (2) 通过变量替换，可以将属于非线性模型的指数增长模型转化成线性模型，并用MA TLAB 函数polyfit 进行计算，请说明转化成线性模型的详细过程，然后写出程序，给出拟合参数和误差平方和的计算结果，并展示拟合效果图. (3) 请分析指数增长模型非线性拟合和线性化拟合的结果有何区别？原因是什么？ (4) 如果用阻滞增长模型00 () 00()()e r t t Nx x t x N x --= +-模拟美国人口从1790年至2000年的变化过程，请用MA TLAB 统计工具箱的函数nlinfit 计算阻滞增长模型的以下三个数据拟合问题： (i) 取定0x =3.9，0t =1790，拟合待定参数r 和N ； (ii) 取定0t =1790，拟合待定参数0x 、r 和N ； (iii) 拟合待定参数0t 、0x 、r 和N . 要求写出程序，给出拟合参数和误差平方和的计算结果，并展示误差平方和最小的拟合效果图. 年份 1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890
关于中国人口增长趋势的研究 【摘要】 本文从中国的实际情况和人口增长的特点出发，针对中国未来人口的老龄化、出生人口性别比以及乡村人口城镇化等，提出了Logistic、灰色预测、动态模拟等方法进行建模预测。 首先，本文建立了Logistic阻滞增长模型，在最简单的假设下，依照中国人口的历史数据，运用线形最小二乘法对其进行拟合，对2007至2020年的人口数目进行了预测，得出在2015年时，中国人口有13.59亿。在此模型中，由于并没有考虑人口的年龄、出生人数男女比例等因素，只是粗略的进行了预测，所以只对中短期人口做了预测，理论上很好，实用性不强，有一定的局限性。 然后，为了减少人口的出生和死亡这些随机事件对预测的影响，本文建立了GM(1,1) 灰色预测模型，对2007至2050年的人口数目进行了预测，同时还用1990至2005年的人口数据对模型进行了误差检验，结果表明，此模型的精度较高，适合中长期的预测，得出2030年时，中国人口有14.135亿。与阻滞增长模型相同，本模型也没有考虑年龄一类的因素，只是做出了人口总数的预测，没有进一步深入。 为了对人口结构、男女比例、人口老龄化等作深入研究，本文利用动态模拟的方法建立模型三，并对数据作了如下处理：取平均消除异常值、对死亡率拟合、求出2001年市镇乡男女各年龄人口数目、城镇化水平拟合。在此基础上，预测出人口的峰值，适婚年龄的男女数量的差值，人口老龄化程度，城镇化水平，人口抚养比以及我国“人口红利”时期。在模型求解的过程中，还对政府部门提出了一些有针对性的建议。此模型可以对未来人口做出细致的预测，但是需要处理的数据量较大，并且对初始数据的准确性要求较高。接着，我们对对模型三进行了改进，考虑人为因素的作用，加入控制因子，使得所预测的结果更具有实际意义。 在灵敏度分析中，首先针对死亡率发展因子θ进行了灵敏度分析，发现人口数量对于θ的灵敏度并不高，然后对男女出生比例进行灵敏度分析得出其灵敏度系数为0.8850，最后对妇女生育率进行了灵敏度分析，发现在生育率在由低到高的变化过程中，其灵敏度在不断增大。 最后，本文对模型进行了评价，特别指出了各个模型的优缺点，同时也对模型进行了合理性分析，针对我国的人口情况给政府提出了建议。 关键字：Logistic模型灰色预测动态模拟 Compertz函数
2011年全国大学生数学建模竞赛测试试题(A) 时量：180分钟满分：150分 院系：专业：学号：姓名： 一、选择题（2分/题×10题=20分） 1、Matlab程序设计中清除当前工作区的变量x,y的命令是( c ) A.clc x,y B.clear(x y) C.clear x y D.remove(x,y) 2、关于Matlab程序设计当中变量名和函数名的描述，下述说法正确的是( B ) A.都不区分大小写 B.都区分大小写 C.变量名区分,函数名不区分 D. 变量名区分,函数名不区分 3、MA TLAB软件中，把二维矩阵按一维方式寻址时的寻址访问是按（B）优先的。 A.行 B.列 C.对角线 D.左上角 4、关于矩阵上下拼接和左右拼接的方式中，下列描述是正确的是（ D ） A．上下拼接的命令为C=[A, B]，要求矩阵A, B的列数相同； B．左右拼接的命令为C=[A; B]，要求矩阵A, B的行数相同； C．上下拼接的命令为C=[A; B]，要求矩阵A, B的行数相同； D．左右拼接的命令为C=[A, B]，要求矩阵A, B的行数相同。 5、Matlab命令a=[65 72 85 93 87 79 62 73 66 75 70];find(a>=70 & a<80)得到的结果为（C ） A.[72 79 73 75] B.[72 79 73 75 70] C.[2 6 8 10 11] D.[0 1 0 0 0 1 0 1 0 1 1] 6、矩阵(或向量)的范数是用来衡量矩阵(或向量)的（A）的一个量 A.维数大小 B.元素的值的绝对值大小 C.元素的值的整体差异程度 D.所有元素的和 7、计算非齐次线性方程组AX=b的解可转化为计算矩阵X=A-1b，可以用Matlab的命令（A）实现 A.左除命令x=A\b B.左除命令x=A/b C.右除命令x=A\b D.右除命令x=A/b 8、关于Matlab的矩阵命令与数组命令，下列说法正确的是（b） A.矩阵乘A*B是指对应位置元素相乘 B.矩阵乘A.*B是指对应位置元素相乘 C.数组乘A.*B是指对应位置元素相乘 D.数组乘A*B是指对应位置元素相乘 9、生成5行4列，并在区间[1:10]内服从均分布的随机矩阵的命令是（d） A.rand(5,4)*10 B.rand(5,4,1,10) C.rand(5,4)+10 D.rand(5,4)*9+1 10、关于Matlab的M文件的描述中，以下错误的是（ d ） A、Matlab的M 文件有脚本M文件和函数M文件两种； B、Matlab的函数M文件中要求首行必须以function顶格开头；
2014年美国数学建模大赛（MCM）试题译文 王景璟大连理工大学 问题A：超车之外靠右行原则 在一些开车必须靠右行驶的国家（比如：美国，中国，以及其他除了英国，澳大利亚，和一些前英国殖民地的国家），行驶在多车道高速路必须遵循一个规则，那就是要求驾驶员在超车之外的情况下，必须在最靠右的车道行驶，超车时，他们向左变道，超车，然后再回到之前行驶的车道。 构建一个数学模型来分析该规则在车流量很少和很大的时候的执行情况。你最好能考察车流量与安全的之间的相互关系，过低或是过量的速度限制的作用（速度设置过低或是过高），以及/或者其他在该问题陈述中没有明确提到的因素。该原则是否能有效促进更好的车流量？如果无效，请建议和分析其他更有助于提高车流量、安全、以及其他你认为重要的因素的其他方案（可以完全不包括该原则）。 在开车靠左行的国家，讨论一下你的方案在经过对方向的简单修改之后或是添加额外的要求之后是否也适用。 最后，以上原则取决于人们遵循交通规则的判断力。如果道路上的车流完全在智能系统（要么是道路体系的一部分，要么是包含在使用道路的所有车辆的设计之中）的控制之下，该改变在多大程度上会影响你先前分析的结果？ 问题B: 大学教练联盟 《体育画报》，一本体育爱好者的杂志，正在寻找上世纪“最好的大学教练”，包括男性和女性。建立一个数学模型以从诸如大学曲棍球，曲棍球，橄榄球，棒球，垒球，篮球，或足球等运动的男性或女性教练中选出最好的一个教练或几个教练（过去的或现在的）。分析中使用的时间分界线是否有影响？即在1913执教和在2013年执教有不同吗？清晰地表达你们模型中的评判标准。讨论你们的模型如何能广泛地应用于两种性别及所有可能的体育运动。分别选出你模型中3种不同运动的前5位教练。 除了MCM格式及要求，准备一篇1-2页的文章给《体育画报》以解释你们的结论并包括一份能让体育迷们看懂的对你们数学模型的非技术性解释。 问题C：使用网络模型测量影响力
网络学院数学建模作业题
数学建模作业题 注意事项： 作业共十题，每题十分，全部是比较简单的建模计算题，题目既是课本上的习题，在课本304~315有参考解答，又是在线题库的题目，在题库里有更详细的解答。学员应该先自己动脑筋解决，然后才参考一下课本及题库的解答。 评分高低主要是看完成作业的态度、独立程度和表达清晰程度。 上传的作业必须是包括全部作业的单独一份word文档，必须自己录入，不允许扫描，不允许直接插入题库答案中的图片。严重违反者，不及格。 请于有效期结束前两周提交上传作业，教师尽快批改，请学员有效期结束前一周查看成绩，不及格的学员可以在课程答疑栏目提出或者课程论坛提出重交申请，教师删除原作业后，这些学员可以在有效期结束前之前重交作业。每人只有一次重交机会。 作业题与考试相关（当然不会一模一样），认真完成作业的学员，必将在考试取得好成绩。 一、教材76页第1章习题1第7题（来自高中数学课本的数学探究问题，满分10分） 表1.17是某地一年中10天的白昼时间（单位：小时），请选择合适的函数模型，并进行数据拟合. 日期1月1日2月28日3月21日4月27日5月6日
白昼时间 5.59 10.23 12.38 16.39 17.26 日期 6月21日 8月14日 9月23日 10月25日 11月21日 白昼时间 19.40 16.34 12.01 8.48 6.13 解：根据地理常识，某地的白昼时间是以一年为周期而变化的，以日期在一年中序号为自变量x ，以白昼时间为因变量y ，则根据表1.17的数据可知在一年（一个周期）内，随着x 的增加，y 大约在6月21日（夏至）达到最大值，在12月21日（冬至）达到最小值，在3月21日（春分） 或9月21日（秋分）达到中间值。选择函数y=(b x A ++)3652sin(?π)作为函数值。根据表1.17的数据，推测A ，b 和?的值， 作非线性拟合得385.123712.13652sin(9022.6+-=x y π，预测该地12月21日的白昼时间为5.49小时。 二、教材100页第2章习题2第1题（满分10分） 继续考虑第2.2节“汽车刹车距离”案例，请问“两秒准则”和“一车长度准则”一样吗？“两秒准则”是否足够安全？对于安全车距，你有没有更好的建议？ 解：“两秒准则”表明前后车距D 与车速v 成正比例关系v K D 2 =，其中s K 22 =，对于小型汽车，“一车长度准则”与“两秒准则”不一致。由)]([1 2 2 K K v K v D d --=-可以计算得到当D d h km K K K v <=-<时有/428.542 12 ，“两秒准则”足够安全，或者把刹车距离实测数据和“两秒准则”都画在同一幅图中，根据图形指出“两秒准则”足够安全的车速范围。用最大刹车距离除以车速，得到最大刹车距离所需的尾随时间，并以尾随时间为依据，提出更安全的准则，如“3秒准则”、“4
对作业题目的说明 1. 本次数学建模周一共提供十五道题目供大家选择。每支队伍（2-3人/队）必须从以下题目中任意选取一题（只须选择一道），并完成一篇论文，对论文的具体要求参阅《论文格式规范》。 2. 题目标注为“A ”的为有一定难度的题目，指导老师会根据题目的难度对论文最后的评分进行调整。 （一）乒乓球赛问题 (A) A 、 B 两乒乓球队进行一场五局三胜制的乒乓球赛，两队各派3名选手上场，并各有3种选手的出场顺序（分别记为123,,ααα 和123,,βββ）。根据过去的比赛记录，可以预测出如果A 队以i α次序出场而B 队以j β次序出场，则打满5局A 队可胜 ij a 局。由此得矩阵()ij R a =如下： 12 3 1232 140345 3 1R βββααα?? = ? ? ??? （1） 根据矩阵R 能看出哪一队的实力较强吗？ （2） 如果两队都采取稳妥的方案，比赛会出现什么结果？ （3） 如果你是A 队的教练，你会采取何种出场顺序？ （4） 比赛为五战三胜制，但矩阵R 中的元素却是在打满五局的情况下得到 的，这样的数据处理和预测方式有何优缺点？ （二）野兔生长问题 在某地区野兔的数量在连续十年的统计数量(单位十万)如下： 分析该数据，得出野兔的生长规律。 并指出在哪些年内野兔的增长有异常现象，
预测T=10 时野兔的数量。 （三）停车场的设计问题 在New England的一个镇上，有一位于街角处面积100 200平方英尺的停车场，场主请你代为设计停车车位的安排方式，即设计在场地上划线的方案。 容易理解，如果将汽车按照与停车线构成直角的方向，一辆紧挨一辆地排列成行，则可以在停车场内塞进最大数量的汽车，但是对于那些缺乏经验的司机来说，按照这种方式停靠车辆是有困难的，它可能造成昂贵的保险费用支出。为了减少因停车造成意外损失的可能性，场主可能不得不雇佣一些技术熟练的司机专门停车；另一方面，如果从通道进入停车位有一个足够大的转弯半径，那么，看来大多数的司机都可以毫无困难地一次停车到位。当然通道越宽，场内所容纳的车辆数目也越少，这将使得场主减少收入。 请你通过建模的计算结果，来给出一个合理的设计方案。 （四）奖学金的评定(A) 背景 A Better Class (ABC)学院的一些院级管理人员被学生成绩的评定问题所困 ）,这使得扰。平均来说，ABC的教员们一向打分较松（现在所给的平均分是A — 无法对好的和中等的学生加以区分。然而,某项十分丰厚的奖学金仅限于资助占总数10%的最优秀学生,因此,需要对学生排定名次。 教务长的想法是在每一课程中将每个学生与其他学生加以比较,运用由此得到的信息构造一个排名顺序。例如,某个学生在一门课程中成绩为A,而在同一课程中所有学生都得A,那么就此课而言这个学生仅仅属于“中等”。反之，如果一个学生得到了课程中唯一的A，那么，他显然处在“中等至上”水平。综合从几门不同课程所得到的信息，使得可以把所有学院的学生按照以10%划分等级顺序（最优秀的10%，其次的10%，等等）排序。 问题 , B+ ,…）这样的方式给出的,教务（1）假设学生成绩是按照（A+，A, A — 长的想法能否实现?
2014年下学期数学实验与数学建模作业习题8 1．轮船的甲板成近似半椭圆面形为了得到甲板的面积。首先测量得到横向最大相间8.534米；然后等间距地测得纵向高度，自左向右分别为：0.914, 5.060, 7.772, 8.717, 9.083, 9.144, 9.083, 8.992, 8.687, 7.376, 2.073，计算甲板的面积。 【1】命令： x=0:0.711:8.534; y2=[0,0.914^2,5.060^2,7.772^2,8.717^2,9.083^2,9.144^2,9.083^2,8.992^2, 8.687^2,7.376^2,2.073^2,0]; %plot(x,y2,'*'); a=polyfit(x,y2,2) 【2】结果： a = -5.2832 46.5248 -16.7465 得y^2=-5.2832*x^2+46.5248*x-16.7465,即y^2/85.68+(x-4.4031)^2/16.2175=1 故面积=0.5*a*b*pi=58.56. 2．物体受水平方向外力作用，在水平直线上运动。测得位移与受力如表8.1 表8.1 X 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 F 20 21 21 20 19 18.5 18.0 13.5 9 4.5 0 求(a) 物体从位移为0到0.4所做的功； (b) 位移为0.4时的速度是多少？ 【1】命令： x=0:0.1:1.0; f=[20,21,21,20,19,18.5,18.0,13.5,9,4.5,0]; plot(x,f,'*');hold on; a=polyfit(x,f,2) f2=-34.4988*x.*x+14.8625*x+19.5979; plot(x,f2); syms t y=-34.4988*t.*t+14.8625*t+19.5979; w=vpa(int(y,t,0,0.4),8) V=diff(y);t=2;v=eval(V)
数学模型课程期末大作业题 要求: 1）选题方式：共53题，每个同学做一题，你要做的题目编号是你的学号mod52所得的值+1。(例如：你的学号为119084157，则你要做的题为mod(119084157,52)+1=50)。 2）该类题目基本为优划问题，要求提交一篇完整格式的建模论文，文字使用小四号宋体，公式用word的公式编辑器编写，正文中不得出现程序以及程序冗长的输出结果，程序以附录形式附在论文的后面，若为规划求解必须用lingo 集合形式编程，其它可用Matlab或Mathmatica编写。 3）论文以纸质文档提交，同时要交一份文章和程序电子文档，由班长统一收上来，我要验证程序。 1、生产安排问题 某厂拥有4台磨床，2台立式钻床，3台卧式钻床，一台镗床和一台刨床，用以生产7种产品，记作p1至p7。工厂收益规定作产品售价减去原材料费用之余。每种产品单件的收益及所需各机床的加工工时（以小时计）列于下表（表1）： 表 到6月底每种产品有存货50件。 工厂每周工作6天，每天2班，每班8小时。 不需要考虑排队等待加工的问题。 在工厂计划问题中，各台机床的停工维修不是规定了月份，而是选择最合
适的月份维修。除了磨床外，每月机床在这6个月中的一个月中必须停工维修；6个月中4台磨床只有2台需要维修。扩展工厂计划模型，以使可作上述灵活安排维修时间的决策。停工时间的这种灵活性价值若何？ 注意，可假设每月仅有24个工作日。 5、生产计划 某厂有4台磨床，2台立钻，3台水平钻，1台镗床和1台刨床，用来生产7种产品，已知生产单位各种产品所需的有关设备台时以及它们的利润如表所示： 台镗床，4月—1台立钻，5月—1台磨床和1台立钻，6月—1台刨床和1台水平钻，被维修的设备在当月内不能安排生产。又知从1月到6月份市场对上述7种产品最大需求量如表所示： 量均不得超过100件。现在无库存，要求6月末各种产品各贮存50件。若该厂每月工作24天，每天两班，每班8小时，假定不考虑产品在各种设备上的加工顺序，要求： （a）该厂如何安排计划，使总利润最大； （b）在什么价格的条件下，该厂可考虑租用或购买有关的设备。 34、瓶颈机器上的任务排序 在工厂车间中，经常会出现整个车间的生产能力取决于一台机器的情况（例如，仅有一台的某型号机床，生产线上速度最慢的机器等）。这台机器就称为关键机器或瓶颈机器。此时很重要的一点就是尽可能地优化此机器将要处理的任务计划。
太阳影子定位问题 摘要 目前，如何确定视频的拍摄地点和拍摄日期是计算机视觉的热点研究问题，是视频数据分析的重要方面，有重要的研究意义。本文通过建立数学模型，给出了通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的方法。 对于问题一，建立空间三维直角坐标系和球面坐标系对直杆投影和地球进行数学抽象，引入地方时、北京时间、太阳赤纬、杆长、太阳高度角等五个参数，建立了太阳光下物体影子的长度变化综合模型。求解过程中，利用问题所给的数据，得到太阳赤纬等变量，将太阳赤纬等参量代入模型，求得了北京地区的9:00至15:00的影子长度变化曲线，当12：09时，影子长度最短；并分析出影长随这些参数的变化规律，利用控制变量法思想，总结了五个参数与影子长度的关系。最后进行模型检验，将该模型运用于东京、西藏两地，得到了这两座城市的影长变化规律曲线，发现变化规律符合实际两地实际情况。 对于问题二，为了消除不同直角坐标系带来的影响，将实际坐标转换为二次曲线的极坐标，建立了极坐标下基于多层优化搜索算法的空间匹配优化模型。求解时，先将未知点的直角坐标系的点转换为极坐标，然后设计了多层优化搜索算法，通过多次不同精度的搜索，最后得出实际观测点的经纬度为东经E115?北纬N25?。同时对模型进行验证，实地测量了现居住地的某个时间段的值，通过模型二来求解出现居住地的经纬度，分析了误差产生的原因：大气层的折射和拟合误差。 对于问题三，将极坐标转换后的基本模型转换为优化模型，建立了基于遗传算法的时空匹配优化模型。将目标函数作为个体的适应度函数，将经度纬度及日期作为待求解变量，用遗传算法进行求解，得到可能的经度纬度及其日期：北纬20度，东经114度，5月21日；北纬20度，东经114度，7月24日；东经94.5度，北纬33.8度，6月19日。最后，将遗传算法与多层优化搜索算法进行对比分析，得出遗传算法的求解效率和求解精度均优于多层次搜索算法。 对于问题四，首先将视频材料以1min为间隔进行采样得到41帧（静态图片），将这些静止图片先利用matlab进行处理，后进行阀值归一化处理，得到这些帧的灰度值矩阵。在图片上建立参考模型，获得影子端点的参考位置。利用投影系统和模型二，建立了基于图形处理的视频拍摄地点搜索模型。利用模型二中多层搜索算法，求得满足精度的最优地点。最优的地点是：东经119,北纬48.7，在内蒙古的呼伦贝尔市。同时假设日期是未知量，将模型四与模型三相结合，得到了可能的地点和时间，并分析了可能出现误差的原因，最后回答了当视频日期未知，也可以确定其位置和日期。 最后，给出了模型的优缺点和改进方案。 关键词：极坐标化，多层优化搜索算法，遗传算法，图像处理，MATLAB
历年全国数学建模试题及解法归纳 赛题解法 93A非线性交调的频率设计拟合、规划 93B足球队排名图论、层次分析、整数规划94A逢山开路图论、插值、动态规划 94B锁具装箱问题图论、组合数学 95A飞行管理问题非线性规划、线性规划 95B天车与冶炼炉的作业调度动态规划、排队论、图论96A最优捕鱼策略微分方程、优化 96B节水洗衣机非线性规划 97A零件的参数设计非线性规划 97B截断切割的最优排列随机模拟、图论 98A一类投资组合问题多目标优化、非线性规划98B灾情巡视的最佳路线图论、组合优化 99A自动化车床管理随机优化、计算机模拟 99B钻井布局0-1规划、图论 00A DNA序列分类模式识别、Fisher判别、人工 神经网络 00B钢管订购和运输组合优化、运输问题 01A血管三维重建曲线拟合、曲面重建
赛题解法 01B 公交车调度问题多目标规划 02A车灯线光源的优化非线性规划 02B彩票问题单目标决策 03A SARS的传播微分方程、差分方程 03B 露天矿生产的车辆安排整数规划、运输问题 04A奥运会临时超市网点设计统计分析、数据处理、优化04B电力市场的输电阻塞管理数据拟合、优化 05A长江水质的评价和预测预测评价、数据处理 05B DVD在线租赁随机规划、整数规划 06A出版社书号问题整数规划、数据处理、优化06B Hiv病毒问题线性规划、回归分析 07A 人口问题微分方程、数据处理、优化07B 最佳交通线路查询多目标规划、图论 08A 照相机问题非线性方程组、优化 08B 大学学费问题数据收集和处理、统计分 析、回归分析 09A制动器试验台的控制方法分析物理模型，计算机仿真09B 眼科病房的合理安排综合评价，决策与预测 10A储油罐的变位识别与罐容标定微积分理论，数值计算10B 2010上海世博会影响力的评价综合评价，统计分析 11A城市表层重金属污染分析综合评价，统计分析
太阳影子定位 （一）摘要 根据影子的形成原理和影子随时间的变化规律，可以建立时间、太阳位置和影子轨迹的数学模型，利用影子轨迹图和时间可以推算出地点等信息，从而进行视频数据分析可以确定视频的拍摄地点。本文根据此模型求解确定时间地点影子的运动轨迹和对于已知运动求解地点或日期。 直立杆的影子的位置在一天中随太阳的位置不断变化，而其自身的所在的经纬度以及时间都会影响到影子的变化。但是影子的变化是一个连续的轨迹，可以用一个连续的函数来表达。我们可以利用这根长直杆顶端的影子的变化轨迹来描述直立杆的影子。众所周知，地球是围绕太阳进行公转的，但是我们可以利用相对运动的原理，将地球围绕太阳的运动看成是太阳围绕地球转动。 我们在解决问题一的时候，利用题目中所给出的日期、经纬度和时间，来解出太阳高度角ｈ，太阳方位角Α，赤纬角δ，时角Ω，直杆高度Ｈ和影子端点位置（x0，y o），从而建立数学模型。影子的端点坐标是属于时间的函数，所以可以借助时间写出参数方程来描述影子轨迹的变化。问题二中给出了日期和随时间影子端点的坐标变化，可以根据坐标变化求出运用软件拟合出曲线找到在正午时纵坐标最小，横坐标最大，影子最短的北京时间，根据时差与经度的关系，求出测量地点的经度。根据太阳方位角Α，赤纬角δ，时角Ω，可以求出太阳高度角ｈ。再结合问题一中的表达式，建立方程求解测量地点的纬度Ф。我们在求解第三问的思路也是沿用之间的模型，但第三问上需要解出日期。 对于问题四的求解，先获取自然图像序列或者视频帧，并对每一帧图像检测出影子的轨迹点；然后确定多个灭点，并拟合出地平线；拟合互相垂直的灭点，计算出仿射纠正和投影纠正矩阵；进而还原出经过度量纠正的世界坐标；在拟合出经过度量纠正世界坐标中的影子点的轨迹，利用前面几问中的关系求出经纬度。 关键字：太阳影子轨迹Matlab 曲线拟合