Matlab 入門篇

Size: px

Start display at page:

Download "Matlab 入門篇"

济闵
7 years ago
Views:

1 Matlab 入門篇 2017/7/23 林崇聖

2 Matlab- 介紹 Matlab 為一種高階的程式語言, 可用於資料分析資料視覺化數值模擬與計算等, 對於科學的理論分析有很大幫助例如 : 現在有量測物體隨時間的位置資料, 此時需要得到該物體的速度加速度, 即可以利用此軟體將位置的資料讀取後, 進行速度加速的計算, 接著將時間對位置速度加速度作圖, 可以快速的得到分析的結果並且以圖形化的方式呈現結果

3 量測的位置座標計算速度加速度, 並畫圖顯示

4 開啟程式 OR 搜尋程式 : matlab 找尋此圖案的程式開啟

6 Matlab- 運算符號可輸入運算方程式或是函數等功能 ( 後面將此處輸入的內容稱為程式碼 ) 數學運算 : 遵守四則運算 4/2+3 4/(2+3) 運算結果不同 * 刮號皆使用小刮號輸入方程式後按下 Enter 將會顯示結果

7 Matlab- 數學運算元

8 Matlab- 運算符號關係運算 : 邏輯比較時 : 正確 (True) 傳回 1 不正確 (False) 傳回 0

9 Matlab- 運算符號邏輯運算元 : 邏輯運算 :

10 Matlab- 關係與邏輯運算元

11 常用運算科學符號 10^3 = 1E3 = 1e3 = 1000 三角函數 π: pi = cosd(90) = cos(pi/2)

12 顯示格式 style

14 符號使用, 若在一行中想要輸入兩段程式碼, 則可以使用逗號隔開, 如 1+1, 2+2 ; 若執行後不想顯示結果則可以在程式碼最後加上分號, 如 1+1; 若一行內無法結束程式碼時, 接在符號後面可以達成分行的動作, 如 1+ 1

15 符號使用

16 函數使用函數為特定名稱的指令, 例如輸入 mod(3,2) 將會得到 3 2 的餘數 1

17 函數查詢 help 函數名可以執行對此函數的說明, 但是說明本身很簡短, 因此有時不易瞭解該函數的使用, 這時候建議到網路上或是官網查詢函數, 將會有詳細的說明和範例

18 函數查詢官網查詢 : 輸入函數名即可搜尋

19 函數查詢

20 開啟 New Script

21 開啟 New Script

22 檔案編輯在這裡輸入要執行的程式碼後, 可以一次全部執行

23 點擊可以存檔副檔名為.m * 無法使用中文且需要英文開頭檔案編輯

24 檔案編輯點擊 Run 可以執行程式執行也可以使用快捷鍵 Ctrl+Enter

25 符號使用 % 程式的註解可以用 % 來分隔程式碼,% 後面的部分將不會被視為程式碼的一部份, 所以不會被執行 %{%} 若想使用不只一行的註解, 則可以用 %{ 想要輸入的註解 %} 來達成

26 符號使用 %% 在編寫一整段程式碼的時候, 可以使用 %% 來分隔不同部分, 在執行程式碼的時候, 如果使用快捷鍵 Ctrl+Enter 就可以只執行被 %% 隔出來的區域, 但是使用 Run 就不會區分

27 變數的使用可以將計算的結果儲存到自己定義名稱的符號中, 但是宣告名稱的時候需要注意 : 1. 第一個字母必須是英文 2. 後續英文數字和底線 _ 可混和使用 3. 字母的大寫和小寫不同 4. 避免與其他函數的名稱相同 5. 最長只能使用 31 個字母 ( 後面忽略 )

28 變數的型態使用的變數型態, 主要有兩種 : 1. 數字型態 (int/double) 2. 文字型態 (char/string) 因此在使用上 1 和 2 不可以視為相同的變數, 在 Matlab 上直接輸入的數字, 將會成為數字的型態, 而字母則會被視為變數的存在所以想要宣告文字型態的變數, 則需要將想要使用的字母用想要儲存的字母隔開來, 如此就可以建立文字的變數

