%雷達系統(tǒng)分析第二個大作業(yè)
%by?hutu
close?all;clear?all;clc
%%?雷達參數(shù)
load（‘IQ_data.mat‘）;%數(shù)據導入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;%脈沖數(shù)
AntenaNumber=20;%天線數(shù)
SampleNumber=3000;%取樣數(shù)
%%?第十一個天線接收的信號
antenna_data=permute（IQ_mid（:11:）[1?3?2]）;%將3000*1*12維變?yōu)?000*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））;%?匹配濾波器系數(shù)
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;%動雜波濾波器階數(shù)
coe_num=amti_rank-1;%需計算的權值個數(shù)＝濾波器階數(shù)-1；第一個權值為1，不需要計算
amti_num=3;%動雜波濾波器個數(shù)
%濾波器零點
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;%第一個常系數(shù)設為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:）））;%系數(shù)歸一化?
????????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（