移動通信蜂窩系統——頻率複用及同信道干擾比

不得不說，移動通信的MATLAB實驗越來越難了，這才第二次試驗，表示上次實驗課上什麼都沒幹，感覺其他同學有做的挺快的，然而自己就畫出了一層六邊形蜂窩，蜂窩不會，複用不會，更別提CCIR了，根本沒時間研究，實在是失落至極╮(╯▽╰)╭

明天上午又是實驗課，新的實驗任務已經出來了，沒辦法，寫了半個上午加半個下午的MATLAB實驗，最終算是大體完成了，過程也是一波三折，尤其頻率複用的標號，模7運算竟然出現了7，見鬼Q_Q弄過來弄去，午休的時候也在想，想了半天，可能因爲六邊形裏的sqrt(3)，然後在計算座標的時候就出現了各種問題，開始是出現接近於0的小數，我就嘗試取整，後來發現一個cluster裏有同頻的，最後就是上面說的模7運算竟然出現了7，真是懵掉了。。最後用fix向0取整，而且在（u，v）座標的時候取整，而不是給標號L取整，終於搞定了，蜂窩圖大體正確沒花很長時間，但是細節上出錯讓我極其地心煩，因爲一直改不對。。

一個cluster裏面有同頻信道Q_Q

mod7出現了7(⊙o⊙)
代碼雖然不是很長，但需要點數學和邏輯理解，迭代加平移變換之類的。附上改正確之後的代碼：

        R = 1;dy = sqrt(3)*R;dx = 1.5*R;N = 7;
        A = pi/3*[1:7];
        yp =[0:dy:100,0:-dy:-100];
        xp = [0:dx:100,0:-dx:-100];
        for t = 1:1:length(yp)
        for m = 1:1:length(xp)
            y = sqrt(3)/3*xp(m)+ yp(t);
            if  -15<xp(m) && xp(m)<15 && -10<y && y<10
            plot(xp(m)+1i*y+R*exp(1i*A),'linewidth',1);
            L = mod(fix(2*xp(m)/R) + fix(yp(t)/(sqrt(3)*R)),N);
            text(xp(m),y,num2str(L),'FontSize',10);
            hold on
            end
        end
        end
        grid on 
        axis equal 
        title('Cell label for 7-cell cluster')

N=7的蜂窩圖
然後是同信道干擾比，由上可發現，N=3 和N=7時第一層的同頻小區數均爲6（N=3時上述程序稍作修改即可），所以，每個座標的夾角爲60°。首先，需要確定第一個同頻小區中心座標的位置，然後再依次旋轉60°即可得到6個同頻小區的座標，然後計算同頻小區到MS的距離，再根據公式即可得到CCIR，注意MS到服務基站的距離爲0時，需要將距離重新賦值爲一個較小的數，否則CCIR爲無窮大。主要代碼如下所示：

pxy = 0:0.2:R;
for k = 1:length(pxy)
    if pxy(k) <exp(-10)
        pxy(k) = 0.1;
    end
end
co_cell(1) = sqrt(3)*R*m1*exp(i*pi/6) + sqrt(3)*R*exp(i*pi/2);
for t = 2:6
    co_cell(t) = co_cell(t - 1)*exp(i*pi/3);
end
sum = 0;
for u = 1:6
    sum = abs(co_cell(u) - pxy).^(-alpha) + sum;
end
gamma = pxy.^(-alpha)./sum;

生成的圖像如下：

CCIR與距離的關係