信源熵实验程序:
clc; close all; clear; linwidd=1
fontt=20
p0=0; pd=1; N=20
p=linspace(p0,pd,N);
I=-log2(p);
plot(p,I,'k');
title('I=-log2(p)函数图'); xlabel('p');ylabel('I');
clc; close all; clear;
linwidd=1
fontt=20
p0=0; pd=1; N=20
p=linspace(p0,pd,N);
H=-p.*log2(p)-(1-p).*log2(1-p);
plot(p,H,'k');
title('H=-p.*log2(p)-(1-p).*log2(1-p)函数图'); xlabel('p');ylabel('H');
信道容量实验程序:
clc; close all; clear;
linwidd=1
fontt=20
p0=0; pd=1; N=20
p=linspace(p0,pd,N);
r=4
c=log2(r)+(1-p).*log2(1-p)+p.*log2(p/(r-1)); plot(p,c,'k');
title('强对称信道容量数值模拟图');