ofdm.rar_OFDM 正交_ofdm上行_spectral efficiency_spectral ofdm_uplink

function SER = ofdm(SP) numSymbols = SP.FFTsize; H_channel = fft(SP.channel,SP.FFTsize); for n = 1:length(SP.SNR), tic; errCount = 0; for k = 1:SP.numRun, tmp = round(rand(2,numSymbols)); tmp = tmp*2 - 1; inputSymbols = (tmp(1,:) + i*tmp(2,:))/sqrt(2); TxSamples = sqrt(SP.FFTsize)*ifft(inputSymbols); ofdmSymbol = [TxSamples(numSymbols-SP.CPsize+1:numSymbols) TxSamples]; RxSamples = filter(SP.channel, 1, ofdmSymbol); % Multipath Channel tmp = randn(2, numSymbols+SP.CPsize); complexNoise = (tmp(1,:) + i*tmp(2,:))/sqrt(2); noisePower = 10^(-SP.SNR(n)/10); RxSamples = RxSamples + sqrt(noisePower)*complexNoise; EstSymbols = RxSamples(SP.CPsize+1:numSymbols+SP.CPsize); Y = fft(EstSymbols, SP.FFTsize); if SP.equalizerType == 'ZERO' Y = Y./H_channel; elseif SP.equalizerType == 'MMSE' C = conj(H_channel)./(conj(H_channel).*H_channel + 10^(-SP.SNR(n)/10)); Y = Y.*C; end EstSymbols = Y; EstSymbols = sign(real(EstSymbols)) + i*sign(imag(EstSymbols)); EstSymbols = EstSymbols/sqrt(2); I = find((inputSymbols-EstSymbols) == 0); errCount = errCount + (numSymbols-length(I)); end SER(n,:) = errCount / (numSymbols*SP.numRun); [SP.SNR(n) SER(n,:)] toc end