数学建模期末作业
《数学建模》期末考试试卷四与参考答案

《数学建模》期末考试试卷 班级 姓名 学号一、(15分)某厂利用甲、乙、丙三种原料生产A 、B 、C 、D 、E 五种产品,单位产品(万件)对原材料的消耗(吨)、原材料的限量(吨)以及单位问五种产品各生产多少才能使总利润达到最大? (1)建立线性规划问题数学模型。
(2)写出用LINGO 软件求解的程序。
二、(15分)用单纯形方法求如下线性规划问题的最优解。
123123123123max 614134248..2460,,0S x x x x x x s t x x x x x x =++++≤⎧⎪++≤⎨⎪≥⎩三、(15分)某厂生产甲、乙、丙三种产品,消耗两种主要原材料A 与B 。
每单位产品生产过程中需要消耗两种资源A 与B 的数量、可供使用的原材料数量以及单位产品利润如下表:设生产甲、乙、丙产品的数量分别为123,,x x x 单位,可以建立线性规划问题的数学模型:123123123123max 4003005006030504500..3040503000,,0S x x x x x x s t x x x x x x =++++≤⎧⎪++≤⎨⎪≥⎩利用LINGO10.0软件进行求解,得求解结果如下:Objective value: 35000.00 Total solver iterations: 2 Variable Value Reduced CostX1 50.00000 0.000000 X2 0.000000 66.66667 X3 30.00000 0.000000 Row Slack or Surplus Dual Price 1 35000.00 1.000000 2 0.000000 3.333333 3 0.000000 6.666667(1)指出问题的最优解并给出原应用问题的答案;(2)写出该线性规划问题的对偶线性规划问题,并指出对偶问题的最优解;(3)灵敏度分析结果如下:Objective Coefficient RangesCurrent Allowable Allowable Variable Coefficient Increase DecreaseX1 400.0000 200.0000 100.0000X2 300.0000 66.66667 INFINITYX3 500.0000 166.6667 66.66667Righthand Side RangesRow Current Allowable AllowableRHS Increase Decrease2 4500.000 1500.000 1500.0003 3000.000 1500.000 750.0000对灵敏度分析结果进行分析四、(10分)一个公司要分派4个推销员去4个地区推销某种产品,4个推销员在各个地区推销这种产品的预期利润(万元)如下表。
数学建模期末考核题目.doc

数学建模期末考核题考题一求出Y对X的回归直线方程,并说明拟合效果的好坏。
(请使用Matlab求解,并附上代码及图形)2据观察,个子高的人一般腿都长,今从16名成年女子测得数据如下表,希望从中得到身高x与腿长y之间的回归关系。
(请使用Matlab求解,并附上代码及图形)身高x与腿长),观测数据3、某人每天由饮食获取10467焦热量,其中5038焦用于新陈代谢,此外每公斤体重需支付69焦热量作为运动消耗,其余热量则转化为脂肪,已知以脂肪形式贮存的热量利用率为100%,每公斤脂肪含热量41868焦,问此人的体重如何随时间而变化?4、在一个巴基斯坦洞穴里,发现了具有古代尼安德特人特征的人骨碎片,科学家们把它们带到实验室,作碳14年代测定。
分析表明C14与C12的比例仅仅是活组织内的6. 24%,此人生活在多少年前?(宇宙射线在大气中能够产生放射性碳一14,并能与氧结合成二氧化碳形后进入所有活组织,先为植物吸收,后为动物纳入.只要植物或动物生存着,它们就会持续不断地吸收碳一14, 在机体内保持一定的水平,这意味着在活体中,C14的数量与稳定的C12的数量成定比。
生物体死亡后,交换过程就停止了,放射性碳便以每年八千分之一的速度减少•并逐渐消失. 对于任何含碳物质,只要测定剩下的放射性碳一14的含量,就可推断其年代.)5、你已经去过几家主要的摩托车商店,基本确定将从三种车型中选购一种。
你选择的标准主要有:价格、耗油量大小、舒适程度和外表美观情况。
经反复思考比较,构造了它们之间的成对比较矩阵13 7 81/3 1 5 51/7 1/5 1 31/8 1/5 1/3 1三种车型(记为a, b , 成对比较矩阵为c)关于价格、耗油量、舒适程度及你对它们表观喜欢程度的(价格)(耗油量)a h c a h ca_123_a_11/51/2_h1/212b517c_l/31/21c_21/71(舒适程度)(外表)a b c a b ca・]35~a・]1/53_b1/314b517c1/51/41c1/31/71(1)根据上述矩阵可以看出四项标准在你心目中的比重是不同的,请按由重到轻的顺序将它们排出。
数学建模期末考查作业

