测试代码:
%%****bin_bp.m*******************%% %% 使用双线性变换法设计带通滤波器 %% 2018年6月13日 16:30:34 %% author:Alimy close all; clear; clc; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %代码正文 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %给定数字滤波器指标 f_sl = 150 ; %阻带下限频率(Hz) f_1 = 200 ; %通带下限频率(Hz) f_3 = 500 ; %通带上限频率(Hz) f_sh = 600 ; %阻带上限频率(Hz) R_p = 0.5 ; %通带允许的最大衰减 R_s = 40 ; %阻带允许的最小衰减 f_s = 2000 ; %采样频率 T_s = 1 / f_s ; %采样间隔 %1.将数字带通滤波器的频率参数变换为归一化的数字角频率参数 omega_sl = 2 * pi * f_sl / f_s; %阻带下限频率 omega_1 = 2 * pi * f_1 / f_s; %通带下限频率 omgea_3 = 2 * pi * f_3 / f_s; %通带上限频率 omega_sh = 2 * pi * f_sh / f_s; %阻带上限频率 %2.预畸变处理,将归一化数字角频率参数变换成模拟带通滤波器的角频率参数 C = 2*f_s ; Omega_sl &#