南京工程学院通信工程学院数字信号处理第9章.doc

1.（1）n1=0;n2=10;a1=0.5; na1=n1:n2; x1=a1.^na1; stem(na1,x1,k); （2）n1=0;n2=10;a1=2; na1=n1:n2; x1=a1.^na1; stem(na1,x1,k); （3）n=0:20; x=sin(n*pi/5); stem(n,x); （4）n=0:20; x=cos(n*pi/10-pi/5); stem(n,x); 2.（1）a=[3,4,1];b=[1,1,0];n=32; hn=impz(b,a,n); stem(hn,k); （2）a=[5/2,6,10];b=[1,0,0];n=32; hn=impz(b,a,n); stem(hn,k); 3. n1=0:15;n01=0; h1=[(n1-n01)=0]; n2=n1;n02=10; h2=[(n2-n02)=0]; h=(7/8).^n1.*(h1-h2); x1=[(n1-n01)=0]; n2=n1;n03=5;x2=[(n2-n03)=0]; x=(x1-x2); y=conv(x,h); stem(y); 4. b=[2,16,44,56,32];a=[3,3,-15,18,-12]; [r p c]=residuez(b,a) N=20;n=0:N-1; x=r(1)*p(1).^n +r(2)*p(2).^n stem(n,x,k); title(用部分分式法求反变换x(n)); axis([0,N-1,1.1*min(x),1.1*max(x)]); 5.（1）b=[0,2,-1.6,-0.9];a=[1,-2.5,1.96,-0.48]; rz=roots(b);rp=roots(a); subplot(1,2,1),zplane(b,a);title(系统的零极点分布图); subplot(1,2,2),impz(b,a,20);title(系统的冲激响应);xlabel(n);ylabel(h(n)); 不是因果稳定系统 （2）b=[0,0,1,-1];a=[1,-0.9,-0.65,0.873]; rz=roots(b);rp=roots(a); subplot(1,2,1),zplane(b,a);title(系统的零极点分布图); subplot(1,2,2),impz(b,a,20);title(系统的冲激响应);xlabel(n);ylabel(h(n)); 因果稳定系统 6. b=[1,0,0];a=[1,-3/4,1/8]; n=(0:1000)*pi/1000; [h,w]=freqz(b,a,n); subplot(2,1,1),plot(n/pi,abs(h));grid axis([0,1,1.1*min(abs(h)),1.1*max(abs(h))]);ylabel(幅度); subplot(2,1,2),plot(n/pi,angle(h));grid axis([0,1,1.1*min(angle(h)),1.1*max(angle(h))]);ylabel(相位);xlabel(以pi为单位的频率); 7. n1=0:10;n2=0:18;N=20; x1=0.8.^n1;x2=0.6.^n2; y=cconv(x1,x2,N) 8. n=0:20;N=20; x=(0.9*exp(j*pi/3)).^n; y=fft(x,20) 9. fp=3400;Omgp=2*pi*fp; fs=4000;Omgs=2*pi*fs; Rp=3;As=40; [n,Omgc]=buttord(Omgp,Omgs,Rp,As,s); [z0,p0,k0]=buttap(n); b0=k0*real(poly(z0)); a0=real(poly(p0)); freqs(b0,a0); 10. T=1;wp=0.2*pi/T;ws=0.3*pi/T;rp=1;rs=15; [N,wc]=cheb1ord(wp,ws,rp,rs,s); [B,A]=cheby1(N,rp,wc,s);[Bz,Az]=impinvar(B,A); [H,w]=freqz(Bz,Az); plot(w/pi,abs(H));xlabel(频率/{\pi});ylabel(幅度); 11. T=1;Fs=1/T;wp=0.25*pi/T;ws=0.4*pi/T;rp=2;rs=20; Omgp=(2/T)*tan(wp/2);Omgs=(2/T)*tan(ws/2); [N,Omgc]=ellipord(Omgp,