MATLAB學習(2)
阿新 • • 發佈:2021-11-18
MATLAB學習第二天
- 向量是數字的一維陣列,矩陣是數字的二維陣列。
- 當引用帶有冒號的向量時,例如v(:),將列出向量的所有成分。
- 引用元素3到7:
rv = [1 2 3 4 5 6 7 8 9]; sub_rv = rv(3:7) % 返回 3 4 5 6 7
- 建立一個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
- zeros()函式建立一個全零的陣列,ones()函式建立一個全1的陣列,eye()函式建立一個單位矩陣,rand()函式在(0,1)上建立一個均勻分佈的隨機數陣列,magic()函式建立一個魔術方陣。
- 通常,要生成多維陣列,我們首先建立一個二維陣列並將其擴充套件,還可以使用cat()函式來構建多維陣列。
- 陣列元素的迴圈移位:
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
- 排序陣列:
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
- 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 }
有兩種方法可以引用單元格陣列的元素-
-
將索引括在第一個方括號()中,以引用單元格集
-
將索引括在大括號{}中,以引用單個單元格中的資料
-
-
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
-
可以通過以下兩種方式垂直組合字串-
-
使用MATLAB連線運算子[],並用分號(;)分隔每一行。請注意,在此方法中,每行必須包含相同數量的字元。對於長度不同的字串,應根據需要使用空格字元填充。
-
使用char函式。如果字串的長度不同,char會用尾隨空格填充較短的字串,以便每行具有相同的字元數。
-
-
可以通過以下兩種方式水平組合字串-
-
使用MATLAB串聯運算子[],並用逗號或空格分隔輸入字串。此方法保留輸入陣列中的所有尾隨空格。
-
使用字串串聯函式 strcat,此方法刪除輸入中的尾隨空格。
-
- cellstr函式將字元陣列轉換為字串的單元格陣列。
- 函式語句語法:
function [out1,out2, ..., outN] = myfun(in1,in2,in3, ..., inN)
從表示式建立匿名函式的語法:
f = @(arglist)expression
- 資料匯入:
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
-
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