南京理工大学数字信号处理matlab上机完美版.pdf

1. 已知3 阶椭圆IIR 数字低通滤波器的性能指标为：通带截止频率 0.4 π，通带波纹为0.6dB，最小阻带衰减为32dB。设计一个6 阶全 通滤波器对其通带的群延时进行均衡。绘制低通滤波器和级联滤波器 的群延时。 %Q1_solution %ellip(N,Ap,Ast,Wp) %NThe order of the filter %Ap--ripple in the passband %Ast-a stopband Rs dB down from the peak value in the passband %Wp--the passband width [be,ae]=ellip(3,0.6,32,0.4); hellip=dfilt.df2(be,ae); f=0:0.001:0.4; g=grpdelay(hellip,f,2); g1=max(g)-g; [b,a,tau]=iirgrpdelay(6,f,[0 0.4],g1); hallpass=dfilt.df2(b,a); hoverall=cascade(hallpass,hellip); hFVT=fvtool([hellip,hoverall]); set(hFVT,Filter,[hellip,hoverall]); legend(hFVT,Lowpass Elliptic filter,Compensated filter); clear; [num1,den1]=ellip(3,0.6,32,0.4); [GdH,w]=grpdelay(num1,den1,512); plot(w/pi,GdH); grid xlabel(\omega/\pi); ylabel(Group delay, samples); F=0:0.001:0.4; g=grpdelay(num1,den1,F,2); % Equalize the passband Gd=max(g)-g; % Design the allpass delay equalizer [num2,den2]=iirgrpdelay(6,F,[0,0.4],Gd); [GdA,w] = grpdelay(num2,den2,512); hold on; plot(w/pi,GdH+GdA,r); legend(Original Filter,Compensated filter); 2．设计巴特沃兹模拟低通滤波器，其滤波器的阶数和3-dB 截止频率 由键盘输入，程序能根据输入的参数，绘制滤波器的增益响应。 clear; N=input(Type in the order N = ); Wn=input(Type in the 3-dB cutoff frequency Wn = ); %模拟频率 [num,den]=butter(N,Wn,s); w=0:2*Wn; h=freqs(num,den,w); plot(w,20*log(abs (h))),grid; 3．已知系统的系统函数为： 1−0.2z−1 +0.5z−2 H (z) = −1 −2 −3 −4 1+3.2z +1.5z −0.8z +1.4z 用MATLAB 进行部分分式展开，并写出展开后的表达式。 % Partial-Fraction Expansion of Rational z-Transform num = [0 0 1 -0.2 0.5]; den = [1 3.2 1.5 -0.8 1.4]; [r,p,k] = residuez(num,den); disp(Residues);disp(r) disp(Poles);disp(p) disp(Constants);disp(k) 4．设计切比雪夫I 型IIR 数字高通滤波器，其性能指标为：通带波 纹p 0.5dB ，最小阻带衰减 = 43dB ，通带和阻带边缘频率 s p 0.75 rad 和s 0.35 rad 绘制所设计的滤波器增益响应。 %a4 disp(prewapping is done,and T=2); Wp