多相信道化接收机

最近在做一个多相信道化接收机,实现的功能是把BW=25MHz@75MHz的中频信号,用f_s=100MSaps的采样频率采样,再对数字中频信号做数字下变频,然后对下变频后的基带信号做信道划分,总共分16个信道。


首先是BW=25MHz@75MHz的模拟中频信号被f_s=100MSaps采样后的数字中频信号为BW=25MHz@25MHz,这是由于带通采样定理决定的,在此不再详述。从这个数字中频信号开始,我开始仔细介绍整个信号处理流程。

为了方便介绍,我假设数字中频信号为BW=25MHz@25MHz,持续时间T=2.6ms的Chirp信号,采样率f_s=100MSaps。首先要做的就是数字下变频,也就是说把中频信号做IQ调制。其信号流程图如下所示。

多相信道化接收机_第1张图片

NCO就是本振,左半部分就是做频移。还记得傅里叶变换的频移性质吧,也就是假设f(t)\leftrightarrow F(w)互为傅里叶变换对,在时域乘以e^{-jw_{c}t}后的傅里叶变换,对应的频域为F(w-w_{c})。对于上述框图中频移部分就是应用了傅里叶变换的这个性质。乘以e^{-jw_{c}t}就等于乘以cos(w_{c}t)-jsin(w_{c}t),因此把复数乘变成两个通道分别乘以cos(w_{c}t)-sin(w_{c}t)。如下图所示,同相通道乘以cos(w_{c}t)等于乘以\frac{1}{2}(e^{jw_{c}t}+e^{-jw_{c}t}),也就是频谱分别向正、负频率方向移动w_{c}。同理正交通道乘以sin(w_{c}t)等于乘以\frac{1}{2j}(e^{jw_{c}t}-e^{-jw_{c}t}),也就是频谱分别向正频率方向移动w_{c},向负频率方向移动w_{c}并且幅度取负,与此同时,相位变换-\pi(虚数j的性质)。相频图就不画了(因为我画不好)。其中向正频率方向频移的认为是镜像频谱,用低通滤波器滤除即可。

 

多相信道化接收机_第2张图片

为了能够恢复信号、降低采样率,我们需要把镜像频谱滤除并做抽取,下图以同相通道为例。设计低通滤波器只滤取0\sim 12.5MHz的信号,再做抽取将采样率降低,由于输入信号BW=25MHz@25MHz,采样率f_s=100MSaps,因此这一步的抽取可以抽4,即采样率降为25MHz。这一阶段处理后的I、Q信号BW=12.5MHz@0MHz,采样率f_{s}=25MHz

多相信道化接收机_第3张图片

下面用Matlab仿真展示数字下变频过程中各个阶段信号的谱、镜像滤波器频率响应。

多相信道化接收机_第4张图片

多相信道化接收机_第5张图片

 

多相信道化接收机_第6张图片

多相信道化接收机_第7张图片

 

我们感兴趣的信号覆盖了第一、二奈奎斯特区域(0\sim 25MHz),信道化就是把这段感兴趣的频段分成K个通道。信道化的本质如下所示。我的案例里25MHz的频段分16个通道,每个通道1.5625MHz,下图中每个分支前的频移项为e^{-jw_kn}w_k=2\pi k\frac{25MHz}{16}。每个分支的低通滤波器都一样,都是滤取一个通道的滤波器,该滤波器保证能滤取一个信道带宽(\frac{25MHz}{16})的频段。由于每个信道的带宽很低,因此可以进一步降低采样率。

 

多相信道化接收机_第8张图片

上述方案为一种低效方案,根据该框图的信号处理流程,我们可以看出每个信道的信号y_k(n)=[x(n)*e^{-jw_kn}]*h_{LP}(n)|_{n=mD},其中x(n)是基带复数信号, h_{LP}(n)是滤波器系数,D是抽取系数。经过下述数学等效变换:y_k(n)=\sum_rx(mD-r)e^{-jw_k(mD-r)}h_{LP}(n)

=\sum_{p=0}^{K-1}\sum_{l=0}^{L-1}x(mD-lK-p)e^{-jw_k(mD-lK-p)}h_{LP}(lK+p)

x_p(m)=x(mD-p)h_p(m)=h_{LP}(mK+p)K=FD,则有:

y_k(n)=\sum_{p=0}^{K-1}\sum_{l=0}^{L-1}x(mD-lFD-p)e^{-jw_k(mD-lFD)}h_{LP}(lK+p)e^{jw_kp}

y_k(n)=\sum_{p=0}^{K-1}[\sum_{l=0}^{L-1}x_p(m-lF)e^{-jw_k(m-lF)D}h_{p}(l)]e^{jw_kp}

y_k(n)=\sum_{p=0}^{K-1}[\sum_{l=0}^{(L-1)F}x_p(m-l)e^{-jw_k(m-l)D}h_{p}(\frac{l}{F})]e^{jw_kp}

h_p^{'}(l)=h_{p}(\frac{l}{F}),也就是说h_p^{'}(l)是由原型滤波器h_{LP}(n)的多相分支经F内插后的结果。

y_k(n)=\sum_{p=0}^{K-1}[x_p(m)e^{-jw_kmD}*h_{p}^{'}(m)]e^{jw_kp}

K=FDK=16,D=8,F=2)代入上式后,信号处理流程相当于下面框图。

多相信道化接收机_第9张图片

原型滤波器频率响应如下图所示。

多相信道化接收机_第10张图片

整个信道化处理流程走完后,各个信道信号的幅频响应如下图所示,

多相信道化接收机_第11张图片

当然我是为了展示方便,把各个信道的信号又频移回了原始的位置。


以上是我的一个小随笔,可能还是有错误希望走过路过的各位大神提意见。

我的邮箱:[email protected]

另外,matlab的代码我会发上来,但是因为我平时想下资料,所以我要收一个积分咯。

 

你可能感兴趣的:(信号与系统,信道化接收机)