matlab 画图 pdf plot

function pdfplot(xin, nbins)
% pdfplot(X, nbins) 
% displays a histogram of the empirical probability density function (PDF)
% for the data in the input array X using nbins number of bins 
% (by default pdfplot sets nbins to 20).
% If input X is a matrix, then pdfplot(X) parses it to the vector and 
% displays PDF of all values.
% For complex input X, pdfplot(X) displays PDF of abs(X).
%
% Example:
% y = randn( 1, 1e5 );
% pdfplot( y );
% pdfplot( y, 100 );

% Version 1.0
% Alex Bur-Guy, September 2003
% alex@wavion.co.il
%
% Revisions:
% Version 1.0 - initial version

if nargin == 1, nbins = 20; end
xin = reshape( xin, numel(xin), 1 );
if ~isreal( xin ), xin = abs( xin ); end
minXin = min(xin); maxXin = max(xin);
if floor( nbins ) ~= nbins, error( 'Number of bins should be integer value' ); end
if nbins < 2, error( 'Number of bins should be positive integer greater than 1 ' ); end
figure;
if minXin == maxXin
bar(minXin,1);
axis([minXin - 10, minXin + 10, 0, 1]);
else
step = (maxXin - minXin) / (nbins-1);
binc = minXin : step : maxXin; 
[N, X] = hist(xin, binc);
bar(X, N/sum(N));
end
xlabel('X', 'FontWeight','b','FontSize',12);
title(['PDF(X) based on ' num2str(length(xin)) ' data samples @ ' num2str(nbins) ' bins'], 'FontWeight','b','FontSize',12);
grid on;


reference: http://www.ilovematlab.cn/viewthread.php?tid=50959