2013年数学建模国赛普通论文new

1、2013高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛

2、规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： （论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。） 日期 年 月 日赛区评阅编号（由赛区组委

3、会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：碎纸片的拼接复原摘要为了将碎纸片拼接在一起，而又无需耗费大量时间和人力。本文利用计算机将图像数字化，用otsu算法确定出图像二值化的最佳阈值，从而便能找出只包含0和1的图像矩阵（0代表黑色，1代表白色）。在第一问中碎片数量小，我们很容易通过人力看出第一个碎片，提取其右边界数据，和别的碎片左边界数据一一比对，找出比对率最高的碎片。为了使结果更

4、加明显，更直接地找出匹配碎片，我们不比对数值为1的点，只找出值为0且一一对应的点的数量k，利用程序计算出两个边界中0值点的数量num，则配对比为：C=第二三问数量大，待比较的边界增多。为了解决图的最优划分问题，我们应用谱聚类算法的思想，对数据进行了处理，减少计算机比对时的运算量。首先将二值化后的图像矩阵提取边界像素点，依0,1点分布特征，将其长分为5块，宽为2块（即均为每36个像素点为一块）。通过软件算出各块0点个数，分别将每个边界的数据作为一个向量。找出相匹配的碎片，即找出相似度最大的向量。向量中带比对的数值均是一一对应的，因此，很自然地想到可以用距离来度量样本点间的相似程度。对于定量变量，

5、常用的是 Minkowski 距离来度量样本点间的相似程度：d（x，y）=当q=1时，我们得到绝对值距离：d（x，y）= 利用计算机通过枚举法比较绝对值距离的大小，便能轻松找出相互匹配的碎片了。关键字：otsu算法 二值化 matlab 枚举法 Minkowski 距离一 问题重述破碎文件的拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。传统的拼接复原工作由人工完成，准确率较高，但由于效率很低，特别是当碎片数量巨大时，人工拼接很难在短时间内完成任务。因此我们试图在计算机高度发展的基础上，开发碎纸片的自动拼接技术，以提高拼接复原效率。1、 建立碎纸片拼接复原模型和算法，并

6、针对附件1、附件2给出的中、英文各一页文件的碎片 （仅纵切)数据进行拼接复原。并写出人工干预的方式和时间点。2、 设计碎纸片拼接复原模型和算法，并针对附件3、附件4给出的中、英文各一页文件的碎片（纵切和横切）数据进行拼接复原。并写出人工干预的方式和时间点。3、 尝试设计相应的碎纸片拼接复原模型与算法，然后将附件5给出的一页英文印刷文字双面打印文件的碎片数据进行复原。并写出人工干预的方式和时间点。二、问题分析2.1问题一的分析问题一要求我们建立拼接复原模型和算法，将其附件1、2给出的中、英文各一页文件的碎纸片（仅纵切）数据进行拼接复原。碎片切痕规则，且给出的每幅图像素均为198072，因此假设两

7、个碎片匹配，则当其对齐拼接时，切断的字符能完全正确地被连接。因此我们考虑从切割处入手。观察可得共享一个切割边的碎片，当像素点无限接近于切割处时，两者的像素点是对称的。因此我们利用计算机软件先将图片进行二值化处理，再选取最边上的一列像素点进行提取比对，利用各种方法进行对比，找出使比对结果最明显的方法。再利用其找出各匹配的碎片。2.2问题二的分析此问中横切与纵切的碎片数量多，面积小和需要考虑的边界数量多。若延用问题一的方法，直接比对像素点，则循环过后会导致预算量太大，matlab无法计算。因此我们将图片二值化处理后，再利用噪声处理的思想，将四个边界上个的像素点分块，算出各块黑色点个数，将各图的数据

