查看: 1663  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

heapyssb

木虫 (著名写手)

[求助] 怎样将txt文件读入到matlab中,谢谢啦

比如有这样的txt文件
>sp|Q9NRA8|4ET_HUMAN E
241  WFSAGPTSQSETIELTGFDDKILEEDHKGRKRTRRRTASVK
>sp|Q9NRA8|4ET_HUMAN E
248  SQSETIELTGFDDKILEEDHKGRKRTRRRTASVKEGIVECN
>sp|Q9NRA8|4ET_HUMAN E
251  ETIELTGFDDKILEEDHKGRKRTRRRTASVKEGIVECNGGV

请问该怎样编写程序来将文件里面的内容读入到matlab中,进行分析?
回复此楼

» 收录本帖的淘帖专辑推荐

source

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

★
xzhdty(金币+1): 欢迎讨论交流 2011-10-28 07:36:10
引用回帖:
1楼: Originally posted by heapyssb at 2011-10-25 14:38:15:
比如有这样的txt文件
>sp|Q9NRA8|4ET_HUMAN E
241  WFSAGPTSQSETIELTGFDDKILEEDHKGRKRTRRRTASVK
>sp|Q9NRA8|4ET_HUMAN E
248  SQSETIELTGFDDKILEEDHKGRKRTRRRTASVKEGIVECN
>sp|Q9NRA8|4ET_HUMAN E
...

格式标准的,可以直接load(test.txt),
你这样的如果load进来格式不正确
还是用fopen打开后fscanf函数读取好点
具体函数用法help会讲的很详细的
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
4楼2011-10-28 00:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

zuoyh09

金虫 (正式写手)

【答案】应助回帖

★ ★
微尘、梦想(金币+2): 欢迎参与应助 2011-10-26 19:19:03
heapyssb(金币+3): 谢谢,希望能再详细点 2011-10-27 17:02:47
转引自:http://ncforest.blog.163.com/blo ... 642009102453056325/

文件内容形式一(假定文件名为test1.txt):

你好

欢迎来到

动力学与控制技术论坛

www.dytrol.com

1 11 111 1111

2 22 222 2222

3 33 333 3333

4 44 444 4444

5 55 555 5555

这样的文件怎么读入数据呢?

方法有多种,现举两个比较简单实用的。

方法一:

file/import data....../next/finish

>> whos

  Name           Size                   Bytes  Class

  data           5x4                      160  double array

  textdata       4x1                      300  cell array

Grand total is 54 elements using 460 bytes

>> data

data =

           1          11         111        1111

           2          22         222        2222

           3          33         333        3333

           4          44         444        4444

           5          55         555        5555

>> textdata

textdata =

    '你好'

    '欢迎来到'

    '动力学与控制技术论坛'

    'www.dytrol.com'

方法二:

[a1,a2,a3,a4]=textread('test1.txt','%s%s%s%s','headerlines',4)

说明:%s可以是其他形式,跟读入的数据类型有关,比如这里也可以用%n,%f等。

这里%s的个数和[a1,a2,a3,a4]对应。

>> [a1,a2,a3,a4]=textread('test1.txt','%s%s%s%s','headerlines',4)

a1 =

    '1'

    '2'

    '3'

    '4'

    '5'

a2 =

    '11'

    '22'

    '33'

    '44'

    '55'

a3 =

    '111'

    '222'

    '333'

    '444'

    '555'

a4 =

    '1111'

    '2222'

    '3333'

    '4444'

    '5555'

因以字符串的形式读入,所以有''。

文件内容形式二(假定文件名为test2.txt):

你好

1 11 111 1111

欢迎来到

2 22 222 2222

动力学与仿真控制论坛

3 33 333 3333

www.dytrol.com

4 44 444 4444

5 55 555 5555

说明:这种内容格式的文件用上面的方法是不行的。

以下是由chinamaker编写的一种方法:

fidin=fopen('test2.txt');                               % 打开test2.txt文件            

fidout=fopen('mkmatlab.txt','w');                       % 创建MKMATLAB.txt文件

while ~feof(fidin)                                      % 判断是否为文件末尾              

    tline=fgetl(fidin);                                 % 从文件读行  

    if double(tline(1))>=48&&double(tline(1))<=57       % 判断首字符是否是数值

       fprintf(fidout,'%s\n\n',tline);                  % 如果是数字行,把此行数据写入文件MKMATLAB.txt

       continue                                         % 如果是非数字继续下一次循环

    end

end

fclose(fidout);

MK=importdata('MKMATLAB.txt');      % 将生成的MKMATLAB.txt文件导入工作空间,变量名为MK,实际上它不显示出来

>> MK

MK =

           1          11         111        1111

           2          22         222        2222

           3          33         333        3333

           4          44         444        4444

           5          55         555        5555
2楼2011-10-25 17:24:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zuoyh09

金虫 (正式写手)

【答案】应助回帖

★
jjdg(金币+1): 感谢参与 2011-10-28 00:45:05
heapyssb(金币+2): 谢谢 2011-10-28 17:34:16
下面的程序可以完成读txt数据,需要注意的是,在txt中不能有空行。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
clc
fidin=fopen('data.txt'); % 打开data.txt文件
fidout=fopen('vr.txt','w'); % 创建vr.txt文件
while ~feof(fidin) % 判断是否为文件末尾
tline=fgetl(fidin); % 从文件读行
% disp(tline)
if double(tline(1))>=48 && double(tline(1))<=57 % 判断首字符是否是数值
fprintf(fidout,'%s\n\n',tline); % 如果是数字行,把此行数据写入文件vr.txt
continue % 如果是非数字继续下一次循环
end
end
fclose(fidout);
MK=importdata('vr.txt'); % 将生成的vr.txt文件导入工作空间,变量名为MK,实际上它不显示出来
3楼2011-10-27 22:36:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦在潇湘

铁杆木虫 (著名写手)

这些在help中都详细讲述的
5楼2012-04-23 13:18:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 7/350 2026-10-12 14:47 by chuans0923
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见