求高斯信道的信道容量 信息论

clear,clc;

%一维二元高斯信道的信道容量

xn=8000;%定义点数

m1=0;

m2=1;%信号参数

sigma1=0.38;

sigma2=0.5;%高斯噪声参数

xmin=min([m1-5*sigma1,m2-5*sigma2]);xmax=max([m1+5*sigma1,m2+5*sigma2]);%定义域边界

x=linspace(xmin,xmax,xn+1);

y1=(((2*pi)^(-0.5))/sigma1)*exp(-(x-m1).^2/(2*(sigma1^2)));

y2=(((2*pi)^(-0.5))/sigma2)*exp(-(x-m2).^2/(2*(sigma2^2)));%信号受干扰后的概率密度分布函数

a=(sigma2^2-sigma1^2)/(2*sigma1^2*sigma2^2);

b=(2*m2*sigma1^2-2*m1*sigma2^2)/(2*sigma1^2*sigma2^2);

c=(sigma2^2*m1^2-sigma1^2*m2^2)/(2*sigma1^2*sigma2^2)-log(sigma2/sigma1);%联立y1,y2解二次方程求阈值

if a==0

    x4=c/b;%二