hist_rate.m

来自「概率统计实验源程序！这是一个学习使用概率中心极限定理的工具箱」· M 代码 · 共 29 行

function [y,n]=hist_rate(x,tag)

%对一列数x1,x1,...,xn从小到大进行排序，并去除多余的重复值，返回处理后向量y，
%当tag==1时，返回累积频率向量n,对应y中每个数字出现的累积频率
%当tag==2时，返回频率向量n,对应y中每个数字出现的频率
%当tag==3时，返回频数向量n,对应y中每个数字出现的频数

if ~isempty(x)
    x=sort(x);
    m=length(x);
    x1=x(2:end)-x(1:end-1);
    x1(end+1)=1;
    x1=find(x1~=0);
    if tag==1
        n=x1/m;
    end
    y=x(x1);
    x1(2:end+1)=x1;
    x1(1)=0;
    if tag==2
        n=(x1(2:end)-x1(1:end-1))/m;
    end
    if tag==3
        n=x1(2:end)-x1(1:end-1);
    end    
else
    warndlg('输入数据不能为空');
end
x=[];x1=[];m=[];