siso_ofdm.m

OFDM系统的仿真代码MATLAB环境下

function siso_ofdm

global sign;
global sim_consts;

set_sim_consts;

for sign = 1:2                       % Comparison between Kalman and LS
    ber(sign,:) = single_packet;
end
semilogy(sim_consts.SNR,ber(1,:),'r',sim_consts.SNR,ber(2,:),'-ro'),
xlabel('SNR (dB)'),ylabel('BER'),
legend('Kalman-FFT','LS-FFT');

% for sign = 1:3
%     if sign == 1
%         sim_consts.NumSubc = 52;
%         sim_consts.UsedSubcIdx = [7:58]';
%         sim_consts.NumDataSubc = 48;
%         sim_consts.NumPilotSubc = 4;
%         sim_consts.DataSubcIdx = [8:19 21:32 34:45 47:58]';
%         sim_consts.PilotSubcIdx = [7 20 33 46]';
%         sim_consts.PilotSubcPatt = [1 14 27 40]';
%         sim_consts.DataSubcPatt = [2:13 15:26 28:39 41:52]';
%         sim_consts.PilotSubcSymbols = [1;1;1;1];
%     elseif sign == 2         sim_consts.NumSubc = 54;
%         sim_consts.UsedSubcIdx = [6:59]';
%         sim_consts.NumDataSubc = 48;
%         sim_consts.NumPilotSubc = 6;
%         sim_consts.DataSubcIdx = [7:14 16:23 25:32 34:41 43:50 52:59]';
%         sim_consts.PilotSubcIdx = [6 15 24 33 42 51]';
%         sim_consts.PilotSubcPatt = [1 10 19 28 37 46]';
%         sim_consts.DataSubcPatt = [2:9 11:18 20:27 29:36 38:45 47:54]';
%         sim_consts.PilotSubcSymbols = [1;1;1;1;1;1];
%     elseif sign == 3
%         sim_consts.NumSubc = 56;
%         sim_consts.UsedSubcIdx = [5:60]';
%         sim_consts.NumDataSubc = 48;
%         sim_consts.NumPilotSubc = 8;
%         sim_consts.DataSubcIdx = [6:11 13:18 20:25 27:32 34:39 41:46 48:53 55:60]';
%         sim_consts.PilotSubcIdx = [5 12 19 26 33 40 47 54]';
%         sim_consts.PilotSubcPatt = [1 8 15 22 29 36 43 50]';
%         sim_consts.DataSubcPatt = [2:7 9:14 16:21 23:28 30:35 37:42 44:49 51:56]';
%         sim_consts.PilotSubcSymbols = [1;1;1;1;1;1;1;1];
%     end
%    ber(sign,:) = single_packet;
% end
% semilogy(sim_consts.SNR,ber(1,:),'r',sim_consts.SNR,ber(2,:),'-rd',sim_consts.SNR,ber(3,:),'-ro'),
% xlabel('SNR (dB)'),ylabel('BER'),
% legend('Pilot Number=4','Pilot Number=6','Pilot Number=8');