1-32是：圖形應用篇

33-66是：界面設計篇

67-84是：圖形處理篇

85-100是：數值分析篇



實例1：三角函數曲線（1）

function shili01
h0=figure('toolbar','none',...
    'position',[198 56 350 300],...
    'name','實例01');
h1=axes('parent',h0,...
   'visible','off');
x=-pi:0.05:pi;
y=sin(x);
plot(x,y);
xlabel('自變量X');
ylabel('函數值Y');
title('SIN( )函數曲線');
grid on

實例2：三角函數曲線（2）

function shili02
h0=figure('toolbar','none',...
    'position',[200 150 450 350],...
    'name','實例02');
x=-pi:0.05:pi;
y=sin(x)+cos(x);
plot(x,y,'-*r','linewidth',1);
grid on
xlabel('自變量X');
ylabel('函數值Y');
title('三角函數');
實例3：圖形的疊加

function shili03
h0=figure('toolbar','none',...
    'position',[200 150 450 350],...
    'name','實例03');
x=-pi:0.05:pi;
y1=sin(x);
y2=cos(x);
plot(x,y1,...
    '-*r',...
    x,y2,...
    '--og');
grid on
xlabel('自變量X');
ylabel('函數值Y');
title('三角函數');
實例4：雙y軸圖形的繪制

function shili04
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例04');
x=0:900;a=1000;b=0.005;
y1=2*x;
y2=cos(b*x);
[haxes,hline1,hline2]=plotyy(x,y1,x,y2,'semilogy','plot');
axes(haxes(1))
ylabel('semilog plot');
axes(haxes(2))
ylabel('linear plot');
實例5：單個軸窗口顯示多個圖形

function shili05
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例05');
t=0:pi/10:2*pi;
[x,y]=meshgrid(t);
subplot(2,2,1)
plot(sin(t),cos(t))
axis equal

subplot(2,2,2)
z=sin(x)-cos(y);
plot(t,z)
axis([0 2*pi -2 2])

subplot(2,2,3)
h=sin(x)+cos(y);
plot(t,h)
axis([0 2*pi -2 2])

subplot(2,2,4)
g=(sin(x).^2)-(cos(y).^2);
plot(t,g)
axis([0 2*pi -1 1])
實例6：圖形標注

function shili06
h0=figure('toolbar','none',...
    'position',[200 150 450 400],...
    'name','實例06');
t=0:pi/10:2*pi;
h=plot(t,sin(t));
xlabel('t=0到2\pi','fontsize',16);
ylabel('sin(t)','fontsize',16);
title('\it{從 0to2\pi 的正弦曲線}','fontsize',16)
x=get(h,'xdata');
y=get(h,'ydata');
imin=find(min(y)==y);
imax=find(max(y)==y);
text(x(imin),y(imin),...
    ['\leftarrow最小值=',num2str(y(imin))],...
    'fontsize',16)
text(x(imax),y(imax),...
    ['\leftarrow最大值=',num2str(y(imax))],...
    'fontsize',16)
    
實例7：條形圖形

function shili07
h0=figure('toolbar','none',...
    'position',[200 150 450 350],...
    'name','實例07');
tiao1=[562 548 224 545 41 445 745 512];
tiao2=[47 48 57 58 54 52 65 48];
t=0:7;
bar(t,tiao1)
xlabel('X軸');
ylabel('TIAO1值');
h1=gca;
h2=axes('position',get(h1,'position'));
plot(t,tiao2,'linewidth',3)
set(h2,'yaxislocation','right','color','none','xticklabel',[])
實例8：區域圖形

function shili08
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例08');
x=91:95;
profits1=[88 75 84 93 77];
profits2=[51 64 54 56 68];
profits3=[42 54 34 25 24];
profits4=[26 38 18 15 4];
area(x,profits1,'facecolor',[0.5 0.9 0.6],...
    'edgecolor','b',...
    'linewidth',3)
hold on
area(x,profits2,'facecolor',[0.9 0.85 0.7],...
    'edgecolor','y',...
    'linewidth',3)
hold on
area(x,profits3,'facecolor',[0.3 0.6 0.7],...
    'edgecolor','r',...
    'linewidth',3)
hold on
area(x,profits4,'facecolor',[0.6 0.5 0.9],...
    'edgecolor','m',...
    'linewidth',3)
hold off
set(gca,'xtick',[91:95])
set(gca,'layer','top')
gtext('\leftarrow第一季度銷量')
gtext('\leftarrow第二季度銷量')
gtext('\leftarrow第三季度銷量')
gtext('\leftarrow第四季度銷量')
xlabel('年','fontsize',16);
ylabel('銷售量','fontsize',16); 
實例9：餅圖的繪制

