第2章 MATLAB語言程序設(shè)計基礎(chǔ)

1、高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)1/1127/16/2021星期六, 2008-9- 6, 13:07:20 第第2章章 MATLAB語言程序設(shè)計基礎(chǔ)語言程序設(shè)計基礎(chǔ)MATLAB 程序設(shè)計語言基礎(chǔ)基本數(shù)學(xué)運算MATLAB語言的流程結(jié)構(gòu)函數(shù)的編寫與調(diào)試二維圖形繪制三維圖形繪制四維圖形繪制高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)2/1127/16/2021星期六, 2008-9- 6, 13:07:20 MATLAB語言有如下的優(yōu)點MATLAB 語言的簡潔高效性MATLAB 語言的科學(xué)運算功能MATLAB 語言的繪

2、圖功能MATLAB 龐大的工具箱與模塊集MATLAB 強大的動態(tài)系統(tǒng)仿真功能MATLAB 語言是當(dāng)前國際上自動控制領(lǐng)域的首選計算機語言，也是很多理工科專業(yè)最適合的計算機數(shù)學(xué)語言，本課程以 MATLAB 為主要計算機語言。高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)3/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.1 MATLAB 程序設(shè)計語言基礎(chǔ)程序設(shè)計語言基礎(chǔ)MATLAB語言的變量與常量數(shù)據(jù)結(jié)構(gòu)MATLAB的基本語句結(jié)構(gòu)冒號表達式與子矩陣提取高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)4/11

3、27/16/2021星期六, 2008-9- 6, 13:07:20 2.1.1 MATLAB語言的變量與常量語言的變量與常量MATLAB語言的變量名規(guī)則由一個字母引導(dǎo)，后面可以為其他字符區(qū)分大小寫Abc ABc有效 MYvar12, MY_Var12 和 MyVar12_錯誤的變量名 12MyVar, _MyVar12 MATLAB的保留常量eps, i, j, pi, NaN, Inf, i=sqrt(-1)lastwarn, lasterr高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)5/1127/16/2021星期六, 2008-9- 6, 13:07:

4、20 2.1.2 數(shù)據(jù)結(jié)構(gòu)數(shù)據(jù)結(jié)構(gòu)數(shù)值型數(shù)型符號型其他數(shù)據(jù)結(jié)構(gòu)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)6/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.1.2.1 2.1.2.1 數(shù)值型數(shù)據(jù)數(shù)值型數(shù)據(jù)雙精度數(shù)值變量IEEE標(biāo)準(zhǔn)，64位 (8字節(jié))，11指數(shù)位，53數(shù)值位和一個符號位 至 double() 函數(shù)的轉(zhuǎn)換其他數(shù)據(jù)類型uint8，常用于圖像表示和處理，8位int8(),int16(),int32(),uint16(),uint32()高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)7/1

5、127/16/2021星期六, 2008-9- 6, 13:07:20 2.1.2.2 符號型符號型符號型，sym(A)，常用于公式推導(dǎo)常用于公式推導(dǎo)和求解析解變量聲明顯示符號變量的任何精度默認(rèn)精度：32位小數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)8/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.1求出 p 的 300 位有效數(shù)字MATLAB代碼可以“任意”增加有效數(shù)字的位數(shù)隨著有效數(shù)字位數(shù)的增加，程序運行的時間變長，如果位數(shù)要求過多可能無法顯示高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 1

6、12)9/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.1.2.3 2.1.2.3 其他數(shù)據(jù)結(jié)構(gòu)其他數(shù)據(jù)結(jié)構(gòu)字符串型數(shù)據(jù)：用單引號括起來 多維數(shù)組：是矩陣的直接擴展，多個下標(biāo)單元數(shù)組：將不同類型數(shù)據(jù)集成到一個變量名下面，用 表示類與對象：可以定義重載函數(shù)，類與對象是“面向?qū)ο蟆本幊痰幕驹馗叩葢?yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)10/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.1.3 MATLAB的基本語句結(jié)構(gòu)直接賦值語句接賦值語句的基本結(jié)構(gòu)賦值語句的末尾加一個分號可以阻止顯示運算結(jié)果保留

7、變量： ans存放最近一次無賦值變量語句的運算結(jié)果高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)11/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.2表示矩陣MATLAB代碼其他命令高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)12/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.3試輸入復(fù)數(shù)矩陣MATLAB代碼需要避免的語句高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)13/1127/16/2021星期六, 2008-9-

