《MATLAB實(shí)踐教程》課件第3章

3.1編程繪圖步驟

3.2二維繪圖

3.3三維繪圖

第3單元MATLAB基本繪圖

編程繪圖步驟如下：

(1)點(diǎn)擊目錄下拉列表框或直接在目錄文本框內(nèi)鍵入文件夾全路徑，選定你的工作文件夾，不選則缺省文件夾是work，如圖3-1所示。3.1編程繪圖步驟圖3-1選定自己的工作文件夾

(2)點(diǎn)擊工具條上的New按鈕或執(zhí)行File_New_M-file菜單命令，啟動(dòng)(打開)Editor窗口，如圖3-2所示。圖3-2點(diǎn)擊File菜單上的New_M-File命令(3)在Editor窗口輸入用戶的繪圖程序，如圖3-3所示。圖3-3在Editor窗口編程

(4)點(diǎn)擊Editor窗口的Save按鈕或執(zhí)行File_Save菜單命令，為程序命名并存盤(如myplot01.m)，如圖3-4所示。圖3-4給程序命名并存盤

(5)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名(如鍵入myplot01，不含擴(kuò)展名)，回車執(zhí)行并觀察程序運(yùn)行結(jié)果，如圖3-5所示。圖3-5在CommandWindow鍵入程序名并回車

(6)繪圖程序運(yùn)行時(shí)自動(dòng)打開Figure窗口并顯示繪圖程序的執(zhí)行結(jié)果，如圖3-6所示。圖3-6Figure窗口被打開并顯示繪圖程序的執(zhí)行結(jié)果

(7)在CommandWindow中直接輸入程序的各個(gè)語句，回車執(zhí)行并觀察運(yùn)行結(jié)果。與Editor窗口編程的操作方式比較。3.2.1用plot函數(shù)繪制散點(diǎn)圖、點(diǎn)線圖、曲線圖和多重點(diǎn)線圖

(1)?plot(x,y)函數(shù)以x為橫坐標(biāo)、以y為縱坐標(biāo)描點(diǎn)繪散點(diǎn)圖。設(shè)置連線的線型和顏色、點(diǎn)標(biāo)記的類型和尺寸的程序如下：3.2二維繪圖clc;closeall;clearall;

x=-pi:pi/10:3*pi;y1=cos(x);

plot(x,y1,'*k','MarkerSize',10);%繪散點(diǎn)圖

%plot(x,y1,'-*k','LineWidth',1.0,'MarkerSize',10);%繪點(diǎn)線圖

axis([-pi3*pi+1-11]);boxoff;

set(gca,'LineWidth',1,'FontSize',16,'FontName','Times');

xlabel('x','FontSize',16,'FontName','Times');

ylabel('cos(x)','FontSize',16,'FontName','Times');

legend('cos(x)');

title('[-\pi~3\pi]上余弦響應(yīng)散點(diǎn)圖');

%title('[-\pi~3\pi]上余弦響應(yīng)點(diǎn)線圖');

plot(x,y1,'*k','MarkerSize',10)的輸出：點(diǎn)標(biāo)記為類型*、尺寸10。如圖3-7所示。圖3-7散點(diǎn)圖

(2)?plot(x,y1,'-*k','LineWidth',1.0,'MarkerSize',10)的輸出：點(diǎn)標(biāo)記為類型*、尺寸10，連線類型為實(shí)線“—”、黑色k。如圖3-8所示。圖3-8點(diǎn)線圖

(3)用plot(x,y)函數(shù)繪曲線圖的程序：縮小x向量的間隔，以x為橫坐標(biāo)y為縱坐標(biāo)繪點(diǎn)，設(shè)置連線的線型“—”和顏色k，不設(shè)點(diǎn)標(biāo)記。程序如下：

clc;closeall;clearall;

x=-pi:pi/100:3*pi;y=cos(x);

plot(x,y,‘-k’,‘LineWidth’,1.0,‘MarkerSize’,10);

axis([-pi3*pi+1-11]);boxoff;

