脉冲压缩程序实验报告.docx

1024点的脉冲压缩程序一 实验原理chirp信号【程序】Tp=2e-6;%脉冲宽度Bandwith = 100e6;%频域带宽Fs = 125e6;%抽样频率r = Bandwith/Tp;%调频率fc = 100;%载波频率N_FFT=1024;%fft抽样点数%定义LFM函数n=round(Tp*Fs);t=linspace(-n/2,n/2,n)/Fs;s_shiyu=exp(1j*pi.*r.*t.^2);% figure;% plot(t,real(s_shiyu));% hold on;plot(t,imag(s_shiyu),r);% title(LFM信号时域-实部和虚部);% xlabel(t/s);s_pinyu=fftshift(fft(s_shiyu,N_FFT));f=linspace(-Fs/2,Fs/2,N_FFT);% figure;% plot(f,real(s_pinyu));% hold on;plot(f,imag(s_pinyu),r);% title(LFM信号pin域-实部和虚部);% xlabel(f/Hz);figureplot(t,abs(s_shiyu));title(LFM信号幅度-shiyu);xlabel(t/s);figure;plot(angle(s_shiyu));title(LFM信号-shiyu);xlabel(t);ylabel(Phase);figureplot(f,abs(s_pinyu));title(LFM信号幅度-pinyu);xlabel(f/Hz);ylabel(幅度);figure;plot(angle(s_pinyu));title(LFM信号幅度-pinyu);xlabel(f);ylabel(Phase);【fftshift函数特点】fftshift是针对频域的，将FFT的DC分量移到频谱中心即对频域的图像，（假设用一条水平线和一条垂直线将频谱图分成四块）对这四块进行对角线的交换与反对角线的交换。记信号f(t)的起始时间t_start,终止时间为t_end，采样频率为Fs，信号的持续时间为duration=t_end-t_start,信号离散化造成的采样点数N=duration*Fs+1；时域的采样化会造成频域的周期化。经过Matlab的fft函数处理后，得到的数据的横坐标为 0：Fs/（N-1）:Fs，如果使用fftshift命令，零频率分量将会移到坐标中心。%Y = fftshift(X) rearranges the outputs of fft, fft2, and fftn % by moving the zero-frequency component to the center of the array. % It is useful for visualizing a Fourier transform with the zero-frequency % component in the middle of the spectrum.【结果】用matlab程序画出 LFM信号时域和频域的图形LFM信号实部和虚部的图像【分析】从LFM信号时域的表达式和傅里叶变换后频域的表达式可以看出，LFM信号的频域也满足LFM信号的特点，即也是LFM信号。时域幅度为常数；图形上表现为矩形窗相位时关于时间t的二次函数，开口朝上时，上式取最大值频域幅度响应若忽略波动，可以近似看作矩形窗函数。相位为关于W的二次函数，开口朝下匹配滤波器匹配滤波器的设计遵循最大信噪比准则，即滤波器的输出信号在某一时刻t0的瞬时功率与噪声平均功率之比达到最大。使得输出信噪比最大的条件是：脉冲压缩体制接收端的匹配滤波器频谱近似为对于LFM信号的压缩网络，其群时延变化与频率变化相反，信号的频率越大，延时越大。chirp信号经过匹配滤波器后的输出信号方法一：u(t)=s(t)h(t)【性质】当 t0=0时是偶函数【程序】syms w tf=1/(2*pi*r).*exp(1i*w*(t-t0));u_shiyu=int(f,w,s1,s2);Tp=2e-6;%脉冲宽度Bandwith = 100e6;%频域带宽Fs = 500e6;%抽样频率r = Bandwith/Tp;%调频率fc = 100;%载波频率N_FFT=1024;t=linspace(-2*N_FFT,2*N_FFT,4*N_FFT)/Fs;%LFM通过匹配滤波器后的输出信号w0=2*pi*fc;%载频为中心频率t0=Tp;%h(t)=s(t0-t)s1=w0-pi*Bandwith;s2=w0+pi*Bandwith;w=linspace(-4*s1,4*s2,4*N_FFT);【结果分析】时域图