Matlab數據導入

在MATLAB中導入數據意味著從外部檔加載數據。importdata函數允許加載不同格式的各種數據檔。它有以下五種形式 -

默認情況下，Octave不支持importdata()函數，因此必須搜索並安裝此軟體包，以使以下示例與Octave安裝配合使用。

示例 - 1

下麵加載並顯示一個圖像檔。 創建腳本檔並在其中鍵入以下代碼：

filename = 'mydog.jpg';
A = importdata(filename);
image(A);

執行上面示例代碼，得到以下結果 -

運行檔時，MATLAB顯示圖像檔。但是，必須將其存儲在當前工作目錄(F:\worksp\matlab)中。

示例 - 2

在這個例子中，導入一個文本檔並使用指定分隔符號和列標題。下麵用列標題創建一個空格分隔的ASCII檔，名稱為weeklydata.txt。

文本檔weeklydata.txt的內容如下 -

週一  週二  週三  週四  週五  週六  周日
95.01   76.21   61.54    40.57       55.79    70.28   81.53
73.11   45.65   79.19    93.55       75.29    69.87   74.68
60.68   41.85   92.18    91.69       81.32    90.38   74.51
48.60   82.14   73.82    41.03       0.99     67.22   93.18
89.13   44.47   57.63    89.36       13.89    19.88   46.60

創建腳本檔並在其中鍵入以下代碼 -

Trial>> filename = 'weeklydata.txt';
delimiterIn = ' ';
headerlinesIn = 1;
A = importdata(filename,delimiterIn,headerlinesIn);
% View data
for k = [1:7]
   disp(A.colheaders{1, k})
   disp(A.data(:, k))
   disp(' ')
end
週一
   95.0100
   73.1100
   60.6800
   48.6000
   89.1300


週二
   76.2100
   45.6500
   41.8500
   82.1400
   44.4700


週三
   61.5400
   79.1900
   92.1800
   73.8200
   57.6300


週四
   40.5700
   93.5500
   91.6900
   41.0300
   89.3600


週五
   55.7900
   75.2900
   81.3200
    0.9900
   13.8900


週六
   70.2800
   69.8700
   90.3800
   67.2200
   19.8800


周日
   81.5300
   74.6800
   74.5100
   93.1800
   46.6000

示例 - 3

在這個例子中，從剪貼板中導入數據。

將以下行複製到剪貼板 -

Matlab is simple，I like it.

創建腳本檔並鍵入以下代碼 -

A = importdata('-pastespecial')

執行上面示例代碼，得到以下結果 -

Trial>> A = importdata('-pastespecial')

A =

  1×1 cell 數組

    {'Matlab is simple，I like it.'}

低級檔I/O

importdata函數是一個高級函數。 MATLAB中的低級檔I/O函數允許對檔讀取或寫入數據控制。但是，這些函數需要更多關於檔的詳細資訊才能有效地工作。

MATLAB為位元組或字元級別的讀寫操作提供以下函數 -

導入具有低級別I/O的文本數據檔

MATLAB提供以下函數，用於低級導入文本數據檔 -

• fscanf函數讀取文本或ASCII檔中的格式化數據。
• fgetl和fgets函數一次讀取檔的一行，換行符分隔每行。
• 讀取函數以位元組或位級別讀取數據流。

示例

假設有一個文本數據檔myfile.txt保存在工作目錄中。此檔存放三個月的降雨資料： 2012年6月，7月和8月降雨量。

myfile.txt中的數據在五個地方包含重複的時間，月份和降雨量測量。標題數據存儲月數M; 所以有M個集合測量值。

檔內容看起來如下 -

Rainfall Data
Months: June, July, August

