管理數(shù)學(xué)實(shí)驗(yàn)- MATLAB在管理運(yùn)籌中的應(yīng)用 2-管理數(shù)學(xué)實(shí)驗(yàn)-MATLAB基礎(chǔ)

第2章MATLAB基礎(chǔ)

1、MATLAB語(yǔ)言概述

2、MATLAB的基本運(yùn)算

3、MATLAB繪圖

4、MATLAB的符號(hào)處理

5、MATLAB數(shù)據(jù)的輸入輸出6、MATLAB數(shù)據(jù)統(tǒng)計(jì)描述和分析1、MATLAB語(yǔ)言概述1.1MATLAB語(yǔ)言的發(fā)展matlab語(yǔ)言是由美國(guó)的CleverMoler博士于1980年開(kāi)發(fā)的MATLAB=

MatrixLaboratory它將一個(gè)優(yōu)秀軟件的易用性與可靠性、通用性與專業(yè)性、一般目的的應(yīng)用與高深的科學(xué)技術(shù)應(yīng)用有機(jī)的結(jié)合MATLAB是一種直譯式的高級(jí)語(yǔ)言，比其它程序設(shè)計(jì)語(yǔ)言簡(jiǎn)單易用2000年，發(fā)布MATLAB6.0；2002年，發(fā)布MATLAB6.5；2004年，發(fā)布MATLAB7.0；2005年，發(fā)布MATLAB7.1；2006年3月，發(fā)布MATLABR2006a（V7.2）；2006年9月，發(fā)布MATLABR2006b（V7.3）；2007年3月，發(fā)布MATLABR2007a（V7.4）；2007年9月，發(fā)布MATLABR2007b（V7.5）；2008年3月，發(fā)布MATLABR2008a（V7.6）；2008年9月，發(fā)布MATLABR2086b（V7.7）；2009年3月，發(fā)布MATLABR2009a（V7.8）；2009年9月，發(fā)布MATLABR2009b（V7.9）；2010年3月，發(fā)布MATLABR2010a（V7.10）；2010年9月，發(fā)布MATLABR2010b（V7.11）；2011年4月，發(fā)布MATLABR2011a（V7.12）；2011年9月，發(fā)布MATLABR2011b（V7.13）；2012年3月，發(fā)布MATLABR2012a（V7.14）；2012年9月，發(fā)布MATLABR2012b（V8.0）；2013年3月，發(fā)布MATLABR2013a（V8.1）。MATLABR2017b。1.2matlab在各領(lǐng)域的應(yīng)用工業(yè)研究與開(kāi)發(fā)數(shù)學(xué)教學(xué)，特別是線性代數(shù)數(shù)值分析和科學(xué)計(jì)算方面的教學(xué)與研究電子學(xué)、控制理論和物理學(xué)等工程和科學(xué)經(jīng)濟(jì)學(xué)、化學(xué)和生物學(xué)等計(jì)算問(wèn)題的所有其他領(lǐng)域中的教學(xué)與研究各種模型的分析與計(jì)算MATLAB語(yǔ)言的簡(jiǎn)潔高效性MATLAB強(qiáng)大的科學(xué)運(yùn)算功能MATLAB豐富的繪圖功能MATLAB龐大的工具箱與模塊集MATLAB強(qiáng)大的動(dòng)態(tài)系統(tǒng)仿真功能MATLAB作為線性系統(tǒng)的一種分析和仿真工具，是理工科學(xué)生應(yīng)該掌握的技術(shù)工具，它作為一種編程語(yǔ)言和可視化工具，可解決工程、科學(xué)計(jì)算和數(shù)學(xué)學(xué)科中許多問(wèn)題。1.3MATLAB語(yǔ)言的功能強(qiáng)大的數(shù)值（矩陣）運(yùn)算廣泛的符號(hào)運(yùn)算功能高級(jí)與低級(jí)兼?zhèn)涞膱D形功能計(jì)算結(jié)果的可視化可靠的容錯(cuò)功能應(yīng)用靈活的兼容與接口功能信息量豐富的聯(lián)機(jī)檢索功能(1)矩陣運(yùn)算功能MATLAB提供了豐富的矩陣運(yùn)算處理功能基于矩陣運(yùn)算的處理工具。變量

矩陣，運(yùn)算

矩陣的運(yùn)算例如C=A+B，A,B,C都是矩陣,是矩陣的加運(yùn)算一個(gè)數(shù)，x=5，MATLAB看做是一個(gè)11的矩陣(2)符號(hào)運(yùn)算功能符號(hào)運(yùn)算即用字符串進(jìn)行數(shù)學(xué)分析允許變量不賦值而參與運(yùn)算用于解代數(shù)方程、微積分、復(fù)合導(dǎo)數(shù)、積分、二重積分、有理函數(shù)、微分方程、泰樂(lè)級(jí)數(shù)展開(kāi)、尋優(yōu)等等，可求得解析符號(hào)解(3)豐富的繪圖功能與計(jì)算結(jié)果的可視化具有高層繪圖功能——兩維、三維繪圖具有底層繪圖功能——句柄繪圖使用plot函數(shù)可隨時(shí)將計(jì)算結(jié)果可視化(4)圖形化程序編制功能動(dòng)態(tài)系統(tǒng)進(jìn)行建模、仿真和分析的軟件包用結(jié)構(gòu)圖編程，而不用程序編程只需拖幾個(gè)方塊、連幾條線，即可實(shí)現(xiàn)編程功能(5)豐富的MATLAB工具箱MATLAB主工具箱符號(hào)數(shù)學(xué)工具箱SIMULINK仿真工具箱控制系統(tǒng)工具箱信號(hào)處理工具箱圖象處理工具箱通訊工具箱系統(tǒng)辨識(shí)工具箱神經(jīng)元網(wǎng)絡(luò)工具箱金融工具箱優(yōu)化工具箱……MATLAB中變量的命名規(guī)則是：（1）變量名必須是不含空格的字符串；（2）變量名以字母開(kāi)頭，之后可以是任意字母、數(shù)字或下劃線，變量名中不允許使用標(biāo)點(diǎn)符號(hào).（3）變量名區(qū)分大小寫；（4）變量名最多不超過(guò)19個(gè)字符；

MATLAB的保留常量eps,pi,NaN,Inf,i,j,i=sqrt(-1)lastwarn,lasterr

