均匀圆阵波束形成Matalb仿真

clc;
clear;
clear all;
c=1500;%声速
f0=1000;%信号频率 
r=4;%均匀圆阵的半径
t=-10:0.001:10;%时间
N=60;%阵元个数
target=40;%目标信号入射角度
theta=-180:180;
for count=1:length(theta)
    for n=1:N
        tao(n)=r*cosd(theta(count)-(n-1)*360/N)/c;%当信号入射方向改变为theta产生的延时
        tao0(n)=r*cosd(target-(n-1)*360/N)/c;%第n路信号添加的延时
        s(n,:)=cos(2*pi*f0*(t+tao(n)-tao0(n)));%第n路形成的信号
    end
    St=sum(s);%对所有阵元的信号进行空间求和
    E(count)=sum(St.^2);%系统的输出
end

E=E/max(E);%能量归一化
figure;
plot(theta,E);xlabel('角度');title('指向性图');
ylabel('指向性归一化值');
xlabel('角度/度');

xi

 相关参数：

阵元数N ：60 圆阵的半径r/m ：4 声速c/m/s ：1500  信号频率f0/Hz  ：1000   目标信号入射角target/° ：    40

当信号的频率升高，主瓣变窄，旁瓣增多，但旁瓣高度基本不变。