2013年数学建模b题

精心整理

碎纸片的拼接复原

【摘要】：碎纸片拼接技术是数字图像处理领域的一个重要研究方向，把计算机视觉和程序识别应用于碎纸片的复原，在考古、司法、古生物学等方面具有广泛的应用，具有重要的现实意义。本文主要结合各种实际应用背景，针对碎纸机绞碎的碎纸片，基于计算机辅助对碎纸片进行自动拼接复原研究。

针对问题1，依据图像预处理理论，通过matlab程序处理图像，将图像转化成适合于计算机处理的数字图像，进行灰度分析，提取灰度矩阵。对于仅纵切的碎纸片，根据矩阵的行提取理论，将

。建

中的任一列与矩阵值，

序列号。将程序进行循环操作，得到最终的碎片自动拼接结果。

、；

分别作为新生成的矩阵、。

，

将矩阵中的任一列分别与矩阵中每一列代入模型，所得p值对应的值即为横排序；将矩

阵中的任一行分别于矩阵中的任一行代入模型，所得q值对应的值即为列排序。循环进行此

程序，得计算机的最终运行结果。所得结果有少许误差，需人工调制，更正排列顺序，得最终拼接结果。

针对问题3，基于碎纸片的文字行列特征，采用遗传算法，将所有的可能性拼接进行比较，进行择优性选择。反面的排序结果用于对正面排序的检验，发现结果有误差，此时，进行人工干预，调换碎纸片的排序。

【关键词】:灰度矩阵欧式距离图像匹配自动拼接人工干预

一、问题重述

破碎文件的拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。传统上，大量的纸质物证复原工作都是以人工的方式完成的，准确率较高，但效率很低。特别是当碎片数量巨大，人工拼接不但耗费大量的人力、物力，而且还可能对物证造成一定的损坏。随着计算机技术的发展，人们试图把计算机视觉和模式识别应用于碎纸片复原，开展对碎纸片自动拼接技术的研究，以提高拼接复原效率。试讨论一下问题，并根据题目要求建立相应的模型和算法：

、附件4

(1)

(2)

(3)

(4)

纸片的自动拼接。

问题1：

根据图像预处理理论，通过程序语言将图像导入matlab程序，对图像进行预处理，将碎纸片转换成适合于计算机处理的数字图像形式，并对数字图像进行灰度分析，提取灰度矩阵。对于仅纵切的碎纸片，根据矩阵的行提取理论，将每个灰度矩阵的第一列提取，作为新矩阵A，提取每个灰度矩阵的最后一列，作为新矩阵B。建立碎纸片匹配模型。根据遗传算法的思想，将矩阵A中的任

一行与矩阵B中的每一行带入模型，由最大相似性原理得s值对应的值，即为所拼接的碎片序列号。将程序进行循环操作，并最终建立对输入的解释。

问题2：

首先对碎纸片图像进行预处理，通过matlab程序将图像信息进行灰度分析，提取灰度矩阵。基于既纵切又横切的碎纸片，根据矩阵的行列提取理论，分别提取每个灰度矩阵的第一列，生成新

矩阵，提取最后一列，结合生成新矩阵；提取所有灰度矩阵的第一行，作为新生成的矩阵，

提取最后一行，作为新矩阵。由于纸质文件边缘空白处的灰度值为常量，通过对灰度矩阵的检

验提取，确定最左列的碎纸片排序。在此基础上，采用从局部到整体，从左到右的方法，建立匹配筛选模型。将矩阵中的任一列分别与矩阵中每一列代入模型，所得p值对应的值即为横排

值对应的值即为列排序。循

运

,采用

D

中第列与矩阵列的欧氏距离

中第中第

四、模型假设

1）同一切线处的两条边界灰度值的差别在误差范围之内。

2）所有的待拼接的碎纸片都是规则的。

3）假设不存在文件缺失、人为因素等意外因素。

4）假设所有的数字图像文件都是已整理好的，不需要再进行旋转、剪裁等操作。

五、模型建立与问题解决

首先，对碎纸片图像进行预处理。由于基于文字的数字图像中文字和背景的灰度值存在明显差异，运用程序把已整理好的数字图像依次导入matlab软件，根据图像显示技术，提取其灰度矩阵，

改变数字图像的视觉效果，将数字图像转化成一种更适合于计算机进行分析处理的数据形式。再者，本文涉及的是无重叠的规则碎纸片拼接，因此只需处理每一张碎纸片的边缘交界处的数据，应用matlab程序对灰度矩阵进行边界数据的提取，减少了数据的处理量，为碎纸片的自动拼接节省了时间。然后，对所提取的每组数据进行相似性度量，采用二维欧式距离，即：

根据遗传算法，将提取的所有待拼接碎纸片的边界数据进行分别左右、上下交叉的两两匹配，根据最大相似性原则，建立匹配准则：

