0001 function [ norm_fact ] = ComputeNormFact( Para )
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021 proto_filt = protofilter_design(Para);
0022 modulated_filters = modulate_proto_filt(proto_filt, Para);
0023 NsubCar = Para.SubBandWidth*Para.NbSubBands;
0024 iSubBand=1;
0025 G_i= toeplitz([modulated_filters(iSubBand,:) zeros(1,Para.nSubcarriers-1)],[modulated_filters(iSubBand,1) zeros(1,Para.nSubcarriers-1)]) ;
0026 B=1/sqrt(Para.nSubcarriers).*fft(G_i');
0027 index=(iSubBand-1)*Para.SubBandWidth+(Para.nSubcarriers-NsubCar)/2+1:(iSubBand-1)*Para.SubBandWidth+(Para.nSubcarriers-NsubCar)/2+Para.SubBandWidth;
0028 norm_fact=sum(diag(B(index,:)*B(index,:)'));
0029
0030 norm_fact=sqrt(Para.SubBandWidth*Para.Es/norm_fact);
0031
0032 end
0033