8、6, 13:07:20 函數(shù)調(diào)用語句函數(shù)調(diào)用舉例函數(shù)可以通過不同的方式被調(diào)用內(nèi)核函數(shù)， *.m函數(shù)匿名函數(shù)、 inline函數(shù)（不建議使用）重載函數(shù)、私有函數(shù)等高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)14/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.1.4 冒號表達式與子矩陣提取冒號表達式與子矩陣提取冒號表達式是定義行向量的有效方法開始于值 ，步長為值 ，終止于值 .默認(rèn)步長值為1高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)15/1127/16/2021星期六, 2008-9- 6, 13

9、:07:20 例例 2.4用不同的步距生 間的向量高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)16/1127/16/2021星期六, 2008-9- 6, 13:07:20 子矩陣提取的基本語句格式 表示子矩陣要保留的行號構(gòu)成的向量 表示要保留的列號構(gòu)成的向量:，表示要提取所有行或列，取決于其位置子矩陣提取子矩陣提取高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)17/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.5矩陣A的不同子矩陣提取MATLAB代碼高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大

10、學(xué)信息學(xué)院Slide (of 112)18/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.2 基本數(shù)學(xué)運算基本數(shù)學(xué)運算矩陣的代數(shù)運算矩陣的邏輯運算矩陣的比較運算解析結(jié)果的化簡與變換基本數(shù)論運算高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)19/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.2.1矩陣的代數(shù)運算矩陣的代數(shù)運算矩陣轉(zhuǎn)置加減法運算矩陣乘法矩陣的除法矩陣翻轉(zhuǎn)矩陣乘方運算點運算高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)20/1127/16/2021星期六,

11、2008-9- 6, 13:07:20 矩陣轉(zhuǎn)置矩陣轉(zhuǎn)置矩陣表示矩陣 A，n 行m列，被稱作 矩陣Hermite轉(zhuǎn)置一般轉(zhuǎn)置高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)21/1127/16/2021星期六, 2008-9- 6, 13:07:20 加減法運算加減法運算數(shù)學(xué)表示用C語言編程難于實現(xiàn)，同A*BMATLAB語法注意：任一個變量可以為標(biāo)量如果矩陣維數(shù)不匹配，系統(tǒng)會報錯高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)22/1127/16/2021星期六, 2008-9- 6, 13:07:20 矩陣乘法矩陣乘法數(shù)學(xué)表示：MA

12、TLAB表示：注意：系統(tǒng)自動檢測矩陣維數(shù)是否匹配高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)23/1127/16/2021星期六, 2008-9- 6, 13:07:20 矩陣的除法矩陣的除法矩陣左除：求解線性方程組：MATLAB解法：最小二乘解若A為非奇異方陣，則高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)24/1127/16/2021星期六, 2008-9- 6, 13:07:20 矩陣右除：求解線性方程組：MATLAB解法：最小二乘解若A為非奇異方陣，則等效的運算高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slid

13、e (of 112)25/1127/16/2021星期六, 2008-9- 6, 13:07:20 矩陣翻轉(zhuǎn)矩陣翻轉(zhuǎn)左右翻轉(zhuǎn) 上下翻轉(zhuǎn)旋轉(zhuǎn) 90o如何旋轉(zhuǎn)180o？高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)26/1127/16/2021星期六, 2008-9- 6, 13:07:20 矩陣乘方矩陣乘方A 為方陣，求矩陣 A 的 x 次冪數(shù)學(xué)描述MATLAB命令x為整數(shù)x為非整數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)27/1127/16/2021星期六, 2008-9- 6, 13:07:20 點運算點運算矩陣對應(yīng)元素的

14、直接運算例如： 為 另一個例子： 為高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)28/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.6求出矩陣A的全部三次方根，并檢驗結(jié)果MATLAB代碼另兩個根高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)29/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.2.2 矩陣的邏輯運算矩陣的邏輯運算邏輯變量：當(dāng)前版本有邏輯變量非 0 表示邏輯 1邏輯運算（相應(yīng)元素間的運算）與運算或運算非運算異或運算高等應(yīng)用數(shù)學(xué)問題的MATLAB求解

15、東北大學(xué)信息學(xué)院Slide (of 112)30/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.2.3 2.2.3 矩陣的比較運算矩陣的比較運算各種允許的比較關(guān)系, =, , =, = =,=, find(), all(), any()實例高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)31/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.2.4 解析結(jié)果的化簡與變換解析結(jié)果的化簡與變換函數(shù) simple() 用于數(shù)學(xué)公式的化簡其他常用化簡函數(shù)numden(), collect(), expand(),

