资源描述：
《离散系统的频率响应分析和零点,极点分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、离散系统的频率响应分析和零点，极点分析电子信息工程四班邱锐煌200530580425一、实验目的加深对离散系统的频率响应分析和零点极点的分布概念理解。零点，极点分布图：>>num=[0.05280.07970.12950.12950.7970.0528];>>den=[1-1.81072.4947-1.88010.9537-0.2336];>>[z,p]=tf2zp(num,den);>>zplane(z,p);幅频响应：num=[0.05280.07970.12950.12950.7970.0528];>>den=
2、[1-1.81072.4947-1.88010.9537-0.2336];>>w=linspace(0,2*pi,512);>>h=freqz(num,den,w);>>plot(w/pi,abs(h));num=[0.05280.07970.12950.12950.7970.0528];>>den=[1-1.81072.4947-1.88010.9537-0.2336];freqz(num,den,128)部分分式展开计算：num=[0.05280.07970.12950.12950.7970.0528];>>de
3、n=[1-1.81072.4947-1.88010.9537-0.2336];>>[r,p,k]=residuez(num,den)r=p=0.8168-0.5527i0.2788+0.8973i0.8168+0.5527i0.2788-0.8973i-2.1440+0.8591i0.3811+0.6274i-2.1440-0.8591i0.3811-0.6274i2.93320.4910k=-0.2260高阶系统分解为2阶系统的串联：>>num=[0.05280.07970.12950.12950.7970.052
4、8];>>den=[1-1.81072.4947-1.88010.9537-0.2336];>>sos=zp2sos(z,p,k)sos=-0.2260-0.015101.0000-0.491001.00003.17404.61251.0000-0.76220.53891.0000-1.73153.23921.0000-0.55750.8829基于MATLAB的数字滤波器设计实验目的：加深对数字数字滤波器的常用指标和设计过程的理解。IIRwp1=0.45*pi;wp2=0.65*pi;Ap=1;ws1=0.3*pi;
5、ws2=0.75*pi;As=40;B=wp2-wp1;w0=sqrt(wp1*wp2);Ws1=(ws2.^2-w0.^2)./(B*ws1);Ws2=(ws2.^2-w0.^2)./(B*ws2);Ws=min(Ws1,Ws2);Ws=Ws/pi;[N,Wc]=buttord(1,Ws,Ap,As,'s');[num,den]=butter(N,Wc,'s');Ws1=tan(ws1/2);Ws2=tan(ws2/2);W0=sqrt(Ws1*Ws2);[numt,dent]=lp2bp(num,den,W0,B
6、);[numt,dent]=bilinear(numt,dent,0.5);w=linspace(0,2*pi,512);h=freqz(numt,dent,w);subplot(2,2,1)plot(w/pi,20*log10(abs(h)));title('IIR')grid;subplot(2,2,2)plot(w/pi,abs(h));title('幅频');gridsubplot(2,2,3)plot(w/pi,angle(h));title('相频');gridsubplot(2,2,4)stem(imp
7、z(numt,dent,80));title('冲激响应')gridFIＲFs1=0.3;Fp1=0.45;Fs2=0.65;Fp2=0.75;dp=0.10875;ds=0.01;f=[Fs1Fp1Fs2Fp2];a=[010];dev=[dsdpds];[M,fo,ao,w]=remezord(f,a,dev);h=remez(M,fo,ao,w);w=linspace(0,pi,1000);mag=freqz(h,[1],w);subplot(2,2,1)hd=plot(w/pi,20*log10(abs(ma
8、g)));title('FIR')grid;subplot(2,2,2)plot(w/pi,abs(mag));title('幅频');gridsubplot(2,2,3)plot(w/pi,angle(mag));title('相频');gridsubplot(2,2,4)stem(impz(h,80));title('冲激响应')grid实
