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

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



%以下為程序
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（180/N）+1）;???%?用來放功率
coef=sqrt（10^（SNR/10））;?
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