set(gca,‘LineWidth’,1,‘FontSize’,16,‘FontName’,‘Times’);

xlabel(‘x’,‘FontSize’,16,‘FontName’,‘Times’);

ylabel(‘cos(x)’,‘FontSize’,16,‘FontName’,‘Times’);

legend(‘cos(x)’);

title(‘[-\pi~3\pi]上余弦響應(yīng)曲線圖’);

程序輸出的曲線圖，如圖3-9所示。圖3-9曲線圖

(4)?plot函數(shù)繪帶標(biāo)記“+”的曲線圖。程序如下：

clc;closeall;clearall;

x1=-pi:pi/100:3*pi;y1=cos(x1);

x2=-pi:pi/10:3*pi;y2=cos(x2);

plot(x1,y1,'-k','LineWidth',2,'MarkerSize',10);

text(x2,y2,'+','FontSize',16);

axis([-pi3*pi+1-11]);boxoff;

set(gca,'LineWidth',1,'FontSize',16,'FontName','Times');

xlabel('x','FontSize',16,'FontName','Times');

ylabel('cos(x)','FontSize',16,'FontName','Times');

legend('cos(x)');

title('[-\pi~3\pi]上余弦響應(yīng)帶標(biāo)記曲線圖');

程序輸出的帶標(biāo)記曲線圖，如圖3-10所示。圖3-10帶標(biāo)記“+”的曲線圖

(5)?plot(x1,y1,’-*k’,x2,y2,’-ob’,…)函數(shù)以(x1,y1)、(x2,y2)，…繪多重點(diǎn)線圖。程序如下：

clc;closeall;clearall;

x=[32.0921.5481.446

52.4241.9421.912

72.7762.3382.238

93.052.6922.662

113.2122.8443.034

133.4583.1463.212

153.6543.4623.378

173.673.53.388

193.6843.5943.404

214.1143.9743.774

234.2284.1063.862

254.34.183.92

27NaN4.59NaN

29NaN4.646NaN];%試驗(yàn)數(shù)據(jù)，第1列為自變量，其余列為因變量

plot(x(:,1),x(:,2),'-*k',x(:,1),x(:,3),'-ok',x(:,1),x(:,4),'->k','MarkerSize',8);

axis([min(x(:,1))max(x(:,1))min(min(x(:,2:4)))max(max(x(:,2:4)))]);boxoff;

set(gca,'LineWidth',1,'FontSize',16,'FontName','Times');

xlabel('growthdays(d)','FontSize',16,'FontName','Times');

ylabel('fruitdiameter(cm)','FontSize',16,'FontName','Times');

title('番茄定株觀測(cè)的果實(shí)直徑動(dòng)態(tài)');

legend('果實(shí)1','果實(shí)2','果實(shí)3');

程序輸出的多重點(diǎn)線圖，如圖3-11所示。圖3-11多重點(diǎn)線圖3.2.2用fplot函數(shù)采樣繪圖

(1)?fplot(‘f(x)’,[lowup],space,‘str’)格式繪圖：采樣函數(shù)f(x)，自變量區(qū)間[lowup]，采樣間隔space，線型、標(biāo)記、顏色的設(shè)置字符串“str”。程序如下：

clc;closeall;clearall;

fplot('[cos(x),sin(x)]',[02*pi],1e-2,'pk');

axis([02*pi-11]);boxoff;gridon;

set(gca,'FontSize',16,'FontName','Times');

xlabel('x','FontSize',16,'FontName','Times');

ylabel('f(x)','FontSize',16,'FontName','Times');

legend('f(x)=cos(x)','f(x)=sin(x)',3);

程序輸出的采樣散點(diǎn)圖，如圖3-12所示。圖3-12采樣散點(diǎn)圖

(2)?fplot(y,[lowup],space,‘str’)格式繪圖：匿名采樣函數(shù)y=@(x)f(x)，自變量區(qū)間[lowup]，采樣間隔space，線型、標(biāo)記、顏色的設(shè)置字符串“str”。程序如下：

