neiji182.m

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

a=[4.58405668130771
4.4130562982672
3.22291941908037
5.62231348485062
5.05860484335573
3.30265363596077
2.47931040755233
4.78198847944689
3.00695184856242
2.38506787232984
6.52926556187061
2.31407284245668
5.22977455268592
4.34416791477895
4.55712883181377
6.29122843134499
3.0076724053085
4.81899494970704
4.16494187793214
4.12315333031638
4.52761114917267
3.32317572131424
2.90151268985269
3.19677059586219
3.65023355657929
3.19501084592912
5.45883623709917
3.39089874719517
5.76496390390767
3.71005219667447
4.03794305659366
5.964444254727
4.09984707488109
2.36130690334176
4.0748988817823
3.04126378504151
3.84895261788849
3.3705517477053
4.01425652991511
2.25753452669035
4.00333392699674
4.04963014177107
2.2064139822269
4.02926474987093
4.99538738797493
3.06299558464004
3.2840274671587
3.59977424021095
4.89929889407317
2.93630960406753
2.77078998943875
3.74504267612992
3.06336779701371
3.51768799672965
6.08255452670413
3.88990205720633
3.67462513910763
4.15763657917519
5.36346566870152
3.8546595129006
4.55002029334105
4.55688350292773
3.14614487760469
4.79547696823041
2.09443449902267
3.92113469092473
5.73196413301273
5.47913758197662
5.93557830546942
7.31816832489603
3.86824736401643
3.71956001490697
4.12235926056381
2.00833048150093
1.8529552725573
3.55676393408745
3.125814457437
4.10851075103693
4.50801000273662
2.5109113761781
8.64136434330229
7.32548113246024
2.66153357644156
2.61654299552428
3.11893264696723
3.02400804116344
2.93700590197382
3.20503461532571
4.64153156278605
6.1267291757938
5.90754396099808
3.18961275249932
3.266728816288
3.03308135729027
3.90081414007225
3.30797058100008
2.5476618727145
2.43897658472602
2.55957165797964
4.53417882560971
4.06722207842827
2.6435642106838
3.34662692028004
5.4635285317982
4.15176897257148
4.15459713044406
3.52586966353127
4.31262670731593
5.93247162214596
2.37441324599765
5.10628229181863
9.4480077662276
4.1473968330111
1.9375620873331
5.34179916554801
3.27830647219412
6.96261119081841
5.66725269375988
3.54422031070653
4.03730322473436
3.14679569525122
3.34155935420827
3.89655470974308
4.72077838275567
2.94649226718013
3.80700435108222
2.59049839660598
2.19075210632255
2.44929148042934
2.32752363830592
2.95710717213425
4.61504011124159
4.2771045489666
4.91445602943651
3.88398869791962
4.42007353219653
5.13391505541756
2.12353027945798
4.16417124849292
3.1396018774963
4.73405166041101
2.47974678742138
2.57468630275072
3.07826289598025
4.09945750852067
2.5924561366911
2.7601488528481
3.8956264102717
4.05053472482713
2.86081301399929
2.36207989645055
3.0299304852429
3.69767205254362
3.90364533605027
5.58088758580974
3.51650007496661
4.64275932432648
5.99863171711491
3.66609444352117
3.40604562307123
2.89733865518666
3.18432147191151
3.38246770797084
3.66966282072018
4.17883598854252
2.79651339038807
3.62914362730552
3.08490146896108
5.24618512232444
2.89520605069221
5.58073588096541
4.11247120283715
3.20286924563218
3.03422981384054
3.00929683860559
5.99836838801188
3.85208840880006
3.62823511139045
3.28603070229023
3.54546426933805
3.02834959343464
3.01107878246098]
t=[3.83965134795695
3.19154286932426
3.0810474381596
4.969582574808
6.15758450024628
3.9559673298717
2.64779721259236
4.42562654007704
3.88214567942917
2.93193202718188
3.42335907588233
3.10160973327135
5.84048716483314
3.91548873465833
6.63700570441949
7.32151011043962
3.81423156970704
3.71275959979076
3.95763300663516
4.40442809952914
4.08551526185859
2.65021845098976
2.88173288319143
2.96774136730766
3.40311930305298
3.01991540016345
7.00533335682937
3.62760662240064
4.67042546703675
3.12198203771221
3.83316951478064
6.55321572329684
5.47066143494438
3.54143234626827
3.85713394360764
4.2226531928457
3.63471710248393
3.4992493728139
3.87261285900389
3.52937629768313
3.80498735201109
4.43404585470354
2.21882959353207
3.68751520049374
3.33054011348779
3.36514272574533
3.69407646502605
3.08068305150721
5.50642273948421
3.29647349070941
3.38971858231334
3.69557721560976
4.16098766070749
3.03671822680508
5.23434258706272
3.07493028149889
3.22273338897735
4.6689740946224
4.99173956472361
3.98834172802169
4.38730911931348
4.79286503254792
3.37715526625945
4.04895223348926
2.73897608901392
4.2640718272654
7.07007091183235
4.49496019934086
5.19100890177864
5.88160595762797
4.83778583436071
3.55298555867153
5.01964794788719
3.69518222671018
2.09741952864311
2.68055325015315
3.85576521853639
3.87143844524227
5.22956220481497
3.48412857387695
3.73219836050727
6.59256070237459
3.4274608902349
2.72614835577257
2.65835759678912
2.86867333805206
5.29709917228992
3.52508483271485
4.41742274699534
4.84385558420997
5.9284843372258
2.69203432263878
4.23844871614142
3.37454289459986
3.72232342054852
3.91147940042796
3.56250347396002
3.15735663173534
3.69921839519755
4.63313271697194
3.86340013261711
3.26879577582907
3.01314517189466
5.86465238636034
3.94766690803492
3.82327669200106
2.77870811160853
4.12459551175543
4.47141130096654
2.85655784988138
4.11663463220692
6.61645631189652
3.52544445227313
3.11098202225215
4.57836756467811
3.03619366943222
6.91830600413641
4.80623718979158
2.74906570315972
5.44320601867352
3.41033861229819
2.71508003625862
2.86208847933307
5.14019875011809
3.55207876838144
3.79674334110451
2.91079632404922
2.86595858377275
3.13184618719751
3.21682844851303
3.14700444926393
3.7858826732662
3.52595834231766
4.39522220134475
3.81235350687063
4.76434311342846
2.91042130871551
3.03280637159502
3.78887830954723
3.25100001106846
5.12502121171406
3.07086758431126
2.87405443949017
2.89234495201885
4.58964781388819
3.07555649375687
4.03536220036336
3.91801253686341
3.64181402880286
4.89210813972318
3.88460817140733
3.79905441761821
3.83824324274797
4.74572160705314
5.92001013722348
2.72456814373099
4.34212178414811
5.1928294593768
2.81763879034074
3.370133197421
2.62362370409363
2.82196023504773
2.59779772660184
3.43858113922808
3.51508911791219
3.23927967587233
3.84930448641949
3.26086732665098
4.05587093290324
2.72498532687408
5.90525148208637
4.40204229658752
2.98149541408761
3.09862378345734
2.97174671819693
4.3032631211791
3.58791184254674
3.52092382602393
3.27599116977845
3.30910575857036
3.16940539557255
2.80375289117727
]
c=[3.44902551645956
4.37761324379531
4.46218273279591
3.76567215875059
1.8338075901723
3.8797006329114
6.30216333437628
2.87842726155699
3.7676247180334
5.82414002864274
4.07424171971592
4.75452910343154
1.70413024638881
2.77720303949234
1.91403868983024
1.85099768617103
4.17043700608662
2.58908304992567
2.87655214045237
4.04529277711046
3.81102351359639
5.73432346471291
5.45861443665953
5.54321091861614
5.00609997364283
4.98505435623316
1.64655348912791
5.34169566654341
3.68071915167738
4.62336773226231
3.83208145734564
1.67820256935377
4.25439476331991
5.87756282193798
3.75796012356858
3.74847100893351
3.55771316317484
3.8906338752097
3.8336103536923
3.65116283703822
3.81924690971644
3.56364317722539
8.11997075421517
3.81589554516637
2.71103554663389
4.96766862921895
3.63474802740121
4.15724083033424
2.10794006397541
4.7251683339886
4.08214178746555
3.51432922493856
3.82368054258796
5.46167751335087
3.81565147196733
4.397306609219
5.22913953124512
3.14005080028017
2.74573569413144
4.30311233445476
2.33144140634741
3.26691267497838
4.04833026575343
3.18671725516145
8.03489726167367
4.27036742384869
2.81365059544277
3.83393266531928
2.64604567954016
2.28438071612892
2.12499836254902
3.71114294497356
2.30465232664636
6.63625692550794
7.33565874117916
4.91595867655016
4.19297472815679
3.90334576513956
2.3177746337197
5.06229734574566
3.18996009546014
3.89872444873409
5.30849297049864
5.66994752317385
6.55137869234888
4.82639986377829
3.69640236727118
3.95343245391644
2.67122186595157
2.50293040287796
3.40523154874578
5.58479856382498
4.34716729376646
4.16913478295078
4.20463213353755
4.52973679265505
5.30977398256973
6.39905193719746
5.35787564716922
2.25444832908353
4.24382419552315
3.87437071396194
5.50191482481884
1.80123244099057
3.83528166924365
3.66062664108899
5.21430706771529
4.40564162978688
1.99072203184704
6.25875212193623
3.72189988778913
1.72904583441985
5.00026290321201
11.1138630346066
3.35498301483925
4.75242178478645
1.8896147354136
3.18680393589698
4.0156497147574
2.19078486692363
5.14132701213927
6.15102131781413
5.06571677328091
3.13605866296735
5.08111206008014
4.79172035795748
6.00369503017464
9.06694852535047
5.55818899499781
5.80159471569576
4.83564020035618
3.79547696727692
4.44890097923884
3.64037064392714
4.12177146272607
2.40154594223993
3.09981432993624
6.56448042589814
3.96805523281086
4.89223440595772
1.9581064039109
5.79151730187122
6.01631801198065
6.21345525752098
4.14184596923812
5.68987930585909
3.57817369630253
2.94215444703402
4.78619935949465
3.791921331406
6.42732356403207
3.93712581425376
4.55051262609209
3.7850047338474
2.35469857865337
5.72306794427777
3.98529711519788
1.95003173675796
5.10437506517866
4.28635450241684
6.10676851148266
4.93954783307925
5.51069959826319
4.58050008108831
4.88066151856979
5.16883401546446
4.35278321173439
6.01971825558614
4.77661077480956
5.74604385771088
2.36009451721527
3.9389469556021
5.31982631801165
5.59081273807287
6.07906306723298
3.85996347778954
4.24391686226955
4.49934103344509
5.14048982266457
5.36780929354124
5.87338444704468
6.15444852716286
]
g=[4.5620857635538
4.43036433356537
3.70912090825776
3.75545079363459
1.94659002327898
4.11125770249527
4.14912488952286
4.03098906151149
3.01918193873448
5.05813626223079
4.00757498622833
4.60497686058731
2.47783183930749
2.29932685304805
2.2664994926372
1.99085075708268
3.25133201543739
2.7399068022911
2.74653418520092
3.61433759175203
2.87409511156452
3.69314126194378
5.08627970070686
3.93364348031753
3.86127225277459
4.81224143297952
2.62492072318273
4.11128804412481
3.9086537913148
4.44395075033754
3.1813933442762
2.57644891118114
3.27416239554183
4.09527108623828
3.18160029897779
3.42939694926035
4.3221621771641
3.30682867178323
3.32708943262992
4.32918822744415
3.35412379973399
2.55014730411195
6.88194981786635
3.29364696882548
2.68268826699893
4.06706692399156
3.30025017819095
5.74320426981588
2.25555728187611
3.71158337231985
4.47450741503766
2.91530560962556
3.578763059274
6.00664237023124
3.3144331911963
4.36517929839789
4.17688380594166
2.71789353373283
2.8939110100738
3.75602447100316
2.53182244718191
3.96862696659689
3.46545455168543
3.33532913042806
4.482396372823
4.72418912511491
2.84370516680693
3.24842556475694
2.94887474404863
2.54343310903865
2.63973356134744
4.49487243609262
2.63149333835876
5.09896375887853
6.83167911021023
5.04497419539195
3.35937565108723
4.20781861156608
2.52322524382984
4.37335509069669
8.69046904999952
5.26088250054987
4.49945504447026
3.95467217595026
4.32685590208085
4.48472152336562
3.78358648445302
3.52888513599272
2.41345578407203
2.84235655370593
3.37517085882163
5.28854722017956
3.48617455821768
5.23638798917963
4.14970845662791
4.92457348680135
4.57363736206569
3.91145788660847
3.58449802093723
2.75757511077479
3.42433289250784
4.82368201475464
4.94422026545281
2.24025397565046
3.81726122839666
2.56546590033973
5.01946051897514
4.64407267001559
2.01260100863833
5.42698459690725
3.76008681436108
1.81467237155086
3.55124363192886
5.45961438367124
3.34073809013046
4.31859757639702
2.03549964089314
3.18027451435545
5.07977234621604
2.51614094297572
4.06178636929318
5.21818371432763
4.30184656458391
3.35894237877829
5.2140052199976
4.85280378528717
4.33058050757073
4.63533125294954
4.5958395606095
4.3092951489229
5.90355659292018
3.66790098150381
4.14570250953731
3.63173859125286
3.91741520328028
2.39434733676411
3.18999434887079
5.01279567236022
3.75555934431885
4.48158948737999
2.61619827113897
4.76955995295965
4.44243867729846
5.35479849334499
4.08103948414671
4.63899862824988
3.58596563473244
2.99306190904264
3.75807413462325
4.08746006225643
5.22450962890135
4.6591835066185
4.80946946622121
3.80373824099358
2.1293763992088
4.0812448352416
4.69127341857472
2.16963010543432
5.00418312838495
4.53945549537302
6.01866010561574
5.21307702667463
5.87375595903832
4.74818439639509
4.54474363895237
5.84896925998765
4.76339530852165
4.55218752158776
5.00625073531939
6.12897939301331
2.12694350717774
4.19957839040588
3.70928616528714
4.087014025835
3.96734210210732
3.47033598102741
4.30826298029895
4.45506048818769
4.86347183432235
5.09239959560266
4.84697456435945
5.72293872922757
]


s=[a,t,c,g]
as=[3.21472881345033 7.73920234155702 5.353333556465 2.55596503346531]
bs=[3.52246934793863 1.85980014068493 12.5705612573442 10.8411699404179]
sas=s*as'

ss=sum((s.*s)')
ass=sum(as.*as)
bss=sum(bs.*bs)
A=sqrt((ss*ass))

neijiA=sas./A'

sbs=s*bs'
B=sqrt(ss*bss)
neijiB=sbs./B'
biji=[neijiA,neijiB]
mmm=[1:182]'
bijjj=neijiA-neijiB


[bijjj,mmm]
%for i = 1:182
  %  if neijiA>neijiB
   %     i,'A'
   % else
   %     i,'B'
        
  %  end
%end