matlab在科學(xué)計(jì)算中的應(yīng)用ppt課件

1、ppt課件.1第二章 MATLAB 語(yǔ)言程序設(shè)計(jì)基礎(chǔ) MATLAB 語(yǔ)言的簡(jiǎn)潔高效性 MATLAB 語(yǔ)言的科學(xué)運(yùn)算功能 MATLAB 語(yǔ)言的繪圖功能 MATLAB 龐大的工具箱與模塊集 MATLAB 強(qiáng)大的動(dòng)態(tài)系統(tǒng)仿真功能ppt課件.2本章主要內(nèi)容 MATLAB 基本命令簡(jiǎn)介 MATLAB 程序設(shè)計(jì)語(yǔ)言基礎(chǔ) 基本數(shù)學(xué)運(yùn)算 MATLAB語(yǔ)言流程控制 MATLAB 函數(shù)的編寫 二維圖形繪制 三維圖形繪制ppt課件.32.1 MATLAB 基本命令簡(jiǎn)介 MATLAB 通用命令ppt課件.4 Matlab常用的窗口幫助命令常用的窗口幫助命令ppt課件.5 help 命令直接輸入 help, MATLA

2、B 將列出所有的幫助主題，每個(gè)幫助主題對(duì)應(yīng)于 MATLAB搜索路徑中的一個(gè)目錄；help 后加幫助主題，可獲得指定幫助主題的幫助信息；help 后加函數(shù)名； help 后加命令名，將得到指定命令的用法；ppt課件.6 demo 命令ppt課件.7 helpwin 命令用于打開(kāi) MATLAB 的幫助文件窗ppt課件.8 lookfor 命令允許用戶通過(guò)完整的或部分關(guān)鍵字來(lái)搜索要查找的內(nèi)容 who 和 whos 的作用是列出在 MATLAB 工作內(nèi)存中駐留的變量名 exist 命令用來(lái)查找或檢查變量和函數(shù)的存在性ppt課件.92.2 MATLAB 程序設(shè)計(jì)語(yǔ)言基礎(chǔ) MATLAB 語(yǔ)言的變量命名規(guī)則

3、命名規(guī)則是：（1）變量名必須是不含空格的單個(gè)詞；（2）變量名區(qū)分大小寫；（3）變量名最多不超過(guò)19個(gè)字符；（4）變量名必須以字母打頭，之后可以是 任意字母、數(shù)字或下劃線，變量名中 不允許使用標(biāo)點(diǎn)符號(hào)ppt課件.10MATLAB 的保留常量特殊變量 取 值 ans 用于結(jié)果的缺省變量名 pi 圓周率 eps 計(jì)算機(jī)的最小數(shù)，當(dāng)和 1 相加就產(chǎn)生一個(gè)比1大的數(shù) flops 浮點(diǎn)運(yùn)算數(shù) inf 無(wú)窮大，如1/0 NaN 不定量，如0/0 i，j i=j=1 nargin 所用函數(shù)的輸入變量數(shù)目 nargout 所用函數(shù)的輸出變量數(shù)目 realmin 最小可用正實(shí)數(shù) realmax 最大可用正實(shí)數(shù) p

4、pt課件.11數(shù)學(xué)運(yùn)算符號(hào)及標(biāo)點(diǎn)符號(hào)數(shù)學(xué)運(yùn)算符號(hào)及標(biāo)點(diǎn)符號(hào)+ 加法運(yùn)算，適用于兩個(gè)數(shù)或兩個(gè)同階矩陣相加. 減法運(yùn)算 * 乘法運(yùn)算 .* 點(diǎn)乘運(yùn)算 / 除法運(yùn)算 ./ 點(diǎn)除運(yùn)算 乘冪運(yùn)算 . 點(diǎn)乘冪運(yùn)算 反斜杠表示左除. ppt課件.12（1）MATLAB的每條命令后，若為逗號(hào)或無(wú)標(biāo)逗號(hào)或無(wú)標(biāo)點(diǎn)點(diǎn)符號(hào)，則顯示命令的結(jié)果；若命令后為分號(hào)分號(hào)，則禁止顯示結(jié)果. （2）“%” 后面所有文字為注釋. （3） “.”表示續(xù)行.數(shù)學(xué)運(yùn)算符號(hào)及標(biāo)點(diǎn)符號(hào)數(shù)學(xué)運(yùn)算符號(hào)及標(biāo)點(diǎn)符號(hào)ppt課件.13 雙精度數(shù)值變量 IEEE標(biāo)準(zhǔn)，64位 (占8字節(jié))，11指數(shù)位，53數(shù)值位和一個(gè)符號(hào)位 double( ) 函數(shù)的轉(zhuǎn)換