clc;closeall;clearall;

y=@(x)[200*sin(x)./xx.^2];

fplot(y,[-2020],1e-3);grid;boxoff;

set(gca,'LineWidth',1,'FontSize',16,'FontName','Times');

xlabel('x','FontSize',16,'FontName','Times');

ylabel('f(x)','FontSize',16,'FontName','Times');

legend('f(x)=200*sin(x)/x','f(x)=x^2',1);

程序輸出的采樣曲線圖，如圖3-13所示。圖3-13采樣曲線圖3.2.3用figure函數(shù)開多個(gè)窗口分別繪圖

用figure函數(shù)打開窗口，用plot函數(shù)在所開窗口中分別繪圖，程序如下：

clc;closeall;clearall;

x=linspace(0,2*pi,60);

y1=sin(x);y2=cos(x);y3=tan(x);y4=cot(x);

figure;plot(x,y1);boxoff;

axis([02*pi-11]);title('sin(x)');

figure;plot(x,y2);boxoff;

axis([02*pi-11]);title('cos(x)');

figure;plot(x,y3);boxoff;

axis([02*pi-4040]);title('tan(x)');

figure;plot(x,y4);boxoff;

axis([02*pi-4040]);title('cot(x)');

上面程序每開1個(gè)窗口繪出1條曲線，4個(gè)窗口分別繪4條曲線，如圖3-14所示。圖3-144個(gè)窗口顯示的4條曲線3.2.4用subplot函數(shù)開多個(gè)子窗口分別繪圖

用subplot函數(shù)打開多個(gè)子窗口，用plot函數(shù)在所開子窗口中分別繪圖，程序如下：

clc;closeall;clearall;

x=linspace(0,2*pi,60);

y1=sin(x);y2=cos(x);y3=tan(x);y4=cot(x);

subplot(2,2,1);

plot(x,y1);boxoff;

axis([02*pi-11]);title('sin(x)');

subplot(2,2,2);

plot(x,y2);boxoff;

axis([02*pi-11]);title('cos(x)');

subplot(2,2,3);

plot(x,y3);boxoff;

axis([02*pi-4040]);title('tan(x)');

subplot(2,2,4);

plot(x,y4);boxoff;

axis([02*pi-4040]);title('cot(x)');

上面程序每開1個(gè)子窗口繪出1條曲線，4個(gè)子窗口分別繪4條曲線，如圖3-15所示。圖3-154個(gè)窗口顯示的4條曲線3.2.5用loglog函數(shù)繪制雙軸對(duì)數(shù)圖

loglog函數(shù)以雙軸坐標(biāo)的對(duì)數(shù)值繪點(diǎn)，刻度顯示仍用原坐標(biāo)值。程序如下：

clc;closeall;clearall;

x=0:pi/180:2*pi;y=abs(1000*sin(4*x))+1;

loglog(x,y,':b','LineWidth',2);boxoff;axis([-1010110^3]);

set(gca,'FontSize',16,'FontName','Times');

xlabel('x','FontSize',16,'FontName','Times');

ylabel('y','FontSize',16,'FontName','Times');

legend('y=|1000sin(4x)|+1',3);

程序輸出的雙軸對(duì)數(shù)圖，如圖3-16所示。圖3-16雙軸對(duì)數(shù)圖3.2.6用semilogx函數(shù)繪制橫軸對(duì)數(shù)圖

semilogx函數(shù)以橫軸坐標(biāo)對(duì)數(shù)值、縱軸坐標(biāo)原值繪點(diǎn)，刻度顯示仍用原坐標(biāo)值。程序如下：

clc;closeall;clearall;

x=0:pi/180:2*pi;

y=abs(1000*sin(4*x))+1;

semilogx(x,y,':b','LineWidth',2);boxoff;

axis([-1010110^3]);

set(gca,'FontSize',16,'FontName','Times');

xlabel('x','FontSize',16,'FontName','Times');

ylabel('y','FontSize',16,'FontName','Times');

legend('y=|1000sin(4x)|+1',3);

