clc;clear?all;close?all;
%%%%%%%%%%%%%%%%%一維DBF%%%%%%%%%%%%%%%%%
K=8;%陣元個數
wavelength=0.1;%波長
d=wavelength/2;%陣元間距
theta0=（-60:60）*pi/180;%波達方向
NFFT=K;%FFT點數
W=chebwin（K40）;%切比雪夫窗
St=zeros（Klength（theta0））;
delta_phase=pi/K;
S=exp（j*2*pi*（0:K-1）‘*（d*sin（theta0）/wavelength-delta_phase/pi/2））;%陣列空域導向矢量
for?ii=1:length（theta0）
????St（:ii）=W.*S（:ii）;
end
B=fftshift（fft（StNFFT1）1）;
for?jj=1:K
????Bn=abs（B（jj:））/max（abs（B（jj:）））;
????plot（theta0*180/pi20*log10（Bn）‘LineWidth‘2）;hold?on;
end
xlabel（‘方位/度‘）;ylabel（‘幅度/dB‘）;
title（‘數字波束形成‘）;axis（[min（theta0）*180/pimax（theta0）*180/pi-50?0]）;

%%%%%%%%%%%%%%%%%二維DBF%%%%%%%%%%%%%%%%%
%?M=8;%陣元行數
%?N=4;%陣元列數
%?wavelength=0.1;%波長
%?d=wavelength/2;%陣元間距
%?theta=（-90:90）*pi/180;%波達方向
%?fai=（-90:90）*pi/180;%波達方向
%?NFFT1=M;%FFT點數
%?NFFT2=N;%FFT點數
%?W1=chebwin（M30）;%切比雪夫窗
%?W2=chebwin（N30）;%切比雪夫窗
%?W=W1*W2.‘;
%?[theta0fai0]=meshgrid（thetafai）;
%?B=zeros（length（theta）length（fai））;
%?for?xx=1: