qamdraw.m

OFDM系统在QPSK与16QAM调制下

%16QAM调制下OFDM系统的误码率
function [ber1] = qamdraw;
para = 52;%并行信道数
fftlen = 64;%FFT长度
noc = 52;%载波数
nd = 3;%每个仿真循环的OFDM数
m1 = 4;%调制水平：QPSK
gilen = 16;%保护间隔长度
%ebn0 = 2;%信噪比
sr = 250000;%OFDM symbol rate (250k syombol/s)
br = sr.*m1;%bit rate per carrier
%**************************主循环部分*********************************
nloop = 200;%仿真的循环数
%noe1 = 0;%信道解码前错误数据数
%nod1 = 0;%信道解码前传输数据数
%*************************发射机**************************************
ebn0 = 0:15;;%信噪比
for ii = 1:length(ebn0);
    noe1 = 0;%信道解码前错误数据数
    nod1 = 0;%信道解码前传输数据数
for iii = 1:nloop
    %****************************数据产生******************************
    xx = ones(1,para*nd*m1);
    %**************************串并转换*********************************
    paradata = reshape(xx,para,nd*m1);
    %************************16qam调制**********************************
    [ich,qch] = qammod2(paradata,para,nd,m1);
    %数据映射
    [ich1,qch1] = cramapping(ich,qch,fftlen,nd);
    %*******************************IFFT********************************
    x = ich1 + qch1.*i;
    y = ifft(x);
    ich2 = real(y);
    qch2 = imag(y);
    %**************************插入保护间隔******************************
    [ich3,qch3] = giins1(ich2,qch2,fftlen,gilen,nd);
    fftlen2 = fftlen + gilen;
    %**************************attenuation calculation*********************
    spow = sum(ich3.^2+qch3.^2)/nd./para;
    attn = 0.5*spow*sr/br*10.^(-ebn0(ii)/10);
    attn = sqrt(attn);
    %**************************AWGN addition******************************
    [ich4,qch4] = comb(ich3,qch3,attn);
    %**********************去除保护间隔***********************************
    [ich5,qch5] = girem1(ich4,qch4,fftlen2,gilen,nd);
    %*************************FFT*****************************************
    rx = ich5 + qch5.*i;
    ry = fft(rx);
    ich6 = real(ry);
    qch6 = imag(ry);
    %载波逆映射
    [ich7,qch7] = crdemapping(ich6,qch6,fftlen,nd);
    %******************************解调**********************************
    [demodata] = qamdemod2(ich7,qch7,para,nd,m1);
    %*****************************并串转换*******************************
    demodata1 = reshape(demodata,1,para*nd*m1);
    %****************************bit error rate(BER)*********************
    noe10 = sum(abs(demodata1 - xx));
    nod10 = length(xx);
    
    noe1 = noe10 + noe1;
    nod1 = nod10 + nod1;
    
    %fprintf('%d\t%e\t\n',iii,noe10/nod10);
end

ber1(ii) = noe1/nod1;
%********************************输出结果**********************************
%fprintf('%f\t%e\t%d\t\n',ebn0(ii),ber1(ii),nloop);
end
%semilogy(ebn0,ber1);
%grid on;
%*******************************end of file********************************