29 而執行的過程則會顯示在此清除顯示和變數使用了變數後, 其結果會顯示於此

30 清除顯示和變數因此想清除儲存的變數時, 可以使用 clear 清除顯示的結果則可以使用 clc 有時候設計程式時, 變數沒有考慮到重複執行的問題, 因此會導致執行結果出錯, 此時可以在最開頭的地方加入 clear 和 clc 確保重複執行程式時不會出錯

31 矩陣的使用想要將同一組的資訊放在同一個變數中, 這時候矩陣就是一個儲存方式儲存的方式為使用 [A,B,C;D,E,F], 其中英文字母表示不同的變數 ( 當然也可以直接輸入數字或文字 ), 而, 表示將資訊存在不同行 ( 也可以使用空白隔開 [A B]), 而 ; 則表示將資訊存在不同列, 因此上面的存放方式為 A B C D E F

32 矩陣的使用但是再儲存矩陣時需要特別注意每個文字型態變數的長度相同, 以及每個矩陣之間如果要存放時兩者之間的長度也需要相同例如說 [ str, s ] 是不能儲存成矩陣只要出現這些紅字就表示矩陣間的長度有不同導致錯誤

33 矩陣的使用變數想要呼叫矩陣中特定位置的資料時, 則是使用 X(1,1) 或是 X(1) 來呼叫第一個儲存的資料, 因此輸入兩個以上的數字時 (A,B)A 是指第幾列,B 是指第幾行只輸入一個數字時, 則是會由列的方向開始排序, 排到第幾個就是對應第幾個數字行增加的方向列增加的方向 1 (1,1) 2 (2,1) 3 (1,2) 4 (2,2) 5 (1,3) 6 (2,3)

34 矩陣的使用若想要呼叫大於一個資訊時, 則可以使用 : 來呼叫如 X(1:2,1) 表示第一列第二列的第一行元素同時被呼叫

35 矩陣的使用此外也可以利用 : 來建立矩陣, 例如說想要創立以 2 為間隔, 數字由 5 到 11 的矩陣, 可以輸入 5:2:11

36 矩陣的使用若不知道此矩陣的大小時, 但是想要呼叫最後一個元素, 則可以利用 end 來達成, 如 X(end,1) 表示第一行最後一列的元素單獨使用 : 則可以表示全部的元素

37 矩陣的使用也可以使用 size(x) 這個函數來得到 X 矩陣的大小, 顯示的結果會是 [ 列行 ] 的矩陣也可以直接 size(x,1) 來得到 X 有幾列, size(x,2) 則是有幾行

38 矩陣的使用也可以直接利用 X(3,4)=2 來建立矩陣的元素, 但是沒有定義元素的部分,Matlab 在數字形式會自動補 0, 而文字形式則會補 null, 表示這個地方沒有文字存在

39 矩陣的使用加減乘除的部分則是符合矩陣運算的方式, 因此需要注意矩陣的大小是否可以符合矩陣運算另一方面, 如果只是想要矩陣間的元素做運算, 則可以使用.* 的方式達成例如想要 X 和 Y 的每一個元素互相相乘, 則 X.*Y 即可 ( 但是 X 和 Y 的大小要相同, 否則會無法相乘 )

40 顯示特定文字若想顯示計算完, 特定格式的文字和數字在指令區, 則需要使用 fprintf fprintf( 需要顯示的格式, 變數 1, 變數 2, ) 當我們不明白計算後的結果是多少時, 沒辦法直接將數字變成文字顯示, 此時可以使用這個函數, 將計算存放的變數指定至函數中, 依照需要的文字格式將變數顯示例如現在想要輸入學生 A 的成績 :