程序輸出的橫軸對(duì)數(shù)圖，如圖3-17所示。圖3-17橫軸對(duì)數(shù)圖3.2.7用semilogy函數(shù)繪制縱軸對(duì)數(shù)圖

semilogy函數(shù)以橫軸坐標(biāo)原值、縱軸坐標(biāo)對(duì)數(shù)值繪點(diǎn)，刻度顯示仍用原坐標(biāo)值。程序如下：

clc;closeall;clearall;

x=0:pi/180:2*pi;

y=abs(1000*sin(4*x))+1;

semilogy(x,y,':b','LineWidth',2);boxoff;

axis([02*pi110^3]);

set(gca,'FontSize',16,'FontName','Times');

xlabel('x','FontSize',16,'FontName','Times');

ylabel('y','FontSize',16,'FontName','Times');

legend('y=|1000sin(4x)|+1',3);

程序輸出的縱軸對(duì)數(shù)圖，如圖3-18所示。圖3-18縱軸對(duì)數(shù)圖3.2.8用plotyy函數(shù)繪制雙縱軸圖

plotyy(x,y1,x,y2)函數(shù)以x為橫坐標(biāo)，以y1和y2為縱坐標(biāo)繪點(diǎn)，顯示兩縱軸。程序如下：

clc;closeall;clearall;

x=0:0.01:20;

y1=200*exp(-0.05*x).*sin(x);

y2=0.8*exp(-0.5*x).*sin(10*x);

%[AX,H1,H2]=plotyy(x,y1,x,y2,'plot');boxoff; %繪線性雙縱軸圖

%[AX,H1,H2]=plotyy(x,y1,x,y2,'semilogx');boxoff;%繪橫軸對(duì)數(shù)雙縱軸圖

[AX,H1,H2]=plotyy(x,y1,x,y2,'semilogy');boxoff; %繪縱軸對(duì)數(shù)雙縱軸圖

set(AX,'LineWidth',1,'FontSize',16,'FontName','Times');

set(H1,'LineStyle','-','LineWidth',2);

set(H2,'LineStyle','--','LineWidth',2);

xlabel('x','FontSize',16,'FontName','Times');

ylabel('f(x)','FontSize',16,'FontName','Times');

legend([H1H2],'f(x)=200e^-^0^.^0^5^xsinx','f(x)=0.8e^-^0^.^5^xsin10x',3);

title('[020]上雙縱軸圖');

%title('[020]上線性雙縱軸圖'); %線性雙縱軸圖圖題

%title('[020]上橫軸對(duì)數(shù)雙縱軸圖'); %橫軸對(duì)數(shù)雙縱軸圖圖題

%title('[020]上縱軸對(duì)數(shù)雙縱軸圖'); %縱軸對(duì)數(shù)雙縱軸圖圖題

函數(shù)plotyy(x,y1,x,y2,'plot')輸出的線性雙縱軸圖，如圖3-19所示。圖3-19線性雙縱軸圖函數(shù)plotyy(x,y1,x,y2,'semilogx')輸出的橫軸對(duì)數(shù)雙縱軸圖，如圖3-20所示。圖3-20橫軸對(duì)數(shù)雙縱軸圖函數(shù)plotyy(x,y1,x,y2,'semilogy')輸出的縱軸對(duì)數(shù)雙縱軸圖，如圖3-21所示。圖3-21縱軸對(duì)數(shù)雙縱軸圖3.2.9用polar函數(shù)繪制極坐標(biāo)圖

polar(theta,rho,'-k')函數(shù)以theta為極角、以rho為極徑繪制極坐標(biāo)圖。程序如下：

clc;closeall;clearall;

theta=0:0.01:2*pi;

rho=sin(2*theta).*cos(2*theta);

h=polar(theta,rho,'-k');

set(h,'LineWidth',2);

title('polarplot','FontSize',16);

程序輸出的極坐標(biāo)圖，如圖3-22所示。圖3-22極坐標(biāo)圖3.2.10用bar函數(shù)繪制垂直柱形圖

