0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014 Para=InitializeChainParameters( 'FBMC-OQAM' )
0015 Para.ActiveSubcarriers= 3:Para.nSubcarriers;
0016 Para.M_PAM=2;
0017 Para.Ns=2000;
0018 Para.N_T=2;
0019 Para.N_R=2;
0020 Para.S=2;
0021 Para.Es_N0_dB=15;
0022
0023
0024 d = GenerateData (Para);
0025
0026 C=GenerateRayleighChannelReal('ITU_VehA', Para);
0027
0028 ParaPreEqualizer.criterion='SVD';
0029 [ d_tilde, A, B ] = PreEqualizer( d, C, Para, 'SingleTap', ParaPreEqualizer );
0030 s = Modulator(d_tilde, Para );
0031
0032
0033 r = Channel_Multipath( s, C );
0034 r=Channel_AWGN( r, Para );
0035
0036
0037 z = Demodulator( r, Para );
0038
0039 ParaEqualizer.criterion='Specific';
0040 ParaEqualizer.B=B;
0041 [ x, B ] = Equalizer( z, C, Para, 'SingleTap', ParaEqualizer );
0042 d_hat=real(x);
0043 MSE = MSEComputes( d, d_hat, Para );
0044
0045
0046 figure
0047 plot(Para.ActiveSubcarriers, 10*log10(MSE),'-xb','markersize',8,'linewidth',1.5)
0048 hold on
0049 xlabel('Subcarrier index')
0050 ylabel('MSE [dB]')
0051 xlim([min(Para.ActiveSubcarriers) max(Para.ActiveSubcarriers)])