双线性映射matlab,[Matlab]双线性变换法设计数字带通滤波器

测试代码：

%%****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 = C * tan( o