16、 factor()高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)32/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.7化簡多項式用不同的函數(shù)求解高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)33/1127/16/2021星期六, 2008-9- 6, 13:07:20 變量替換變量替換變量替換語句該函數(shù)執(zhí)行點運算轉(zhuǎn)換成LaTeX表示高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)34/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例

17、2.8試用 對下式進行雙線性變換，得出其LaTeX排版格式MATLAB代碼高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)35/1127/16/2021星期六, 2008-9- 6, 13:07:20 結(jié)果：LaTeX排版顯示高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)36/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.2.5 基本數(shù)論運算基本數(shù)論運算高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)37/1127/16/2021星期六, 2008-9- 6, 13:07

18、:20 例例 2.9運用各種函數(shù)，對下面的數(shù)據(jù)進行取整運算-0.2765,0.5772,1.4597,2.1091,1.191,-1.6187MATLAB代碼：高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)38/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.10假設(shè) 的Hilbert矩陣可以由 定義，試對其進行有理數(shù)變換。結(jié)果：高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)39/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.11試求1856120和1483

19、720的最大公約數(shù)與最小公倍數(shù)，并求出所得出的最小公倍數(shù)的質(zhì)因數(shù)分解。MATLAB求解語句：高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)40/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.12試求出11000間的全部質(zhì)數(shù)得到的所有質(zhì)數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)41/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.3 MATLAB語言的流程結(jié)構(gòu)語言的流程結(jié)構(gòu)循環(huán)結(jié)構(gòu)轉(zhuǎn)移結(jié)構(gòu)開關(guān)結(jié)構(gòu)試探結(jié)構(gòu)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Sli

20、de (of 112)42/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.3.1 循環(huán)結(jié)構(gòu)循環(huán)結(jié)構(gòu)for循環(huán)結(jié)構(gòu)v 為一個向量，循環(huán)變量 每次從 v 向量中取一個數(shù)值，執(zhí)行一次循環(huán)體的內(nèi)容，如此下去，直至執(zhí)行完 v 向量中所有的分量高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)43/1127/16/2021星期六, 2008-9- 6, 13:07:20 while循環(huán)結(jié)構(gòu)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)44/1127/16/2021星期六, 2008-9- 6, 13:07:20 例

21、例 2.13用循環(huán)求解更簡單的語句高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)45/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.14用循環(huán)求解求最小的 m，使下式成立 使用while結(jié)構(gòu) 這里不宜應(yīng)用for結(jié)構(gòu) 高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)46/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.15求如下值使用MATLAB的循環(huán)結(jié)構(gòu)和向量化運算方法高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)47/1127

22、/16/2021星期六, 2008-9- 6, 13:07:20 2.3.2 轉(zhuǎn)移結(jié)構(gòu)轉(zhuǎn)移結(jié)構(gòu)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)48/1127/16/2021星期六, 2008-9- 6, 13:07:20 高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)49/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.16用for循環(huán)求解最大的 m，使下式成立MATLAB代碼遠(yuǎn)比while結(jié)構(gòu)繁瑣高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)50/1127/16/

23、2021星期六, 2008-9- 6, 13:07:20 2.3.3 2.3.3 開關(guān)結(jié)構(gòu)開關(guān)結(jié)構(gòu)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)51/1127/16/2021星期六, 2008-9- 6, 13:07:20 高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)52/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.3.4 2.3.4 試探結(jié)構(gòu)試探結(jié)構(gòu)全新結(jié)構(gòu)優(yōu)點：設(shè)置錯誤陷阱提高算法的執(zhí)行效率高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)53/1127/16/2

24、021星期六, 2008-9- 6, 13:07:20 2.4 函數(shù)編寫與調(diào)試函數(shù)編寫與調(diào)試MATLAB 語言函數(shù)的基本結(jié)構(gòu)可變輸入輸出個數(shù)的處理inline函數(shù)與匿名函數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)54/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.4.1 MATLAB語言函數(shù)的語言函數(shù)的基本結(jié)構(gòu)基本結(jié)構(gòu) nargin, nargout, varargin, varargout高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)55/1127/16/2021星期六, 2008-9- 6

