雷达系统仿真大作业

%雷达系统分析第二个大作业
%by hutu
close all;clear all;clc
%% 雷达参数
load（‘IQ_data.mat‘）;%数据导入3000*20*12
C=3e8;%光速
Lambda=1.25;
Pfa=10e-6;%
Ts=1e-6;%采样周期
BandWidth=8e5;%带宽800KHz
taoup=420e-6;%时宽420us
mu=BandWidth/taoup;%μ
T1=4114e-6;%三变T，时间
T2=4400e-6;
T3=4686e-6;
d=0.65;%天线间距
PulseNumber=12;%脉冲数
AntenaNumber=20;%天线数
SampleNumber=3000;%取样数
%% 第十一个天线接收的信号
antenna_data=permute（IQ_mid（:11:）[1 3 2]）;%将3000*1*12维变为3000*12*1维，方便画图
figure（1）;plot（20*log10（abs（antenna_data）））;
title（‘第十一根天线接收的信号‘）;
xlabel（‘距离单元‘）;ylabel（‘信号幅度/dB‘）;
axis tight;zoom on;grid on;

%% DBF
theta0=0;%天线指向
a0=exp（j*2*pi*d/Lambda*（0:AntenaNumber-1）‘*sin（theta0））;%导向矢量
win=taylorwin（AntenaNumber5-25）;%-25dB 泰勒窗
ww=a0.*win;
IQ_DBF=zeros（SampleNumberPulseNumber）;
for i=1:SampleNumber
    for k=1:PulseNumber
        IQ_DBF（ik）=IQ_mid（i:k）*ww;
    end
end
figure（2）;plot（20*log10（abs（IQ_DBF）））;
xlabel（‘距离单元‘）;ylabel（‘信号幅度/dB‘）;title（‘DBF‘）;
axis（[0 3000 60 170]）;zoom on;grid on;
%% 脉压
number=fix（taoup/Ts）;
t=linspace（-0.5*taoup0.5*taoupnumber）;
tay=taylorwin（number9-35）;%-35dB泰勒窗
ht=exp（sqrt（-1）*pi*mu*（t.^2））;% 匹配滤波器系数
ht=ht.*tay.‘;%加窗
ht_fft=fft（ht4096）;
buff=zeros（14096）;
IQ_ifft=zeros（4096PulseNumber）;
for i=1:PulseNumber
    buff=ifft（fft（IQ_DBF（:i）4096）.*ht_fft.‘）;
    IQ_ifft（:i）=buff;
end
figure（3）;
plot（20*log10（abs（IQ_ifft（::））））;
xlabel（‘距离单元‘）;ylabel（‘信号幅度/dB‘）;title（‘脉压‘）;
axis（[0 3000 80 180]）;zoom on;grid on;
%% MTI
bianT=[4.114e-3  4.4e-3  4.686e-3];%变T
amti_rank=6;%动杂波滤波器阶数
coe_num=amti_rank-1;%需计算的权值个数＝滤波器阶数-1；第一个权值为1，不需要计算
amti_num=3;%动杂波滤波器个数
%滤波器零点
fc（1）=0-1;
fc（2）=0-0.5;
fc（3）=0;
fc（4）=0+0.5;
fc（5）=0+1;
f=[-100:0.1:100 101:1:2000];%频率范围
    TT=repmat（bianT1amti_num）;
    %产生A矩阵 计算权值
    W（1）=1;%第一个常系数设为1
    U=ones（coe_num1）;
    for n=1:amti_num;
        T1=TT（n:n+coe_num-1）;
        for k=1:coe_num Ti（k）=sum（T1（1:k））; 
        end;
        for p=1:coe_num        
            A1（p1:coe_num）=exp（-1*i*2*pi*fc（p）.*Ti（1:coe_num））;            
        end;    
        W（2:amti_rank）=-W（1）*inv（A1）*U;
        WW（n:）=W;
        WW（n:）=WW（n:）/max（abs（WW（n:）））;%系数归一化 
        hd（n:）=WW（n1）+WW（n2）*exp（-j*2*pi*Ti（1）*f）+WW（n3）*exp（-j*2*pi*Ti（2）*f）+WW（n4）*exp（-j*2*pi*Ti（3）*f）+WW（n5）*exp（-j*2*pi*Ti（4）*f）++WW（n6）*exp（-j*2*pi*Ti（5）*f）;     
    end
hd（n+1:）=mean（abs（hd（1:n:）））; 
Hd=10*log10（abs（hd））;
figure（4）;
plot（fHd）;xlim（[-100 2000]）;ylim（[-100 20]）;grid on;
title（‘MTI地杂波滤波器特性‘）;xlabel（‘Frequency/Hz‘）;ylabel（‘H（f）/dB‘）;legend（‘T1:T2:T3‘‘T2:T3:T1‘‘T3:T2:T1‘‘mean‘）;    
IQ_mti=zeros（SampleNumberPulseNumber-amti_rank+1）;
for i=1:PulseNumber-amti_rank+1
    n=mod（（i-1）3）+1;
    IQ_mti（:i）=IQ_ifft（1:SampleNumberi+（1:amti_rank）-1）*WW（n:）.‘;
end
figure（5）;
plot（20*log10（abs（IQ_mti）））;
title（‘MTI后原始视频‘）;xlabel（‘距离单元‘）;ylabel（