www.gusucode.com > 瑞利信道下噪声能量的估计源码程序 > 瑞利信道下噪声能量的估计源码程序/MSP_estimate_version1/Histogram/CalculateHistogram.m

    % Histogram calculation procedure with specified granularity
% 
% Tomasz Pieciak
% AGH university of Science and Technology, Krakow, Poland
% pieciak@agh.edu.pl, http://home.agh.edu.pl/pieciak/
%
% ARGUMENTS
%   data - input data
%   granularity - histogram granularity
%
% FUNCTION RETURNS
%   histogram_x - histogram abscissa
%   histogram_p - histogram frequencies
%
% USAGE
%   [histogram_x, histogram_p] = CalculateHistogram(data, 5);

function [histogram_x, histogram_data] = CalculateHistogram(data, granularity)

data = round(data ./ granularity) .* granularity;
tabulate_results = tabulate(data);

histogram_x = tabulate_results(:, 1);
histogram_data = tabulate_results(:, 2);
histogram_data = (histogram_data ./ sum(histogram_data)) ./ granularity;