值即为拼接序列。

，生成新矩阵

建立碎纸片匹配模型:

将矩阵中的任一列与矩阵中的每一列带入模型，值，

纸片，根据矩阵的行列提取理论，分别提取每个灰度矩阵的第一列，生成新矩阵，提取最后一

列，结合生成新矩阵；提取所有灰度矩阵的第一行，作为新生成的矩阵，提取最后一行，作为

新矩阵。由于纸质文件边缘空白处的灰度值为常量，通过对灰度矩阵的检验提取，确定最左列的碎纸片排序。在此基础上，采用从局部到整体，从左到右的方法，建立匹配筛选模型：

，

将矩阵中的任一列分别与矩阵中每一列代入模型，所得p值对应的值即为横排序；将矩

阵中的任一行分别于矩阵中的任一行代入模型，所得q值对应的值即为列排序。循环进行此

程序，得计算机的最终运行结果。所得结果有少许误差，需人工调制，更正排列顺序，得最终拼接结果。

针对双面打印文件的碎纸片拼接处理，首先要对碎纸片两面的图像信息都进行图像预处理，运用matlab程序提取图像的灰度矩阵。基于碎纸机绞碎的碎纸片，每个碎纸片的大小是一样的,采用差距度量法。取矩阵中特征值从第一行中开始到限制条件取i行，限制条件是当取到i行时满足每行都是特征值的矩阵还有2m个矩阵，则这2m个矩阵为同一行碎片所产生的矩阵。包括正反面碎纸片。然后取其中矩阵左边列向量满足，此行矩阵所剩两个矩阵，此过程虽然误差较小，针对其他问题可能要产生误差，需要人工干预，剔除多余矩阵。然后取出两个中的一个矩阵，与产生的这一行

矩阵进行匹配，获取一行排序较优的碎片链接。运用嵌套循环实现每行的排序连接。

附件5的碎纸片复原结果为：

正面结果：

从实验分析可看出，边缘灰度值的精度对减少匹配的计算量和提高匹配的准确性都非常重要。它在很大程度上决定了纸片匹配的排序问题，有时甚至会导致某些匹配的失败。

2）数据库的开发管理

在实际应用中，由于设备等条件的限制，一次的处理往往不能解决全部的拼接问题。因此，必须开发合理的数据库，以充分利用每一次的拼接结果。

3）程序优化

虽然在进行简单的拼接复原工作时，该系统可以满足要求。但是在实际的应用中，通常会涉及对大量不规则纸片或彩色碎片数据进行管理和处理工作。因此，必须继续优化程序结构，提高程序的交互能力，真正实现快速有效的计算机辅助碎纸片自动拼接复原。

本文基于计算机辅助对碎纸机绞碎的碎纸片实现自动快速拼接，从图形显示技术中提取图形的

灰度信息，并提出欧氏距离的配准方法，对图像中的误匹配点产生极大的抑制，求解的变换矩阵精度得到提高，大大缩短了匹配时间，在现实应用中具有重要意义，类似的研究可广泛应用到文物碎片的自动修复、司法物证修复、计算机辅助设计等领域。

七、参考文献

[1]WolfsonH，KalvinA，SchonbergE，LambdanY，Solvingjigsawpuzzlesbycomputer，AnnalesofOperationResearch，1988，12:51-64

[2]FreemanH，GarderL，Apictorialjigsawpuzzles:thecomputersolutionofaprobleminpatternrecognition，IEEE，1964,13:118-127

[3]余宏生，金伟其，散字图像拼接方法研究进展口[J]，红外技术，2009,31（6）：348-353

附件

附件

s(:,:,i)=imread(files(i).name);

end;%将所有碎片导入一个三维数组

forj=1:19

n=0;

fori=1:1980

ifs(i,1,j)==255

n=n+1;

end;

end;

ifn==1980

end;

end;

clc

cell{7}=imread('006.bmp'); cell{8}=imread('007.bmp'); cell{9}=imread('008.bmp'); cell{10}=imread('009.bmp'); cell{11}=imread('010.bmp');

cell{12}=imread('011.bmp');

cell{13}=imread('012.bmp');

cell{14}=imread('013.bmp');

cell{15}=imread('014.bmp');

cell{16}=imread('015.bmp');

end

t=8

j=1;

end

