使用matlab实现s={'s1','s2','s3','s4','s5'}; p=[0.4,0.2,0.2,0.1,0.1];针对二进制编码，对离散信源进行霍夫曼编码，能够打印出编码的结果

在MATLAB中，你可以使用霍夫曼编码算法对离散信源进行编码，其基本步骤包括计算频率、构造霍夫曼树和生成编码规则。对于给定的字符串数组`s`和概率向量`p`，可以按照以下步骤操作：

1. **计算频率**：

   freq = p; % 字符串概率直接对应于频率

2. **构建霍夫曼树**：
首先，将频率作为叶子节点创建哈希表（结构体），然后通过合并频率最小的两个节点不断构建非叶节点，直到所有元素都被纳入。可以手动编写函数或者使用MATLAB的并查集工具。

3. **生成编码规则**：
- 对于每个字符，从根节点开始遍历霍夫曼树，遇到分叉就向左走代表0，向右走代表1，记录路径上的字符序列即为其编码。
- 可以使用递归的方式来实现这个过程。

function [codes, code_book] = huffman_encode(s, p)
    % ... (此处省略具体的霍夫曼树构建代码)

    % 初始化空编码映射和代码书
    codes = cell(length(s), 1);
    code_book = containers.Map;

    % 递归生成霍夫曼编码
    function encode(codes, current_node, parent_code)
        if ismember(current_node, leaves)
            % 如果当前节点是叶子，添加到代码书和对应的编码
            codes{current_node} = [parent_code char(current_node)];
            code_book(s{current_node}) = codes{current_node};
        else
            left_child = find(current_node == children(1));
            right_child = find(current_node == children(2));
            
            % 继续编码左右子节点
            codes{current_node} = [encode(codes, left_child, [parent_code '0'])];
            codes{current_node} = [encode(codes, right_child, [parent_code '1'])];
        end
    end

    leaves = find(freq > 0); % 获取频率非零的叶子节点
    roots = ones(size(leaves)); % 根节点初始化为1
    for i = 2:length(leaves) % 构建霍夫曼树
        [i1, i2] = min([freq(leaves).*roots]); % 找到频率最小的两个节点
        children = [leaves(i1) leaves(i2)]; % 子节点索引
        freq(children) = freq(children) + roots(i1); % 更新频率
        roots = roots - roots(i1); % 更新剩余节点的权重
        current_node = roots'*[1; 2]; % 新节点由频率较小的两个节点组成
        children = [children(2) children(1)]; % 按从小到大排序
        [codes, ~] = huffman_encode(codes, current_node, '');
    end

    % 输出编码结果
    codes
end

运行上述代码后，会返回一个字符串向量，表示每个字符对应的霍夫曼编码。例如：

s = {'s1', 's2', 's3', 's4', 's5'};
p = [0.4, 0.2, 0.2, 0.1, 0.1];
[codes, code_book] = huffman_encode(s, p);
disp(code_book); % 显示编码字典