Matlab编程100例

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1. 实例1：三角函数曲线（1）
   
2. 
   
3. function shili01
   
4. h0=figure('toolbar','none',...
   
5.     'position',[198 56 350 300],...
   
6.     'name','实例01');
   
7. h1=axes('parent',h0,...
   
8.    'visible','off');
   
9. x=-pi:0.05:pi;
   
10. y=sin(x);
    
11. plot(x,y);
    
12. xlabel('自变量X');
    
13. ylabel('函数值Y');
    
14. title('SIN( )函数曲线');
    
15. grid on

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：MATLAB编程 MATLAB atlab matla Mat function position visible figure parent

1. 实例2：三角函数曲线（2）
   
2. 
   
3. function shili02
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 350],...
   
6.     'name','实例02');
   
7. x=-pi:0.05:pi;
   
8. y=sin(x)+cos(x);
   
9. plot(x,y,'-*r','linewidth',1);
   
10. grid on
    
11. xlabel('自变量X');
    
12. ylabel('函数值Y');
    
13. title('三角函数');

复制代码

1. 实例3：图形的叠加
   
2. 
   
3. function shili03
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 350],...
   
6.     'name','实例03');
   
7. x=-pi:0.05:pi;
   
8. y1=sin(x);
   
9. y2=cos(x);
   
10. plot(x,y1,...
    
11.     '-*r',...
    
12.     x,y2,...
    
13.     '--og');
    
14. grid on
    
15. xlabel('自变量X');
    
16. ylabel('函数值Y');
    
17. title('三角函数');

复制代码

1. 实例4：双y轴图形的绘制
   
2. 
   
3. function shili04
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 250],...
   
6.     'name','实例04');
   
7. x=0:900;a=1000;b=0.005;
   
8. y1=2*x;
   
9. y2=cos(b*x);
   
10. [haxes,hline1,hline2]=plotyy(x,y1,x,y2,'semilogy','plot');
    
11. axes(haxes(1))
    
12. ylabel('semilog plot');
    
13. axes(haxes(2))
    
14. ylabel('linear plot');

复制代码

1. 实例5：单个轴窗口显示多个图形
   
2. 
   
3. function shili05
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 250],...
   
6.     'name','实例05');
   
7. t=0:pi/10:2*pi;
   
8. [x,y]=meshgrid(t);
   
9. subplot(2,2,1)
   
10. plot(sin(t),cos(t))
    
11. axis equal
    
12. 
    
13. subplot(2,2,2)
    
14. z=sin(x)-cos(y);
    
15. plot(t,z)
    
16. axis([0 2*pi -2 2])
    
17. 
    
18. subplot(2,2,3)
    
19. h=sin(x)+cos(y);
    
20. plot(t,h)
    
21. axis([0 2*pi -2 2])
    
22. 
    
23. subplot(2,2,4)
    
24. g=(sin(x).^2)-(cos(y).^2);
    
25. plot(t,g)
    
26. axis([0 2*pi -1 1])

复制代码

1. 实例6：图形标注
   
2. 
   
3. function shili06
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 400],...
   
6.     'name','实例06');
   
7. t=0:pi/10:2*pi;
   
8. h=plot(t,sin(t));
   
9. xlabel('t=0到2\pi','fontsize',16);
   
10. ylabel('sin(t)','fontsize',16);
    
11. title('\it{从 0to2\pi 的正弦曲线}','fontsize',16)
    
12. x=get(h,'xdata');
    
13. y=get(h,'ydata');
    
14. imin=find(min(y)==y);
    
15. imax=find(max(y)==y);
    
16. text(x(imin),y(imin),...
    
17.     ['\leftarrow最小值=',num2str(y(imin))],...
    
18.     'fontsize',16)
    
19. text(x(imax),y(imax),...
    
20.     ['\leftarrow最大值=',num2str(y(imax))],...
    
21.     'fontsize',16)

复制代码

1. 实例7：条形图形
   
2. 
   
