[Matlab]双线性变换法设计数字低通滤波器

测试代码：

%%****bin_lp.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