PD.rar_PD_PD检测功能_PD雷达_radar_什么是PD测速

clc; clear all; close all; %参数设置 N=4096;%采样点 f0=10.^10;%载频 Tp=10*(10.^(-6));%脉冲宽度 B=10*(10.^6);%信号带宽 fs=100*(10.^6);%采样率 c=3*(10.^8); r0=3*(10.^3); a=2*r0/c; K=B./Tp; t=(0:N-1)/fs; %回波信号 Sb=rectpuls((t-a)-Tp/2,Tp).*exp(j*pi*K*((t-Tp/2-a).^2)).*exp(-2*j*pi*f0*a) r=c*t; figure; plot(r-3000,real(Sb)) figure; plot(r-3000,imag(Sb)) %回波信号频谱 N=4096; SF1=fft(Sb,N); fr=(0:N-1)*fs/N; figure; plot(fr,abs(SF1)); %参考信号 Sr=rectpuls(t-Tp/2,Tp).*exp(j*pi*K*((t-Tp/2).^2)); r=c*t; figure; plot(r,real(Sr)); figure; plot(r,imag(Sr)); %参考信号频谱 N=4096; SF2=fft(Sr,N); fr=(0:N-1)*fs/N; figure; plot(fr,abs(SF2)); %脉压处理 N=4096; S1=fft(Sb,N).*conj(fft(Sr,N)); So=ifft(S1,N); r=c*t/2; figure; plot(r,abs(So)) %真值转换为分贝值 N=4096; S1=fft(Sb,N).*conj(fft(Sr,N)); So=ifft(S1,N); r=c*t/2; figure; plot(r,10.*log10(abs(So))) %纵轴归一化db表示 N=4096; S1=fft(Sb,N).*conj(fft(Sr,N)); So=ifft(S1,N); r=c*t/2; figure; plot(r,10.*log10(abs(So)/1000))