x1=sum(sqrt((cell2{t}(:,72)-cell2{i}(:,1)).'*(cell2{t}(:,72)-cell2{i}(:,1)))); ifx1

x=x1;

j=i;

end j

t=j; end 结果

t= 8

j= 18 j= 1

j= 7

j= 9

j= 15 j= 13 j= 16 j=

j=

11

j=

3

j=

17

j=

2

j=

5

j=

6

j=

10

j=

14

j=

1

j=

7

figure2

IV1=imread('008.bmp');

IV2=imread('014.bmp');

IV3=imread('012.bmp');

IV4=imread('015.bmp');

IV5=imread('003.bmp');

IV6=imread('010.bmp');

PicData=[IV1,IV2,IV3,IV4,IV5,IV6,IV7,IV8,IV9,IV10,IV11,IV12,IV13,IV14,IV15,IV16, IV17,IV18,IV19];

imshow(PicData);

附件2figure1

clearall

cell{1}=imread('000.bmp'); cell{2}=imread('001.bmp'); cell{3}=imread('002.bmp'); cell{4}=imread('003.bmp');

cell{18}=imread('017.bmp'); cell{19}=imread('018.bmp'); fori=1:19

cell2{i}=im2bw(cell{i},0.5); end

forx=1:18

x=sum(abs(cell2{t}(:,72)-cell2{1}(:,1))); j=1;

fori=1:19

end

j=i;

end

end

j

t=j;

end

IV1=imread('003.bmp');

IV2=imread('006.bmp');

IV3=imread('002.bmp');

IV4=imread('007.bmp');

IV5=imread('015.bmp');

IV6=imread('018.bmp');

IV7=imread('011.bmp');

IV8=imread('000.bmp');

IV9=imread('005.bmp');

IV10=imread('001.bmp');

附件

clear

>>filename=dir('C:\Users\127\Desktop\附件3\*.bmp');%文件路径t=209;

cell(1,t);

I=cell(1,t);

name=imread(strcat('C:\Users\127\Desktop\附件3\',filename(i).name));%文件路径I{i}=name;

end

c=cell2mat(I);

结果

A=

clear

>>filename=dir('C:\Users\127\Desktop\附件3\*.bmp');%文件路径

t=209;

cell(1,t);

I=cell(1,t);

name=imread(strcat('C:\Users\127\Desktop\附件3\',filename(i).name));%文件路径I{i}=name;

end

c=cell2mat(I);

end

end

b=C'

HIMCM 2014美国中学生数学建模竞赛试题

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

2017全国数学建模竞赛B题

2017年高教社杯全国大学生数学建模竞赛题目 （请先阅读“全国大学生数学建模竞赛论文格式规范”） B题“拍照赚钱”的任务定价 “拍照赚钱”是移动互联网下的一种自助式服务模式。用户下载APP，注册成为APP的会员，然后从APP上领取需要拍照的任务（比如上超市去检查某种商品的上架情况），赚取APP对任务所标定的酬金。这种基于移动互联网的自助式劳务众包平台，为企业提供各种商业检查和信息搜集，相比传统的市场调查方式可以大大节省调查成本，而且有效地保证了调查数据真实性，缩短了调查的周期。因此APP成为该平台运行的核心，而APP中的任务定价又是其核心要素。如果定价不合理，有的任务就会无人问津，而导致商品检查的失败。 附件一是一个已结束项目的任务数据，包含了每个任务的位置、定价和完成情况（“1”表示完成，“0”表示未完成）；附件二是会员信息数据，包含了会员的位置、信誉值、参考其信誉给出的任务开始预订时间和预订限额，原则上会员信誉越高，越优先开始挑选任务，其配额也就越大（任务分配时实际上是根据预订限额所占比例进行配发）；附件三是一个新的检查项目任务数据，只有任务的位置信息。请完成下面的问题： 1.研究附件一中项目的任务定价规律，分析任务未完成的原因。 2.为附件一中的项目设计新的任务定价方案，并和原方案进行比较。 3.实际情况下，多个任务可能因为位置比较集中，导致用户会争相选择，一种 考虑是将这些任务联合在一起打包发布。在这种考虑下，如何修改前面的定价模型，对最终的任务完成情况又有什么影响？ 4.对附件三中的新项目给出你的任务定价方案，并评价该方案的实施效果。 附件一：已结束项目任务数据 附件二：会员信息数据 附件三：新项目任务数据

数学建模练习试题

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万元；

2016年数学建模大赛试题B题

2016高教社杯全国大学生数学建模竞赛题目（请先阅读“全国大学生数学建模竞赛论文格式规范”） B题小区开放对道路通行的影响 2016年2月21日，国务院发布《关于进一步加强城市规划建设管理工作的若干意见》，其中第十六条关于推广街区制，原则上不再建设封闭住宅小区，已建成的住宅小区和单位大院要逐步开放等意见，引起了广泛的关注和讨论。 除了开放小区可能引发的安保等问题外，议论的焦点之一是：开放小区能否达到优化路网结构，提高道路通行能力，改善交通状况的目的，以及改善效果如何。一种观点认为封闭式小区破坏了城市路网结构，堵塞了城市“毛细血管”，容易造成交通阻塞。小区开放后，路网密度提高，道路面积增加，通行能力自然会有提升。也有人认为这与小区面积、位置、外部及内部道路状况等诸多因素有关，不能一概而论。还有人认为小区开放后，虽然可通行道路增多了，相应地，小区周边主路上进出小区的交叉路口的车辆也会增多，也可能会影响主路的通行速度。 城市规划和交通管理部门希望你们建立数学模型，就小区开放对周边道路通行的影响进行研究，为科学决策提供定量依据，为此请你们尝试解决以下问题： 1. 请选取合适的评价指标体系，用以评价小区开放对周边道路通行的影响。 2. 请建立关于车辆通行的数学模型，用以研究小区开放对周边道路通行的影响。交通流分配模型 3. 小区开放产生的效果，可能会与小区结构及周边道路结构、车流量有关。请选取或构建不同类型的小区，应用你们建立的模型，定量比较各类型小区开放前后对道路通行的影响。 4. 根据你们的研究结果，从交通通行的角度，向城市规划和交通管理部门提出你们关于小区开放的合理化建议。

西南大学2016年春《数学建模》作业及答案(已整理)(共5次)

西南大学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．概

2015年全国数学建模B题论文思路

B题“互联网+”时代的出租车资源配置 出租车是市民出行的重要交通工具之一，“打车难”是人们关注的一个社会热点问题。随着“互联网+”时代的到来，有多家公司依托移动互联网建立了打车软件服务平台，实现了乘客与出租车司机之间的信息互通，同时推出了多种出租车的补贴方案。 请你们搜集相关数据，建立数学模型研究如下问题： (1)试建立合理的指标，并分析不同时空出租车资源的“供求匹配”程度。指标：里程利用率，车辆满载率，车辆拥有量（万人）等，从这些指标去按以下步骤收集数据并分析 1分别收集一线（比如北上广），二线（比如西安），三线（比如拉萨）城市各一个的出租车数据来分析，这样就能代表全国了。这就是第一问中的“空” 2主要分析各个城市早（7:00——8:30） 中（11:30——2:30） 晚（17:30——18:30）上班高峰 和平时时段的打车的供求情况这就是第一问中的“时” 3最后总结哈供求匹配程度

(2)分析各公司的出租车补贴方案是否对“缓解打车难”有帮助？ 1选取几个打车平台的补贴方案去分析，比如： 快的打车补贴变化 2014年1月20日快的打车乘客车费返现10元，司机奖励10元 2014年2月17日快的打车乘客返现11元，司机返5-11元[10] 2014年2月18日快的打车乘客返现13元[11] 2014年3月4日快的打车乘客返现10元/单，司机端补贴不变[6] 2014年3月5日快的打车乘客补贴金额变为5元 2014年3月22日快的打车乘客返现3—5元 2014年5月17日软件乘客补贴“归零” 2014年7月9日，将司机端补贴降为2元/单。[12] 2014年8月9日，滴滴、快的两大打车软件再出新规，全面取消司机端现金补贴。 滴滴打车 1月10日，滴滴打车乘客车费立减10元、司机立奖10元 2月17日，滴滴打车乘客返现10-15元，新司机首单立奖50元 2月18日，滴滴打车乘客返现12至20元 3月7日，滴滴打车乘客每单减免随机“6-15元” 3月23日，滴滴打车乘客返现3-5元 5月17日，打车软件乘客补贴“归零” 7月9日，软件司机端补贴降为2元/单 8月12日，滴滴打车取消对司机接单的常规补贴 2分析传统出租车公司的补贴方案 3最后一定要联系到是否对“缓解打车难”有帮助上，结论是：有一定帮助，但并未完全解决问题（），同时产生了新的问题。 注意要用数据和案例论证，不能自己在那空口说。这样就为下

数学建模作业

习 题 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

2011年全国大学生数学建模竞赛B题

承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。 我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员(打印并签名) ：1. 韩晓峰 2. 杨晓帆 3. 李弘倩 指导教师或指导教师组负责人(打印并签名)： 日期： 2011 年 9 月 11 日赛区评阅编号（由赛区组委会评阅前进行编号）：

编号专用页 赛区评阅编号（由赛区组委会评阅前进行编号）： 赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号（由赛区组委会送交全国前编号）： 全国评阅编号（由全国组委会评阅前进行编号）：

交巡警服务平台的设置与调度 摘要 在（1）第一问中，我们根据附表1所给各路口坐标算出A图中每条路线的长度，然后通过floyd算法找出了两点之间的最短路程，得出矩阵D，通过使用matlab圈出各服务平台到周围路口小于3min（即3km）的点，再根据就近原则，将各路口划分到这个圈中离此路口最近的交巡警平台。对于任意到交巡警平台路程大于3min（即3km）序号为28，29，38，39，61，92的五个路口，则采用就近原则人工划入距离其最近的交巡警平台辖区，这样就在保证出警时间基本都小于3min的条件下，划分出各警务平台合理的管辖范围。 对于（1）第二问中，我们采用指派模型，用lingo软件对20个巡警服务平台对17个城市出入口进行封锁的方法进行了优化，得到初步的调度方案。在这个方案的基础上，如果在某条巡警服务平台调度路线中经过其他的的调度点，则与所经过的调度点互换目标路口，由此得到最佳调度方案，即最快8分钟可以实现快速封锁路口。 对于（1）第三问，我们按照工作量均衡和出警时间尽可能短的原则考虑增加交巡警平台。首先，利用excel算出了各辖区内交巡警服务平台的工作量总和以及出警时间的平均值，求出了所有辖区的工作总量与出警时间的平均值，取出工作总量明显高于平均值且平均出警时间超过3min的四个辖区来增设新的交巡警平台。我们通过（1）第一问中的最短路程矩阵得到各辖区内种增设平台的所有可能的方法，通过比较每种方法的工作量总和及出警时间综合考虑得到一个最优的设置。最后通过spss软件求得优化前后两组数据的标准差，比较后发现优化后的数据标准差明显下降，达到了优化的目的。 在（2）第一问中，首先参照（1）中A区的处理方法分别求出了城区B，C，D，E，F中路口到最近交巡警服务平台的路程，出警时间以及工作量。通过spss软件计算出标准差，并与优化后的A城区进行比较，找出了交巡警平台明显分配不合理（标准差远高于A区优化后方案的标准差）的三个区，在其距离较远工作量较大的路口增设服务平台。 对于（2）中第二问搜捕嫌疑犯问题，我们采用时间圈法，以事发地点为圆心，以一分钟所行使的路程为单位半径，画出等间隔的数个同心圆，这样就可以确定每个时间段嫌疑犯的逃逸范围（因为任意两点间的路程大于等于两点间的距离，所以在某时间内嫌疑犯必定在对应的时间圈内），根据每个路口在这些时间圈上的位置，要求警察到这些路口的时间小于罪犯到达的时间（即可以围堵住罪犯），找出所有这些路口，构成闭圈，最小的闭圈便是围堵罪犯的最佳方案。 关键词：floyd算法，指派模型，spss分析数据，时间圈法

2014年美国数学建模大赛(MCM)试题译文

2014年美国数学建模大赛（MCM）试题译文 王景璟大连理工大学 问题A：超车之外靠右行原则 在一些开车必须靠右行驶的国家（比如：美国，中国，以及其他除了英国，澳大利亚，和一些前英国殖民地的国家），行驶在多车道高速路必须遵循一个规则，那就是要求驾驶员在超车之外的情况下，必须在最靠右的车道行驶，超车时，他们向左变道，超车，然后再回到之前行驶的车道。 构建一个数学模型来分析该规则在车流量很少和很大的时候的执行情况。你最好能考察车流量与安全的之间的相互关系，过低或是过量的速度限制的作用（速度设置过低或是过高），以及/或者其他在该问题陈述中没有明确提到的因素。该原则是否能有效促进更好的车流量？如果无效，请建议和分析其他更有助于提高车流量、安全、以及其他你认为重要的因素的其他方案（可以完全不包括该原则）。 在开车靠左行的国家，讨论一下你的方案在经过对方向的简单修改之后或是添加额外的要求之后是否也适用。 最后，以上原则取决于人们遵循交通规则的判断力。如果道路上的车流完全在智能系统（要么是道路体系的一部分，要么是包含在使用道路的所有车辆的设计之中）的控制之下，该改变在多大程度上会影响你先前分析的结果？ 问题B: 大学教练联盟 《体育画报》，一本体育爱好者的杂志，正在寻找上世纪“最好的大学教练”，包括男性和女性。建立一个数学模型以从诸如大学曲棍球，曲棍球，橄榄球，棒球，垒球，篮球，或足球等运动的男性或女性教练中选出最好的一个教练或几个教练（过去的或现在的）。分析中使用的时间分界线是否有影响？即在1913执教和在2013年执教有不同吗？清晰地表达你们模型中的评判标准。讨论你们的模型如何能广泛地应用于两种性别及所有可能的体育运动。分别选出你模型中3种不同运动的前5位教练。 除了MCM格式及要求，准备一篇1-2页的文章给《体育画报》以解释你们的结论并包括一份能让体育迷们看懂的对你们数学模型的非技术性解释。 问题C：使用网络模型测量影响力

2016年高教社杯全国大学生数学建模竞赛题目A.B

2016年高教社杯全国大学生数学建模竞赛题目 （请先阅读“全国大学生数学建模竞赛论文格式规范”） A题系泊系统的设计 近浅海观测网的传输节点由浮标系统、系泊系统和水声通讯系统组成（如图1所示）。某型传输节点的浮标系统可简化为底面直径2m、高2m的圆柱体，浮标的质量为1000kg。系泊系统由钢管、钢桶、重物球、电焊锚链和特制的抗拖移锚组成。锚的质量为600kg，锚链选用无档普通链环，近浅海观测网的常用型号及其参数在附表中列出。钢管共4节，每节长度1m，直径为50mm，每节钢管的质量为10kg。要求锚链末端与锚的链接处的切线方向与海床的夹角不超过16度，否则锚会被拖行，致使节点移位丢失。水声通讯系统安装在一个长1m、外径30cm的密封圆柱形钢桶内，设备和钢桶总质量为100kg。钢桶上接第4节钢管，下接电焊锚链。钢桶竖直时，水声通讯设备的工作效果最佳。若钢桶倾斜，则影响设备的工作效果。钢桶的倾斜角度（钢桶与竖直线的夹角）超过5度时，设备的工作效果较差。为了控制钢桶的倾斜角度，钢桶与电焊锚链链接处可悬挂重物球。 图1 传输节点示意图（仅为结构模块示意图，未考虑尺寸比例）系泊系统的设计问题就是确定锚链的型号、长度和重物球的质量，使得浮标的吃水深度和游动区域及钢桶的倾斜角度尽可能小。 问题1某型传输节点选用II型电焊锚链22.05m，选用的重物球的质量为1200kg。

现将该型传输节点布放在水深18m、海床平坦、海水密度为1.025×103kg/m3的海域。若海水静止，分别计算海面风速为12m/s和24m/s时钢桶和各节钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。 问题2在问题1的假设下，计算海面风速为36m/s时钢桶和各节钢管的倾斜角度、锚链形状和浮标的游动区域。请调节重物球的质量，使得钢桶的倾斜角度不超过5度，锚链在锚点与海床的夹角不超过16度。 问题3 由于潮汐等因素的影响，布放海域的实测水深介于16m~20m之间。布放点的海水速度最大可达到1.5m/s、风速最大可达到36m/s。请给出考虑风力、水流力和水深情况下的系泊系统设计，分析不同情况下钢桶、钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。 说明近海风荷载可通过近似公式F=0.625×Sv2(N)计算，其中S为物体在风向法平面的投影面积(m2)，v为风速(m/s）。近海水流力可通过近似公式F=374×Sv2(N)计算，其中S为物体在水流速度法平面的投影面积(m2)，v为水流速度(m/s）。

2016年全国大学生数学建模B题思路

2016 高教社杯全国大学生数学建模竞赛题目 、 B 题分析初稿，旨在交流，有各种做题思路，大家自由发挥！ 注意：这只是看了3篇文章，找到的思路，请大家多看文献，思路会很多！ 我们后续会整理更多的思路！ 关键词： 1.评价指标体系，评价开放对周边道路通行的效果。 2.车辆通行的数学模型，研究小区开放对周边道路通行的影响。 3.小区开放产生的效果，可能会与小区结构及周边道路结构、车流量有关。请选取或构建不同类型的小区，应用你们建立的模型，定量比较各类型小区开放前后对道路通行的影响。 4. 根据你们的研究结果，从交通通行的角度，向城市规划和交通管理部 门提出你们关于小区开放的合理化建议。 相关资料整理： 1.评价指标体系，评价开放对周边道路通行的效果。

参考文献《居住小区开发交通影响分析研究_商仲华》第 48 页，有 5 个指标，并用层次分析 AHP 进行了研究。 我们要做的可能是强调类似哪些指标是针对开放对周边道路通行的效果，不属于这类的指标可以删除。 2.车辆通行的数学模型，研究小区开放对周边道路通行的影响。 参考文献《城市交通拥堵对策_封闭型小区交通开放研究_李向朋》第 11 页，图 6 上面，给出一句话，关于开放小区的定义。 是不是建模就是选取小区附件的某些范围研究，这就是理论依据。 参考文献《城市交通拥堵对策_封闭型小区交通开放研究_李向朋》第 26 页，图 3.2，了解道路系统的简图，用简图做分析。

简单的车辆模型，可以化个节点，图，权重。分析流量。类似文献《城市应急车辆优先通行关键问题研究_毕煦东》第 23 页， 用其中的符号定义等，后面的应急什么别管，太复杂。利用这里模型分析第一个问题中指标系统的指标。 3.小区开放产生的效果，可能会与小区结构及周边道路结构、车流量有关。请选取或构建不同类型的小区，应用你们建立的模型，定量比较各类型小区开放前后对道路通行的影响。 小区结构：参考文献《城市交通拥堵对策_封闭型小区交通开放研究_李向朋》第10页，

2013数学建模B题国家一等奖Matlab程序

附录3：程序源文件 1.duqu_image.m文件 %数据读取预处理文件 %将附件中的图片读取到matlab矩阵中，并保存为image_1，image_2，image_3，image_4，image_5a，image_5b %所有附件均放在文件夹 D:\B 中%%%%%%%%%%%%%%%%%%%%%%%%图片名序列 %图像名称序号 b = [ones(1,10);0 1 2 3 4 5 6 7 8 9]'; image_num= [ strcat( num2str(0*b(:,1)),num2str(0*b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(2*b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(3*b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(4*b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(5*b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(6*b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(7*b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(8*b(:,1)),num2str(b(:,2)) ), strcat( num2str(0*b(:,1)),num2str(9*b(:,1)),num2str(b(:,2)) ), strcat( num2str(b(:,1)),num2str(0*b(:,1)),num2str(b(:,2)) ), strcat( num2str(b(:,1)),num2str(b(:,1)),num2str(b(:,2)) ), strcat( num2str(b(:,1)),num2str(2*b(:,1)),num2str(b(:,2)) ), strcat( num2str(b(:,1)),num2str(3*b(:,1)),num2str(b(:,2)) ), strcat( num2str(b(:,1)),num2str(4*b(:,1)),num2str(b(:,2)) ), strcat( num2str(b(:,1)),num2str(5*b(:,1)),num2str(b(:,2)) ), strcat( num2str(b(:,1)),num2str(6*b(:,1)),num2str(b(:,2)) ),

2014年下学期数学实验与数学建模作业习题8

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)

全国数学建模2004年B题

2004高教社杯全国大学生数学建模竞赛题目 B题电力市场的输电阻塞管理 我国电力系统的市场化改革正在积极、稳步地进行。2003年3月国家电力监管委员会成立，2003年6月该委员会发文列出了组建东北区域电力市场和进行华东区域电力市场试点的时间表，标志着电力市场化改革已经进入实质性阶段。可以预计，随着我国用电紧张的缓解，电力市场化将进入新一轮的发展，这给有关产业和研究部门带来了可预期的机遇和挑战。 电力从生产到使用的四大环节——发电、输电、配电和用电是瞬间完成的。我国电力市场初期是发电侧电力市场，采取交易与调度一体化的模式。电网公司在组织交易、调度和配送时，必须遵循电网“安全第一”的原则，同时要制订一个电力市场交易规则，按照购电费用最小的经济目标来运作。市场交易-调度中心根据负荷预报和交易规则制订满足电网安全运行的调度计划――各发电机组的出力（发电功率）分配方案；在执行调度计划的过程中，还需实时调度承担AGC（自动发电控制）辅助服务的机组出力，以跟踪电网中实时变化的负荷。 设某电网有若干台发电机组和若干条主要线路，每条线路上的有功潮流（输电功率和方向）取决于电网结构和各发电机组的出力。电网每条线路上的有功潮流的绝对值有一安全限值，限值还具有一定的相对安全裕度（即在应急情况下潮流绝对值可以超过限值的百分比的上限）。如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值，称为输电阻塞。当发生输电阻塞时，需要研究如何制订既安全又经济的调度计划。 ●电力市场交易规则： 1. 以15分钟为一个时段组织交易，每台机组在当前时段开始时刻前给出下一个时段的报价。各机组将可用出力由低到高分成至多10段报价，每个段的长度称为段容量，每个段容量报一个价（称为段价），段价按段序数单调不减。在最低技术出力以下的报价一般为负值，表示愿意付费维持发电以避免停机带来更大的损失。 2. 在当前时段内，市场交易-调度中心根据下一个时段的负荷预报,每台机组的报价、当前出力和出力改变速率，按段价从低到高选取各机组的段容量或其部分（见下面注释），直到它们之和等于预报的负荷，这时每个机组被选入的段容量或其部分之和形成该时段该机组的出力分配预案（初始交易结果）。最后一个被选入的段价（最高段价）称为该时段的清算价，该时段全部机组的所有出力均按清算价结算。 注释： （a）每个时段的负荷预报和机组出力分配计划的参照时刻均为该时段结束时刻。 （b）机组当前出力是对机组在当前时段结束时刻实际出力的预测值。 （c）假设每台机组单位时间内能增加或减少的出力相同，该出力值称为该机组的爬坡速率。由于机组爬坡速率的约束，可能导致选取它的某个段容量的部分。 （d）为了使得各机组计划出力之和等于预报的负荷需求，清算价对应的段容量可能只选取部分。 市场交易-调度中心在当前时段内要完成的具体操作过程如下： 1、监控当前时段各机组出力分配方案的执行，调度AGC辅助服务，在此基础上给出各机组的当前出力值。 2、作出下一个时段的负荷需求预报。 3、根据电力市场交易规则得到下一个时段各机组出力分配预案。 4、计算当执行各机组出力分配预案时电网各主要线路上的有功潮流，判断是否会出现输电阻塞。如果不 出现，接受各机组出力分配预案；否则，按照如下原则实施阻塞管理： ●输电阻塞管理原则： （1）调整各机组出力分配方案使得输电阻塞消除。 （2）如果（1）做不到，还可以使用线路的安全裕度输电，以避免拉闸限电（强制减少负荷需求），但要使每条线路上潮流的绝对值超过限值的百分比尽量小。 （3）如果无论怎样分配机组出力都无法使每条线路上的潮流绝对值超过限值的百分比小于相对安全裕度，则必须在用电侧拉闸限电。

2013数学建模国赛B题

3v2013高教社杯全国大学生数学建模竞赛题目 （请先阅读“全国大学生数学建模竞赛论文格式规范”） B题碎纸片的拼接复原 破碎文件的拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。传统上，拼接复原工作需由人工完成，准确率较高，但效率很低。特别是当碎片数量巨大，人工拼接很难在短时间内完成任务。随着计算机技术的发展，人们试图开发碎纸片的自动拼接技术，以提高拼接复原效率。请讨论以下问题： 1. 对于给定的来自同一页印刷文字文件的碎纸机破碎纸片（仅纵切），建立碎纸片拼接复原模型和算法，并针对附件1、附件2给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结果以图片形式及表格形式表达（见【结果表达格式说明】）。 2. 对于碎纸机既纵切又横切的情形，请设计碎纸片拼接复原模型和算法，并针对附件3、附件4给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结果表达要求同上。 3. 上述所给碎片数据均为单面打印文件，从现实情形出发，还可能有双面打印文件的碎纸片拼接复原问题需要解决。附件5给出的是一页英文印刷文字双面打印文件的碎片数据。请尝试设计相应的碎纸片拼接复原模型与算法，并就附件5的碎片数据给出拼接复原结果，结果表达要求同上。 【数据文件说明】 (1)每一附件为同一页纸的碎片数据。 (2)附件1、附件2为纵切碎片数据，每页纸被切为19条碎片。 (3)附件3、附件4为纵横切碎片数据，每页纸被切为11×19个碎片。 (4)附件5为纵横切碎片数据，每页纸被切为11×19个碎片，每个碎片有正反两面。该附 件中每一碎片对应两个文件，共有2×11×19个文件，例如，第一个碎片的两面分别对应文件000a、000b。 【结果表达格式说明】 复原图片放入附录中，表格表达格式如下： (1)附件1、附件2的结果：将碎片序号按复原后顺序填入1×19的表格； (2)附件3、附件4的结果：将碎片序号按复原后顺序填入11×19的表格； (3)附件5的结果：将碎片序号按复原后顺序填入两个11×19的表格； (4)不能确定复原位置的碎片，可不填入上述表格，单独列表。

2016年全国研究生数学建模竞赛B题

2016年全国研究生数学建模竞赛B题 具有遗传性疾病和性状的遗传位点分析 人体的每条染色体携带一个DNA分子，人的遗传密码由人体中的DNA携带。DNA是由分别带有A,T,C,G四种碱基的脱氧核苷酸链接组成的双螺旋长链分子。在这条双螺旋的长链中，共有约30亿个碱基对，而基因则是DNA长链中有遗传效应的一些片段。在组成DNA 的数量浩瀚的碱基对（或对应的脱氧核苷酸）中，有一些特定位置的单个核苷酸经常发生变异引起DNA的多态性，我们称之为位点。染色体、基因和位点的结构关系见图1. 在DNA长链中，位点个数约为碱基对个数的1/1000。由于位点在DNA长链中出现频繁，多态性丰富，近年来成为人们研究DNA遗传信息的重要载体，被称为人类研究遗传学的第三类遗传标记。 大量研究表明，人体的许多表型性状差异以及对药物和疾病的易感性等都可能与某些位点相关联，或和包含有多个位点的基因相关联。因此，定位与性状或疾病相关联的位点在染色体或基因中的位置，能帮助研究人员了解性状和一些疾病的遗传机理，也能使人们对致病位点加以干预，防止一些遗传病的发生。 近年来，研究人员大都采用全基因组的方法来确定致病位点或致病基因，具体做法是：招募大量志愿者（样本），包括具有某种遗传病的人和健康的人，通常用1表示病人，0表示健康者。对每个样本，采用碱基(A,T,C,G)的编码方式来获取每个位点的信息(因为染色体具有双螺旋结构，所以用两个碱基的组合表示一个位点的信息）；如表1中，在位点rs100015位置，不同样本的编码都是T和C的组合，有三种不同编码方式TT,TC和CC。类似地其他的位点虽然碱基的组合不同，但也只有三种不同编码。研究人员可以通过对样本的健康状况和位点编码的对比分析来确定致病位点，从而发现遗传病或性状的遗传机理。 表1. 在对每个样本采集完全基因组信息后，一般有以下的数据信息 rs