数字信号处理基础实验报告

1、中南大学数字信号处理实验报告课程名称数字信号处理指导教师去逐学 院 信息科学与工程学院专业班级学号 姓名实验一 常见离散时间信号的产生和频谱分析一、实验目的(1) 熟悉matlab应用环境，常用窗口的功能和使用方法；(2) 加深对常用离散时间信号的理解；(3) 掌握简单的绘图命令；(4) 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离 散信号进行频域分析。二、实验内容及要求(1)复习常用离散时间信号的有关内容；常用离散时间信号a) 单位抽样序列5(")=e如果55)在时间轴上延迟了 k个单位，得到- k)即:n-k) = l h = k0心0b) 单位阶跃序列0 n &

2、lt; 0、加f 宀fl 0 < n < a-lc) 矩形序列心%其他d) 正弦序列x(n) = a sin( wn + (p)e) 实指数序列x(n) = anu(n)f) 复指数序列x(n) = d"”(2)用matlab编程产生上述任意3种序列(长度可输入确定，对(d) (e) (f)中的参数可自行选择)，并绘出其图形；程序如下：1) 单位阶跃序列：n 二20:20;xn=heaviside(n);xn(n=0)=l;plot(n,xn);stem(n,xn);axis(-20 20 0 1.2);title('单位阶跃序列,);xlabel('n,

3、);ylabel('u(n),);box on得到图像如2单位阶跃序列2) 单位抽样序列：n 二20:20;xn=heaviside(n)-heaviside(n+1);xn(n=0)=l;plot(n,xn);stem(n,xn);axis(-20 20 0 1.2);titlec 单位抽样序列');xlabelcn');ylabel('delta(n)');box on得到图像如下:0.8单位抽样序列0.60.40.23) 矩阵序列:n=-20:20;n=5;xn=heaviside(n)-heaviside(n-n);xn(n=0)=l ;xn(n

4、=n)=0;plot(n,xn);stem(n,xn);axis(-20 20 0 1.2);titled矩阵序列');xlabel(h);ylabel('r_n(n)');box on得到图像如下：矩阵序列qqqqq0.80.60.40.210154) 正弦序列： n=-40:40;a=2;w=pi/8;f=pi/4;xn 二 a*sin(w.*n+f);plot(n,xn);stem(n,xn);axis(-40 40 -4.2 4.2)title('正弦序列,);xlabel('n,);ylabel('x(n),);box on得到图像如下