41 顯示特定文字 fprintf 特殊符號 : 如同前面的文字格式變數, 內容表示的即是文字 % 後面接序格式, 對應變數的形式, 詳細的對照請參照下一頁 \ 後面接續特定的字母,\n 表示換行 \t 表示 Tab \r\n 則是另外一種換行所以如果需要使用 \ % 這幾個特殊符號到文字中, 則同時輸入兩次即可未使用也不影響結果

42 顯示特定文字 fprintf 特殊符號 : %s 將變數顯示文字的時候 %d 將變數顯示整數的時候 %f 將變數顯示浮點數的時候 ( 小數點 ) %g 比 %f 或 %e 顯示簡潔 ( 可能沒有差異 ) %e 科學記號表示

43 格式參考

44 建立文字的相似函數 sprintf : 與 fprintf 相似, 但是 sprintf 是建立文字格式的字串輸出, 因此可以將 sprintf 建立的結果儲存至變數中

45 函數使用 Matlab 彩蛋, 可以執行結果 : vibes teapotdemo logo fifteen xpbombs why

46 Matlab 程式流程控制篇 2017/7/7 林崇聖

47 條件選擇一般在執行的過程中, 會需要利用條件去選擇要執行的動作, 此時就需要使用條件選擇的功能第一種選擇的方法 : if( 第一個判斷的條件 )%True or False 若第一個判斷成立則執行這裡的程式碼 elseif( 第二個判斷的條件 ) 若第二個判斷成立則執行這裡的程式碼 elseif( ) %elseif 可以使用很多次 else % 不需要輸入判斷條件上述的判斷皆未成立時則執行這裡的程式碼 end % 判斷結束的地方需要加入 end * 可以只使用 if-end 或 if-elseif-end 或 if-else-end

48 條件選擇判斷的條件一旦成功, 則後面的條件不會判斷就會結束這個結構

49 條件選擇另一種選擇的方法 : switch-case-otherwise-end switch 變數 case 數值 / 文字 * 大於一種選擇時用 { 選擇 1, 選擇 2, } otherwise end 比對變數存放的數值或文字, 尋找是否有符合的 case 存在, 若都不存在時則執行 otherwise, 且與 if 相同由上而下執行, 若已找到符合的 case 則下面的 case otherwise 皆不會執行 * 與 if 相似可以只放一個 case,otherwise 並非必要

50 條件選擇

51 條件選擇條件選擇需注意存放的變數是否皆存在, 由於執行至某個符合的項目就會跳出, 後面選項中的變數可能不會被執行, 後續若還是需要使用到此變數, 就會出現變數未使用的錯誤, 因此注意每個選項的變數使用, 是否符合後續的程式執行

52 迴圈有時候我們會需要執行重複的動作, 或是每次動作與上次動作類似, 此時若需要不斷重複寫入近乎相同的程式碼, 則會使整體的程式碼變長變複雜且耗費時間, 迴圈即可處理這類行的動作例如 : 計算若是要自行輸入 1+~+50, 則會使程式碼變的很長, 而且很花時間若是要加到更多的數字, 或是減少數字, 也是非常耗費時間去修改

53 迴圈 for-end 迴圈 : for 變數名稱 = 起使數值 :( 間隔 ): 結束數值執行程式 end 如同前面建立矩陣 1:2:8 的用法, 由 1 開始計算, 以 2 為間隔, 加到不大於 8 為止所以對應迴圈會執行的數值為 : 這四個數值會對應到定義的變數中, 每次迴圈依序改變變數的數值

54 迴圈 while-end 迴圈 : while 變數執行程式 end 當變數是數值 1(True) 時,while 將會執行, 因此可用來判別當某數值發生之前仍然繼續執行

55 迴圈 * 需注意給予的迴圈條件, 有可能導致無法離開 while 迴圈, 造成程式卡住的情形, 此時只要在 Matlab 的指令區按下 Ctrl+C 即可強制停止程式執行 * 使用 break 的函數可以強制跳出迴圈,continue 則可以跳過此次迴圈

56 迴圈 - 練習不同數值之加總, 可以快速更改

57 迴圈 - 練習不同數值之加總, 可以快速更改

