切比雪夫I型滤波器频域特性
filter_design_cbi_001.png
-
时，
在1和
之间等幅振荡
-
大于
-
的值为：
切比雪夫I型滤波器设计步骤
步骤一：由通带截频
确定
步骤二：由阻带衰减As确定
:
$\varepsilon = \sqrt{10^{0.1Ap - 1}}$
步骤三：由通带、阻带指标确定N :
步骤四：由模方%7C%5E2)
求模拟低通滤波器的极点
filter_design_cbi_002.png
步骤5:由极点确定系统函数H(s)
,
filter_design_cbi_003.png
切比雪夫II型滤波器频域特性
filter_design_cbi_004.png
-
时，
- 对任意N,
,|H(0)| = 1
- 在通带
切比雪夫II型滤波器设计步骤
步骤一：由通带截频确定
步骤二：由阻带衰减As确定 :
$\varepsilon=\frac{1}{\sqrt(10^{0.1A_s} - 1)}$
步骤三：由通带、阻带指标确定N :
步骤四：由模方求模拟低通滤波器的极点
步骤5:由极点确定系统函数H(s)
使用matlab求模拟低通滤波器的系统函数以及频率响应
- 设计一个满足下列指标CBI型模拟低通滤波器fp=1KHz , fs=2KHz,Ap ≤1dB, As ≥40dB,求出其系统函数，并画出其频率响应特性图
Wp=2*pi*1000;%将数字频率转换成模拟频率
Ws=2*pi*2000;%将数字频率转换成模拟频率
Ap=1;%最大通带衰耗
As=40; % %最小阻带衰耗
%求阶数和Wc
[N,Wc]=cheb1ord(Wp,Ws,Ap,As,'s'); %Computer analog filter order
fprintf('Order of the filter=%.0f\n',N)
%求分之和分母系数
[num,den] = cheby1(N,Ap,Wc,'s'); %Compute AF coefficients
disp('Numerator polynomial');
fprintf('%.4e\n',num);
disp('Denominator polynomial');
fprintf('%.4e\n',den);
% 求Wp和Ws两点对应的频率响应
omega=[Wp Ws];
h = freqs(num,den,omega); %Compute Ap and As of the AF
ap = -20*log10(abs(h(1)));% 将通带频率响应转换成衰耗值，
as = -20*log10(abs(h(2)));% 将阻带频率响应转换成衰耗值，
hp = abs(h(1));
hs = abs(h(2));
gain_p = 20*log10(abs(h(1)));
gain_s = 20*log10(abs(h(2)));
fprintf('Ap= %.4f\n',-20*log10(abs(h(1))));
fprintf('As= %.4f\n',-20*log10(abs(h(2))));
%求0～6000pi之间的频率响应
omega = [0: 200: 3000*2*pi];
h = freqs(num,den,omega); %Compute the frequency response of the AF
subplot(211);
plot(omega/(2*pi),abs(h),'LineWidth',3);%画出频率响应曲线
hold on;%在画完函数虚线之后保持曲线图
plot([0 Wp/(2*pi)],[hp hp],'r--'); %画两个虚线
plot([Wp/(2*pi) Wp/(2*pi)],[0 hp],'r--');
plot([0 Ws/(2*pi)],[hs hs],'r--'); %画两个虚线
plot([Ws/(2*pi) Ws/(2*pi)],[0 hs],'r--');
hold off;
xlabel('Frequency in Hz');
ylabel('H(jw)');
gain=20*log10(abs(h));
subplot(212);
plot(omega/(2*pi),gain,'LineWidth',3);
hold on;%在画完函数虚线之后保持曲线图
plot([0 Wp/(2*pi)],[gain_p gain_p],'r--'); %画两个虚线
plot([Wp/(2*pi) Wp/(2*pi)],[0 gain_p],'r--');
plot([0 Ws/(2*pi)],[gain_s gain_s],'r--'); %画两个虚线
plot([Ws/(2*pi) Ws/(2*pi)],[0 gain_s],'r--');
hold off;
xlabel('Frequency in Hz');
ylabel('Gain in dB');
-
得到其图像如下
filter_design_cbi_005.png
-
其系数如下：
Order of the filter=5 %5阶滤波器
Numerator polynomial
0.0000e+00
0.0000e+00
0.0000e+00
0.0000e+00
0.0000e+00
1.2028e+18
Denominator polynomial
1.0000e+00
5.8862e+03
6.6672e+07
2.4170e+11
9.0479e+14
1.2028e+18
Ap= 1.0000
As= 45.3060
- 使用matlab分析其零极点分布情况
num=[1.2028e+18];
den=[1.0000e+00 5.8862e+03 6.6672e+07 2.4170e+11 9.0479e+14 1.2028e+18];
subplot(211);
sys=tf(num,den);%求出系统的零极点
pzmap(sys,'r');%画出系统零极点
grid on
subplot(212);
omega = [0: 200: 3000*2*pi];
H=freqs(num,den,omega);%得到频率响应
plot(omega/(2*pi),abs(H),'LineWidth',3);%画出频率响应
xlabel('Frequency Hz');
title('Magnitude Respone');
filter_design_cbi_006.png
- 由图可知极点的分布都在S平面的左半平面，说明是系统稳定
- N=5为奇数，一共有5个极点，和上面理论分析互相匹配
网友评论