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

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

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

3、27/16/2021星期六, 2008-9- 6, 13:07:20 2.1.1 MATLAB語(yǔ)言的變量與常量語(yǔ)言的變量與常量MATLAB語(yǔ)言的變量名規(guī)則由一個(gè)字母引導(dǎo)，后面可以為其他字符區(qū)分大小寫Abc ABc有效 MYvar12, MY_Var12 和 MyVar12_錯(cuò)誤的變量名 12MyVar, _MyVar12 MATLAB的保留常量eps, i, j, pi, NaN, Inf, i=sqrt(-1)lastwarn, lasterr高等應(yīng)用數(shù)學(xué)問(wèn)題的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ù)型符號(hào)型其他數(shù)據(jù)結(jié)構(gòu)高等應(yīng)用數(shù)學(xué)問(wèn)題的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ù)值位和一個(gè)符號(hào)位 至 double() 函數(shù)的轉(zhuǎn)換其他數(shù)據(jù)類型uint8，常用于圖像表示和處理，8位int8(),int16(),int32(),uint16(),uint32()高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)7/1

5、127/16/2021星期六, 2008-9- 6, 13:07:20 2.1.2.2 符號(hào)型符號(hào)型符號(hào)型，sym(A)，常用于公式推導(dǎo)常用于公式推導(dǎo)和求解析解變量聲明顯示符號(hào)變量的任何精度默認(rèn)精度：32位小數(shù)高等應(yīng)用數(shù)學(xué)問(wèn)題的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ù)的增加，程序運(yùn)行的時(shí)間變長(zhǎng)，如果位數(shù)要求過(guò)多可能無(wú)法顯示高等應(yīng)用數(shù)學(xué)問(wèn)題的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ù)：用單引號(hào)括起來(lái) 多維數(shù)組：是矩陣的直接擴(kuò)展，多個(gè)下標(biāo)單元數(shù)組：將不同類型數(shù)據(jù)集成到一個(gè)變量名下面，用 表示類與對(duì)象：可以定義重載函數(shù)，類與對(duì)象是“面向?qū)ο蟆本幊痰幕驹馗叩葢?yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)10/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.1.3 MATLAB的基本語(yǔ)句結(jié)構(gòu)直接賦值語(yǔ)句接賦值語(yǔ)句的基本結(jié)構(gòu)賦值語(yǔ)句的末尾加一個(gè)分號(hào)可以阻止顯示運(yùn)算結(jié)果保留

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

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

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

10、學(xué)信息學(xué)院Slide (of 112)18/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.2 基本數(shù)學(xué)運(yùn)算基本數(shù)學(xué)運(yùn)算矩陣的代數(shù)運(yùn)算矩陣的邏輯運(yùn)算矩陣的比較運(yùn)算解析結(jié)果的化簡(jiǎn)與變換基本數(shù)論運(yùn)算高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)19/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.2.1矩陣的代數(shù)運(yùn)算矩陣的代數(shù)運(yùn)算矩陣轉(zhuǎn)置加減法運(yùn)算矩陣乘法矩陣的除法矩陣翻轉(zhuǎn)矩陣乘方運(yùn)算點(diǎn)運(yùn)算高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)21/1127/16/2021星期六, 2008-9- 6, 13:07:20 加減法運(yùn)算加減法運(yùn)算數(shù)學(xué)表示用C語(yǔ)言編程難于實(shí)現(xiàn)，同A*BMATLAB語(yǔ)法注意：任一個(gè)變量可以為標(biāo)量如果矩陣維數(shù)不匹配，系統(tǒng)會(huì)報(bào)錯(cuò)高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)22/1127/16/2021星期六, 2008-9- 6, 13:07:20 矩陣乘法矩陣乘法數(shù)學(xué)表示：MA

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

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

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

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

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

18、:20 例例 2.9運(yùn)用各種函數(shù)，對(duì)下面的數(shù)據(jù)進(jìn)行取整運(yùn)算-0.2765,0.5772,1.4597,2.1091,1.191,-1.6187MATLAB代碼：高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)38/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.10假設(shè) 的Hilbert矩陣可以由 定義，試對(duì)其進(jìn)行有理數(shù)變換。結(jié)果：高等應(yīng)用數(shù)學(xué)問(wèn)題的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ǔ)句：高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)41/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.3 MATLAB語(yǔ)言的流程結(jié)構(gòu)語(yǔ)言的流程結(jié)構(gòu)循環(huán)結(jié)構(gòu)轉(zhuǎn)移結(jié)構(gòu)開(kāi)關(guān)結(jié)構(gòu)試探結(jié)構(gòu)高等應(yīng)用數(shù)學(xué)問(wèn)題的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 為一個(gè)向量，循環(huán)變量 每次從 v 向量中取一個(gè)數(shù)值，執(zhí)行一次循環(huán)體的內(nèi)容，如此下去，直至執(zhí)行完 v 向量中所有的分量高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)44/1127/16/2021星期六, 2008-9- 6, 13:07:20 例

21、例 2.13用循環(huán)求解更簡(jiǎn)單的語(yǔ)句高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的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ùn)算方法高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)48/1127/16/2021星期六, 2008-9- 6, 13:07:20 高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)50/1127/16/

