huffman,m

%huffman源程序
function CODE=huffman(p)
%huffman为源符号建立一个可变长度的Huffman编码
%对于符号概率向量p在单元数组CODE中的二进制字符串返回一个Huffman编码
%在CODE中的每个元素对应一个字符，它的概率与指数p符合。

%检查输入参数
error(nargchk(1,1,nargin));
if (ndims(p)~=2) | (min(size(p))>1) | ~isreal(p) | ~isnumeric(p)
    error('P must be a real numeric vector.');
end

%全局变量存在于所有递归函数'makecode'
global CODE
CODE=cell(length(p),1);%初始的全局单位数组

if length(p)>1%当多余一个符号时
    p=p/sum(p);%标准化的输入概率
    s=reduce(p);%简约的源符号
    makecode(s,[]);%递归产生编码
else
    CODE={'1'};%平凡的一个符号例子
end