function shili09
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例09');
t=[54 21 35;
    68 54 35;
    45 25 12;
    48 68 45;
    68 54 69];
x=sum(t);
h=pie(x);
textobjs=findobj(h,'type','text');
str1=get(textobjs,{'string'});
val1=get(textobjs,{'extent'});
oldext=cat(1,val1{:});
names={'商品一：';'商品二：';'商品三：'};
str2=strcat(names,str1);
set(textobjs,{'string'},str2)
val2=get(textobjs,{'extent'});
newext=cat(1,val2{:});
offset=sign(oldext(:,1)).*(newext(:,3)-oldext(:,3))/2;
pos=get(textobjs,{'position'});
textpos=cat(1,pos{:});
textpos(:,1)=textpos(:,1)+offset;
set(textobjs,{'position'},num2cell(textpos,[3,2]))
實例10：階梯圖

function shili10
h0=figure('toolbar','none',...
    'position',[200 150 450 400],...
    'name','實例10');
a=0.01;
b=0.5;
t=0:10;
f=exp(-a*t).*sin(b*t);
stairs(t,f)
hold on
plot(t,f,':*')
hold off
glabel='函數e^{-(\alpha*t)}sin\beta*t的階梯圖';
gtext(glabel,'fontsize',16)
xlabel('t=0:10','fontsize',16)
axis([0 10 -1.2 1.2])
實例11：枝干圖

function shili11
h0=figure('toolbar','none',...
    'position',[200 150 450 350],...
    'name','實例11');
x=0:pi/20:2*pi;
y1=sin(x);
y2=cos(x);
h1=stem(x,y1+y2);
hold on
h2=plot(x,y1,'^r',x,y2,'*g');
hold off
h3=[h1(1);h2];
legend(h3,'y1+y2','y1=sin(x)','y2=cos(x)')
xlabel('自變量X');
ylabel('函數值Y');
title('正弦函數與余弦函數的線性組合');
實例12：羅盤圖

function shili12
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例12');
winddirection=[54 24 65 84 
    256 12 235 62 
    125 324 34 254];
windpower=[2 5 5 3 
    6 8 12 7 
    6 14 10 8];
rdirection=winddirection*pi/180;
[x,y]=pol2cart(rdirection,windpower);
compass(x,y);
desc={'風向和風力',
    '北京氣象臺',
    '10月1日0：00到',
    '10月1日12：00'};
gtext(desc)
實例13：輪廓圖

function shili13
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例13');
[th,r]=meshgrid((0:10:360)*pi/180,0:0.05:1);
[x,y]=pol2cart(th,r);
z=x+i*y;
f=(z.^4-1).^(0.25);
contour(x,y,abs(f),20)
axis equal
xlabel('實部','fontsize',16);
ylabel('虛部','fontsize',16);
h=polar([0 2*pi],[0 1]);
delete(h)
hold on
contour(x,y,abs(f),20)
實例14：交互式圖形

function shili14
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例14');
axis([0 10 0 10]);
hold on
x=[];
y=[];
n=0;
disp('單擊鼠標左鍵點取需要的點');
disp('單擊鼠標右鍵點取最后一個點');
but=1;
while but==1
    [xi,yi,but]=ginput(1);
    plot(xi,yi,'bo')
    n=n+1;
    disp('單擊鼠標左鍵點取下一個點');
    x(n,1)=xi;
    y(n,1)=yi;
end
t=1:n;
ts=1:0.1:n;
xs=spline(t,x,ts);
ys=spline(t,y,ts);
plot(xs,ys,'r-');
hold off
實例14：交互式圖形

function shili14
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例14');
axis([0 10 0 10]);
hold on
x=[];
y=[];
n=0;
disp('單擊鼠標左鍵點取需要的點');
disp('單擊鼠標右鍵點取最后一個點');
but=1;
while but==1
    [xi,yi,but]=ginput(1);
    plot(xi,yi,'bo')
    n=n+1;
    disp('單擊鼠標左鍵點取下一個點');
    x(n,1)=xi;
    y(n,1)=yi;
end
t=1:n;
ts=1:0.1:n;
xs=spline(t,x,ts);
ys=spline(t,y,ts);
plot(xs,ys,'r-');
hold off
實例15：變換的傅立葉函數曲線

function shili15
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例15');
axis equal
m=moviein(20,gcf);
set(gca,'nextplot','replacechildren')
h=uicontrol('style','slider','position',...
    [100 10 500 20],'min',1,'max',20)
for j=1:20
    plot(fft(eye(j+16)))
    set(h,'value',j)
    m(:,j)=getframe(gcf);
end
clf;
axes('position',[0 0 1 1]);
movie(m,30)
實例16：勞倫茲非線形方程的無序活動

