数字信号处理课程设计报告

1、现代信号处理课程设计任务实验报告 班级：通信080*姓名：*学号：*指导老师：*目 录第一部分：实验设计目的第二部分：实验设计原理第三部分：实验设计方案第四部分：实验设计过程第五部分：个人小结参考文献第一部分：实验设计目的1全面复习课程所学理论知识，巩固所学知识重点和难点，将理论与实践很好地结合起来。2提高综合运用所学知识独立分析和解决问题的能力；3熟练使用一种高级语言进行编程实现。第二部分：实验设计原理（一）信号分析1、编制信号生成程序，产生下述各序列，绘出它们的时域波形1） 三角波序列2）反三角波序列3） gaussian序列4） 衰减正弦序列 对连续信号 进行采样，可得到测试序列令a=5

2、0,采样周期t=1ms,即fs=1000hz,f0=62.5,a=1005）已调信号 ，其中调制信号，载波为。6）无限长序列2. 对上述信号完成下列信号分析1）对三角波序列和反三角波序列，作n=8点的fft，观察比较它们的幅频特性，说明它们有什么异同？绘出两序列及其它们的幅频特性曲线。 在和的尾部补零，作n=16点的fft，观察它们的幅频特性发生了什么变化？分析说明原因。2)、观察高斯序列，固定信号中的参数p=8，令q分别等于2，4，8，观察它们的时域和幅频特性，了解当q取不同值时，对信号序列的时域幅频特性的影响；固定q=8，令p分别等于8，13，14，观察参数p变化对信号序列的时域及幅频特性

3、的影响，观察p等于多少时，会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。3）观察衰减正弦序列的时域和幅频特性，绘出幅频特性曲线，改变采样频率fs,使fs=300hz， 观察此时的频谱的形状和谱峰出现位置？说明产生现象的原因。4）观察已调信号的时域与频域特性，绘出其幅频特性曲线，并对它进行解调，以恢复原调制信号。5）对序列 ，截取序列长度m=10点，计算其频谱，画出幅频特性曲线；截取序列长度m=20点，计算其频谱，画出幅频特性曲线；（二）、数字滤波器的设计与实现设计和实现iir或fir数字滤波器：要求：输入数字滤波器的滤波指标,包括通带截止频

4、率,通带最大衰减,阻带截止频率,阻带最小衰减。设计出相应的数字滤波器。显示得到的滤波器的系数和阶数，求出系统的单位脉冲响应h(n)。并画出设计得到的滤波器的幅频特性图(要有坐标标度)和相频特性。1对于iir数字滤波器：使用双线性变换法。可选用bw型、cb型和椭圆型等滤波器类型来进行设计2对于fir数字滤波器：要对多种窗口（三种及以上）、不同阶数设计结果进行分析比较。（三）、数字滤波1）根据设计的滤波器，选择2到3个测试信号进行滤波。展示滤波后信号的幅频与相频特性，分析是否满足滤波要求。对同一滤波要求，根据输出信号频谱，对比分析各类滤波器的差异。2）将一段语音作为测试信号，通过频谱展示和语音播放

5、，对比分析滤波前后语音信号的变化，进一步加深对数字信号处理的理解。第三部分：实验设计方案模拟带通滤波器在信号的检测和传输过程中起着很重要的作用，但在传统模拟带通滤波器的设计过程中需要大量繁琐的数值计算，如果手工计算则费时费力。也有许多公司开发了相关的滤波器设计软件，但是它们大多数使用起来不容易上手，而且价格比较昂贵。matlab中的图形界面编程环境guide为用户开发软件界面提供了强有力的工具。在图形界面编程环境guide下，设计好滤波器的设计工具界面后，对界面上各个控件的回调函数进行编写就能设计出交互式滤波器的分析设计软件。本文在matlab gui基础上设计了针对chebyshev型模拟带

