查看: 2234|回复: 7

关于matlab如何读取“*.text”文件的问题

[复制链接]

字体大小: 正常放大

欣赏世界

1 主题	4 听众	35 积分

升级 31.58%

该用户从未签到

群组: 数学天下

电梯直达

1^#

发表于 2009-7-8 19:00 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

怎么读取很多字母的个数，比如“tgacctcttgtcctgtatagcaacctatttggtaatgattccagcactcacagaaaagct
tgcacacatacacacacaccccacccctccccactaacaaatgcaagttggtaaacaaat
tccaaaaaggcataacaaaccttatatatatagacaaatatatattaaagttttttagtc
tgtactagaaagagcttcagacagaactgaccaccattccattgctcatcaatttcctgg
gacagcacctgagcgtgcgcttacgcgcgtacacacacatagacacgcactgcgatacaa
gtcctgatttgggagtccgtccttttaaaaacagccacatgctttcacgctctgagaccc
acccgtttctgtgagcagggggagggcaaggaaagccctggcctcagtccagccttttct
ctgcttccacctgctcaggctgtgtgctcttggttctgtcctgcacttgtgtgaattcca
aaactgtttttttaaaaaatggcccgcaccccaaatgtctccctgccccatactttgcaa”，怎么样去读取其中的字母t,g,a,c,的个数，并求出各占总数的比例？谢谢！

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

zhaoying

1 主题	4 听众	68 积分

升级 66.32%

该用户从未签到

2^#

发表于 2009-7-8 20:17 |只看该作者 |招呼Ta 关注Ta

好像还不能读吧~~只能读一些数字~~

我也说一句

发表

使用道具举报

mcm-dlu-edu

98 主题	230 听众	1万积分

mcm

升级 0%

TA的每日心情

	怒 2015-11-22 11:00

签到天数: 679 天

[LV.9]以坛为家II

网络挑战赛参赛者

群组: 2015美赛优秀论文解析

群组: 第六届国赛赛前冲刺培

群组: 2014年地区赛数学建模

群组: 2013年美赛优秀论文解

群组: 2014美国中学生数学建

3^#

发表于 2009-7-9 10:06 |只看该作者 |招呼Ta 关注Ta

可以读!!用fopen!!
把你上面的字段按放入abc.text中,再用以下程序,就可能得到每个字母的个数num和比例rate

alp={'t';'g';'a';'c'};
num=zeros(1,length(alp));
fid = fopen('abc.txt','r');
tline=1;
total=0;
i=0;
while tline~=-1
i=i+1;
tline = fgetl(fid);
for j=1:length(tline)
      for k=1:length(alp)
         if strcmp(tline(j),alp(k));
            num(k)=num(k)+1;
         end
      end
end
if tline~=-1
   total=total+j;
end
end

num
rate=num/total

我也说一句

发表

使用道具举报

mcm-dlu-edu

98 主题	230 听众	1万积分

mcm

升级 0%

TA的每日心情

	怒 2015-11-22 11:00

签到天数: 679 天

[LV.9]以坛为家II

网络挑战赛参赛者

群组: 2015美赛优秀论文解析

群组: 第六届国赛赛前冲刺培

群组: 2014年地区赛数学建模

群组: 2013年美赛优秀论文解

群组: 2014美国中学生数学建

4^#

发表于 2009-7-9 10:08 |只看该作者 |招呼Ta 关注Ta

1# 欣赏世界

可以读!!用fopen!!
把你上面的字段按放入abc.text中,再用以下程序,就可能得到每个字母的个数num和比例rate

alp={'t';'g';'a';'c'};
num=zeros(1,length(alp));
fid = fopen('abc.txt','r');
tline=1;
total=0;
i=0;
while tline~=-1
i=i+1;
tline = fgetl(fid);
for j=1:length(tline)
      for k=1:length(alp)
         if strcmp(tline(j),alp(k));
            num(k)=num(k)+1;
         end
      end
end
if tline~=-1
   total=total+j;
end
end

num
rate=num/total