无码间干扰基带系统误码率程序程序求助
 

程序是仿真不同信噪比下的误码率情况，其中sigma=sqrt(N0/2);这句是有什么作用呢？下面并没有用到sigma的值啊。。向大家请教。
EbN0dB=0:0.5:10;
N0=10.^(-EbN0dB/10);
sigma=sqrt(N0/2);
%理论计算的误码率
Pb=0.5*erfc(sqrt(1./N0));
%仿真误码率
for n=1:length(EbN0dB)
%产生等概信源
a=sign(randn(1,100000));
%离散等效接收模型
rk=a+sigma(n)*randn(1,100000);
dec_a=sign(rk); %判决
%计算误码率
ber(n)=sum(abs(a-dec_a)/2)/length(a);
end
semilogy(EbN0dB,Pb);
hold;
semilogy(EbN0dB,ber,'rd-');
legend('理论值','仿真结果');
xlabel('Eb/N0(dB)');
ylabel('Pb');