D2D-HS算法

HS算法解决D2D资源分配优化问题，基于和声搜索算法提出了一种集中式的联合频谱资源分配和功率控制（JRAPC）算法，针对所考虑的系统场景设计了相应的算法编码方案，通过对不可行解的处理很大程度上避免了算法陷入局部最优的问题。

function [SINR_cue_dBRate]=cal_sinr_rate_ee（alpha_d_newP_d_newPr_k_ue_mwC_ue_posD_ue_pos_txD_ue_pos_rxpower_noise）
 d_dx_rx = sqrt（ （D_ue_pos_tx（:1）-D_ue_pos_rx（:1））.^2 + （D_ue_pos_tx（:2）-D_ue_pos_rx（:2））.^2）; 
 Pr_drx = P_d_new - 148-40*log10（0.001*d_dx_rx）;
 Pr_drx_mw = power（10Pr_drx/10）;
 
 
 for j_ = 1:length（C_ue_pos）
     temp_matric = find（alpha_d_new == j_）;
    
     if isempty（temp_matric） ~= 1
         inter_ = 0;
         for i_ = 1:length（temp_matric）
             d_dtx_C_ue = sqrt（ （D_ue_pos_tx（temp_matric（i_）1）-C_ue_pos（j_1））^2 + （D_ue_pos_tx（temp_matric（i_）2）-C_ue_pos（j_2））^2）; 
             Pr_k_inter_ue =P_d_new（temp_matric（i_））- 148-40*log10（d_dtx_C_ue）;
             inter_ = inter_ + power（10Pr_k_inter_ue/10）;
         end
         SINR_cue（j_） = Pr_k_ue_mw（j_）/（inter_+power_noise）;
        
     else
         SINR_cue（j_） = Pr_k_ue_mw（j_）/power_noise;
     end
 end
 SINR_cue_dB = 10 * log10（SINR_cue）;
 %% 计算DUE SINR 
 P_K = 23;
 for k_ = 1:length（D_ue_pos_tx）
     alph_temp = alpha_d_new;
     alph_temp（k_） = [];
     temp_matr = find（alph_temp == alpha_d_new（k_））;
     d_cue_drx = sqrt（ （D_ue_pos_rx（ k_1）-C_ue_pos（alpha_d_new（k_）1））^2 + （D_ue_pos_rx（ k_2）-C_ue_pos（alpha_d_new（k_）2））^2）; 
     Pr_cue_inter = P_K-48-40*log10（d_cue_drx）;
     Pr_cue_inter = power（10Pr_cue_inter/10）;
     if isempty（temp_matr）~=1
         inter_1 = 0;
         for i_ = 1:length（temp_matr）
            d_rx_inter_tx = sqrt（ （D_ue_pos_tx（temp_matr（i_）1）-D_ue_pos_rx（k_1））^2 + （D_ue_pos_tx（temp_matr（i_）2）-D_ue_pos_rx（k_2））^2）; 
            Pr_d_inter = P_d_new（temp_matr（i_）） - 148-40*log10（d_rx_inter_tx）;
            inter_1 = inter_1 + power（10Pr_d_inter/10）;
         end
         SINR_d（k_） = Pr_drx_mw（k_）/（Pr_cue_inter+inter_1+power_noise）;
     else
         SINR_d（k_） = Pr_drx_mw（k_）/（Pr_cue_inter+power_noise）;
     end
     
 end
   Rate = log（1+SINR_d）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1968  2018-07-14 09:43  HS\cal_sinr_rate_ee.m

     文件       5392  2018-07-14 09:12  HS\main1.asv

     文件       5474  2018-07-14 09:54  HS\main1.m

     文件       5497  2018-07-14 09:27  HS\main2.asv

     文件       5782  2018-07-14 09:49  HS\main2.m

     文件        972  2018-07-12 20:03  HS\topology.m

     目录          0  2018-07-14 09:56  HS

----------- ---------  ---------- -----  ----

                25085                    7