DBF的matlab仿真程序

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: