一种声波滤波器带通域直接综合设计方法与流程

1.本发明涉及滤波器综合设计技术，特别是涉及一种声波滤波器带通域直接综合设计方法，属于基本电气元件的技术领域。


背景技术：

2.通信系统的快速发展对滤波器的尺寸和性能提出了越来越高的要求。薄膜体声波滤波器因其具有体积小、插损低、矩形系数高的优势在各通信系统中得到广泛应用，如何高效快捷地设计薄膜体声滤波器成为关键。
3.目前，针对传统腔体、介质、微带等滤波器综合方法的研究已经相当成熟，但针对声波滤波器的综合研究还较少。alfredgimenez等人于2018年在文章“general synthesis methodology for the design of acoustic wave ladder filters and duplexers”(声波阶梯滤波器和双工器的综合设计方法)中提出一种低通到带通的声波滤波器综合方法，其主要设计思路是将串联谐振器和并联谐振器等效为非谐振节点，并给出对应的等效电路，基于切比雪夫综合计算等效电路参数再将其转换为低通bvd模型电路参数，最后再转换为带通bvd(butterworth
‑
van dyke，巴特沃斯范戴克)模型电路参数，然而，该方法中非谐振节点对应的等效电路中的电抗为非频变元件，因此其综合结果只在中心频率f0附近准确度较高，远带外差异较大。iuliiaevdokimova等人2017年在文章“synthesis of ladder
‑
type acoustic filters in theband
‑
pass domain”(梯形声波滤波器的带通域综合设计方法)中提出一种带通域直接综合声波滤波器的方法，其核心是根据滤波器特性构造适合带通域直接综合的切比雪夫函数，这种方法避免了低通域综合采用非频变电抗的缺点，但文中的方法仅适合于对称结构的声波滤波器综合。
4.综上，本发明旨在提出一种带通域直接综合设计声波滤波器的方法以适应任意阶数对称和非对称拓扑结构的电路参数设计需求。


技术实现要素：

