Диссертация (Снижение пик-фактора неортогональных многочастотных сигналов путем добавления корректирующих поднесущих), страница 13

B. S. Krongold and D. L. Jones, "An Active-Set approach for OFDM PARreduction via tone reservation," IEEE Trans. on Signal Processing, vol. 52, no. 2, pp.495–509, Feb. 2004.13. Феер К. Беспроводная цифровая связь. Методы модуляции и расширения спектра. – М.: Радио и связь, 2000. – 502 с.14. X.

Yang, W. Ai, T. Shuai, D. Li, “A fast decoding algorithm for non-orthogonal frequency division multiplexing signals,” International Conference on Communications and Networking in China (CHINACOM), pp. 595–598, Aug. 2007.15. I. Kanaras, A. Chorti, M. Rodrigues, I. Darwazeh, “An overview of optimaland sub-optimal detection techniques for a non-orthogonal spectrally efficient FDM,”International Symposium on Communication and Information Technologies, pp. 460–465, Sept.

2009.16. S. Bharadwaj, N. Krishna, P. Sudheesh, M. Jayakumar, “Low complexity detection scheme for NOFDM systems based on ML detection over hyperspheres”, International Conference on Devices and Communications (ICDeCom), pp. 1–5, Feb.2011.17. IEEE Standard for Information technology – Telecommunications and information exchange between systems Local and metropolitan area networks – Specificrequirements, IEEE Standard 802.11, 2012.18. Telesystem Innovations Inc., LTE in a Nutshell: The Physical Layer, 2010.19. Digital Video Broadcasting (DVB); Frame structure channel coding and modulation for a second generation digital terrestrial television broadcasting system (DVBT2), ETSI Standard 302 755 V1.4.1, 2015.9620.

Safa Isam A Ahmed. "Spectrally Efficient FDM Communication Signals andTransceivers: Design, Mathematical Modelling and System Optimization", A thesissubmitted for the degree of Doctor of Philosophy, 2011, pp 182–201.21. L. Yang, R. Chen, K. Soo, and Y. Siu, “An efficient sphere decoding approachfor PTS assisted PAPR reduction of OFDM signals,” AEU - International Journal ofElectronics and Communications, vol. 61, pp. 684–688, 2007.22.

J. Cimini, L.J. and N. Sollenberger, “Peak-to-average power ratio reductionof an OFDM signal using partial transmit sequences,” IEEE Communications Letters,vol. 4, no. 3, pp. 86–88, 2000.23. P. Boonsrimuang and T. Paungma, “Proposal of improved PTS method forOFDM signal in the multi-path fading channel,” in Proc. 5th International Conferenceon List of References 306 Electrical Engineering/Electronics, Computer, Telecommunications and Informa- tion Technology ECTI-CON 2008, vol. 1, May 14–17, 2008,pp.

