基于matlab的数据处理与统计作图汇总.doc

Smooth函数:;c=smooth(count(:));C1=reshape(c,24,3);subplot(3,1,1);plot(count,':');holdon;plot(C1,'-');C2=zeros(24,3);forI=1:3C2(:,I)=smooth(count(:,I));end;subplot(3,1,2);plot(count,':');holdon;plot(C2,'-');subplot(3,1,3);plot(C2-C1,'o-');>>x=15*rand(150,1);y=sin(x)+*(rand(size(x))-);y(ceil(length(x)*rand(2,1)))=3;noise=normrnd(0,15,150,1);y=y+noise;>>yy1=smooth(x,y,,'loess');>>yy2=smooth(x,y,,'rloess');>>yy3=smooth(x,y,,'moving');>>yy4=smooth(x,y,,'lowess');>>yy5=smooth(x,y,,'sgolay');>>yy6=smooth(x,y,,'rlowess');>>[xx,ind]=sort(x);subplot(3,2,1);plot(xx,y(ind),'b-.',xx,yy1(ind),'r-');subplot(3,2,2);plot(xx,y(ind),'b-.',xx,yy2(ind),'r-');subplot(3,2,3);plot(xx,y(ind),'b-.',xx,yy3(ind),'r-');subplot(3,2,4);plot(xx,y(ind),'b-.',xx,yy4(ind),'r-');subplot(3,2,5);plot(xx,y(ind),'b-.',xx,yy5(ind),'r-');subplot(3,2,6);plot(xx,y(ind),'b-.',xx,yy6(ind),'r-');Smoothts函数:>>x=122+rand(500,4);p=x(:,4)';out1=smoothts(p,'b',30);out2=smoothts(p,'b',100);out3=smoothts(p,'g',30);out4=smoothts(p,'g',100,100);out5=smoothts(p,'e',30);out6=smoothts(p,'e',100);subplot(2,2,1);plot(p);subplot(2,2,2);plot(out1,'k');holdon;plot(out2,'m.');subplot(2,2,3);plot(out3,'k');holdon;plot(out4,'m.');subplot(2,2,4);plot(out5,'k');holdon;plot(out6,'m.');Medfilt1函数:x=linspace(0,2*pi,250)';y=sin(x)*150;y(ceil(length(x)*rand(2,1)))=3;noise=normrnd(0,15,250,1);y=y+noise;subplot(1,2,1);plot(x,y);yy=medfilt1(y,50);subplot(1,2,2);plot(x,y,'r-.');holdon;plot(x,yy,'.');直方图:Hist函数:x=randn(499,1);y=randn(499,3);subplot(3,1,1);hist(x);subplot(3,1,2);hist(x,100);subplot(3,1,3);hist(y,25);Histc函数:x=-::;y=randn(10000,1);subplot(1,2,1);hist(y,x);n=histc(y,x);c=cumsum(n);subplot(1,2,2);bar(x,c);Histfit函数:>>r=normrnd(10,1,10,1);>>histfit(r);>>h=get(gca,'Children');盒子图:>>N=1024;x1=normrnd(5,1,N,1);>>x2=normrnd(6,1,N,1);x=[x1x2];subplot(2,2,1);sym1='*';notch1=1;boxplot(x,notch1,sym1);>>subplot(2,2,2);notch2=0;boxplot(x,notch2);x=randn(100,25);subplot(2,1,1);boxplot(x);误差条图:>>x=0: