Matlab具体例子分析

Matlab网站用两个例子非常详尽地介绍了各个命令的基本用法，实际中，面对手头上的数据，如何选用合适的命令呢？以下结合几个示例给出一些总结，大家举一反三就可以了：

1. 纯数据（列数相同）：
源文件：

[Copy to clipboard] [-]

CODE:

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.字段名（中、英文字段均可）+数据：
源文件：

[Copy to clipboard] [-]

CODE:

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文件

源文件：

[Copy to clipboard] [-]

CODE:

Group 2   12.02.2006 Limei
Samples CH1   CH2   CH3
0.000123   0.000325 0.000378 0.000598
0.000986   0.000256 0.000245 0.000698

目标文件：

[Copy to clipboard] [-]

CODE:

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); ）。一个简单的、非通用的包含注释的读取方法如下：
-------------------------------------转 and 改(by eight)---------------------------------------------------------------------------------------

[Copy to clipboard] [-]

CODE:

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);

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

4. 注释（不含独立的数字串）+数据（列数相同）：
源文件：

[Copy to clipboard] [-]

CODE:

你好 abc
欢迎来到 我们
振动论坛
vib.hit.edu.cn
1 11 111 1111
2 22 222 2222
3 33 333 3333
4 44 444 4444
5 55 555 5555

解答：直接用 importdata 便可

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

5. 注释与数据混排：
对此当然只能自己编程，举例：

源文件：

[Copy to clipboard] [-]

CODE:

1 11 111 1111
你好
2 22 222 2222
欢迎来到
3 33 333 3333
振动论坛
4 44 444 4444
vib.hit.edu.cn
5 55 555 5555

解答：
--------------------------------------------转 and 改(by eight)-------------------------------------------------------------

[Copy to clipboard] [-]

CODE:

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);

 

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

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

6.各列数据的分离：
源文件：

[Copy to clipboard] [-]

CODE:

          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'); 便可

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