yanzhengfangchaneiji20.m

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

a=[2.23516615533062
3.18594511964954
2.46552585344185
2.02734923669708
2.99276906337035
2.33411217781108
2.52723430433715
3.460117598161
4.80011724243539
5.25212293474148
2.15211034739588
2.56675579167899
3.92921702617028
3.05262974942864
5.68529512384021
2.31873202820496
2.2045149584413
3.30037469361883
4.0646942234853
4.08014939035558
]
t=[5.5137191532666
6.07590871118606
17.9962959151784
4.29478303180994
9.04935959562383
13.6836979557278
4.81117121094777
5.19994939246629
5.2563453717071
5.03395066927154
1.57815272953407
1.54695887502102
2.03316655337585
1.57205238691265
1.30487222032932
1.50319387862468
3.41700540326262
1.92304789528165
1.0630402837187
1.23008729531621
]
c=[5.21091112678405
3.58612522336394
3.01314642581458
10.3677822648293
2.19241115365405
5.75348278024399
5.82713956890991
3.95191688409859
4.59743011217828
3.15565018321387
28.3313724811677
50.9116882454314
16.0350310951436
8.26243651369601
0
3.81491434092188
4.18072500153521
5.94099777356081
5.14735253263637
4.68797655529617
]
g=[2.03947749154979
2.1468793419964
1.8058389800225
4.9689091840697
2.51536342107788
2.17378929784062
2.32222146889725
2.18620382844034
1.64071688489886
2.19359050952251
8.8950672972284
8.09291283131346
7.15756727043269
6.18588324454572
11.2190314495801
17.0473197632681
8.08777397755946
9.77383127426383
12.876704695283
20.1270961641266
]


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]