5、:-40-30-2010010203040n(3)混叠现象对连续信号x(/) = sin(2*pzvoi*o其中，/oi = 500hz进行采样，分别取采样频率2000hz, 1200hz,800hz,观察|x(/)|的变化，并做记录(打印曲线), 观察随着采样频率降低频谱混叠是否明显存在，说明原因。程序如下：_(l) fs=2000hz时的釆样 fs=2000;t=o:l/fs:o；x二sin(2*pi*500*t)subplot(2,l,l);plot(t,x)；title('x=sin(2*pi*500*t)');xlabel(f);ylabel('幅度);n=l

6、ength(t);subplot(2,l,2);plot(-n/2:n/2-1 )*fs/n,abs(fftshift(fft(x,n);axis(o 1000 0 100)title(,fft(x=sin(2*pi*500*t)');xlabel(,f/hz,);ylabel('|x(eajomega)r);(2) fs二1200hz时的采样fs=1200;t=0:l/fs:0.1;x=sin(2*pi*500*t)subplot(2,l,l);plot(t,x);title(,x=sin(2*pi*500*t),);xlabel(,t,);ylabel(幅度');n

7、=length(t);subplot(2,l,2);plot(-n/2:n/2-l)*fs/n,abs(fftshift(fft(x,n);axis(ro 1000 0 100)title(,fft(x=sin(2*pi*500*t)');xlabel(,f/hz,);ylabel('|x(eajomega)r);得到结果图如下:x=sin(2*pi*500*t)fft(x=sin(2*pi*500*t)100500i 01002003004005006007008009001000f/hz(3) fs=800iiz时的釆样fs=800;t=0:l/fs:0.1;x=sin(2

8、 *pi *500* t)subplot(2,l,l);plot(t,x);title(,x=sin(2*pi*500*t),);xlabel(,t');ylabel('幅度');n=length(t);subplot(2,l,2);plot(-n/2:n/2-1 )*fs/n,abs(fftshift(fft(x,n);axis(o 1000 0 100)title(,fft(x=sin(2*pi*500*t)');xlabel(,f/hz,);ylabel('|x(eajomega)r);得到结果图如下:x=sin(2*pi*500*t)结论：当采样

9、频率fs小于两倍信号(这里指是信号)最大频率500hz时，经过 采样就会发生频谱混叠，这使得采样后的信号序列频谱不能真实地反映原信号的 频谱。fs=2000ilz时和fs=1200hz时未出现混叠，当fs=800ilz<2*500hz时，出 现明显混叠现象。(4)截断效应给定x() = cos(/?),截取一定长度的信号y(n) = x(/i)iv(/t),)为窗函数,4长度为n, w(n) = rn(n).做2n点dft变换，分析当n逐渐增大时，分析是否有频谱泄露现象、主瓣的宽度变化？如何减小泄露？ 程序如下:%n=20 时：n=20;x=cos(pi/4*n);w=heaviside

10、(n)-heaviside( n-n);w(n=0)=1 ;w(n=n)=0;y=x.*w;subplot(2,l,l);stem(n,y);box on;title(n=20 截断j;xlabel(h);ylabel(幅度)； subplot(2,l,2);n1=2*n;plot(-n:n-l)*2*pi/nl/pi,abs(fft(y,nl);titlec 频谱,);xlabel(,omega/pi,);ylabel(,|x(eajomegaj)r);n=20截断-0 o0.5得到结果图如下：-0.5x-1i 03iiiiiiii0102030405060708090100n%n=30dt

11、：n=30;x=cos(pi/4*n);w=heaviside(n)-heaviside(n-n);w(n=0)=1 ;w(n=n)=0;y=x.*w;subplot(2,l,l);stem(n,y);box on;title(n=30 截断);xlabel(,n,);ylabel(幅度')；subplot(2,l,2);n1=2*n;plot(-n:n-1 )*2邙i/n l/pi,abs(fft(y,n 1);title('频谱');xlabel('omega/pi,);ylabel('|x(eajomega)r);n=30截断%n=50: n=0:

12、1:100;n=50;x=cos(pi/4*n);w=heaviside (n)heaviside(nn);w(n=0)=l ;w(n=n)=0;y=x.*w;subplot(2,l,l);stem(n,y);box on;title(fn=5 0 截断');xlabel(h);ylabel('幅度);subplot(2,l,2);n1=2n;plot(-n:n-l)*2*pi/nl/pi,abs(fft(y,n 1);title('频谱');xlabel(,omega/pi');ylabel('|x(eajomega)r); 得到结果图如下：n

13、=50截断%n=100 时：n=0:l:100;n=100;x=cos(pi/4*n);w=heaviside (n)heaviside(nn);w(n=0)=l ;w(n=n)=0;y=x.*w;subplot(2,l,l);stem(n,y);box on;title(fn=100 截断');xlabel(h);ylabel('幅度);subplot(2,l,2);n1=2n;plot(-n:n-l)*2*pi/nl/pi,abs(fft(y,n 1);title('频谱');xlabel('omega/pi');ylabel(lx(ea j

14、omega )|');结论：dft变换频谱泄漏的根本原因是信号的截断。即时域加窗，对应 为频域卷积，因此，窗函数的主瓣宽度等就会影响到频谱。当n取20、30、50 截断时，频谱都有一定程度的泄漏，当n逐渐取值增大至100,泄漏减小，主瓣 宽度逐渐变窄。(5)栅栏效应给定x(z?) = r4(n),分别计算x(ejw)在频率区间0,2龙上的16点、32点、64点等间隔采样，绘制x(/)采样的幅频特性图，分析栅栏效应，如何减小栅 栏效应？程序如下： xn=l 1 1 11; yl=fft(xn,16); y2=fft(xn,32); y3=fft(xn,64);誓n=100截断60频谱40

15、3x 20-1-0.8-0.6-0.4-0.200.20.40.60.81ovpisubplot(3 丄 1);stem(0:15)*2*pi/16/pi,abs(y l);box on;titlec16 点 dft 幅频图');xlabel('omega/pi');ylabel('幅度); subplot(3,l,2);stem(0:31 )*2*pi/32/pi,abs(y2);box on;title(*32 点 dft 幅频图');xlabelcomega/pi');ylabel('幅度); subplot(3,l,3)；stem

16、(0:63)*2*pi/64/pi,abs(y3);box on;titlec64 点 dft 幅频图');xlabel('omega/pi');ylabelc幅度);得到结果图如下：16点dft幅频图civpi32点dft幅频图lgovpi64点dft幅频图ovpi栅栏效应，就一定意义上看，dft來观看频谱只能在离散点上看到真实频 谱，这样就可能发生一些频谱的峰点或谷点被“尖桩的栅栏”所挡住，不能被观 察到。减小栅栏效应的一个方法就是借助在原序列的末端添补一些零值，从而变 动dft的点数。这一方法实际上是人为地改变了对真实谱采样的点数和位置，相 当于搬动了每一根“尖桩

17、栅栏”的位置，从而使得频谱的峰点或者谷点暴露出来。 我们在这里所做的就是为了减小栅栏效应，使频域抽样更密，增加频域抽样点数 n,使一个周期内的点数增加，但不改变原有的记录数据。(6)提高题：给定信号 xjn) = 5sin(/?), x2(n) = cos( n), xj(n) = x1(n) + x2(w) ® 分别计算x(e>)在频率区间0,2刃上的频谱，观察其幅频特性图，分析是否存在 谱间干扰，如何减小谱间干扰？ 程序如下:l.n=40 时n 二 20:19;n=length(n)xlu5*sin(pv5n)_ x2hcos(pi/4 诧 n)； x3ux1+x2八sub

18、pot(221)pot(on 1 )诜2 诜 pvn、plabs(ffshif(b(x lz)；hta-xl 遵谪)eab2-(omeg 巴 pi)；ylab2-(_x(？9omegas; axiao 2 0 400dsubplot(222)potxon i)*2*pi/n、pipbs(ffshift(fft(x2z)；tit?r(x2 遊e)；xlabe5omega/pi)；ypb2-(_x(e>9omega)_)； axis(6 2 000)subpoi(223)plot(pn 1 )*2*pi/n、pipbs(fftshiftat(x3z)；s-12x3<)babe5omeg

19、 巴 pi)abel(-x(e>50mega)-)； axisao 2 0 400)400300x200soxis.0500ovpi0 05 1 5gyp-2n"100 巨.is50少50ix00 05 1 15ovpif1h5p49 八nilength(n)xl=5*sin(pi/5j)；x2"cos(pv4 诧 n)八x3hx1+xj?subp-£(22-)pot(3n 1)诜 2 诜 pi/n/pibbs(fftshift(fft(x 12)；5-a-x 1 遵蒔)eabe5omeg a、pi_)八 ylab2-(-x(e> iomega 一 )

20、-)； axis(o2 0 400)subplot(222)p-o(pn 一 )*2兴pi/n、pibbs(ffshif(ff(x2z)r ti二 2x2li«)eabe 三omega/piy_ab2.(_x(e>jsmeg 巴二)；axis(©2 0 loo) subplot(223)ploson 一)诧 2 诧 pi/n、pipbs(fftshift(fft(x3z)rs-a.x3w)；xlabel(omega/pi-)；ysrbel(-x(e>jomega7)-)； axis(o2 0 402)q.o ou1o|x（曲）1ooo14i|x（如）1goood

21、i3nh500三n 卩 25p24tnhongrh(n)xlh5*sin(pv5-x-n)八x2hcos(pv4 诜 ro-xjhxl+xn subplor(221) plo(n*2*2*2*pi/n、pibbs(fkshift(s(xln)； tita.xl<)eabe三omega、pi-)-yfbe5x(e>9omega)-)；axis(l2 1.2 0 2000) subpor(222) pots诜2诜2诜2诜pvn、plabs(sshift(fft(x2z)八 tito(x2<)eabe三omega、pi-)-yfbe5x(e>9omega)-)；axis(l2

22、 1.2 0 400)cjg00000 l00200000寸a(oxicopxi 抄 co迪w.sx(oooeoz一 z一 丄)sxe o(e 斜 luo三 vgx_)oqewad、e 和 eocaqe 一汽(.班巔 ex. )3-2 _(ndh)uj)sqs£)sq5d、n=d£ied (ezz)o-dqnsx(n7x)yj)5qs)sq5d、§r(i n0)oid(ezz)lokqns (ooo0 z odsxeo(e 斜 luo三 vgx_)oqewad、e 和 eocaqe 一汽(.班巔 ix. )3-2"(nj =)£)主 qse)sq5d、n/-d 沃尖(i noved (i77)o-dqns “x+ixh咚 xu* 寸二 d)sofzx (ul二d)ws:上 x(u)q导hn "6660021 fc 000址.寸title(!x2 频谱,);xlabel(,omega/pi,);ylabel(,|x(eajomega)r); axis(o 2 0 1000)subplot(2,2,3)plot(0:n-l )* 2*pi/n/pi,abs(fftshift(fft(x3,n);titlecx3 频谱,);xlabel(,omega/pi,);ylabelc|x(eajome