基于ADS\HFSS\Matlab的平面带通发夹滤波器设计

秋水离月

已于 2024-12-03 23:49:14 修改

阅读量2.1k

点赞数 24

文章标签：平面

于 2024-11-14 03:02:00 首次发布

本文链接：https://blog.csdn.net/likaifeng12138/article/details/143755940

版权

写在最前面，前段时间完成了一个小项目，实现13.5~14GHz的带通滤波器，算是第一个微波项目，在老师、学长、各种大神的指点下，最终完成了指标，想来这段时间收获良多，想到老师教导的，有能力的话要尽可能帮助别人，以及网络上很少面向于小白的教程，于是决定把自己做这个小项目的经过，一些关键步骤，还有软件操作和公式讲解都发出来，供大家指正借鉴。由于我目前大三，刚学习电磁场、微波工程等等，技术实在有限，同时这样是我第一次发帖子，只能尽可能的将自己确认过正确的内容，步骤以及一些经验写出来，供大家借鉴。也请各位大神不吝赐教。

微带发夹型滤波器是一种经典半波长耦合微带线滤波器，使用半波长耦合谐振器折叠来实现的。该设计继承了原有半波长耦合微带线滤波器的特性，通过开路式对称耦合微带单元的级联构造，并保持耦合微带线的电长度为90°。巴拉巴拉~

总之这个东西就是靠一个个的谐振单元，让所需要通过的频率通过，将不要的频率分量去除。

微带发夹滤波器的性能主要由谐振器臂长L、微带线宽度W，谐振器之间的间距s以及抽头距离谐振器的位置t等参数定．其中，谐振器臂长L决定了滤波器中心频率的大小，其长度计算公式

其中: $\lambda_{0}$ 为滤波器的中心频率在空间传播的波长; $\lambda_{g}$ 为滤波器的中心频率在介质与自由空间传播的相对波长； $\varepsilon _{eff}$ 为有效介电常数，其表达式为

谐振器之间的间距s与抽头位置t对滤波器的性能有着重要的影响，其中谐振器的间距与耦合系数有关，它影响着滤波器的带宽值。

其中:B为相对带宽,i为谐振器的序号,n为滤波器的阶数,gi为滤波器低通原型中第i个归一化元件值.

本项目中归一化参数选择的是0.5dB纹波Chebyshev原型，通过查询0.5dB纹波Chebyshev归一化表可得归一化参数，再由归一化参数计算计算奇模阻抗和偶模阻抗，再通过ADS的计算工具LineCalc计算线宽，耦合器间距，耦合系数，等等，是不是感觉很麻烦，不用担心，我已经通过Matlab写好了计算的代码，你只需要修改参数就可以很简单的计算出上面的数据，如果有人想看的话，以后可以和大家一起推导公式

我这边设计的是

通带：13.5~14GHz；

带内插损IL：＜4dB；

带内反射系数 RL：＞-10dB；

上边带抑制：14.25GHz以上至少抑制13dB、15.5-17.5 GHz 抑制25dB以上；

下边带抑制：13.25GHz 以下至少抑制 13dB、12GHz 以下至少抑制35dB、9 -11GHz 以下至少抑制 45dB；

板材RO4350B，相对介质 $\varepsilon _{r}$ =3.66（对应代码中er）,板材厚度h=0.254mm,中心频率f=13.75G

直接复制到Matlab中就可以运行


% 本程序用于确定平面微带发夹滤波器的
%     阶数n，耦合系数k(j,j+1)，间距s(j,j+1),偶模阻抗Zoe,奇模阻抗Zoo,抽头位置t
% 配合ADS的LineCale功能，可以计算出线宽，线距，线长


%% 计算ws以在低通模型中查表计算阶数
%设置阶数
n=5;
%设置中心频率（GHz）
f=13.75;
%上边频和下边频(GHz)
f1=13.5;
f2=14;
%带外抑制频率,频率有显著下降的点,单位GHz
dw=14.25;



%相对带宽
FBW=(f2-f1)/f;
%带宽
BW=f2-f1;

%光速,频率单位是G，波长单位mm,所以光速设置为300
v0=300;
%线宽
w=0.5;
%板材厚度
h=0.254;
%相对介质
er=3.66;
%有效介质常数
eeff=(er+1)/2+((er-1)/2)/(sqrt(1+(12*h/w)));


%特征阻抗，在ADS计算工具LineCale中输入奇模阻抗，和偶模阻抗之后计算得来
Zr=50.049451;


%计算1/4波长(cm)
L=v0/f/(sqrt(eeff))/4;


%确定阶数所用的数据，w'-w'1=ws
ws=1/FBW*(dw/f-f/dw);



%% 依次输入归一化参数

g(1)=1;%g0
g(2)=1.1468;
g(3)=1.3712;
g(4)=1.9750;
g(5)=1.3712;
g(6)=1.1468;
g(7)=1;


%% 计算奇模阻抗和偶模阻抗

%输入端口阻抗
z0=50;
%应该是z0*j(c,c+1)的，但是为了方便，将z0*j(j,j+1)写为z0j(c+1)
for c=0:n
    if c==0
        z0j(c+1)=sqrt(pi*FBW/2/g(1)/g(2));
    elseif c==n
        z0j(c+1)=sqrt(pi*FBW/2/g(1)/g(2));
    else
        z0j(c+1)=pi*FBW/(sqrt(g(c+1)*g(c+2)))/2;
    end

end


%偶模阻抗Zoe，奇模阻抗zoo
for c=0:n
    if c==0
        zoe(c+1)=z0*(1+(z0j(c+1)))+((z0j(c+1))^2);
        zoo(c+1)=z0*(1-(z0j(c+1)))+((z0j(c+1))^2);
    elseif c==n
        zoe(c+1)=z0*(1+(z0j(c+1)))+((z0j(c+1))^2);
        zoo(c+1)=z0*(1-(z0j(c+1)))+((z0j(c+1))^2);
    else
        zoe(c+1)=z0*(1+(z0j(c+1)))+((z0j(c+1))^2);
        zoo(c+1)=z0*(1-(z0j(c+1)))+((z0j(c+1))^2);
    end

end

%% 计算耦合系数,输入使用馈电线的情况

for c=1:n-1

    k(c)=BW/f/(sqrt(g(c+1)*g(c+2)));

end

% %% 计算耦合系数,输入使用耦合线的情况
% 
% for c=0:n-1
% 
%     k(c+1)=BW/f/(sqrt(g(c+1)*g(c+2)));
% 
% end

%外部品质因数
Qe=f/BW*g(1);

%计算输入输出抽头位置
t=2*L/pi*asin(sqrt(pi*Zr/2/Qe/z0));

下面就教大家如何使用这个代码，所有需要我们手动输入的参数我全部在后面加了一个符号‘%’，仔细一点的话应该没有问题

1.计算阶数

先随便写一个阶数，设置中心频率，上下边频，相对介质，dw(带外抑制频率,频率有显著下降的点)，然后其他的不用管，计算，在旁边工作区里找到ws