% beehive
clc,close all
r=1;
rc=2/3;
K=20;
alpha=0:2*pi/K:2*pi;
figure;
axis square
hold on;
A=pi/3*[0:6];
% aa=linspace(0,pi*2,80);
plot(r*exp(i*A),'k','linewidth',2);
plot(rc*exp(i*alpha),'ro','linewidth',2);
plot(0,0,'k+','linewidth',2);
% Author email: zjliu2001@163.com
% g1=fill(real(r*exp(i*A)),imag(r*exp(i*A)),'k');
% set(g1,'FaceColor',[1,0.5,0])
% g2=fill(real(rc*exp(i*aa)),imag(rc*exp(i*aa)),'k');
% set(g2,'FaceColor',[1,0.5,0],'edgecolor',[1,0.5,0],'EraseMode','xor')
text(0,0,'1','fontsize',10);
Z=0;
At=pi/6;
RA=-pi/2;
N=1;At=-pi/2-pi/3*[0:6];
k=1;
Z=Z+sqrt(3)*r*exp(i*pi/6);
plot(Z,'k+','linewidth',2);
for pp=1:6;
for p=1:k;
N=N+1;
zp=Z+r*exp(i*A);
zr=Z+rc*exp(i*alpha);
plot(zp,'k','linewidth',2);
if pp==1
plot(zr,'bo','linewidth',2);
elseif pp==2
plot(zr,'ko','linewidth',2);
elseif pp==3
plot(zr,'co','linewidth',2);
elseif pp==4
plot(zr,'go','linewidth',2);
elseif pp==5
plot(zr,'yo','linewidth',2);
else
plot(zr,'mo','linewidth',2);
end
% g1=fill(real(zp),imag(zp),'k');
% set(g1,'FaceColor',[1,0.5,0],'edgecolor',[1,0,0]);
% g2=fill(real(zr),imag(zr),'k');
% set(g2,'FaceColor',[1,0.5,0],'edgecolor',[1,0.5,0],'EraseMode','xor');
text(real(Z),imag(Z),num2str(N),'fontsize',10);
Z=Z+sqrt(3)*r*exp(i*At(pp));
plot(Z,'k+','linewidth',2);
end
end
% xlim([-6,6]*r)
% ylim([-6.1,6.1]*r)
legend('Cell','User','Base');
axis on;