6、通滤波器的设计分析软件。用户在使用时只需输入要设计的滤波器设计指标值，选择要设计的滤波器类型，即可在分析滤波器特性时只需选择相应菜单就能得出对应的特性。该工具解决了用户在设计chebyshe型模拟带通滤波器时大量繁琐的数值计算，不仅可以作为研究设计chebyshev型模拟带通滤波器的辅助工具，还可以作为设计滤波器等相关课程的辅助教学软件。1 模拟带通滤波器的设计理论目前，模拟高通滤波器、带通滤波器和带阻滤波器的设计方法都是先将要设计的滤波器技术指标通过某种频率转换关系转换成模拟低通滤波器的技术指标，并依据这些技术指标设计出低通滤波器的转移函数，然后再依据频率转换关系变成所要设计的滤波器的转移函

7、数。工程实际中，设计高通、带通、带阻滤波器的常用方法是借助低通原型滤波器，经频率变换和元件变换得到。在matlab中，图形用户界面(graphical user interface，gui)是matlab中一个专用于gui程序设计的向导设计器，而gui是由各种图形对象，如图形窗口、图轴、菜单、按钮、文本框等构建的用户界面，是人机交互的有效工具和方法。通过guide可以很方便地设计出各种符合要求的图形用户界面。用户通过一定的方法(如鼠标或键盘)选择、激活这些图形对象，使计算机产生某种动作或变化，比如实现计算、绘图等。gui设计既能以基本的matlab程序设计为主，也能以鼠标为主，利用guide工

8、具进行设计，也可综合以上两种方法进行设计。guide主要是一个界面设计工具集。matlab将所有gul支持的用户控件都集成在这个环境中并提供界面外观、属性和行为相应方式的设置方法。guide将用户保存设计好的图形用户界面保存在一个fig资源文件中，同时自动生成包含图形用户界面初始化和组件界面布局控制代码的m文件，这个m文件为实现回调函数的编写提供了一个参考框架。fig文件是一个二进制文件，包含系列化的图形窗口对象。所有对象的属性都是用户创建图形窗口时保存的属性。该文件最主要的功能是对象句柄的保存。m文件包含gui设计、控制函数及控件的回调函数，主要用来控制gui展开时的各种特征。该文件基本上可

