1. 程式人生 > 其它 >MATLAB學習(2)

MATLAB學習(2)

MATLAB學習第二天
  1. 向量是數字的一維陣列,矩陣是數字的二維陣列。
  2. 當引用帶有冒號的向量時,例如v(:),將列出向量的所有成分。
  3. 引用元素3到7:
    rv = [1 2 3 4 5 6 7 8 9];
    sub_rv = rv(3:7)  % 返回 3 4 5 6 7
  4. 建立一個3×3矩陣m,然後將複製該矩陣的第二行和第三行兩次以建立4×3矩陣:
    a = [ 1 2 3 ; 4 5 6; 7 8 9];
    new_mat = a([2,3,2,3],:)

    顯示如下結果:

    new_mat =
          4     5     6
          7     8     9
          4     5     6
          7
    8 9
  5. zeros()函式建立一個全零的陣列,ones()函式建立一個全1的陣列,eye()函式建立一個單位矩陣,rand()函式在(0,1)上建立一個均勻分佈的隨機數陣列,magic()函式建立一個魔術方陣。
  6. 通常,要生成多維陣列,我們首先建立一個二維陣列並將其擴充套件,還可以使用cat()函式來構建多維陣列。
  7. 陣列元素的迴圈移位:
    a = [1 2 3; 4 5 6; 7 8 9]  % the original array a
    b = circshift(a,1)         %  circular shift first dimension values down by 1
    . c = circshift(a,[1 -1]) % circular shift first dimension values % down by 1 % and second dimension values to the left % by 1.
  8. 排序陣列:
    v = [ 23 45 12 9 5 0 19 17]  % horizontal vector
    sort(v)                      % sorting v
    m = [2 6 4; 5 3 9; 2 0 1]    % two dimensional array
    sort(m, 1)                   % sorting m along the row
    
    sort(m, 2) % sorting m along the column
  9. cell函式用於創造單元格陣列:
    c = cell(2, 5);
    c = {'Red', 'Blue', 'Green', 'Yellow', 'White'; 1 2 3 4 5}

    輸出:

    c = 
    {
       [1,1] = Red
       [2,1] =  1
       [1,2] = Blue
       [2,2] =  2
       [1,3] = Green
       [2,3] =  3
       [1,4] = Yellow
       [2,4] =  4
       [1,5] = White
       [2,5] =  5
    }

    有兩種方法可以引用單元格陣列的元素-

    • 將索引括在第一個方括號()中,以引用單元格集

    • 將索引括在大括號{}中,以引用單個單元格中的資料

  10. A = [1 2 3 4; 4 5 6 7; 7 8 9 10]
    A(2:3,2:3)  %第二行和第三行以及第二列和第三列

    輸出:

    A =
          1     2     3     4
          4     5     6     7
          7     8     9    10
    ans =
          5     6
          8     9
  11. 可以通過以下兩種方式垂直組合字串-

    • 使用MATLAB連線運算子[],並用分號(;)分隔每一行。請注意,在此方法中,每行必須包含相同數量的字元。對於長度不同的字串,應根據需要使用空格字元填充。

    • 使用char函式。如果字串的長度不同,char會用尾隨空格填充較短的字串,以便每行具有相同的字元數。

  12. 可以通過以下兩種方式水平組合字串-

    • 使用MATLAB串聯運算子[],並用逗號或空格分隔輸入字串。此方法保留輸入陣列中的所有尾隨空格。

    • 使用字串串聯函式 strcat,此方法刪除輸入中的尾隨空格。

  13. cellstr函式將字元陣列轉換為字串的單元格陣列。
  14. 函式語句語法:
    function [out1,out2, ..., outN] = myfun(in1,in2,in3, ..., inN)

    從表示式建立匿名函式的語法:

    f = @(arglist)expression
  15. 資料匯入:
    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

    weeklydata.txt:

    SunDay  MonDay  TuesDay  WednesDay  ThursDay  FriDay  SaturDay
    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

    執行結果:

    SunDay
       95.0100
       73.1100
       60.6800
       48.6000
       89.1300
     
    MonDay
       76.2100
       45.6500
       41.8500
       82.1400
       44.4700
     
    TuesDay
       61.5400
       79.1900
       92.1800
       73.8200
       57.6300
    
    WednesDay
       40.5700
       93.5500
       91.6900
       41.0300
       89.3600
     
    ThursDay
       55.7900
       75.2900
       81.3200
       0.9900
       13.8900
     
    FriDay
       70.2800
       69.8700
       90.3800
       67.2200
       19.8800
    
    SaturDay
       81.5300
       74.6800
       74.5100
       93.1800
       46.6000
  16. filename = '/data/myfile.txt';
    rows = 7;
    cols = 5;
     
    %開啟檔案
    fid = fopen(filename);
     
    %讀取檔案頭,找到M(月數)
    M = fscanf(fid, '%*s %*s\n%*s %*s %*s %*s\nM=%d\n\n', 1);
     
    %讀取每一組測量資料
    for n = 1:M
       mydata(n).time = fscanf(fid, '%s', 1);
       mydata(n).month = fscanf(fid, '%s', 1);
     
       %fscanf按列順序填充陣列,
       %把結果轉置
       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
     
    %關閉檔案
    fclose(fid);

    myfile.txt:

    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

    執行結果:

    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