基于MATLAB的数据处理与统计作图.doc

Smooth函数:
load ;
c=smooth(count(:));
C1=reshape(c,24,3);
subplot(3,1,1);plot(count,':');
hold on;
plot(C1,'-');
C2=zeros(24,3);
for I=1:3
C2(:,I)=smooth(count(:,I));
end;
subplot(3,1,2);plot(count,':');hold on;
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');hold on;plot(out2,'m.');
subplot(2,2,3);plot(out3,'k');hold on;plot(out4,'m.');
subplot(2,2,4);plot(out5,'k');hold on;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-.');
hold on;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