雷达原理---直接数字中频正交采样

一、模拟正交鉴相

模拟正交鉴相又称为零中频处理。这里，零中频是指由于相干振荡器的频率与中频信号的中心频率相等，当不考虑多普勒频移时，其差频为零。

1. 原理

传统雷达对接收信号经过模拟混频、滤波得到中频信号，再经过模拟正交相干检波器得到基带I、Q信号。

中频信号分为同相的两路信号：其中一路与本振源输出的信号进行混频、滤波后，得到同相支路基带信号I(t)；另一路与经过90°相移后的本振源信号进行混频、滤波后，得到正交支路基带信号Q(t)。I(t)和Q(t)分别经过A/D变换器后，输出同相的数字基带信号I(n)和正交的数字基带信号Q(n)。

2. 优点

可以处理较宽的基带信号，对A/D变换器的转换速率要求也相对较低。

3. 缺点

由于模拟正交相干检波器需要两路完全正交的本振源、两个混频器和滤波器，如果这两路模拟器件的幅度和相位特性不一致，将导致I、Q不平衡，产生镜频分量。

• 若中频输入信号模型为$s(t)=cos2\pi (f_0+f_d)t$，则在理想情况下，正交两路混频器的参考信号和输出的基带信号为：
  

• 若两个本振信号存在幅度相对误差${\epsilon }_A$和正交相位误差${\epsilon }_{\phi }$（即相位差不等于90°），正交两路混频器的参考信号和输出的基带信号为：
  
  
  因此，模拟正交鉴相器在高精度、高性能的现代全相参雷达接收机的应用中受到了一定的限制。

二、数字正交鉴相

数字正交鉴相（又称为数字下变频）的基本原理是首先对中频模拟信号直接进行A/D变换，然后实现I/Q数字分离。数字鉴相的最大优点是可以得到很高的I/Q正交精度和稳定度。

实现数字鉴相的方法很多。两种常用的方法：

1. 带通采样数字混频低通滤波法

这种方法的原理类似于零中频模拟正交鉴相法。首先用A/D变换器对中频信号进行直接采样，采样数据分别与来自数字压控振荡器（NCO）的两路相位相差90°的数字本振$cos(n{w}_{0}T)$和$sin(n{w}_{0}T)$进行数字混频，混频后的数据分别经过数字低通滤波器输出基带同相数字信号I(n)和基带正交数字信号Q(n)。

由于两路相干振荡信号、混频和低通滤波器等都是用数字方法实现的，这就消除了采用模拟方法引入的直流偏移以及I、Q通道的增益、相位失配等问题，从而保证了I、Q通道在幅度一致性和相位正交性方面的精度远高于传统的模拟方法，这使得采用数字化接收机的信号处理精度和稳定性都明显优于传统模拟接收机的信号处理的处理精度和稳定性。

2. 带通采样数字插值滤波法

通过选用适当的采样频率对中频信号进行A/D变换，可以交替得到I(n)和Q(n)，然后通过数字内插滤波器进行内插运算，从而分别输出基带同相数字信号I(n)和基带正交数字信号Q(n)。

已知输入的中频信号表达式为
$$s(t)=A(t)\cdot cos[{\omega }_{0}t+\varphi (t)]$$
式中，$A(t)$和$\varphi (t)$分别为信号的幅度调制信号和相位调制函数。

