画波束程序

clear all;clc;
n=7;% 7 satellite beam
i=1;
beam_radius=36000*1e3*tan(0.15*pi/180); %satellite beam width

% method 1
t=linspace(0,2*pi,n);
x=cos(t+pi/6)*sqrt(3)*beam_radius;
y=sin(t+pi/6)*sqrt(3)*beam_radius;

 
% method 2 contains center beam point in [0 0]
beam0_pos=[0,0];
beam1_pos=[0,sqrt(3)*beam_radius];

c=cosd(60);
s=sind(60);

temp=zeros(6,2);
temp(1,:)=beam1_pos;          
for i=2:6
    temp(i,1)=temp(i-1,1)*c-temp(i-1,2)*s;
    temp(i,2)=temp(i-1,1)*s+temp(i-1,2)*c;
end

beam_pos(1,:)=beam0_pos;
beam_pos(2:7,:)=temp;
c=temp(:,1);
d=temp(:,2);


r = beam_radius;
%rectangle('Position',[x-r,y-r,2*r,2*r],'Curvature',[1,1],'EdgeColor','m');
plot([0,x],[0,y],'o','Markersize',200) 
xlim(-2*e5, 2*e5);
axis equal;