【matlab编程代做】MATLAB实现基于无线宽带系统的正交频分复用(OFDM)的链路级和系统级仿真

% this function load the run-time data, save in a single mat file% called by script_load_part2% after finishing simulating each case, this function need to be called% all the analysis is based on the output of this functionfunctionload_part2(test_resource_allocation_idx,UE,MIMO,loaddir,savedir,a,b,overwrite_old,MaxSaveInd ex,save_info_factor,Comm_conf)filename = ['RA',num2str(test_resource_allocation_idx),'_UE',num2str(UE),'_FLA',num2str(Comm_conf.F_LA ),'_FPA',num2str(Comm_conf.F_PA),MIMO,'_trms',num2str(Comm_conf.t_rms),'f',num2str(Com m_conf.fd),'_s',num2str(Comm_conf.state),'N',num2str(Comm_conf.N),'subN',num2str(Comm_co nf.subN),'_E',num2str(Comm_conf.flg_csi_err),'c',num2str(Comm_conf.flg_offset),'D',num2str(Co mm_conf.DELAY),'BO',num2str(Comm_conf.BO_dB),'HPA',num2str(Comm_e_HPA_flg),'_B LER',num2str(Comm_conf.targetFER),'.mat'];if exist([savedir,'/',filename],'file') && ~overwrite_oldreturn;elsetest_possible_user_velocities = 30;test_possible_scenarios = {'micro1','macro1','macro3'};test_pissible_resource_allocations = {'RR','PF','MAXCI','PF2'};test_nu_active_users = UE;test_scenario = test_possible_scenarios{1};test_resource_allocation = test_pissible_resource_allocations{test_resource_allocation_idx};test_avg_velo_kmph = test_possible_user_velocities(1);%test_frame_period_factor = LARArate;test_case=['testCaseResult-Scn',test_scenario,'-ResAlloc',test_resource_allocation,'-PowCon','no', '-G','no','-NoUE',num2str(test_nu_active_users),'-kmph',num2str(test_avg_velo_kmph),'F_LA',nu m2str(Comm_conf.F_LA),'-F_PA',num2str(Comm_conf.F_PA)];save_user_fair = zeros(UE,save_info_factor * (b-a + 1) * MaxSaveIndex);save_user_ActiveCnt = uint16(zeros(UE,save_info_factor * (b-a + 1) * MaxSaveIndex));save_user_distance = save_user_ActiveCnt;save_user_mobID = save_user_ActiveCnt;save_user_velocity = uint8(save_user_ActiveCnt);if Comm_conf.N == 512,Comm_ed = 301;elseif Comm_conf.N == 1024,Comm_ed = 601;endNu_groups = floor(Comm_ed / Comm_conf.subN);save_user_sinr = int8(zeros(UE,save_info_factor * (b-a + 1) * MaxSaveIndex,Nu_groups));save_bs_sinr = int8(zeros(Nu_groups,save_info_factor * (b-a + 1) * MaxSaveIndex));save_bs_mobID = uint16(zeros(Nu_groups,save_info_factor * (b-a + 1) * MaxSaveIndex));save_bs_mod = int8(zeros(Nu_groups,save_info_factor * (b-a + 1) * MaxSaveIndex));save_bs_fec = int8(zeros(Nu_groups,save_info_factor * (b-a + 1) * MaxSaveIndex));save_bs_poweradjust = (zeros(Nu_groups,save_info_factor * (b-a + 1) * MaxSaveIndex));save_bs_FrminErr = int8(zeros(Nu_groups,save_info_factor * (b-a + 1) * MaxSaveIndex));save_bs_dist = uint16(zeros(Nu_groups,save_info_factor * (b-a + 1) * MaxSaveIndex));save_outageCnt = 0;save_totalBitsSent = 0;save_totalFrmErr = 0;nuUser = 0;for lp_nuDrops=a:bdir_name=[loaddir,'/',test_case,'Drp',num2str(lp_nuDrops),'_trms',num2str(Comm_conf.t_rms),'f' ,num2str(Comm_conf.fd),'_s',num2str(Comm_conf.state),'N',num2str(Comm_conf.N),'subN',num 2str(Comm_conf.subN),'_E',num2str(Comm_conf.flg_csi_err),'c',num2str(Comm_conf.flg_offset),' D',num2str(Comm_conf.DELAY),'BO',num2str(Comm_conf.BO_dB),'HPA',num2str(Comm_ e_HPA_flg),'_BLER',num2str(Comm_conf.targetFER)];lp_nuDropsfor k = 1:MaxSaveIndexload([dir_name,'/','user_gain_info_and_distance/saveIndex',num2str(k)]);idx = (lp_nuDrops-1) * MaxSaveIndex*save_info_factor + (k-1)*save_info_factor + 1: (lp_nuDrops-1) * MaxSaveIndex*save_info_factor + k*save_info_factor;save_user_ActiveCnt(:,idx) = user_ActiveCnt;save_user_distance(:,idx) = user_distance;save_user_mobID(:,idx) = user_mobID + nuUser;save_user_velocity(:,idx) = user_velocity;save_user_sinr(:,idx,:) = user_sinr;save_user_fair(:,idx) = user_fair;%save_user_sinr_noErr(:,idx,:) = user_sinr_noFeedBackErr;clear user*load([dir_name,'/','measurement/dataVals',num2str(k)]);%save_user_remFileSz(:,idx) = st_measurement.remFileSzStore;save_bs_sinr(:,idx) = st_measurement.snrFeedBack;save_bs_mobID(:,idx) = st_measurement.mobIdxStore + nuUser;%save_bs_sinr_noErr(:,idx) = st_measurement.snrFeedBackNoChErr;save_bs_mod(:,idx) = st_measurement.modStore;save_bs_fec(:,idx) = st_measurement.fecStore;save_bs_dist(:,idx) = st_measurement.distanceStore;save_bs_poweradjust(:,idx) = st_measurement.poweradjust;save_bs_FrminErr(:,idx) = st_measurement.FrminError;%save_bs_instTP(1,idx) = st_measurement.instThpt;endnuUser = max(max(save_user_mobID(:,lp_nuDrops * MaxSaveIndex*save_info_factor)));save_outageCnt=save_outageCnt+st_measurement.outageCnt;save_totalBitsSent=save_totalBitsSent+st_measurement.totalBitsSent;save_totalFrmErr = save_totalFrmErr + st_measurement.totalFramesinError;clear st_measurementendfec1 = [0 3 5 7]; fec2 = [0 1/3 1/2 2/3];fecs = zeros(size(save_bs_fec));for k = 2:4fecs = fecs + fec2(k)*(save_bs_fec == fec1(k));endsave_bs_fec = fecs;if ~exist(savedir,'dir'),mkdir(savedir);endsave([savedir,'/',filename],'save_user_ActiveCnt','save_user_distance','save_user_mobID','save_u ser_velocity',...'save_user_sinr','save_bs_sinr',...'save_bs_mod','save_bs_fec','save_bs_dist','save_bs_poweradjust','save_bs_FrminErr','save_bs_ mobID',...'save_outageCnt','save_totalBitsSent','save_user_fair');End% load the data file for each testing case% output is the BLER, Total UEs and Dropped UEs,% the served UE number can be found by total UEs - Dropped UEsfunction [result,fileidx] = res_bler_NoUEOUT(test_resource_allocation_idx,UE,MIMO,loaddir,result,fileidx,Comm_conf) filename = ['RA',num2str(test_resource_allocation_idx),'_UE',num2str(UE),'_FLA',num2str(Comm_conf.F_LA ),'_FPA',num2str(Comm_conf.F_PA),MIMO,'_trms',num2str(Comm_conf.t_rms),'f',num2str(Com m_conf.fd),'_s',num2str(Comm_conf.state),'N',num2str(Comm_conf.N),'subN',num2str(Comm_co nf.subN),'_E',num2str(Comm_conf.flg_csi_err),'c',num2str(Comm_conf.flg_offset),'D',num2str(Co mm_conf.DELAY),'BO',num2str(Comm_conf.BO_dB),'HPA',num2str(Comm_e_HPA_flg),'_B LER',num2str(Comm_conf.targetFER),'.mat'];load([loaddir,'/',filename],'save_outageCnt','save_user_mobID','save_bs_FrminErr','save_bs_mod ');fileidx = fileidx + 1result.bler(fileidx)=sum(sum(save_bs_FrminErr))/sum(sum(save_bs_mod~=0));erPersecond(fileidx) = double(max(max(save_user_mobID)))/(double(size(save_user_mobID,2))*5e-4);result.outagePersecond(fileidx) = double(save_outageCnt)/(size(save_user_mobID,2)*5e-4);% the power utilization figure% generates two figures% figure 1: distance vs. utilized power (per block)% figure 2: CDF of utilized power (for each OFDM frame)function [result,fileidx] = res_bs_powerad(test_resource_allocation_idx,UE,MIMO,loaddir,result,fileidx,av_factor,Comm_c onf)filename = ['RA',num2str(test_resource_allocation_idx),'_UE',num2str(UE),'_FLA',num2str(Comm_conf.F_LA ),'_FPA',num2str(Comm_conf.F_PA),MIMO,'_trms',num2str(Comm_conf.t_rms),'f',num2str(Com m_conf.fd),'_s',num2str(Comm_conf.state),'N',num2str(Comm_conf.N),'subN',num2str(Comm_co nf.subN),'_E',num2str(Comm_conf.flg_csi_err),'c',num2str(Comm_conf.flg_offset),'D',num2str(Co mm_conf.DELAY),'BO',num2str(Comm_conf.BO_dB),'HPA',num2str(Comm_e_HPA_flg),'_B LER',num2str(Comm_conf.targetFER),'.mat'];load([loaddir,'/',filename],'save_bs_poweradjust','save_bs_dist');colors = {[0.5 0.5 0.5],[0 0 0],[0 0.5 1],[0 0 1],[1 0 0.5],[1 0 0]};lst = {'-','-.','-','-.','-','-.'};RAs = {'RR','PF','MX','PF2'};para.RA = RAs{test_resource_allocation_idx};UEs = [4 10];para.UE = UE;para.ra_ix = test_resource_allocation_idx;if Comm_conf.F_LA == 2 && Comm_conf.F_PA == 2,para.ra_ix = 4;endif MIMO == 2,para.ra_ix = 5;endpara.ue_ix = find(UEs == para.UE);para.nuSymbol = 6;para.subN = Comm_conf.subN;para.Tfrm = 5e-4;fileidx = fileidx + 1;%%result.ra_idx(fileidx) = test_resource_allocation_idx;%%fig_ix = 0;dist = reshape(save_bs_dist,[],1);[dist_sort,idx] = sort(dist);tmpV = save_bs_poweradjust;yy=reshape(tmpV,[],1);yy_sort = yy(idx);fig_ix = fig_ix + 1;cfigure(fig_ix);plot(average2(double(dist_sort),av_factor),10*log10(average2(double(yy_sort),av_factor)),'color' ,colors{fileidx},'linestyle',lst{fileidx},'linewidth',2);hold on;result.LGD{fileidx} = [para.RA,',UE',num2str(para.UE),',F',num2str(Comm_conf.F_LA),'P',num2str(Comm_conf.F_PA),',', MIMO];result.figname{fig_ix}='bs_powerad_vs_dist';xlabel('Distance in meter','FontSize',14,'fontname','times','FontWeight','b');ylabel('Utilized Power per Block in dB','FontSize',14,'fontname','times','FontWeight','b');fig_ix = fig_ix + 1;cfigure(fig_ix);yy=mean(tmpV,1);mns = sprintf('%3.2f',10*log10(mean(yy)));result.LGDmn{fileidx} = [para.RA,',UE',num2str(para.UE),',F',num2str(Comm_conf.F_LA),'P',num2str(Comm_conf.F_PA),',', MIMO,',M:',mns];handCDF = cdfplot(10*log10(double(yy)));hold on;set(handCDF,'color',colors{fileidx},'linestyle',lst{fileidx},'linewidth',2);result.figname{fig_ix}='CDF_of_bs_powerad';xlabel('Averaged Utilized Power at NB in dB','FontSize',14,'fontname','times','FontWeight','b'); ylabel('CDF of Utilized Power','FontSize',14,'fontname','times','FontWeight','b');% two figures for the SINR at Base Station% figure 1: distance vs. SINR (per block)% figure 2: CDF of SINR (per block)function [result,fileidx] = res_bs_sinr(test_resource_allocation_idx,UE,MIMO,loaddir,result,fileidx,av_factor,Comm_conf) filename = ['RA',num2str(test_resource_allocation_idx),'_UE',num2str(UE),'_FLA',num2str(Comm_conf.F_LA ),'_FPA',num2str(Comm_conf.F_PA),MIMO,'_trms',num2str(Comm_conf.t_rms),'f',num2str(Com m_conf.fd),'_s',num2str(Comm_conf.state),'N',num2str(Comm_conf.N),'subN',num2str(Comm_co nf.subN),'_E',num2str(Comm_conf.flg_csi_err),'c',num2str(Comm_conf.flg_offset),'D',num2str(Co mm_conf.DELAY),'BO',num2str(Comm_conf.BO_dB),'HPA',num2str(Comm_e_HPA_flg),'_B LER',num2str(Comm_conf.targetFER),'.mat'];load([loaddir,'/',filename],'save_bs_sinr','save_bs_dist');colors = {[0.5 0.5 0.5],[0 0 0],[0 0.5 1],[0 0 1],[1 0 0.5],[1 0 0]};lst = {'-','-.','-','-.','-','-.'};RAs = {'RR','PF','MX','PF2'};para.RA = RAs{test_resource_allocation_idx};UEs = [4 10];para.UE = UE;para.ra_ix = test_resource_allocation_idx;if Comm_conf.F_LA == 2 && Comm_conf.F_PA == 2,para.ra_ix = 4;endif MIMO == 2,para.ra_ix = 5;endpara.ue_ix = find(UEs == para.UE);para.nuSymbol = 6;para.subN = Comm_conf.subN;para.Tfrm = 5e-4;fileidx = fileidx + 1;%%result.ra_idx(fileidx) = test_resource_allocation_idx;result.LGD{fileidx} = [para.RA,',UE',num2str(para.UE),',F',num2str(Comm_conf.F_LA),'P',num2str(Comm_conf.F_PA),',', MIMO];%%fig_ix = 0;dist = reshape(save_bs_dist,[],1);[dist_sort,idx] = sort(dist);tmpV = save_bs_sinr;yy=reshape(tmpV,[],1);yy_sort = yy(idx);fig_ix = fig_ix + 1;cfigure(fig_ix);plot(average2(double(dist_sort),av_factor),average3(double(yy_sort),av_factor),'color',colors{filei dx},'linestyle',lst{fileidx},'linewidth',2);hold on;result.figname{fig_ix}='bs_sinr_vs_dist';xlabel('Distance in Meter','FontSize',14,'fontname','times','FontWeight','b');ylabel('Averaged per Block SINR at NB','FontSize',14,'fontname','times','FontWeight','b');fig_ix = fig_ix + 1;cfigure(fig_ix);handCDF = cdfplot(double(yy));hold on;set(handCDF,'color',colors{fileidx},'linestyle',lst{fileidx},'linewidth',2);result.figname{fig_ix}='CDF_of_bs_sinr';xlabel('SINR per Block at NB in dB','FontSize',14,'fontname','times','FontWeight','b');ylabel('CDF of SINR','FontSize',14,'fontname','times','FontWeight','b');% CDF of instantaneous cell throughput% calculated by B * R * * 6 * 12 (~Err) / T;function [result,fileidx] = res_cdf_instTP(test_resource_allocation_idx,UE,MIMO,loaddir,result,fileidx,av_factor,Comm_con f)filename = ['RA',num2str(test_resource_allocation_idx),'_UE',num2str(UE),'_FLA',num2str(Comm_conf.F_LA ),'_FPA',num2str(Comm_conf.F_PA),MIMO,'_trms',num2str(Comm_conf.t_rms),'f',num2str(Com m_conf.fd),'_s',num2str(Comm_conf.state),'N',num2str(Comm_conf.N),'subN',num2str(Comm_co nf.subN),'_E',num2str(Comm_conf.flg_csi_err),'c',num2str(Comm_conf.flg_offset),'D',num2str(Co mm_conf.DELAY),'BO',num2str(Comm_conf.BO_dB),'HPA',num2str(Comm_e_HPA_flg),'_B LER',num2str(Comm_conf.targetFER),'.mat'];load([loaddir,'/',filename],'save_bs_mod','save_bs_fec','save_bs_FrminErr');colors = {[0.5 0.5 0.5],[0 0 0],[0 0.5 1],[0 0 1],[1 0 0.5],[1 0 0]};lst = {'-','-.','-','-.','-','-.'};RAs = {'RR','PF','MX','PF2'};para.RA = RAs{test_resource_allocation_idx};UEs = [4 10];para.UE = UE;para.ra_ix = test_resource_allocation_idx;if Comm_conf.F_LA == 2 && Comm_conf.F_PA == 2,para.ra_ix = 4;endif MIMO == 2,para.ra_ix = 5;endpara.ue_ix = find(UEs == para.UE);para.nuSymbol = 6;para.subN = 12;para.Tfrm = 5e-4;save_bs_instTP = 1e-6 * sum(double(save_bs_mod) .* save_bs_fec .* double(~save_bs_FrminErr),1) * para.nuSymbol * para.subN / 5e-4;mns = sprintf('%3.2f',mean(save_bs_instTP));fileidx = fileidx + 1;%%result.ra_idx(fileidx) = test_resource_allocation_idx;result.LGD{fileidx} = [para.RA,',UE',num2str(para.UE),',F',num2str(Comm_conf.F_LA),'P',num2str(Comm_conf.F_PA),',', MIMO];result.LGDmn{fileidx} =[para.RA,',UE',num2str(para.UE),',F',num2str(Comm_conf.F_LA),'P',num2str(Comm_conf.F_PA),',', MIMO,',M:',mns];%%cfigure;handCDF = cdfplot(save_bs_instTP);hold on;set(handCDF,'color',colors{fileidx},'linestyle',lst{fileidx},'linewidth',2);result.figname='CDF_of_bs_instTP';xlabel('Instantaneous Cell Throughput in Mbps','FontSize',14,'fontname','times','FontWeight','b'); ylabel('CDF of Instantaneous Cell Throughput','FontSize',14,'fontname','times','FontWeight','b');。