matlab实程序百例1.doc

频谱分析: fs=500; t=0:1/fs:1; x1=5+10*sin(2*pi*100*t)+sin(2*pi*50*t); Y1=fft(x1,512); Y2=abs(Y1); f=(0:256)*250/256; subplot(3,1,1); plot(f,Y2(1:257)) 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