频域宽带波束形成（fft和LFMsource）代码

%fft
xx = zeros(J,NN,M); % x的第三维是M个阵元，第一维是要做的FFT的点数，第二维是要做的J点FFT的次数
for n=1:NN
for m=1:M
xx(:,n,m) = x(m,(n-1)*J+1:(n-1)*J+J);
end
end
X = fft(xx,J,1); % FFT变换完成 第一维是子频点的个数
% figure
% plot((0:fs/J:fs-fs/J),abs(X(1:J,1,1)))

%LFM
function [x] = LFMsource(M,f0,fl,B,fs,T,seta,snr,K)
 
c = 310^8;                       % 光速
d = 0.5c/f0;                     % 阵元间距
ft = 0:1/fs:T-1/fs;               % 时间变量
N = length(ft);                   % 快拍数
k = B/T;                          % 调频速率
for m=1:M
    for n=1:N
        x(m,n) = 10^{(snr(K)/20)*exp(i*(2*pi*fl*(ft(n)-(m-1)*d/c*sin(seta(K)))+pi*k*(ft(n)-(m-1)*d/c*sin(seta(K))).}2));
    end
end
end