不同调制比较.m

不同调制方法,和不同编码方式的信道仿真比较

% 比较了不同调制情况时的误比特率
%  BPSK,QPSK,.....
%  QAM,16QAM.....
clear 
n_frame=10;
L_frame=10000;
M=4;    %1--BPSK    2--QPSK or QAM    3---8PSK   4---16PSK or 16QAM
          % psk        psk     qask        psk         psk     qask
Es=10;   % 数字映射后的符号能量 psk--1
        %                      16QAM--10
start=-16;
step=2;
finish=20;

Eb=Es/M; % 每比特能量
plot_pe=[];
Q=1;
axis_EbN0=start:step:finish;
%*********************************************************************
for EbN0=start:step:finish
    Liner_EbN0=10^(EbN0/10);
    pe_number=0;
    variance=0.5*(Eb/Liner_EbN0);  % 噪声方差   variance是单独加在同相（正交相）的噪声功率普密度
    for i=1:1:n_frame
        symbol=randint(L_frame/M,1,2^M); % 生成随机符号，化为0 1 二进制后比特个数应为L_frame
        x_msg=symbol2bit(symbol,M); % 化为0 1 二进制
        
        modulation_sig=modmap(symbol,1,1,'qask',2^M); % 数字映射
        
        x_ch=g_ch(modulation_sig,variance);  % 经过高斯信道，噪声方差var
        
        x_de_symbol=demodmap(x_ch,1,1,'qask',2^M); %逆映射
        
        x_de_msg=symbol2bit(x_de_symbol,M);
        
        pe_number=pe_number+sum(x_msg~=x_de_msg);
        
        current_time=fix(clock);
       fprintf('i am working %g, total(%g)     %g年  %g月  %g日  %g时  %g分  %g秒\n\n',Q,(((finish-start)/step)+1)*n_frame,current_time(1),current_time(2),current_time(3),current_time(4),current_time(5),current_time(6))
        Q=Q+1;
        fprintf('\n\n')
        
    end
    pe=pe_number/(L_frame*n_frame);
    plot_pe=[plot_pe,pe];
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
semilogy(axis_EbN0,plot_pe,'b*-')

xlabel('Eb/N0 信号功率/噪声功率(dB)')
ylabel('BER')
title('误码率')
grid on