fangchaneiji20.asv

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

a=[2.23516615533062
4.92442890089805
4.49205648100861
3.15524255098646
5.02908208987154
3.74887370473507
3.79188322319904
5.03778823937922
6.84870512613603
7.46100976186646
]
t=[2.26172160430602
5.63154381269344
7.66197639992066
3.37967782816975
5.53344807683669
3.55567956132961
6.06030302272689
2.75317997958724
7.46363734066314
4.65617869072913
3.78709720804179
3.66450152525162
3.48055824750653
10.6957935656968
5.66150435652463
6.28490254498827
7.11188774998741
2.97266457809783
5.54428809088981
3.20806275915336
]
c=[2.70801280154532
4.67183978283949
5.8130600088651
6.47688196588451
5.60438140289282
5.30437212461227
6.28280862437543
8.767363723872
6.84835746730557
4.94553670831926
6.557438524302
5.67986527703056
6.57747351912091
5.81534021517072
6.77126280689208
4.33205109234259
5.00416493200081
7.56453714527348
4.02225949294317
7.19126454208755

]
g=[3.71756122181041
4.48144321991625
3.36917665778323
6.00416522091123
4.09471203708718
6.06543268730356
3.94277244403663
7.75806032459145
3
5.01901147542783
5.46928180644048
5.02217305777312
5.40282113256363
3.37638860322683
3.55050145839438
4.39545219516718
3.72491610643783
7.0789201938651
5.09901951359278
7.30296743340222]


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]