波束形成MATLAB代码

常规的波束形成方法

clc;
clear;
close all;

fs = 1000;
c = 150;
N = 128;
f = 100;
lambda = c / f;
d = 0.5 * lambda;
theta = 1: 1: 180;
t = (0:1:1000-1) / fs;

A = zeros(1, length(theta));
A(5) = 5;
A(20) = 4;
A(25) = 5;
A(50) = 3;   % 在这四个方向上有目标
S = zeros(length(theta), length(t));
for theta_i = 1: 1: 180
    S(theta_i, :) = sin(2*pi*f.*t) * A(theta_i);
end
W = zeros(N, length(theta));
for i = 1: 1: 180
    for j = 1: N
        W(j, i) = exp(-2j*pi*f*(j-1)*d*sind(i)/c);  % 阵列流形
    end 
end
p = real(W * S);  % 各阵元的接收信号

figure(1)
for i = 1: N
    i
   plot(p(i, :));
   ylim([-20, 20]);
   pause(0.01)
end
figure(2)
plot(sum(p))

y = (W.' * p);  % 波束形成
y2 = real(y);
y3 = sum(y2, 2);

figure(3)
plot(sum(S, 2));
figure(4)
plot(y3)