相控阵雷达的波位编排..docxVIP

%本程序用于雷达球坐标系到正弦空间坐标系的转换,以及相控阵雷达的波位编排。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