资源描述：
《matlab仿真程序》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、窗型选择仿真程序:clear,clcbw=3e6;%信号带宽T=le-4;%信号脉冲宽度A=2;%信号幅度fs=4氺bw;lfft=round(T*fs);%采样点数lfft=2nextpow2(lfft);dt=l/fs;%采样间隔f0=le6;(0:1fft-1);%时域采样点s=A*cxp(j*2*pi*fO*t+j*pi*bw*t.*t/T);%产生线性调频信号S=(fft(s));%线性调频信号的傅立叶变换fftH=conj(S);%闪配滤波器的频率响极Y=S.*H;%线性调频信号的频域匹配滤波输出y=fftshift(ifft(Y))线性
2、调频信号的时域匹配滤波输出%对chirp信号进行时域加权hl=(triang(lfft))’;%三角窗函数sl=s.氺hl;Sl=fft(si);Hl=conj(SI);Y1=S1.*H1;yl=fftshift(ifft(Yl))加三角窗后的线性调频信号的吋域匹配滤波输出h2=(harming(lfft))’;%汉宁窗函数s2=s.*h2;S2=fft(s2);H2=conj(S2);Y2=S2.*H2;y2=fftshift(ifft(Y2));%加汉宁窗后的线性调频信号的时域匹配滤波输出h3=(hamming(lfft))’;%海明窗函数s3=
3、s.賦S3=fft(s3);H3=conj(S3);Y3=S3.*H3;y3=fftshift(ifft(Y3));%加海明窗后的线性调频信号的吋域匹配滤波输出figure;subplot(3,1,1),plot(t,real(s)),titleCchirpsignal’)；subplot(3,1,2),plot(q,abs⑸),title(’线性调频信号幅度谱’);subplot(3,1,3),plot(q,angle(S)),title(’线性调频信号相位谱’)；figure;subplot(2,1,1),plot(q,abs(H));title
4、(’MF的幅度谱’);subplot(2,1,2),plot(q,angle(H));title(’MF的相位谱’)；figure;subplot(3,1,1),plot(t,real(y)),title(’脉压信号’);subplot(3,1,2),plot(q,abs(Y)),title(’脉压信号幅度谱’)；subplot(3,1,3),plot(q,angle(Y)),title(’脉压信号相位谱’);figure;subplot(2,2,1),plot(t,20*logl0(abs(y)/max(abs(y))));title(’未加窗时的
5、时域输出’)；subplot(2,2,2),plot(t,20*logl0(abs(y1)/max(abs(yl))));title(’加三角窗时的时域输出’)；subplot(2,2,3),plot(t,20*logl0(abs(y2)/max(abs(y2))));title(’加汉宁窗吋的吋域输出’);subplot(2,2,4),plot(t,20*logl0(abs(y3)/max(abs(y3))));title(’力口海明窗时的时域输出’);叠加3个频移多普勒干扰程序：BandWidth=10e6;%发射信号带宽TimeWidth=20
6、e-6;%发射信号时宽mu=BandVidth/TimeWidth%调频率Fs=2*BcindWidth;%采样频率Ts=l/Fs;Ns=fix(Fs*TimeWidth);%计算一个脉冲周期的采样点数400；N=1024;%FFT点数t=0:Ts:TimeWidth-Ts;y=exp(j*pi*mu*t，2);%产生LFM信号h二zeros(1，Ns);fori=l:Nsh(i)=conj(y(Ns~i+1));endfd=6e6;yl=exp(j*2*pi*(fd*t+0.2));%频移干扰信号fdl=+le6;y2=exp(j*2*pi*(f
7、dl*t+O.5*mu*t.2));%频移干扰信fd2=+10e6;y3=exp(j*2*pi*(fd2*t+0.2));%频移干扰信y=y1+exp(j*pi*mu*t."2)+y2+y3;%产生叠加了干扰的LFM信号yfft=fft(y,1024);win二hamming(Ns)’；h_w=h.*win;hfft_w=fft(h_w,1024);ycomp=abs(ifft(yfft.*hfftw));%脉冲压缩maxvall=max(ycomp);ycomp_w=eps+ycomp./maxvall;%利用ycomp的最大值归一化tt=0:Ts
8、:2*TimeWidth~Ts;plot(tt,ycompw(l:2*Ns),’b’)xlabel(’t-s
