清除全部;
参数定义
c=3e 8; % speed of light
pi=3.1415926;
j=sqrt(-1;
Tp=1e-6; % transmitted pulse width
fc=1e 9; %运营公司frequency
Br=50.e 6; %传输的bandwidth
Fs=200.e 6; % a/d样本速率
kr=Br/Tp; % range chirp rate
Nr=Tp*Fs;
Ni=1:Nr;
tr=(Ni-NR/2 ) *Tp/Nr;
%Generate a Chirp pulse
%===========================
SIG=exp(j*pi*kr* ) tr ).^2);
图形;
subplot(2、3、1 ); 打印(实际) SIG );
subplot(2、3、2 ); 打印(imag ) SIG );
% Spectrum of this Chirp pulse
%===========================
SIG_spec=FFTshift(FFT(SIG );
subplot(2、3、3 ); plot(ABS ) SIG_spec );
%Match filtering in time domain
%===========================
MF_SIG=conj(fliplr ) SIG );
subplot(2、3、4 ); 打印(real ) MF_SIG );
subplot(2、3、5 ); plot(imag ) MF_SIG );
MF_out=conv(SIG,mf_sig );
subplot(2、3、6 ); plot(ABS(MF_out ) NR/21:NR*3/2) );
%===================================
% matchfilteringinfrequencydomain
%===================================
fr=(Ni-NR/2 ) *Fs/Nr;
MF_spec=exp(j*pi*fr.^2/kr );
MF_spec_out=IFFT(FFTshift ) SIG_spec.*MF_spec );
图形;
subplot(1、2、1 ); plot(ABS ) MF_spec_out );
SIG _ un=20 * log 10 (ABS (MF_spec_out )/max (ABS ) MF _ spec _ out ) );
for i=1 :Nr
if(SIG_un(I )-40 ) ) ) ) ) )。
SIG_un(I )=-40;
结束
结束
subplot(1、2、2 ); plot(SIG_un );
D86