程序學習前篇學習好一門程序語言的建議：掌握基本的語法結構，熟悉軟件使用環境；養成養好的編程習慣，如文件的命名，函數的命名，程序的注釋；培養良好的“邏輯結構”；數學多動手，不怕出問題，重點學會尋找問題的方法Matlab技術與應用2012-2013學年課程安排Matlab的基本知識基于Matlab語言進行實例開發Matlab技術與應用2012-2013學年第一講Matlab概述Matlab的由來NewMexico大學計算機系，系主任CleverMoler博士，1980年Matlab：MatrixLaboratory，Fortran1984年，MathWorks公司，CMatlab技術與應用2012-2013學年第一講Matlab概述Matlab的特點“智能化”的語言，一種直譯式的高級語言，語法寬松，設計自由大；運算符、函數庫豐富；圖像功能強大；功能全面，應用廣泛，現已成為大學工科類專業的必修課。Matlab技術與應用2012-2013學年第一講Matlab概述Matlab的功能數據運算：矩陣為基本單位，豐富的數值計算函數，高質量的數值計算功能符號運算：解析解，結合符號計算語言Maple繪圖功能：二維、三維圖形語言體系：程序結構控制、函數調用、數據結構、輸入輸出等工具箱：基本和可選工具箱Matlab技術與應用2012-2013學年第一講Matlab概述Matlab工具箱Matlab技術與應用2012-2013學年MATLAB主工具箱符號數學工具箱SIMULINK仿真工具箱控制系統工具箱信號處理工具箱圖象處理工具箱通訊工具箱系統辨識工具箱神經元網絡工具箱金融工具箱Matlab在工程中的應用Matlab技術與應用2012-2013學年應用系統的開發流程：算法設計與驗證算法在硬件系統的實現Matlab是首選匯編語言，C語言，C++Matlab也可與C,C++混合編程Matlab在工程中的應用Matlab技術與應用2012-2013學年Matlab在工程中的應用Matlab技術與應用2012-2013學年第一講Matlab概述Matlab的開放性與兼容性Matlab技術與應用2012-2013學年Matlab中的所有函數均對用戶開放，用戶可以根據自己的需要對函數進行查看、更改；Matlab語言可以與其它語言兼容，如C語言，C++，跨平臺兼容。(此處用Amigo平臺舉例說明)。Matlab安裝Matlab技術與應用2012-2013學年MatlabR2008a，Matlab7.6.0第一步：第二步：Matlab使用Matlab技術與應用2012-2013學年Matlab使用Matlab技術與應用2012-2013學年Matlab使用Matlab技術與應用2012-2013學年Matlab使用Matlab技術與應用2012-2013學年Matlab使用Matlab技術與應用2012-2013學年Matlab使用Matlab技術與應用2012-2013學年Matlab使用Matlab技術與應用2012-2013學年Matlab使用Matlab技術與應用2012-2013學年Matlab編程實例Matlab技術與應用2012-2013學年兩種編程方式舉一個加減法的例子說明。在調試過程中，介紹whos,who,help,clc,clear等命令符。示范matlab的安裝與打開。第二講Matlab數據及其運算-2.1數據的特點Matlab技術與應用2012-2013學年Matlab的基本數據對象：矩陣向量：一行或一列的矩陣單個數據標量：單個元素的矩陣Matlab技術與應用2012-2013學年Matlab的數據類型：雙精度，Double，64位，8個字節

單精度，Single，32位，4個字節帶符號整數，int無符號整數，uint字符數據，Char邏輯數據：沒有專門的邏輯數據，‘1’為‘真’，‘0’為‘假’第二講Matlab數據及其運算-2.1數據的特點Matlab技術與應用2012-2013學年矩陣中元素的數據類型應該一致結構體（Structure）和單元（Cell）數據類型第二講Matlab數據及其運算-2.1數據的特點Matlab技術與應用2012-2013學年變量命名變量-內存單元-命名以訪問變量以字母開頭，后接字母、數字、下劃線，最多63字符Speed_12,Var1合法字符12Speed,_Var1非法字符區分大小寫：Var與var不同Matlab中提供的函數名和命令都是小寫clc,clear,who,whos第二講Matlab數據及其運算-2.2變量及其操作Matlab技術與應用2012-2013學年預定義變量ans計算結果的默認賦值變量pi圓周率的近似值i,j虛數單位inf,Inf無窮大1/0NaN,nan非數，0/0,inf/inf函數命名時，應盡量避開第二講Matlab數據及其運算-2.2變量及其操作Matlab技術與應用2012-2013學年數據的保存：方法一save——