401–404.24. B. S. Krongold and D. L. Jones, "An Active-Set approach for OFDM PARreduction via tone reservation," IEEE Trans. on Signal Processing, vol. 52, no. 2, pp.495–509, Feb. 2004.25. Xilinx LogiCORE IPFast Fourier Transform v8.0 (http://www.xilinx.com/support/documentation/ip_documentation/ds808_xfft.pdf).26. Digital Video Broadcasting (DVB); Frame structure channel coding and modulation for a second generation digital terrestrial television broadcasting system (DVBT2), ETSI Standart 302 755 V1.1.1, 2009.27.

Нгок Нуен Тан, Антонов Е.О., Фадеев Д.К., Рашич А.В. Снижение пикфактора многочастотных сигналов с неортогональным частотным уплотнениемметодом резервирования поднесущих частот // Радиотехника, №1, 2018 г., с.1623.28. Нгок Н. Тан, Рашич А. В., Сальников В.А. Расчет пик-фактора многочастотных сигналов с ортогональным и неортогональным частотным уплотнением// НТВ СПбГПУ. Информатика. Телекоммуникации.

Управление №3, 2018. Принята в печать9729. Ngok Nuen Tan, Alexandr B. Kislitsyn, Andrey V. Rashic, Generation ofSEFDM-Signals Using FFT/IFFT // 14th International Conference, NEW2AN 2014and 7th Conference, ruSMART 2014, Proceedings, 8638 LNCS, pp. 488-501.30. Ngoc Tan, Antonov E.O., Fadeev D.K., Rashich A.V. Reduced ComplexityTone Reservation Peak-to-Average Power Ratio Reduction Algorithm for SEFDMSignals // 39th International conference on telecommunications and signal processing,2016.31. Ngoc Tan, Kislitsyn A.B., Fadeev D.K., Rashich A.V. FFT-based trellis receiver for SEFDM signals // 2016 IEEE Global Communications Conference(GLOBECOM); Pages: 1 - 6, DOI: 10.1109/GLOCOM.2016.784184132. Нгок Нуен Тан, Антонов Е.О., Фадеев Д.К., Рашич А.В. Пик-факторSEFDM-сигналов // 18-я Международная конференция.

Цифровая обработка сигналов и ее применение. Доклады Выпуск: XVIII; Том-1. – М.: ООО «БРИС-М»,2016. – 356-362 с.33. Antognetti P. Power Integrated Circuits: Physics, Design, and Applications.New York, McGraw-Hill, 1986. 544 p.34. Макаров С.Б., Марков А.М.

Влияние амплитудного ограничения квадратурных составляющих сигналов с OFDM в модуляторе радиопередатчика на помехоустойчивость приема // Научно-технические ведомости СПбГПУ. 2012, №2(145).35. Thompson S. C., Proakis J. G., Zeidler J.

R. The Effectiveness of Signal Clipping for PAPR and Total Degradation Reduction in OFDM Systems // Global Telecommunications Conference, 2005. GLOBECOM '05. IEEE.36. Guel, D. Palicot, J. Analysis and comparison of clipping techniques forOFDM Peak-to-Average Power Ratio reduction // Digital Signal Processing, 2009 16thInternational Conference.98Приложение 1. Исходный код программы разработанногоалгоритма снижения пик-фактора SEFDM-сигналов%% ============================================================================function PAPR_stat = PAPR_calc( init_params )% exp_par - number of concurrent experiments (1 to X), 1 - one concurrentexperiment% exp_seq - number of concurrent attempts (1 to Nattempt). The less thisnumber the less%needed RAM for calculation, but more time.%% InitializationMod_info= init_params.Mod_info;% AlphabetMod_res= init_params.Mod_res;% Reserved subcarriers alphaabetNexp= init_params.Nexp;% Experiment numalpha= init_params.alpha;% Normalized subcarier frequency spac-Nattempt= init_params.Nattempt;% Attempt numgpu_calc= init_params.gpu_calc;% calc using gpuexp_par= init_params.exp_par;% Number of concurrent experimentsexp_seq= init_params.exp_seq;% Number of concurrent attemptsNinfo= init_params.Ninfo;% Info subarriers numNfft= init_params.Nfft;% IFFT volNfft_res= init_params.Nfft_res;% IFFT vol reserved subcarriers (2^p)ing% --rng( 'shuffle' );% Reset the CPU random number generatorparallel.gpu.rng( 'shuffle' ); % Reset the GPU random number generator%-------------------------------------------------------------------------------% Additinal init parametersgamma_apsk = [];constell_info= get_constell_v2(Mod_info,gamma_apsk); % Get constellationvectorlen_const_info= length(constell_info); % alphabet volumeconstell_res= get_constell_v2(Mod_res,gamma_apsk); % Get constellationvectorlen_const_res= length(constell_res);% Energy normalization to E = 199E_info= sum(abs(constell_info).^2)/length(constell_info);constell_info= constell_info/sqrt(E_info);E_res= sum(abs(constell_res).^2)/length(constell_res);constell_res= constell_res/sqrt(E_res);%----------------------frag_num = Nfft/Nfft_res;% replication num IFFT reservedaddr_res = 1:frag_num:Nfft;% addresses of reserved subcarriersinfo_sub_places = true(Nfft,1);info_sub_places(addr_res) = false; % true only on info placesif( mod(Ninfo,2) ~= 0 )% odd or even info subcarriers numNright= (Nfft-1-Ninfo)/2;% left guard intervalNleft= (Nfft-1-Ninfo)/2+1;% right guard intervalnew_g_l = (Nfft-1-Nfft_res-Ninfo)/2;new_g_r = (Nfft-1-Nfft_res-Ninfo)/2+1;elseNright= (Nfft-Ninfo)/2;Nleft= Nright;new_g_l = (Nfft-Nfft_res-Ninfo)/2; % new guard interval left (only infosubcarrier places)new_g_r = new_g_l;end%----------------------% zero res subcarrier places% Example: ooo:|:::|:::|:::|:::|:::|:::|ooo, o -- zeros, | -- reserved subcarriers,% : -- info subcarriers (Nfft = 32, Nfft_res = 8)cut_pl = false(Nfft,1);cut_pl(info_sub_places) = [false(new_g_l,1); true(Ninfo,1);false(new_g_r,1)]; % true on positions of info subcarriersNres_cut_l= sum(~info_sub_places(1:find(cut_pl,1,'first')-1));% num ofzero res subcarriers on the right from info subNres_cut_r= sum(~info_sub_places(1+find(cut_pl,1,'last'):end)); % on theleftNres_interl = Nfft_res - (Nres_cut_l + Nres_cut_r);% num ofres interleaving with info subcarriers% --PAPR_stat_w_o_red= zeros(Nexp,1); % w/o reductionPAPR_stat_zero_dec = zeros(Nexp,1); % zero decimatedPAPR_stat_res= zeros(Nexp,1); % with reserved subcarriersfprintf(['----------------\n'....'Nfft= %d\n'....'Ninfo= %d\n'....100'Nexp= %.1e\n'....'Nattempt = %d\n'....'Nfft_res = %d\n'....'Mod_info = %s\n'....'Mod_res= %s\n'....'alpha= %.2f\n'....'exp_par= %d\n'],...Nfft,...Ninfo,...Nexp,...Nattempt,...Nfft_res,...Mod_info,...Mod_res,...alpha,...exp_par);if( gpu_calc )constell_info= gpuArray(constell_info);constell_res= gpuArray(constell_res);enditer_show = 1;%% Calcif( gpu_calc )for i = 0:ceil(Nexp/exp_par)-1word_info= gpuArray.randi( [0 len_const_info-1], Ninfo, exp_par);word_mod_info = constell_mapping_nums_in( word_info, constell_info);% PAPR calc w/o reduction and decimationword_mod_gi= [zeros(Nleft,exp_par); word_mod_info; ze-ros(Nright,exp_par)]; % with guard intervalst_symb_no= ifft(word_mod_gi,Nfft); % OFDM time domaint_symb_no= t_symb_no(1:ceil(Nfft*alpha),:); % SEFDMP= t_symb_no.*conj(t_symb_no);PAPR_stat_w_o_red(i*exp_par+1 : i*exp_par+exp_par) =gather(max(P,[],1)./mean(P,1));% Zero decimated info vectorinfo_wz101= gpuArray(zeros(Nfft,exp_par));info_wz(repmat(info_sub_places,1,exp_par)) = [zeros(new_g_l,exp_par); word_mod_info; zeros(new_g_r,exp_par)]; % info subcarrierswith zerosinfo_wz_time= ifft(info_wz,Nfft); % with zeros in timeti_symb_wz= info_wz_time(1:ceil(Nfft*alpha),:); %P= ti_symb_wz.*conj(ti_symb_wz);domainSEFDMPAPR_stat_zero_dec(i*exp_par+1 : i*exp_par+exp_par) =gather(max(P,[],1)./mean(P,1)); % zero decimated info vectorfor exp_seq_iter = 0:ceil(Nattempt/exp_seq)-1word_res= gpuArray.randi( [0 len_const_res-1], Nres_in-terl, ceil(Nattempt*exp_par/exp_seq) );word_mod_res= constell_mapping_nums_in( word_res, constell_res);% Addition of reserved subcarriersinfo_wz_time_att = reshape(repmat(info_wz_time.',1,exp_seq).',Nfft,[]); % replicate according to exp_seq, i.e.cols: 1 2 => 1 1 2 2word_mod_res_gi= [zeros(Nres_cut_l, exp_seq*exp_par);word_mod_res; zeros(Nres_cut_r, exp_seq*exp_par)];word_res_time= ifft(word_mod_res_gi,Nfft_res);word_res_time= 1/frag_num*repmat(word_res_time,frag_num,1);t_symb_ins= info_wz_time_att + word_res_time;t_symb_ins= t_symb_ins(1:ceil(Nfft*alpha),:); % SEFDMP= t_symb_ins.*conj(t_symb_ins);PAPR_res_exp_par(exp_seq_iter*exp_seq+1 :exp_seq_iter*exp_seq+exp_seq)....= max(P,[],1)./mean(P,1); % calc papr for all coloumnsendfor j = 0:exp_par-1PAPR_stat_res(i*exp_par+1+j) =gather(min(PAPR_res_exp_par(j*Nattempt+1:j*Nattempt+Nattempt))); % with reservedendfprintf('%d ',i);if mod(iter_show,100) == 0102fprintf('\n');endif mod(iter_show,100e3) == 0save('temp.mat','PAPR_stat_w_o_red','PAPR_stat_zero_dec','PAPR_stat_res');fprintf('SAVED\n');enditer_show = iter_show + 1;endelsefor i = 1:exp_par:Nexpword_info= randi( [0 len_const_info-1], Ninfo, exp_par );word_mod_info = constell_mapping_nums_in( word_info, constell_info);% PAPR calc w/o reduction and decimationword_mod_gi= [zeros(Nleft,exp_par); word_mod_info; ze-ros(Nright,exp_par)]; % with guard intervalst_symb_no= ifft(word_mod_gi,Nfft); % OFDM time domaint_symb_no= t_symb_no(1:ceil(Nfft*alpha),:); % SEFDMP= t_symb_no.*conj(t_symb_no);PAPR_stat_w_o_red(i*exp_par+1 : i*exp_par+exp_par) =gather(max(P,[],1)./mean(P,1));% Zero decimated info vectorinfo_wz= zeros(Nfft,exp_par);info_wz(repmat(info_sub_places,1,exp_par)) = [zeros(new_g_l,exp_par); word_mod_info; zeros(new_g_r,exp_par)]; % info subcarrierswith zerosinfo_wz_time= ifft(info_wz,Nfft); % with zeros in timeti_symb_wz= info_wz_time(1:ceil(Nfft*alpha),:); %P= ti_symb_wz.*conj(ti_symb_wz);domainSEFDMPAPR_stat_zero_dec(i*exp_par+1 : i*exp_par+exp_par) =max(P,[],1)./mean(P,1); % zero decimated info vectorfor exp_seq_iter = 0:ceil(Nattempt/exp_seq)-1word_res= randi( [0 len_const_res-1], Nres_interl,ceil(Nattempt*exp_par/exp_seq) );word_mod_res);103= constell_mapping_nums_in( word_res, constell_res% Addition of reserved subcarriersinfo_wz_time_att = reshape(repmat(info_wz_time.',1,exp_seq).',Nfft,[]); % replicate according to exp_seq, i.e.cols: 1 2 => 1 1 2 2word_mod_res_gi= [zeros(Nres_cut_l, exp_seq*exp_par);word_mod_res; zeros(Nres_cut_r, exp_seq*exp_par)];word_res_time= ifft(word_mod_res_gi,Nfft_res);word_res_time= 1/frag_num*repmat(word_res_time,frag_num,1);t_symb_ins= info_wz_time_att + word_res_time;t_symb_ins= t_symb_ins(1:ceil(Nfft*alpha),:); % SEFDMP= t_symb_ins.*conj(t_symb_ins);PAPR_res_exp_par(exp_seq_iter*exp_seq+1 :exp_seq_iter*exp_seq+exp_seq)....= max(P,[],1)./mean(P,1); % calc papr for all coloumnsendfor j = 0:exp_par-1PAPR_stat_res(i*exp_par+1+j) = min(PAPR_res_exp_par(j*Nattempt+1:j*Nattempt+Nattempt)); % with reservedendfprintf('%d ',i);if mod(iter_show,100) == 0fprintf('\n');endif mod(iter_show,100e3) == 0save('temp.mat','PAPR_stat_w_o_red','PAPR_stat_zero_dec','PAPR_stat_res');fprintf('SAVED\n');enditer_show = iter_show + 1;endendfprintf('\n');PAPR_stat.w_o_red= PAPR_stat_w_o_red;PAPR_stat.zero_dec= PAPR_stat_zero_dec;PAPR_stat.res= PAPR_stat_res;end104.