smldpe.m

通原实验课的作业

function [p]=smldPe(snr_in_dB)
% [p]=smldPe(snr_in_dB)
%	SMLDPE仿真求出以dB为单位的特定信噪比的误码率
E=1;
SNR=exp(snr_in_dB*log(10)/10); %信噪比
sgma=E/sqrt(2*SNR); %噪声的标准偏差
if (snr_in_dB<=6),
   N=5000;
elseif (snr_in_dB<=7)
   N=30000;
elseif (snr_in_dB<=8)
   N=60000;
end;
%产生信源
for i=1:N,
	temp=rand;	%在（0，1）区间内的均匀随机变量
	if (temp<0.5),
		dsource(i)=0;  %信源输出为0的概率为1/2
	else
		dsource(i)=1;  %信源输出为1的概率为1/2
	end
end
%随后为判决和误码概率的计算
numoferr=0;
for i=1:N,
	%匹配滤波器的输出
	if (dsource(i)==0),
		r=-E+sgma*randn; %如果信源输出为“0”
	else
		r=E+sgma*randn;  %如果信源输出为"1"
	end;
	%随后为判决
	if (r<0),
		decis=0;		 %判决为“0”
	else 
		decis=1;	  	%判决为“1”
	end;
	if(decis~=dsource(i)),
		numoferr=numoferr+1; %如果出错，计数器加1
	end;
end;
p=numoferr/N;	  %误差概率估算