8、储存在矩阵中，则只用比对各矩阵的数据，利用模糊聚类的方法来比较，配对的碎片。2.3问题三的分析问题三在问题二的基础上加上了正反面，加大了运算量，但我们问题二的解法已经尽量对数据进行了处理简化，此问依旧能延用，只需改动数据便能达到解题的目的。三、模型假设1、该文字文件被碎纸机破碎后，没有受到其它因素影响而有其它的损坏。2、该文字文件在被破碎之前文字数据正常，没有缺印、重印或字迹模糊、数据错位等现象。3、利用的电脑软件计算误差小，可以忽略不计。 四、符号说明dMinkowski 距离c改进后的陪对比c改进前的配对比num边界中0值点的数量n前景的点数csum前景的质量矩阵m前景的平均灰度n当前阈值

9、下的背景的点数sum-csum当前阈值下的质量矩m当前阈值下的平均灰度t前景与背景的分割阈值w前景点数占图像比例u平均灰度w点数占图像比例u。平均灰度u图像的总平均灰度g前景和背景图象的方差k位置与值均相同的点的数量p二值化处理后的灰度矩阵元素五、模型的建立及求解5.1 问题一的模型建立及求解首先将图片输入matlab软件，便自动将图片处理为数据矩阵，此时便可取出边界点进行匹配，通过程序，算出位置与值均相同的点的数量num，我们能得出直接匹配率c：c=100%例如：匹配附件1中的“008.bmp”与“001.bmp”运行程序（程序见附录）得：c = 0.7712num= 1527再用枚举法匹配

10、出合适的碎片时，数据如下：num= 1862c = 0.9404我们发现是否匹配所反映出的比率差异并不大，不难发现，空白处均被识别为匹配，导致匹配率都很高，因此我们需要进行改进，然后我们发现当灰度为40的点和灰度为50的点匹配时，肉眼看来没有区别，而计算机却判别为不同，此时必然也会参数误差。要解决这个问题，我们需要将图像数据处理为只有0和1（0代表黑色，1代表白色），至于判定是0还是1的准则，通过确定图像二值化最佳阈值即能实现。阈值将原图象分成前景，背景两个图象。前景：用n,csum,m来表示在当前阈值下的前景的点数，质量矩，平均灰度后景：用n, sum-csum,m来表示在当前阈值下的背景的

11、点数，质量矩，平均灰度当取最佳阈值时，背景应该与前景差别最大，关键在于如何选择衡量差别的标准，在otsu算法中这个衡量差别的标准就是最大类间方差。记t为前景与背景的分割阈值，前景点数占图像比例为w， 平均灰度为u；点数占图像比例为w，平均灰度为u。则图像的总平均灰度为：u= wu+ w u。前景和背景图象的方差：g= w (u)(u-u)+ w(u-u) (u-u)=w0w(u- u) (u- u)我们得到阈值后便能实现二值化了：p=得出二值化处理后的灰度矩阵：p=提取灰度矩阵左右边界构成新的匹配向量q:q=(q，qq)或q=(q，qq)接下来再来对二值化后的图片矩阵进行比对处理。改写后的程序

12、不再比对值为1的点了，此时利用以下公式计算：C=此时我们再来看与合适碎片与不合适碎片的比对率：不合适：K= 35C= 0.14056合适：K= 176C= 0.90683由此看出差距明显，实现起来误差小，具有较高可行性。通过编程实现循环查找，便能快速找出配对碎片了（代码见附录）。得出的结论和人工干预方式和时间见附录。5.2问题二、三的模型建立及求解首先处理数据，简化比对运算。二值化处理方法同上，此问中图像像素为：18072。我们先从宏观观察的左右边界一般有3个字符密集处，2个空白处；上下1个密集处，1个空白处。由此应左右边界分为5块，上下2块。再从微观，即像素点的分配考虑，发现如此刚好每块有3

13、6个像素点。此分法合理。然后再由计算机算出每块0（即黑点）的个数，做好编号，依次储存在向量中。由此，每个图像的数据我们均有记录之后，左右上下的边界都能通过比对找到匹配碎片了。对于向量的匹配，我们通过计算出的绝对值距离来判断是否匹配：d（x，y）= 如我们在matlab软件中输入附件4的图片“000.bmp”和“001.bmp”对000的上下左右边界处理后得出的向量a，b，c，d：a= 0 0b= 0 0c= 0 5 0 32 0d = 0 2 0 33 0对001的上下左右边界处理后得出的向量a，b，c，d：a= 0 0b= 32 4 c= 0 3 32 22 5d = 0 23 30 9 0

