通信原理仿真实验系统设计与实现 第14页

通信原理仿真实验系统设计与实现 第14页
mapping=[-3*d 3*d; 
    -d 3*d; d 3*d; 3*d 3*d;  -3*d d;
    -d d; d d; 3*d d; -3*d -d;
    -d -d; d -d;3*d -d; -3*d -3*d;
    -d -3*d; d -3*d; 3*d -3*d];
for i=1:N
    qam_sig(i,: )=mapping(dsource(i),: );
end;
for i=1:N,
    n=gngauss(sgma); %产生高斯随机噪声
    r(i,: )=qam_sig(i,: )+n;%在信号上參加噪声
end;
numoferr=0; %误比特数初始值置为０
for i=1:N,
    for j=1:M,
        metrics(j)=(r(i,1)-mapping(j,1))^2+(r(i,2)-mapping(j,2))^2;
    end;
    [min_metric decis]=min(metrics);
    if(decis~=dsource(i)), %若出现错的情况，误比特数加一
        numoferr=numoferr+1;
    end;
end;
p=numoferr/(N);

%qammotocalo.m
echo on
SNRindB1=0:2:15;   SNRindB2=0:0.1:15;
M=16;  k=log2(M);
for i=1:length(SNRindB1),
    smld_err_prb(i)=qammoto(SNRindB1(i));
end;
for i=1:length(SNRindB2),
    SNR=exp(SNRindB2(i)*log(10)/10);
    theo_err_prb(i)=4*Qfunct(sqrt(3*k*SNR/(M-1)));
end;
semilogy(SNRindB1,smld_err_prb,'*'); %用对数坐标做出实际信噪比
hold
semilogy(SNRindB2,theo_err_prb); %画出理论信噪比
%函数文件gugauss.m产生高斯分布的随机函数
function grsv=gngauss(m,sgma)

if(nargin==0),
    m=0;   sgma=1;
elseif nargin==1
    sgma=m; m=0;
end;
u=rand;
z=sgma*(sqrt(2*log(1/(1-u)))); %利用上面的Ｕ产生一个瑞利公布随机数
u=rand;
grsv(1)=m+z*cos(2*pi*u);   grsv(2)=m+z*sin(2*pi*u);
%函数文件Qfunct.m用于求出理论误码率
function[y]=Qfunct(x)
y=(1/2)*erfc(x/sqrt(2));

 << 上一页  [11] [12] [13] [14] 

通信原理仿真实验系统设计与实现 第14页下载如图片无法显示或论文不完整，请联系qq752018766