数学建模期末考查作业一、某化工厂生产A,B,C,D 四种化工产品,每种产品生产1吨消耗的工时,能该厂明年的总利润最高的数学模型,并利用MATLAB 写出简单的求解程序。
解:设该厂明年生产1A ,2A ,3A ,四种产品的数量分别为1x ,2x ,3x ,4x (单位:t ),总利润为z 。
约束条件 :工时限额:18480753802501004321≤+++x x x x能耗限额:1001.05.03.02.04321≤+++x x x x确定目标函数:4321852x x x x Z +++=4321852m ax x x x x Z +++=()⎪⎩⎪⎨⎧=∈≥≤+++≤+++4,3,2,1,01001.05.03.02.01848075380250100..43214321i N x x x x x x x x x x t s i i 且 求解:model:max=2*x1+5*x2+8*x3+x4;100*x1+250*x2+380*x3+75*x4<=18480; 0.2*x1+0.3*x2+0.5*x3+0.1*x4<=100; @gin(x1); @gin(x2); @gin(x3); @gin(x4); endGlobal optimal solution found.Objective value: 388.0000 Objective bound: 388.0000 Infeasibilities: 0.000000 Extended solver steps: 0 Total solver iterations: 0Variable Value Reduced Cost X1 2.000000 -2.000000 X2 0.000000 -5.000000 X3 48.00000 -8.000000 X4 0.000000 -1.000000Row Slack or Surplus Dual Price 1 388.0000 1.000000 2 40.00000 0.000000 3 75.60000 0.000000分析:由程序及结果可知,当四种化工产品生产数量分别为1x =2,2x =0,3x =48,4x =0时,该厂利润取最大值,最大值为388万元。
数学建模期末作业(嘉兴)

