1、数 字 信 号 处 理 实 验 报 告 学 部：信息科学与技术学部 学 号：090102031126 姓 名：肖 伟 指导老师：刘怡 专业班级：09电子信息工程（3）班 2011年12月20日试验一 n=-4:6; x=1,0,2,-1,2,1,2,4,0.5,0,2; subplot(3,2,1);stem(n,x,.);line(-6,6,0,0) 第二题 n=0:7; xn=1 1 1 1 1 1 1 1; subplot(3,2,1);stem(n,xn,.);line(0,8,0,0) Xk8=fft(xn,8); fft(xn,8); subplot(3,2,1);stem(n,X
2、k8,.);line(0,8,0,0) abs(Xk8); angle(Xk5); subplot(3,2,1);stem(n,angle(Xk8),.);line(0,8,0,0) U(n)图 DFT 图相位图 第三题 实部 y1=real(Y); plot(y1);stem(y1,.);line(0,16,0,0) 虚部 y2=imag(Y); plot(y2);stem(y2,.);line(0,16,0,0)幅值 Y=ifft(K1,12); y=abs(Y); plot(y);stem(y,.);line(0,16,0,0)实验二 一、在 MATLAB 产生窗函数的方法： 1矩形窗（
3、Rectangle Window） 代码： w=boxcar(8); b=fir1(7,0.4,w); freqz(b,1) 2 三角窗（Triangular Window） w= triang (8); b=fir1(7,0.4,w); freqz(b,1)3 汉宁窗（Hanning Window） w= hanning (8); b=fir1(7,0.4,w); freqz(b,1) 4 海明窗 Hamming Window） w= hamming (8); b=fir1(7,0.4,w); freqz(b,1)5布莱克曼窗（Blackman Window） w= blackman (8)
4、; b=fir1(7,0.4,w); freqz(b,1) 6 凯撒窗（Kaiser Window） w= kaiser (8,6); b=fir1(7,0.4,w); freqz(b,1)二、 Q1: ， ( ) 0.7 ( 1) 0.45 ( 2) 0.6 ( 3) 0.8 ( ) 0.44 ( 1) 0.36 ( 2) 0.02 ( 3) y n y n y n y n x n x n x n x n 其中 ，求 y（n） 。 ( ) ( ),0 1 x n n n N 用函数 y=filter(p,d,x)计算输出 y，p表示 x前面的系数矩阵，d表示 y前面的系 数矩阵。要求绘制。要
5、求绘制 y的图。 a=1,0.7,-0.45,-0.6; b=0.8,-0.44,0.36,0.02 ; x=1,zeros(1,20); yfilter=filter(b,a,x); stem(n,yfilter); title(yfilter); xlabel(x); ylabel(y);Q2：设计指标为：wp=0.2 Rp=0.25dB wa=0.3 As=50dB的低通数字FIR 滤波器 wp=0.2*pi; wa=0.3*pi; w=(wp+wa)/2; Window=Blackman(47) ; b=fir1(46,w,window) ; freqz(b,1) ;Q3：采样频率为1Hz，通带临界频率fp=0.2Hz,通带内衰减小于1dB（ap=1） ； 阻带临界频率fs=0.3Hz，阻带内衰减大于25dB(as=25)。设计一个数字滤波器 满足以上参数。 （绘制幅频响应和相频响应曲线函数为freqz（b，a，N） ） FS=1; n,Wn=buttord(0.2/(FS/2) ,0.3/(FS/2) ,1,25) ; b,a=butter(n,Wn) ;freqz(b,a,512,FS) ;
