comp40juli.m

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

a2=[9.00900900900901E-02 .036036036036036 9.00900900900901E-02 8.10810810810811E-02 .045045045045045 .036036036036036 9.00900900900901E-03 .045045045045045 .036036036036036 .036036036036036 .018018018018018 8.10810810810811E-02 .117117117117117 .027027027027027 5.40540540540541E-02 .189189189189189 
9.90990990990991E-02 .036036036036036 7.20720720720721E-02 5.40540540540541E-02 .027027027027027 5.40540540540541E-02 .018018018018018 5.40540540540541E-02 .045045045045045 .018018018018018 9.00900900900901E-03 9.00900900900901E-02 9.90990990990991E-02 .045045045045045 5.40540540540541E-02 .216216216216216 
5.40540540540541E-02 .036036036036036 .117117117117117 5.40540540540541E-02 .027027027027027 9.00900900900901E-03 .018018018018018 9.00900900900901E-03 5.40540540540541E-02 9.00900900900901E-03 9.00900900900901E-03 .144144144144144 .135135135135135 9.00900900900901E-03 7.20720720720721E-02 .234234234234234 
.189189189189189 .117117117117117 5.40540540540541E-02 5.40540540540541E-02 .108108108108108 .108108108108108 .018018018018018 5.40540540540541E-02 5.40540540540541E-02 .018018018018018 9.00900900900901E-03 .027027027027027 6.30630630630631E-02 .045045045045045 .027027027027027 .045045045045045 
6.30630630630631E-02 .018018018018018 8.10810810810811E-02 7.20720720720721E-02 .018018018018018 .036036036036036 .027027027027027 .027027027027027 5.40540540540541E-02 .045045045045045 .027027027027027 .108108108108108 9.90990990990991E-02 9.00900900900901E-03 9.00900900900901E-02 .216216216216216 
.153153153153153 6.30630630630631E-02 .027027027027027 9.90990990990991E-02 .036036036036036 5.40540540540541E-02 .018018018018018 .018018018018018 .045045045045045 0 0 8.10810810810811E-02 .108108108108108 9.00900900900901E-03 8.10810810810811E-02 .198198198198198 
.153153153153153 .108108108108108 .018018018018018 7.20720720720721E-02 .045045045045045 5.40540540540541E-02 .027027027027027 6.30630630630631E-02 9.00900900900901E-03 .018018018018018 9.00900900900901E-03 6.30630630630631E-02 .135135135135135 9.00900900900901E-03 .045045045045045 .162162162162162 
8.10810810810811E-02 6.30630630630631E-02 .036036036036036 9.90990990990991E-02 5.40540540540541E-02 7.20720720720721E-02 .036036036036036 .027027027027027 .027027027027027 .036036036036036 .018018018018018 8.10810810810811E-02 .108108108108108 .018018018018018 7.20720720720721E-02 .162162162162162 
9.00900900900901E-02 .045045045045045 9.00900900900901E-03 6.30630630630631E-02 0 .045045045045045 .036036036036036 7.20720720720721E-02 .036036036036036 .027027027027027 .027027027027027 .117117117117117 7.20720720720721E-02 .036036036036036 .135135135135135 .18018018018018 
6.36363636363636E-02 1.81818181818182E-02 3.63636363636364E-02 6.36363636363636E-02 1.81818181818182E-02 3.63636363636364E-02 5.45454545454545E-02 2.72727272727273E-02 5.45454545454545E-02 3.63636363636364E-02 4.54545454545455E-02 .136363636363636 4.54545454545455E-02 3.63636363636364E-02 .136363636363636 .181818181818182]
b2=[.154545454545455 .145454545454545 2.72727272727273E-02 2.72727272727273E-02 .163636363636364 .3 9.09090909090909E-03 1.81818181818182E-02 9.09090909090909E-03 9.09090909090909E-03 9.09090909090909E-03 1.81818181818182E-02 2.72727272727273E-02 4.54545454545455E-02 0 2.72727272727273E-02 
.136363636363636 .109090909090909 9.09090909090909E-03 6.36363636363636E-02 .154545454545455 .309090909090909 1.81818181818182E-02 1.81818181818182E-02 9.09090909090909E-03 9.09090909090909E-03 0 9.09090909090909E-03 2.72727272727273E-02 7.27272727272727E-02 0 4.54545454545455E-02 
6.36363636363636E-02 .1 4.54545454545455E-02 4.54545454545455E-02 .127272727272727 .345454545454545 1.81818181818182E-02 2.72727272727273E-02 2.72727272727273E-02 2.72727272727273E-02 1.81818181818182E-02 1.81818181818182E-02 3.63636363636364E-02 4.54545454545455E-02 1.81818181818182E-02 2.72727272727273E-02 
8.18181818181818E-02 .127272727272727 9.09090909090909E-03 7.27272727272727E-02 .136363636363636 .281818181818182 6.36363636363636E-02 1.81818181818182E-02 2.72727272727273E-02 4.54545454545455E-02 0 9.09090909090909E-03 5.45454545454545E-02 4.54545454545455E-02 9.09090909090909E-03 9.09090909090909E-03 
.136363636363636 .127272727272727 0 1.81818181818182E-02 .136363636363636 .481818181818182 0 2.72727272727273E-02 0 0 0 0 1.81818181818182E-02 3.63636363636364E-02 0 9.09090909090909E-03 
.163636363636364 .154545454545455 3.63636363636364E-02 9.09090909090909E-03 .136363636363636 .227272727272727 4.54545454545455E-02 4.54545454545455E-02 1.81818181818182E-02 5.45454545454545E-02 0 9.09090909090909E-03 4.54545454545455E-02 2.72727272727273E-02 0 1.81818181818182E-02 
.172727272727273 .109090909090909 5.45454545454545E-02 1.81818181818182E-02 .1 5.45454545454545E-02 6.36363636363636E-02 4.54545454545455E-02 4.54545454545455E-02 7.27272727272727E-02 9.09090909090909E-02 2.72727272727273E-02 3.63636363636364E-02 2.72727272727273E-02 3.63636363636364E-02 3.63636363636364E-02 
8.18181818181818E-02 .118181818181818 7.27272727272727E-02 1.81818181818182E-02 .154545454545455 .309090909090909 1.81818181818182E-02 9.09090909090909E-03 3.63636363636364E-02 3.63636363636364E-02 1.81818181818182E-02 2.72727272727273E-02 1.81818181818182E-02 3.63636363636364E-02 9.09090909090909E-03 2.72727272727273E-02 
2.72727272727273E-02 .136363636363636 2.72727272727273E-02 1.81818181818182E-02 .145454545454545 .318181818181818 9.09090909090909E-02 9.09090909090909E-03 1.81818181818182E-02 8.18181818181818E-02 1.81818181818182E-02 2.72727272727273E-02 2.72727272727273E-02 2.72727272727273E-02 9.09090909090909E-03 9.09090909090909E-03 
6.36363636363636E-02 6.36363636363636E-02 6.36363636363636E-02 9.09090909090909E-03 9.09090909090909E-02 .327272727272727 .1 3.63636363636364E-02 2.72727272727273E-02 .136363636363636 9.09090909090909E-03 0 1.81818181818182E-02 3.63636363636364E-02 0 9.09090909090909E-03]
c2=[5.30973451327434E-02 .079646017699115 4.42477876106195E-02 8.84955752212389E-02 9.73451327433628E-02 .185840707964602 6.19469026548673E-02 1.76991150442478E-02 6.19469026548673E-02 4.42477876106195E-02 4.42477876106195E-02 4.42477876106195E-02 6.19469026548673E-02 4.42477876106195E-02 4.42477876106195E-02 1.76991150442478E-02 
8.57142857142857E-02 3.80952380952381E-02 9.52380952380952E-02 7.61904761904762E-02 .104761904761905 6.66666666666667E-02 3.80952380952381E-02 9.52380952380952E-03 2.85714285714286E-02 .019047619047619 7.61904761904762E-02 .114285714285714 7.61904761904762E-02 8.57142857142857E-02 2.85714285714286E-02 4.76190476190476E-02 
2.94117647058824E-02 5.88235294117647E-02 3.92156862745098E-02 4.90196078431373E-02 3.92156862745098E-02 9.80392156862745E-02 2.94117647058824E-02 1.96078431372549E-02 0 1.96078431372549E-02 .127450980392157 9.80392156862745E-02 .107843137254902 9.80392156862745E-03 4.90196078431373E-02 .215686274509804 
1.73913043478261E-02 3.47826086956522E-02 4.34782608695652E-02 .11304347826087 .130434782608696 .226086956521739 1.73913043478261E-02 2.60869565217391E-02 2.60869565217391E-02 9.56521739130435E-02 4.34782608695652E-02 2.60869565217391E-02 3.47826086956522E-02 4.34782608695652E-02 8.69565217391304E-02 2.60869565217391E-02 
6.66666666666667E-02 3.80952380952381E-02 3.80952380952381E-02 9.52380952380952E-02 5.71428571428571E-02 9.52380952380952E-02 .019047619047619 4.76190476190476E-02 7.61904761904762E-02 4.76190476190476E-02 7.61904761904762E-02 2.85714285714286E-02 4.76190476190476E-02 3.80952380952381E-02 9.52380952380952E-02 .123809523809524 
3.50877192982456E-02 5.26315789473684E-02 3.50877192982456E-02 9.64912280701754E-02 7.89473684210526E-02 .245614035087719 .043859649122807 1.75438596491228E-02 7.89473684210526E-02 6.14035087719298E-02 1.75438596491228E-02 .043859649122807 2.63157894736842E-02 2.63157894736842E-02 .114035087719298 1.75438596491228E-02 
5.76923076923077E-02 4.80769230769231E-02 4.80769230769231E-02 7.69230769230769E-02 6.73076923076923E-02 .105769230769231 2.88461538461538E-02 2.88461538461538E-02 2.88461538461538E-02 2.88461538461538E-02 7.69230769230769E-02 6.73076923076923E-02 7.69230769230769E-02 4.80769230769231E-02 4.80769230769231E-02 .153846153846154 
3.41880341880342E-02 .094017094017094 5.12820512820513E-02 6.83760683760684E-02 .11965811965812 .239316239316239 5.12820512820513E-02 3.41880341880342E-02 2.56410256410256E-02 6.83760683760684E-02 2.56410256410256E-02 2.56410256410256E-02 7.69230769230769E-02 3.41880341880342E-02 1.70940170940171E-02 2.56410256410256E-02 
1.98019801980198E-02 3.96039603960396E-02 1.98019801980198E-02 6.93069306930693E-02 3.96039603960396E-02 8.91089108910891E-02 2.97029702970297E-02 2.97029702970297E-02 1.98019801980198E-02 9.9009900990099E-03 8.91089108910891E-02 8.91089108910891E-02 6.93069306930693E-02 4.95049504950495E-02 7.92079207920792E-02 .247524752475248 
9.34579439252336E-02 2.80373831775701E-02 5.60747663551402E-02 .102803738317757 7.47663551401869E-02 6.54205607476635E-02 5.60747663551402E-02 5.60747663551402E-02 8.41121495327103E-02 7.47663551401869E-02 2.80373831775701E-02 5.60747663551402E-02 3.73831775700935E-02 8.41121495327103E-02 9.34579439252336E-02 0 
2.67857142857143E-02 4.46428571428571E-02 5.35714285714286E-02 .116071428571429 .151785714285714 .169642857142857 1.78571428571429E-02 8.92857142857143E-03 3.57142857142857E-02 .0625 3.57142857142857E-02 4.46428571428571E-02 2.67857142857143E-02 7.14285714285714E-02 7.14285714285714E-02 5.35714285714286E-02 
5.45454545454545E-02 2.72727272727273E-02 2.72727272727273E-02 6.36363636363636E-02 6.36363636363636E-02 .154545454545455 7.27272727272727E-02 4.54545454545455E-02 4.54545454545455E-02 5.45454545454545E-02 6.36363636363636E-02 6.36363636363636E-02 9.09090909090909E-03 .1 6.36363636363636E-02 8.18181818181818E-02 
5.40540540540541E-02 7.20720720720721E-02 7.20720720720721E-02 7.20720720720721E-02 .108108108108108 .153153153153153 .018018018018018 5.40540540540541E-02 .036036036036036 .045045045045045 .027027027027027 7.20720720720721E-02 7.20720720720721E-02 6.30630630630631E-02 6.30630630630631E-02 9.00900900900901E-03 
7.84313725490196E-02 9.80392156862745E-03 4.90196078431373E-02 8.82352941176471E-02 4.90196078431373E-02 7.84313725490196E-02 9.80392156862745E-03 2.94117647058824E-02 2.94117647058824E-02 3.92156862745098E-02 9.80392156862745E-02 6.86274509803922E-02 7.84313725490196E-02 3.92156862745098E-02 6.86274509803922E-02 .176470588235294 
.058252427184466 3.88349514563107E-02 4.85436893203883E-02 9.70873786407767E-02 7.76699029126214E-02 6.79611650485437E-02 3.88349514563107E-02 1.94174757281553E-02 3.88349514563107E-02 .029126213592233 3.88349514563107E-02 9.70873786407767E-02 6.79611650485437E-02 6.79611650485437E-02 .087378640776699 .116504854368932 
4.76190476190476E-02 .019047619047619 3.80952380952381E-02 .123809523809524 8.57142857142857E-02 6.66666666666667E-02 5.71428571428571E-02 0 5.71428571428571E-02 3.80952380952381E-02 .104761904761905 .104761904761905 3.80952380952381E-02 8.57142857142857E-02 9.52380952380952E-02 2.85714285714286E-02 
5.76923076923077E-02 2.88461538461538E-02 2.88461538461538E-02 .105769230769231 5.76923076923077E-02 5.76923076923077E-02 9.61538461538462E-03 6.73076923076923E-02 5.76923076923077E-02 5.76923076923077E-02 9.61538461538462E-02 3.84615384615385E-02 4.80769230769231E-02 5.76923076923077E-02 .115384615384615 .105769230769231 
3.41880341880342E-02 5.98290598290598E-02 .094017094017094 3.41880341880342E-02 .102564102564103 .273504273504274 1.70940170940171E-02 4.27350427350427E-02 5.12820512820513E-02 3.41880341880342E-02 4.27350427350427E-02 3.41880341880342E-02 2.56410256410256E-02 6.83760683760684E-02 1.70940170940171E-02 5.98290598290598E-02 
8.49056603773585E-02 4.71698113207547E-02 5.66037735849057E-02 8.49056603773585E-02 4.71698113207547E-02 6.60377358490566E-02 8.49056603773585E-02 2.83018867924528E-02 .113207547169811 1.88679245283019E-02 9.43396226415094E-02 5.66037735849057E-02 2.83018867924528E-02 9.43396226415094E-02 4.71698113207547E-02 3.77358490566038E-02 
3.44827586206897E-02 4.31034482758621E-02 7.75862068965517E-02 4.31034482758621E-02 .103448275862069 .275862068965517 8.62068965517241E-03 3.44827586206897E-02 1.72413793103448E-02 6.03448275862069E-02 8.62068965517241E-02 3.44827586206897E-02 4.31034482758621E-02 5.17241379310345E-02 1.72413793103448E-02 6.03448275862069E-02]



juna2=mean(a2)
junb2=mean(b2)
for m=1:20
sb2=c2(m,:)-junb2;
sa2=c2(m,:)-juna2;
   if (sum(sb2.*sb2)/4)^(1/2)>(sum(sa2.*sa2)/4)^(1/2)
 
         m,out='A'
   else
        m, out='B'
   end
end