bar(x,y)函數(shù)以x為橫坐標(biāo)繪柱的位置、以y為縱坐標(biāo)繪柱的高度。程序如下：

clc;closeall;clearall;

x=[0123456];

y1=[2715171383];

y2=[3917161152];

bar(x,y1);boxoff; %彩色填充bar圖

%bar(x,y1,0.8,'w','LineWidth',2);boxoff; %白色填充bar圖

%bar(x,[y1'y2'],'LineWidth',2);boxoff; %分組彩色填充bar圖

set(gca,'LineWidth',1,'FontSize',16,'FontName','Times');

axis([min(x)-1max(x)+10max(y1)+1])

xlabel('x','FontSize',16,'FontName','Times');

ylabel('n(x)','FontSize',16,'FontName','Times');

title('frequencedistribution');

bar(x,y1)函數(shù)輸出橫坐標(biāo)x、柱高y1、柱寬默認(rèn)、柱面彩色填充的垂直柱形圖，如圖3-23所示。圖3-23垂直柱形圖

bar(x,y1,0.8,'w','LineWidth',2)函數(shù)輸出柱寬因子0.8、柱邊線寬2、柱面白色填充的柱形圖，如圖3-24所示。圖3-24柱面白色填充的柱形圖

bar(x,[y1'y2'],'LineWidth',2)函數(shù)輸出橫坐標(biāo)x、柱邊線寬2、柱高y1和y2的二重柱形圖，如圖3-25所示。圖3-25二重柱形圖3.2.11用hist函數(shù)繪制直方圖

將下面播種機(jī)試驗(yàn)測(cè)定的粒距樣本輸入Excel并存盤為lijusample.xls文件。注意，數(shù)據(jù)按一列輸進(jìn)Excel，第一行是變量名(如起名liju)，第二行以后是數(shù)據(jù)，參見第1單元表1-1。

9810211093961309412011095999810511484121

1089599103103113111979910810211294105103114

871011017711312998103961049799117929510099

11512191941069712011710610697108112102120107

9910411197100124961101091161067692106103120

1151051018511910285105102951059387839794

10585

(1)調(diào)用lijusample.xls數(shù)據(jù)，用hist函數(shù)統(tǒng)計(jì)粒距頻數(shù)和繪頻數(shù)分布直方圖。程序如下：

clc;closeall;clearall;

file='D:\Users\MyMATLABFiles\lijusample.xls';

y=xlsread(file,'Sheet1');

N=length(y);zushu=floor(1+3.322*log10(N));

x=linspace(min(y),max(y),zushu);

zuju=(max(y)-min(y))/(zushu-1);

hist(y,x,'w');boxoff;

axis([min(x)-zujumax(x)+zuju036]);

xlabel('粒距(mm)','FontSize',16,'FontName','Times');

ylabel('頻數(shù)','FontSize',16,'FontName','Times');

set(gca,'LineWidth',1,'FontSize',16,'FontName','Times');

h=findobj(gca,'Type','patch');

set(h,'FaceColor','w','EdgeColor','k','LineWidth',2)

title('粒距的頻數(shù)分布');

程序輸出的粒距頻數(shù)分布直方圖，如圖3-26所示。圖3-26粒距頻數(shù)分布直方圖

(2)調(diào)用lijusample.xls數(shù)據(jù)，用hist函數(shù)統(tǒng)計(jì)粒距頻數(shù)count，以組中值x為橫坐標(biāo)、以所統(tǒng)計(jì)頻數(shù)count計(jì)算頻率freq，以該頻率freq為縱坐標(biāo)，用bar(x,freq)函數(shù)繪制頻率分布直方圖。程序如下：

clc;closeall;clearall;

file='D:\Users\MyMATLABFiles\lijusample.xls';%Excel數(shù)據(jù)文件的路徑

y=xlsread(file,'Sheet1'); %讀入Excel數(shù)據(jù)文件并賦值給y

N=length(y); %計(jì)算樣本容量