25、, 13:07:20 例例 2.17為什么需要函數(shù)？問題M-腳本以M-文件的形式被存取若最大值不為 10000， M-文件就需要被修改則需要一種新的文件形式高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)56/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.18按例2.17的要求編寫一個MATLAB函數(shù)M-文件如下舉例優(yōu)點：無需修改程序高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)57/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.19編寫一個函數(shù)生成 Hil

26、bert矩陣要求若只給出一個輸入?yún)?shù)，則自動生成一個方陣在函數(shù)中給出合適的幫助信息檢測輸入和返回變量的個數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)58/1127/16/2021星期六, 2008-9- 6, 13:07:20 MATLAB函數(shù)如下高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)59/1127/16/2021星期六, 2008-9- 6, 13:07:20 在線幫助命令生成Hilbert矩陣高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)60/1127/16/2021星期六, 200

27、8-9- 6, 13:07:20 例例 2.20利用函數(shù)的遞歸調(diào)用，計算階乘遞歸調(diào)用函數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)61/1127/16/2021星期六, 2008-9- 6, 13:07:20 計算11！其他函數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)62/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.21試比較遞歸算法和循環(huán)算法在 Fibonacci 數(shù)列中應(yīng)用的優(yōu)劣Fibonacci數(shù)列MATLAB代碼高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide

28、(of 112)63/1127/16/2021星期六, 2008-9- 6, 13:07:20 第25項用 for循環(huán)結(jié)構(gòu)計算第100項注意：本例不宜采用遞歸結(jié)構(gòu)符號運算高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)64/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.4.2 可變輸入輸出可變輸入輸出個數(shù)的處理個數(shù)的處理vararginvarargoutvarargin1, varargin2, vararginn高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)65/1127/16/2021星期六, 2

29、008-9- 6, 13:07:20 例例 2.22conv()可以計算兩個多項式的積，時使用varargin實現(xiàn)任意多個多項式的積MATLAB代碼高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)66/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.4.3 inline 函數(shù)與匿名函數(shù)函數(shù)與匿名函數(shù)inline 函數(shù)（功能重疊，不建議使用）可以免去編輯M-文件MATLAB7.0以上版本提供的匿名函數(shù)在MATLAB工作空間內(nèi)已經(jīng)定義的變量可以直接使用高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)67/11

30、27/16/2021星期六, 2008-9- 6, 13:07:20 2.4.4 偽代碼與代碼保密處理偽代碼與代碼保密處理偽代碼化的目的提高程序的執(zhí)行速度保密：把ASCII的 .m文件轉(zhuǎn)換成二進制代碼偽代碼語句用戶一定要在安全的位置保留.m源文件，不能輕易刪除，因為.p文件是不可逆的高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)68/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5 二維圖形繪制二維圖形繪制二維圖形繪制基本語句多縱軸曲線繪制其他二維圖形繪制語句隱函數(shù)繪制及應(yīng)用圖形修飾數(shù)據(jù)文件的讀取與存儲高等應(yīng)用數(shù)學(xué)問題的MATLA

31、B求解東北大學(xué)信息學(xué)院Slide (of 112)69/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.1 二維圖形繪制基本語句二維圖形繪制基本語句兩個序列 和 構(gòu)造向量:利用這些數(shù)據(jù)繪圖高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)70/1127/16/2021星期六, 2008-9- 6, 13:07:20 其他方法 仍為向量，而 為矩陣，亦即 和 均為矩陣，且 和 的行和列數(shù)均相同假設(shè)有多對這樣的向量或矩陣高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)71/1127/16/2021星期六,

32、2008-9- 6, 13:07:20 改變曲線性質(zhì)選項列表高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)72/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.23繪制函數(shù)MATLAB代碼問題：如何檢查曲線？高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)73/1127/16/2021星期六, 2008-9- 6, 13:07:20 改用不同的步長，觀察曲線是否相同換用較小步長在 附近用變步長算法高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)74/1127/16/2

33、021星期六, 2008-9- 6, 13:07:20 例例 2.24繪制飽和函數(shù)方程MATLAB繪圖語句更簡單的命令高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)75/1127/16/2021星期六, 2008-9- 6, 13:07:20 圖像修飾圖像修飾每一個窗口、曲線和坐標(biāo)軸都是一個對象，對象的屬性可以通過函數(shù)set()來設(shè)置，也可以通過函數(shù)get()來獲取。高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)76/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.2 多縱軸曲線的繪制多縱軸曲線的