1.4MATLAB語(yǔ)言基礎(chǔ)特殊變量表數(shù)值型數(shù)據(jù)結(jié)構(gòu)(1)雙精度數(shù)值變量IEEE標(biāo)準(zhǔn)，64位(8字節(jié))，11指數(shù)位，53數(shù)值位和一個(gè)符號(hào)位-1.7*10308

至1.7*10308double()函數(shù)的轉(zhuǎn)換(2)其他數(shù)據(jù)類型uint8，常用于圖像表示和處理，8位int8(),int16(),int32(),uint16(),uint32()符號(hào)型變量數(shù)據(jù)類型(3)符號(hào)型，sym(A),常用于公式推導(dǎo)變量聲明symsABCpositivesym(sin(x))采用變精度函數(shù)求值vpa(A),或vpa(A,n)例求出p

的300位有效數(shù)字

vpa(pi,300)MATLAB支持的其他數(shù)據(jù)結(jié)構(gòu)(4)基本數(shù)值變量類型：雙精度復(fù)數(shù)矩陣(5)字符串型數(shù)據(jù)：用單引號(hào)括起來(lái)(6)多維數(shù)組：是矩陣的直接擴(kuò)展，多個(gè)下標(biāo)(7)單元數(shù)組：將不同類型數(shù)據(jù)集成到一個(gè)變量名下面，用{}表示。例。(8)結(jié)構(gòu)體：A.b，引用也用A.b，不是A->b(9)類與對(duì)象：可以定義重載函數(shù)返回例：?jiǎn)卧?元胞)數(shù)組的創(chuàng)建C_str=char('北京理工大學(xué)','管理科學(xué)與工程系');R=reshape(1:9,3,3);Cn=[1+2i,1-2i];

S_sym=sym('sin(-3*t)*exp(-t)');%符號(hào)函數(shù)量（1）直接創(chuàng)建法之一：“外標(biāo)識(shí)單元元素賦值法”（2）直接創(chuàng)建法之二：“編址單元元素內(nèi)涵的直接賦值法”A(1,1)={C_str};A(1,2)={R};A(2,1)={Cn};A(2,2)={S_sym};A%顯示單元數(shù)組類型B{1,1}=C_str;B{1,2}=R;B{2,1}=Cn;B{2,2}=S_sym;celldisp(B)%顯示單元數(shù)組內(nèi)容結(jié)構(gòu)體示例：通過(guò)溫室數(shù)據(jù)（包括溫室名、容積、溫度、濕度等）演示：結(jié)構(gòu)體green_house的創(chuàng)建和顯示。(它有三個(gè)域：name，volume，parameter。而parameter又有三個(gè)子域：temperature，humidity。)（1）直接對(duì)域賦值產(chǎn)生“結(jié)構(gòu)體”，即結(jié)構(gòu)數(shù)組。green_=‘一號(hào)房’；結(jié)構(gòu)的域由（結(jié)構(gòu)名）.（域名）標(biāo)識(shí)。green_house.volume='2000立方米';green_house.parameter.temperature=…[31.230.431.628.7

29.731.130.929.6];green_house.parameter.humidity…=[62.159.557.761.562.061.959.257.5];（2）顯示“結(jié)構(gòu)體”結(jié)構(gòu)和內(nèi)容：green_house，……返回MATLAB的基本語(yǔ)句結(jié)構(gòu)