58 Matlab 繪圖篇 2017/7/7 林崇聖

59 繪圖指令 Matlab 有許多的繪圖指令, 這裡將只會挑選一小部分的繪圖方式來介紹使用到的指令 : plot plotyy semilogy semilogx loglog plot3

60 繪圖指令首先是最主要的 plot, 大部分的繪圖需求, 使用這個指令即可完成 plot(x 座標,y 座標, 繪圖形式, 選擇格式設定, 格式設定輸入值, ) 想要畫在 x 軸 y 軸的哪一點就分別輸入數字形式的矩陣至座標中, 若只想畫一點, 也可只輸入兩個數值接著後面的繪圖形式格式設定若無給予, 則會使用預設值, 若要設定則可以針對線的形式顏色大小等做調整

61 繪圖指令繪圖形式可分為 : -*b - 第一項表示線條形式 (line style) * 第二項表示點的形式 (marker) b 第三項表示顏色的形式 (color) 所以 - 是實線, 點用 * 的形式,b 表示使用藍色

62 其他線的形式可參考 : 繪圖指令

63 其他點的形式可參考 : 繪圖指令

64 其他顏色的形式可參考 : 繪圖指令

65 繪圖指令選擇的格式 : 若要調整線的寬度,0.5 是預設值 LineWidth, 大小如 plot(x, y, linewidth,0.5) 若要調整點的的大小,6 是預設值 : MarkerSize,6 其他可參考 :

66 繪圖指令另外可以輸入多組 x,y 座標來繪圖 :

67 繪圖指令並且可以同時調整多組的繪圖形式 :

68 繪圖指令 plotyy 則是用來建立兩個 y 軸的圖, 因此輸入的形式必須是兩組 x,y plotyy(x1,y1,x2,y2,funciton1,function2) 但是這裡要控制圖的線寬或顏色等, 則不能直接使用 plotyy 達成, 後面所輸入的是不同 y 軸的圖案類別, 所以 funcion 是指可以達成 function(x,y) 的函數, 如 plot(x,y), 這樣一來第一組 x1,y1 就是使用這個 function 來繪圖 plotyy(x1,y1,x2,y2, plot, loglog )

69 plotyy 的繪圖 : 繪圖指令

70 繪圖指令 semilogy semilogx loglog 屬於對數繪圖, 也就是座標軸的刻度是使用過 log() 處理, 刻度分布不是等比例的變化, 因此適用於數值變化極大的情況, 例如數值一開始在 10^-3 附近的變化, 可是卻突然攀升到 10^5 次方時, 要將在小範圍的變化表現出來,plot 會有困難, 這時候使用對數圖就有幫助 semilogy/x 是指對數的座標軸只有 x 軸或 y 軸 loglog 則是兩個座標軸皆為對數做圖使用的形式幾乎等同於 plot 的用法 loglog(x,y, 繪圖形式, ) semilogy(x,y, 繪圖形式, )

71 對數的繪圖 : 繪圖指令

72 繪圖指令 plot3 則是用來繪製 3D 線圖, 如果想要使用 x,y,z 的座標軸, 就需要用這個函數, 其餘的用法等同於 plot *grid on 則是開啟背景網格 ( 虛線 )

73 其他繪圖指令如果要調整的不是繪圖的線, 而是標題或是座標軸, 這時候需要用到不同的指令 xlabel: 可以給予 x 座標軸名稱 ylabel: 可以給予 y 座標軸名稱 title: 可以給予圖的名稱 legend: 可以給予圖例 axis: 可以調整座標軸的範圍

74 其他繪圖指令

75 其他繪圖指令若想要將每次使用的 plot 都繪製在同一張圖上時可以使用 hold on, 輸入 hold off 就可以關閉這個功能

76 其他繪圖指令若需要繪製兩張以上的圖在同一張影像上時, 則使用 subplot subplot( 數字 1, 數字 2, 擺放位置 ), 將會由數字 1 2 依照矩陣的切割方式, 切割出可以擺放的區域, 接著對應擺放位置的編號, 放入對應的區域中例如使用 subplot(2,3,5), 則繪圖後會被放置在 (2,3) 切割的區域中編號第五格的位置

