compdan.m

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

c=[.274336 .362832 .19469 .168142 
.295238 .219048 .238095 .247619 
.176471 .186275 .254902 .382353 
.208696 .408696 .191304 .191304 
.247619 .219048 .228571 .304762 
.219298 .385965 .210526 .184211 
.230769 .230769 .201923 .336538 
.25641 .444444 .145299 .153846 
.148515 .188119 .217822 .445545 
.28972 .252336 .242991 .214953 
.241071 .357143 .178571 .223214 
.172727 .336364 .227273 .263636 
.27027 .333333 .189189 .207207 
.235294 .166667 .235294 .362745 
.242718 .203883 .213592 .339806 
.228571 .209524 .304762 .257143 
.221154 .201923 .25 .326923 
.222222 .435897 .17094 .17094 
.273585 .235849 .283019 .207547 
.198276 .431034 .198276 .172414]
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:20  
sb=c(m,:)-junb;
sa=c(m,:)-juna;
m,(sum(sb.*sb)/4)^(1/2)-(sum(sa.*sa)/4)^(1/2)
end