comp182dan.m

数学建模竞赛：DNA序列的分析

c=[.241748438893845 .257805530776093 .280107047279215 .220338983050847 
.265832681782643 .28068803752932 .220484753713839 .232994526974199 
.266900790166813 .269534679543459 .2071992976295 .256365232660228 
.20655487804878 .254573170731707 .270579268292683 .268292682926829 
.17032967032967 .154761904761905 .332417582417582 .342490842490843 
.27332361516035 .240524781341108 .241253644314869 .244897959183673 
.342105263157895 .320175438596491 .133333333333333 .204385964912281 
.183127572016461 .184499314128944 .338134430727023 .294238683127572 
.273840769903762 .241469816272966 .220472440944882 .26421697287839 
.350649350649351 .296969696969697 .165367965367965 .187012987012987 
.206033376123235 .274711168164313 .25609756097561 .263157894736842 
.363157894736842 .279824561403509 .156140350877193 .200877192982456 
.153577661431065 .136998254799302 .381326352530541 .328097731239092 
.2 .222881355932203 .26271186440678 .314406779661017 
.205128205128205 .129807692307692 .342948717948718 .322115384615385 
.13203125 .115625 .4015625 .35078125 
.284920634920635 .251587301587302 .208730158730159 .254761904761905 
.193573068094874 .208110175975516 .299923488905891 .298393267023718 
.213053263315829 .22430607651913 .27006751687922 .292573143285821 
.229771463866584 .209388511426807 .258801729462631 .302038295243978 
.232592592592593 .251851851851852 .232592592592593 .282962962962963 
.281517747858017 .323133414932681 .163402692778458 .231946144430845 
.306317044100119 .32181168057211 .16865315852205 .203218116805721 
.289178001182732 .301596688350089 .179775280898876 .229450029568303 
.263343108504399 .270381231671554 .201759530791789 .264516129032258 
.282952548330404 .304042179261863 .196250732278852 .216754540128881 
.163146394756009 .124544792425346 .399854333576111 .312454479242535 
.3068669527897 .284692417739628 .188841201716738 .219599427753934 
.19988577955454 .274700171330668 .263278126784694 .262135922330097 
.279545454545455 .302272727272727 .210795454545455 .207386363636364 
.216506500847937 .233465234595817 .248728094968909 .301300169587337 
.154545454545455 .137762237762238 .396503496503496 .311188811188811 
.24759751271905 .231769361221029 .233465234595817 .287167891464104 
.366108786610879 .283124128312413 .152719665271967 .198047419804742 
.21682665160926 .232072275550536 .250141163184641 .300959909655562 
.28088336783989 .243616287094548 .209799861973775 .265700483091787 
.257014590347924 .254769921436588 .261503928170595 .226711560044893 
.263945578231293 .254421768707483 .225850340136054 .25578231292517 
.220872446162341 .232468249585864 .252346769740475 .29431253451132 
.339325842696629 .258426966292135 .21498127340824 .187265917602996 
.219834710743802 .235261707988981 .249035812672176 .295867768595041 
.227272727272727 .221256684491979 .250668449197861 .300802139037433 
.372496662216288 .367156208277704 .113484646194927 .146862483311081 
.216171617161716 .240924092409241 .245324532453245 .297579757975798 
.184 .233333333333333 .282666666666667 .3 
.298750678978816 .278109722976643 .19663226507333 .226507332971211 
.277358490566038 .251572327044025 .219496855345912 .251572327044025 
.258264462809917 .300103305785124 .238119834710744 .203512396694215 
.172350791717418 .15468940316687 .343483556638246 .329476248477466 
.308372641509434 .295400943396226 .183372641509434 .212853773584906 
.314035087719298 .246198830409357 .213450292397661 .226315789473684 
.236658932714617 .240139211136891 .250580046403712 .27262180974478 
.28137310073157 .25267304445695 .205965109735509 .259988745075971 
.276671771311894 .312404287901991 .194997447677386 .215926493108729 
.18910741301059 .225416036308623 .275844679778114 .309631870902673 
.25 .28369140625 .22607421875 .240234375 
.263996138996139 .285714285714286 .211389961389961 .238899613899614 
.243888888888889 .200555555555556 .275555555555556 .28 
.173686701372456 .173686701372456 .345480359678183 .307146237576905 
.269631031220435 .236991485335856 .225638599810785 .267738883632923 
.186674008810573 .177312775330396 .325991189427313 .310022026431718 
.223880597014925 .202809482001756 .298507462686567 .274802458296752 
.281284004352557 .273667029379761 .216539717083787 .228509249183896 
.189371573175875 .232813159004639 .28469000421763 .293125263601856 
.388559093362115 .294657312466271 .116567728008635 .200215866162979 
.255389718076285 .254145936981758 .242951907131012 .247512437810945 
.185414091470952 .159456118665019 .325504738360115 .329625051503914 
.204481792717087 .211134453781513 .284663865546218 .299719887955182 
.167444963308873 .174449633088726 .343895930620414 .314209472981988 
.14 .149677419354839 .367741935483871 .34258064516129 
.217668972857903 .174028738690793 .325705162320383 .282064928153273 
.260126582278481 .252848101265823 .261075949367089 .225949367088608 
.200954400848356 .183457051961824 .316012725344645 .299575821845175 
.423498694516971 .254830287206266 .145691906005222 .17597911227154 
.377310061601643 .343429158110883 .124743326488706 .154517453798768 
.247717972930438 .338684293358514 .202706956248033 .210890777463015 
.286511156186613 .248985801217039 .201825557809331 .262677484787018 
.24510966212211 .25489033787789 .259336099585062 .240663900414938 
.188622754491018 .164670658682635 .335828343313373 .310878243512974 
.347222222222222 .287202380952381 .179563492063492 .186011904761905 
.177586206896552 .254022988505747 .310919540229885 .257471264367816 
.194603903559127 .219575200918485 .29735935706085 .288461538461538 
.322124756335283 .272904483430799 .181286549707602 .223684210526316 
.30952380952381 .314761904761905 .151428571428571 .224285714285714 
.296545914069082 .337545633249087 .151923616961528 .213984835720303 
.29331450094162 .322975517890772 .175141242937853 .208568738229755 
.319236016371078 .220100045475216 .23556161891769 .225102319236016 
.278314310798946 .273924495171203 .21027216856892 .237489025460931 
.195964125560538 .186995515695067 .283408071748879 .333632286995516 
.155622766015947 .185317569425351 .354687929612318 .304371734946384 
.182880434782609 .18179347826087 .314945652173913 .320380434782609 
.290023815824292 .349298756284731 .181264884890183 .179412543000794 
.303474257011302 .217664294683968 .22059439095856 .25826705734617 
.300487554529125 .262509622786759 .232486528098537 .204516294585579 
.242439144332432 .252766166707647 .253503811163019 .251290877796902 
.288512323943662 .235255281690141 .25 .226232394366197 
.331037277147488 .267828200972447 .177471636952998 .223662884927066 
.336965294592413 .296610169491525 .137610976594027 .228813559322034 
.325381306218224 .257723895189675 .173249902229175 .243644896362925 
.197714285714286 .192380952380952 .303238095238095 .306666666666667 
.235141560406311 .239464015560839 .244867084503998 .279230602982494 
.320872274143302 .283489096573209 .212227414330218 .183411214953271 
.291166380789022 .315608919382504 .18417667238422 .209048027444254 
.149962034927866 .140091116173121 .37243735763098 .337509491268033 
.237944162436548 .243231810490694 .259729272419628 .256979695431472 
.221008107155446 .23863235812478 .235107507930913 .305252026788861 
.282777428151878 .33606041535557 .179777637927418 .201384518565135 
.246646646646647 .248448448448448 .259459459459459 .245445445445445 
.147523709167545 .17351598173516 .332982086406744 .345978222690551 
.36744966442953 .304026845637584 .156375838926174 .172147651006711 
.19577824170726 .235908141962422 .299466481094874 .268847135235444 
.142068965517241 .141379310344828 .359310344827586 .357241379310345 
.24811320754717 .273396226415094 .211509433962264 .266981132075472 
.419063545150502 .297658862876254 .110702341137124 .17257525083612 
.184827073261435 .207326143547787 .305317962067683 .302528821123094 
.279499011206328 .295319709953856 .206328279499011 .218852999340804 
.151 .144 .339666666666667 .365333333333333 
.172644927536232 .192028985507246 .314311594202899 .321014492753623 
.276493545481837 .332933053137196 .221855298709096 .16871810267187 
.211875 .1646875 .3275 .2959375 
.29728982300885 .288440265486726 .195796460176991 .218473451327434 
.296564195298373 .331464737793852 .173417721518987 .198553345388788 
.265427643450018 .312342114760014 .200288704438831 .221941537351137 
.207160725453408 .197310819262039 .305972482801751 .289555972482802 
.320067869813358 .2776492364646 .195434212555915 .206848681166127 
.268901132537496 .276553412917049 .225283134374043 .229262320171411 
.331598125976054 .298542425819885 .152785007808433 .217074440395627 
.385416666666667 .307738095238095 .116964285714286 .189880952380952 
.329233922623623 .292083013066872 .168332052267487 .210351012042019 
.332564693825263 .289008455034589 .161926723033564 .216500128106585 
.318154311649017 .302571860816944 .200756429652042 .178366111951588 
.21213026590977 .238273080370481 .273976695548252 .275619958171497 
.244734931009441 .26957153231663 .232244008714597 .253449527959332 
.204913678618858 .216467463479416 .294289508632138 .284329349269588 
.241325127667932 .253109859892628 .253895508707608 .251669503731832 
.200402313301484 .177772190093035 .314307266784008 .307518229821473 
.193037974683544 .289435248296008 .256572541382668 .26095423563778 
.361388611388611 .296453546453546 .151348651348651 .190809190809191 
.227329347284504 .241155954160438 .266691579471849 .264823119083209 
.294330518697226 .288781664656212 .199517490952955 .217370325693607 
.195104562737643 .182271863117871 .33745247148289 .285171102661597 
.341290475064996 .291656818718979 .16048215551879 .206570550697235 
.335288596902863 .302205537306429 .151102768653214 .211403097137494 
.312419373753958 .323208631406122 .167937140846722 .196434853993198 
.243127743127743 .242434742434742 .256987756987757 .257449757449757 
.330045871559633 .29197247706422 .169495412844037 .20848623853211 
.30042194092827 .238396624472574 .228481012658228 .232489451476793 
.219447287615148 .233162743091095 .277584442169908 .269805527123849 
.254760897164621 .281950909860347 .208738891239949 .254549301735083 
.330874820733456 .224134398688793 .225568531038722 .219422249539029 
.384115029099623 .279869907565902 .154228004108182 .181787059226292 
.308529945553539 .23956442831216 .231479953803003 .220425672331298 
.264883155334673 .256529554827112 .233477847097388 .245109442740827 
.255074219933354 .208724628900333 .267191760072705 .269009391093608 
.178233193603478 .166278528178854 .303058531283962 .352429746933706 
.27423096318709 .324457892082703 .177407967725668 .223903177004539 
.246254588748884 .255680126996726 .259152693719615 .238912590534775 
.146286990508096 .171040387120789 .341522426949563 .341150195421552 
.272520518144962 .323939483832691 .1981607831504 .205379214871947 
.283414059431863 .27529781384998 .220709516952481 .220578609765676 
.32694911975237 .331785645192494 .172954149738828 .168311085316309 
.292726587525909 .321179574147353 .197569248162804 .188524590163934 
.287150095522884 .346789600465155 .189301436996428 .176758867015533 
.261509221455723 .277066028438688 .234971138955371 .226453611150218 
.250998934469899 .280367607884923 .235082578582845 .233550879062334 
.324059191930389 .292911378373885 .198913706146428 .184115723549299 
.26317761010328 .262160860491251 .24268207845026 .231979450955209 
.303043155550843 .295037641819531 .179196267628035 .22272293500159 
.251426586612415 .275790204932268 .237929836748871 .228998164045055 
.321848818601148 .327685211333898 .181869528381813 .16859644168314 
.178333074654664 .166537327331988 .30296445755083 .352165140462517 
.236490272996558 .232287246817709 .271635577615883 .25958690256985 
.286812536960378 .298639858072147 .182732111176818 .231815493790656 
.289752650176678 .29093050647821 .187868080094229 .231448763250883 
.293113596233078 .294879340788699 .182460270747499 .229546792230724 
.192566148677026 .220915581688366 .294834103317934 .291684166316674 
.262682161962738 .261206419479801 .23464305478694 .241468363770522 
.277489626556017 .272475795297372 .227394536652835 .222640041493776 
.296981926653799 .290665028952448 .19415686962625 .218196174767503 
.294440093970243 .290872705124859 .195858348559993 .218828852344906 
.3064344589205 .300008862891075 .179739431002393 .213817247186032 
.320910818767935 .330451309954161 .168151157157232 .180486714120672]
b=[.354545 .5 .045455 .1 
.327273 .5 .027273 .145455 
.254545 .518182 .1 .127273 
.3 .5 .081818 .118182 
.290909 .645455 0 .063636 
.363636 .463636 .081818 .090909 
.354545 .263636 .245455 .136364 
.290909 .5 .118182 .090909 
.218182 .563636 .145455 .072727 
.2 .563636 .172727 .063636 
]
a=[.297297 .135135 .171171 .396396 
.27027 .153153 .162162 .414414 
.27027 .063063 .216216 .45045 
.423423 .288288 .108108 .18018 
.234234 .108108 .234234 .423423 
.351351 .126126 .126126 .396396 
.351351 .189189 .099099 .36036 
.279279 .189189 .162162 .369369 
.207207 .153153 .207207 .432432 
.181818 .136364 .272727 .409091 
]
juna=mean(a);
junb=mean(b);
for m=1:182
sb=c(m,:)-junb;
sa=c(m,:)-juna;
   if (sum(sb.*sb)/4)^(1/2)>(sum(sa.*sa)/4)^(1/2)
       
         out=m,  'A'
   else
         out=m,'B'
       
   end
end