77 其他繪圖指令

78 繪圖視窗放大縮小移動 * 不改變圖旋轉視角資料點存檔移動 * 會改變圖

79 繪圖 - 練習

80 繪圖 - 練習

81 Matlab 寫檔與讀檔篇 2017/7/7 林崇聖

82 寫讀檔指令 Matlab 同樣有許多的讀寫檔指令, 因此不會全部做介紹一般在處理資料時, 不可能自己手動不斷輸入資料, 因此資料會以檔案的方式儲存, 接著才用程式讀取檔案內的資料, 讀取以後再進行處理, 處理完畢後又需要將這些資料再次儲存, 因此就需要使用寫檔的方式建立處理完資料的檔案主要的指令有 fopen fclose fprintf fscanf fseek ftell dlmwrite dlmread

83 寫讀檔指令不論要讀取檔案或是寫入檔案, 一開始必須先指定路徑開啟 / 建立目標檔案, 才能使用後續部分的指令開啟 / 建立的指令為 fopen( 儲存的路徑, 開啟檔案的權限設定 ) 儲存路徑如絕對路徑 C:\matlab\data\Data.txt 開啟檔案的權限如 w+

84 寫讀檔指令 * 可注意若非絕對路徑, 單純輸入檔案名稱時如 Data.txt 則此時 Matlab 會搜尋程式本身已經選擇的路徑做搜尋程式已經選擇的路徑在此

85 寫讀檔指令開啟檔案的權限設定 : 如果未設定則預設為 r 常使用的有 : r 僅讀取檔案, w 僅寫入檔案 ( 若無存在檔案會自動建立檔案, 若有存在檔案則會覆蓋 ), a 與 w 只差別在不會覆蓋檔案, 而是會將新資料接著寫入檔案最後面, r+ w+ a+ 與原本的功能多加讀取 / 寫入的功能,

86 寫讀檔指令 [ 開啟檔案的編號, 錯誤訊息 ]=fopen( 檔案路徑, 權限設定 ) fopen 還有同時輸出兩個資訊的形式, 每個開啟的檔案皆會建立一個指定的數值, 後面再使用時即是這個編號若開啟 / 建立失敗時則檔案編號固定會回傳 -1, 且會說明為何失敗

87 寫讀檔指令且使用完檔案後, 必須要使用 fclose( all ) 的指令, 釋放讀取檔案的指令, 若一直開啟檔案, 則會佔用系統資源, 若使用大量的讀寫檔案時, 容易造成應用程式終止等問題因此完整的讀檔方式為 : fid=fopen( 檔案路徑, 存取權限 ); 利用 fid 使用其他指令達成檔案處理的過程 fclose( all ); * 也可針對特定開啟的檔案執行關閉 fclaose( fid )

88 寫讀檔指令寫入檔案的指令 : fprintf 這個指令除了能將特定文字顯示於指令區以外, 同樣也能寫入特定文字至檔案中, 並且使用方式與顯示於指令區時相同, 僅需要多加一個開啟檔案編號的變數輸入即可 fprintf( 開啟檔案之編號, 需要顯示的格式, 變數 1, 變數 2, ) 開啟檔案之編號, 為前面使用 fopen 後輸出的變數開啟檔案之編號 =fopen( 檔案路徑, 存取權限 )

89 寫讀檔指令 fprintf: * 若讓 fprintf 輸出時, 則會說明特定文字的大小為多少, 這裡的 5 表示 1+1=2 這串文字大小為 5 位元組 (bytes)

90 寫讀檔指令讀取檔案的指令 : fscanf( 開啟檔案之編號, 讀取資料之特定格式 ) 將檔案內容依照特定的格式讀出, 特定格式與 fprintf 使用上相同

