adaptive_blind.zip_Blind OFDM_adaptive-blind_blind_ofdm_ofdm bli

% % SISO-OFDM Subcarrier-Adaptive % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [ber, per, throughput, mean_SNRdB, Mode_err] = ... total_sim_adaptive_blind(info_byte, Cycle, Coding_Rate, SNR, Model, Doppler_Speed, SNR_Th, sub_in_group, Blind_Detection) %%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%/ %%%%%%%%%%%%%%%%%%%%% rand('state',0); sub_per = 0; sub_mean_SNRdB = zeros(1,Cycle); sub_ber = zeros(1,Cycle); sub_CBPS = zeros(1,Cycle); sub_DBPP = zeros(1,Cycle); packet_Mode_err = 0; Num_trans_packet = Cycle; act_sym = 6; MAC_header = MAC_head; % LLC v0_kmh = Doppler_Speed; CarrierFrequency_Hz = 5.25e9; Wavelength_m = 3e8/CarrierFrequency_Hz; v0_ms = v0_kmh/3.6; f_D_Hz = v0_ms/Wavelength_m; fid = fopen('mode_est_error.txt','w'); fclose(fid); switch(Model) case 0 Multipath_Delay = 0; Multipath_Gain = 0; case 1 Multipath_Delay = [0 10 20 30 40 50 60 70 80 90 110 140 170 200 240 290 340 390] * 10^(-9); Multipath_Gain = [0.0 -0.9 -1.7 -2.6 -3.5 -4.3 -5.2 -6.1 -6.9 -7.8 -4.7 -7.3 -9.9 -12.5 -13.7 -18.0 -22.4 -26.7]; case 2 Multipath_Delay = [0 10 20 30 50 80 110 140 180 230 280 330 400 490 600 730 880 1050] * 10^(-9); Multipath_Gain = [-3.3 -3.6 -3.9 -4.2 0.0 -0.9 -1.7 -2.6 -1.5 -3.0 -4.4 -5.9 -5.3 -7.9 -9.4 -13.2 -16.3 -21.2]; case 3 Multipath_Delay = [0 10 20 40 70 100 140 190 240 320 430 560 710 880 1070 1280 1510 1760] * 10^(-9); Multipath_Gain = [-4.9 -5.1 -5.2 -0.8 -1.3 -1.9 -0.3 -1.2 -2.1 0.0 -1.9 -2.8 -5.4 -7.3 -10.6 -13.4 -17.4 -20.9]; otherwise error('Not prepared for this channel model.'); end ch = rayleighchan(1/(80*10^6), f_D_Hz, Multipath_Delay, Multipath_Gain); ch.ResetBeforeFiltering = 0; % ------------ Mode -------------- for Z=0:Cycle %Mapping decision if Z == 0 Mode = 3*ones(1,480); % Num_ofdm_sym = cal_symbol_adaptive(info_byte, Mode, Coding_Rate); else [Mode, SNRdB_each] = criterion2_group(feed_back_CSI, pwr_n, SNR_Th, sub_in_group); sub_mean_SNRdB(Z) = mean(SNRdB_each); Num_ofdm_sym = cal_symbol_adaptive(info_byte, Mode, Coding_Rate); end if ((Num_ofdm_sym >= 1) && (Num_ofdm_sym <= 70)) source_dat = source_generator_adaptive_code(info_byte, MAC_header, Num_ofdm_sym, Mode, Coding_Rate); else source_dat(1,1) = 'n'; end if char(source_dat(1,1)) == 'n' demap_dat = 110*ones(1,Num_ofdm_sym); trans_dat = []; Num_trans_packet = Num_trans_packet - 1; else coded_dat = conv_encord_adaptive(source_dat, Coding_Rate); inter_dat = interleave_adaptive_new(coded_dat, Mode);%cmdadd for I=1:Num_ofdm_sym map_dat(:,I) = mapping_adaptive(inter_dat(:,I), Mode); end pilot_add=[zeros(1,Num_ofdm_sym); ... map_dat(1:22,:); ... ones(1,Num_ofdm_sym); ... map_dat(23:44,:); ... ones(1,Num_ofdm_sym); ... map_dat(45:66,:); ... ones(1,Num_ofdm_sym); ... map_dat(67:88,:); ... -ones(1,Num_ofdm_sym); ... map_dat(89:110,:); ... ones(1,Num_ofdm_sym); ... map_dat(111:132,:); ... -ones(1,Num_ofdm_sym); ... map_dat(133:154,:); ... ones(1,Num_ofdm_sym); ... map_dat(155:176,:); ... ones(1,Num_ofdm_sym); ... map_dat(177:198,:); ... -ones(1,Num_ofdm_sym); ... map_dat(199:220,:); ... ones(1,Num_ofdm_sym); ... map_dat(221:240,:); ... zeros(11,Num_ofdm_sym); ... map_dat(241:260,:); ... ones(1,Num_ofdm_sym); ... map_dat(261:282,:); ... ones(1,Num_ofdm_sym); ... map_dat(283:304,:); ... -ones(1,Num_ofdm_sym); ... map_dat(305:326,:); ... ones(1,Num_ofdm_sym); ... map_dat(327:348,:); ... ones(1,Num_ofdm_sym); ... map_dat(349:370,:); ... ones(1,Num_ofdm_sym); ... map_dat(371:392,:); ... -ones(1,Num_ofdm_sym); ... map_dat(393:414,:); ... -ones(1,Num_ofdm_sym); ... map_dat(415:436,:); ... ones(1,Num_ofdm_sym); ... map_dat(437:458,:); ... ones(1,Num_ofdm_sym); ... map_dat(459:480,:)]; ifft_dat = ifft(pilot_add, 512); gi_add = [ifft_dat(449:512,:); ifft_dat]; trans_dat = reshape(gi_add,576*Num_ofdm_sym,1); end if (Multipath_Delay==0) multipath_add = trans_dat; else multipath_add = filter(ch, trans_dat); end signal_power = mean(abs(trans_dat) .* abs(trans_dat)); SNR0 = SNR -10*log10(512/480); noise_add = awgn(multipath_add, SNR0, 10*log10(signal_power), 'dB'); gi_rem = synchro_dat(65:576,:); fft_dat = fft(gi_rem, 512); pilot_rem = rem_pilot(fft_dat); rec_pilot = ex_pilot(fft_dat); end rec_pream_fft1=fft(noise_add(288+1:800), 512); rec_pream_fft2=fft(noise_add(800+1:1312), 512); rec_pream_fft_ave = (rec_pream_fft1 + rec_pream_fft2) ./2; source_pream_fft=fft(long_pre(128+1:640), 512); data_shift = rec_pream_fft_ave ./ (source_pream_fft + 10^(-12)); feed_back_CSI = rem_pilot(data_shift); pilot_h = ex_pilot(data_shift); h = feed_back_CSI; pwr_n = mean(abs(rec_pream_fft1-rec_pream_fft2).^2)/2; if Z>=1 && char(source_dat(1,1)) ~= 'n' Mode_est = mode_estimator_rev2(pilot_rem, h, sub_in_group, act_sym, Mode, SNR_Th, Blind_Detection); end if Z>=1 && char(source_dat(1,1)) ~= 'n' temp_Mode_err = 0; for I=1:480 if Mode_est(I) ~= Mode(I) temp_Mode_err = temp_Mode_err + 1; end end if sum(temp_Mode_err) == 0 for I=1:Num_ofdm_sym equ_dat(:,I) = pilot_rem(:,I)./h; end for I = 1:Num_ofdm_sym demap_dat(:,I) = soft_demap_adaptive_fast_CSI(equ_dat(:,I), Mode, h, pwr_n); end deinter_dat = deinterleave_adaptive_new(demap_dat,Mode); for I=1:Num_ofdm_sym viterbi_dat(:,I)=soft_viterbi_adaptive_fast(deinter_dat(:,I), Mode, Coding_Rate); end sub_ber(Z) = ber_calc(viterbi_dat, source_dat); if(sub_ber(Z) ~= 0) sub_per = sub_per + 1; end else % In case of Mode Estimation Error packet_Mode_err = packet_Mode_err + 1; sub_ber(Z) = 0.5; sub_per = sub_per + 1; end end y_ = filter(ch, ones(1,2000)); clear source_dat coded_dat inter_dat map_dat synchro_dat equ_dat demap_dat fading_add deinter_dat viterbi_dat; end if Num_trans_packet >= 1 Mode_err = packet_Mode_err/Num_trans_packet; ber = sum(sub_ber)/Num_trans_packet; per = sub_per/Num_trans_packet; CBPS = sum(sub_CBPS)/Num_trans_packet; DBPS = Coding_Rate * CBPS; throughput = sum(sub_DBPP)/((Num_ofdm_sym*Num_trans_packet*7.2+Cycle*16.4)*10^(-6)); % Mbps ブラインドのためシグナルシンボルなし mean_SNRdB = sum(sub_mean_SNRdB)/Cycle; else Mode_err = 1; ber = 1; per = 1; DBPS = 0; throughput = 0; mean_SNRdB = sum(sub_mean_SNRdB)/Cycle; end