香农码字matlab,香农编码的MATLAB实现.ppt

* * 实验二 香农编码的MATLAB实现 一、实验目的 熟悉MATLAB平台及工作环境 熟练使用条件判断语句、循环语句 掌握香农编码的步骤 掌握利用MATLAB实现香农编码 二、实验原理 1. 香农编码步骤 (1)将信源符号按概率从大到小的顺序排序： (2)令 p(a0) = 0，用 pa (aj)(j = i+1)表示第 i 个码字的累加概率，则： 二、实验原理 1. 香农编码步骤 (3)令 ki 为第 i 个码字的码长，按照下式计算出各概率对应的码字长度 ki ： (4)把各个累加概率pa (aj)由十进制小数转化为二进制小数，并取小数点后的前 ki 位作为对应信源符号ai的码字。 二、实验原理 【例 3.3.1】 设有一单符号离散无记忆信源： 对该信源编二进制香农码。 2. 以教材P39例题3.3.1信源为例 二、实验原理 i ai p(ai) j pa(aj) ki 码 字 a0 0 1 a1 0.25 1 0.000 2 00 (0.000)2 2 a2 0.25 2 0.250 2 01 (0.010)2 3 a3 0.20 3 0.500 3 100 (0.1000)2 4 a4 0.15 4 0.700 3 101 (0.1011)2 5 a5 0.10 5 0.850 4 1101 (0.11011)2 6 a6 0.05 6 0.950 5 11110 (0.111100)2 2. 以教材P39例题3.3.1信源为例 二、实验原理 平均码长： 信源熵： 信息率： 编码效率： 2. 以教材P39例题3.3.1信源为例 三、实验内容 1. 程序设计 (1)代码初始化 clc; % 清除命令窗口的内容 clear all; % 清除工作空间的内容 p = [0.25 0.25 0.20 0.15 0.10 0.05]; % 输入信源分布 n = length(p); % 计算信源符号的个数 三、实验内容 (2)判断信源是否合法，计算信源熵 if min(p)<0 || max(p)>1 % 判断概率值是否介于0到1之间 disp(['信源分布 P(X)=[',num2str(p),']']); disp('概率值必须介于0到1之间，请重新输入信源分布'); return elseif sum(p) ~= 1 % 判断概率累加和是否等于1 disp(['信源分布 P(X)=[',num2str(p),']']); disp('概率累加和必须等于1，请重新输入信源分布'); return 三、实验内容 (2)判断信源是否合法，计算信源熵 else H = 0; % 熵值初始化为0 for i = 1:n H = -p(i)*log2(p(i)) + H; % 累加各个符号的概率加权 end end 三、实验内容 (3)按照概率大小排序，并计算累加概率 p = sort(p, 'descend'); % 降序排列 pa = zeros(1,n); % 对累加概率矩阵初始化 for j = 2:n pa(j) = pa(j-1) + p(j-1); % 计算累加概率 end 三、实验内容 (4)计算每个码字的长度 k = zeros(1,n); % 对码长矩阵初始化 for i = 1:n k(i) = ceil(-log2(p(i))); % 计算每个码字长度，向上取整 end 三、实验内容 (5)十进制小数转换成二进制小数 k_max = max(k); bin = zeros(n,k_max); % 对二进制小数矩阵初始化 for i = 1:n for j = 1:k_max bin(i,j)=floor(pa(i)*2); % 将累加概率转换成二进制 pa(i)=pa(i)*2-floor(pa(i)*2); end end 三、实验内容 (6)取小数点后 ki 位作为码字 codeword = zeros(n,k_max); % 对码字矩阵初始化 for i = 1:n for j = 1:k_max if j > k(i) codeword(i,j) = blanks(1); % 创建空白串 else