5、其他數(shù)據(jù)類型 uint8( )，無(wú)符號(hào)8位整形數(shù)據(jù)類型，值域?yàn)?至255，常用于圖像表示和處理。（節(jié)省存儲(chǔ)空間，提高處理速度） int8( ), int16( ), int32( ),uint16( ), uint32( )ppt課件.14 符號(hào)型，sym(A), 常用于公式推導(dǎo)、解析解解法 符號(hào)變量聲明 syms var_list var_props 例：syms a b real syms c positiveppt課件.15 符號(hào)型數(shù)值可采用變精度函數(shù)求值 vpa(A), 或 vap(A,n) vpa(pi)ans = 3.1415926535897932384626433832795

6、vpa(pi,60)ans =3.14159265358979323846264338327950288419716939937510582097494ppt課件.16 字符串型數(shù)據(jù)：用單引號(hào)括起來(lái) 。 多維數(shù)組：是矩陣的直接擴(kuò)展，多個(gè)下標(biāo)。 單元數(shù)組：將不同類型數(shù)據(jù)集成到一個(gè)變量名下面，用表示；例：用Ai,j可表示單元數(shù)組A的第i行,第j列的內(nèi)容。 類與對(duì)象：允許用戶自己編寫包含各種復(fù)雜詳細(xì)的變量，可以定義傳遞函數(shù)。MATLAB支持的其它數(shù)據(jù)結(jié)構(gòu)ppt課件.17 直接賦值語(yǔ)句 賦值變量賦值表達(dá)式 例： a=pi2 a = 9.8696 例：表示矩陣 B=1+9i,2+8i,3+7j;4+6j

7、 5+5i,6+4i;7+3i,8+2j 1iB = 1.0000 + 9.0000i 2.0000 + 8.0000i 3.0000 + 7.0000i 4.0000 + 6.0000i 5.0000 + 5.0000i 6.0000 + 4.0000i 7.0000 + 3.0000i 8.0000 + 2.0000i 0 + 1.0000iMATLAB 的基本語(yǔ)句結(jié)構(gòu)ppt課件.18 函數(shù)調(diào)用語(yǔ)句返回變量列表函數(shù)名（輸入變量列表） 例：a,b,c=my_fun(d,e,f,c) 冒號(hào)表達(dá)式 v=s1:s2:s3 該函數(shù)生成一個(gè)行向量v，其中s1是起始值， s2是步長(zhǎng)（若省略步長(zhǎng)為1），

8、s3是最大值。 例：用不同的步距生成 (0,p) 間向量。 v1=0:0.2:piv1 = Columns 1 through 9 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 Columns 10 through 16 1.8000 2.0000 2.2000 2.4000 2.6000 2.8000 3.0000ppt課件.19 v2=0:-0.1:pi %步距為負(fù)，不能生成向量，得出空矩陣v2 = Empty matrix: 1-by-0 v3=0:piv3 = 0 1 2 3 v4=pi:-1:0 逆序排列構(gòu)成新向量

