实验目的
学习利用MATLAB提取DNA序列特征建立向量的方法，掌握利用FCM命令进行DNA分类的方法，学会做出分类图形直接给出分类结果的MATLAB编程。 知识扩展
DNA序列分类
DNA(Deoxyribonucleic acid)，中文译名为脱氧核苷酸，是染体的主要化学成分，同时也是基因组成的，有时被称为“遗传微粒”。DNA是一种分子，可组成遗传指令，以引导生物发育与生命机能运作。主要功能是长期性的资讯储存，可比喻为“蓝图”或“食谱”。DNA分子是由两条核苷酸链以互补配对原则所构成的双螺旋结构的分子化合物。其中两条DNA链中对应的碱基A-T以双键形式连接，C-G以三键形式连接，糖-磷酸-糖 形成的主链在螺旋外侧，配对碱基在螺旋内侧。
FCM算法中样本点隶属于某一类的程度是用隶属度来反映的，不同的样本点以不同的隶属度属于每一类；但是算法中的概率约束∑uij=1使得样本的典型性反映不出来，不适用于有噪音，样本分布不均衡，存在两个或者两个以上样本分别距两个类的距离相等的样本等等。 欧氏距离（ Euclidean distance）也称欧几里得距离，它是一个通常采用的距离定义，它是在m维空间中两个点之间的真实距离。
公式gal芯片
在二维和三维空间中的欧式距离的就是两点之间的距离，二维的公式是
d = sqrt((x1-x2)^+(y1-y2)^)
三维的公式是
d=sqrt(x1-x2)^+(y1-y2)^+(z1-z2)^)
推广到n维空间，欧式距离的公式是
d=sqrt( ∑(xi1-xi2)^ ) 这里i=
xi1表示第一个点的第i维坐标,xi2表示第二个点的第i维坐标
n维欧氏空间是一个点集,它的每个点可以表示为(x(1),x(2),...x(n)),其中x(i)(i=)是实数,称为x的第i个坐标,两个点x和y=(y(1),y(2)...y(n))之间的距离d(x,y)定义为上面的公式.
欧氏距离判别准则如下：
若dA<dB,则将Xi点判为A类
若dA>dB,则将Xi点判为B类
若dA=dB,则将Xi点判为不可判别点。
欧氏距离看作信号的相似程度。 距离越近
昂立口译
1.问题的提出m序列
2000年6月，人类基因组计划中DNA全序列草图完成，预计2001以完精确的全序列图，此
后人类将拥有一本记录着自身生老病死及遗的全部信息的“天书”，这本大自然写成的“天书淡马锡模式”，是由4个字符A,T,C,G按一定顺序排成的长约30亿的序列，其中没有“断句”也没有标点符号，除了这4种碱基以外，人们对它包含的“内容”知之甚少，难以读懂，破译这部世界上最巨量信息的“天书”是21实际最重要的任务之一。在这个目标中，研究DNA全序列具有什么结构，由这4个字符排成的看似随机的序列中隐藏着什么规律，又是解读这部天书的基础，是生物信息学最重要的课题之一。
2．问题的分析
这是一个比较典型的分类问题，为了表述的严格和方便，我们用数学的方法来重述这个问题。在这里问题的关键就是要从已知的20个字母序列中提取用于分类的特征。知道了这些特征，我们就可以比较容易的，对那些未标明类型的序列进行分类，下面我们将首先对用于分类的标准问题进行必要的讨论。
3．分类的方法
为了在众多可能的分类中寻求合理的分类结果，为此，就要确定合理的聚类准则。定义目标函数为
显然，J(U,V)表示了各类中样本到聚类中心的加权距离平方和，权重是样本XK对第i类隶
度的m次方，聚类准则取为求J(U,V)的极小值(min){J(U,V)}。
其中，U=[uik]为模糊分类矩阵，i=1，2；k=1，2，···，20；且满足0≤ uik ≤1和
若uik=max{uik}>0.5，则xk∈第j类。
在MATLAB中，我们只要直接调用如下程序即可：
[Center,U,obj_fcm]=fcm(data,cluster_n)
data:要聚类的数据函数，每一行为一个样本
cluster_n:聚类数（大于1）
Center:最终的聚类中心矩阵，其每一行为聚类中心的坐标值
U:最终的模糊分区矩阵
obj_fcm:在迭代过程中的目标函数值
4. 对DNA序列组合分类的分析
（1）提取DNA序列特征建立两类序列的特征向量
（2）确定两类序列的中心
（3）分类方法
（4）回代误判率
（5）未知的20个序列判别结果
5.提取DNA序列特征建立两类序列的特征向量
为了对DNA序列进行分类，我们首先对已知的两类DNA序列进行研究，从中到两类序列的特征。由于在不用于编码蛋白质的序列片段中，A和T的含量特别多些，于是我们利用MATLAB软件，通过编程计算出A，B两类序列中4中碱基对含量的百分比，对每个序列构造四维向量
xk =( xk1，xk2，xk3，xk4 ) (k=1,2···,20)
其中， xk1，xk2，xk3，xk4分别表示第k个序列所含有的碱基对A,T,C,G含量的百分比，利用MATLAB软件，我们得到A,B两类序列的特征矩阵A=(xkj)20×4
6. 确定两类序列的中心
为第k个序列到第i类中心的欧式距离，实际计算时要对取定的初始值进行迭代计算直至max｛|uikt-uikt-1|｝<ɛ ，ɛ为事先指定的精度。
回代误判率
用欧氏距离作为判据虽然简洁直观，但是存在着明显的缺陷：从概率统计的角度来看，用欧氏距离描述随机点之间的距离并不是很好。因此对待分类样本是随机样本，具有一定的统计性质时，这个模型并不是很能很好的描述两个随机点之间的接近程度。
我们对于已知的A，B两类序列利用上述方宣传法进行判别，结果如3.15表所示
未知的20个序列判别结果 利用模糊均值聚类方法，对于未知的20个序列进行判别，结果如表3.16
MATLAB编程
(1) 提取特征建立特征矩阵
A1='aggcacggaaaaacgggaataacggaggaggacttggcacggcattacacggaggacgaggtaaaggaggcttgtctacggccggaagtgaagggggatatgaccgcttgg';
A2='cggaggacaaacgggatggcggtattggaggtggcggactgttcggggaattattcggtttaaacgggacaaggaaggcggctggaacaaccggacggtggcagcaaagga';
A3='gggacggatacggattctggccacggacggaaaggaggacacggcggacatacacggcggcaacggacggaacggaggaaggagggcggcaatcggtacggaggcggcgga';
A4='atggataacggaaacaaaccagacaaacttcggtagaaatacagaagcttagatgcatatgttttttaaataaaatttgtattattatggtatcataaaaaaaggttgcga';
A5='cggctggcggacaacggactggcggattccaaaaacggaggaggcggacggaggctacaccaccgtttcggcggaaaggcggagggctggcaggaggctcattacggggag';
A6='atggaaaattttcggaaaggcggcaggcaggaggcaaaggcggaaaggaaggaaacggcggatatttcggaagtggatattaggagggcggaataaaggaacggcggcaca';
A7='atgggattattgaatggcggaggaagatccggaataaaatatggcggaaagaacttgttttcggaaatggaaaaaggactaggaatcggcggcaggaaggatatggaggcg';
A8='atggccgatcggcttaggctggaaggaacaaataggcggaattaaggaaggcgttctcgcttttcgacaaggaggcggaccataggaggcggattaggaacggttatgagg';
A9='atggcggaaaaaggaaatgtttggcatcggcgggctccggcaactggaggttcggccatggaggcgaaaatcgtgggcggcggcagcgctggccggagtttgaggagcgcg';明报
A10='tggccgcggaggggcccgtcgggcgcggatttctacaagggcttcctgttaaggaggtggcatccaggcgtcgcacgctcggcgcggcaggaggcacgcgggaaaaaacg';
诚实的商人A11='gttagatttaacgttttttatggaatttatggaattataaatttaaaaatttatattttttaggtaagtaatccaacgtttttattactttttaaaattaaatatttatt';
A12='gtttaattactttatcatttaatttaggttttaattttaaatttaatttaggtaagatgaatttggttttttttaaggtagttatttaattatcgttaaggaaagttaaa';
A13='gtattacaggcagaccttatttaggttattattattatttggattttttttttttttttttttaagttaaccgaattattttctttaaagacgttacttaatgtcaatgc';
A14='gttagtcttttttagattaaattattagattatgcagtttttttacataagaaaatttttttttcggagttcatattctaatctgtctttattaaatcttagagatatta';
A15='gtattatatttttttatttttattattttagaatataatttgaggtatgtgtttaaaaaaaatttttttttttttttttttttttttttttttaaaatttataaatttaa';
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议