雷达实验.zip_相控阵雷达_编码雷达_脉冲压缩回波_距离分辨力_雷达发射机

clear all; clc; close all; MHz=1e6; us=1e-6; %%产生线性调频信号 %产生线性调频信号参数 T=10*us; %带宽 B=10*MHz; %调频斜率 k=B/T; %采样间隔 fs=2*B; %产生标准的线性调频信号 t=-T/2:1/fs:T/2-1/fs; st=exp(1i*pi*k*t.^2); figure(1);plot(t,real(st)); hold on; plot(t,imag(st),'r'); L=T*fs; N=2^nextpow2(L); fst=fft(st,N); f=(0:N-1).*fs/N; figure(2); plot(f,abs(fst)); %%时域脉冲压缩处理 %不加窗处理 ht=conj(st(L:-1:1)); yt=conv(ht,st); %画出脉冲压缩的输出的幅度(纵轴用dB) tt=1:2*L; figure(3); plot(tt,20*log10(real(yt)));axis([0,400,-50,50]); fst1=fft(st); N1=length(fst1); f1=(0:N1-1).*fs/N1; figure(4); plot(f1,abs(fst1)); zhufubi=rms(abs(yt)) % %加窗处理频域 % wn=hamming(N)'; % figure(5);plot(wn); % %对fst加窗处理 % wh(1:N/2)=wn(1:N/2).*fst(N/2+1:N); % wh(N/2+1:N)=wn(N/2+1:N).*fst(1:N/2); % %作IFFT,得到wht % wht=ifft(wh); % %进行匹配滤波处理 % yt2=conv(wht,ht); % figure(6); % plot(abs(yt2)); w=zeros(N,1); w1=hamming(B/(fs/N)); w(1:B/(fs/N))=w1; w=circshift(w,-N/4); figure(8); plot(w); wHt=w'.*fst; wht=ifft(wHt); yt2=conv(wht,ht); figure(9); plot(20*log10(abs(yt2))); zhufubi2=rms(abs(yt2)) %分析主副比 %分析脉冲宽度的变化 % %%频域加窗脉冲压缩 % % %确定DFT的点数 % % %对st作M点DFT % fst1=?? % %对st的频域求共轭 % Hf=???; % %对Hf进行加窗 % w1=hamming(M); % wHf=w1.*Hf; % %在频域做频谱的乘积 % Yf=fst1.*wHf; % %作IFFT % yt=ifft(Yf); % % %加噪声 % %作脉冲压缩处理