MATLAB仿真IIR带通滤波器

%巴特沃斯带通滤波器
clear all
n=0:180;fs=2000;%采样频率2000hz
k=0:fs-1;f=2*k/fs;
x=cos（0.1*n*2*pi）;%[（n-17）>=0]-[n-234>=0];
y=sin（0.015*n*2*pi）;
w=x+y;
figure（1）;
subplot（621）;stem（x‘.‘）;title（‘输入信号x（n）‘）;ylabel（‘x（n）‘）;
subplot（623）;stem（y‘.‘）;title（‘输入信号y（n）‘）;ylabel（‘y（n）‘）;
subplot（625）;stem（w‘.‘）;title（‘输入信号w（n）‘）;ylabel（‘w（n）‘）
X=fft（xfs）;
subplot（622）;plot（fabs（X））;title（‘输入信号x（n）幅频曲线‘）;ylabel（‘|FT[x（n）]|‘）;
Y=fft（yfs）;
W=fft（wfs）;
subplot（626）;plot（fabs（W））;title（‘输入信号w（n）幅频曲线‘）;ylabel（‘|FT[w（n）]|‘）;
subplot（624）;plot（fabs（Y））;title（‘输入信号y（n）幅频曲线‘）;ylabel（‘|FT[y（n）]|‘）;
Wp=[2*pi*160/fs2*pi*310/fs];Ws=[2*pi*120/