將工作空間中所有的變量存到matlab.mat文件中第二講Matlab數據及其運算-2.2變量及其操作Matlab技術與應用2012-2013學年數據的保存方法二save文件名存儲路徑：默認/指定saved:\data2aascii省略時以二進制格式存儲把所有數據存到指定的文件中第二講Matlab數據及其運算-2.2變量及其操作Matlab技術與應用2012-2013學年數據的保存方法三save文件名變量名ascii存儲路徑：默認/指定saved:\data2aascii省略時以二進制格式存儲文件名可以不帶擴展名，默認‘.mat’格式load文件名變量名ascii第二講Matlab數據及其運算-2.2變量及其操作Matlab技術與應用2012-2013學年數據的讀入loadload文件名load文件名變量名ascii第二講Matlab數據及其運算-2.2變量及其操作Matlab技術與應用2012-2013學年矩陣的建立A=[123;456],直接建立法利用M文件建立大矩陣C=[A;B]C=[AB]矩陣的修改A(,)=來修改第二講Matlab數據及其運算-2.3矩陣的表示

Matlab技術與應用2012-2013學年冒號表達式e1:e2:e3,e1為初始值，e2為步長（省略時步長為1），e3為終止值矩陣的拆分A(:,j),A(i,:),A(i,j)A(:,K:K+M),A(K:K+M,:)矩陣中元素的刪除A(:,[12])=[]第二講Matlab數據及其運算-2.3矩陣的表示

Matlab技術與應用2012-2013學年矩陣加減運算A+B,兩者的大小要相同；矩陣也可以和標量加減,標量與矩陣的所有元素分別進行加減操作。舉例第二講Matlab數據及其運算-2.3數據的運算

Matlab技術與應用2012-2013學年矩陣乘法C=A*B，A為m*n矩陣，B為n*p矩陣C為m*p矩陣矩陣也可以和標量相乘,標量與矩陣的所有元素分別進行乘法操作。第二講Matlab數據及其運算-2.3數據的運算

Matlab技術與應用2012-2013學年矩陣除法\和/表示左除和右除標量運算：兩者一樣1\2,2/1矩陣不同A=[1,2,3;4,2,6;7,4,9];B=[4,3,2;7,5,1;12,7,92]A\B不等于B/A第二講Matlab數據及其運算-2.3數據的運算

Matlab技術與應用2012-2013學年矩陣乘方A^x每個元素的乘方矩陣點運算點運算表示矩陣的元素的進行相關的運算兩者的維數要相同A=[12;34],B=[11;-1-1]A.*B不等于A*B第二講Matlab數據及其運算-2.3數據的運算

Matlab技術與應用2012-2013學年矩陣的其他運算inv——矩陣求逆det——行列式的值eig——矩陣的特征值diag——對角矩陣’——矩陣轉置sqrt——矩陣開方第二講Matlab數據及其運算-2.3數據的運算

Matlab技術與應用2012-2013學年Name='ShandongUniversity‘Name(1:3)第二講Matlab數據及其運算-2.4字符串

Matlab技術與應用2012-2013學年結構矩陣的建立a(1).x1='Tom';a(1).x2='[001,99]'a(2).x1='Devi';a(2).x2='[002,69]'結構成員的修改a(1).x2='[001,80]'第二講Matlab數據及其運算-2.5結構數據

Matlab技術與應用2012-2013學年單元數據的建立b={10,'liu';12,'wang'}單元數據的引用b{1,2}第二講Matlab數據及其運算-2.6單元數據

Matlab技術與應用2012-2013學年結構數據和單元數據的比較相同點：兩者均可把不同屬性的數據放在一個變量中不同點：兩者對成員的引用方式不同

結構變量名.成員名第二講Matlab數據及其運算Matlab技術與應用2012-2013學年Matlab語言把多項式表達成一個行向量，該向量中的元素是按多項式降冪排列的。

f(x)=anxn+an-1xn-1+……+loa0可用行向量p=[anan-1

……a1+a0]表示第二講Matlab數據及其運算-2.7多項式運算Matlab技術與應用2012-2013學年conv例子：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)c=4.0013.0028.0027.0018.00p=poly2str(c,'x')p=4x^4+13x^3+28x^2+27x+18第二講Matlab數據及其運算-2.7多項式乘法Matlab技術與應用2012-2013學年deconv例子：a=[123];c=[4.0013.0028.0027.0018.00][d,r]=deconv(c,a)d:c除a后的整數r:余數第二講Matlab數據及其運算-2.7多項式除法Matlab技術與應用2012-2013學年ployder例子：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+4第二講Matlab數據及其運算-2.7多項式微分Matlab技術與應用2012-2013學年[]：產生空矩陣zeros(m,n):產生全0矩陣zeros(size(A))ones(m,n);產生全1矩陣eye:單位矩陣rand:產生0~1間均勻分布的隨機矩陣randn:產生均值為0，方差為1的標準分布隨機矩陣第三講Matlab矩陣分析與處理-3.1特殊矩陣

