以下是我对照 森荣二的《LC滤波器设计与制作》

这本书的电容耦合滤波式带通滤波器章节编写的一个Matlab程序。

电容耦合滤波式带通滤波器适合设计高Q的窄带滤波器，N阶谐振器耦合式BPF由N个谐振器和N－1个耦合元件K构成。

n=3; % the order of the butterworth LPF
for i=1:n
g(i)=2*sin((2*i-1)*pi/(2*n)); % butterworth LPF's g_i computation
end
disp('the normilized value of Butterworth LPF is: ')
g
BW=5e6; % bandwidth,unit Hz
F_ss=100e6; % unit Hz and means "suan shu ping jun "
Fl=F_ss-BW/2;
Fh=F_ss+BW/2;
F_jh=sqrt(Fl*Fh); % unit Hz and means "ji he ping jun "
for i=1:n-1
k(1:n-1)=1/(sqrt(g(i)*g(i+1)));
K(1:n-1)=BW*k(i)/F_jh;
end
K
L_resonator=10e-9; % unit is F
Z1=2*pi*F_jh^2*L_resonator*g(1)/BW
Z2=2*pi*F_jh^2*L_resonator*g(n)/BW % here Z2 maybe is not equal to Z1
C_resonator=1/[(2*pi*F_jh)^2*L_resonator];
L=input('Please input the L of the prototype filter:\n') % xu yao shou dong shu ru de L, for example, L=100e-9 H
C=input('Please input the C of the prototype filter:\n') % xu yao shou dong shu ru de C, for example, C=253.461e-12F
for i=1:n-1
C_couple(1:n-1)=K(1:n-1)*C_resonator;
C_new(1)=C_resonator-C_couple(1);
C_new(n)=C_resonator-C_couple(n-1);
C_new(2:n-1)=C_resonator-C_couple(n-2)-C_couple(n-1); % here i am not sure
end
C_new
K=50/Z1;
L=L*K % unit H
C=C_new/K % unit F

由于编程很少，里面的数组可能编的不规范，不过对应3阶的计算是正确的，还未验证其通用性，如果也有朋友看 森荣二的《LC滤波器设计与制作》这本书，欢迎讨论。

按照标准做出来的一定是不对称的,因为呢,低端的零点比较近,高端零点在无穷远;当你的中心频率比较高的时候,高低端看起来就比较对称了