移动通信蜂窝系统——频率复用及同信道干扰比

版权声明:本文为博主原创文章,遵循 CC 4.0 by-sa 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/zhzxlcc/article/details/70214518

不得不说,移动通信的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与距离的关系

展开阅读全文

没有更多推荐了,返回首页