>>A=[1,2,3;45,6;7,8,0]>>A=[1,2,3;45,6;7,8,0];>>A=[[A;[1,2,3]][4;5;6;7];函數(shù)調(diào)用語(yǔ)句冒號(hào)表達(dá)式[返回變量列表]=函數(shù)名（輸入變量列表）[a,b,c]=myfun(d,e,f,g)v=s1:s2:s3例用不同的步距生成(0,p)間向量

>>v1=0:0.2:pi

>>v2=0:-0.1:pi

>>v3=0:pi

>>v4=pi:-1:0

>>v5=[0:0.2:pi,pi]

>>v6=[pi:-1:00]子矩陣提取基本語(yǔ)句格式例子矩陣提取提取A

矩陣全部奇數(shù)行，所有列提取A

矩陣3,2,1行、2,3,4列構(gòu)成子矩陣將A

矩陣左右翻轉(zhuǎn)B3=A(:,end:-1:1)B1=A(1:2:end,:)B2=A([3,2,1],[2,3,4])B=A(v1,v2)2.1創(chuàng)建矩陣2.2矩陣運(yùn)算2.3多項(xiàng)式運(yùn)算2.4線性方程組2.5MATLAB語(yǔ)言流程控制2.6數(shù)學(xué)函數(shù)2.7.m文件2、MATLAB的基本運(yùn)算2.1創(chuàng)建矩陣規(guī)則：矩陣元素必須用[]括住矩陣元素必須用逗號(hào)或空格分隔

[]內(nèi)矩陣的行與行之間必須用分號(hào)分隔矩陣元素任何matlab表達(dá)式可以是實(shí)數(shù)，也可以是復(fù)數(shù)，復(fù)數(shù)可用i,j輸入

a=[123;456]x=[2pi/2;sqrt(3)3+5i]

矩陣的建立逗號(hào)或空格作用：分隔某一行的元素。分號(hào)作用：區(qū)分不同的行。Enter鍵作用之一：輸入矩陣時(shí)，按Enter鍵開(kāi)始一新行。注意：輸入矩陣時(shí)，嚴(yán)格要求所有行有相同的列。

例m=[1234；5678；9101112]p=[111122223333]創(chuàng)建簡(jiǎn)單的數(shù)組x=[a,b,

c,

d,

e,

f]創(chuàng)建包含指定元素的行向量x=first：last

創(chuàng)建從first開(kāi)始，加1計(jì)數(shù)，到last結(jié)束的行向量x=first：increment：last創(chuàng)建從first開(kāi)始，加increment計(jì)數(shù)，last結(jié)束的行向量x=linspace(first，last，n）創(chuàng)建從first開(kāi)始，到last結(jié)束，有n個(gè)元素的行向量x=logspace(first，last，n）創(chuàng)建從10first開(kāi)始，到10last結(jié)束，n個(gè)元素為10.^linspace(first，last，n）的對(duì)數(shù)分隔行向量數(shù)組元素的訪問(wèn)（1）訪問(wèn)一個(gè)元素：x(i)表示訪問(wèn)數(shù)組x的第i個(gè)元素。（2）訪問(wèn)一塊元素：x(b：s：e)表示訪問(wèn)數(shù)組x的從第b個(gè)元素開(kāi)始，以步長(zhǎng)為s到第e個(gè)元素（但不超過(guò)e），s可以為負(fù)數(shù)，s缺省時(shí)為1。（3）直接使用元素編址序號(hào)：x([ijkl])表示提取數(shù)組x的第i、i、k、l個(gè)元素構(gòu)成一個(gè)新的數(shù)組[x(i)，x(j)，x(k)，x(l)]。符號(hào)的作用逗號(hào)和分號(hào)的作用逗號(hào)和分號(hào)可作為指令間的分隔符，matlab允許多條語(yǔ)句在同一行出現(xiàn)。分號(hào)如果出現(xiàn)在指令后，屏幕上將不顯示結(jié)果。當(dāng)一個(gè)指令或矩陣太長(zhǎng)時(shí)，可用”…”續(xù)行冒號(hào)的作用用于生成等間隔的向量，默認(rèn)間隔為1。用于選出矩陣指定行、列及元素。循環(huán)語(yǔ)句用matlab函數(shù)創(chuàng)建矩陣a=[]——空陣[]，大小為零，Matlab允許輸入空陣，當(dāng)一項(xiàng)操作無(wú)結(jié)果時(shí)，返回空陣。rand——隨機(jī)矩陣，b=rand([m,]n)eye——單位矩陣，c=eye(n)zeros——全部元素都為0的矩陣，d=zeros([m,]n)ones——全部元素都為1的矩陣，e=ones(m，n)還有伴隨矩陣、稀疏矩陣、魔方矩陣、對(duì)角矩陣、范德蒙等矩陣的創(chuàng)建，就不一一介紹了。注意：matlab嚴(yán)格區(qū)分大小寫字母。

matlab函數(shù)名必須小寫。矩陣中元素的操作（1）矩陣A的第r行：A（r，：）（2）矩陣A的第j列：A（：，j）（3）依次提取矩陣A的每一列，將A拉伸為一列向量：

A（：）（4）取矩陣A的第i1~i2行、第j1~j2列構(gòu)成新矩陣：

A（i1：i2，j1：j2）（5）逆序提取矩陣A的第i1~i2行，構(gòu)成新矩陣：

A（i2：-1：i1，：）（6）逆序提取矩陣A的第j1~j2列，構(gòu)成新矩陣：

A（：，j2：-1：j1）（7）刪除A的第i1~i2行，構(gòu)成新矩陣:A（i1：i2，：）=[]（8）刪除A的第j1~j2列，構(gòu)成新矩陣:A（：，j1：j2）=[]（9）將矩陣A和B拼接成新矩陣：[A，B]；[A；B]2.2矩陣運(yùn)算1）矩陣加、減（＋,－）運(yùn)算規(guī)則：相加、減的兩矩陣必須有相同的行和列；允許參與運(yùn)算的兩矩陣之一是標(biāo)量：標(biāo)量與矩陣的所有元素分別進(jìn)行加減操作。2）矩陣乘（）運(yùn)算規(guī)則：A矩陣的列數(shù)必須等于B矩陣的行數(shù)標(biāo)量可與任何矩陣相乘。a=[123;456;780];b=[1;2;3];c=a*bc=143223

3）矩陣乘方——a^n,a^p,p^a

a^p——a自乘p次冪

方陣>1的整數(shù)對(duì)于p的其它值，計(jì)算將涉及特征值和特征向量如果p是矩陣，a是標(biāo)量a^p使用特征值和特征向量自乘到p次冪；如a,p都是矩陣，a^p則無(wú)意義。當(dāng)一個(gè)方陣有復(fù)數(shù)特征值或負(fù)實(shí)特征值時(shí)，非整數(shù)冪是復(fù)數(shù)陣。if[V,D]=eig(a),thena^p=V*D.^p/V4）矩陣的其它運(yùn)算inv——矩陣求逆，inv（A）det——方陣行列式的值，det（A）eig——矩陣的特征值，[V，D]=eig[A]diag——對(duì)角矩陣，

'——矩陣轉(zhuǎn)置sqrt——矩陣開(kāi)方5）借助cat,repmat,reshape等函數(shù)構(gòu)作高維數(shù)組

cat(3,ones(2,3),ones(2,3)*2,ones(2,3)*3)

ans(:,:,1)=ans(:,:,2)=ans(:,:,3)=111222333111222333reshape(1:12,2,2,3)

ans(:,:,1)=ans(:,:,2)=ans(:,:,3)=135791124681012repmat(ones(2,3)*2,[1,1,3])ans(:,:,1)=ans(:,:,2)=ans(:,:,3)=2222222222222222226）矩陣的一些特殊操作矩陣的變維

a=[1:12];b=reshape(a,3,4)c=zeros(3,4);c(:)=a(:)矩陣的變向

rot90:旋轉(zhuǎn);fliplr:左右翻;flipud:上下翻矩陣的抽取

diag:抽取主對(duì)角線;tril:抽取主下三角;

triu:抽取主上三角維數(shù)、大小、長(zhǎng)度和降維

ndims(A);size(A);length(A);squeeze(A)

7）矩陣的數(shù)組運(yùn)算

數(shù)組運(yùn)算指元素對(duì)元素的算術(shù)運(yùn)算，與通常意義上的由符號(hào)表示的線性代數(shù)矩陣運(yùn)算不同數(shù)組加減(.+,.-)（與矩陣加減等效）

a.+b=a+ba.-b=a–b數(shù)組乘除(，./，.\）ab——a，b兩數(shù)組必須有相同的行和列兩數(shù)組對(duì)應(yīng)元素相乘。a./b=b.\aa.\b=b./a——

給出a,b對(duì)應(yīng)元素間的商例:a=[123];b=[456];c1=a.\b;c2=b./ac1=4.00002.50002.0000c2=4.00002.50002.0000數(shù)組乘方(.^)—元素對(duì)元素的冪例:a=[123];b=[456];z=a.^2z=1.004.009.00z=a.^bz=1.0032.00729.008）矩陣的邏輯運(yùn)算邏輯變量：當(dāng)前版本有邏輯變量對(duì)double變量來(lái)說(shuō)，非0表示邏輯1邏輯運(yùn)算（相應(yīng)元素間的運(yùn)算）與運(yùn)算A&B或運(yùn)算A|B非運(yùn)算~A異或運(yùn)算xor(A,B)9）矩陣的比較運(yùn)算各種允許的比較關(guān)系

>,>=,<,<=,==,~=,find(),all(),any()實(shí)例Linear_index=find(R>0)

>>A=[1,2,3;4,5,6;7,8,9];find(A>=5)’,any(A>=5)>>[i,j]=find(A>=5)>>all(A>=5),A>=510）數(shù)據(jù)類型轉(zhuǎn)換函數(shù)int2str；num2str；mat2str；str2mat；str2numeval(Expression)

eval(expression,catch_expr)[a1,a2,a3,...]=eval(‘function(b1,b2,b3,...)’)strcat(s1,s2,...),strvcat(s1,s2,...)fprintf(fileID,formatSpec,A1,…,An),str=sprintf(formatSpec,A1,…,An),A=sscanf(str,format)A=eye(2,4);A_str1=int2str(A)

rand('state',0);B=rand(2,4);B3=num2str(B,3)B_str=mat2str(B,4)fprintf,sprintf,sscanfs_s=sprintf('%.10e\n',B)

[s,errmsg]=sprintf(format,A,...)11）基本數(shù)論運(yùn)算RemainderGreatestcommondivisorLeastcommonmultiplePrimefactorsRationalfractionapproximation函數(shù)名調(diào)用格式floor()ceil()round()fix()rat()rem()gcd()lcm()factor()isprime()n=floor(x)n=ceil(x)n=round(x)n=fix(x)[n,d]=rat(x)B=rem(A,C)k=gcd(n,m)k=lcm(n,m)k=factor(n)v1=isprime(v)例對(duì)下面的數(shù)據(jù)進(jìn)行取整運(yùn)算

-0.2765,0.5772,1.4597,2.1091,1.191,-1.6187例3×3Hilbert矩陣，試用rat()函數(shù)變換例1856120，1483720，最大公約數(shù)、最小公倍數(shù)，質(zhì)因數(shù)分解例1-1000間質(zhì)數(shù)N./DapproximatingAmatlab語(yǔ)言把多項(xiàng)式表達(dá)成一個(gè)行向量，該向量中的元素是按多項(xiàng)式降冪排列的。

f(x)=anxn+an-1xn-1+……+a0

可用行向量p=[anan-1……a1a0]表示poly——

產(chǎn)生特征多項(xiàng)式系數(shù)向量特征多項(xiàng)式一定是n+1維的特征多項(xiàng)式第一個(gè)元素一定是12.3多項(xiàng)式運(yùn)算例:a=[123;456;780];p=poly(a)p=1.00-6.00-72.00-27.00p是多項(xiàng)式p(x)=x3-6x2-72x-27的matlab描述方法，可用：p1=poly2str(p,‘x’)—函數(shù)文件，顯示數(shù)學(xué)多項(xiàng)式的形式p1=x^3-6x^2-72x-27roots——

求多項(xiàng)式的根a=[123;456;780];p=poly(a)p=1.00-6.00-72.00-27.00r=roots(p)r=12.12-5.73——

顯然r是矩陣a的特征值

-0.39可用poly返回多項(xiàng)式形式p2=poly(r)p2=1.00-6.00-72.00-27.00conv,convs多項(xiàng)式乘運(yùn)算例:a(x)=x2+2x+3;b(x)=4x2+5x+6;c=(x2+2x+3)(4x2+5x+6)a=[123];b=[456];c=conv(a,b)=conv([123],[456])c=4.0013.0028.0027.0018.00p=poly2str(c,'x')p=4x^4+13x^3+28x^2+27x+18deconv多項(xiàng)式除運(yùn)算a=[123];c=[4.0013.0028.0027.0018.00]d=deconv(c,a)d=4.005.006.00[d,r]=deconv(c,a)余數(shù)c除a后的整數(shù)多項(xiàng)式微分matlab提供了polyder函數(shù)多項(xiàng)式的微分。命令格式：polyder(p):求p的微分polyder(a,b):求多項(xiàng)式a,b乘積的微分[p,q]=polyder(a,b):求多項(xiàng)式a,b商的微分

deconv(p,q)例：a=[12345];poly2str(a,'x')ans=x^4+2x^3+3x^2+4x+5b=polyder(a)

b=4664poly2str(b,'x')

ans=4x^3+6x^2+6x+42.4代數(shù)方程組求解Matlab中有兩種除運(yùn)算左除和右除。對(duì)于方程Ax=b，A為m×n矩陣Am×n，有三種情況：當(dāng)n=m時(shí)，此方程成為“恰定”方程當(dāng)n<m時(shí)，此方程成為“超定”方程當(dāng)n>m時(shí)，此方程成為“欠定”方程matlab定義的除運(yùn)算可解上述三種方程1）恰定方程組的解方程Ax=b(A為非奇異)

x=A-1

b兩種解語(yǔ)句:x=inv(A)b

——

采用求逆運(yùn)算解方程x=A\b——

采用左除運(yùn)算解方程2）超定方程組的解方程Ax=b,n<m時(shí)此時(shí)不存在唯一解。方程解(A'A)x=A'bx=(A'A)-1

A'b——

求逆法

x=A\b——matlab用最小二乘法找一個(gè)近似的解。A=[12;23;34];b=[1;2;3];

解1

x=A\b

解2x=inv(A'A)

A'b

x=x=1.001.0000.00

=Ax=b3）欠定方程組的解當(dāng)方程數(shù)少于未知量個(gè)數(shù)時(shí),即不定情況,有無(wú)窮多個(gè)解存在。matlab可求出兩個(gè)解：用除法求的解x是具有最多零元素的解基于偽逆pinv求得具有最小長(zhǎng)度或范數(shù)的解A=[123;234];b=[1;2];

x=A\bx=pinv(A)bx=x=1.000.8300.330-0.17=Ax=b2.4MATLAB語(yǔ)言流程控制

2.4.1循環(huán)結(jié)構(gòu)for結(jié)構(gòu)while結(jié)構(gòu)fori=V

循環(huán)結(jié)構(gòu)體endwhile（條件式）

循環(huán)結(jié)構(gòu)體end例用循環(huán)求解例用循環(huán)求解求最小的m例求>>tic,s=0;fori=1:100000,s=s+1/2^i+1/3^i;end;toc>>tic,i=1:100000;s=sum(1./2.^i+1./3.^i);toc

2.5.2條件結(jié)構(gòu)if(條件表達(dá)式1)%若條件1滿足，則執(zhí)行語(yǔ)句組1

語(yǔ)句組1%可以嵌套下一層的

if結(jié)構(gòu)elseif(條件表達(dá)式2)%若條件2滿足，則執(zhí)行語(yǔ)句組2

語(yǔ)句組2

┆┆%可設(shè)置n種轉(zhuǎn)移條件

else%上面n種條件均不滿足，則執(zhí)行語(yǔ)句組n+1

語(yǔ)句組n+1end

s=0;form=1:100000s=s+i;

ifs>10000

break;end;end;2.5.3開(kāi)關(guān)結(jié)構(gòu)switch(開(kāi)關(guān)表達(dá)式

)case表達(dá)式1%若表達(dá)式1與開(kāi)關(guān)表達(dá)式一致，則執(zhí)行語(yǔ)句組1

語(yǔ)句組1%可以嵌套下一層的

case結(jié)構(gòu)case{表達(dá)式2.1,表達(dá)式2.2,…

表達(dá)式2-m}

語(yǔ)句組2%若條件2之一滿足表達(dá)式

，則執(zhí)行語(yǔ)句組2

┆┆%可設(shè)置n種開(kāi)關(guān)表達(dá)式otherwise%上面n種表達(dá)式均不滿足，則執(zhí)行語(yǔ)句組n+1

語(yǔ)句組n+1end和C語(yǔ)言的區(qū)別當(dāng)開(kāi)關(guān)表達(dá)式的值等于某表達(dá)式，執(zhí)行該語(yǔ)句后結(jié)束該結(jié)構(gòu)，不用break同時(shí)滿足若干個(gè)條件之一，則用單元形式otherwise語(yǔ)句，不是default程序的執(zhí)行結(jié)果和各個(gè)case順序無(wú)關(guān)case語(yǔ)句中條件不能重復(fù)，否則列在后面的條件將不能執(zhí)行2.5.4試探結(jié)構(gòu)全新結(jié)構(gòu)try,

語(yǔ)句段1，

catch,

語(yǔ)句段2，end2.6數(shù)學(xué)函數(shù)函數(shù)是MATLAB編程的主流方法除了函數(shù)外，還可以采用M-script文件M-script適合于小規(guī)模例前例中，若最大值不為10000，需修改程序，如M-script更適合變化的K。

MATLAB的內(nèi)部函數(shù)是有限的，有時(shí)為了研究某一個(gè)函數(shù)的各種性態(tài)，需要為MATLAB定義新函數(shù)，為此必須編寫函數(shù)文件.函數(shù)文件是文件名后綴為M的文件，這類文件的第一行必須是一特殊字符function開(kāi)始，格式為：

function因變量名=函數(shù)名（自變量名）函數(shù)值的獲得必須通過(guò)具體的運(yùn)算實(shí)現(xiàn)，并賦給因變量.

2.7.M文件M文件建立方法：(1)在Matlab中，點(diǎn):File->New->M-file(2)在編輯窗口中輸入程序內(nèi)容

(3)點(diǎn)：File->Save，存盤，M文件名必須與函數(shù)名一致。Matlab的應(yīng)用程序也以M文件保存。例：定義函數(shù)f(x1,x2)=100(x2-x12)2+(1-x1)2function

f=fun(x)f=100*(x(2)-x(1)^2)^2+(1-x(1))^21)建立M文件：fun.m2)可以直接使用函數(shù)fun.m例如：計(jì)算f(1,2),只需在Matlab命令窗口鍵入命令：x=[12]fun(x)2.7.1MATLAB語(yǔ)言函數(shù)的基本結(jié)構(gòu)

nargin，nargout,varargin,varargoutfunction[返回變量列表]=函數(shù)名(輸入變量列表)

注釋說(shuō)明段，有%引導(dǎo)

輸入、返回變量格式的檢測(cè)

函數(shù)體語(yǔ)句段例對(duì)于變化K的解決方案無(wú)需修改程序function[m,s]=findsum(k)s=0;m=0;

while(s<=k),m=m+1;m=m+1;s=s+m;end>>[m,s]=findsum(145323)例若只給出一個(gè)輸入?yún)?shù)，則會(huì)自動(dòng)生成一個(gè)方陣在函數(shù)中給出合適的幫助信息檢測(cè)輸入和返回變量的個(gè)數(shù)editmyhilbn×m階的Hilbert矩陣

functionA=myhilbert(n,m)

ifnargin==1,m=n;elseifnargin==0|nargin>2error(‘Wrongnumberofinputarguments.’);else

fori=1:nforj=1:mA(i,j)=1/(i+j-1);endendendhelpmyhilbertA=myhilbert(3,4)A=myhilbert(4)例函數(shù)的遞歸調(diào)用：階乘functionk=my_fact(n)

ifn>1%如果n>1,則進(jìn)行遞歸調(diào)用

k=n*my_fact(n-1);

elseifany([01]==

n%0!=1!=1

k=1;

endendn!=n*(n-1)!2.7.2可變輸入輸出個(gè)數(shù)例可以計(jì)算兩個(gè)多項(xiàng)式的積用varargin實(shí)現(xiàn)任意多個(gè)多項(xiàng)式的積p=conv(p1,p2)functionp=convs(varargin)p=1;for1:length(varargin),p=conv(p,varargin{i});

endendP=[12405];Q=[12];F=[123]D=convs(P,Q,F)2.7.3inline函數(shù)和匿名函數(shù)inline函數(shù)，可以免去文件MATLAB7.0之后fun=inline(‘函數(shù)內(nèi)容’,自變量列表)fun=@(自變量列表)函數(shù)內(nèi)容,例如，f=@(x,y)sin(x.^2+y.^2)3、Matlab繪圖——Matlab語(yǔ)言豐富的圖形表現(xiàn)方法，使得數(shù)學(xué)計(jì)算結(jié)果可以方便地、多樣性地實(shí)現(xiàn)了可視化，這是其它語(yǔ)言所不能比擬的。matlab語(yǔ)言的繪圖功能

不僅能繪制幾乎所有的標(biāo)準(zhǔn)圖形，而且其表現(xiàn)形式也是豐富多樣的。

matlab語(yǔ)言不僅具有高層繪圖能力，而且還具有底層繪圖能力——句柄繪圖方法。在面向?qū)ο蟮膱D形設(shè)計(jì)基礎(chǔ)上，使得用戶可以用來(lái)開(kāi)發(fā)各專業(yè)的專用圖形。3.1二維繪圖3.1.1plot——

最基本的二維圖形指令plot的功能：

plot命令自動(dòng)打開(kāi)一個(gè)圖形窗口Figure

用直線連接相鄰兩數(shù)據(jù)點(diǎn)來(lái)繪制圖形根據(jù)圖形坐標(biāo)大小自動(dòng)縮擴(kuò)坐標(biāo)軸，將數(shù)據(jù)標(biāo)尺及單位標(biāo)注自動(dòng)加到兩個(gè)坐標(biāo)軸上，可自定坐標(biāo)軸，可把x,y軸用對(duì)數(shù)坐標(biāo)表示如果已經(jīng)存在一個(gè)圖形窗口，plot命令則清除當(dāng)前圖形，繪制新圖形可單窗口單曲線繪圖；可單窗口多曲線繪圖；可單窗口多曲線分圖繪圖；可多窗口繪圖可任意設(shè)定曲線顏色和線型可給圖形加坐標(biāo)網(wǎng)線和圖形加注功能plot的調(diào)用格式

plot(x)——

缺省自變量繪圖格式，x為向量,以x元素值為縱坐標(biāo)，以相應(yīng)元素下標(biāo)為橫坐標(biāo)繪圖

plot(x,y)——

基本格式，以y(x)的函數(shù)關(guān)系作出直角坐標(biāo)圖，如果y為n×m的矩陣，則以x為自變量，作出m條曲線二維圖形繪制基本語(yǔ)句構(gòu)造向量:x=x1,x2,…,xny=y(x1),y(x2),…,y(xn)x=[x1,x2,…,xn]y=[y(x1),y(x2),…,y(xn)]plot(x,y)xx(x1,y1)，(x2,y2)…,(xn,yn)plot(x1,y1,x2,y2,…,xn,yn)多條曲線繪圖格式x選項(xiàng)取值指定曲線的線型和顏色表plot(x1,y1,選項(xiàng)1,x2,y2,選項(xiàng)2,…,xn,yn,選項(xiàng)n)plot(x,y,’s’)——

開(kāi)關(guān)格式，開(kāi)關(guān)量字符串s設(shè)定曲線顏色和繪圖方式，使用顏色字符串的前1~3個(gè)字母，如

yellow—yel表示等?；騪lot(x1,y1,’s1’,x2,y2,’s2’，…)S的標(biāo)準(zhǔn)設(shè)定值如下：

字母顏色標(biāo)點(diǎn)線型

y黃色·點(diǎn)線

m粉紅○圈線

c亮藍(lán)××線

r大紅＋＋字線

g綠色－實(shí)線

b藍(lán)色星形線

w白色：虛線

k黑色－·(--)點(diǎn)劃線Matlab線形：[+|o|*|.|x|square|diamond|v|^|>|<|pentagram|hexagram]square正方形

diamond菱形

pentagram五角星

hexagram六角星(1)單窗口單曲線繪圖例：x=[0,0.48,0.84,1,0.91,0.6,0.14][x1,x2,x3,x4,x5,x6,x7,]plot(x)(2)單窗口多曲線繪圖例：x=0:pi/100:2*pi;y=sin(x);y1=sin(x+0.25);y2=sin(x+0.5);plot(x,y,x,y1,x,y2)yy1y2例：y=sin(x);y1=sin(x+0.25);y2=sin(x+0.5);

y3=cos(x);y4=cos(x+0.25);y5=cos(x+0.5);

plot(x,[y',y1',y2',y3',y4',y5'])y3=cos(x);y4=cos(x+0.25);y5=cos(x+0.5);plot(x,y3);holdon;plot(x,y4);plot(x,y5);x=peaks;plot(x)x=1:length(peaks);y=peaks;plot(x,y)(3)單窗口多曲線分圖繪圖subplot——

子圖分割命令調(diào)用格式：

subplot(m,n,p)

——

按從左至右，從上至下排列行列繪圖序號(hào)subplot(1,3,1);plot(x,y)

subplot(1,3,2);plot(x,y3)

subplot(1,3,3);plot(x,y2)subplot(3,1,1);

plot(x,y)

subplot(3,1,2);

plot(x,y3)

subplot(3,1,3);

plot(x,y2)

subplot('position',[leftbottomwidthheight])Left,bottom,width,height

∈[0,1]clf;x=(pi*(0:1000)/1000)';y1=sin(x);y2=sin(10*x);y12=sin(x).*sin(10*x);subplot(2,2,1),plot(x,y1);axis([0,pi,-1,1])subplot(2,2,2),plot(x,y2);axis([0,pi,-1,1])subplot('position',[0.2,0.05,0.6,0.45]) plot(x,y12,'b-',x,[y1,-y1],'r:');axis([0,pi,-1,1])

subplot('position',[leftbottomwidthheight])Left,bottom,width,height

∈[0,1]clf;x=(pi*(0:1000)/1000)';y1=sin(x);y2=sin(10*x);y12=sin(x).*sin(10*x);subplot(2,2,1),plot(x,y1);axis([0,pi,-1,1])subplot(2,2,2),plot(x,y2);axis([0,pi,-1,1])subplot('position',[0.2,0.05,0.6,0.45]) plot(x,y12,'b-',x,[y1,-y1],'r:');axis([0,pi,-1,1])

(4)多窗口繪圖figure(n)——?jiǎng)?chuàng)建窗口函數(shù)，n為窗口順序號(hào)。x=0:pi/100:2*pi;y=sin(x);y1=sin(x+0.25);y2=sin(x+0.5);plot(x,y)——自動(dòng)出現(xiàn)第一個(gè)窗口figure(2)plot(x,y1)——

在第二窗口繪圖figure(3)plot(x,y2)——在第三窗口繪圖(5)可任意設(shè)置顏色與線型例：plot(x,y,'r-',x,y1,'g:',x,y2,'b*')(6)圖形加注功能將標(biāo)題、坐標(biāo)軸標(biāo)記、網(wǎng)格線及文字注釋加注到圖形上，這些函數(shù)為：title('string','PropertyName',PropertyValue,...)

——

給圖形加標(biāo)題

xlabel('string')

——

給x軸加標(biāo)注

ylabel('string')

——

給y軸加標(biāo)注

text(x,y,'string')

——

在圖形指定位置加標(biāo)注

gtext('string')

——

將標(biāo)注加到圖形任意位置

gridon(off)

——

打開(kāi)、關(guān)閉坐標(biāo)網(wǎng)格線

legend

——

添加圖例

axis

——

控制坐標(biāo)軸的刻度例：t=0:0.1:10y1=sin(t);y2=cos(t);plot(t,y1,'r',t,y2,'b--');x=[1.7*pi;1.6*pi];y=[-0.3;0.8];s=['sin(t)';'cos(t)'];text(x,y,s);title('正弦和余弦曲線');legend('正弦','余弦')xlabel('時(shí)間t'),ylabel('正弦、余弦')gridaxissquare

axis的用法還有：

axis([xminxmaxyminymax])——用行向量中給出的值設(shè)定坐標(biāo)軸的最大和最小值如axis([-2

205])axis(equal)——

將兩坐標(biāo)軸設(shè)為相等

axison(off)——

顯示和關(guān)閉坐標(biāo)軸的標(biāo)記、標(biāo)志axisauto——

將坐標(biāo)軸設(shè)置返回自動(dòng)缺省值(7)fplot——

繪制函數(shù)圖函數(shù)fplot的調(diào)用格式:fplot(fun,lims)—

繪制函數(shù)fun在x區(qū)間lims=[xminxmax]的函數(shù)圖。fplot(fun,lims,'corline')—以指定線形繪圖。[x,y]=fplot(fun,lims)—

只返回繪圖點(diǎn)的值,而不繪圖。用plot(x,y)來(lái)繪圖。fplot('[sin(x),tan(x),cos(x)]',2*pi*[-11-11])

fplot('humps',[01],'rp')(8)ezplot—符號(hào)函數(shù)的簡(jiǎn)易繪圖函數(shù)ezplot的調(diào)用格式：ezplot(f)—這里f為包含單個(gè)符號(hào)變量x的符號(hào)表達(dá)式，在x軸的默認(rèn)范圍[-2*pi2*pi]內(nèi)繪制f(x)的函數(shù)圖ezplot(f,xmin,xmax)—給定區(qū)間ezplot(f,[xmin,xmax],figure(n))—指定繪圖窗口繪圖。ezplot('sin(x)')

ezplot('sin(x)','cos(y)',[-4*pi4*pi],figure(2))繪圖：>>ezplot('sin(x)^2+cos(y)^2-1')例sin2(x)+con2(y)=1隱函數(shù)繪制及應(yīng)用隱函數(shù)：例3.1.2其他二維圖形繪制語(yǔ)句bar––––繪制直方圖compass––––復(fù)數(shù)向量圖(羅盤圖)feather––––復(fù)數(shù)向量投影圖(羽毛圖)hist––––繪制統(tǒng)計(jì)直方圖polar––––繪制極坐標(biāo)圖stairs––––繪制階梯圖comet––––繪制彗星曲線

errorbar––––繪制誤差棒圖

fill––––繪制二維多邊形圖

loglog––––繪制對(duì)數(shù)坐標(biāo)圖

quiver––––向量場(chǎng)圖

stem––––繪制火柴桿圖area––––區(qū)域圖

rose––––繪制統(tǒng)計(jì)扇形圖(玫瑰)pie––––餅圖

convhull––––凸殼圖

scatter––––離散點(diǎn)圖fill：繪制二維多邊形并填充顏色例：x=[12345];y=[41514];fill(x,y,'r')例，繪制階梯曲線x=0:pi/20:2*pi;y=sin(x);stairs(x,y)例：階梯繪圖h2=[11;1-1];h4=[h2h2;h2-h2];h8=[h4h4;h4-h4];t=1:8;subplot(8,1,1);stairs(t,h8(1,:));axis('off')subplot(8,1,2);stairs(t,h8(2,:));axis('off')subplot(8,1,3);stairs(t,h8(3,:));axis('off')subplot(8,1,4);stairs(t,h8(4,:));axis('off')subplot(8,1,5);stairs(t,h8(5,:));axis('off')subplot(8,1,6);stairs(t,h8(6,:));axis('off')subplot(8,1,7);stairs(t,h8(7,:));axis('off')subplot(8,1,8);stairs(t,h8(8,:));axis('off')h2=[11;1-1];h4=[h2h2;h2-h2];h8=[h4h4;h4-h4];t=1:8;fori=1:8subplot(8,1,i);stairs(t,h8(i,:))axis('off')end例：繪制極坐標(biāo)繪圖t=0:2*pi/90:2*pi;y=cos(4*t);polar(t,y)例：繪制火柴桿繪圖

t=0:0.2:2*pi;y=cos(t);stem(y)例：繪制直方圖

t=0:0.2:2*pi;y=cos(t);bar(y)例：繪制彗星曲線圖t=-pi:pi/500:pi;y=tan(sin(t))-sin(tan(t));comet(t,y)x=magic(6);area(x)x=[1234567];y=[0001000];pie(x,y)pie(x,y,{'North','South','East','West','middle','fa','white'})loadseamount

scatter(x,y,50,z)a=rand(200,1);b=rand(200,1);c=rand(200,1);

scatter(a,b,100,c,'p')3.2三維繪圖三維繪圖的主要功能：繪制三維線圖繪制等高線圖繪制偽彩色圖繪制三維網(wǎng)線圖繪制三維曲面圖、柱面圖和球面圖繪制三維多面體并填充顏色3.2.1三維線圖plot3——

基本的三維圖形指令調(diào)用格式：plot3(x,y,z)——x,y,z是長(zhǎng)度相同的向量plot3(X,Y,Z)——X,Y,Z是維數(shù)相同的矩陣plot3(x,y,z,s)——

帶開(kāi)關(guān)量plot3(x1,y1,z1,'s1',x2,y2,z2,'s2',…)二維圖形的所有基本特性對(duì)三維圖形全都適用。定義三維坐標(biāo)軸大小

axis([xminxmaxyminymaxzminzmax

])gridon(off)繪制三維網(wǎng)格

text(x,y,z,‘string’)三維圖形標(biāo)注子圖和多窗口也可以用到三維圖形中例：繪制三維線圖t=0:pi/50:10*pi;plot3(t,sin(t),cos(t),'r:')3.2.2三維餅圖

pie3([43689])3.2.3三維多邊形fill3=fill——

三維多邊形的繪制和填色與二維多邊形完全相同調(diào)用格式：

fill3(x,y,z,‘s’)——

與二維相同例：用隨機(jī)頂點(diǎn)坐標(biāo)畫出5個(gè)粉色的三角形，并用黃色的○表示頂點(diǎn)y1=rand(3,5);y2=rand(3,5);y3=rand(3,5);fill3(y1,y2,y3,'m');holdon;plot3(y1,y2,y3,'yo')3.2.4三維網(wǎng)格圖mesh——

三維網(wǎng)線繪圖函數(shù)調(diào)用格式：

mesh(z)

——

z為n×m的矩陣，x與y

坐標(biāo)為元素的下標(biāo)

mesh(x,y,z)——x,y,z分別為三維空間的坐標(biāo)位置例，矩陣的三維網(wǎng)線圖z=rand(6);0.99610.77490.80010.13610.62210.12330.07820.81730.43140.86930.35100.18390.44270.86870.91060.57970.51320.24000.10670.08440.18180.54990.40180.41730.96190.39980.26380.14500.07600.04970.00460.25990.14550.85300.23990.9027z=round(z)

111010010100011110000100100000000101mesh(z)例：8階hadamard矩陣的網(wǎng)線圖h2=[11;1-1];h4=[h2h2;h2-h2]h8=[h4h4;h4-h4]111111111-11-11-11-111-1-111-1-11-1-111-1-111111-1-1-1-11-11-1-11-1111-1-1-1-1111-1-11-111-1mesh(h8)三維網(wǎng)線圖作圖要領(lǐng)生成坐標(biāo)——[X,Y]=meshgrid(x,y)表達(dá)式點(diǎn)運(yùn)算——Z=X.^2+Y.^2x——1×n的矩陣y——1×m的矩陣X，Y是n×m的矩陣，維數(shù)可任定X——nY——mZ——默認(rèn)方位角：37。5o，俯角30o。m×n維meshgrid——

網(wǎng)線坐標(biāo)值計(jì)算函數(shù)

z=f(x,y)—根據(jù)x,y坐標(biāo)找出z的高度例：繪制z=x2+y2的三維網(wǎng)線圖形

x=-5:5;y=x;[X,Y]=meshgrid(x,y)Z=X.^2+Y.^2

mesh(X,Y,Z)坐標(biāo)矩陣坐標(biāo)向量縱坐標(biāo)矩陣?yán)L圖函數(shù)colormap([R,G,B])——

色圖設(shè)定函數(shù)matlab的顏色數(shù)據(jù)集合為紅、綠、藍(lán)三顏色矩陣[R，G，B]，維數(shù)m×3r,g,b在[01]區(qū)間連續(xù)取值，理論上顏色種類可達(dá)無(wú)窮多種matlab使用三維向量表示一種顏色，常用顏色數(shù)據(jù)見(jiàn)下表飽和色

[000]—黑色

[001]—蘭色

[010]—綠色

[011]—淺蘭

[100]—紅色

[101]—粉紅

[110]—黃色

[111]—白色調(diào)和色[]—

灰色[0.500]—

暗紅色[10.620.4]—

銅色[0.4910.8]—

淺綠[0.4910.83]—寶石蘭3.2.5三維曲面圖surf——

三維曲面繪圖函數(shù)，與網(wǎng)格圖看起來(lái)一樣與三維網(wǎng)線圖的區(qū)別：網(wǎng)線圖：線條有顏色，空擋是黑色的（無(wú)顏色）曲面圖：線條是黑色的，空擋有顏色（把線條之間的空擋填充顏色，沿z軸按每一網(wǎng)格變化）調(diào)用格式：surf(x,y,z)——

繪制三維曲面圖，x,y,z

為圖形坐標(biāo)向量例：[X,Y,Z]=peaks(30)※peaks為matlab自動(dòng)生成的三維測(cè)試圖形surf(X,Y,Z)

surfc(X,Y,Z)—帶等高線的曲面圖

[X,Y,Z]=peaks(30);surfc(X,Y,Z)surfl(X,Y,Z)——

被光照射帶陰影的曲面圖[X,Y,Z]=peaks(30);surfl(X,Y,Z)

cylinde(r,n)—三維柱面繪圖函數(shù)

r為半徑；n為柱面圓周等分?jǐn)?shù)例：繪制三維陀螺錐面t1=0:0.1:0.9;t2=1:0.1:2;r=[t1-t2+2];[x,y,z]=cylinder(r,30);surf(x,y,z);grid