Matlab技術與應用2012-2013學年例子1：在區間[20,50]內均勻分布的5階隨機矩陣x=20+(50-20)*rand(5)例子2:均值為0.6，方差為0.1的5階正態分布隨機矩陣。y=0.6+sqrt(0.1)*randn(5)第三講Matlab矩陣分析與處理-3.1特殊矩陣

Matlab技術與應用2012-2013學年矩陣求逆運算：對于一個方陣A，如果存在一個與其同階的方陣B，使得：A*B=B*A=I(I為單位矩陣)則稱B為A的逆矩陣，當然A也為B的逆矩陣Matlab中的調用函數：inv(A)實例：A=[1,-1,1;5,-4,3;2,1,1],求A的逆矩陣B第三講Matlab矩陣分析與處理

3.2線性方程求解

Matlab技術與應用2012-2013學年廣義逆矩陣：若矩陣A不是一個方陣，或者A是一個非滿秩的方陣時，矩陣A沒有逆矩陣，此時可以找到一個與A的轉置矩陣同型的矩陣B，使得：A*B*A=AB*A*B=B此時，矩陣B為矩陣A的偽逆，也稱為廣義逆矩陣。Matlab中的調用函數：pinv(A)實例：A=[3,1,1,1;1,3,1,1;1,1,3,1]求A的廣義逆矩陣第三講Matlab矩陣分析與處理

3.2線性方程求解

Matlab技術與應用2012-2013學年線性方程組的矩陣表示第三講Matlab矩陣分析與處理

3.2線性方程求解

Matlab技術與應用2012-2013學年線性方程組的矩陣表示第三講Matlab矩陣分析與處理

3.2線性方程求解

a為m*n的系數矩陣。m=n,該方程為“恰定”方程m>n,該方程為“超定”方程m<n,該方程為“欠定”方程Matlab技術與應用2012-2013學年第三講Matlab矩陣分析與處理

3.2線性方程求解

恰定方程求解方程ax+b(a為非奇異)x=a-1

b兩種解:x=inv(a)b—采用求逆運算解方程x=a\b—采用左除運算解方程方程ax=ba=[12;23];b=[8;13]x=inv(a)*b

x=a\b

=

ax=b例:x1+2x2=82x1+3x2=13Matlab技術與應用2012-2013學年第三講Matlab矩陣分析與處理

3.2線性方程求解

