scfde.m

单载波频域均衡SCFDE的MATLAB源码

% FDE 
function SER=scfde(originalsymbol,sb,SNR)

cprsb=sb(65:length(sb));

load hf
hf=ans;
hf=hf(1:206);
% hf=[1 0 10^(1/20) 0 10^(-9/20) 10^(-10/20) 0 0 0 10^(-15/20) 0 0 0 10^(-20/20)];
pulse=pulseshaping(50e9,2e-9,1);
hf=equichannel(hf,pulse,100);
E=sum(abs(hf).^2);
hf=hf./E;
H_channel=fft(hf,256);

Y = fft(cprsb, 256);

equalizerType='ZERO';

if equalizerType == 'ZERO'
    Y = Y./H_channel;
elseif equalizerType == 'MMSE'
    C = conj(H_channel)./(conj(H_channel).*H_channel + 10^(-SP.SNR/10));
    Y = Y.*C;
end

EstSymbols = ifft(Y);

EstSymbols = sign(EstSymbols);

I = find((originalsymbol(65:length(sb))-EstSymbols) == 0);
errCount = length(sb)-length(I)-64;

SER= errCount / (length(sb)-64);