為球面等分?jǐn)?shù)，缺省為20例：繪制三維球面[x,y,z]=sphere(30);surf(x,y,z);3.2.6

圖形修飾方法圖形顏色的修飾matlab有極好的顏色表現(xiàn)功能，其顏色數(shù)據(jù)又構(gòu)成了一維新的數(shù)據(jù)集合，也可稱為四維圖形colormap(MAP)——

色圖設(shè)定函數(shù)，MAP為m×3維色圖矩陣圖形顏色可根據(jù)需要任意生成，也可用matlab配備的色圖函數(shù)matlab的色圖函數(shù)：

hsv——飽和值色圖

gray——線性灰度色圖

hot——暖色色圖

cool——冷色色圖

bone——蘭色調(diào)灰色圖

copper——銅色色圖

pink——粉紅色圖

prism——光譜色圖

jet——飽和值色圖IIflag——紅、白、藍(lán)交替色圖shadingfaceted—網(wǎng)格修飾，缺省方式shadingflat——

去掉黑色線條，根據(jù)小方塊的值確定顏色shadinginterp——

顏色整體改變，根據(jù)小方塊四角的值差補(bǔ)過(guò)度點(diǎn)的值確定顏色peaks(30);shadinginterp;colormap(hot)[X,Y,Z]=peaks(30);surfl(X,Y,Z)

shadinginterp;colormap(cool);axisoffpeaks(30);colormap(hot);colorbar('horiz')figure(2);colormap(cool);圖形效果修飾透視與消隱——

用于網(wǎng)線圖裁剪修飾——用于網(wǎng)線圖、曲面圖視角修飾——

觀察不同角度的三維視圖其它修飾：

a.水線修飾

b.等高線修飾透視與消隱p=peaks(30);mesh(p);hiddenonp=peaks(30);mesh(p);hiddenoff裁減修飾p=peaks;p(30:40,20:30)=nan*p(30:40,20:30);surf(p)裁減修飾p=peaks;p(30:40,20:30)=nan*p(30:40,20:30);surf(p)p=peaks;p(30:40,20:30)=nan*p(30:40,20:30);mesh(peaks,p)視角修飾(函數(shù)view(az,el))az----方位角；el----俯視角省缺值為：az=-37.5;el=30例：觀察不同視角的波峰圖形z=peaks(40);subplot(2,2,1);mesh(z);subplot(2,2,2);mesh(z);view(-15,6