interference.m

来自「这是一个基于matlab/RUNE环境的移动通信网络仿真程序包」· M 代码 · 共 100 行

function  [CIR] = interference(complex_dist,logn4,logn4a,logn4b,logn4c,logn4d,logn4e,logn4f,mode) 


%positions of first tier interference base stations
a= 0 - 1500i;
b= 1300 - 750i;
c= 1300 + 750i;
d= 0 + 1500i;
e= -1300 + 750i;
f= -1300 - 750i;

[m,n]=size(complex_dist);
CIR=zeros(m,n);
temp=zeros(m,n);


%logn4   is the location-based correlated shadow fading for donwlink channel from BS to MS
%gain antenna = 0 dB (omnidirectional)
%transmitted power = 1 watt = 30 dBm

%logn4a  is the location-based correlated shadow fading for interference from interfere-BS1 to MS
%logn4b  is the location-based correlated shadow fading for interference from interfere-BS2 to MS
%logn4c  is the location-based correlated shadow fading for interference from interfere-BS3 to MS
%logn4d  is the location-based correlated shadow fading for interference from interfere-BS4 to MS
%logn4e  is the location-based correlated shadow fading for interference from interfere-BS5 to MS
%logn4f  is the location-based correlated shadow fading for interference from interfere-BS6 to MS




% distance dependence interference
if (mode==1)
    
for i=1:m
    
    pathloss(i,:)=32.4 - 40.*log10(abs(complex_dist(i,:))./1000) - 20.*log(900);
    power_rcvd(i,:)=1.*10.^(pathloss(i,:)./10).*(10.^(logn4(i,:)./10));
    
  
    pathloss(i,:)=32.4 - 40.*log10(abs((complex_dist(i,:)-a))./1000) - 20.*log(900);
    interference(i,:)=1.*10.^(pathloss(i,:)./10).*(10.^(logn4a(i,:)./10));
    
    pathloss(i,:)=32.4 - 40.*log10(abs((complex_dist(i,:)-b))./1000) - 20.*log(900);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10).*(10.^(logn4b(i,:)./10));
    
    pathloss(i,:)=32.4 - 40.*log10(abs((complex_dist(i,:)-c))./1000) - 20.*log(900);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10).*(10.^(logn4c(i,:)./10));
    
    pathloss(i,:)=32.4 - 40.*log10(abs((complex_dist(i,:)-d))./1000) - 20.*log(900);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10).*(10.^(logn4d(i,:)./10));
    
     pathloss(i,:)=32.4 - 40.*log10(abs((complex_dist(i,:)-e))./1000) - 20.*log(900);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10).*(10.^(logn4e(i,:)./10));
    
    pathloss(i,:)=32.4 - 40.*log10(abs((complex_dist(i,:)-f))./1000) - 20.*log(900);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10).*(10.^(logn4f(i,:)./10));
    
    CIR(i,:)=power_rcvd(i,:)./interference(i,:);
    CIR(i,:)=10.*log10(CIR(i,:));
        
end
end



%COST231 ikegami

if (mode==2)
    for i=1:m
        
    pathloss(i,:)=cost231wi(abs(complex_dist(i,:))./1000);
    power_rcvd(i,:)=1.*10.^(pathloss(i,:)./10);%.*(10.^(logn4(i,:)./10));
        
    pathloss(i,:)=cost231wi(abs((complex_dist(i,:)-a))./1000);
    interference(i,:)=1.*10.^(pathloss(i,:)./10);%.*(10.^(logn4a(i,:)./10));
    
    pathloss(i,:)=cost231wi(abs((complex_dist(i,:)-b))./1000);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10);%.*(10.^(logn4b(i,:)./10));
    
    pathloss(i,:)=cost231wi(abs((complex_dist(i,:)-c))./1000);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10);%.*(10.^(logn4c(i,:)./10));
    
    pathloss(i,:)=cost231wi(abs((complex_dist(i,:)-d))./1000);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10);%.*(10.^(logn4d(i,:)./10));
    
    pathloss(i,:)=cost231wi(abs((complex_dist(i,:)-e))./1000);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10);%.*(10.^(logn4e(i,:)./10));
    
    pathloss(i,:)=cost231wi(abs((complex_dist(i,:)-f))./1000);
    interference(i,:)=interference(i,:) + 1.*10.^(pathloss(i,:)./10);%.*(10.^(logn4f(i,:)./10));
    
    CIR(i,:)=power_rcvd(i,:)./interference(i,:);
    CIR(i,:)=10.*log10(CIR(i,:));
    
    
    end
end