clear?all
clc
close?all
%?n基元的數量
%?b波束的指向，單位為角度
%?fs抽樣頻率
%?c波速
%?f信號頻率
%?an為陣元點的位置

%以下為設置可以修改
N=1?????????????????????????%?N波束圖的角度分辨力
T=0.01??????????????????????%?T信號的長度
f=100;c=1500;b=-30;fs=5*f;an=[0?0;1?1;?-1?2;?4?5;6??7;-3?-4;-7?-5;-2?-5]???%?信號初始信息



%以下為程序
t=0:1/fs:T;?????????????????%?抽樣時間
s0=exp（j*2*pi*f*t）;?????????%?原始信號
number=length（t）;???????????%?采樣點的長度
[nv]=size（an）;?????????????%?n為陣元個數
S=zeros（nnumber）;??????????%?建立一個N行length（t）列的零矩陣用來表示收到的信號???????????????????????
P=zeros（nnumber）;??????????%?建立一個N行length（t）列的零矩陣用來表示處理后的信號
power=zeros（1（360/N）+1）;???%?用來放功率
for?m=1:n
????bn（m:）=an（m:）-an（1:）;
end
cn（:1）=-sin（b*pi/180）*bn（:1）;
cn（:2）=-cos（b*pi/180）*bn（:2）;
dn=cn（:1）+cn（:2）;
tao=dn/c;?????????????????????%?tao為各陣元點與第一個陣元點的時延
for?m=1:n
????kn=tao（m1）;
????S（m:）=exp（j*2*pi*f*（t+kn））;??