function?y=maiya（PulseNumberrangeFsTargetVelocitySNR）

C=3e8;
Tr=800e-6;%單脈沖最大可測(cè)距離0.5*Tr*C=120000m
Tp=160e-6;
fr=1/Tr;
lamda=0.03;
B=1e6;
%%?產(chǎn)生線性調(diào)頻信號(hào)?
number=fix（Fs*Tp）;%單脈沖發(fā)射信號(hào)內(nèi)采樣點(diǎn)數(shù)
SampleNumber=fix（Fs*Tr）;%一個(gè)脈沖采樣點(diǎn)數(shù)800；則每個(gè)采樣間距代表120000/799m
TotalNumber=SampleNumber*PulseNumber;
for?i=-fix（number/2）:fix（number/2）-1
????Chirp（i+fix（number/2）+1）=exp（j*（pi*（B/Tp）*（i/（Fs））^2））;?
end??
figure
plot（real（Chirp））;
title（‘線性調(diào)頻信號(hào)‘）
coeff=conj（fliplr（Chirp））;%匹配濾波函數(shù)h
figure
plot（real（coeff））
title（‘匹配濾波函數(shù)h‘）


W?=?linspace（-pipi4096）;????
[h1w0]?=?freqz（coeff1W）;??
figure
plot（w0/pi20*log10（abs（h1/max（h1））））;?
title（‘匹配濾波函數(shù)H頻率響應(yīng)‘）
%-------
%%?產(chǎn)生目標(biāo)回波
TargetDistance（1:2）=range;
DelayNumber（1:2）=fix（Fs*2*TargetDistance（1:2）/C）;

TargetFd（1:2）=2*TargetVelocity（1:2）/lamda;?%計(jì)算目標(biāo)多卜勒

SignalAll=zeros（1TotalNumber）;%所有脈沖的信號(hào)先填0
for?k=1:2%?依次產(chǎn)生各個(gè)目標(biāo)1?2?
????SignalTemp=zeros（1SampleNumber）;%?一個(gè)脈沖
???
????SignalTemp（DelayNumber（k）+1:DelayNumber（k）+number）=（10^（SNR（k）/10））*Chirp;??%一個(gè)脈沖的1個(gè)目標(biāo)（未加多普勒速度）
????Signal=zeros（1TotalNumber）;
????for?i=1:PulseNumber
????????Signal（（i-1）*SampleNumber+1:i*SampleNumber）=SignalTemp;
????end
????FreqMove=exp（j*2*pi*TargetFd（k）*（0:TotalNumber-1）/Fs）;%目標(biāo)的多普勒速度*時(shí)間=目標(biāo)的多普勒相移
????Signal=Signal.*FreqMove;
????SignalAll=SignalAll+Signal;
end
figure
subplot（211）;plot（real（SignalAll）‘r-‘）;title（‘目標(biāo)信號(hào)的實(shí)部‘）;grid?on;zoom?on;
subplot（212）;plot（imag（SignalAll））;title（‘目標(biāo)信號(hào)的虛部‘）;
grid?on;zoom?on;
%%?產(chǎn)生噪聲
SystemNoise=normrnd（010^（-12/10）1TotalNumber）+j*normrnd（010^（-12/10）1TotalNumber）;
%%?回波+噪聲
Echo=SignalAll+SystemNoise;???%?+SeaClutter+TerraClutter;
for?i=1:PulseNumber???%在接收機(jī)閉鎖期接收的回波為0
????Echo（（i-1）*SampleNumber+1:（i-1）*SampleNumber+number）=0;
end
figure
plot（real（Echo））;title（‘總回波信號(hào)的實(shí)部閉鎖期為0‘）;
%subplot（212）;plot（imag（Echo））;title（‘總回波信號(hào)的虛部閉鎖期為0‘）;

%%?回波積累
Echoplus=zeros（1SampleNumber）;
for?i=1:PulseNumber
???Echoplus1（1:SampleNumber）=Echo（（i-1）*SampleNumber+1:i*SampleNumber）;
???Echoplus=Echoplus+Echoplus1（1:SampleNumber）;
end
Echoplus=Echoplus/PulseNumber;
figure
plot（real（Echoplus））;title（‘回波脈沖積累信號(hào)的實(shí)部閉鎖期為0‘）;


%%?時(shí)域脈壓
pc_time0=conv（Echocoeff）;
%figure
%plot（abs（pc_time0））;title（‘時(shí)域脈壓結(jié)果的幅度有暫態(tài)點(diǎn)‘）;
%%?頻域脈壓
Echo_fft=fft（EchoTotalNumber+number-1）;%進(jìn)行TotalNumber+number-1點(diǎn)FFT
coeff_fft=fft（coeffTotalNumber+number-1）;
pc_fft=Echo_fft.*coeff_fft;
pc_freq0=ifft（pc_fft）;
pc_freq_abs=abs（pc_freq0）;
figure
plot（pc_freq_abs）;title（‘（不加窗）頻域脈壓結(jié)果的幅度（有暫態(tài)點(diǎn)）‘）;

%%??數(shù)據(jù)重排到距離門(mén)
pc_freq1=pc_freq0（number:TotalNumber+number-1）;
for?i=1:PulseNumber
???pc（i1:SampleNumber）=pc_freq1（（i-1）*SampleNumber+1:i*SampleNumber）;
end


figure

%subplot（211）
%plot（（0:SampleNumber-1）*（120000/（SampleNumber-1））?abs（pc（5:））?）;
%?title（‘不加窗脈壓結(jié)果‘）;xlabel（‘距離‘）ylabel（‘脈壓輸出‘）;
?
plot（（0:SampleNumber-1）*（120000/（SampleNumber-1））?20*log10（abs（pc（1:）））?）;
xlabel（‘距離‘）ylabel（‘脈壓輸出/db‘）;
axis（[-infinf?-1050]）;
title（‘不加窗脈壓結(jié)果‘）
%%?加窗處理?
win=chebw