相控阵雷达的波位编排..docx
%本程序用于雷达球坐标系到正弦空间坐标系的转换,以及相控阵雷达的波位编排。070129,bygoodtony.
%070228添加波位数据保存
clearall
closeall
clc
%波束宽度,单位：度
BeamWidth=3.6;
theta05=round(sin(BeamWidth*pi/180)*1000)/1000;
%雷达作用空域
AZ=[-1515];%方位角范围(雷达球坐标系)，单位：度
EL=[030];%俯仰角范围(雷达球坐标系)，单位：度
%阵面倾角，单位：度
thetaT=atan(-(cos(EL(1)*pi/180)-cos(EL(2)*pi/180))/(sin(EL(1)*pi/180)-sin(EL(2)*pi/180))*cos(AZ(2)*pi/180))*180/pi;
%左边界
az=AZ(1);
el=EL(1):EL(2);
alpha1=cos(el*pi/180)*sin(az*pi/180);
beta1=sin(el*pi/180)*cos(thetaT*pi/180)-cos(el*pi/180)*cos(az*pi/180)*sin(thetaT*pi/180);
%右边界
az=AZ(2);
el=EL(1):EL(2);
alpha2=cos(el*pi/180)*sin(az*pi/180);
beta2=sin(el*pi/180)*cos(thetaT*pi/180)-cos(el*pi/180)*cos(az*pi/180)*sin(thetaT*pi/180);
%下边界
az=AZ(1):AZ(2);
el=EL(1);
alpha3=cos(el*pi/180)*sin(az*pi/180);
beta3=sin(el*pi/180)*cos(thetaT*pi/180)-cos(el*pi/180)*cos(az*pi/180)*sin(thetaT*pi/180);
%上边界
az=AZ(1):AZ(2);
el=EL(2);
alpha4=cos(el*pi/180)*sin(az*pi/180);
beta4=sin(el*pi/180)*cos(thetaT*pi/180)-cos(el*pi/180)*cos(az*pi/180)*sin(thetaT*pi/180);
lowup=max(beta3);
lowest=min(beta3);
rightest=max(alpha3)+theta05/2;
leftest=min(alpha3)-theta05/2;
upest=max(beta4);
uplow=min(beta4);
figure(1),plot(alpha1,beta1,r);axis([leftest-0.1rightest+0.05lowest-0.05upest+0.1]);holdon;gridon;
title(雷达波束在正弦空间的编排);xlabel(alpha(rad));ylabel(beta(rad));
plot(alpha2,beta2,r);
plot(alpha3,beta3,r);
plot(alpha4,beta4,r);
%pause(0.1)
%==========================================================================
BoweiNumber=0;
Mx=2*round(rightest/theta05)-1;
My=round((upest-lowest)/(0.866*theta05));
arr_x0=zeros(My,Mx);%波位坐标（方位角）
arr_y0=zeros(My,Mx);%波位坐标（俯仰角）
arr_NboweiRow=zeros(My,1);%每一行的波位数目
%---------------奇数行------------------------
m=0;
fory0=lowest+theta05/2:2*0.866*theta05:upest
m=m+1;
n=0;
NboweiRow=0;
forx0=0:theta05:rightest