91 寫讀檔指令 fseek ftell 指令的使用, 與寫入後的游標位置有關, 當 fprintf 寫入完畢後, 檔案若未使用 fopen 重新開啟, 也未使用 fclose( all ), 則此時游標位置將停留在檔案最後若想要繼續使用其他指定, 如 fscanf 就會發生讀不到東西的狀況 ftell 則會回傳現在的游標位置, 因此使用 fseek 移動位置後 ftell 就改變了

92 寫讀檔指令結果 =fseek( 開啟檔案之編號, 游標移動的距離, 游標起使位置 ) 若指令執行成功, 則結果會得到 0, 反之失敗時, 會得到 -1 開啟檔案之編號為 fopen 的結果游標移動的距離則是根據游標起使位置改變負的方向 -1-2 游標位置正的方向 1 2 游標起使位置則以 -1 為檔案開頭,0 為當前位置,-1 為檔案最後前面使用的 fseek(fid,0,-1) 表示將游標移動至 : fid 開啟的檔案的開頭位置且額外移動 0 的距離所以才會使游標回到檔案開頭的位置

93 寫讀檔指令讀寫檔案時需要注意的地方是特殊符號也會包含在一個字中, 所以移動幾個位置是以幾個字計算, 若存在空白 Tab( \t ) 換行 ( \n ), 都會計算為一個字 * 要特別注意換行的使用, 在 windows 的作業系統時, 是以 \r\n 兩個字來完成的這是使用 \n 換行的結果, 利用文字文件開啟後並沒有換行使用 \r\n 後, 才有真正的使之間換行

94 寫讀檔指令 dlmwrite( 檔案路徑, 資料, 間隔方式,R,C) 為利用特定間隔方式將資料儲存的指令, 可以不使用 fopen 的指令來達成寫檔案的動作 * 也因此不需要 fclose 檔案路徑為需要存檔的檔案位置資料則是將要寫入的數值或文字間隔方式則是決定以矩陣儲存的每個元素之間的間隔方式, 若未填入時則預設為, ( 逗號為間隔方式 ) R C 則是決定開始寫檔的位置要平移多少距離, 並且 R 是以列的方式移動,C 則是以行的方式移動

95 寫讀檔指令 * 若是以行區隔的數值會以間隔方式隔開, 但是列區隔的數值會用 \n 的換行方式間隔

$將換行的方式改為 windows 的格式 coffset 0 1 0 1 2 3 roffset 間隔為 Tab( \t ) * 有確實換行 -apeend$

96 寫讀檔指令若是需要使用其他的設定時, 則形式與前面將不同每種設定前需要用文字來指定 dlmwrite( 檔案路徑, 資料, 指定格式敘述, 指定格式 ) 將換行的方式改為 windows 的格式 coffset roffset 間隔為 Tab( \t ) * 有確實換行 -apeend 所以不會覆蓋掉原本的檔案

97 寫讀檔指令 dlmread( 檔案路徑, 間隔方式,R,C) 為利用特定間隔方式將資料讀入, 一樣不需 fopen 的指令

98 寫讀檔指令 dlmread( 檔案路徑, 間隔方式,[R1,C1,R2,C2]) 如果讀取的檔案也要限制結束的位置, 則要改由矩陣的形式輸入

99 迴圈條件選擇 - 作業費式數列 : 請問為第幾項? 2. 請問加到此項目時總和為多少?

100 繪圖 - 作業 1. 將下列位置座標輸入至矩陣中, 並繪製成圖 ( 每一個位置點間隔時間為 0.5 秒 ) 後來的位置原來的位置 2. 計算 = 後來的時間原來的時間並繪製成圖後來的時間原來的時間的速度 2 3. 利用相同的方式計算加速度並繪製成圖 4. 將三張圖合併至一起如下一頁的圖 * 以上皆須加入圖名圖例座標軸之說明 * 速度加速度繪圖時需取絕對值

101 繪圖 - 作業

102 讀檔寫檔 - 作業 1. 將繪圖題之位置寫入檔案中 2. 將繪圖題之速度寫入檔案中 3 將繪圖題之加速度寫入檔案中

