MATLAB 二脉冲对消、重频参差滤波器设计仿真

雷达原理与系统作业，仅供参考，欢迎讨论，杜绝完全照搬。

一、 问题提出

二、问题解决

1.二脉冲对消原理

2.重频参差滤波器原理：

3.仿真实验

Matlab代码：

c=3e8;%光速
PRF=2000;%脉冲重复频率
Tr=1/PRF;%脉冲重复间隔
two_pulseH=@(w) 1-exp(-1j.*w.*Tr);%二脉冲对消系统函数
Doppler_w=@(v,lambda) 2.*v./lambda;%多普勒频率
lambda=c/(3e9);%S 波段雷达 2-4Ghz
v=0:0.1:800;%速度刻度
w=2.*pi.*Doppler_w(v,lambda);%速度对应的多普勒角频率
plot(v,20.*log10(abs(two_pulseH(w))))%幅度取对数
ylim([-30,7])
xlabel("径向速度/(m/s)")
ylabel("|H(j\omega)|/dB")
grid on

%枚举法，可能要5分钟左右
for i=1:300
for j=1:300
c=3e8;%光速
lambda=c/(3e9);%S 波段雷达 2-4Ghz
% PRFn=[2000,1754];%脉冲重复频率
% PRFn=[2910,2560];%脉冲重复频率
PRFn=[i*10,j*10];%脉冲重复频率
wn=[1,-2,1];%滤波器系数
Tn=1./PRFn;%脉冲重复间隔
tn=cumsum(Tn);%累计求和
tn=[0,tn];%初始为 0
H=@(w) sum(wn'.*exp(-1j.*w.*tn'));%系统函数
Doppler_w=@(v,lambda) 2.*v./lambda;%多普勒频率
v=0:0.1:800;%速度刻度
w=2.*pi.*Doppler_w(v,lambda);%速度对应的多普勒角频率
result=abs(H(w));%幅度
dBresult=20.*log10(result);%取对数
% plot(v,dBresult)
% ylim([-30,20])
% xlabel("径向速度/(m/s)")
% ylabel("|H(j\omega)|/dB")
% grid on
RMSE(i,j)=sqrt(sum((result-mean(result)).^2)./length(result));
end
end
mesh(RMSE)%画图
xlabel("PRF1")
ylabel("PRF2")
zlabel("RMSE")
colorbar
[m,n]=find(RMSE==min(min(RMSE)));%找到最小的序号，频率时在此基础上乘以 2