yanzhengdangezifu.m

来自「数学建模竞赛：DNA序列的分析」· M 代码 · 共 56 行

c=[.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
.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]
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
   if (sum(sb.*sb)/4)^(1/2)>(sum(sa.*sa)/4)^(1/2)
 
         out='A'
   else
         out='B'
   end
end