14、得出无论怎么摆放，均无法匹配。通过程序枚举算出000与102左右匹配。 102的数据：a= 0 0b= 0 0c= 0 10 0 0 0d = 0 5 0 32 0102的左边与000的右边绝对距离为：d（x，y）= =（0-0）+（5-5）+（0-0）+（32-32）+（0-0）=0由此例我们能看出匹配的两幅图的绝对值距离可以非常的小，我们通过多次试验发现此法可行性很高。人工调入图片数据后，再整合检查，便能得出匹配后的完整图像了。结果见附录。六、模型的评价及推广通过建立数学模型，以及调试完善的程序，对于碎纸机纵向剪碎纸片，我们能很好的迅速找出相匹配的图片，并能通过求匹配率的程序进行检验。所以

15、可以推广到各种纵向、横向剪切的碎纸片，运用此方法，均能复原，并能查看复原图。而既被纵切又被横切的碎纸片，只需少量人工工作，便能将十分零碎的碎纸片复原。对于我们模型的评价，优点是我们运用MATLAB将图片二值化，处理灰度，准确的找出特征点和与之匹配的点。当碎纸片较少时，我们在适时进行人工干预，输入碎纸片的编号，逐一匹配，这样匹配率较高，并且运用图像拼接的程序，观察匹配图，这样就可以及时找出错误；并且当碎纸片较多时，为提高准确率，我们运用了图像导入的程序，交给软件找出相匹配的碎纸片。缺点是在MATLAB中导入图片时，不能太多，增加了人工工作量。八、参考文献【1】 杨丹，赵海滨，MATLAB从入门到

16、精通：中国铁道出版社（第一版），2013年2月1日。【2】 周开利，邓春晖，MATLAB基础及其应用教程：北京大学出版社，2007年3月1日。【3】 张强，精通MATLAB图像处理：电子工业出版社，2013年1月。【4】 李明哲，金俊，石端银，图论及其算法：机械工业出版社，2010年10月。【5】 王海英，黄强，李传涛，褚宝增，图论算法及其MATLAB实现，北京市：北京航空航天大学出版社，2010年1月。【6】 邓敏，刘启亮，李光强，黄健柏，空间聚类分析及应用：科学出版社，2011年9月。附录一、两张图像匹配度的程序：clearclose allI=imread(099.bmp);figure

