wildgoat的个人空间 https://blog.eetop.cn/wildgoat [收藏] [复制] [分享] [RSS]

空间首页 动态 记录 日志 相册 主题 分享 留言板 个人资料

日志

MATLAB小技:从文本数据第N行开始读取

热度 1已有 31437 次阅读| 2011-8-8 01:16 |个人分类:Matlab

因为用ocnprint从Cadence导出的文本数据中第一行会有留白,很多人介绍的方法都是手工删除前面的空白或空白和文字。个人觉得比较麻烦。

其实如果在MATLAB中可以用textread命令来控制从第N行开始读取数据。比如输出的原始数据,第1行是空白数据,第2行是数据名称,第3行是空白行,而有用数据是从第4行开始,如下所示:
---------------------------------------------------------------------------------------------------------------

vgs               getData("P2:id    getData("P2:gm    getData("P2:gd    getData("P2:vt 
 0.00000e+00      -2.02811e-14       2.35844e-13       1.23904e-14      -8.09185e-01  

----------------------------------------------------------------------------------------------------------------

可用下面语句读取这些数据:
d=textread('nmos_vgs.dat', '' , 'headerlines', 2)  %意味着读取数据的时候跳过前2行。
或者:
d=textread('nmos_vgs.dat', '' , 'headerlines', 3)  %意味着读取数据的时候跳过前3行。

(这说明textread命令似能自动跳过数据段落的空白行?)

顺手转一段网上的关于textread命令的详细解释。
---------------------------------------------------------------------------------------------------------------

实际应用中经常要读取txt文件,这个时候就需要用到强大的textread函数。
它的基本语法是:
[A,B,C,...] = textread(filename,format)
[A,B,C,...] = textread(filename,format,N)
其中
filename就是文件名, format就是要读取的格式,A,B,C就是从文件中读取到的数据。中括号里面变量的个数必须和format中定义的个数相同。 如果每N行相同格式的数据,可采用[A,B,C,...] = textread(filename,format,N)的语法,读取N次。
举例如下:
例子1:myfile.txt 中的内容如下:
1 2 3 4
5 6 7 8
9 10 11 12
相应的语句为:
filename = '.\myfile.txt';
[data1,data2,data3,data4]=textread(filename,'%n%n%n%n');
或者为
filename = '.\myfile.txt';
[data1,data2,data3,data4]=textread(filename,'%n%n%n%n',3);


例子2:myfile.txt 中的内容如下:
1, 2, 3, 4
5, 6, 7, 8
9, 10, 11, 12
相应的语句为:
filename = '.\myfile.txt';
[data1,data2,data3,data4]=textread(filename,'%n%n%n%n',
'
delimiter', ',');
这里
delimiter指出分隔符,读数据的时候会自动跳过分隔符。


例子3:myfile.txt 中的内容如下:
% this a comment
1, 2, 3, 4
5, 6, 7, 8
9, 10, 11, 12
相应的语句为:
filename = '.\myfile.txt';
[data1,data2,data3,data4]=textread(filename,'%n%n%n%n',
'
delimiter', ',','headerlines'1);
这里告诉textread跳过一开始的1行,1可以替换为任意你要跳过的行数。

例子4:
myfile.txt 中的内容如下:

Sally Level1 12.34 45 Yes

相应语句为:

filename = '.\myfile.txt';

[names, types, x, y, answer] = textread(filename , '%s %s %f %d %s', 1)

4.1如果要忽略12.34这个浮点数。
[names, types, y, answer] = textread(filename , '%s %s %*f %d %s', 1)

%*f 告诉textread跳过一个浮点数。
4.2 如果要忽略Level,指读取后面的数字,

[names, levelnum, x, y, answer] = textread(filename , ... '%s Level%d %f %d %s', 1)


例子5:myfile.txt 中的内容如下
1,2,3,4,,6 
7,8,9,,11,12

想用nan替代为空的部分。
语句如下:
data = textread('myfile.txt', 'delimiter', ',', 'emptyvalue', NaN);


例子6: myfile.txt 中的内容如下 
 Sally    Type1 12.34 45 Yes
 Joe      Type2 23.54 60 No
 Bill     Type1 34.90 12 No

如果只想读第一列,其余的跳过
filename = '.\myfile.txt';
[names]=textread(filename,'%s%*[^\n]')

 %[^...] - reads characters not matching characters between the brackets until first matching character 
%[^\n] 就是一直读到行尾。
%*[^\n] 就是从当前直接跳到行尾。
可对比
[names,rest]=textread(filename,'%s%[^\n]')


例子7: myfile.txt 中的内容如下                 
Location;date;discharge
Lobith;1989-01-01;00:00;2801
Lobith;1989-01-02;00:00;2619

语法如下:
[Location,... % 1
    yyyy              ,... % 2a datestring year
    mm                ,... % 2b datestring month
    dd                ,... % 2c datestring day
    HH                ,... % 3a timestring hour
    MM                ,... % 3b timestring minute
    discharge] =  ...
    textread(filename,'%s%f-%f-%f%f:%f%f',...
        'headerlines',1,...
        'delimiter'  ,';');

注意:
textread不用先fopen那个文件,适用于格式统一的txt文件的一次性大批量读取。textread读取某个文件后,下次再用textread读取这个文件时,还是会从文件头开始读取。
  

----------------------------------------------------------------------------------------------------------------

 

1

点赞

刚表态过的朋友 (1 人)

发表评论 评论 (1 个评论)

回复 Reteaa 2023-4-18 18:01
This article helps me a lot! Greatest appreciation!

facelist

您需要登录后才可以评论 登录 | 注册

  • 关注TA
  • 加好友
  • 联系TA
  • 0

    周排名
  • 0

    月排名
  • 0

    总排名
  • 0

    关注
  • 124

    粉丝
  • 43

    好友
  • 304

    获赞
  • 118

    评论
  • 22143

    访问数
关闭

站长推荐 上一条 /1 下一条


小黑屋| 手机版| 关于我们| 联系我们| 在线咨询| 隐私声明| EETOP 创芯网
( 京ICP备:10050787号 京公网安备:11010502037710 )

GMT+8, 2024-11-13 07:09 , Processed in 0.015443 second(s), 8 queries , Gzip On, Redis On.

eetop公众号 创芯大讲堂 创芯人才网
返回顶部