直接用A/D变换器对中频带通信号采样，为了避免混叠，采样频率$f_s$应满足$f_s\ge 2B$。采样频率$f_s$与中频信号的中心频率$f_0$以及信号带宽B的关系为
$$f_s=\frac{1}{T}=\frac{4}{M+1}{f}_0$$
$$f_s\ge 2B$$
式中，T为采样周期，$M=0,1,2,3,4,\cdots$，则有
$${\omega }_{0}T=2\pi f_0\frac{M+1}{4f_0}=\frac{\pi }{2}(M+1)$$
如果取$n=M+1$，则可以得到
$$\begin{array}{rl}s(nT)& =A(nT)\cdot cos[n{\omega }_{0}T+\varphi (nT)]\\ & =A(nT)\cdot cos(nT)cos(\frac{n\pi }{2})-A(nT)\cdot sin(nT)sin(\frac{n\pi }{2})\\ & =I(n)\cdot cos(\frac{n\pi }{2})-Q(n)\cdot sin(\frac{n\pi }{2})\\ & =\{\begin{array}{ll}(-1{)}^{\frac{n}{2}}I(n)& n为偶数\\ (-1{)}^{\frac{n-1}{2}}Q(n)& n为奇数\end{array}\end{array}$$
式中，$I(n)=A(nT)\cdot cos(nT)；Q(n)=A(nT)\cdot sin(nT)。$因此，$n$为偶数时可得到$I(n)$，$n$为奇数时可得到$Q(n)$，其输出顺序为$I(0),Q(1),-I(2),-Q(3),I(4),Q(5),\cdots$。由此可见，$I(n)和Q(n)$交替出现，它们在时间上差1个采样周期T。采用数字内插滤波方法，将两序列符号统一变正，群延时对齐，即可分离出$I(n)和Q(n)$信号。

这种数字内插滤波器法的一个突出优点是不需要数字压控振荡器（NCO）即可保证较高的精度和稳定性。在某些应用中，为了方便起见，将经过A/D变换器中频采样的数字信号$s(nT)$直接送至接收机后续的数字信号处理器，由数字信号处理器完成内插滤波和$I(n)、Q(n)$分离。

三、MATLAB仿真

1. 带通采样数字混频低通滤波法的时域实现

clear all;
close all;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%% 时域数字正交相干检波 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%********** 线性调频信号 ***********%%
T=10e-6;            %信号时宽
B=30e6;             %信号带宽
fs=4*B;             %采样频率
f0=3/4*fs;          %载频频率
N=fix(T*fs);        %采样点数

t=(-N/2:N/2-1)/fs;
sr=cos(2*pi*f0*t+pi*(B/T).*t.^2);  %sr是以fs=4B的采样频率采样后的数字信号

figure; plot(t*1e6,sr); title('线性调频信号波形'); xlabel('t/us'); ylabel('sr');
%%********** 1.混频 **********%%
x_I=sr.*cos(2*pi*f0*t);
x_Q=-sr.*sin(2*pi*f0*t);
%%********** 2.低通滤波 **********%%
%使用Parks McClellan算法设计FIR滤波器:Fs=4B,截止频率为B/2，过渡带宽度为B/10.
fl=200;
fbe=[0 (B/2)/(fs/2) (B/2+B/10)/(fs/2) 1];
Damps=[1 1 0 0];
b=firpm(fl,fbe,Damps);
figure;
freqz(b);
title('FIR滤波器的幅频和相频特性');
%使用由分子和分母系数 b 和 a=1 定义的有理传递函数对输入数据 x 进行滤波。
x_I_filter=filter(b,1, x_I );
figure;subplot(2,1,1); plot(t*1e6,x_I_filter);xlabel('t/us');ylabel('x_I_ filter');title('低通滤波后的I路信号');
x_Q_filter=filter(b,1, x_Q );
subplot(2,1,2); plot(t*1e6,x_Q_filter);xlabel('t/us');ylabel('x_Q_ filter');title('低通滤波后的Q路信号');
%%********** 3.抽样 **********%%
x_I_ex=x_I_filter(1:2:end);
figure; subplot(2,1,1);plot(t(1:2:end)*1e6,x_I_ex);xlabel('t/us');ylabel('x_I_ ex');title('抽样后的I路信号');
x_Q_ex=x_Q_filter(1:2:end);
subplot(2,1,2); plot(t(1:2:end)*1e6,x_Q_ex);xlabel('t/us');ylabel('x_Q_ ex');title('抽样后的Q路信号');
%%************ 4.合成 ***********%%
x_ex= x_I_ex+1j*x_Q_ex;
figure;plot(t(1:2:end)*1e6,x_ex);xlabel('t/us');ylabel('x_ ex');title('正交下变频的最终合成信号');

2. 带通采样数字混频低通滤波法的频域实现

clear all;
close all;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%% 频域数字正交相干检波 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%********** 线性调频信号波形 **********%%
T=10e-6;            %信号时宽
B=30e6;             %信号带宽
fs=4*B;             %采样频率
f0=3/4*fs;          %载频频率
N=fix(T*fs);        %采样点数

t=(-N/2:N/2-1)/fs;
sr=cos(2*pi*f0*t+pi*(B/T).*t.^2);
figure;
subplot(2,1,1);plot(t*1e6,sr); title('线性调频信号波形'); xlabel('t/us'); ylabel('sr');
%%********** 线性调频信号的频谱 **********%%
L=2^nextpow2(N);
Sf1=fftshift(abs(fft(sr,L)));
f=linspace(-fs/2,fs/2,L);
subplot(2,1,2);plot(f*1e-6,Sf1); title('线性调频信号频谱'); xlabel('f/MHz'); ylabel('Sf1');
%%********** 1.频谱搬移 **********%%
Sf=fftshift(abs(fft(sr.*exp(-1j*2*pi*f0.*t),L)));
figure;plot(f*1e-6,Sf); title('频谱搬移后的信号'); xlabel('f/MHz'); ylabel('Sf');
%%********** 2.频率滤波+抽取 **********%%
Sf_ex=[Sf(1:L/4),Sf(3/4*L+1:end)];
figure;plot(f(1:2:end)*1e-6,fftshift(Sf_ex)); title('滤波抽取后的信号'); xlabel('f/MHz'); ylabel('Sf_ ex');