Matlab 基本操作练习
及初识 Bioinformatics 工具箱
生物医学工程 2 班
20090260230
资艳格
2011-11-18
实验一、Matlab 基本操作练习及初识 Bioinformatics 工具
箱
一、实验目的:
1、熟悉 Matlab 实验环境,掌握其基本操作方法
2、练习 Matlab 常用指令的使用
3、了解 Matlab 中 BioInformatics 工具箱的功能及用法
4、进一步认识生物信息学
二、实验设备:
PC 机,Matlab 实验环境,Internet 网络
三、 实验内容:
比较以下两个碱基序列,分别计算各碱基的百分含量,并作出直
方图。要求在实验报告中写出完整的 MATLAB 实验程序。
四、程序代码:
function ziyange
close all
A='ATGTATCCGGATTTAAAAGGAAAAGTCGTCGCTATTAC
AGGAGCTGCTTCAGGGCTCGGAAAGGCGATGGCCATTCGC
TTCGGCAAGGAGCAGGCAAAAGTGGTTATCAACTATTATA
GTAATAAACAAGATCCGAACGAGGTAAAAGAAGAGGTCAT
CAAGGCGGGCGGTGAAGCTGTTGTCGTCCAAGGAGATGT
CACGAAAGAGGAAGATGTAAAAAATATCGTGCAAACGGCA
ATTAAGGAGTTCGGCACACTCGATATTATGATTAATAATGC
CGGTCTTGAAAATCCTGTGCCATCTCACGAAATGCCGCTCA
AGGATTGGGATAAAGTCATCGGCACGAACTT';
B='AACGGGTGCCTTTTTAGGAAGCCGTGAAGCGATTAAAT
ATTTCGTAGAAAACGATATCAAGGGAAATGTCATTAACATG
TCCAGTGTGCACGAAGTGATTCCTTGGCCGTTATTTGTCCA
CTATGCGGCAAGTAAAGGCGGGATAAAGCTGATGACAGAA
ACATTAGCGTTGGAATACGCGCCGAAGGGCATTCGCGTCA
ATAATATTGGGCCAGGTGCGATCAACACGCCAATCAATGCT
GAAAAATTCGCTGACCCTAAACAGAAAGCTGATGTAGAAA
GCATGATTCCAATGGGATATATCGGCGAACCGGAGGAGAT
CGCCGCAGTAGCAGCCTGGCTTGCTTCGA';
m=length(A)
n=length(B)
M=zeros(1,4);
N=zeros(1,4);
for i=1:m
if A(i)=='A'
M (1)= M(1)+1;
elseif A(i)=='T'
M(2)= M(2)+1;
elseif A(i)=='G'