34、繪制有時兩條曲線的幅值相差甚多，需要調(diào)用雙縱軸繪制函數(shù)plotyy()例2-25：試?yán)L制曲線三、四縱軸圖形可以下載相應(yīng)函數(shù)繪制plotyyy()、plot4y()，從MathWorks File Exchange下載，還可以使用plotxx()函數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)77/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.3 其他二維圖形繪制語句其他二維圖形繪制語句高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)78/1127/16/2021星期六, 2008-9- 6, 1

35、3:07:20 例例 2.26繪制極坐標(biāo)函數(shù) 和繪制曲線高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)79/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.27以正弦數(shù)據(jù)為例，試在同一窗口的不同區(qū)域用不同的繪圖方式繪制出相應(yīng)的曲線MATLAB求解語句：高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)80/1127/16/2021星期六, 2008-9- 6, 13:07:20 下述函數(shù)可以用于將一個繪圖窗口分割成幾個不同的區(qū)域n和m分別是行數(shù)和列數(shù)， k則指示區(qū)域的序號高等應(yīng)用數(shù)學(xué)問題的MATLA

36、B求解東北大學(xué)信息學(xué)院Slide (of 112)81/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.4 隱函數(shù)繪制及應(yīng)用隱函數(shù)繪制及應(yīng)用隱函數(shù)繪圖語句默認(rèn)區(qū)域是-2p,2p其他語法高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)82/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.28試?yán)L制隱函數(shù)MATLAB語句擴大范圍高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)83/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.5 圖形

37、修飾圖形修飾直接采用工具欄文字修飾 特殊符號表下標(biāo)分別用 和 _ 表示 優(yōu)越性MATLAB 7.0 的新功能高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)84/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.6 數(shù)據(jù)文件的讀取與存儲數(shù)據(jù)文件的讀取與存儲可以采用save和load命令存儲和讀取數(shù)據(jù)MATLAB和Excel交互數(shù)據(jù)寫文件 xlswrite()高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)85/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2-29已知E

38、xcel文件 census.xls給出某省人口數(shù)第5-67行存儲數(shù)據(jù)B列存儲年份，C列存儲人口數(shù)先讀入MATLAB再繪圖高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)86/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.6 三維圖形表示三維圖形表示三維曲線繪制三維曲面繪制三維圖形視角設(shè)置三維隱函數(shù)圖形繪制三維圖形的旋轉(zhuǎn)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)87/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.6.1 三維曲線繪制三維曲線繪制三維曲線繪制stem3,

39、 fill3, bar3等高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)88/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.30試?yán)L制參數(shù)方程其中，MATLAB繪圖高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)89/1127/16/2021星期六, 2008-9- 6, 13:07:20 使用stem3() 函數(shù)圖形窗口的工具欄3D繪圖和視角變換讀取坐標(biāo)值局部放大高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)90/1127/16/2021星期六, 2008-9-

40、 6, 13:07:20 2.6.2 三維曲面繪制三維曲面繪制一般曲面繪制其他函數(shù)surfl(), surfc()等高線繪制 contour(), contours()高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)91/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.31給出二元函數(shù)如下，繪制3D圖像MATLAB繪圖表面圖高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)92/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.32試?yán)L制出二元函數(shù)繪制3D圖像高等應(yīng)用

41、數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)93/1127/16/2021星期六, 2008-9- 6, 13:07:20 變步長方法高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)94/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.33繪制如下分段函數(shù)的三維圖像分段函數(shù)求值互斥的不等式、點運算循環(huán)，但麻煩高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)95/1127/16/2021星期六, 2008-9- 6, 13:07:20 分段函數(shù)的數(shù)學(xué)形式MATLAB繪圖高等應(yīng)

42、用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)96/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.6.3 等高線繪制等高線繪制各種等高線高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)97/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例2-34例2-33的分段函數(shù)，繪制各種等高線三維等高線高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)98/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.6.4 三維隱函數(shù)圖繪制三維

43、隱函數(shù)圖繪制下載ezimplot3()函數(shù)fun可以為隱函數(shù)字符串匿名函數(shù)M函數(shù)高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)99/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例2-35已知三元隱函數(shù)繪圖語句高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)100/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.6.5 三維圖形視角設(shè)置三維圖形視角設(shè)置兩種方法改變圖形的視角直接采用工具欄命令語句view() 定義為方位角， 定義為仰角高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)101/1127/16/2021星期六, 2008-9- 6, 13:07:20 視點定義高等應(yīng)用數(shù)學(xué)問題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)102/1127/16/2021星期六, 2008-9-