function shili15
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例15');
axis equal
m=moviein(20,gcf);
set(gca,'nextplot','replacechildren')
h=uicontrol('style','slider','position',...
    [100 10 500 20],'min',1,'max',20)
for j=1:20
    plot(fft(eye(j+16)))
    set(h,'value',j)
    m(:,j)=getframe(gcf);
end
clf;
axes('position',[0 0 1 1]);
movie(m,30)
實例17：填充圖

function shili17
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例17');
t=(1:2:15)*pi/8;
x=sin(t);
y=cos(t);
fill(x,y,'r')
axis square off
text(0,0,'STOP',...
    'color',[1 1 1],...
    'fontsize',50,...
    'horizontalalignment','center')
例18：條形圖和階梯形圖

function shili18
h0=figure('toolbar','none',...
    'position',[200 150 450 250],...
    'name','實例18');
subplot(2,2,1)
x=-3:0.2:3;
y=exp(-x.*x);
bar(x,y)
title('2-D Bar Chart')

subplot(2,2,2)
x=-3:0.2:3;
y=exp(-x.*x);
bar3(x,y,'r')
title('3-D Bar Chart')

subplot(2,2,3)
x=-3:0.2:3;
y=exp(-x.*x);
stairs(x,y)
title('Stair Chart')

subplot(2,2,4)
x=-3:0.2:3;
y=exp(-x.*x);
barh(x,y)
title('Horizontal Bar Chart') 
實例19：三維曲線圖

function shili19
h0=figure('toolbar','none',...
    'position',[200 150 450 400],...
    'name','實例19');
subplot(2,1,1)
x=linspace(0,2*pi);
y1=sin(x);
y2=cos(x);
y3=sin(x)+cos(x);
z1=zeros(size(x));
z2=0.5*z1;
z3=z1;
plot3(x,y1,z1,x,y2,z2,x,y3,z3)
grid on
xlabel('X軸');
ylabel('Y軸');
zlabel('Z軸');
title('Figure1:3-D Plot')

subplot(2,1,2)
x=linspace(0,2*pi);
y1=sin(x);
y2=cos(x);
y3=sin(x)+cos(x);
z1=zeros(size(x));
z2=0.5*z1;
z3=z1;
plot3(x,z1,y1,x,z2,y2,x,z3,y3)
grid on
xlabel('X軸');
ylabel('Y軸');
zlabel('Z軸');
title('Figure2:3-D Plot')
實例20：圖形的隱藏屬性

function shili20
h0=figure('toolbar','none',...
    'position',[200 150 450 300],...
    'name','實例20');
subplot(1,2,1)
[x,y,z]=sphere(10);
mesh(x,y,z)
axis off
title('Figure1:Opaque')
hidden on

subplot(1,2,2)
[x,y,z]=sphere(10);
mesh(x,y,z)
axis off
title('Figure2:Transparent')
hidden off
實例21PEAKS函數曲線

function shili21
h0=figure('toolbar','none',...
    'position',[200 100 450 450],...
    'name','實例21');
[x,y,z]=peaks(30);
subplot(2,1,1)
x=x(1,:);
y=y(:,1);
i=find(y>0.8&y<1.2);
j=find(x>-0.6&x<0.5);
z(i,j)=nan*z(i,j);
surfc(x,y,z)
xlabel('X軸');
ylabel('Y軸');
zlabel('Z軸');
title('Figure1:surfc函數形成的曲面')

subplot(2,1,2)
x=x(1,:);
y=y(:,1);
i=find(y>0.8&y<1.2);
j=find(x>-0.6&x<0.5);
z(i,j)=nan*z(i,j);
surfl(x,y,z)
xlabel('X軸');
ylabel('Y軸');
zlabel('Z軸');
title('Figure2:surfl函數形成的曲面') 
實例22：片狀圖

function shili22
h0=figure('toolbar','none',...
    'position',[200 150 550 350],...
    'name','實例22');
subplot(1,2,1)
x=rand(1,20);
y=rand(1,20);
z=peaks(x,y*pi);
t=delaunay(x,y);
trimesh(t,x,y,z)
hidden off
title('Figure1:Triangular Surface Plot');

subplot(1,2,2)
x=rand(1,20);
y=rand(1,20);
z=peaks(x,y*pi);
t=delaunay(x,y);
trisurf(t,x,y,z)
title('Figure1:Triangular Surface Plot');
實例23：視角的調整

function shili23
h0=figure('toolbar','none',...
    'position',[200 150 450 350],...
    'name','實例23');
x=-5:0.5:5;
[x,y]=meshgrid(x);
r=sqrt(x.^2+y.^2)+eps;
z=sin(r)./r;