ulaqierxuefufxt.m

天线方向图可以看出此算法的性能

%线阵　指向30度的-30dB的切比雪夫方向图
clc
clear all
close all
%………………………………基本参数………………………………………………………
c=1500;                                                 
f=9000; 
lamda=c/f;                                             %波长λ
M=32;                                                  %阵元数
d=lamda/2;                                             %阵元间距
Q1=-90;
Q2=90;                                                 %扫描范围
theta=Q1:0.1:Q2;                                       %方向角向量（１x1800）
k=length(theta);                                       %扫描点数（1800）
%…………………………指向30度………………………………………
Bd=pi*sin(30*pi/180);                                           %指向30度
for i=1:k
    B(i)=pi*(sin((Q1+(i-1)*0.1)*pi/180));        %空间频率β=(2πdsinθ)/λ
    g(i)=sin(M*(B(i)-Bd)/2)./(M*sin((B(i)-Bd)/2));      %阵列输出
end
G=abs(g);
dB=20*log10(G);
%－－－－－－－－－－－－绘制图形－－－－－－－－－－－－－
subplot(211)
plot(theta,dB);xlabel('方位角(度)');ylabel('G(\theta)(dB)')
axis([-50,50,-50,1]);
title('30度方向图')
grid on
%…………………………加-30dB切比雪夫权………………………………………
y2=zeros(1,k);
w=chebwin(M,30);                                 %w是-30dB的切比雪夫权
for i=1:M
   y1=exp(j*(i-1)*Bd)*exp(-j*(i-1)*pi*sin(theta*pi/180))*w(i);
   y2=y2+y1;
end
ymax=max(abs(y2));
yy=20*log10(abs(y2)/ymax); 
subplot(212)
 plot(theta,yy,'m')
    axis([-50,50,-50,0])
     xlabel('方位角(度)');
  ylabel('G(\theta)(dB)');
    title('加30dB切比雪夫权阵的30度方向图');