9、以分为gui初始化和回调函数2个部分，控件的回调函数可根据用户与gui的具体交互行为分别调用。第四部分：实验设计过程（一）信号分析1、利用matlab编制信号生成程序，产生下述各序列，绘出它们的时域波形，并对各信号进行相关分析：（1）三角波序列、反三角波序列及其fft三角波序列matlab源程序：subplot(1,1,1); n1=0:3;x1=n1+1;n2=4:7;x2=8-n2; stem(n1,n2,x1,x2,.);axis(-1,9,0,7);反三角波序列matlab源程序：subplot(1,1,1); n1=0:3;x1=4-n1;n2=4:7;x2=n2-3; stem(n

10、1,n2,x1,x2,.);axis(-1,9,0,10);它们的fft变换的matlab源程序：n1=8,n2=16;n1=0:n1-1;x1=1,2,3,4,4,3,2,1;subplot(4,2,1);stem(n1,x1,.);axis(0,10,0,5);title(x3(n);y1=fft(x1,n1);subplot(4,2,2);stem(n1,abs(y1),.);axis(0,10,0,8);title(fft,n=8);x2=4,3,2,1,1,2,3,4;subplot(4,2,3);stem(n1,x2,.);axis(0,10,0,5);title(x4(n);y2

11、=fft(x2,n1);subplot(4,2,4);stem(n1,abs(y2),.);axis(0,10,0,8);title(fft,n=8);n2=0:n2-1;x3=x1,zeros(1,8);subplot(4,2,5);stem(n2,x3,.);axis(0,20,0,5);title(x3(n);y3=fft(x3,n2);subplot(4,2,6);stem(n2,abs(y3),.);axis(0,20,0,20);title(fft,n=16);x4=x2,zeros(1,8);subplot(4,2,7);stem(n2,x4,.);axis(0,20,0,5);

12、title(x4(n);y4=fft(x4,n2);subplot(4,2,8);stem(n2,abs(y4),.);axis(0,20,0,20);title(fft,n=16)生成波形为：（2） gaussian序列及其fftp为定值：matlab源程序：n=0:15;x1=exp(-(n-8).2/2);subplot(3,2,1);stem(n,x1,.);axis(0,20,0,1);title(x5(n),p=8,q=2); y1=fft(x1,16);subplot(3,2,2);stem(n,abs(y1),.);axis(0,20,0,3);title(fft,p=8,q=

13、2);x2=exp(-(n-8).2/4);subplot(3,2,3);stem(n,x2,.);axis(0,20,0,1);title(x5(n),p=8,q=4);y2=fft(x2,16);subplot(3,2,4);stem(n,abs(y2),.);axis(0,20,0,4);title(fft,p=8,q=4); x3=exp(-(n-8).2/8);subplot(3,2,5);stem(n,x3,.);axis(0,20,0,1);title(x5(n),p=8,q=8);y3=fft(x3,16);subplot(3,2,6);stem(n,abs(y3),.);ax

14、is(0,20,0,6);title(x5(n),p=8,q=8);生成波形为：q为定值：n=0:15;x1=exp(-(n-8).2/8);subplot(4,2,1);stem(n,x1,.);axis(0,20,0,1);title(x5(n),p=8,q=8);y1=fft(x1,16);subplot(4,2,2);stem(n,abs(y1),.);axis(0,20,0,6);title(fft,p=8,q=8);x2=exp(-(n-13).2/8);subplot(4,2,3);stem(n,x2,.);axis(0,20,0,1);title(x5(n),p=13,q=8)

15、;y2=fft(x2,16);subplot(4,2,4);stem(n,abs(y2),.);axis(0,20,0,6);title(fft,p=13,q=8);x3=exp(-(n-14).2/8);subplot(4,2,5);stem(n,x3,.);axis(0,20,0,1);title(x5(n),p=14,q=8);y3=fft(x3,16);subplot(4,2,6);stem(n,abs(y3),.);axis(0,20,0,4);title(x5(n),p=14,q=8);x4=exp(-(n-15).2/8);subplot(4,2,7);stem(n,x4,.);

16、axis(0,20,0,1);title(x5(n),p=15,q=8);y4=fft(x4,16);subplot(4,2,8);stem(n,abs(y4),.);axis(0,20,0,4);title(fft,p=15,q=8);生成波形为：（3）衰减正弦序列及其不同采样频率时的fft：n=0:50;x1=50*exp(-0.1.*n).*sin(62.5*pi/500*n);subplot(2,2,1);stem(n,x1,.);axis(0,50,-20,40);title(x7(n),fs=1000hz);y1=fft(x1,51);subplot(2,2,2);stem(n,a

17、bs(y1),.);axis(0,50,0,250);title(fft,fs=1000hz);x2=50*exp(-1/3.*n).*sin(62.5*pi/150*n);subplot(2,2,3);stem(n,x2,.);axis(0,50,-20,40);title(x7(n),fs=300hz);y2=fft(x2,51);subplot(2,2,4);stem(n,abs(y2),.);axis(0,50,0,150);title(fft,fs=300hz生成波形为：（4）载波调制解调信号t=0:0.01:3; x=cos(pi*t).*cos(9*pi*t); plot(t,x

18、); axis(0 3 -2 2); grid; title(x5(t)=cos(pi*t)*cos(9*pi*t) t=0:0.01:3; x=cos(pi*t); plot(t,x); axis(0 3 -1 1); grid; title(解调信号)（5）无限长序列n=0:100; x=0.8.n; subplot(1,1,1);stem(n,x,.);axis(0 100 0 1);title(x6(n)=0.8n)n=0:9 x1=0.8.n; subplot(2,2,1);stem(n,x1,.);axis(0,20,0,1);title(序列x6(n),m=10); y1=fft

19、(x1,10); subplot(2,2,2);stem(n,abs(y1),.);axis(0,20,0,3);title(m=10,fft,幅频特性);m=0:19 x2=0.8.m; subplot(2,2,3);stem(m,x2,.);axis(0,20,0,1);title(序列x6(n),m=20); y2=fft(x2,20); subplot(2,2,4);stem(m,abs(y2),.);axis(0,20,0,3);title(m=20,fft,幅频特性)（二）、数字滤波器的设计与实现采用双线性变换法设计iir数字滤波器：源程序：低通：global fpglobal f

20、sglobal apglobal asglobal fs wp=2*pi*fp/fs;ws=2*pi*fs/fs;ts=1/fs;wp1=2/ts*tan(wp/2); ws1=2/ts*tan(ws/2); n,wn=buttord(wp1,ws1,ap,as,s); z,p,k=buttap(n); bap,aap=zp2tf(z,p,k);b,a=lp2lp(bap,aap,wn); bz,az=bilinear(b,a,fs); figure(1)freqz(bz,az); y,t=impz(bz,az,50);figure(2)stem(t,y,.);高通：global fpglob

21、al fsglobal apglobal asglobal fs wp=2*pi*fp/fs;ws=2*pi*fs/fs;ts=1/fs;wp1=2/ts*tan(wp/2); ws1=2/ts*tan(ws/2); n,wn=buttord(wp1,ws1,ap,as,s); z,p,k=buttap(n); bap,aap=zp2tf(z,p,k);b,a=lp2hp(bap,aap,wn); bz,az=bilinear(b,a,fs);figure(1)freqz(bz,az); y,t=impz(bz,az,50);figure(2)stem(t,y,.);带通：global fp1

22、global fp2global fs1global fs2global apglobal asglobal fs wp1=2*pi*fp1/fs;wp2=2*pi*fp2/fs;ws1=2*pi*fs1/fs;ws2=2*pi*fs2/fs;b=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2);ws=1;n,wc=buttord(wp,ws,ap,as,s);num,den=butter(n,wc,s);numt,dent=lp2bp(num,den,w0,b);numd,dend=bilinear(numt,dent,1);w=linspa

23、ce(0,pi,512);figure(1)freqz(numd,dend,w);y,t=impz(numd,dend,50);figure(2)stem(t,y,.)；带阻：global fp1global fp2global fs1global fs2global apglobal asglobal fs wp1=2*pi*fp1/fs;wp2=2*pi*fp2/fs;ws1=2*pi*fs1/fs;ws2=2*pi*fs2/fs;b=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2);ws=1;n,wc=buttord(wp,ws,ap,a

24、s,s);num,den=butter(n,wc,s);numt,dent=lp2bs(num,den,w0,b);numd,dend=bilinear(numt,dent,1);w=linspace(0,pi,512);figure(1)freqz(numd,dend,w);y,t=impz(numd,dend,50);figure(2)stem(t,y,.);生成波形为：低通及其单位脉冲响应：高通及其单位脉冲响应：带通及其单位脉冲序列：带阻及其单位脉冲响应：（三）、综合应用语音信号的分析处理源程序：读入语音信号：x1,fs=wavread(c:tada.wav);sound(x1,fs);

25、fs y1=fft(x1,1024); f=fs*(0:511)/1024;figure(1)plot(x1) title(原始语音信号);xlabel(time n);ylabel(fuzhi n);figure(2)freqz(x1) title(频率响应图)figure(3)subplot(2,1,1);plot(abs(y1(1:512) title(原始语音信号fft频谱)subplot(2,1,2);plot(f,abs(y1(1:512);title(原始语音信号频谱)xlabel(hz);ylabel(fuzhi);语音信号低通处理：global fpglobal fsglob

26、al apglobal asx1,fs=wavread(c:tada.wav);t=0:1/fs:(size(x1)-1)/fs;fswp=2*pi*fp/fs;ws=2*pi*fs/fs;ts=1/fs;wp1=2/ts*tan(wp/2); ws1=2/ts*tan(ws/2); n,wn=buttord(wp1,ws1,ap,as,s); z,p,k=buttap(n); bap,aap=zp2tf(z,p,k);b,a=lp2lp(bap,aap,wn); bz,az=bilinear(b,a,fs); h,w=freqz(bz,az); figure(1)plot(w*fs/(2*p

27、i),abs(h)gridxlabel(频率hz)ylabel(频率响应幅度)title(butterworth)f1=filter(bz,az,x1);figure(2)subplot(2,1,1)plot(t,x1); title(滤波前的时域波形);subplot(2,1,2)plot(t,f1); title(滤波后的时域波形);sound(f1,22050); f0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512); title(滤波前的频谱)xlab

28、el(hz);ylabel(fuzhi);subplot(2,1,2)f1=plot(f,abs(f0(1:512); title(滤波后的频谱)xlabel(hz);ylabel(fuzhi);语音信号高通处理：global fpglobal fsglobal apglobal asx1,fs=wavread(c:tada.wav);t=0:1/fs:(size(x1)-1)/fs;fswp=2*pi*fp/fs;ws=2*pi*fs/fs;ts=1/fs;wp1=2/ts*tan(wp/2); ws1=2/ts*tan(ws/2); n,wn=buttord(wp1,ws1,ap,as,s

29、); z,p,k=buttap(n); bap,aap=zp2tf(z,p,k);b,a=lp2hp(bap,aap,wn); bz,az=bilinear(b,a,fs); h,w=freqz(bz,az); figure(1)plot(w*fs/(2*pi),abs(h)gridxlabel(频率hz)ylabel(频率响应幅度)title(butterworth)f1=filter(bz,az,x1);figure(2)subplot(2,1,1)plot(t,x1); title(滤波前的时域波形);subplot(2,1,2)plot(t,f1); title(滤波后的时域波形);s

30、ound(f1,22050); f0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512); title(滤波前的频谱)xlabel(hz);ylabel(fuzhi);subplot(2,1,2)f1=plot(f,abs(f0(1:512); title(滤波后的频谱)xlabel(hz);ylabel(fuzhi);语音信号带通处理：global fp1global fp2global fs1global fs2global apglobal asx1,fs=

31、wavread(c:tada.wav);t=0:1/fs:(size(x1)-1)/fs;wp1=2*pi*fp1/fs;wp2=2*pi*fp2/fs;ws1=2*pi*fs1/fs;ws2=2*pi*fs2/fs;b=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2);ws=1;n,wc=buttord(wp,ws,ap,as,s);num,den=butter(n,wc,s);numt,dent=lp2bp(num,den,w0,b);numd,dend=bilinear(numt,dent,1);w=linspace(0,pi,512);

32、figure(3)h,w=freqz(numd,dend,w); figure(1)plot(w*fs/(2*pi),abs(h)gridxlabel(频率hz)ylabel(频率响应幅度)title(butterworth)f1=filter(numd,dend,x1);figure(2)subplot(2,1,1)plot(t,x1) title(滤波前的时域波形);subplot(2,1,2)plot(t,f1); title(滤波后的时域波形);sound(f1,22050); f0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1

33、024);subplot(2,1,1);plot(f,abs(y2(1:512); title(滤波前的频谱)xlabel(hz);ylabel(fuzhi);subplot(2,1,2)f1=plot(f,abs(f0(1:512); title(滤波后的频谱)xlabel(hz);ylabel(fuzhi);语音信号带阻处理：global fp1global fp2global fs1global fs2global apglobal asx1,fs=wavread(c:tada.wav);t=0:1/fs:(size(x1)-1)/fs;wp1=2*pi*fp1/fs;wp2=2*pi*

34、fp2/fs;ws1=2*pi*fs1/fs;ws2=2*pi*fs2/fs;b=ws2-ws1;w0=sqrt(ws1*ws2);wp=max(abs(wp1),abs(wp2);ws=1;n,wc=buttord(wp,ws,ap,as,s);num,den=butter(n,wc,s);numt,dent=lp2bs(num,den,w0,b);numd,dend=bilinear(numt,dent,1);w=linspace(0,pi,512);figure(3)h,w=freqz(numd,dend,w); figure(1)plot(w*fs/(2*pi),abs(h)gridx

35、label(频率hz)ylabel(频率响应幅度)title(butterworth)f1=filter(numd,dend,x1);figure(2)subplot(2,1,1)plot(t,x1) title(滤波前的时域波形);subplot(2,1,2)plot(t,f1); title(滤波后的时域波形);sound(f1,22050); f0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512); title(滤波前的频谱)xlabel(hz);ylabel(fuzhi);subplot(2,1,2)f1=plot(f,abs(f0(1:512); title(滤波后的频谱)xlabel(hz);ylabel(fuzhi);生成