103 進階挑戰題 - 寫檔將讀寫檔案的位置速度加速度對時間的關係寫入檔案中, 若無對應時間的資料時, 寫入 NaN, 如下一頁之結果 * 建立矩陣時, 可將要儲存的值輸入 NaN

104 進階挑戰題 - 寫檔

105 進階挑戰 -IV Curve 將各位已經儲存好的檔案 ( 使用 LabVIEW 儲存的檔案 ) 將一次完整掃描電壓的資料繪製成圖並將多次掃描電壓的圖綜合成 GIF 檔案儲存

106 進階挑戰 -IV Curve 參考程式碼 : 加入圖標說明即時繪圖且將 plot 得到的圖案轉換成可以輸出圖案檔案的資料輸出 GIF 檔

進階挑戰題 - 資料處理讀取附件的檔案 ( 請至網站上下載 ): 進階挑戰資料.txt 檔案標題時間第三項資料第三項資料 0.23 選擇第三項資料處理, 0.

107 進階挑戰題 - 資料處理讀取附件的檔案 ( 請至網站上下載 ): 進階挑戰資料.txt 檔案標題時間第三項資料第三項資料 0.23 選擇第三項資料處理, 得到新的資料後, 將新的資料寫入, 寫入時檔案標題時間保留, 並將第二項改為新資料, 第三項以後刪除時間與資料間隔使用 Tab 且檔案標題與下方資料要再增加標籤說明, 時間目標溫度, 之間間隔需要與下方資料對齊

108 進階挑戰題 - 資料處理寫入的檔案就如下面所示 : 加入標籤檔案標題第三項資料計算後結果時間保留紀錄

109 進階挑戰題 - 資料處理除了寫入檔案以外還需要繪製圖, 將 x 軸為時間變化,y 軸為量測的第二項資料, 以及經過計算後的第三項資料, 利用 subplot 繪製成上下兩張圖第二項資料需要使用 semilogy 繪製, 皆須加入圖名圖例座標說明並調整繪製範圍且繪製虛線標示特定數值第二項資料圖名 :Pressure, Start time=star_time *star_time 為進階挑戰資料內的檔案標題 x 座標 :time (s) y 座標 :Pressure (pa) 圖例 :pressure 繪製範圍 :x 軸範圍 0~ y 軸範圍 10 5 ~2 10 5

110 進階挑戰題 - 資料處理第三項資料圖名 :Thermocouple Temparature, Start time=star_time *star_time 為進階挑戰資料的檔案標題 x 座標 :time (s) y 座標 :temparature ( ) 圖例 : 目標溫度繪製範圍 :x 軸範圍 0~ y 軸範圍 -60~90 虛線繪製第二項資料 :y 軸繪製虛線第三項資料 :y 軸繪製虛線

111 進階挑戰題 - 資料處理最後必須在程式前後加上 tic toc 指令這兩個指令可以達成計時的功能, 因此最後的挑戰題會需要加上這兩個指令來得到執行完整體程式的時間如果需要建立相同值特定大小的矩陣, 可以使用 ones( 列數量, 行數量 ) 可建立數值皆 1 的矩陣 zeros( 列數量, 行數量 ) 可建立數值皆 0 的矩陣

112 結果圖 : 進階挑戰題 - 資料處理

113 若加入座標軸調整則可完成此圖進階挑戰題 - 資料處理

0 0 = 1 0 = 0 1 = = 1 1 = 0 0 = 1

0 0 = 1 0 = 0 1 = = 1 1 = 0 0 = 1 0 0 = 1 0 = 0 1 = 0 1 1 = 1 1 = 0 0 = 1 : = {0, 1} : 3 (,, ) = + (,, ) = + + (, ) = + (,,, ) = ( + )( + ) + ( + )( + ) + = + = = + + = + = ( + ) + = + ( + ) () = () ( + ) = + + = ( + )( + ) + = = + 0