9、v4 = 3.1416 2.1416 1.1416 0.1416 v5=0:0.4:pi,piv5 = 0 0.4000 0.8000 1.2000 1.6000 2.0000 2.4000 2.8000 3.1416ppt課件.20 基本語(yǔ)句格式 B=A(v1,v2) v1、 v2分別表示提取行（列）號(hào)構(gòu)成的向量。例： A=1,2,3,4;3,4,5,6;5,6,7,8;7,8,9,0A = 1 2 3 4 3 4 5 6 5 6 7 8 7 8 9 0 B1=A(1:2:end,:) 提取全部奇數(shù)行、所有列。B1 = 1 2 3 4 5 6 7 8子矩陣提取ppt課件.21 B2=A(3,

10、2,1,2,3,4) 提取3，2，1行、2，3，4列構(gòu)成子矩陣。 A =B2 = 1 2 3 4 6 7 8 3 4 5 6 4 5 6 5 6 7 8 2 3 4 7 8 9 0 B3=A(:,end:-1:1) 將A矩陣左右翻轉(zhuǎn)，即最后一列排在最前面。B3 = 4 3 2 1 6 5 4 3 8 7 6 5 0 9 8 7ppt課件.22 矩陣表示 矩陣轉(zhuǎn)置 數(shù)學(xué)表示 (若A有復(fù)數(shù)元素，先轉(zhuǎn)置再取各元素共軛復(fù)數(shù)值，Hermit轉(zhuǎn)置） MATLAB 求解 BA. C=A2.3 基本數(shù)學(xué)運(yùn)算矩陣的代數(shù)運(yùn)算ppt課件.23 矩陣加減法 C=A+B D=A-B 注意維數(shù)是否相等 注意其一為標(biāo)量的情

11、形 矩陣乘法 數(shù)學(xué)表示 MATLAB 表示 C=A*B 注意兩個(gè)矩陣相容性 ppt課件.24 矩陣除法 矩陣左除：AX = B，求 X MATLAB 求解：X=AB 若A為非奇異方陣，則 X=A-1B 最小二乘解（若A不是方陣） 矩陣右除：XA = B，求 X MATLAB求解：X=B/A 若A為非奇異方陣，則 X=BA-1 最小二乘解（若A不是方陣）ppt課件.25 矩陣翻轉(zhuǎn) 左右翻轉(zhuǎn) B=fliplr(A) 上下翻轉(zhuǎn) C=flipud(A) 旋轉(zhuǎn) 90o (逆時(shí)針） D=rot90(A) 如何旋轉(zhuǎn)180o？ D=rot180(A) ? Undefined function or varia

12、ble rot180. D=rot90(rot90(A)矩陣乘方 A 為方陣，求 MATLAB 實(shí)現(xiàn)： F=Ax,1iji mjba 1,ijni jca ppt課件.26 點(diǎn)運(yùn)算-矩陣對(duì)應(yīng)元素的直接運(yùn)算數(shù)學(xué)表示 : MATLAB 實(shí)現(xiàn)： C=A.*B例: A=1,2,3;4,5,6;7,8,0; B=A.AB = 1 4 27 256 3125 46656 823543 16777216 1 C=A.*AC = 1 4 9 16 25 36 49 64 0ppt課件.27 邏輯變量： 當(dāng)前版本有邏輯變量 對(duì) double 變量來(lái)說(shuō)，非 0 表示邏輯 1 邏輯運(yùn)算（相應(yīng)元素間的運(yùn)算） 與運(yùn)算

13、A&C 或運(yùn)算 A|C 非運(yùn)算 A 異或運(yùn)算 xor(A,C)矩陣的邏輯運(yùn)算ppt課件.28 各種允許的比較關(guān)系 , =, , AA = 1 2 3 4 5 6 7 8 0 find(A=5), 大于或等于5元素的下標(biāo) ans = 3 5 6 8矩陣的比較運(yùn)算ppt課件.29 i,j=find(A=5);i,j 顯示行標(biāo)，列標(biāo)ans = A= 3 1 1 2 3 2 2 4 5 6 3 2 7 8 0 2 3 all(A=5) 某列元素全大于或等于5時(shí)，相應(yīng)元素為1，否則為0。ans = 0 0 0 any(A=5) 某列元素中含有大于或等于5時(shí)，相應(yīng)元素為1，否則為0。ans = 1

14、 1 1ppt課件.30解析結(jié)果的化簡(jiǎn)與變換MATLAB 實(shí)現(xiàn)： s1=simple(s) 從各種方法中自動(dòng)選擇最簡(jiǎn)格式 s1,how=simple(s) 化簡(jiǎn)并返回實(shí)際采用的化簡(jiǎn)方法 其中，s為原始表達(dá)式，s1為化簡(jiǎn)后表達(dá)式，how為采用的化簡(jiǎn)方法。 其他常用化簡(jiǎn)函數(shù)（信息與格式可用 help命令得出） collect( ) 合并同類項(xiàng) expand( ) 展開(kāi)多項(xiàng)式 factor( ) 因式分解 numden( ) 提取多項(xiàng)式的分子和分母 sincos( ) 三角函數(shù)的化簡(jiǎn)ppt課件.31例： syms s; P=(s+3)2*(s2+3*s+2)*(s3+12*s2+48*s+64)P

15、=(s+3)2*(s2+3*s+2)*(s3+12*s2+48*s+64) simple(P) % 一系列化簡(jiǎn)嘗試，得出計(jì)算機(jī)認(rèn)為的最簡(jiǎn)形式ans = (s+3)2*(s+2)*(s+1)*(s+4)3 ppt課件.32 a,m=simple(P) % 返回化簡(jiǎn)方法為因式分解方法，用 factor( ) 函數(shù)將得同樣結(jié)果 a = (s+3)2*(s+2)*(s+1)*(s+4)3m =factor expand(P) ans = s7+21*s6+185*s5+883*s4+2454*s3+3944*s2+3360*s+1152ppt課件.33 變量替換 其中，f為原表達(dá)式，用x*替換x得出新

16、的。例：求其 Taylor 冪級(jí)數(shù)展開(kāi) syms a b c d t; % 假設(shè)這些變量均為符號(hào)變量 f=cos(a*t+b)+sin(c*t)*sin(d*t); % 定義給定函數(shù) f(t) f1=subs(f,a,b,c,d,t,0.5*pi,pi,0.25*pi,0.125*pi,4)f1 = -1.0000*111(,)fsubs f xx*11212(,)nnfsubs fx xxxxxppt課件.34基本數(shù)論運(yùn)算下取整、上取整、四舍五入、離0近方向取整、最簡(jiǎn)有理數(shù)、求模的余數(shù)、最大公約數(shù)、最小公倍數(shù)、質(zhì)因數(shù)分解、判定是否為質(zhì)數(shù)ppt課件.35例：對(duì)下面的數(shù)據(jù)進(jìn)行取整運(yùn)算 -0.27

17、65, 0.5772, 1.4597, 2.1091, 1.191, -1.6187 A=-0.2765,0.5772,1.4597,2.1091,1.191,-1.6187; floor(A) % 向 -inf 方向取整ans = -1 0 1 2 1 -2 ceil(A) % 向 +inf 方向取整ans = 0 1 2 3 2 -1 round(A) %取最近的整數(shù)ans = 0 1 1 2 1 -2 fix(A) %向 0 的方向取整ans = 0 0 1 2 1 -1ppt課件.36 例：3x3 Hilbert 矩陣，試用 rat() 函數(shù)變換 A=hilb(3); n,d=rat(

18、A) 將元素變換成最小有理數(shù)，n,d分別為分子、分母矩陣。n = 1 1 1 1 1 1 1 1 1d = 1 2 3 2 3 4 3 4 5ppt課件.37 例:1856120，1483720，最大公約數(shù)、最小公倍數(shù)，質(zhì)因數(shù)分解。 format long m=1856120; n=1483720; gcd(m,n), lcm(m,n)求m,n的最大公約數(shù)、最小公倍數(shù)。ans = 1.0e+009 * 0.00000196000000 1.40508284000000 factor(lcm(n,m) 對(duì)lcm(n,m)進(jìn)行質(zhì)因數(shù)分解。ans = 2 2 2 5 7 7 757 947ppt課件

19、.38 例：1-100間質(zhì)數(shù) A=1:10; isprime(A) 若向量A中某個(gè)整數(shù)值為質(zhì)數(shù)，則相應(yīng)位置為1，其他為零。ans = 0 1 1 0 1 0 1 0 0 0 A=1:100; B=A(isprime(A)B = Columns 1 through 16 2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 Columns 17 through 25 59 61 67 71 73 79 83 89 97 rem(A,C) %A中元素對(duì)C中元素求模得出的余數(shù)。ppt課件.39 2.4.1 循環(huán)結(jié)構(gòu) for 結(jié)構(gòu) while 結(jié)構(gòu)2.4 MATLAB

20、 語(yǔ)言流程控制ppt課件.40例：用循環(huán)求解 s=0;for i=1:100 s=s+i;end s=0; i=1;while (i sum(1:100)ans = 5050例：用循環(huán)求解求最小的 m s=0; m=0;while (s tic, s=0; for i=1:100000, s=s+1/2i+1/3i; end; tocelapsed_time = 1.1820 tic, i=1:100000; s=sum(1./2.i+1./3.i); toc 向量化所需時(shí)間少elapsed_time = 0.3010 i=1:10;s=1./2.i+1./3.i, ss=sum(1./2.i

21、+1./3.i)s = 0.8333 0.3611 0.1620 0.0748 0.0354 0.0170 0.0083 0.0041 0.0020 0.0010ss = 1.4990ppt課件.422.4.2 轉(zhuǎn)移結(jié)構(gòu)ppt課件.43ppt課件.44例：用循環(huán)求解求最大的 m s=0;for i=1:10000 s=s+i; if s10000, break; endend ii = 141ppt課件.452.4.3 開(kāi)關(guān)結(jié)構(gòu)ppt課件.46ppt課件.47 和 C 語(yǔ)言的區(qū)別當(dāng)開(kāi)關(guān)表達(dá)式的值等于某表達(dá)式，執(zhí)行該語(yǔ)句后結(jié)束該結(jié)構(gòu)，不用 break當(dāng)需要在開(kāi)關(guān)表達(dá)式滿足若干個(gè)表達(dá)式之一時(shí)執(zhí)行某

22、一程序段，則用單元形式 （用大括號(hào)把這些表達(dá)式括起來(lái)，用逗號(hào)分隔）otherwise 語(yǔ)句，不是C語(yǔ)言中的 default（但與之等價(jià)）程序的執(zhí)行結(jié)果和各個(gè)case順序無(wú)關(guān)case 語(yǔ)句中條件不能重復(fù)，否則列在后面的條件將不能執(zhí)行ppt課件.48 全新結(jié)構(gòu)（首先試探性執(zhí)行語(yǔ)句1，若執(zhí)行過(guò)程中有錯(cuò)，將錯(cuò)誤信息賦給保留的lasterr變量，并終止這段語(yǔ)句的執(zhí)行，轉(zhuǎn)而執(zhí)行語(yǔ)句2。） 應(yīng)將不保險(xiǎn)但快的算法放在語(yǔ)句1，保險(xiǎn)的放在語(yǔ)句2；或語(yǔ)句2說(shuō)明語(yǔ)句1失效原因。2.4.4 試探結(jié)構(gòu)ppt課件.49 函數(shù)是 MATLAB 編程的主流方法 除了函數(shù)外，還可以采用 M-script（M腳本文件） 文件 M-

23、script 適合于小規(guī)模運(yùn)算例：若最大值不為 10000，需修改程序?qū)?m 和 10000 值的設(shè)置，不適合于M-script2.5 MATLAB 函數(shù)的編寫ppt課件.502.5.1 MATLAB 語(yǔ)言的函數(shù)的基本結(jié)構(gòu) nargin， nargout 分別表示輸入和返回變量的實(shí)際個(gè)數(shù)，此為MATLAB保留變量，只要進(jìn)入該函數(shù)， MATLAB就將自動(dòng)生成這兩個(gè)變量。varargin, varargout 輸入、輸出變量列表（可變輸入輸出個(gè)數(shù)）。ppt課件.51例：前面的要求，m, 10000 function m,s=findsum(k) s=0; m=0; while (s m1,s1=f

24、indsum(145323)m1 = 539s1 = 145530 無(wú)需修改程序ppt課件.52例： 若只給出一個(gè)輸入?yún)?shù)，則會(huì)自動(dòng)生成一個(gè)方陣 在函數(shù)中給出合適的幫助信息 檢測(cè)輸入和返回變量的個(gè)數(shù) edit myhilb,1/(1)i jhijppt課件.53function A=myhilb(n, m)%產(chǎn)生A=MYHILB(N,M)或A=MYHILB(N)；if nargout1, error(Too many output arguments.); endif nargin=1, m=n; elseif nargin=0 | nargin2 error(Wrong number of

25、input arguments.); endA1=zeros(n,m);for i=1: n for j=1:m A1(i,j)=1/(i+j-1);end, endif nargout=1, A=A1; elseif nargout=0, disp(A1); endppt課件.54 help myhilb 產(chǎn)生A=MYHILB(N,M)或A=MYHILB(N)； A=myhilb(3,4)A = 1.0000 0.5000 0.3333 0.2500 0.5000 0.3333 0.2500 0.2000 0.3333 0.2500 0.2000 0.1667 A=myhilb(4)A =

26、1.0000 0.5000 0.3333 0.2500 0.5000 0.3333 0.2500 0.2000 0.3333 0.2500 0.2000 0.1667 0.2500 0.2000 0.1667 0.1429 A=myhilb(3,4,5)? Error using = myhilbToo many input arguments.ppt課件.55 例：函數(shù)的遞歸調(diào)用：階乘function k=my_fact(n)if nargin=1, error(輸入變量個(gè)數(shù)錯(cuò)誤，只能有一個(gè)輸入變量); endif nargout1, error(輸出變量個(gè)數(shù)過(guò)多); endif abs(n

27、-floor(n)eps | n1 % 如果 n1, 進(jìn)行遞歸調(diào)用 k=n*my_fact(n-1);elseif any(0 1=n) % 0!=1!=1 k=1;endppt課件.56 my_fact(11)ans = 39916800 其實(shí)MATLAB提供了求取階乘的函數(shù)factorial(),其核心算法為 prod(1:n),從結(jié)構(gòu)上更簡(jiǎn)單、直觀，速度也更快。 prod(1:11)ans = 39916800 prod(1:3:11)ans = 280ppt課件.57例： conv( ) 可以計(jì)算兩個(gè)多項(xiàng)式的積用 varargin 實(shí)現(xiàn)任意多個(gè)多項(xiàng)式的積function a=convs

28、(varargin) a=1; for i=1:length(varargin), a=conv(a,varargini); end P=1 2 4 0 5; Q=1 2; F=1 2 3; D=convs(P,Q,F)D = 1 6 19 36 45 44 35 30 poly2sym(D)ans =x7+6*x6+19*x5+36*x4+45*x3+44*x2+35*x+302.5.2 可變輸入輸出個(gè)數(shù)ppt課件.58 E=conv(conv(P,Q),F) % 若采用 conv() 函數(shù)，則需要嵌套調(diào)用E = 1 6 19 36 45 44 35 30 poly2sym(E)ans =x

29、7+6*x6+19*x5+36*x4+45*x3+44*x2+35*x+30 G=convs(P,Q,F,1,1,1,3,1,1)G = 1 11 56 176 376 578 678 648 527 315 90ppt課件.592.5.3 inline 函數(shù)和匿名函數(shù) f=inline(sin(x.2+y.2),x,y)MATLAB 7.022( , )sin()f x yxyppt課件.602.6 二維圖形繪制2.6.1 二維圖形繪制基本語(yǔ)句構(gòu)造向量:ppt課件.61ppt課件.62例：選項(xiàng)為紅色點(diǎn)劃線且每個(gè)轉(zhuǎn)折點(diǎn)用五角星表示r-.pentagramppt課件.63例： x=-pi : 0

30、.05: pi; % 以 0.05 為步距構(gòu)造自變量向量 y=sin(tan(x)-tan(sin(x); % 求出各個(gè)點(diǎn)上的函數(shù)值 plot(x,y) plot(x,y,r-.pentagram)ppt課件.64 x=-pi:0.05:-1.8,-1.801:.001:-1.2, -1.2:0.05:1.2,. 1.201:0.001:1.8, 1.81:0.05:pi; % 以變步距方式構(gòu)造自變量向量 y=sin(tan(x)-tan(sin(x); % 求出各個(gè)點(diǎn)上的函數(shù)值 plot(x,y) % 繪制曲線ppt課件.65例： x=-2:0.02:2; % 生成自變量向量 y=1.1*s

31、ign(x).*(abs(x)1.1) + x.*(abs(x) plot(-2,-1.1,1.1,2,-1.1,-1.1,1.1,1.1)ppt課件.66圖形元素屬性獲取與修改 圖形中，每條曲線、坐標(biāo)軸、圖形窗口分別是一個(gè)對(duì)象?？捎胹et( )函數(shù)設(shè)置對(duì)象的屬性，用get( )函數(shù)獲得對(duì)象的某個(gè)屬性。這兩個(gè)語(yǔ)句在界面編程中特別有用。ppt課件.672.6.2 其他二維圖形繪制語(yǔ)句 二維條形圖、羅盤圖、羽毛狀圖、直方圖、極坐標(biāo)圖、階梯圖形、x-半對(duì)數(shù)圖、彗星狀軌跡圖、誤差限圖形、二維填充圖、對(duì)數(shù)圖、磁力線圖、火柴桿圖、y-半對(duì)數(shù)圖。ppt課件.68例：繪制極坐標(biāo)曲線 theta=0:0.01:

32、6*pi; rho=5*sin(4*theta/3); polar(theta,rho) rho=5*sin(theta/3); polar(theta,rho) 周期確定，可以采用試湊方法ppt課件.69例：用不同曲線繪制函數(shù)表示正弦曲線 t=0:.2:2*pi; y=sin(t); % 先生成繪圖用數(shù)據(jù) subplot(2,2,1), stairs(t,y) % 分割窗口，在左上角繪制階梯曲線 subplot(2,2,2), stem(t,y) % 火柴桿曲線繪制 subplot(2,2,3), bar(t,y) % 條型圖繪制 subplot(2,2,4), semilogx(t,y)

33、% 橫坐標(biāo)為對(duì)數(shù)的曲線ppt課件.702.6.3 隱函數(shù)繪制及應(yīng)用 隱函數(shù) 例： ezplot(x2 *sin(x+y2) +y2*exp(x+y)+5*cos(x2+y) x自選 ezplot(x2 *sin(x+y2) +y2*exp(x+y)+5*cos(x2+y),-10 10)ppt課件.712.7 三維圖形繪制 2.7.1 三維曲線繪制 stem3（三維火柴桿型曲線）, fill3（三維填充圖形）, bar3（三維直方圖）等。ppt課件.72例：參數(shù)方程 t=0:.1:2*pi; % 構(gòu)造 t 向量，注意下面的點(diǎn)運(yùn)算 x=t.3.*sin(3*t).*exp(-t); y=t.3.

34、*cos(3*t).*exp(-t); z=t.2; plot3(x,y,z), grid % 三維曲線繪制 stem3(x,y,z); hold on; plot3(x,y,z), gridppt課件.732.7.2 三維曲面繪制 一般曲面繪制 mesh( )繪制網(wǎng)格圖,surf( )繪制表面圖。 其他函數(shù)，光照下 surfl( ),等高線surfc( )，瀑布型waterfall( ) 等高線繪制 contour( ), contour3( )ppt課件.74例:Butterworth 濾波器 x,y=meshgrid(0:31); n=2; D0=200; D=sqrt(x-16).2+(y-16).2); % 求距離 z=1./(1+D.(2*n)/D0); mesh(x,y,z), % 計(jì)算并繪制濾波器 axis(0,31,0,31,0,1) % 重新設(shè)置坐標(biāo)系，增大可讀性 surf(x,y,z) % 繪制三維表面圖 contour3(x,y,z,30) 三維等高線圖，30等高線條數(shù)ppt課件.75ppt課件.76 例：試?yán)L制出二元函數(shù) x,y=meshgrid(-2:.1:2); z=1./(sqrt(1-x).2+y.2)+1./(sqrt(