zushu=floor(1+3.322*log10(N)); %計(jì)算統(tǒng)計(jì)分組個(gè)數(shù)

x=linspace(min(y),max(y),zushu); %生成組中值

zuju=(max(y)-min(y))/(zushu-1); %計(jì)算組距

count=hist(y,x);freq=count/N;%統(tǒng)計(jì)或計(jì)算所有組區(qū)間的組頻數(shù)和組頻率

bar(x,freq,1);boxoff;

axis([min(x)-zujumax(x)+zuju00.4]);

xlabel('粒距(mm)','FontSize',16,'FontName','Times');

ylabel('頻率','FontSize',16,'FontName','Times');

set(gca,'LineWidth',1,'FontSize',16,'FontName','Times');

h=findobj(gca,'Type','patch');

set(h,'FaceColor','w','EdgeColor','k','LineWidth',2)

title('粒距的頻率分布');

程序輸出的粒距頻率分布直方圖，如圖3-27所示。圖3-27粒距頻率分布直方圖3.2.12用stem函數(shù)繪制離散序列散點(diǎn)圖

(1)?stem(t,y)函數(shù)以離散時(shí)間t為橫坐標(biāo)、對(duì)應(yīng)離散序列y為縱坐標(biāo)繪制附高度線的散點(diǎn)圖。程序如下：

clc;closeall;clearall;

t=0:5:150;

y=exp(-0.02*t).*cos(0.5*t);

h=stem(t,y,'--ob');%繪制不填充離散序列散點(diǎn)圖并為圖形句柄h賦值

