一种基于DFT滤波器组结构的原型滤波器及设计方法

一种基于dft滤波器组结构的原型滤波器及设计方法
技术领域
1.本发明属于数字信号处理领域，适用于完全重建滤波器组中的原型滤波器设计，具体涉及一种基于dft滤波器组结构的原型滤波器及设计方法。


背景技术：

2.基于高效信道化结构的并行解调原理框图如图1所示，数字中频信号经过下变频处理后，利用多相滤波器组技术进行信道化分解，将宽带信号分解为多路并行窄带信号，然后进行子带并行匹配滤波处理。完成匹配滤波后进入信道化合成模块，将同步与均衡模块移至综合滤波和内插中间进行，此时各子带间数据仍处于并行低速状态。
3.1、多相分析滤波器组实现结构
4.m通道，d倍抽取的级联复指数滤波器组结构如图2所示，其中hm(z)与fm(z)，m＝0，
…
，m-1为带通滤波器，宽带信号经过不同选通频率的带通滤波器滤波后，被d倍抽取，为避免混叠，必须满足d≤m，然后通过复数旋转器将子带信号变频到以零频率为中心，即可得到m路子带信号。暂不考虑子带信号处理模块，综合部分是分析部分的逆过程，如果滤波器hm(z)与fm(z)能够满足功率互补条件，且完全滤除带外频谱，则该模型实现了对信号的完全重建。
5.下面推导各节点信号的频域表达式，对于不同频段的带通滤波器可以用其低通原型滤波器h(z)表示为
[0006][0007]
其中，wm＝e-j2π/m
，所以
[0008][0009]
对带通滤波后的信号进行d倍抽取，得到
[0010][0011]
然后进行下变频处理
[0012][0013]
将式(3)代入式(4)，即可得到各子带信号的频域表达式
[0014][0015]
不考虑子带信号处理模块，则有
[0016][0017][0018]
对上式进一步整理，可得
[0019][0020]
下面推导分析滤波器组的多相实现形式，首先对原型滤波器h(z)进行多相位分解：
[0021][0022]
其中
[0023][0024]
将式(9)代入式(3)，可得
[0025][0026]
当m是d的整数倍时，则上式可进一步变形为
[0027][0028]
其中[
·
]|
↓d表示序列进行d倍抽取，根据离散傅里叶变换的定义，式(11)中的累加可以用idft运算来实现，因此可得：
[0029][0030]
式(13)的实现框图如图3所示，其中idft运算采用更加高效的ifft实现，可以发现通过对原型滤波器进行多相分解，将抽取操作调整到了滤波计算之前，从而大大减小了运算速率。
[0031]
2、多相综合滤波器组实现结构
[0032]
多相综合滤波器组的推导方法与多相分析滤波器组大致相同，综合原型滤波器的多相形式为
[0033][0034]
其中
[0035][0036]
将式(14)代入式(7)中，可得
[0037][0038]
式(16)的最后一行等式表明，可以将子带信号yk(z)经过idft变换后，直接进行多相综合滤波，然后执行内插操作，由此实现了降速处理，其实现过程是多相分析滤波器组结构的逆过程。
[0039]
图4所示为m通道，d倍抽取的多相滤波器组总体实现结构，信号先抽取再滤波后，进行ifft变换，即可得到m路并行数据。综合模块首先对子带信号进行ifft变换，然后先滤波再内插，最后延迟累加，即可得到恢复的信号。级联的多相滤波器组结构为信号并行降速处理提供了较好的方案，且抽取倍数可以根据实际情况调整，只需满足d整除m的条件即可。
[0040]
3、滤波器组的完全重建条件
[0041]
根据公式(8)，可以得到y(z)与x(z)的关系，下面用矩阵形式描述式(8)，定义下列矩阵：
[0042][0043][0044][0045][0046][0047][0048]
则有
[0049][0050]
其中是分析及综合滤波器组中m个路径的总传递函数，f
th(m
×
1)
是所需的信号传递函数，而是混叠干扰，式(23)可以改写为
[0051][0052]
混叠能量被完全消除的条件为传递函数即
[0053][0054]
由于分析及综合滤波器组本质上是原型滤波器在频域的移位和叠加，所以可以去除式(25)中的频移因子和累加部分不会影响结论，因此混叠消除的条件变为
[0055][0056]
系统对输入信号产生无失真响应的条件为传递函数f
th(m
×
1)
＝d，即
[0057][0058]
其中滤波器c(z)满足功率互补条件。
[0059]
综上所述，式(26)和式(27)构成了分析及综合滤波器组的完全重建条件。
[0060]
当前的完全重建滤波器组设计中，常采用两个完全相同的根升余弦滤波器作为分析及综合原型滤波器，这是因为升余弦滤波器具有良好的功率互补特性，能够满足式(27)的要求，且设计简单易行，但是根升余弦滤波器的阻带衰减性能不好，仿真发现对于不同的滚降因子β，其阻带抑制均小于-100db，无法良好地满足式(26)要求的混叠消除条件，会给重建信号带来一定的误差。