3. function shili07
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 350],...
   
6.     'name','实例07');
   
7. tiao1=[562 548 224 545 41 445 745 512];
   
8. tiao2=[47 48 57 58 54 52 65 48];
   
9. t=0:7;
   
10. bar(t,tiao1)
    
11. xlabel('X轴');
    
12. ylabel('TIAO1值');
    
13. h1=gca;
    
14. h2=axes('position',get(h1,'position'));
    
15. plot(t,tiao2,'linewidth',3)
    
16. set(h2,'yaxislocation','right','color','none','xticklabel',[])

复制代码

1. 实例8：区域图形
   
2. 
   
3. function shili08
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 250],...
   
6.     'name','实例08');
   
7. x=91:95;
   
8. profits1=[88 75 84 93 77];
   
9. profits2=[51 64 54 56 68];
   
10. profits3=[42 54 34 25 24];
    
11. profits4=[26 38 18 15 4];
    
12. area(x,profits1,'facecolor',[0.5 0.9 0.6],...
    
13.     'edgecolor','b',...
    
14.     'linewidth',3)
    
15. hold on
    
16. area(x,profits2,'facecolor',[0.9 0.85 0.7],...
    
17.     'edgecolor','y',...
    
18.     'linewidth',3)
    
19. hold on
    
20. area(x,profits3,'facecolor',[0.3 0.6 0.7],...
    
21.     'edgecolor','r',...
    
22.     'linewidth',3)
    
23. hold on
    
24. area(x,profits4,'facecolor',[0.6 0.5 0.9],...
    
25.     'edgecolor','m',...
    
26.     'linewidth',3)
    
27. hold off
    
28. set(gca,'xtick',[91:95])
    
29. set(gca,'layer','top')
    
30. gtext('\leftarrow第一季度销量')
    
31. gtext('\leftarrow第二季度销量')
    
32. gtext('\leftarrow第三季度销量')
    
33. gtext('\leftarrow第四季度销量')
    
34. xlabel('年','fontsize',16);
    
35. ylabel('销售量','fontsize',16);

复制代码

1. 实例9：饼图的绘制
   
2. 
   
3. function shili09
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 250],...
   
6.     'name','实例09');
   
7. t=[54 21 35;
   
8.     68 54 35;
   
9.     45 25 12;
   
10.     48 68 45;
    
11.     68 54 69];
    
12. x=sum(t);
    
13. h=pie(x);
    
14. textobjs=findobj(h,'type','text');
    
15. str1=get(textobjs,{'string'});
    
16. val1=get(textobjs,{'extent'});
    
17. oldext=cat(1,val1{:});
    
18. names={'商品一：';'商品二：';'商品三：'};
    
19. str2=strcat(names,str1);
    
20. set(textobjs,{'string'},str2)
    
21. val2=get(textobjs,{'extent'});
    
22. newext=cat(1,val2{:});
    
23. offset=sign(oldext(:,1)).*(newext(:,3)-oldext(:,3))/2;
    
24. pos=get(textobjs,{'position'});
    
25. textpos=cat(1,pos{:});
    
26. textpos(:,1)=textpos(:,1)+offset;
    
27. set(textobjs,{'position'},num2cell(textpos,[3,2]))

复制代码

1. 实例10：阶梯图
   
2. 
   
3. function shili10
   
4. h0=figure('toolbar','none',...
   
5.     'position',[200 150 450 400],...
   
6.     'name','实例10');
   
7. a=0.01;
   
8. b=0.5;
   
9. t=0:10;
   
10. f=exp(-a*t).*sin(b*t);
    
11. stairs(t,f)
    
12. hold on
    
13. plot(t,f,':*')
    
14. hold off
    
15. glabel='函数e^{-(\alpha*t)}sin\beta*t的阶梯图';
    
16. gtext(glabel,'fontsize',16)
    
17. xlabel('t=0:10','fontsize',16)
    
18. axis([0 10 -1.2 1.2])

复制代码