武汉理工大学信息理论编码实验.docVIP

实验1 1.p=0.00001:0.001:0.99999; h=-p.*log2(p)-(1-p).*log2(1-p); plot(p,h); title(二进熵函数曲线); ylabel(H(p,1-p)) 2. p=linspace(eps,1-eps,100); q=linspace(eps,1-eps,100); [P,Q]=meshgrid(p,q); P_Q=P+Q; for n=1:100 for m=1:100 if P_Q(n,m)=1 Q(n,m)=nan; end end end H=-P.*log2(P)-Q.*log2(Q)-(1-P-Q).*log2(1-P-Q); mesh(P,Q,H) title(三维熵函数的图形) 3.function H=entropy(p) %该函数用来计算包含任意多个符号的信源熵 %p为DMS的概率分布，行向量 if sum(p)~=1 %判断是否满足概率完备性 error(!!!不满足概率完备性，重新输入信源分布) return else L=length(p) %得到信源符号的个数 H=0; %熵值初始化为零 for i=1:L H=H-p(i)*log2(p(i)); %累加熵函数各个子项 end end 3.[p,q]=meshg