clear all
close all
num=16;        %% 陣元數
k=1000;          %% 快拍數
as=0;            %波束指向角
C=340;
fc=3400;
f0=1700; 
D=(C/f0)/2;
 
e_position=[0:num-1]';
aa=100*[-1:0.001:1];                                            %% scan vector
vs=exp(j*2*pi*fc*e_position*sin(as/180*pi)/C);
v=exp(j*pi*e_position*sin(aa/180*pi));  %% 
w=1/num*vs';

figure(1),plot(aa,20*log10(abs(w*v)));
axis([-90 90 -50 0]);
hold on