%?===========================================================================================%
%????????????該程序完成16個脈沖信號的【脈壓、動目標顯示/動目標檢測（MTI/MTD）】????
%?===========================================================================================%
%??程序中根據每個學生學號的末尾三位（依次為XYZ）來決定仿真參數，034
%??目標距離為[3000?8025?9000+（Y*10+Z）*200?8025]，4個目標
%??目標速度為[50?0?（Y*10+X+Z）*6?100]
%?===========================================================================================%
close?all;?%關閉所有圖形
clear?all;?%清除所有變量
clc;
%?===================================================================================%
%????????????????????????????????????雷達參數???????????????????????????????????????%
%?===================================================================================%
C=3.0e8;??%光速（m/s）
RF=3.140e9/2;??%雷達射頻?1.57GHz
Lambda=C/RF;%雷達工作波長
PulseNumber=16;???%回波脈沖數
BandWidth=2.0e6;??%發射信號帶寬?帶寬B=1/τ，τ是脈沖寬度?
TimeWidth=42.0e-6;?%發射信號時寬
PRT=240e-6;???%?雷達發射脈沖重復周期（s）240us對應1/2*240*300=36000米最大無模糊距離
PRF=1/PRT;
Fs=2.0e6;??%采樣頻率
NoisePower=-12;%（dB）;%噪聲功率（目標為0dB）
%?---------------------------------------------------------------%
SampleNumber=fix（Fs*PRT）;%計算一個脈沖周期的采樣點數480；
TotalNumber=SampleNumber*PulseNumber;%總的采樣點數480*16=；
BlindNumber=fix（Fs*TimeWidth）;%計算一個脈沖周期的盲區-遮擋樣點數；
%===================================================================================%
%????????????????????????????????????目標參數???????????????????????????????????????%
%===================================================================================%
TargetNumber=4;%目標個數
SigPower（1:TargetNumber）=[1?1?1?0.25];%目標功率無量綱
TargetDistance（1:TargetNumber）=[3000?8025?15800?8025];%目標距離單位m???距離參數為[3000?8025?9000+（Y*10+Z）*200?8025]
DelayNumber（1:TargetNumber）=fix（Fs*2*TargetDistance（1:TargetNumber）/C）;%?把目標距離換算成采樣點（距離門）?fix函數向0靠攏取整
TargetVelocity（1:TargetNumber）=[50?0?204?100];%目標徑向速度?單位m/s???速度參數為[50?0?（Y*10+X+Z）*6?100]
TargetFd（1:TargetNumber）=2*TargetVelocity（1:TargetNumber）/Lambda;?%計算目標多卜勒頻移2v/λ
?
%====================================================================================%
%???????????????????????????????????產生線性調頻信號??????????????????????????????????%
%====================================================================================%
number=fix（Fs*TimeWidth）;%回波的采樣點數=脈壓系數長度=暫態點數目+1
if?rem（number2）~=0??%rem求余
???number=number+1;
end???%把number變為偶數
?
for?i=-fix（number/2）:fix（number/2）-1
???Chirp（i+fix（number/2）+1）=exp（j*（pi*（BandWidth/TimeWidth）*（i/Fs）^2））;%exp（j*fi）*，產生復數矩陣Chirp
end
coeff=conj（fliplr（Chirp））;%把Chirp矩陣翻轉并把復數共軛，產生脈壓系數
figure（1）;%脈壓系數的實部
plot（real（Chirp））;axis（[0?90?-1.5?1.5]）;title（‘脈壓系數實部‘）;
%-------------------------產生目標回波串-----------------------------------------------------------------------------------------%
%-------------------------產生前3個目標的回波串-------%
SignalAll=zeros（1TotalNumber）;%所有脈沖的信號先填0
for?k=1:TargetNumber-1?%?依次產生各個目標
???SignalTemp=zeros（1SampleNumber）;%?一個PRT
???SignalTemp（DelayNumber（k）+1:DelayNumber（k）+numbe