M=3
12:00:00
June-2012
17.21  28.52  39.78  16.55 23.67
19.15  0.35   17.57  NaN   12.01
17.92  28.49  17.40  17.06 11.09
9.59   9.33   NaN    0.31  0.23
10.46  13.17  NaN    14.89 19.33
20.97  19.50  17.65  14.45 14.00
18.23  10.34  17.95  16.46 19.34
09:10:02
July-2012
12.76  16.94  14.38  11.86 16.89
20.46  23.17  NaN    24.89 19.33
30.97  49.50  47.65  24.45 34.00
18.23  30.34  27.95  16.46 19.34
30.46  33.17  NaN    34.89  29.33
30.97  49.50  47.65  24.45 34.00
28.67  30.34  27.95  36.46 29.34
15:03:40
August-2012
17.09  16.55  19.59  17.25 19.22
17.54  11.45  13.48  22.55 24.01
NaN    21.19  25.85  25.05 27.21
26.79  24.98  12.23  16.99 18.67
17.54  11.45  13.48  22.55 24.01
NaN    21.19  25.85  25.05 27.21
26.79  24.98  12.23  16.99 18.67

從此檔導入數據並顯示此數據。請按照以下步驟 -

• 使用fopen函數打開檔並獲取檔識別字。
• 使用格式說明符描述檔中的數據，例如字串的%s，整數%d或浮點數%f。
• 要跳過檔中的文字，請將其包含在格式說明中。要跳過數據字段，請在說明符中使用星號(*)。例如，要讀取頭檔並返回M的單個值，可以書寫為：

  M = fscanf(fid, '%*s %*s\n%*s %*s %*s %*s\nM=%d\n\n', 1);
  

• 默認情況下，fscanf會根據格式描述讀取數據，直到找不到數據的匹配，或到達檔的末尾。這裏使用for迴圈讀取3組數據，每次讀取7行和5列。
• 將在工作空間中創建一個名為mydata的結構來存儲從檔讀取的數據。這個結構有三個字段 - 時間，月份和雨天數組。

創建腳本檔(readraindata.m)並在其中鍵入以下代碼 -

filename = 'myfile.txt';
rows = 7;
cols = 5;

% open the file
fid = fopen(filename);

% read the file headers, find M (number of months)
M = fscanf(fid, '%*s %*s\n%*s %*s %*s %*s\nM=%d\n\n', 1);

% read each set of measurements
for n = 1:M
   mydata(n).time = fscanf(fid, '%s', 1);
   mydata(n).month = fscanf(fid, '%s', 1);

   % fscanf fills the array in column order,
   % so transpose the results
   mydata(n).raindata  = ...
      fscanf(fid, '%f', [rows, cols]);
end
for n = 1:M
   disp(mydata(n).time), disp(mydata(n).month)
   disp(mydata(n).raindata)
end

% close the file
fclose(fid);

執行上面示例代碼，得到以下結果 -

Trial>> readraindata
12:00:00
June-2012
   17.2100   17.5700   11.0900   13.1700   14.4500
   28.5200       NaN    9.5900       NaN   14.0000
   39.7800   12.0100    9.3300   14.8900   18.2300
   16.5500   17.9200       NaN   19.3300   10.3400
   23.6700   28.4900    0.3100   20.9700   17.9500
   19.1500   17.4000    0.2300   19.5000   16.4600
    0.3500   17.0600   10.4600   17.6500   19.3400

09:10:02
July-2012
   12.7600       NaN   34.0000   33.1700   24.4500
   16.9400   24.8900   18.2300       NaN   34.0000
   14.3800   19.3300   30.3400   34.8900   28.6700
   11.8600   30.9700   27.9500   29.3300   30.3400
   16.8900   49.5000   16.4600   30.9700   27.9500
   20.4600   47.6500   19.3400   49.5000   36.4600
   23.1700   24.4500   30.4600   47.6500   29.3400

15:03:40
August-2012
   17.0900   13.4800   27.2100   11.4500   25.0500
   16.5500   22.5500   26.7900   13.4800   27.2100
   19.5900   24.0100   24.9800   22.5500   26.7900
   17.2500       NaN   12.2300   24.0100   24.9800
   19.2200   21.1900   16.9900       NaN   12.2300
   17.5400   25.8500   18.6700   21.1900   16.9900
   11.4500   25.0500   17.5400   25.8500   18.6700