数字信号处理实验报告(自己的实验报告)

1、数字信号处理实验报告西南交通大学信息科学与技术学院姓名：伍先春学号：20092487班级：自动化1班指导老师：张翠芳实验一序列的傅立叶变换实验目的进一步加深理解dfs,dft算法的原理;研究补零问题;快速傅立叶 变换(fft)的应用。实验步骤1. 复习dfs和dft的定义，性质和应用；2. 熟悉matlab语言的命令窗口、编程窗口和图形窗口的使用；利用提 供的程序例子编写实验用程序；按实验内容上机实验，并进行实验结 果分析；写出完整的实验报告，并将程序附在后面。实验内容1.周期方波序列的频谱试画出下面四种情况下的的幅度频谱,并分析 补零后，对信号频谱的影响。xn) = cos(0.48n) +

2、 cos(0.5271)1, mn < n < mn + £ - 1x(77)= 0, mn + l < n< (m + 1)n 1 m - o,± 1,±2, 文伙)=dfs x(n)(1)厶=5,n = 20;(2)l = 5,n 二 40;(3)厶=5,n = 60; (4) l = 7,7v = 602.有限长序列x(n)的dft _(1) 取x(n)(n=0:10)时，画出x(n)的频谱x(k)的幅度；(2) 将(1)中的x(n)以补零的方式，使x(n)加长到(n:0100)时, 画出x(n)的频谱x(k)的幅度；(3) 取x(n

3、)(n:0-100)时，画出x(n)的频谱x(k)的幅度。利用fft 进行谱分析已知：模拟信号x(t) = 2 sin(4加)+ 5 cos(8 加)以t=0.01n(n=0:n-l)进行采样，求n点dft的幅值谱。请分别画出n=45; n=50;n=55;n=60时的幅值曲线。数字信号处理实验一1. (1)l=5;n=20;n=l:n;xn=ones(1 l)zeros (1a n-l);xk=dfs(xnz n);magxk=abs(xk(n/2 + 1:n) xk (1:n/2 + 1); k=-n/2:n/2;figure (1)subplot(2a1a1);stem(n,xn);xl

4、abel(1n1);ylabel(1xtide(n) 1);title ( 'dfs of sq.wave:l=5z n=20');subplot (2,1,2);stem(k,magxk);axis ( -n/2小/2。16);xlabel ();ylabel ( 1 xtide (k) 1);l=5;n=40;n=l:n;xn=ones(1,l),zeros(1,n-l);xk=dfs(xn,n);magxk=abs(xk(n/2+l：n) xk(l:n/2+l)；k=-n/2:n/2;figure (2)subplot(2,lz1);stem (n, xn) ; xlab

5、el ( 1 n 1 ) ; ylabel ( 1 xtide (n) 1 );title ( 'dfs of sq.wave:l=5a n=40');subplot(2flf2);stem(k magxk);axis(-n/2,n/2,0,16);xlabel(1k1);ylabel(1xtide(k) 1);l=5;n=60;n=l:n;xn=ones(lfl)fzeros (1,n-l);xk=dfs(xnfn);magxk=abs ( xk(n/2 + 1:n) xk(1:n/2 + 1)； k=-n/2:n/2;figure(3)subplot(2flfl);stem

6、 (nz xn) ; xlabel ( 1 n 1) ; ylabel ( 1 xtide (n) 1); title ( 1dfs of sq.wave:l=5a n=601);subplot (2,lz 2);stem(k,magxk);axis(-n/2,n/2,0,16);xlabel ( 1 k 1) ; ylabel ( 1 xtide (k) 1 );l=7;n=60;n=l:n;xn=ones(lz l)9 zeros (1,n-l);xk=dfs(xn,n);magxk=abs(xk(n/2+1:n) xk(1:n/2+1)； k=-n/2:n/2;figure(4)sub

7、plot (2 w);stem(na xn);xlabel(1n1);ylabel(1xtide(n) 1);title ( 1dfs of sqwave:l=7 fn=601);subplot(2,1,2);stem(kz magxk);axis(-n/2,n/2,0,16);q figure 4l cd 回file editview insert tools desktop window help d d k吳的相0 口xlabel ( 1 k 1) ; ylabel ( 1 xtide (k) 1 );dfs of sq wave丄=7.n=601? (专 p=x2. (1)m=10;n

8、=10;n=l:m;xn=cos(048*pi*n)+cos(052*pi*n);nl=0:1:n-l;yl=xn(1:1:m)zeros(1n-m); figure(1)subplot(2,1,1);stem(nl,yl);xlabel(1n1 );title ( 1 signal x (n),0<=n<=101);axis(0,n,-25,2 5);yl=fft(yl) ;magyl=abs(yl (1:1:n/2 + 1);kl=0:l:n/2;wl=2*pi/n*kl;subplot (2/1,2);title(1 samples of dtft magnitude 1);

9、 stem(wl/piz magyl); axis ( 0zlz 0z10);xlabe丄('frequency in pi units 1);m=10;n=100;n=l:m;xn=cos(0.48*pi*n)+cos(0.52*pi*n);nl=0:1:n-l;yl=xn(1:1:m), zeros(1,n-m); figure(2)subplot(2alz1);stem(nl yl);xlabel(1n1); title ( 1 signal x (n), 0<=n<=101);axis(0,nz-2.5,2.5);yl = f f t (yl) ; rnagyl=a

10、bs (yl (1:1:n/2 + 1);kl=0:1:n/2;wl=2*pi/n*kl;subplot (2,1/2);title(1 samples of dtft magnitude 1); stem(wl/pi,magyl); axis (0,1,0,10); xlabel(1 frequency in pi units 1);(3)m=100;n=100;n=l:m;xn=cos(048*pi*n)+cos(052*pi*n);nl=o:l:n-l;yl=xn(l:l:m),zeros(1,n-m); figure(3)subplot(2z1,1);stem(nl,yl);xlabe

11、l(1n1); title ( 1 signal x (n)0<=n<=1001);axis(0,n,-2.5,2.5);yl=fft(yl);magyl=abs(yl(1:1:n/2+1);kl=0:l:n/2;wl=2*pi/n*kl;subplot (2,1,2);title(1 samples of dtft magnitude 1); stem(wl/pi,magyl); axis (0,1,0,10);xlabel(1 frequency in pi units 1);3.figurc(l)subplot(292j)n=45;n=0:n-l;t=0.01*n;q=n*2

12、*pi/n;x=2*sin(4*pi*t)+5*cos(8*pi*t);y=fft(x,n);plot(q,abs(y)stem(q,abs(y)titlecfft n=45*)%subplot(2,2,2)n=50;n=0:n-l;t=0.01*n;q=n*2*pi/n;x=2*sin(4*pi*t)+5*cos(8*pi*t);y=fft(x,n);plot(q,abs(y)titlecfft n=5(t) subplot(2923)n=55;n=0:n-l;t=0>01*n;q=n*2*pi/n;x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,n);p

13、lot(q,abs(y)title(ffft n=55r)%subplot(2,2,4)n=16;n=0:n-l;t=0.01*n;q=n*2*pi/n;x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,n);plot(q,abs(y)title(：fft n=16*)functionxk=dfs(xnn) n=o:l:n-l;k=o:l:n-l;wn=exp(-j*2*pi/n);nk=n 怙 k;wnnk=wn.ank;xk=xn*vvnnk;fj figure 1实验二用双线性变换法设计iir数字滤波器一、实验目的1. 熟悉用双线性变换法设计iir数字滤波器的

14、原理与方法；2. 掌握数字滤波器的计算机仿真方法；3. 通过观察对实际心电图的滤波作用，获得数字滤波器的感性知 识。二、实验内容1. 用双线性变换法设计一个巴特沃斯低通iir滤波器，设计指标 参数为：在通带内频率低于0.2n时，最大衰减小于ldb；在阻 带内0.3ji,兀频率区间上，最小衰减大于1 5db。2. 以0.2n为采样间隔，打印出数字滤波器在频率区间0, 0.2兀 上的幅值响应曲线。3. 用所设计的滤波器对实际的心电图信号采样序列x (n)(见教科 书上254页实验二的第六部分)进行仿真滤波处理，并分别打 印出滤波前后的心电图信号波形图，观察总结滤波作用与效果。三、实验步骤1. (参

15、考数字信号处理教材实验二)；2. (参考数字信号处理教材实验二)；3. 用matlab语言编写仿真滤波程序，完成实验内容13；4. 写出完整的实验报告并回答教科书上的思考题。四、iir滤波器设计的常用函数1. butter功能：butterwoeth (巴特沃斯)数字滤波器设计。格式：b, a=butter (n, wn)可设计出截止频率为wn的n阶低通butterwoeth滤波器，其滤波器为hq 二 b(z)二 b(l) + b(2)z+ + bg + l)z”a(z) 14- a(2)zl + + a(n + l)z_,j截止频率是滤波器幅度下降至i/j2处的频率.2. buttordhr

16、 (巴特沃斯)滤波器阶的选择格式n,wn=buttord(wp,ws,rp,rs)数字域n,wn=buttord(wp,ws,rp,rs/s,)模拟域说明：buttord可在给定滤波器性能的情况下，选择模拟或数字butterword 滤波器的最小阶，其中wp和ws分别是通带和阻带的截止频率，其值 0<()<1，当其值为1时，表示05 fs ,rp,rs分别是通带和阻带区的波纹系数。例如：设计一低通滤波器，通带范围0100血,通带波纹小于3db,阻带为30db, 并利用最小的阶来实现，其程序如下wp=l00/500；ws=200/500；n, wn=buttord (wp,ws,3,

17、 30)；b, a=butter (n, wn)；freqz (b, a, 512, 1000)mp) pmcbew(s ajb p) sejzd21、t=pi/2;fs=l/t;wpz=0.2;wsz=0.3;wp=2*tan(wpz*pi/2);ws=2*tan(wsz*pi/2);rp=l;rs=15;n,wc=buttord(wp,ws,rp,rs,s,)b,a=butter(n,wc,'s')k=0:511;fk=0:0.2*pi:pi/2;wk=2*pi*fk;hk=freqs(b,a,wk);subplot(l；l 91);plot(fk,20*logl0(abs

18、(hk);grid on xlabelc频率(hz) lylabel('幅度(db) f)n 二wc -0. 76621.0000(ap) §140000.20242.96054.3r224. 11242. 57281. 02050. 2024020-4060-80-100-1200.20.40.6 0.8 1频率(hz)2-2b, a =butterw(0. 2, 0. 3, 1, 15); figure(l)freqz (b, a, 0:0. 02*pi:0 2*pi)b, a=butterw(0. 2, 0. 3, 1, 15);xn二-4, -2, 0, -4, -6

19、, -4, -2, -4, -6, -6, -4, -4, -6, -6, -2, 6, 12, & 0, -16, -38 ,-60, -84, -90, -66, -32, -4, -2, -4, 8, 12, 12, 10, 6, 6, 6, 4, 0, 0, 0, 0, 0, -2, -4 ,0, 0, 0, -2, -2, 0, 0, 一2, 一2, 一2, 一2, 0;n=length(xn);figure (2)subplot (2, 1, 1) ;plot (xn)title (' before，)yn=f 订tcr (b, a, xn);subplot (2,

20、 1,2); plot (yn)title (' aftei?)7276780.020 040 060 080 10 12 0 140 16 0 180 2normalized frequency (xg< rad/sample)0.020.040.060.080.10.12 0 140.16 0.180.2normalized frequency (xx rad/sample) _ o 0500050-1d (s8aj68p) asmq准备函数：butterwfunction b,a=butterw(wp,ws,ap,as)t=0.2;wp=(2*tan(wp/2)/t;ws=

21、(2*tan(ws/2)/t;11=0. pap;t2=0.1 *as;k=(10atl-l)/(10at2-l)a0.5;r=ws/wp;n=-loglo(k)/loglo(r);n=floor(n)+l;wn=wp*( 10a0.1 *ap-1 )a(-l/(2*n);z,p,k=butter(n,wn；s,);lp=length(p);z=-ones(l,lp);z=z'p=l/(l-p)；k=k/(2/t)alp);b,aj=zp2tf(z,p,k);实验三用窗函数法设计fir数字滤波器1. 实验目的及原理方法参考(301-302)页2. 实验内容 用matlab产生各种窗函数

22、利用窗函数设计fir滤波器参考程序:(1)用matlab实现各种窗函数m=200;a=boxcar (m);矩形窗b=hanning (m);汉宁窗 c=hamming (m);海明窗 d=blackman (m);布莱克曼窗 e=kaiser (m, 7.865);凯泽窗 m=l:200;plot(m,af *r*1zm,b,'g+',m,c, fy*!d, fb.1)(2)利用窗函数设计fir滤波器设计具有下列指标(op=0.257t,rp=0.25db, =03兀,rp=50db的低通数 字滤波器。要求：选择合适的窗函数；画岀滤波器的频率特性将有理分式变换成部分分式表示

23、将部分分式换成有理分式变表示几个主要的matlab函数(1) r,p,k=residuez(b,a)(2) lb,aj=residuez(r,p,k)式中各参数意义如下：h(z)=仇+b|z+ + $/"1畀+饥1) +£(加一卅+ 1)广(旳)z + + q“z厂,厂)d 1一0(1比一1 - /?(/2)z(3) freqz(b,a)画连续系统的频率特性曲线(4) impz(b,a, 101)画系统的脉冲响应 三、实验程序1m=200;a=boxcar(m);% 矩形窗 b=hanning(m);% 汉宁窗 c=haniming(m); % 海明窗 d=blackman(m); %布莱克曼窗 e=kaiser(m,7.865); % 凯泽窗 m=l:200;2)wp=0.25*pi;ws=0.3*pi;tr_wid th=w s-wp; % 求过渡带宽n=ceil(6.6*pi/tr_width)+l % 求窗的宽度n=0:l:n-l;wc=(ws+wp)/2;%理想低通滤波器的截止频率 hd=ideal