5.本发明的发明目的是针对上述背景技术的不足，提供一种声波滤波器带通域直接综合设计方法，先将切比雪夫多项式转换到带通域推导滤波器输入阻抗和导纳，再根据推导滤波器输入阻抗和导纳逐阶提取谐振器bvd模型参数，实现带通域直接综合设计声波滤波器电路参数的发明目的，解决先低通域设计电路参数后转换为带通域电路参数的设计方法只在窄带内较准确的技术问题。
6.本发明为实现上述发明目的采用如下技术方案：
7.一种声波滤波器综合方法，具体包括如下步骤：
8.(1)选定滤波器拓扑，根据滤波器指标，确定滤波器阶数n，通带频率范围[ω
l
,ω
h
]，带外零点位置ω
n
(n＝1、2、3、
…
、n,n为滤波器的阶数)，以及回波损耗rl；
[0009]
(2)将低通域综合的切比雪夫多项式按照映射
关系转换到带通域，得到带通域切比雪夫多项式：
[0010][0011]
(3)根据带通滤波器零点位置，由步骤(2)得到的带通域切比雪夫多项式g(ω)计算传输多项式p(ω)、反射多项式f(ω)、以及分母多项式e(ω)；进一步计算声波滤波器的输入阻抗z
in
和输入导纳y
in
；
[0012]
(4)根据z
in
和y
in
直接逐级计算声波滤波器对应的串联谐振器和并联谐振器bvd模型参数，保证了本发明提出的方法可以适用于任意阶对称和非对称的声波滤波器综合。
[0013]
本发明采用上述技术方案，具有以下有益效果：发明提出的方法是在带通域直接综合，避免在低通综合过程中采用非频变电抗导致综合只在窄带内较准确的问题；且本发明是根据传输多项式、反射多项式、以及分母多项式计算出滤波器输入阻抗，然后逐阶计算各级谐振器bvd模型参数，这样计算的优点是该方法可以适合于各种对称的或非对称的滤波器拓扑。本发明提出的方法计算过程简洁明了，计算结果准确度高，且适用范围广。
附图说明
[0014]
图1为本发明带通域直接综合设计声波滤波器的流程图。
[0015]
图2为本发明实施例中七阶声波滤波器等效电路图。
[0016]
图3为本发明实施例中八阶声波滤波器等效电路图。
[0017]
图4为本发明实施例中bvd电路示意图。
[0018]
图5为本发明实施例一的声波滤波器频率响应曲线。
[0019]
图6为本发明实施例二的声波滤波器频率响应曲线。
具体实施方式
[0020]
下面结合附图对发明的声波滤波器带通域直接综合的方法进行详细说明。
[0021]
本发明所使用的变量符号并不限定于本发明之中所采用的符号。如果采用其它变量符号来替代本发明声波滤波器带通域直接综合方法中所采用的变量符号，皆可认为替换为其它变量符号的技术方案落入本发明的保护范围。
[0022]
本发明的说明书和附图示例性地描述了本发明的具体实施方式。本发明只考虑无损耗的情况，图1为声波滤波器带通域直接综合方法的流程图，具体包括如下6个步骤。
[0023]
步骤1)根据滤波器指标要求选定滤波器阶数n及拓扑结构
[0024]
图2为一种对称拓扑结构示意图，图3为一种非对称拓扑结构示意图，明确滤波器通带频率范围[ω
l
,ω
h
]，ω
l
和ω
h
分别为声波滤波器的上边带和下边带，滤波器中心频率带外零点位置ω
n
(n＝1、2、3、
…
、n,n为滤波器的阶数)，通带回波rl。
[0025]
步骤2)将切比雪夫函数由低通域变换到带通域
[0026]
切比雪夫函数为：
[0027][0028]
其中，式中，ω
n
(n＝1、2、3、
…
、n)为带通滤波器各阶零点转换到
低通对应的零点，且函数x
n
(ω)满足条件：
[0029]
1)当ω＝ω
n
时，x
n
(ω)＝
±
∞，其中，ω
n
指带通率波器零点变换到低通对应的零点；
[0030]
2)在通带有x
n
(
±
ω
h
)＝
‑
x
n
(
±
ω
l
)＝
‑
1，其中，ω
h
和ω
l
为带通滤波器通带边界变换到低通的值；
[0031]
3)当频率满足ω＝
±
ω
h
和ω＝
±
ω
l
时，1≥x
n
(ω)≥
‑
1。
[0032]
根据低通域与带通域间的映射关系，将切比雪夫函数转换得到带通域切比雪夫多项式，有：
[0033][0034]
步骤3)计算传输多项式、反射多项式、分母多项式
[0035]
根据带通滤波器零点位置，由步骤(2)修正得到的带通域切比雪夫多项式计算出传输多项式p(ω)、反射多项式f(ω)，并运用交替极点法构造出分母多项式e(ω)。
[0036]
步骤4)计算滤波器输入阻抗z
in
[0037]
根据传输多项式p(ω)、反射多项式f(ω)、及分母多项式e(ω)，计算滤波器的输入阻抗z
in
(ω)和输入导纳y
in
(ω)：
[0038][0039]
y
in
(ω)＝1/z
in
(ω)。
[0040]
其中，s
11
(ω)为滤波器的反射系数，θ
11
代表s
11
的相角。
[0041]
步骤5)根据输入阻抗z
in
(ω)和输入导纳y
in
(ω)递归计算bvd模型参数
[0042]5‑
a)首先提取第一个串联谐振器bvd模型的相关参数：根据第一个串联谐振器在ω1的位置有传输零点，则可以将滤波器输入阻抗z
in
(ω)表示为：
[0043][0044]
其中，z1(ω)为从第二个谐振器输入端位置等效的输入阻抗，即第二谐振器的输入阻抗。
[0045]
根据图4，单个串联谐振器resseries1的bvd模型等效阻抗可以由l
a1
、c
a1
、c
01
表示为：则有：
[0046][0047]
则根据串联谐振器1在ω1的位置有传输零点，可以得到l
a1
＝1/k1，c
a1
＝ω1/k1。对于电路参数c
01
，根据阻抗的保留部分z1(ω)在频率ω2处有传输零点，且该零点位于滤波器的右带外，则c
01
为z
in
(ω)在ω2处的虚部，即c
01
＝imag(z
in
(ω2))。l
a1
、c
a1
为bvd模型中与机械相关的动态电感、动态电容，c
01
为bvd模型中的静态电容。
[0048]5‑
b)再提取第一个并联谐振器bvd模型的相关参数：根据并联谐振器在ω2的位置有传输零点，则可以将从第二个谐振器的输入导纳y1(ω)表示为：
[0049][0050]
其中，y2(ω)为从第二个谐振器电路元件c
02
输入端位置等效的输入导纳，即第三个谐振器的输入导纳。
[0051]
根据图4，l
a2
和c
a2
串联电路的导纳为则有：
[0052][0053]
根据第一个并联谐振器在ω2处有传输零点，可以得到l
a2
＝1/k2，c
a2
＝ω2/k2。
[0054]
对于电路参数c
02
，根据阻抗的保留部分y2(ω)在ω3处有传输零点，且该零点位于滤波器的左带外。则c
02
为y2(ω)在ω3处的虚部，即c
02
＝imag(y2(ω3))。
[0055]
至此，对于阶数大于2阶的滤波器，重复第5)步的a)、b)步骤，可以实现任意阶数及任意拓扑的串联谐振器和并联谐振器的bvd模型参数的依次提取。
[0056]
同样地，当滤波器第一个谐振器为并联谐振器，第二个谐振器为并联谐振器时，调换步骤b)和步骤a)，逐阶递归计算每个谐振器的bvd模型参数。
[0057]
步骤6)根据提取的bvd参数仿真得声波滤波器频率响应曲线
[0058]
以下通过具体的声波滤波器综合实例进一步说明本发明。
[0059]
实施例1：
[0060]
采用图2所示拓扑结构，滤波器阶数n＝7，通带频率范围取[2402mhz，2482mhz]，带外零点位置分别为[2510mhz，2350mhz，2525mhz，2359mhz，2525mhz，2350mhz，2510mhz],回波损耗il＝18db。
[0061]
根据带通域切比雪夫多项式计算传输多项式p(ω)、反射多项式f(ω)、以及分母多项式e(ω)有：
[0062]
f(ω)＝ω7‑
0.508ω6 1.78ω5‑
0.81jω4 0.91ω3‑
0.34jω2 0.1163ω
‑
0.0236j，
[0063]
p(ω)＝ω7‑
0.798ω6 11.8ω5‑
11.04jω4 46.8322ω3‑
49.29jω2 61.75ω
‑
71.524j，
[0064]
e(ω)＝ω7 (1.958
‑
0.508j)ω6 (3.699
‑
1.011j)ω5 (4.04
‑
1.82j)ω4 (3.569
‑
1.921j)ω3 (2.1149
‑
1.5558j)ω2 (0.848
‑
0.803j)ω (0.1602
‑
j0.2385)。
[0065]
进一步，根据f(ω)、p(ω)、e(ω)计算滤波器输入阻抗再计算各串并联谐振器bvd模型参数：
[0066][0067]
[0068]
最后，根据表中的bvd模型参数仿真的滤波器频率响应曲线如图5所示。
[0069]
实施例2：
[0070]
采用图3所示拓扑结构，滤波器阶数n＝8，通带频率范围取[2496mhz，2690mhz]，带外零点位置分别为[2450mhz，2750mhz，2460mhz，2740mhz，2465mhz，2745mhz，2455mhz，2755mhz],回波损耗il＝18db。
[0071]
根据带通域切比雪夫多项式计算传输多项式p(ω)、反射多项式f(ω)、以及分母多项式e(ω)有：
[0072]
f(ω)＝ω9 0.724ω8‑
2.365jω7‑
1.6387jω6 1.8781ω5 1.2127jω4‑
0.5531ω3‑
0.3114jω2 0.042ω 0.0152，
[0073]
p(ω)＝ω9 0.6432ω8‑
9.2767jω7‑
6.7157jω6 32.3605ω5 25.8206jω4‑
50.2854jω3‑
43.4729jω2 29.3552ω 27.1307，
[0074]
e(ω)＝ω9 (1.0237
‑
0.037j)ω8 (0.9425 0.128j)ω7 (0.739
‑
0.265j)ω6 (0.326
‑
0.427j)ω5 (0.246
‑
0.426j)ω4 (0.658
‑
0.268j)ω3 (0.868
‑
0.141j)ω2 (1.012
‑
0.019j)ω (0.970
‑
0.066j)。
[0075]
进一步，根据f(ω)、p(ω)、e(ω)计算滤波器输入阻抗再计算各串并联谐振器bvd模型参数：
[0076][0077]
最后，根据表中的bvd模型参数仿真的滤波器频率响应曲线如图6所示。
[0078]
上述仅为本发明的优选实施例，并不对本发明起到任何限制作用。任何所属技术领域的技术人员，在不脱离本发明的发明宗旨，对本发明揭露的技术方案和技术内容做任何形式的等同替换或修改等变动均落入本发明的保护范围之内。