23、2021星期六, 2008-9- 6, 13:07:20 2.3.3 2.3.3 開(kāi)關(guān)結(jié)構(gòu)開(kāi)關(guān)結(jié)構(gòu)高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)51/1127/16/2021星期六, 2008-9- 6, 13:07:20 高等應(yīng)用數(shù)學(xué)問(wèn)題的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)點(diǎn)：設(shè)置錯(cuò)誤陷阱提高算法的執(zhí)行效率高等應(yīng)用數(shù)學(xué)問(wèn)題的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 語(yǔ)言函數(shù)的基本結(jié)構(gòu)可變輸入輸出個(gè)數(shù)的處理inline函數(shù)與匿名函數(shù)高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)54/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.4.1 MATLAB語(yǔ)言函數(shù)的語(yǔ)言函數(shù)的基本結(jié)構(gòu)基本結(jié)構(gòu) nargin, nargout, varargin, varargout高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)55/1127/16/2021星期六, 2008-9- 6

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

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

27、8-9- 6, 13:07:20 例例 2.20利用函數(shù)的遞歸調(diào)用，計(jì)算階乘遞歸調(diào)用函數(shù)高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)61/1127/16/2021星期六, 2008-9- 6, 13:07:20 計(jì)算11！其他函數(shù)高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide

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

29、008-9- 6, 13:07:20 例例 2.22conv()可以計(jì)算兩個(gè)多項(xiàng)式的積，時(shí)使用varargin實(shí)現(xiàn)任意多個(gè)多項(xiàng)式的積MATLAB代碼高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的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)換成二進(jìn)制代碼偽代碼語(yǔ)句用戶一定要在安全的位置保留.m源文件，不能輕易刪除，因?yàn)?p文件是不可逆的高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)68/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5 二維圖形繪制二維圖形繪制二維圖形繪制基本語(yǔ)句多縱軸曲線繪制其他二維圖形繪制語(yǔ)句隱函數(shù)繪制及應(yīng)用圖形修飾數(shù)據(jù)文件的讀取與存儲(chǔ)高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLA

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

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

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

34、繪制有時(shí)兩條曲線的幅值相差甚多，需要調(diào)用雙縱軸繪制函數(shù)plotyy()例2-25：試?yán)L制曲線三、四縱軸圖形可以下載相應(yīng)函數(shù)繪制plotyyy()、plot4y()，從MathWorks File Exchange下載，還可以使用plotxx()函數(shù)高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)77/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.3 其他二維圖形繪制語(yǔ)句其他二維圖形繪制語(yǔ)句高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的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ǔ)句：高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)80/1127/16/2021星期六, 2008-9- 6, 13:07:20 下述函數(shù)可以用于將一個(gè)繪圖窗口分割成幾個(gè)不同的區(qū)域n和m分別是行數(shù)和列數(shù)， k則指示區(qū)域的序號(hào)高等應(yīng)用數(shù)學(xué)問(wèn)題的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ù)繪圖語(yǔ)句默認(rèn)區(qū)域是-2p,2p其他語(yǔ)法高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)82/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.28試?yán)L制隱函數(shù)MATLAB語(yǔ)句擴(kuò)大范圍高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)83/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.5 圖形

37、修飾圖形修飾直接采用工具欄文字修飾 特殊符號(hào)表下標(biāo)分別用 和 _ 表示 優(yōu)越性MATLAB 7.0 的新功能高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)84/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.5.6 數(shù)據(jù)文件的讀取與存儲(chǔ)數(shù)據(jù)文件的讀取與存儲(chǔ)可以采用save和load命令存儲(chǔ)和讀取數(shù)據(jù)MATLAB和Excel交互數(shù)據(jù)寫文件 xlswrite()高等應(yīng)用數(shù)學(xué)問(wèn)題的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行存儲(chǔ)數(shù)據(jù)B列存儲(chǔ)年份，C列存儲(chǔ)人口數(shù)先讀入MATLAB再繪圖高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的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é)問(wèn)題的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é)問(wèn)題的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é)問(wèn)題的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é)問(wèn)題的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é)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)93/1127/16/2021星期六, 2008-9- 6, 13:07:20 變步長(zhǎng)方法高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)94/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例 2.33繪制如下分段函數(shù)的三維圖像分段函數(shù)求值互斥的不等式、點(diǎn)運(yùn)算循環(huán)，但麻煩高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)96/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.6.3 等高線繪制等高線繪制各種等高線高等應(yīng)用數(shù)學(xué)問(wèn)題的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é)問(wèn)題的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é)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)99/1127/16/2021星期六, 2008-9- 6, 13:07:20 例例2-35已知三元隱函數(shù)繪圖語(yǔ)句高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)100/1127/16/2021星期六, 2008-9- 6, 13:07:20 2.6.5 三維圖形視角設(shè)置三維圖形視角設(shè)置兩種方法改變圖形的視角直接采用工具欄命令語(yǔ)句view() 定義為方位角， 定義為仰角高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)101/1127/16/2021星期六, 2008-9- 6, 13:07:20 視點(diǎn)定義高等應(yīng)用數(shù)學(xué)問(wèn)題的MATLAB求解東北大學(xué)信息學(xué)院Slide (of 112)102/1127/16/2021星期六, 2008-9-