cm_sm32.m

本程序实现对多进制PSK数字调制系统的仿真。程序运行需要几分钟

function[pb,ps]=cm_sm32(snr_in_dB)
%[pb,ps]=cm_sm32(snr_in_dB)
%        CM_SM32求出QPSK以dB为单位的给定信噪比的
%        比特误码率和符号误码率
N=10000;
E=1;                      %每符号能量
snr=10^(snr_in_dB/10);    %信噪比
sgma=sqrt(E/snr)/2;       %噪声标准方差
%信号映射
s00=[1  0];
s01=[0  1];
s11=[-1 0];
s10=[0 -1];
%数据源的产生
for i=1:N,
	temp=rand;         %在区间（0，1 )间的一个均匀随机变量
	if  (temp<0.25),        %信源输出为“00”的概率为1/4
	    dsource1(i)=0;    
	    dsource2(i)=0;
	elseif (temp<0.5),      %信源输出为“01”的概率为1/4
	    dsource1(i)=0;
	    dsource2(i)=1;
	elseif (temp<0.75),     %信源输出为“10”的概率为1/4
	    dsource1(i)=1;
	    dsource2(i)=0;
	else                    %信源输出为“11”的概率为1/4
    	dsource1(i)=1; 
	    dsource2(i)=1;
	end;
end;
%判决、误码率的计算
numofsymbolerror=0;
numofbiterror=0;
for i=1:N,
    %在判决器处的接收信号，对于第i个符号为：
	n(1)=gngauss(sgma);
	n(2)=gngauss(sgma);
	if ((dsource1(i)==0) & (dsource2(i)==0)),
		r=s00+n;
	elseif ((dsource1(i)==0) & (dsource2(i)==1)),
		r=s01+n;
	elseif ((dsource1(i)==1) & (dsource2(i)==0)),
		r=s10+n;
	else 
		r=s11+n;
	end;
    %以下为计算互相关量度
	c00=dot(r , s00);
	c01=dot(r , s01);	
	c10=dot(r , s10);
	c11=dot(r , s11);
    %第i个符号的判决按如下方式进行
	c_max=max([c00 c01 c10 c11]);
	if (c00==c_max)
		decis1=0;decis2=0;
	elseif (c01==c_max)
		decis1=0;decis2=1;
	elseif (c10==c_max)
		decis1=1;decis2=0;
	else
		decis1=1;decis2=1;
	end;
	%如果判决不对，计错器加1
	symbolerror=0;
	if (decis1~=dsource1(i)),
		numofbiterror=numofbiterror+1;
		symbolerror=1;
	end;
	if (decis2~=dsource2(i)),
		numofbiterror=numofbiterror+1;
		symbolerror=1;
	end;
	if (symbolerror==1),
		numofsymbolerror=numofsymbolerror+1;
	end;
end;
ps=numofsymbolerror/N;     %因为总共有N个符号
pb=numofbiterror/(2*N);    %因为发送了2N个比特