fengwo.txt

matlab画的峰窝图

% beehive
clc,close all
r=1;
rc=0.7;
figure;
axis square
hold on;
A=pi/3*[0:6];
aa=linspace(0,pi*2,80);
plot(r*exp(i*A),'k','linewidth',2);
% Author email: [e-mail]
% From: 萝卜驿站
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];
for k=1:3;
    Z=Z+sqrt(3)*r*exp(i*pi/6);
    for pp=1:6;
        for p=1:k;
            N=N+1;
            zp=Z+r*exp(i*A);
            zr=Z+rc*exp(i*aa);
            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));
        end
    end
end
xlim([-6,6]*r)
ylim([-6.1,6.1]*r)
axis off;