%h=stem(t,y,'--ob',‘fill'); %繪制填充離散序列散點(diǎn)圖并為圖形句柄h賦值

axis([0150-11]);boxoff;

set(gca,'FontSize',16);

set(h,'MarkerSize',10,'LineWidth',2);

xlabel('Timein\musecs','FontSize',16)

ylabel('Magnitude','FontSize',16)

程序輸出的不填充附高度線的離散序列散點(diǎn)圖，如圖3-28所示。圖3-28不填充附高度線的離散序列散點(diǎn)圖程序輸出的填充附高度線的離散序列散點(diǎn)圖，如圖3-29所示。圖3-29填充附高度線的離散序列散點(diǎn)圖

(2)?stem函數(shù)以離散時(shí)間向量t為橫坐標(biāo)、以對(duì)應(yīng)離散序列矩陣y為縱坐標(biāo)繪制附高度線的散點(diǎn)圖，其中y矩陣的每一列是一個(gè)離散時(shí)間序列。程序如下：

clc;closeall;clearall;

t=0:5:150; %離散時(shí)間向量

y=[sin(t);exp(-t/50)]'; %離散序列矩陣

h=stem(t,y,'--o'); %繪制不填充二重離散序列散點(diǎn)圖，并為圖形句柄h賦值

axis([0150-11]);boxoff;

set(gca,'FontSize',16);

set(h,'MarkerSize',10,'LineWidth',2);

xlabel('Timein\musecs','FontSize',16)

ylabel('Magnitude','FontSize',16)

程序輸出的二重離散序列散點(diǎn)圖，如圖3-30所示。圖3-30二重離散序列散點(diǎn)圖3.2.13用errorbar函數(shù)繪制誤差圖

根據(jù)試驗(yàn)或調(diào)查數(shù)據(jù)計(jì)算誤差，用errorbar函數(shù)繪制以0.95置信區(qū)間為誤差帶的誤差圖。程序如下：

clc;closeall;clearall;

p=[9.147.7410.53;108.5511.46;2.711.923.49;

0.2400.47;10.458.9711.94;1.550.952.15;5.364.276.45;

11.229.6912.75;3.852.924.78]; %矩陣p的3列分別為如下數(shù)據(jù)：

%手機(jī)品牌喜好百分率，喜好百分率置信下限，喜好百分率置信上限

pinpai=1:size(p,1); %為品牌編號(hào)，自1開始順序編號(hào)

E=(p(:,3)-p(:,2))/2; %計(jì)算喜好百分率偏差，即置信區(qū)間半長(zhǎng)度

h=errorbar(pinpai,p(:,1),E,'--o');boxoff;

axis([0.5max(pinpai)+0.5min(p(:,1)-E)max(p(:,1)+E)]);

set(gca,'FontSize',16);set(h,'LineWidth',2);

xlabel('品牌','FontSize',16);

ylabel('喜好百分率(%)','FontSize',16);

errorbar(pinpai,p(:,1),E,'--o')函數(shù)繪出以手機(jī)品牌編號(hào)為橫坐標(biāo)、喜好百分率為縱坐標(biāo)、0.95置信區(qū)間為誤差帶的誤差圖(點(diǎn)標(biāo)記“o”，誤差條標(biāo)記“工”，虛線線型連線“--”)如圖3-31所示。圖3-31誤差圖3.2.14用pie函數(shù)繪制二維餅圖

pie(x,explode)函數(shù)以參數(shù)x中各個(gè)元素值占總和的百分率決定扇形面積、以參數(shù)explode中非0元素的序號(hào)指定炸開的扇形塊繪餅圖。程序如下：

clc;closeall;clearall;

x=[194.8266.5230.9119.7];

explode=[0010]; %變量explode指定第3號(hào)扇形塊炸開

h=pie(x,explode);

%計(jì)算x中各個(gè)元素值占總和的百分率及繪餅圖

textObjs=findobj(h,'Type','text')

set(textObjs,'FontSize',16);

程序輸出第3號(hào)扇形塊炸開的餅圖，如圖3-32所示。圖3-32餅圖3.2.15用contour函數(shù)繪制二維等高線圖

contour(x,y,z)函數(shù)在xoy坐標(biāo)系中繪制z等間隔所取各個(gè)定值下的y=f(x)圖形。程序如下：

clc;closeall;clearall;

x=[0:0.15:2*pi];

y=[0:0.35:2*pi];

[x,y]=meshgrid(x,y);

z=(x+pi/5).^2+(y+pi/3).^2-5*x.*y+115;

[C,h]=contour(x,y,z,7,'--k','LineWidth',2);boxoff;

clabel(C,h,'FontSize',11);

set(gca,'FontSize',16);

xlabel('x','FontSize',16);

ylabel('y','FontSize',16);

title('Contourplotofmatrix')

程序輸出的等高線圖，如圖3-33所示。圖3-33等高線圖3.3.1用plot3函數(shù)繪制三維散點(diǎn)圖、點(diǎn)線圖和曲線圖

plot3(x,y,z)函數(shù)以x為橫坐標(biāo)、y為縱坐標(biāo)、z為垂直坐標(biāo)描點(diǎn)繪制z=f(x,y)圖形，通過連線線型和顏色、點(diǎn)標(biāo)記類型和尺寸等設(shè)置或不設(shè)置實(shí)現(xiàn)散點(diǎn)圖、點(diǎn)線圖和曲線圖等。3.3三維繪圖clc;closeall;clearall;

t=0:pi/30:10*pi;y1=sin(t);y2=cos(t);

plot3(y1,y2,t,'*b','MarkerSize',7,'LineWidth',2);%不設(shè)置線型繪散點(diǎn)圖

%plot3(y1,y2,t,'--b','MarkerSize',7,'LineWidth',2);%不設(shè)置點(diǎn)標(biāo)記繪曲線圖

%plot3(y1,y2,t,'-b','LineWidth',2); %不設(shè)置點(diǎn)標(biāo)記繪曲線圖

%plot3(y1,y2,t,'-*b','MarkerSize',7,'LineWidth',2);%設(shè)置線型點(diǎn)標(biāo)記繪點(diǎn)線圖

axis([-11-11010*pi]);boxoff;

set(gca,'FontSize',16,'FontName','Times');

xlabel('sin(t)','FontSize',16,'FontName','Times');

ylabel('cos(t)','FontSize',16,'FontName','Times');

zlabel('t','FontSize',16,'FontName','Times');

程序中分別采用4種格式的plot3函數(shù)，輸出圖形如圖3-34所示。(a)三維散點(diǎn)圖(b)三維虛線圖(c)三維實(shí)線圖(d)三維點(diǎn)線圖圖3-34三維圖形3.3.2用bar3函數(shù)繪制三維垂直柱形圖

bar3(x,y)函數(shù)以x為柱中心橫坐標(biāo)、矩陣y每一列為一個(gè)樣本、矩陣y元素值為柱高縱坐標(biāo)繪制分組比較的三維柱形圖。程序如下：

clc;closeall;clearall;

x=[12345];%成績(jī)值

y=[2332;139107;15161818;912911;2113];%每列為一個(gè)班級(jí)的頻數(shù)

bar3(x,y);boxoff;

xlabel('class','FontSize',16,'FontName','Times');

ylabel('record','FontSize',16,'FontName','Times');

zlabel('freq','FontSize',16,'FontName','Times');

函數(shù)輸出的分組三維柱形圖，如圖3-35所示。圖3-35分組三維柱形圖3.3.3用stem3函數(shù)繪制三維離散序列圖

stem3(x,y,z)函數(shù)以二維離散序列(x,y)為橫、縱坐標(biāo)，以離散序列z=f(x,y)為垂直坐標(biāo)繪制三維離散序列圖。程序如下：

clc;closeall;clearall;

th=(0:127)/128*2*pi; %單位循環(huán)的轉(zhuǎn)角序列，128個(gè)點(diǎn)，周期2*pi

x=cos(th); %復(fù)數(shù)實(shí)部

y=sin(th); %復(fù)數(shù)虛部

f=abs(fft(ones(10,1),128)); %傅立葉變換的振幅

stem3(x,y,f','d','fill');boxoff;

set(gca,'FontSize',16);

view([-6530])

xlabel('Real','FontSize',16)

ylabel('Imaginary','FontSize',16)

zlabel('Amplitude','FontSize',16)

title('MagnitudeFrequencyResponse')

程序輸出的三維離散序列圖，如圖3-36所示。圖3-36三維離散序列圖3.3.4用pie3函數(shù)繪制三維餅圖

pie3(x,explode)函數(shù)以參數(shù)x中各個(gè)元素值占總和的百分率決定扇形面積，以參數(shù)explode中非0元素的序號(hào)指定炸開的扇形塊繪制三維餅圖。程序如下：

clc;closeall;clearall;

x=[194.8266.5230.9119.7];

explode=[0010]; %變量explode指定第3號(hào)扇形塊炸開

h=pie3(x,explode); %計(jì)算x中各個(gè)元素值占總和的百分率及繪餅圖

textObjs=findobj(h,'Type','text')

set(textObjs,'FontSize',16);

程序輸出第3號(hào)扇形塊炸開的三維餅圖，如圖3-37所示。圖3-37三維餅圖3.3.5用mesh函數(shù)、meshc函數(shù)和meshz函數(shù)繪制三維網(wǎng)格圖

mesh函數(shù)以x(長(zhǎng)度m)為橫坐標(biāo)、以y(長(zhǎng)度n)為縱坐標(biāo)，以x與y的m?×?n個(gè)組合構(gòu)建網(wǎng)格點(diǎn)，計(jì)算網(wǎng)格點(diǎn)(x,y)上的z值，再以z為垂直坐標(biāo)在網(wǎng)格點(diǎn)上方描點(diǎn)，直線連接相鄰4個(gè)點(diǎn)形成曲面網(wǎng)格，最終繪出表達(dá)二元函數(shù)z=f(x,y)的三維網(wǎng)格圖。與mesh函數(shù)的輸入?yún)?shù)和運(yùn)行格式相同，meshc函數(shù)繪制附等高線三維網(wǎng)格圖，meshz函數(shù)繪制附高度線三維網(wǎng)格圖。程序如下：clc;closeall;clearall;

x=[0:0.15:2*pi];y=[0:0.35