数字信号处理.doc

数字信号处理实验报告学 部：信息科学与技术学部学 号：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,Xk8,.);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矩形窗（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);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: ，其中，求y（n）。用函数y=filter(p,d,x)计算输出y，p表示x前面的系数矩阵，d表示y前面的系数矩阵。要求绘制。要求绘制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=b