频域宽带波束形成.m

% 宽带信号频域处理用一致聚焦算法进行波束形成
% 阵元数为8，间距为半波长，信号中心频率为3.1G，带宽50M
% 干扰为窄带频率为3.1G
%RSS算法
% 2008.3.10



%%%添加注释研究
clear all;
clc;
f0 = 3.1*10^9;                  % 信号中心频率
fr = 3.1*10^9;                 % 干扰信号频率
B = 300*10^6;                    % 信号带宽
fl = f0-B/2;                    % 信号起始频率
fh = f0+B/2;                    %信号截止频率
Tr = 50*10^（-6）;                % 工作周期，脉冲重复周期？？
T = 10*10^（-6）;                 % 信号持续时间，要求大于孔径渡越时间
                                %T和Tr的关系
fs=400*10^6;                    %采样频率？
snr = [40;0;40;40];                   % 信噪比？？？
theta0 = 0/180*pi;              % 指向俯仰角，方位角为0
% phi0 = 0/180*pi;                % 指向方位角
M = 16;                          % 行阵元数为M 
% Q = 8;                          % 列阵元数为M 
Ks = 2;                         % 信号数目
Theta =[35/180*pi0/180*pi-40/180*pi60/180*pi];     % 干扰信号方向
% Phi =[0/180*pi0/180*pi];     % 干扰信号方向
c=3*10^8;
d=0.5*c/f0;
N=T*fs;                       %一个周期采样点数，4000
Nr=Tr*fs;                     %%五个周期采样点数
Nm = 5;                        % 驻留周期个数
J = 200;                       % FFT的点数？？？FFT点数在频域宽带波束形成中起到什么作用？？
                               %当信号带宽较小时，相对于固定采样率，则信号所占的频率点数越少，则无法较好的左频带划分。
                               %将fft点数增多，则对同样带宽信号具有更详细的表示，同时更容易做频带划分。
fw=linspace（0fsJ+1）;              %前nfft点是fft后的频域频率
                                    %%%即将信号对应的频率点数，实际中这些点值有少部分点有值，对应带通信号
                                    %%%还是变换到基带进行处理，所以其频带
                                %%应当好奇的是这些点所对应的频带范围始终在一个频率周期以内即0—fs之内
                                    
fwn=[];                                %fft后落在带宽内频率
kn=[];
m=fix（f0/fs）;                  %%此处m=7;
for k=1:J                               %%%这里到底是用的带通采样，根据带通采样定理
    if fw（k）>=f0-m*fs-B/2 && fw（k）<=f0-m*fs+B/2  %%判断带通信号在一个频率周期内的数据是那哪些点，信号所占频带点数与信号带宽跟采样频率之比相等。
        fwn=[fwn fw（k）];              %fft后落在带宽内频点，保存这些数据，fwn对应频带的横轴
        kn=[kn k];                    %%kn为索引
    end
end
%%%%%对信号频带做子带划分和处理应仅对这部分完全的信号频带点做处理

G=length（kn）;                         % 信号频率带宽所对应的点数，频带划分，划分了多少个子带
NN = Nm*Nr/J;                   % 频域快拍数500个，
tic                             %%%%倒数计时器
% 产生阵列接收的宽带数据
x=[];
for nm=1:Nm
%     phase_rand=randn（Ks1）;
    x2 = LFMsource（Mf0flBfsTThetasnr2）; % 期望信号的脉冲宽度为T，出现时间不定，跟噪声功率相同，方向为0
    % 产生干扰信号
    x1 = LFMsource（Mf0flBfsTrThetasnr1）;%干扰信号1，整个重复周期的数据，功率大，方向为40
%     ft = 0:1/fs:Tr-1/fs;    
%     N1 = length（ft）;                   % 快拍数
%                         % 调频速率
%     for m=1:M
%     for n=1:N1
%         x3（mn） = 10^（snr（3）/20）*exp（i*（2*pi*fl*（ft（n）-（m-1）*d/c*sin（Theta（3）））））;
%                             %%%%fl是否该改成f0
%     end
%    end
%                          % 调频速率

    x3 = LFMsource（Mf0flBfsTrThetasnr3）;
    x4 = LFMsource（Mf0flBfsTrThetasnr4）;
%     nn=1:Nr;
%     for n=1:Nr
%         fai（1n） = rand;
%     end
%     s（1:） = exp（i*2*pi*（fr*nn/fs+k*fai））;
%     for k=1:Kr
%         for m=1:M
%