%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