IIR巴特沃斯低通滤波器源程序

function [ output_args ] = Untitled（ input_args ）
%UNtitleD Summary of this function goes here
%   Detailed explanation goes here
wp=0.35*pi;ws=0.4*pi;
Rp=1;Rs=60;                         %% 通带\阻带波动（dB）
Fs=1;Ts=1/Fs;
wp1=2/Ts*tan（wp/2）;                 %将模拟指标转换成数字指标
ws1=2/Ts*tan（ws/2）; 
[NWn]=buttord（wp1ws1RpRs‘s‘）;  %选择滤波器的最小阶数
[ZPK]=buttap（N）;                  %求N阶的butterworth模拟滤波器的零、极点和增益因子
[BapAap]=zp2tf（ZPK）;             %根据系统的零、极点、增益还原出H（z）的分子分母系数
[ba]=lp2lp（