方程ax=b,m>n時此時不存在唯一解。方程解(a'a)x=a'bx=(a'

a)-1

a'b——求逆法x=a\bmatlab用最小二乘法找一個準確地基本解。第三講Matlab矩陣分析與處理

3.2線性方程求解

超定方程求解Matlab技術與應用2012-2013學年例:x1+2x2=12x1+3x2=23x1+4x2=3a=[12;23;34];b=[1;2;3];

解1x=a\b

解2x=inv(a'*a)*a'*b

=

ax=bMatlab技術與應用2012-2013學年第三講Matlab矩陣分析與處理

3.2線性方程求解

3.欠定方程組的解當方程數少于未知量個數時,即不定情況,有無窮多個解存在。matlab可求出兩個解：（1）用除法求的解x是具有最多零元素的解（2）是具有最小長度或范數的解，這個解是基于偽逆pinv求得的。第三講Matlab矩陣分析與處理

3.2線性方程求解

欠定方程求解Matlab技術與應用2012-2013學年x1+2x2+3x3=12x1+3x2+4x3=2a=[123;234];b=[1;2];x=a\bx=pinv(a)bx=x=1.000.8300.330-0.17=ax=b第三講Matlab矩陣分析與處理

3.2線性方程求解

Matlab技術與應用2012-2013學年六、微分方程求解微分方程求解的仿真算法有多種，常用的有Euler（歐拉法）、RungeKutta(龍格-庫塔法。Euler法稱一步法，用于一階微分方程。龍格-庫塔法：實際上取兩點斜率的平均斜率來計算的，其精度高于歐拉算法。第三講Matlab矩陣分析與處理

3.3微分方程求解

Matlab技術與應用2012-2013學年六、微分方程求解實例：第三講Matlab矩陣分析與處理

3.3微分方程求解

Matlab技術與應用2012-2013學年六、微分方程求解第三講Matlab矩陣分析與處理

3.4多項式擬合多項式擬合x0=0:0.1:1；y0=[-.4471.9783.115.255.024.664.014.583.455.359.22];p=polyfit(x0,y0,3)p=56.6915-87.117440.0070-0.9043xx=0:0.01:1;yy=polyval(p,xx);figure;plot(xx,yy,'-b',x0,y0,'or')Matlab技術與應用2012-2013學年謝謝大家！Matlab技術與應用2012-2013學年六、微分方程求解第四講Matlab繪圖Matlab語言豐富的圖形表現方法，使得數計算結果可以方便的、多樣性地實現了可視化，這是其它語言所不能比擬的。Matlab技術與應用2012-2013學年一、二維繪圖（一）plot——最基本的二維圖形指令plot的功能：用直線連接相鄰兩數據點來繪制圖形根據圖形坐標大小自動縮擴坐標軸，將數據標尺及單位標注自動加到兩個坐標軸上，可自定坐標軸，可把x,y軸用對數坐標表示Matlab技術與應用2012-2013學年如果已經存在一個圖形窗口，plot命令則清除當前圖形，繪制新圖形可單窗口單曲線繪圖；可單窗口多曲線繪圖；可單窗口多曲線分圖繪圖；可多窗口繪圖可任意設定曲線顏色和線型可給圖形加坐標網線和圖形加注功能一、二維繪圖Matlab技術與應用2012-2013學年plot的調用格式

plot(x)——缺省自變量繪圖格式，x為向量,以x元素值為縱坐標，以相應元素下標為橫坐標繪圖plot(x,y)——基本格式，以y(x)的函數關系作出直角坐標圖，如果y為n×m的矩陣，則以x為自變量，作出m條曲線plot(x1,y1,x2,y2)——多條曲線繪圖格式Matlab技術與應用2012-2013學年plot(x,y,’s’)——開關格式，開關量字符串s設定曲線顏色和繪圖方式，使用顏色字符串的前1~3個字母，如yellow—yel表示等。或plot(x1,y1,’s1’,x2,y2,’s2’，…)Matlab技術與應用2012-2013學年S的標準設定值如下：

字母顏色標點線型

y黃色·

點線

m粉紅○圈線

c亮藍××線

r大紅＋＋字線

g綠色－實線

b藍色星形線

w白色：虛線

k黑色－·(--)點劃線Matlab技術與應用2012-2013學年1.單窗口單曲線繪圖例1：x=[0,0.48,0.84,1,0.91,0.6,0.14][x1,x2,x3,x4,x5,x6,x7,]plot(x)Matlab技術與應用2012-2013學年2.單窗口多曲線繪圖例2：t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y,t,y1,t,y2)yy1y2Matlab技術與應用2012-2013學年3.單窗口多曲線分圖繪圖subplot——子圖分割命令調用格式：

subplot(m,n,p)——按從左至右，從上至下排列

Matlab技術與應用2012-2013學年行列繪圖序號subplot(1,3,1);plot(t,y)

subplot(1,3,2);plot(t,y1)

subplot(1,3,3);plot(t,y2)Matlab技術與應用2012-2013學年subplot(3,1,1);

plot(t,y)

subplot(3,1,2);

plot(t,y1)

subplot(3,1,3);

plot(t,y2)

Matlab技術與應用2012-2013學年4.多窗口繪圖

figure(n)——

創建窗口函數，n為窗口順序號。t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y)——

自動出現第一個窗口figure(2)plot(t,y1)——

在第二窗口繪圖figure(3)plot(t,y2)——在第三窗口繪圖Matlab技術與應用2012-2013學年5.可任意設置顏色與線型例4(例2)：plot(t,y,'r-',t,y1,'g:',t,y2,'b*')Matlab技術與應用2012-2013學年6.圖形加注功能將標題、坐標軸標記、網格線及文字注釋加注到圖形上，這些函數為：

title——

給圖形加標題xlable——

給x軸加標注ylable——

給y軸加標注text——

在圖形指定位置加標注gtext——

將標注加到圖形任意位置gridon(off)——

打開、關閉坐標網格線legend——

添加圖例axis——

控制坐標軸的刻度Matlab技術與應用2012-2013學年例：t=0:0.1:10;y1=sin(t);y2=cos(t);figure;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('時間t'),ylabel('正弦、余弦')gridMatlab技術與應用2012-2013學年

axis的用法還有：

axis([xminxmaxyminymax])——

用行向量中給出的值設定坐標軸的最大和最小值。如axis([-2205])axis(equal)——

將兩坐標軸設為相等axison(off)——

顯示和關閉坐標軸的標記、標志axisauto——

將坐標軸設置返回自動缺省值Matlab技術與應用2012-2013學年（二）fill––––

基本二維繪圖函數

fill的功能：繪制二維多邊形并填充顏色例：x=[12345];y=[41514];fill(x,y,'r')Matlab技術與應用2012-2013學年（三）特殊二維繪圖函數bar––––

繪制直方圖polar––––

繪制極坐標圖

hist––––

繪制統計直方圖

stairs––––

繪制階梯圖

stem––––

繪制火柴桿圖

rose––––

繪制統計扇形圖

comet––––

繪制彗星曲線Matlab技術與應用2012-2013學年errorbar––––

繪制誤差棒圖compass––––

復數向量圖(羅盤圖)feather––––

復數向量投影圖(羽毛圖)

quiver––––

向量場圖

area––––

區域圖

pie––––

餅圖

convhull––––

凸殼圖

scatter––––

離散點圖Matlab技術與應用2012-2013學年例，繪制階梯曲線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例：繪制極坐標繪圖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')二、三維繪圖三維繪圖的主要功能：繪制三維線圖繪制等高線圖繪制偽彩色圖繪制三維網線圖繪制三維曲面圖、柱面圖和球面圖繪制三維多面體并填充顏色（一）三維線圖plot3——

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

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

axis([xminxmaxyminymaxzminzmax])gridon(off)繪制三維網格

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

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

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

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

與二維相同例：用隨機頂點坐標畫出5個粉色的三角形，并用黃色的○表示頂點y1=rand(3,5);y2=rand(3,5);y3=rand(3,5);fill3(y1,y2,y3,'m');holdon;plot3(y1,y2,y3,'yo')（四）三維網格圖mesh——

三維網線繪圖函數調用格式：

mesh(z)——z為n×m的矩陣，x與y

坐標為元素的下標

mesh(x,y,z)——x,,y,z分別為三維空間的坐標位置例，矩陣的三維網線圖z=rand(6);

0.88080.07290.41680.76940.37750.47760.33810.71010.09640.63520.38260.70860.18950.87910.67470.89650.68760.23800.74310.35940.56260.87840.12170.39100.71890.08990.81300.48650.07680.97590.87920.16100.87820.11310.14330.6288z=round(z)

100100010101011110101100101001101001mesh(z)；mesh(z)例：8階hadamard矩陣的網線圖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)三維網線圖作圖要領生成坐標——[X,Y]=meshgrid(x,y)表達式點運算——Z=X.^2+Y.^2X，Y是n×m的矩陣，維數可任定X—nY—mZ—默認方位角：37。5o，俯角30o。n×m維meshgrid——