17、(1)imshow(I)J=find(I=150);I(J)=255;title( 图像二值化( 域值为150 );figure(2)imshow(I)clc;aw=im2bw(I,0.5);%选取阈值为0.5figure(3);imshow(aw) %显示二值图象k=imread(128.bmp);figure(1)imshow(k)J=find(k=150);I(J)=255;title( 图像二值化( 域值为150 );figure(2)imshow(k)clc;bw=im2bw(I,0.5);%选取阈值为0.5figure(3);imshow(bw) %显示二值图象I1=aw(:,72

18、);imshow(I1);k1=bw(:,1);imshow(k1);sum=0;for i=1:1980if I1(i)=k1(i)sum=sum+1;endc=sum/180end二、图像匹配、二值化程序：clearclose allfilename,pathname = uigetfile(*.bmp, Select all the bmp-files , MultiSelect, on);I=imread(filename1);figure(1)imshow(I)J=find(I=150);I(J)=255;title( 图像二值化( 域值为150 );figure(2)imshow(

19、I)clc;aw=im2bw(I,0.5);%选取阈值为0.5figure(3);imshow(aw) %显示二值图象max=0;for p=1:110k=imread(filenamep);figure(1)imshow(k)J=find(k=150);I(J)=255;title( 图像二值化( 域值为150 );figure(2)imshow(k)clc;bw=im2bw(I,0.5);%选取阈值为0.5figure(3);imshow(bw) %显示二值图象I1=aw(:,72);imshow(I1);k1=bw(:,1);imshow(k1);sum=0;for i=1:180if

20、I1(i)=k1(i)sum=sum+1;endendnum=0;for j=1:180;if (I1(j)=1)num=num+1;endendc=sum/180;if cmaxq=p;max=c;endend三、图像组合程序：附件一：I1=imread(008.bmp);I2=imread(014.bmp);I3=imread(012.bmp);I4=imread(015.bmp);I5=imread(003.bmp);I6=imread(010.bmp);I7=imread(002.bmp);I8=imread(016.bmp);I9=imread(001.bmp);I10=imread

21、(004.bmp);I11=imread(005.bmp);I12=imread(009.bmp);I13=imread(013.bmp);I14=imread(018.bmp);I15=imread(011.bmp);I16=imread(007.bmp);I17=imread(017.bmp);I18=imread(000.bmp);I19=imread(006.bmp);I=I1,I2,I3,I4,I5,I6,I7,I8,I9,I10,I11,I12,I13,I14,I15,I16,I17,I18,I19;imshow(I)附件二：I1=imread(003.bmp);I2=imread

22、(006.bmp);I3=imread(002.bmp);I4=imread(007.bmp);I5=imread(015.bmp);I6=imread(018.bmp);I7=imread(011.bmp);I8=imread(000.bmp);I9=imread(005.bmp);I10=imread(001.bmp);I11=imread(009.bmp);I12=imread(013.bmp);I13=imread(010.bmp);I14=imread(008.bmp);I15=imread(012.bmp);I16=imread(014.bmp);I17=imread(017.bm

23、p);I18=imread(016.bmp);I19=imread(004.bmp);I=I1,I2,I3,I4,I5,I6,I7,I8,I9,I10,I11,I12,I13,I14,I15,I16,I17,I18,I19;imshow(I)附件三：I001=imread(049.bmp);I002=imread(054.bmp);I003=imread(065.bmp);I004=imread(143.bmp);I005=imread(186.bmp);I006=imread(002.bmp);I007=imread(057.bmp);I008=imread(192.bmp);I009=im

24、read(178.bmp);I010=imread(118.bmp);I011=imread(190.bmp);I012=imread(095.bmp);I013=imread(011.bmp);I014=imread(022.bmp);I015=imread(129.bmp);I016=imread(028.bmp);I017=imread(091.bmp);I018=imread(188.bmp);I019=imread(141.bmp);I020=imread(061.bmp);I021=imread(019.bmp);I022=imread(078.bmp);I023=imread(0

25、67.bmp);I024=imread(069.bmp);I025=imread(099.bmp);I026=imread(162.bmp);I027=imread(096.bmp);I028=imread(131.bmp);I029=imread(079.bmp);I030=imread(063.bmp);I031=imread(116.bmp);I032=imread(163.bmp);I033=imread(072.bmp);I034=imread(006.bmp);I035=imread(177.bmp);I036=imread(020.bmp);I037=imread(052.bmp

26、);I038=imread(036.bmp);I039=imread(168.bmp);I040=imread(100.bmp);I041=imread(076.bmp);I042=imread(062.bmp);I043=imread(142.bmp);I044=imread(030.bmp);I045=imread(041.bmp);I046=imread(023.bmp);I047=imread(147.bmp);I048=imread(191.bmp);I049=imread(050.bmp);I050=imread(179.bmp);I051=imread(120.bmp);I052

27、=imread(086.bmp);I053=imread(195.bmp);I054=imread(026.bmp);I055=imread(001.bmp);I056=imread(087.bmp);I057=imread(018.bmp);I058=imread(038.bmp);I059=imread(148.bmp);I060=imread(046.bmp);I061=imread(161.bmp);I062=imread(024.bmp);I063=imread(035.bmp);I064=imread(081.bmp);I065=imread(189.bmp);I066=imrea

28、d(122.bmp);I067=imread(103.bmp);I068=imread(130.bmp);I069=imread(193.bmp);I070=imread(088.bmp);I071=imread(167.bmp);I072=imread(025.bmp);I073=imread(008.bmp);I074=imread(009.bmp);I075=imread(105.bmp);I076=imread(074.bmp);I077=imread(014.bmp);I078=imread(128.bmp);I079=imread(003.bmp);I080=imread(159.

29、bmp);I081=imread(082.bmp);I082=imread(199.bmp);I083=imread(135.bmp);I084=imread(012.bmp);I085=imread(073.bmp);I086=imread(160.bmp);I087=imread(203.bmp);I088=imread(169.bmp);I089=imread(134.bmp);I090=imread(039.bmp);I091=imread(031.bmp);I092=imread(051.bmp);I093=imread(107.bmp);I094=imread(115.bmp);I

30、095=imread(176.bmp);I096=imread(094.bmp);I097=imread(034.bmp);I098=imread(084.bmp);I099=imread(183.bmp);I100=imread(090.bmp);I101=imread(047.bmp);I102=imread(121.bmp);I103=imread(042.bmp);I104=imread(124.bmp);I105=imread(144.bmp);I106=imread(077.bmp);I107=imread(112.bmp);I108=imread(149.bmp);I109=im

31、read(097.bmp);I110=imread(136.bmp);I111=imread(164.bmp);I112=imread(127.bmp);I113=imread(058.bmp);I114=imread(043.bmp);I115=imread(125.bmp);I116=imread(013.bmp);I117=imread(182.bmp);I118=imread(109.bmp);I119=imread(197.bmp);I120=imread(016.bmp);I121=imread(184.bmp);I122=imread(110.bmp);I123=imread(1

32、87.bmp);I124=imread(066.bmp);I125=imread(106.bmp);I126=imread(150.bmp);I127=imread(021.bmp);I128=imread(173.bmp);I129=imread(157.bmp);I130=imread(181.bmp);I131=imread(204.bmp);I132=imread(139.bmp);I133=imread(145.bmp);I134=imread(029.bmp);I135=imread(064.bmp);I136=imread(111.bmp);I137=imread(201.bmp

33、);I138=imread(005.bmp);I139=imread(092.bmp);I140=imread(180.bmp);I141=imread(048.bmp);I142=imread(037.bmp);I143=imread(075.bmp);I144=imread(055.bmp);I145=imread(044.bmp);I146=imread(206.bmp);I147=imread(010.bmp);I148=imread(104.bmp);I149=imread(098.bmp);I150=imread(172.bmp);I151=imread(171.bmp);I152

34、=imread(059.bmp);I153=imread(007.bmp);I154=imread(208.bmp);I155=imread(138.bmp);I156=imread(158.bmp);I157=imread(126.bmp);I158=imread(068.bmp);I159=imread(175.bmp);I160=imread(045.bmp);I161=imread(174.bmp);I162=imread(000.bmp);I163=imread(137.bmp);I164=imread(053.bmp);I165=imread(056.bmp);I166=imrea

35、d(093.bmp);I167=imread(153.bmp);I168=imread(070.bmp);I169=imread(166.bmp);I170=imread(032.bmp);I171=imread(196.bmp);I172=imread(071.bmp);I173=imread(156.bmp);I174=imread(083.bmp);I175=imread(132.bmp);I176=imread(200.bmp);I177=imread(017.bmp);I178=imread(080.bmp);I179=imread(033.bmp);I180=imread(202.

36、bmp);I181=imread(198.bmp);I182=imread(015.bmp);I183=imread(133.bmp);I184=imread(170.bmp);I185=imread(205.bmp);I186=imread(085.bmp);I187=imread(152.bmp);I188=imread(165.bmp);I189=imread(027.bmp);I190=imread(060.bmp);I191=imread(089.bmp);I192=imread(146.bmp);I193=imread(102.bmp);I194=imread(154.bmp);I195=imread(114.bmp);I196=imread(040.bmp);I197=imread(151.bmp);I198=imread(207.bmp);I199=imread(155.bmp);I200=imread(140.bmp);I201=imread(185.bmp);I202=imread(108.bmp);I203=imread(117.bmp);I204=imread(004.bmp);I205=imread(101.bmp);I206=imread(113.bmp);I207=imread(194.bmp);I20