因为最近在做宽带数字信号的数字下变频，所以重点看了这一方面的论文。

下文重点对IEEE 上的一篇论文Real-time wideband DDC based on parallel architecture in synthetic instrument 做了个读书笔记，并对自己在这方面的工作做了一个小结。

主要内容：这篇论文主要解决宽带数字中频信号处理过程中的一些问题，并行处理结构可以有效的降低采样率，多通道并行处理则可以有效降低处理成本。

该论文介绍了数字下变频技术，分析了多种数控振荡器导致输出相位的不同，并总结了减少信噪比的误差因素，提供了特殊工程设计的指数分布依据。

图1DDC总体结构
原理：数字下变频器(D DC)是接收机A／D变换后，首先要完成的处理工作，一般的DDC 由本地振荡器(NCO)、混频器、低通滤波器和抽取器组成.主要作用：其一是把中频信号变为零中频信号；其二是降低采样率。

从频谱上看，数字下变频将A/D采样后信号从中频变换,到基带。

这样的处理由两步完成：首先是将输入信号与正交载波相乘，然后进行数字滤波滤除不需要的频率分量。

NCO，混频器，数字滤波器速率要等于采样率，采样率低于600MHz，很难实时的在FPGA中进行处理。

结构：
NCO：产生正余弦序列，即I/Q两路信号。

CIC滤波器，适用于系统中的第一级抽取和进行大的抽取因子的抽取工作，并降低速率。

FIR滤波器，完成对整个信道的整形滤波。

DDC主要有三种实现途径：采用专用芯片、自制专用芯片、基于DSP或FPGA等通用芯片。

论文主体：
1.并行处理结构:抽取滤波器模块通过多相滤波器结构降低采样率和实现低通滤波。

图2 DDC并行处理结构
2.通过合成，所有的i I 和i Q 信号将混合成最后的I 和Q 路信号。

图3 多路信号的合成
3.NCO 输出序列的相位差分析
NCO 输出多路正余弦信号，各有不同的相位差，但在相位分离后有相同的采样率。

以正弦信号为例，等同于,...2,1,0),2sin()(=∙=n nT f n d s π分析此式，很容易发现多路信号间存在相位差。

4.量化影响和误差分析
离散信号的精确性与寄存器长度相关，寄存器长度越长，精确度越高，硬件实现却越复杂。

多路相位合成滤波器在DDC 并行处理过程中不会产生其他噪声，因此在并行处理过程中信噪比的分析等同于传统的DDC 结构。

我做的工作：
通过matlab 程序了解数字下变频几个过程，重点在CIC 滤波器，HB 滤波器和FIR 滤波器。

基于matlab 的数字下变频设计步骤。

步骤一：确定数字频率合成器输出频率 步骤二：确定总抽取率D
步骤三：确定各级滤波器设计参数 CIC 滤波器（所谓的积分梳妆滤波器）：抽取率、级数； 补偿滤波器和FIR 滤波器：抽取率、阶数及系数。

数控振荡器
NCO 采用的直接数字频率合成技术(DDS)是一种实用的频率合成技术，DDS 由相位概念出发直接合成所需波形的一种新的频率合成技术。

DDS 合成技术采用了简便和有效的查表法， 即根据各个正弦波相位计算好相位对应的正弦值，存在ROM 中。

随着时钟周期变换，相位累加器不断频率控制字累加，送出相应的相位累加值；按照产生的相位累加值在正弦数据存储单元中查出对应的正弦幅度值并输出。

当一个周期相位满后，查找表回到初始值。

由matlab 程序得到NCO 信号如图1。

0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
x 10
-6
-200
-1000100
200NCO 信号
0.1
0.2
0.3
0.4
0.50.6
0.7
0.8
0.9
1
x 10
-6
-200
-1000100
200时间/s
图1
CIC 滤波器
CIC 滤波器是一种重要的用于采样率变换的高效率滤波器，结构整齐，运算简捷，硬件资源利用率很高，常用于高采样率数据端口和高变换比率的场合，即靠近模拟部分的一端。

由matlab 的CIC 滤波器程序可以得到CIC 滤波器幅频特性、相频特性：
0.5
1
1.5
2
-2002040
60H1(z)的幅频特性曲线0
0.5
1
1.5
2
-2-101
2H1(z)的相频特性曲线
0.5
1
1.5
2
-40
-20
20H2(z)的幅频特性曲线0
0.51 1.52
-2-1
1
2H2(z)的相频特性曲线
0.5
1
1.5
2
-60-40-200
20CIC 滤波器的幅频特性曲线00.51 1.52
-4
-202
4CIC 滤波器的相频特性曲线
0.5
1
1.5
2
-200
-100
1004级CIC 滤波器的幅频特性曲线0
0.5
1
1.5
2
-4-2
2
44级CIC 滤波器的相频特性曲线
图2
CIC 算法模型：
CIC 滤波器的传递函数为:
N
N RM N c
N
l z z H H z )1()1()(H 1----==N
RM k K z ][1
∑-=-= N 阶CIC 抽取滤波器的幅频响应为:
N
jw ]/2)
Sa(R/2)Sa([
|)H(e |ωω=
HB （半带滤波器)
HB 滤波器非常多的用于实现D=2M 倍的抽取，计算效率比普通FIR 效率高出近一倍，时间也更快。

HB 性质：
5
.0)()(1)(H 2
/)(=-=-πωπj j jw e
H e H e
HB 的冲激响应偶数点全是零(零点除外)，所以HB 只有一般FIR 运算量的二分之一。

由matlab 的HB 滤波器程序可以得到HB 滤波器幅频特性、相频特性：
00.5
1
-100
-500
50HB0滤波器的幅度曲线00.5
1
-3000
-2000
-1000
0HB0滤波器的相位曲线
0.5
1
-300
-200-1000
100HB4滤波器的幅度曲线
00.51
-6000
-4000
-2000
0HB4滤波器的相位曲线
FIR 滤波器
FIR 滤波器设计方法以直接逼近所需离散时间系统的频率响应为基础，主要是选择有限长度的h(n)，尽可能逼近传输函数。

线性相位FIR 滤波器单位冲激响应h(n)的设计常用等纹波法、窗函数法或频率抽样法。

窗函数法是指定连续的理想频率响应，然后用积分方法求出理想滤波器的单位抽样响应。

容易做到线性相位、稳定，可以设计各种特殊类型的滤波器，设计方法简单。

但是幅频指标不理想，不易控制边缘频率。

其中采用汉宁(hanli)窗设计的滤波器，主瓣和旁瓣有大约50dB 的差距，并且过渡带宽、通带波动比较小，而且其旁瓣下降比较快，基本满足设计要求。

等波纹算法就是使滤波器的频率响应在所感兴趣的频率范围内与理想滤波器的频率响应之间的最大逼近误差最小。

MA TLAB 中的函数REMEZ 是用于设计最佳滤波器的，而REMEZORD 函数则用来计算所需的滤波器阶数。

下图给出的matlab 程序是用汉宁窗实现31阶和127阶FIR 低通滤波器的频响应：
0123456
7
x 10
5
-150
-100-50031阶FIR 低通滤波器的幅频响应
1
2
345
6
7
x 10
5
-2000
-1500-1000-500031阶FIR 低通滤波器的相频响应
频率/Hz
0123456
7
x 10
5
-150
-100-50031阶FIR 低通滤波器的幅频响应
1
2
345
6
7
x 10
5
-2000
-1500-1000-500031阶FIR 低通滤波器的相频响应
频率/Hz。