MTD雷达信号处理

MTD雷达信号相参处理，有详尽的仿真程序和信号源。MTD雷达信号相参处理，有详尽的仿真程序和信号源MTD雷达信号相参处理，有详尽的仿真程序和信号源

%MTD滤波器组设计
%
%tic;
clear all;
clc;
N=32;%阶数可变
Tr=0.001;Fr=1/Tr;%可修改的脉冲重复周期，%数值由自己填写
f=0:2:Fr;  
delta2_clutter=16;  %杂波方差确定零陷凹口的大小
sigma2_clutter=10^6;  %杂波功率确定凹口深度
sigma2_noise=1;    % 噪声功率
%for fi=200:65:4035
Df=40;   %主瓣频率和副瓣衰减
interval=20;  %interval_former=20
f_interrupt=2:interval:Fr; %干扰频率    
L=length（f_interrupt）;  %干扰频率个数
%%%%%%%%%%%%%%%%%%%%%%基本参数设定%%%%%%%%%%%%%%%%%%%%%%
a_fi=a_f（NfiTr）; %导频矢量函数
%
for j=1:N          %杂波协方差矩阵Rc
    for i=j:N
     Rc（ij）=sigma2_clutter*exp（（-2）*（pi^2）...
         *delta2_clutter*（j-i）^2*Tr^2）;
     Rc（ji）=Rc（ij）;
    end
end
%%
for i=1:length（f）
    a_f_fd（::i）=a_f（Nf（i）Tr）;    %导频矢量组成的矩阵并且降维过后的函数
end
[a_f_fd_rowa_f_fd_colz]=size（a_f_fd）;
for i=1:a_f_fd_row
    for j=1:z
        a_f_fd2（ij）=a_f_fd（i1j）;
    end
end
%初始化干扰功率
for i=1:L
    sigma2_interrupt（1i）=0;
end
%%%%%%%%%%%%%%%%%以上不需改动%%%%%%%%%%%%%%%%%%%%%%%%%%%
k=0;
%%%%%%%%%%%%%%%%%修改参数部分%%%%%%%%%%%%%%%%%%%%%%%%%%%
K=0.08;                               %迭代增益0.1
K1=4;                                 %K1=3
K2=8;                                 %K2=5
K3=1.05;
K4=1.03;
for k=0:20                            %迭代次数20
Rc_inter=zeros（size（Rc））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%加权和
for j=1:L
a_f2=a_f（Nf_interrupt（j）Tr）*ctranspose（a_f（Nf_interrupt（j）Tr））;
Rc_interrupt（::j）=sigma2_interrupt（k+1j）.*a_f2;
Rc_inter=Rc_inter+Rc_interrupt（::j）;
end
%构造权矢量函数
Ra（::k+1）=Rc+sigma2_noise.*eye（size（Rc））+Rc_inter;   %协方差矩阵
Wi（::k+1）=inv（Ra（::k+1））*a_fi;%权矢量函数

%%%%%%%%%%频响应函数%%%%%%%%%%%%%%%%%
Hi（k+1:）=ctranspose（Wi（::k+1））*a_f_fd2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h（k+1）=abs（max（Hi（k+1:）））;
for i=1:L
    d（ik+1）=h（k+1）/（10^（Df/20））;
end
for i=1:L
    eta（i）=sigma2_interrupt（k+1i）+K*（abs（Hi（k+1（（i-1）*interval+2）/2+1））...
                                       -d（ik+1））;
    sigma2_interrupt（k+2i）=max（0eta（i））;
    if（i<=1|i>=L-1）
        sigma2_interrupt（k+2i）=0;
    end
  if（（i>（fi-K1*interval）/interval）&（i<（fi+K2*interval）/interval））  
      sigma2_interrupt（k+2i）=0;
  end
 if（i>=L-2）
       sigma2_interrupt（k+2i）=K3*sigma2_interrupt（k+2i）;
   end
  if（i>=8&i<=10）
     sigma2_interrupt（k+2i）=K4*sigma2_interrupt（k+2i）;  
   end
    %%%%%%%%%%
end
end
W（::1）=Wi（::k+1）;
plot（f20*log10（abs（Hi（k+1:）/max（abs（Hi（k+1:））））））;
% hold on;
% plot（fabs（Hi（k+1:）））;
% title（‘修改后第一个滤波器的增益函数‘）;xlabel（‘频率/Hz‘）;ylabel（‘增益值‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2678  2011-01-24 10:41  mtd.m

----------- ---------  ---------- -----  ----

                 2678                    1