数学建模期末作业题一 河水污染问题如图是一个容量为32000m 的小湖,小河A 以310.1/m s -的速率向小湖注入河水,而小湖又以同样的速率通过小河B 流出,在上午9:20, 该地区发生交通事故,一个装有有毒化学物质的容器倾翻,在图中点X 处注入湖中,采取紧急措施后,于9:50分得到控制。
但数量不祥的有毒化学物质Z 流入湖中。
据估计Z 的量在35m 与320m 之间。
建立相应的模型来估计湖水受污染的程度随时间的变化函数关系并估计⑴湖水何时达到污染高峰;⑵何时污染可降至安全水平?(≤题二 选址问题考虑A,B,C 三地,每地都生产一定数量的原料,也消耗一定数量的产品(见下表),已知制成每吨产品需3吨原料,各地之间的距离为A —B :150,km B —C :200,km A —B :100.km 又每万吨原料运输1km 的运价是5000元,每万吨产品运输1km 的运价是6000元,由于地区条件的差异,在不同地点设厂的生产费用也不同,问怎样在何处设厂,规模多大,才能使总费用最小,由于其它条件限制,在B 处建厂的规模不能超过5万吨。
题三 雇员的聘用问题某服务部门一周中每天需要不同数目的雇员:周一到周四每天至少50人,周五和周日每天至少需要80人,周六至少需要90人,现规定应聘者每周需连续工作5天,试确定聘用方案,即周一到周日每天需要聘用多少人,使在满足需要的条件下,所聘用的总人数最少。
如果周日的需要量由80增至90人,方案应该怎样改变?若全时雇员(一天工作8小时)可以通过临时聘用的半时雇员(一天工作4小时,且无需连续工作)来代替,但规定半时雇员的工作量不得超过工作总量的四分之一,又设全时雇员和半时雇员每小时的酬金分别为5元和3元,试确定聘用方案,使在满足需要的前提下,所付的酬金为最小。
题四 肿瘤问题肿瘤大小V 生长的速率与V 的a 次方成正比,其中a 为形状参数,01;a ≤≤而其比例系数K 随时间减小,减小速率又与当时的K 值成正比,比例系数为环境系数.b 设某肿瘤参数1,0.1,a b ==K 的初始值为2,V 的初始值为1,问⑴此肿瘤生长不会超过多大?⑵过多长时间肿瘤大小翻一倍?⑶何时肿瘤生长速率由递增转为递减?⑷若参数2/3a =呢?题五 油气田的开发问题油气田开发试验表明:准确预测油气产量和可开采储量,对石油工作者来说,始终是一项既重要又困难的工作. 1995年,有人通过对国内外一些油气田的开发资料,得出结论:油气田的产量与累积产量之比()r t 与其开发时间存在着半指数关系:()lg .r t A Bt =-根据某气田1957~1976年总共20个年度的产气量数据(如下表),建立该气田的产量预测模型,并将预测与实际值进行比较.10m.注:产量单位83要求:每位学生在上面五题中可以任选一题,最迟于17周的周二前上交作业.。
《数学建模》期末试卷A

《数学建模》期末试卷A一、填空题(每题2分,共20分)1、在数学建模中,我们将所要研究的问题________化。
2、在解决实际问题时,我们常常需要收集大量的数据,这些数据通常是不________的。
3、在建立数学模型时,我们通常需要对变量进行假设,这些假设通常是对________的描述。
4、在解决实际问题时,我们通常需要对多个因素进行________,以确定哪些因素对所要研究的问题有显著影响。
5、在建立数学模型时,我们通常需要对数据进行________,以发现数据之间的规律和关系。
6、在解决实际问题时,我们通常需要将复杂的问题________化,以方便我们更好地理解和解决它们。
7、在建立数学模型时,我们通常需要将实际问题________化,以将其转化为数学问题。
8、在解决实际问题时,我们通常需要考虑实际情况的________性,以避免我们的解决方案过于理想化。
9、在建立数学模型时,我们通常需要使用数学语言来________模型,以方便我们更好地描述和解决它。
10、在解决实际问题时,我们通常需要使用计算机来帮助我们进行________和计算。
二、选择题(每题3分,共30分)11、在下列选项中,不属于数学建模步骤的是()。
A.确定变量和参数B.建立模型C.进行实验D.验证模型12、在下列选项中,不属于数学建模方法的是()。
A.归纳法B.演绎法C.类比法D.反证法13、在下列选项中,不属于数学建模应用领域的是()。
A.物理学B.工程学C.经济学D.政治学14、在下列选项中,不属于数学建模语言的是()。
A.文字语言B.符号语言C.图形语言D.自然语言15、在下列选项中,不属于数学建模原则的是()。
A.简洁性原则B.一致性原则C.可行性原则D.可重复性原则16、在下列选项中,不属于数学建模步骤的是()。
A.对数据进行分析和处理B.对模型进行假设和定义C.对模型进行检验和修正D.对结果进行解释和应用17、在下列选项中,不属于数学建模应用领域的是()。
数学建模期末大作业

数学建模承诺书
我们仔细阅读了数学建模作业的对应规则。
我们完全明白,在开始做题后不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反规则的。
如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守规则,以保证公正、公平性。
如有违反规则的行为,我们将受到严肃处理。
我们选择的题号是(A/B/C/D题): D
参赛队员:
1. 专业年级软件工程姓名段永春学号201410413112 成绩
2. 专业年级软件工程姓名殷福贵学号201410413113 成绩
3. 专业年级软件工程姓名高培富学号201410413107 成绩
日期: 2015 年 6 月 15 日。
数学建模期末作业题

数学建模期末作业题1、数学规划有三种物品:A、B和C。
它们的重量、体积和价值如下表所示:A、B和C重量(单位:公斤)体积(单位:l)123213价值(单位:100元)357当有人旅行时,选择10件物品陪伴他。
根据情况,个人物品的总重量不得超过18kg,体积不得超过100L。
在这三件物品中,你能选择多少件来最大限度地提高你物品的价值?2、谣言的传播假设一个城市有n+1个人。
其中一人出于某种目的编造了一个谣言,所以他利用他认识的人来传播谣言。
该市初中及以上文化程度的人口比例为p。
只有1%的人相信这个谣言,而其他人中约有B%会相信。
还假设单位时间内每个人相信谣言的平均人数与当时没有听到谣言的人数成正比,而不相信谣言的人不会传播谣言。
试图建立一个数学模型来反映谣言的传播,并简单分析其规律。
假设1第一个人仍将参与第二次谣言传播。
也就是说,第一个人和相信谣言的人会继续传播谣言。
假设2相信此谣言的人每人在单位时间内传播的平均人数正比于当时尚未听说此谣言的人数这个比恒定不变假设3在传播的同时,它也会传播给那些传播谣言和听到谣言的人设第i个单位时间开始时相信谣言总人数xyz(i)没有听说过MT(I)的人数受传播人数中没听过的人数占总人数比例(共有n+1个人,出去自己就有n个人)t(i)=mt(i)/n;受传播人数如果k为定植scb(i)=k*mt(i)*xyz(i);没有听到谣言的人数(考虑到他们也会传播给那些传播谣言和听到谣言的人)sch_mt(i)=scb(i)*t(i);其中相信的有scb_uumt_uxx(i)=sch_mt(i)*p*a/100+sch_mt(i)*(1-p)*b/100;有些人不相信scb_mt_bxx(i)=sch_mt(i)-scb_xx(i);在时间I+1的单位时间开始时相信谣言的总人数xyz(i+1)=xyz(i)+scb_mt_xx(i);没听过人数mt(i+1)=mt(i)-sch_umt(i);受传播人数中没听过的人数占总人数比例t(i+1)=mt(i+1)/n;如果K为定殖,则SCB(I+1)=K*MT(I+1)*XYZ(I+1);受传播人数中没听过谣言的人数(考虑到传播的时候也会传给传播谣和听过谣言的人)sch_mt(i+1)=scb(i+1)*t(i+1);其中包括scb_mt_xx(i+1)=sch_mt(i+1)*p*a/100+sch_mt(i+1)*(1-p)*b/100;其中不相信的有scb_umt_bxx(i+1)=sch_mt(i+1)-scb_xx(i+1);可以看到各种数构成了一个循环,这样就可以无限迭代下去根据由1单位时刻相信谣言总人数xyz(1)=1没听过人数mt(1)=n然后重复。
数学建模期末大作业

数学建模期末大作业论文题目:A题美好的一天组长:何曦(2014112739)组员:李颖(2014112747)张楚良(2014112740)班级:交通工程三班指导老师:陈崇双美好的一天摘要关键字:Dijkstra算法多目标规划有向赋权图 MATLAB SPSS1 问题的重述Hello!大家好,我是没头脑,住在西南宇宙大学巨偏远的新校区(节点22)。
明天我一个外地同学来找我玩,TA叫不高兴,是个镁铝\帅锅,期待ing。
我想陪TA在城里转转,当然是去些不怎么花钱的地方啦~~。
目前想到的有林湾步行街(节点76)、郫郫公园(节点91),大川博物院(节点72)。
交通嘛,只坐公交车好了,反正公交比较发达,你能想出来的路线都有车啊。
另外,进城顺便办两件事,去老校区财务处一趟(节点50),还要去新东方(节点34)找我们宿舍老三,他抽奖中了两张电影票,我要霸占过来明晚吃了饭跟TA一起看。
电影院嘛,TASHIWODE电影院(节点54)不错,比较便宜哈。
我攒了很久的钱,订了明晚开心面馆(节点63)的烛光晚餐,额哈哈,为了TA,破费一下也是可以的哈。
哦,对了,老三说了,他明天一整天都上课,只有中午休息的时候能接见我给我票。
我主要是想请教一下各位大神:1)明天我应该怎么安排路线才能够让花在坐车上的时间最少?2)考虑到可能堵车啊,TA比较没耐心啊,因为TA叫不高兴嘛。
尤其是堵车啊,等车啊,这种事,万一影响了气氛就悲剧了。
我感觉路口越密的地方越容易堵,如果考虑这个,又应该怎么安排路线呢?3)我们城比较挫啊,连地图也没有,Z老师搞地图测绘的,他有地图,跟他要他不给,只给了我一个破表格(见附件,一个文件有两页啊),说“你自己画吧”。
帮我画一张地图吧,最好能标明我们要去的那几个地方和比较省时的路线啊,拜托了~2 问题的分析2.1 对问题一的分析问题一要求安排路线使得坐车花费的时间最少。
对于问题一,假设公交车的速度维持不变,要使花费的时间最少,则将问题转化为对最短路径的求解。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
指纹识别技术研究摘要目前大部分研究认为指纹具有唯一性,可以作为个人身份验证的凭证,因此指纹识别技术一直以来就会研究热点,而随着移动设备对指纹识别的使用,指纹识别技术更需要进一步研究与改进。
本文主要研究了指纹识别从采集到匹配所采用的算法。
指纹识别的主要步骤是图像预处理、特征点提取、匹配。
在指纹预处理过程上,本文对预处理方法进行分析筛选,选择了最优的方案处理;在提取上,采用3*3矩阵遍历提取的算法,并且对于提取的特征点进行筛选,保证了特征点的准确度与代表性;在匹配部分,采取了使用距离与角度的匹配,保证了在指纹图像在旋转一定角度后的匹配鲁棒性。
关键词指纹识别;预处理;特征提取;匹配一、问题的重述指纹自古便被作为人身份的鉴别工具,而且随着指纹学的研究基本成熟,指纹识别越来越多的使用,但是靠人工对比指纹卡鉴别速度慢、效率低。
随着计算机技术的发展,人们开始研究使用计算机识别指纹的技术,以提高效率。
1.对于采集到的指纹图像进行预处理,将图形中的噪点,冗余信息剔除,并将图像转化成为二值图像以便于处理,最后要将图形中指纹纹线细化为单像素宽的轨迹,以便于提取特征点。
2.提取特征点,主要是提取出指纹中具有代表性的端点和分叉点作为特征点提取出来。
3.匹配,就是将已经提取到的两幅图的特征点集进行匹配。
二、模型假设1.不考虑提取过程中未提取到指纹中心点的的指纹图像。
2.假设指纹图像是没有巨大伤痕的。
3.假设指纹图像是基本清晰可是别的。
三、变量说明四、指纹图像预处理直接采集到的原始图像并不能满足提取特征点的要求,会存在很多的噪声。
这些噪声点由采集仪器,环境因素,人为操作等所产生的各类噪声所组成。
如采样和量化产生的高频散粒噪声、光照不均引起的低频噪声以及采集头上的污渍所引入的噪声;手指被弄脏,手指有刀伤、疤、痕、干燥、湿润或撕破等。
【1】这些噪声严重影响了指纹图像的特性,致使指纹图像包含的部分细节特征不清晰甚至出现丢失,同时引入许多虚假的特征信息,如果直接对原始指纹图像进行特征提取,势必会影响指纹识别结果的准确性。
所以,在提取特征之前需要对采集到的指纹图像进行预处理,以去除噪声信号,使其变成一幅纹线结构清晰特征信息明显的二值点线图。
主要流程如下图4-1所示:图4-1 图像预处理流程归一化首先需要将指纹图像与背景图分离开,所以首先进行归一化处理。
归一化利用公式如下4-1所示:【2】200200((,))I(,);(,)((,))(,);i ii iV I i j M M i j M V G i j V I i j M M I i j M V ⎧-+〉⎪⎪=⎨-⎪-〈⎪⎩若若 (4-1) 如果i i i M V M j i I V M 〈--200)),((,则把灰度值ii V M j i I V M 200)),((--归一化为255背景处理,其中0M 和0V 为期望的均值和方差,根据实际情况而定,i M 和i V 为指纹图像的均值和方差。
分割分割图像基于块特征的指纹图像分割,本文处理将采用均值方差法【3】。
该算法基于背景区灰度方差小,而指纹区方差大的思想,将指纹图像分成块,计算每一块的方差,如果该块的方差小于阈值为背景,否则为前景。
具体步骤分以下三步:(1)将图像分块,将其分为N ×N 的方块。
(2)计算出每一块的均值和方差。
原图归一化分割二值化 细化(3)如果方差很小,接近于0就认为是背景,对于方差不为零的区域在进行阈值分割算法。
二值化二值化是对分割后的图像的进一步处理,使图像灰度值只有两个值(黑、白)。
以便于进一步处理。
二值化即选择某一阈值,对于大于阈值的灰度置为1,小于阈值的灰度值都置为0,如下4-2公式:【4】()1,0,x Tf x x T ≥⎧=⎨<⎩公式4-2 二值化公式 但是由于原始指纹图像不同区域深浅不一,如对整幅图像用同一阈值进行二值分割,会造成大量有用信息的丢失。
所以我们可以选用局部阈对图像进行二值化。
局部阈值法即选取N ×N 的块,求该区域的阈值并对该区域二值化,可以有效地保证信息的可靠性。
二值化结果如下图4-3所示:图4-2 二值化细化二值化以后还需要将纹线宽度细化为单像素才便于提取特征点,本文细化采用MATLAB 里的bwmorph 函数对图像进行细化。
代码为:w=bwmorph(u,'thin',Inf);细化图如下图4-4所示:图4-3 细化图五、特征点的提取特征点提取概述特征点提取是指纹识别的关键,特征点指能代表指纹特异性的指纹纹路信息。
一般说来,这种特征应有以下性质:【5】(1)单一性:要求这种特征能够充分体现指纹的唯一性。
(2)可测试性:适用于指纹匹配算法,便于在匹配算法中应用。
(3)紧凑性:要求提取的特征不应包含指纹唯一性以外的冗余信息,并且信息量要尽量小,便于存储、管理和计算。
(4)鲁棒性:要求这种特征对噪声的存在与指纹形变不敏感。
为了比较两个指纹是否相同,需要从指纹图像中提取出能表示指纹唯一性的特征。
Galton提出的指纹细节点是人工指纹匹配中最常用的特征。
指纹由脊线和谷线交替构成,在大多数地方纹线连续且相互平行,而某些局部不连续的地方构成了细节点。
Galton定义了4种细节点类型:分叉点,端点,环、岛,并指出细节点具有唯一性,可以用于指纹匹配。
但是基本上的特征点都是由两类最基本的特征点构成的,即端点和分叉点。
本文提取的特征点即最基本的特征点,端点和分叉点【6】。
提取到特征点后并不可直接使用,因为由于指纹纹路的复杂性,往往产生很多特征点,其中很多都是伪特征点,需要我们剔除,所以提取到端点和分叉点以后需要去除为特征点。
端点提取概述特征点提取去采用3×3的块进行遍历提取,端点的判断条件为:周围的8邻域两两相邻当且仅当存在2个不同值。
如下图5-1所示:(1) (2)(3) (4)图 5-1 端点提取的四种情况分叉点提取概述特征点提取去采用3×3的块进行遍历提取,分叉点的判断条件为:周围的8邻域两两相邻当且仅当存在6个不同值。
如下图5-2所示:(1) (2)(3) (4)图 5-2 分叉点提取的四种情况伪特征点剔除因为采集、预处理过程中的空因素和提取算法的原因,提取到的特征点中会存在很多为特征点。
为了不影响后续的识别,必须进行剔除。
文本中剔除的伪特征点主要有三种:指纹范围外的点、指纹边沿的端点、断点,并根据这三种伪特征点的特征进行剔除。
在伪特征点剔除算法中,主要使用到了特征点距离的概念,特征点1p 和2p 间距离用欧式距离计算,见式(5-1):()12,D p p =(5-1)下面是三类伪特征点的筛选方法:第一类:指纹范围外的点,本文中范围用了17。
此时把它和周围8点记为0,并在上述基础上减去不是特征点的个数。
第二类:指纹边沿的端点,【7】观察细化图可以发现,真正的端点的在它周围半径()2R R λλ<<,λ为平均纹线宽度,本文中8λ=。
在R 范围内没有其它特征点存在,则是真的端点。
而周围有其他特征点的,则可以认为是指纹的边沿端点,即第二类特征点。
如下图所示,展示了局部的端点标记,在图中绿色标记的是真正的端点,而红色标记的则是指纹边沿的端点,即第二类伪特征点。
图5-3 第二类伪特征点对比第三类为断点:断点,【8】即被认定为端点且在12的上半范围或下半范围,具有被确定为端点的特征点的点,则认为是第三类为特征点,即端点。
如下图所示,展示了局部的端点标记,其中绿色的为较长纹线的端点,判断为真特征点;而很短的纹线(一般是有假信号)的端点则根据12的范围,被判断为断点,即第三类伪特征点。
图5-4 第三类伪特征点对比六、匹配指纹匹配算法综述对于指纹识别的算法人们已经有了很多研究,也提出了很多匹配的算法,目前来看主要有两类:一类是基于图形的匹配方式,包括点模式匹配和基于图论的方法;另一类是采用人工神经网络的方法。
图形匹配是针对纹线几何形状及其特征点拓扑结构的匹配方式,它的原理是基于相似变换的方法把两个特征点集中的相对应点匹配起来,这些相似变换可以是平移变换、旋转变换、伸缩变换等线性变换,可以在一定程度内允许少量伪特征点的存在、真正特征点的丢失以及轻微的特征点定位偏差,且对图像的平移和旋转也不敏感。
但这种方法有两点不足:一是匹配速度比较慢;二是对指纹图像的质量要求比较高,低质量的图像匹配效果不佳,下面对这些算法进行一些简要的介绍【8】。
Ranade和R.Seinfeld提出了点模式匹配的松弛算法,其思想是寻找一对匹配点,使得反映匹配程度的相似变换最大,则将该点对作为基准点对,然后根据相似变换的计算结果调整待识别图像的位置,统计最终的匹配点对数,给出匹配结果。
Stookinan等提出的基于Hough变换的方法,把点模式匹配转化成对转换参数的Hough空间中的峰值检测。
这种方法的缺点在于当特征点数目较少(少于30个)时,很难在Hough空间里积累起足够大的证据来保证一个可靠的匹配,另外,该方法有计算量较大的缺点。
Sparrows与A.K.Hrechak等都提出了基于结构信息的特征匹配方法,而D.K.ISenor与S.G.zaky使用图来表示指纹特征,并用图匹配的方法来匹配指纹图像。
这类方法利用了指纹图像的拓扑结构,允许一般的图像平移旋转、特征点丢失以及伪特征点的存在,但是这类方法的准确性在很大程度上依赖于所提取的指纹特征信息及分类信息的准确性。
采用人工神经网络的指纹匹配方法也有很多。
Vinod将非对称神经网络应用于指纹匹配中,提出了一种基于非对称神经网络的点模式匹配算法,而田捷等人将遗传算法应用于指纹匹配中,提出了基于遗传算法的指纹图匹配方法,利用指纹图像的结构信息进行初匹配,缩小搜索空间,然后采用遗传算法和补偿算法匹配指纹图像,有较强的抗噪声与非线性形变的能力。
但由于神经网络固有的反复处理特性,速度难以得到提高,计算量偏大,因此不适合用于对实时性要求较高的在线指纹识别系统【9】。
指纹中心点的提取因为中心点一般位置固定,以指纹中心点定位具有鲁棒性强,准确性高的有点,所以本文采用中心点为参考等进行匹配。
以中心点为参考进行匹配,首先需要寻找中心点。
本文采用了基于Sobel算子的指纹中心点定位,这里选取的中心是指纹中心的一个小区域,先求出指纹图像的点方向,相邻8个灰度值之和的平均值,再求这8个灰度值与平均值之差的和,最小和所在的方向即此点所在指纹脊线的方向,如此得到点方向图。
把点方向图分为若干块16×16大小的小块,对每块计算直方图,其峰值方向即为块方向,即每块中点的主导方向。
然后在这个粗的块方向图上按照以下原则去搜索中心区域,逐行检查块方向数组。
然后再根据求出各个方向的角度以及相邻8个灰度值之和的平均值,再求这8个灰度值与平均值之差的和,最小和所在的方向即此点所在指纹脊线的方向,如此得到点方向图。