Matlab实现SC-FDMA 上传版本.zip

%* File Name: detectit.m %* Function Name: detectit() %* Description: one detection algorithm (hard decision at this version) %* %* Current Version: 0.1 %* Time: 27.3.2008 %* Programmer: Neo function detected = detectit(bits_modu,modu) % variables Eucli_dis10 = zeros(1,2^(modu-1)); % vector for Euclidian distance Eucli_dis11 = zeros(1,2^(modu-1)); Eucli_dis20 = zeros(1,2^(modu-1)); Eucli_dis21 = zeros(1,2^(modu-1)); Eucli_dis30 = zeros(1,2^(modu-1)); Eucli_dis31 = zeros(1,2^(modu-1)); Eucli_dis40 = zeros(1,2^(modu-1)); Eucli_dis41 = zeros(1,2^(modu-1)); Eucli_dis50 = zeros(1,2^(modu-1)); Eucli_dis51 = zeros(1,2^(modu-1)); Eucli_dis60 = zeros(1,2^(modu-1)); Eucli_dis61 = zeros(1,2^(modu-1)); constel = zeros(1,2^modu); % constellations counter1 = 1; % counters counter2 = 1; counter3 = 1; counter4 = 1; counter5 = 1; counter6 = 1; counter7 = 1; counter8 = 1; counter9 = 1; counter10 = 1; counter11 = 1; counter12 = 1; if modu == 1 constel = exp(j*[0 -pi]); constel = constel([1 0]+1); elseif modu == 2 constel = exp(j*[-3/4*pi 3/4*pi 1/4*pi -1/4*pi]); constel = constel([2 3 1 0]+1); elseif modu == 4 counter_1 = 1; for counter_2 = -3:2:3 for counter_3 = -3:2:3 constel(counter_1) = (counter_2+j*counter_3)/sqrt(10); % power normalization counter_1 = counter_1 + 1; end; end; constel = constel([0 1 3 2 4 5 7 6 12 13 15 14 8 9 11 10]+1); elseif modu == 6 counter_1 = 1; for counter_2 = -7:2:7 for counter_3 = -7:2:7 constel(counter_1) = (counter_2+j*counter_3)/sqrt(42); % power normalization counter_1 = counter_1 + 1; end; end; constel = constel([[ 0 1 3 2 7 6 4 5]... 8+[ 0 1 3 2 7 6 4 5]... 24+[ 0 1 3 2 7 6 4 5]... 16+[ 0 1 3 2 7 6 4 5]... 56+[ 0 1 3 2 7 6 4 5]... 48+[ 0 1 3 2 7 6 4 5]... 32+[ 0 1 3 2 7 6 4 5]... 40+[ 0 1 3 2 7 6 4 5]]+1); else error('Unimplemented modulation'); end %-------------------------------computer first bit LLR value--------------- if modu == 2 for counter0 = 1:2^modu if real(constel(1,counter0)) > 0 Eucli_dis10(1,counter1) = abs(bits_modu-constel(counter0)); counter1 = counter1+1; else Eucli_dis11(1,counter2) = abs(bits_modu-constel(counter0)); counter2 = counter2+1; end end %-------------------------------computer second bit LLR value-------------- for counter0 = 1:2^modu if imag(constel(1,counter0)) > 0 Eucli_dis20(1,counter3) = abs(bits_modu-constel(counter0)); counter3 = counter3+1; else Eucli_dis21(1,counter4) = abs(bits_modu-constel(counter0)); counter4 = counter4+1; end end LLR0 = min( Eucli_dis10)/min( Eucli_dis11); LLR1 = min( Eucli_dis20)/min( Eucli_dis21); detected = [log(LLR0) log(LLR1)] ; elseif modu == 4 %16QAM---LLR for counter0 = 1:2^modu %第一个比特 %--------------------find the first bit is 0 in the constel---------------------------------- if real(constel(1,counter0)*sqrt(10)) == -1 Eucli_dis10(1,counter1) = abs(bits_modu-constel(counter0)); counter1 = counter1+1; end if real(constel(1,counter0)*sqrt(10)) == -3 Eucli_dis10(1,counter1) = abs(bits_modu-constel(counter0)); counter1 = counter1+1; end %------------------find the first bit is 1 in the constel--------------------------------------- if real(constel(1,counter0)*sqrt(10)) == 1 Eucli_dis11(1,counter2) = abs(bits_modu-constel(counter0)); counter2 = counter2+1; end if real(constel(1,counter0)*sqrt(10)) == 3 Eucli_dis11(1,counter2) = abs(bits_modu-constel(counter0)); counter2 = counter2+1; end end %----------------------computer the first bit------------------------------------------------ LLR1 = min(Eucli_dis10)/min(Eucli_dis11); %-------------------------------------------------------------------------------------------- %第二个比特 for counter0 = 1:2^modu %--------------------find the second bit is 0 in the constel-------------------------------------- if real(constel(1,counter0)*sqrt(10)) == 3 Eucli_dis20(1,counter3) = abs(bits_modu-constel(counter0)); counter3 = counter3+1; end if real(constel(1,counter0)*sqrt(10)) == -3 Eucli_dis20(1,counter3) = abs(bits_modu-constel(counter0)); counter3 = counter3+1; end %--------------------find the second bit is 1 in the constel------------------------------------ if real(constel(1,counter0)*sqrt(10)) == 1 Eucli_dis21(1,counter4) = abs(bits_modu-constel(counter0)); counter4 = counter4+1; end if real(constel(1,counter0)*sqrt(10)) == -1 Eucli_dis21(1,counter4) = abs(bits_modu-constel(counter0)); counter4 = counter4+1; end end %----------------------computer the second bit------------------------------------------- LLR2 = min(Eucli_dis20)/min(Eucli_dis21); %---------------------------------------------------------------------------------------- %第三个比特 for counter0 = 1:2^modu %--------------------find the third bit is 0 in the constel--------------------------------- if imag(constel(1,counter0)*sqrt(10)) == -1 Eucli_dis30(1,counter5) = abs(bits_modu-constel(counter0)); counter5 = counter5+1; end if imag(constel(1,counter0)*sqrt(10)) == -3 Eucli_dis30(1,counter5) = abs(bits_modu-constel(counter0)); counter5 = counter5+1; end %--------------------find the third bit is 1 in the constel-------------------------------- if imag(constel(1,counter0)*sqrt(10)) == 1 Eucli_dis31(1,counter6) = abs(bits_modu-constel(counter0)); counter6 = counter6+1; end if imag(constel(1,counter0)*sqrt(10)) == 3 Eucli_dis31(1,counter6) = abs(bits_modu-constel(counter0)); counter6 = counter6+1; end end %----------------------computer the third bit-------------------------------------------- LLR3 = min(Eucli_dis30)/min(Eucli_dis31); %---------------------------------------------------------------------------------------- %第四个比特 for counter0 = 1:2^modu %--------------------find the fouth bit is 0 in the constel------------------------------ if imag(constel(1,counter0)*sqrt(10)) == 3 Eucli_dis40(1,counter7) = abs(bits_modu-constel(counter0)); counter7 = counter7+1; end if imag(constel(1,counter0)*sqrt(10)) == -3 Eucli_dis40(1,counter7) = abs(bits_modu-constel(counter0)); counter7 = counter7+1; end %--------------------find the fouth bit is 0 in the constel--------------------------- if imag(constel(1,counter0)*sqrt(10)) == 1 Eucli_dis41(1,counter8) = abs(bits_modu-constel(counter0)); counter8 = counter8+1; end if imag(constel(1,counter0)*sqrt(10)) == -1 Eucli_dis41(1,counter8) = abs(bits_modu-constel(counter0)); counter8 = counter8+1; end end LLR4 = min(Eucli_dis40)/min(Eucli_dis41); detected = [log(LLR1) log(LLR2) log(LLR3) log(LLR4)]; else %64QAM软比特信息 %------------------the first bit------------------------------ for counter0 = 1:2^modu if real(constel(counter0)*sqrt(42)) < 0 Eucli_dis10(1,counter1) = abs(bits_modu-constel(counter0)); counter1 = counter1+1; else Eucli_dis11(1,counte