Matlab数据存取(二):各类常见数据文件应用实例

为了方便学习,部分内容略有修改,侵删。转载,来自@声振之家公众号。


Matlab网站用两个例子非常详尽地介绍了各个命令的基本用法,实际中,面对手头上的数据,如何选用合适的命令呢?

以下结合几个示例给出一些总结,大家举一反三就可以了:

1. 纯数据(列数相同)

源文件:

0 3866.162 2198.938 141.140
1 3741.139 2208.475 141.252
2 3866.200 2198.936 141.156
3 3678.048 2199.191 141.230
4 3685.453 2213.726 141.261
5 3728.769 2212.433 141.277
6 3738.785 2214.381 141.256
7 3728.759 2214.261 141.228
8 3748.886 2214.299 141.243
9 3748.935 2212.417 141.253
10 3733.612 2226.653 141.236
11 3733.583 2229.248 141.223
12 3729.229 2229.118 141.186

对于这个txt文件,由于各行列数相同,故简单地使用load,importdata均可。

2. 字段名(中、英文字段均可)+数据

源文件:

CH0 CH1 CH2 CH3
0.000123 0.000325 0.000378 0.000598
0.000986 0.000256 0.000245 0.000698

由于是记录的形式,因此各行列数必相同(缺少部分列时请自行在文件中补上 Inf 或 NaN),故直接使用 importdata 便可。

3. 注释(含有独立的数字串)+数据(列数相同)

这个文件有4列,但前6行是文字说明,4列数字是从第8行开始的。

现在我想把这个文件的前2列文字说明提出来组成一个新的dat文件

源文件:

Group 2 12.02.2006 Limei
Samples of datas: 50000
CH0 CH1 CH2 CH3
0.000123 0.000325 0.000378 0.000598
0.000986 0.000256 0.000245 0.000698

目标文件:

Group 2 12.02.2006 Limei
Samples of datas: 50000
CH0 CH1
0.000123 0.000325
0.000986 0.000256

由于注释中含有独立的数字串,且注释部分没有明显的格式,这时候用importdata, load等高级命令直接读取会失败,用textread,dlmwrite 等格式化命令也不太合适,因此只能使用低级命令进行读取。(当然了,可以跳过注释部分直接用高级命令读取数据,即:[a b c d] = textread(filename,'%f %f %f %f','headerlines',4); )。

一个简单的、非通用的包含注释的读取方法如下

clc;clear;
fid = fopen('exp.txt', 'r');
fid_n=fopen('ex.dat','w');
while ~feof(fid)
    tline=fgetl(fid);
    if ~isempty(tline)
        if double(tline(1))>=48 && double(tline(1))<=57 
            %数值开始
            a=strread(tline);
            a(3:4)=[];
            fprintf(fid_n,'%f %f\n',a);
            clear a;
        elseif double(tline(1))==67 
            %字母C开始
            [b1,b2,b3,b4]=strread(tline,'%s %s %s %s');
            b=[b1{1},' ',b2{1}];
            fprintf(fid_n,'%s\n',b);
            clear b b1 b2 b3 b4;
        else
            fprintf(fid_n,'%s\n',tline);
        end
    else
        fprintf(fid_n,'%s\n',tline);
    end
end
fclose(fid);
fclose(fid_n);
fclose(fid);
fclose(fid_n);

4. 注释(不含独立的数字串)+数据(列数相同)

源文件:

你好 Evryone
欢迎来到 我们
CSDN
1 11 111 1111
2 22 222 2222
3 33 333 3333
4 44 444 4444
5 55 555 5555

直接用importdata即可。

注:有时候注释中含有独立的数字串也可以 importdata 成功,不过得到的结果有可能不正确,建议这时候使用第3种情形的读取方式

5. 注释与数据混排

只能自己编程,举例:

源文件:

1 11 111 1111
你好
2 22 222 2222
欢迎来到
3 33 333 3333
CSDN
4 44 444 4444
大家庭
5 55 555 5555

编程:

function [data]=distilldata(infile)
    %功能说明:
    %将保存数据的原始文件中的数值数据读入到一个data变量中
    %使用说明:
    % infile——原始数据文件名;
    % data=数据变量
    tmpfile='tmp2.mat';
    fidin=fopen(infile,'r');
    %打开原始数据文件(.list)
    fidtmp=fopen(tmpfile,'w'); 
    %创建保存数据文件(不含说明文字)
    while ~feof(fidin) 
        %判断是否为文件末尾
        tline=fgetl(fidin); 
        %从文件读入一行文本(不含回车键)
        if ~isempty(tline) 
            %判断是否空行
            [m,n]=size(tline);
            flag=1;
            for i=1:n 
                %判断一行中有没有字符(+-.Ee和空格键除外)
                if ~(tline(i)==' '|tline(i)=='-'|tline(i)=='.'|tline(i)=='E'...
                    |tline(i)=='e'|tline(i)=='+'...
                    |(double(tline(i))>=48&&double(tline(i))<=57))
                    flag=0;
                    break;
                end
             end
            if flag==1 
                %如果是数字行,把此行数据写入文件
                fprintf(fidtmp,'%s\n',tline);
            end
        end
    end
    fclose(fidin);
    fclose(fidtmp);
    data=textread(tmpfile);
    delete(tmpfile);
end

如果要求不高,也可以使用 textread 函数跳过注释部分进行读取,不过前提是需要事先知道文件内容的结构(即哪行是数据、哪行是注释)

6.各列数据的分离

源文件:

0 + 47038.7 1.05 09:26:07 C
2 + 46477.7 1.03 09:28:38 C
4 + 44865.7 1.04 09:28:48 C
6 + 41786.4 1.03 09:28:56 C
8 + 39896.0 0.97 09:29:03 C
10 + 37518.4 0.93 09:29:15 C
12 + 35858.5 0.92 09:29:30 C
14 + 46105.0 1.03 09:30:21 C
16 + 46168.6 6.89 09:30:30 C
18 + 48672.3 4.33 09:30:40 C
20 + 49565.7 0.49 09:30:48 C
22 + 49580.7 0.53 09:30:55 C
24 + 49602.3 0.84 09:31:03 C
26 + 49582.5 1.51 09:31:11 C
28 + 49577.0 1.39 09:31:19 C
30 + 49589.3 0.61 09:31:27 C
32 + 49578.3 1.06 09:31:29 C
34 + 49512.5 1.77 09:31:38 C

直接用 [a,b,c,d,e,f]=textread(yourfilename,'%d %c %f %f %s %c');

 

 

 

 

 

 

 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值