黃聰明 國立臺灣師範大學數學系 min@math.ntnu.edu.tw http://math.ntnu.edu.tw/~min MATLAB 基本功能介紹 黃聰明 國立臺灣師範大學數學系 min@math.ntnu.edu.tw http://math.ntnu.edu.tw/~min MATLAB 基本功能介紹 T.-M.Huang T.-M.Huang
2-1 變數和陣列 變數 變數名稱 常用變數型態 double char 使用者設定名稱的陣列,在實體上是由一塊記憶體區域所組成。 大小寫字母有別 不可超過63個字元 第一個字必須是文字,其後可用文字、數字及底線任意組合。 佔64位元 擁有15到16個有效位數 可處理實數、虛數和複數 內定 i 和 j 代表虛數 Var = 10 + 10i ; MATLAB 基本功能介紹 T.-M.Huang
2-2 MATLAB 變數的初始化 利用宣告的方式,指定資料給變數 三種用來初始化變數的方式 從鍵盤輸入資料給變數 從檔案讀取資料 var = 40i; var2 = var/5; array = [ 1 2 3 4 ]; x = 1; y = 2; var = expression MATLAB 基本功能介紹 T.-M.Huang
特殊符號(I) 分號(;) 冒號(:) 轉置transpose( ’) first:incr:last MATLAB 基本功能介紹 T.-M.Huang
特殊符號(II) 百分比符號(%) 程式中加入註解 (Comments) >> y = (5*2+3.5)/5; % 將運算結果儲存在變數 y,但不用顯示於螢幕 >> z = y^2 % 將運算結果儲存在變數 z,並顯示於螢幕 z = 7.2900 MATLAB 基本功能介紹 T.-M.Huang
特殊用途矩陣 指令 說明 zeros(m, n) 產生維度為 m×n ,構成元素全為 0 的矩陣 ones(m, n) eye(n) 產生維度為 n×n ,對角線的各元素全為 1 ,其他各元素全為 0 的單位矩陣 pascal(m, n) 產生維度為 m×n 的 Pascal 矩陣 vander(m, n) 產生維度為 m×n 的 Vandermonde 矩陣 hilb(n) 產生維度為 n×n 的 Hilbert 矩陣 rand(m, n) 產生 [0, 1] 均勻分佈的亂數矩陣,其維度為 m×n randn(m, n) 產生 µ = 0, σ= 1 的正規分佈亂數矩陣,其維度為 m×n magic(n) 產生維度為 n×n 的魔方陣,其各個直行、橫列及兩對角線的元素和都相等 MATLAB 基本功能介紹 T.-M.Huang
從鍵盤輸入初始化變數 浮點數資料 字元字串型態 30 T.M. Huang MATLAB 基本功能介紹 T.-M.Huang
2-3 多維陣列與子陣列 在 MATLAB 的資料型態中,向量可視為一維陣列,矩陣可視二維陣列,對於維度(Dimensions)超過 1 的陣列則均可視為「多維陣列」(Multidimesional Arrays,簡稱 N-D Arrays)。 MATLAB 基本功能介紹 T.-M.Huang
二維陣列(I) 二維陣列又簡稱矩陣,具有兩個維度(Dimensions) 「橫列」(Row) 「直行」(Column) 直行 橫列 (1,1) (1,2) (1,3) (1,4) (2,1) (2,2) (2,3) (2,4) (3,1) (3,2) (3,3,) (3,4) 橫列 MATLAB 基本功能介紹 T.-M.Huang
二維陣列(II) 二維陣列,可對應至一個 X - Y 二維平面座標,圖示如下: Y(行) X(列) MATLAB 基本功能介紹 T.-M.Huang
矩陣的索引或下標 (I) 所有矩陣的內部 A( i , j ) 表示法都是以直行 為主的一維向量 第i橫列、第j直行的元素 A(i,j)和A(i+(j-1)*m) 是完全一樣的。 ★ m為矩陣A的列數 i 與 j 即是此元素的 下標(Subscript) 或索引(Index) 可以使用一維或 二維下標來存取 矩陣 MATLAB 基本功能介紹 T.-M.Huang
矩陣的索引或下標 (II) A= A(1:5,5) A(:,5) A(21;25) A(1:5,end) A(:,end) 4 1 10 6 11 16 2 21 8 7 9 12 17 22 3 5 13 18 23 14 19 24 15 20 25 A(2,3) A(12) A= A(1:5,end) A(:,end) A(21;end) A(4:5,2:3) A([ 9 14; 10 15 ]) MATLAB 基本功能介紹 T.-M.Huang
三維陣列(I) 將兩個二維陣列疊在一起,就形成第三個維度,此第三個維度稱為「頁」(Page),圖示如下: 行 頁 列 列 頁 行 (1,1) (1,2) (1,3) (1,4) (2,1) (2,2) (2,3) (2,4) (3,1) (3,2) (3,3) (3,4) 列 (1,1) (1,2) (1,3) (1,4) (2,1) (2,2) (2,3) (2,4) (3,1) (3,2) (3,3) (3,4) 列 頁 行 MATLAB 基本功能介紹 T.-M.Huang
三維陣列(II) 三維陣列,可對應至一個 X - Y - Z 三維立體座標,圖示如下: Z(頁) Y(行) X(列) MATLAB 基本功能介紹 T.-M.Huang
三維陣列(III) 三維陣列定址,可以(列,行,頁)定之。 以維度為 3×4×2 的三維陣列為例,其定址方式可圖示如下: 陣列 A 是三維陣列,其中 A(:,:,1)代表第一頁的二維陣列,A(:,:,2)代表第二頁的二維陣列。 MATLAB 基本功能介紹 T.-M.Huang
2-5 特殊的數值 函式 目的 pi 代表 到15位有效數字的值 i , j 代表 的值 inf 代表無窮大,通常是除以0的結果 代表 到15位有效數字的值 i , j 代表 的值 inf 代表無窮大,通常是除以0的結果 Nan `不是數字’,是由未定義的數學運算而來 eps 電腦上兩個數字間的最小差異 ans 用來儲存一個敘述式的結果 MATLAB 基本功能介紹 T.-M.Huang
2-6 顯示輸出資料 指令格式 結 果 舉 例 format short 顯示4位小數(預設值) 12.3457 format long 結 果 舉 例 format short 顯示4位小數(預設值) 12.3457 format long 顯示14位小數 12.345678901234567 format short e 顯示5個數字加冪次方 1.2346e+001 format short g 總共顯示5個數字(可加或不加冪次方) 12.346 format long e 顯示15位小數字加冪次方 1.234567890123457e+001 format long g 總共顯示15個數字(可加或不加冪次方) 12.3456789012346 format hex 16位元進位格式 4028b0fcd32f707a MATLAB 基本功能介紹 T.-M.Huang
disp num2str int2str 轉換數字成字串 轉換整數成字串 >> str = [‘ The value of pi = ’ , num2str(pi)]; >> disp(str); The value of pi = 3.1416 轉換數字成字串 int2str >> FileName1 = 'rslt_w'; >> FileName = strcat(FileName1, int2str(10)) FileName = rslt_w10 轉換整數成字串 MATLAB 基本功能介紹 T.-M.Huang
fprintf格式化輸出 >> fprintf('The value of pi is %6.2f \n', pi) fprintf(format, data) 一個或多個陣列變數 字串:描述輸出資料的方式 >> fprintf('The value of pi is %f \n', pi) The value of pi is 3.141593 >> fprintf('The value of pi is %6.2f \n', pi) The value of pi is 3.14 MATLAB 基本功能介紹 T.-M.Huang
Example/conv_rate/conv_ratio.m 格式字串 說 明 %d 以指數格式顯示數值 %e %f 以浮點數格式顯示數值 %g 以浮點數或指數格式顯示數值,由何者較短為優先顯示 \n 跳到新的一行 MATLAB 基本功能介紹 T.-M.Huang
save filename var1 var2 var3 2-7 資料檔案 save 把工作區的資料 存進一個磁碟 檔案中 save filename var1 var2 var3 檔案延伸檔名為“.mat” load 把磁碟檔案中 的資料存進 工作區 load filename MATLAB 基本功能介紹 T.-M.Huang
variable_name = expression; 2-8 純量與陣列運算 variable_name = expression; ii=ii+1; 計算等號右邊 敘述式的結果 並將其儲存於 左邊變數內 運算方法 MATLAB 形式 加法 a + b 減法 a - b 乘法 a * b 除法 a / b 冪次方 a^b MATLAB 基本功能介紹 T.-M.Huang
陣列與矩陣運算 矩陣的加減與一般純量(Scalar)的加減類似 相加或相減的矩陣必需具有相同的維度 矩陣與純量可以直接進行加減,MATLAB 會直接將加減應用到每一個元素 >> A = [12 34 56 20]; >> B = [1 3 2 4]; >> C = A + B C = 13 37 58 24 >> A = [1 2 3 2 1] + 5 A = 6 7 8 7 6 MATLAB 基本功能介紹 T.-M.Huang
矩陣的乘法與除法 純量對矩陣的乘或除,可比照一般寫法 欲進行矩陣相乘,必需確認第一個矩陣的直行數目( Column Dimension) 必須等於第二個矩陣的橫列數目(Row Dimension) 矩陣的除法,常藉由反矩陣或解線性方程式來達成 >> A = [123 , 442]; >> C = A/3 >> B = 2*A C = B = 41.0000 147.3333 246 884 >> A = [1; 2]; >> B = [3, 4, 5]; >> C = A * B C = 3 4 5 6 8 10 MATLAB 基本功能介紹 T.-M.Huang
矩陣的左、右除法 A\b:矩陣除法由 inv(A)*b 來定義,其中inv(A)為A的反矩陣 a\B:矩陣除法由 a*inv(B) 來定義 >> y = b' / A y = -0.0333 0.4667 -0.0333 >> A = magic(3) A = 8 1 6 3 5 7 4 9 2 >> b = [1; 2; 3]; >> e = A * x - b e = 1.0e-015 * -0.4441 >> x = A \ b x = 0.0500 0.3000 >> f = y * A - b' f = 1.0e-015 * 0 0.4441 0 MATLAB 基本功能介紹 T.-M.Huang
矩陣的次方運算 矩陣的次方運算,可由「^」來達成,但矩陣必需是方陣,其次方運算才有意義 在「*」,「/」及「^」之前加上一個句點,MATLAB 將會執行矩陣內「元素對元素」(Element-by-element) 的運算 >> A = magic(3); >> B = A^2 B = 91 67 67 67 91 67 67 67 91 >> A = [12; 45]; >> B = [2; 3]; >> C = A.*B C = 24 135 >> C = A. * B ??? C = A. * B | Error: Unexpected MATLAB operator. >> E = A.^2 E = 144 2025 >> D = A./B D = 6 15 *前後不能有任何空格 MATLAB 基本功能介紹 T.-M.Huang
轉置和「共軛轉置」矩陣 複數矩陣 z,其「共軛轉置」矩陣 (Conjugate Transpose) 可表示成 矩陣 z′ 想得到任何矩陣 z 的轉置(Transpose),則可表示成矩陣 z.′ 若 z 為實數,則 z′和 z.′的結果是一 樣的 MATLAB 基本功能介紹 T.-M.Huang
2-11 繪圖功能簡介 ylabel('y'); xlabel('x'); title('Plot of y = x.^2-10.*x+15') x = 0:0.1:10; y = x.^2–10.*x+15; plot(x,y); ylabel('y'); grid on; grid off; xlabel('x'); MATLAB 基本功能介紹 T.-M.Huang
print <options> <filename> 功能敘述 -deps 產生一個灰階的 eps 圖形 -depsc 產生一個彩色的 eps 圖形 -dpng 產生一個 png 圖形 -dtiff 產生一個壓縮的 tiff 圖形 x = 0:0.1:10; y = x.^2–10.*x+15; plot(x,y); print –djpeg fig_ch2.jpeg 產生一個 jpeg 圖形, 檔案名稱fig_ch2.jpeg MATLAB 基本功能介紹 T.-M.Huang
圖形化影像輸出(II) MATLAB 基本功能介紹 T.-M.Huang
圖形化影像輸出(II) MATLAB 基本功能介紹 T.-M.Huang
圖形化影像輸出(II) MATLAB 基本功能介紹 T.-M.Huang
圖形化影像輸出(II) MATLAB 基本功能介紹 T.-M.Huang
多重線條繪圖 x = 0:pi/100:2*pi; y1 = sin(2*x); y2 = 2*cos(2*x); plot( x, y1, x, y2 ); plot(x, y1); hold on; plot(x, y2); hold off MATLAB 基本功能介紹 T.-M.Huang
線條顏色、形式、 資料標記及說明文字(I) plot(x, y1); hold on; plot(x, y2); hold off plot(x, y1,'ro-.' ); hold on; plot(x, y2,'gx:'); hold off MATLAB 基本功能介紹 T.-M.Huang
線條顏色、形式、 資料標記及說明文字(I) legend('sin(2x)','2cos(2x)','Location','NorthWest'); text(3.3, 1.5, ‘2cos(2x)') text(1.5, 0.5, 'sin(2x)') MATLAB 基本功能介紹 T.-M.Huang
線條顏色、形式、 資料標記及說明文字(II) x = 0:pi/100:2*pi; y1 = sin(2*x); y2 = 2*cos(2*x); plot(x, y1,'ro-.' ); hold on; plot(x, y2,'gx:'); hold off legend('sin(2x)','2cos(2x)','Location','NorthWest'); text(1.5, 0.5, 'sin(2x)') text(3.3, 1.5, '2cos(2x)') set(gca,'xtick',[0 1/2*pi pi 3/2*pi 2*pi]) set(gca,'xticklabel',{'0','1/2 pi','pi','3/2 pi','2 pi'}) MATLAB 基本功能介紹 T.-M.Huang