matlab实验——求信源可靠传输的ε范围

实验要求：
概率为P(X=0)=1/4的二进制信源，在交叉概率为ε的二进制对称信道中传输，确定使信源可靠传输的ε的范围。假定信道对每个输入信源能够使用一次。附代码和结果图。
代码如下：

Se=-(1/4.*log2(1/4)+3/4.*log2(3/4));  %求信源熵
p=0:0.005:1;   %在0-1区间上以0.005为步长
Hb=-(p.*log2(p)+(1-p).*log2(1-p));   %p为传输错误率
C=1-Hb;
k=find(abs(C-Se)<=0.01);   %返回满足条件的值
plot(p,Se,'b+',p,C,'r',p(k),Se,'ro','MarkerSize',5);  %绘图

结果图如下：

由上面的运行结果图可以看出，在两个交点之间的区域就是信源可靠传输的范围，故信源可靠传输的ε的范围是（0.03，0.97）。