main.m

OFDM的simulink仿真

clear
clc
Es_No = 0:15;    % 信噪比  dB
BER = Es_No;   % 误比特率，长度与Es_No相同

%   循环执行仿真程序
for ii = 1:length(Es_No)
    EsNo = Es_No(ii);
    sim('OFDMModulate');
    sim('MC_CDMASim');
    load rxerror;
    rxerror(1,:) = 0;
    error=sum(sum(rxerror));
    BER(ii) = error/(1.2e5);
end

hold on
grid on
semilogy(Es_No,BER,'r+:');
save data_for_plot BER Es_No


% %%%计算MQAM的理论误码率曲线
% %%%MQAM的理论误码率公式——P=[2*（L-1）/L]*Q[sqrt（[6*log2(L)/(L^2-1)]*(Eb/n0)）]  
% M = 16;
% L = sqrt(M);
% for x = 1:16
%     y(x) = x-1;
%     EbNo = 10^(y(x)/10);
%     P(x) = 2*(L-1)/L*q(sqrt(6*log2(L)/(L^2-1)*EbNo));
% end
% semilogy(y,P,'k:');

%%%计算BPSK的理论误码率曲线
%%%P=Q（sqrt(2*Eb/no)）
for ii = 1:11
    Ebno_indB(ii) = ii-1;
    Ebno = 10^(Ebno_indB(ii)/10);
    Pe(ii) = Q(sqrt(2*Ebno));
end
semilogy(Ebno_indB,Pe,'k-');