網線坐標值計算函數

z=f(x,y)—

根據x,y坐標找出z的高度例：繪制z=x2+y2的三維網線圖形

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

mesh(X,Y,Z)坐標矩陣坐標向量縱坐標矩陣繪圖函數colormap([R,G,B])——

色圖設定函數matlab的顏色數據集合為紅、綠、蘭三顏色矩陣[R，G，B]，維數m×3r,g,b在[01]區間連續取值，理論上顏色種類可達無窮多種matlab使用三維向量表示一種顏色，常用顏色數據見下表飽和色

[000]—

黑色

[001]—

蘭色

[010]—

綠色

[011]—

淺蘭

[100]—

紅色

[101]—

粉紅

[110]—

黃色

[111]—

白色

調和色[0.50.505]—

灰色[0.500]—

暗紅色[10.620.4]—

銅色[0.4910.8]—

淺綠[0.4910.83]—寶石蘭（五）三維曲面圖surf——

三維曲面繪圖函數，與網格圖看起來一樣與三維網線圖的區別：網線圖：線條有顏色，空擋是黑色的（無顏色）曲面圖：線條是黑色的，空擋有顏色（把線條之間的空擋填充顏色，沿z軸按每一網格變化調用格式：surf(x,y,z)——

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

為圖形坐標向量例：[X,Y,Z]=peaks(30)※peaks為matlab自動生成的三維測試圖形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)—

三維柱面繪圖函數

r為半徑；n為柱面圓周等分數例：繪制三維陀螺錐面t1=0:0.1:0.9;