fp=1000; fst=1500; Rp=1; As=15; Ap=Rp;
epc=sqrt(10^(Ap/10)-1); Wp=2*pi*fp; Wc=Wp;
Wst=2*pi*fst;
N=ceil(acosh(sqrt(10^(0.1*As)-1)/epc)/acosh(Wst/Wc)); [b a]=cheby1(N,Rp,Wc,'s'); [bz az]=impinvar(b,a,fs); w=[0:500]*pi/500; [H w]=freqz(bz,az); subplot 221; plot(w/pi,abs(H)); xlabel('\\omega(\\pi)');
ylabel('|H(e^j^\\omega)|(dB)'); grid on; subplot 222;
plot(w/pi,20*log10(abs(H))); xlabel('\\omega(\\pi)');
ylabel('|H(e^j^\\omega)|(dB)'); grid on; subplot 223;
plot(w/pi,angle(H)/pi); xlabel('\\omega(\\pi)'); ylabel('|H(e^j^\\omega)|'); grid on; subplot 224;
grd=grpdelay(bz,az,w); plot(w/pi,grd);
xlabel('\\omega(\\pi)'); ylabel('|H(e^j^\\omega)|'); grid on;
1.550|H(ej?)|(dB)1|H(ej?)|(dB)00.5?(?)100.5-500-10000.5?(?)110.51510|H(ej?)|0-0.5-100.5?(?)1|H(e)|j?5000.5?(?)1
（3）切比雪夫II型 fs=10000; fp=1000; fst=1500; Rp=1; As=15; Ap=Rp;
epc=sqrt(10^(Ap/10)-1); Wp=2*pi*fp; Wst=2*pi*fst; Wc=Wst;
N=ceil(acosh(sqrt(10^(0.1*As)-1)/epc)/acosh(Wst/Wp)); [b a]=cheby2(N,As,Wc,'s'); [bz az]=impinvar(b,a,fs) w=[0:500]*pi/500; [H w]=freqz(bz,az); subplot 221; plot(w/pi,abs(H)); xlabel('\\omega(\\pi)');
ylabel('|H(e^j^\\omega)|(dB)'); grid on; subplot 222;
plot(w/pi,20*log10(abs(H)));
xlabel('\\omega(\\pi)');
ylabel('|H(e^j^\\omega)|(dB)' ); grid on;subplot 223;
plot(w/pi,angle(H)/pi); xlabel('\\omega(\\pi)'); ylabel('|H(e^j^\\omega)|'); grid on; subplot 224;
grd=grpdelay(bz,az,w); plot(w/pi,grd);
xlabel('\\omega(\\pi)'); ylabel('|H(e^j^\\omega)|'); grid on;
1.5)B1(d|)j?(eH|0.5000.51?(?)10.5|)j?(e0H|-0.5-100.51?(?)
5)B0(d|)(eH|-5-1000.51?(?)5|)j?(e0H|-500.51?(?)j?
（4）椭圆模拟原型滤波器 fs=10000; fp=1000; fst=1500; Rp=1; As=15; Ap=Rp;
epc=sqrt(10^(Ap/10)-1); Wp=2*pi*fp; Wst=2*pi*fst; Wc=Wp;
A=10^(As/20);
k1=epc/(sqrt(A^2-1)); k=Wp/Wst;
N=ceil(ellipke(k)*ellipke(sqrt(1-k1^2))/ellipke(k1)/ellipke(sqrt(1-k^2)));
[bz az]=ellip(N,Rp,As,fp/fs*pi) w=[0:500]*pi/500; [H w]=freqz(bz,az); subplot 221; plot(w/pi,abs(H)); xlabel('\\omega(\\pi)');
ylabel('|H(e^j^\\omega)|(dB)'); grid on; subplot 222;
plot(w/pi,20*log10(abs(H))); xlabel('\\omega(\\pi)');
ylabel('|H(e^j^\\omega)|(dB)'); grid on; subplot 223;
plot(w/pi,angle(H)/pi); xlabel('\\omega(\\pi)'); ylabel('|H(e^j^\\omega)|'); grid on; subplot 224;
grd=grpdelay(bz,az,w); plot(w/pi,grd);
xlabel('\\omega(\\pi)'); ylabel('|H(e^j^\\omega)|'); grid on;
1.5|H(ej?)|(dB)50|H(ej?)|(dB)100.5-50000.5?(?)1-10000.5?(?)110.5|H(ej?)|15100-0.5-100.5?(?)1|H(e)|j?5000.5?(?)1
用脉冲响应不变法设计IIR数字滤波器的优缺点及适用范围： 优点：运算简单，容易实现。可利用模拟滤波器设计技术。
缺点：由于发生混叠效应是阻带增益降低，容易不满足最初要求。且不能用FFT技术。不能够设计高通和带阻滤波器。
适用范围：不能够设计高通和带阻滤波器。用于设计规格化选频滤波器，h(n)无限长,用于模拟滤波器数字化。低阶次更容易实现。
五、实验心得与体会
通过这次实验，我掌握了利用MATLAB实现脉冲响应不变法设计IIR数字滤波器的具体方法。加深了对其设计原理的理解和认识，在实验过程中使我对数字滤波器和模拟滤波器之间的技术指标转化有了深刻理解，并且由直观的实验结果和理论课知识相结合使我掌握了脉冲响应不变法设计IIR数字滤波器的整个过程以及这种方法的优缺点及适用范围。所以，这次实验使我受益匪浅。