测试代码:
%%****bin_hp.m*******************%% %% 使用双线性变换法设计高通滤波器 %% 2018年6月13日 14:27:37 %% author:Alimy close all; clear; clc; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %代码正文 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %给定数字滤波器指标 f_p = 200 ; %阻带上截止频率 f_st = 210 ; %通带下截止频率 R_p = 3 ; %通带允许的最大衰减 R_st = 30 ; %阻带允许的最小衰减 f_s = 1000 ; %采样频率 T_s = 1 / f_s ; %采样间隔 %1.将数字高通滤波器的频率参数变换为归一化的数字角频率参数 omega_p = 2 * pi * f_p / f_s; omega_st = 2 * pi * f_st / f_s; %2.预畸变处理,将归一化数字角频率参数变换成模拟高通滤波器的角频率参数 C = 2 / T_s ; Omega_p = C * tan( omega_p / 2 ); Omega_st = C * tan( omega_st / 2 ); %3.对模拟高通滤波器的角频率参数做归一化处理 lamda_p = Omega_p / Omega_p; lamda_st = Omega_st / Omega_p; %4.设计归一化模拟滤波器,得到归一化模拟高通滤波器的转移函数 [ N , Wn ] = buttord( Omega_p