技术实现要素：

[0061]
有鉴于此，本发明的目的是提供一种基于dft滤波器组结构的原型滤波器及设计方法，可以增大滤波器的阻带抑制，提高完全重建性能。
[0062]
一种基于dft滤波器组结构的原型滤波器，原型综合滤波器采用凯撒窗滤波器，原型分析滤波器采用升余弦函数滤波器。
[0063]
一种原型滤波器的设计方法，包括如下步骤：
[0064]
第一步，根据需要设置以下参数的值：滤波器组的信道个数k，原型分析滤波器滚降因子β和原型综合滤波器阻带最小衰减as；
[0065]
第二步，根据如下公式，计算原型综合滤波器的形状参数βk：
[0066]
βk＝0.1102(a
s-8.7)；
[0067]
第三步，设计原型分析滤波器和原型综合滤波器的阶数nf；
[0068]
第四步，根据参数β和nf设计升余弦函数的原型分析滤波器；
[0069]
根据参数βk和nf设计原型综合滤波器w(n)：
[0070][0071]
其中，n表示滤波器的采样点；i0(
·
)是第一类零阶变型贝塞尔函数。
[0072]
较佳的，当滤波器阶数的允许设计范围超过设定值时，采用最优化设计方法，根据如下公式计算得出阶数nf：
[0073][0074]
较佳的，当原型滤波器阶数的允许设计范围受限时，采用次优化设计方法，根据如下公式计算得出阶数nf：
[0075][0076]
本发明具有如下有益效果：
[0077]
本发明基于非最大抽取dft滤波器组的框架，提供了一种原型分析及综合滤波器及设计方法，该滤波器增大滤波器的阻带抑制，提高完全重建性能，预期比传统的采用两个根升余弦滤波器作为原型滤波器的方法在重建误差方面降低约2个数量级；设计方法分析了滤波器参数的约束优化问题，并分别设计了最优化设计方法和次优化设计方法，两种方法实现简单且性能良好。
附图说明
[0078]
图1为高效信道化结构并行解调原理框图；
[0079]
图2为m通道，d倍抽取复指数滤波器组模型；
[0080]
图3为分析滤波器组的多相实现结构；
[0081]
图4为m通道，d倍抽取多相滤波器组结构；
[0082]
图5为原型分析与综合滤波器的频率响应；
[0083]
图6为不同子带滤波器的频域分布；
[0084]
图7为传统根升余弦滤波器组的总功率谱；
[0085]
图8为改进方法滤波器组的总功率谱；
[0086]
图9为输入单脉冲信号与输出信号；
[0087]
图10为输入高斯白噪声信号与重建信号；
[0088]
图11为传统根升余弦滤波器组的重建误差；
[0089]
图12为改进滤波器组的重建误差。
具体实施方式
[0090]
下面结合附图并举实施例，对本发明进行详细描述。
[0091]
为了增大滤波器的阻带抑制，提高完全重建性能，本发明重新设计了原型分析及综合滤波器，其中综合滤波器利用凯撒窗滤波器取代根升余弦滤波器，从而具备了较强的阻带抑制能力，而分析端直接使用升余弦滤波器，根据升余弦滤波器的功率互补特性可得：
[0092][0093]
其中，h
rc
(z)表示采用升余弦函数设计的原型分析滤波器，为了使级联的完整滤波器组仍能满足式(27)的无失真响应条件，综合滤波器需要满足以下等式：
[0094]hrc
(z)f
ks
(z)＝h
rc
(z)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(29)
[0095]
其中，f
ks
(z)表示采用凯撒窗设计的原型综合滤波器。式(29)表明综合滤波器的通带宽度要包含分析滤波器的通带和过渡带，且其通带波纹要尽可能小，这显然是容易实现
的。
[0096]
下面讨论原型滤波器设计的参数优化问题，凯撒窗函数表示为
[0097][0098]
其中i0(
·
)是第一类零阶变型贝塞尔函数，βk是可自由选择的形状参数。
[0099]
根据凯撒窗函数性质可得其形状参数βk：
[0100]
βk＝0.1102(a
s-8.7)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(31)
[0101]
归一化过渡带宽δω表示为：
[0102][0103]
其中nf是原型滤波器阶数。式(29)要求综合滤波器的过渡带处在分析滤波器及其镜像的阻带范围内，因此有：
[0104][0105]
结合式(32)和式(33)可得参数设计的约束表达式为：
[0106][0107]
在实际工程设计中，需要预先设定的指标包括信道个数k，分析滤波器的滚降因子β和综合滤波器的阻带衰减as。本发明提供最优化与次优化两种参数设计方法，当滤波器阶数的允许设计范围较大时，采用最优化设计方法，此时：
[0108][0109]
结合式(31)和式(35)即可设计出最优化原型分析与综合滤波器。
[0110]
当滤波器阶数的允许设计范围受限时，采用次优化设计方法。考虑到与多相分析滤波器相邻的两个镜像已经受到升余弦滤波器一定程度的抑制，造成的频谱混叠能量是有限的，而多相综合滤波器的主要作用是最大程度地抑制远处的镜像叠加造成的混叠，因此可以适当放宽原型综合滤波器过渡带宽的限制，达到降低滤波器阶数的目的。
[0111]
将式(33)的过渡带约束变为：
[0112][0113]
此时设计的原型滤波器阶数nf′
为：
[0114][0115]
结合式(31)和式(37)即可设计出次优化原型分析与综合滤波器。
[0116]
比较两种设计方法的原型滤波器阶数，可得：
[0117]
[0118]
从式(38)可以看出，原型分析滤波器的滚降因子β越大，次优化设计的滤波器阶数相对越小，从而具有更高的实用价值。
[0119]
下面总结设计本原型滤波器的实现步骤：
[0120]
第一步，根据需要设置以下参数的值：滤波器组的信道个数k，原型分析滤波器滚降因子β和原型综合滤波器阻带最小衰减as；
[0121]
第二步，根据公式(31)，计算凯撒窗滤波器的形状参数βk；
[0122]
第三步，当原型滤波器阶数的允许设计范围较大时，采用最优化设计方法，根据公式(35)计算得出原型滤波器的阶数nf；当原型滤波器阶数的允许设计范围受限时，采用次优化设计方法，根据公式(37)计算得出原型滤波器的阶数nf′
。
[0123]
第四步，根据参数β和nf设计升余弦函数的原型分析滤波器；根据参数βk和nf，结合式(30)设计凯撒窗原型综合滤波器。
[0124]
实施例：
[0125]
按照上述步骤即可设计出性能良好的原型分析及综合滤波器，例如在信道个数为64，过采样倍数为2的滤波器组框架下，采用最优化方法设计原型分析及综合滤波器阶数均为2048。其中分析滤波器利用升余弦函数实现，滚降因子为0.65，综合滤波器采用理想滤波器加凯撒窗的方法实现，窗函数的形状参数βk＝12.2653。仿真得到原型滤波器的频率响应如图图5所示，其中h0(ω)和f0(ω)分别表示原型分析滤波器和原型综合滤波器，表示分析滤波器经过抽取和内插后的相邻频谱镜像，其间隔为2个子带宽度。图中还标注了在分析滤波器的镜像通带频率处，升余弦滤波器和所设计的原型综合滤波器的阻带衰减幅度，可以看出，本发明所采用的综合滤波器在镜像抑制方面具有明显优势。
[0126]
仿真验证：
[0127]
下面通过仿真验证改进原型滤波器的完全重建性能，部分仿真参数如表所示。
[0128]
表1仿真参数
[0129][0130]
计算分析滤波器组与综合滤波器组的联合频率响应，可得不同子带滤波器的频域分布如图6所示，为便于观察，图中仅画出了8个子带频率范围的滤波器组分布情况，可以看出，相邻滤波器频率响应交点的幅度为-3db，在该频率点处满足功率互补条件。
[0131]
下面分别对传统根升余弦原型滤波器和改进的原型滤波器在全频带内的功率互补条件进行验证，两种方法中的滤波器阶数均为2048，图7和图8分别表示使用两种方法对各个子带滤波器的功率谱相加后的结果，由图中可以看出，传统方法功率合成后误差在10-3
数量级，而改进方法的合成误差在10-5
数量级，因此改进方法在功率互补性能方面更加优越。
[0132]
下面使用单脉冲信号考察所设计滤波器组的总处理时延和重建性能，如图9所示为输入单脉冲信号和输出信号，可以看出，改进的滤波器组能够对单脉冲信号进行无失真传输，系统的总处理时延为分析与综合滤波器组的总时延，即2048个采样时刻。
[0133]
为了验证改进方法在输入宽带随机信号情况下的重建性能，将功率为0db的高斯白噪声作为输入信号，图10表示输入信号与重建信号的时域波形，可以发现，本发明所采用的滤波器组结构在宽带信号无失真传输方面仍具有较好的性能。
[0134]
对上图中的重建信号和输入信号做差即可得到重建误差，图11和图12分别表示采用传统根升余弦滤波器组和改进滤波器组的重建误差，两种方法的滤波器设计阶数均为2048。为了便于观察，图中截取了0～256个采样数据的重建误差，可以发现，传统方法对宽带信号的重建误差为10-3
数量级，而改进方法的重加误差为10-5
数量级。这是因为传统方法中仅用两个根升余弦滤波器作为原型滤波器，对信号抽取后的频谱镜像抑制能力较弱，而改进方法中的综合滤波器具有较好的带外抑制能力，降低了混叠干扰。
[0135]
综合以上分析，本发明所采用的非最大抽取dft滤波器组，通过设置2倍过采样并重新设计原型分析及综合滤波器，对输入信号具有良好的重建性能。当信道个数为64时，采用2048阶原型滤波器，所组成的高效信道化结构对0db功率高斯白噪声信号的重建误差在3＝10-5
以内，对信号后续处理的影响可以忽略不计。
[0136]
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种基于dft滤波器组结构的原型滤波器，其特征在于，原型综合滤波器采用凯撒窗滤波器，原型分析滤波器采用升余弦函数滤波器。2.一种如权利要求1所述的原型滤波器的设计方法，其特征在于，包括如下步骤：第一步，根据需要设置以下参数的值：滤波器组的信道个数k，原型分析滤波器滚降因子β和原型综合滤波器阻带最小衰减a
s
；第二步，根据如下公式，计算原型综合滤波器的形状参数β
k
：β
k
＝0.1102(a
s-8.7)；第三步，设计原型分析滤波器和原型综合滤波器的阶数n
f
；第四步，根据参数β和n
f
设计升余弦函数的原型分析滤波器；根据参数β
k
和n
f
设计原型综合滤波器w(n)：其中，n表示滤波器的采样点；i0(
·
)是第一类零阶变型贝塞尔函数。3.如权利要求2所述的设计方法，其特征在于，当滤波器阶数的允许设计范围超过设定值时，采用最优化设计方法，根据如下公式计算得出阶数n
f
：4.如权利要求2所述的设计方法，其特征在于，当原型滤波器阶数的允许设计范围受限时，采用次优化设计方法，根据如下公式计算得出阶数n
f
：

技术总结
本发明基于非最大抽取DFT滤波器组的框架，提供了一种原型分析及综合滤波器及设计方法，该滤波器增大滤波器的阻带抑制，提高完全重建性能，预期比传统的采用两个根升余弦滤波器作为原型滤波器的方法在重建误差方面降低约2个数量级；设计方法分析了滤波器参数的约束优化问题，并分别设计了最优化设计方法和次优化设计方法，两种方法实现简单且性能良好。两种方法实现简单且性能良好。两种方法实现简单且性能良好。


技术研发人员：滕飞 焦义文 杨文革 马宏 吴涛 李超 陈雨迪 卢志伟 高泽夫 李雪健 王育新
受保护的技术使用者：中国人民解放军战略支援部队航天工程大学
技术研发日：2023.03.21
技术公布日：2023/7/7