FBMC 的 AWGN 信道仿真（matlab）

FBMC 的 AWGN 信道仿真(matlab实现) The simulation of FBMC with AWGN channel.

%% FBMC/OQAM
%% FBMC-AWGN: BER verus SNR for FBMC/OQAM using polyphase structure  

%%------------------------------------------------------------------------------------------------------------%%
%
%             “Prototype filter and filter bank structure“
% FBMC/OQAM with noise
%%------------------------------------------------------------------------------------------------------------%%
clc
clear all;
%%------------------------------------------------------------------------------------------------------------%%
%                                               Intialization
%%------------------------------------------------------------------------------------------------------------%%
nd=2; %number of OQAM sample per symbols （not change）
M=512;% number of sub-channels
K=4; % overlapping factor（good K=4） 
lp=K*M-1;% prototype filter length
SNR= 0:5; % signal to noise ratio
nloop=10;  % Number of simulation loops
Rate0= zeros（1 length（SNR））;
%%------------------------------------------------------------------------------------------------------------%%
%                                           Design of Prototype Filter 
%%------------------------------------------------------------------------------------------------------------%%
%  low pass filter FIR
y=[1 0.97195983 sqrt（2）/2 0.23514695 ]; % coefficient of filter when K=4
yy=[1 0.911438 0.411438];% coefficient of filter when K=3
yyy=[1 sqrt（2）/2];% coefficient of filter when K=2
u=lp-1;
s=2*pi/（K*M）;
for m=0:u-1
    r=m+1;
p（m+1）=y（11）-2*y（12）*cos（r*s）+2*y（13）*cos（2*r*s）-2*y（14）*cos（3*r*s）;% prototype filter equation K=4
%p（m+1）=yy（11）-2*yy（12）*cos（r*s）+2*yy（13）*cos（2*r*s）; % prototype filter equation K=3
%p（m+1）=yyy（11）-2*yyy（12）*cos（r*s）; % prototype filter equation K=2
%p（m+1）=1; % prototype filter equation for OFDM
end
p=[0 p 0]; % extra delay sample （z^-2）
ty=p; % 
%%------------------------------------------------------------------------------------------------------------%
%%------------------------------------------------------------------------------------------------------------%%
%%------------------------------------------------------------------------------------------------------------%%
%                                               Transmitter 
%%------------------------------------------------------------------------------------------------------------%%
for ebn0=1:length（SNR）
ebn0
for iii=1:nloop  
% Data generation
q=2;% binary data
m=4;% modulation level
In_Data=randi（[0 q-1]M1）;% random input data M x 1
% OQAM Modulation
[OQAM_In_Data]=OQAM_preprocessing（In_DatamM）; % OQAM modulation M x nd
%%------------------------------------------------------------------------------------------------------------%%
%                                               Transform Block  
%%------------------------------------------------------------------------------------------------------------%%
%% beta carrier 
for k

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        7575  2013-10-14 08:57  FMBC_OQAM_AWGN.m
     文件         510  2013-10-14 08:57  OQAM_postprocessing.m
     文件         490  2013-10-13 22:54  OQAM_preprocessing.m