《数字信号处理》课程设计(2013)

1、 数字信号处理课程设计 指导书适用专业： 通信工程、电子信息工程 教师姓名： 余尤好 莆田学院电信系电信教研室2011.02目 录设计一 正余弦信号的谱分析1设计二 数字滤波器的设计及实现3设计三 语音信号滤波处理6设计四 调制解调系统的设计及实现7设计五 倒频系统的软件实现及应用11数字信号处理课程设计 莆田学院电信系设计一 正余弦信号的谱分析【一】 设计目的1. 用DFT实现对正余弦信号的谱分析；2. 观察DFT长度和窗函数长度对频谱的影响；3. 对DFT进行谱分析中的误差现象获得感性认识。【二】 设计原理一、谱分析原理数字信号处理方法的一个重要用途是在离散时间域中确定一个连续时间信号的频

2、谱，通常称为频谱分析，更具体地说，它也包括确定能量谱和功率谱。数字频谱分析可以应用在很广阔领域，频谱分析方法是基于以下的观测：如果连续时间信号是频带有限的，那么对其离散时间等效信号的DFT进行谱分析。然而，在大多数情况下，是在范围内定义的，因此也就定义在的无线范围内，要估计一个无限长信号的频谱是不可能的。实用的方法是：先让模拟连续信号通过一个抗混叠的模拟滤波器，然后把它采样成一个离散序列。假定反混叠滤波器的设计是正确的，则混叠效应可以忽略，又假设A/D变换器的字长足够长，则A/D变换中的量化噪声也可忽略。假定表征正余弦信号的基本参数，如振幅、频率和相位不随时间改变，则此信号的傅立叶变换可以用计

3、算它的DTFT得到 （1.1）实际上无限长序列首先乘以一个长度为M的窗函数，使它变成一个长为M的有限长序列，对求出的DTFT 应该可以作为原连续模拟信号的频谱估计，然后求出在区间等分为N点的离散傅立叶变换DFT。为保证足够的分辨率，DFT的长度N选的比窗长度M大，其方法是在截断了的序列后面补上NM个零。计算采用FFT算法。二、MATLAB函数介绍1. 输入函数input( ) 格式：R=input(string) 功能：在屏幕上显示input括号后的string内容，提示用户从键盘输入某值，并将输入的值赋给R。例如，在命令窗口输入R=input(How many apples) 会显示How

4、many apples 从键盘输入 3 会显示 R=32. 一维快速傅里叶变换函数fft( ) 格式： y=fft(x) y=fft(x,n) 说明：fft函数用于计算矢量或矩阵的傅里叶变换。幅度函数abs() 格式：abs（X） 功能：对X取绝对值，当X是复数时，得到X的复模值。 例1-1 对余弦序列进行8点谱分析，绘出序列和幅频特性曲线。 MATLAB命令如下：N=input('输入谱分析的长度')n=0:N-1;x=cos(pi*n/4);subplot(1,2,1),stem(n,x) %绘制序列的波形xlabel('n');ylabel('x1

5、(n)'); title('余弦序列');X=abs(fft(x,N); %求余弦序列的N点FFTsubplot(1,2,2)k=0:N-1;stem(k,X) %绘制序列的幅频特性曲线xlabel('k');ylabel('X(k)'); string=num2str(N),'点FFT幅频曲线'title(string);程序运行结果如图1.1所示。 图1.1 例1-1的运行结果【三】 设计内容1. 对一个频率为10Hz，采样频率为64Hz的32点余弦序列进行谱分析，画出其频谱图；若将频率改为11Hz，其他参数不变，重新

6、画出该序列的频谱图，观察频率泄漏现象，分析原因。2. 考察DFT的长度对双频率信号频谱分析的影响。设待分析的信号为 （1.2） 令两个长度为16的正余弦序列的数字频率为及。取N为四个不同值16，32，64，128。画出四个DFT幅频图，分析DFT长度对频谱分辨率的影响。3. 在上题中若把两个正弦波的频率取得较近，令，试问怎样选择FFT参数才能在频谱分析中分辨出这两个分量？【四】课程设计报告要求完成课程设计任务后，应按要求提交课程设计报告。设计报告应包含如下几个方面的内容。 1. 课程设计目的2. 课程设计要求。3. 详细设计过程。4. 调试分析。5. 结果分析与体会。6. 附录或参考资料。设计

7、二 数字滤波器的设计及实现【一】设计目的1. 熟悉IIR数字滤波器和FIR数字滤波器的设计原理和方法；2. 学会调用MATLAB信号处理工具箱中的滤波器设计函数设计各种IIR和FIR数字滤波器，学会根据滤波要求确定滤波器指标参数；3. 掌握用IIR和FIR数字滤波器的MATLAB实现方法，并能绘制滤波器的幅频特性、相频特性；4. 通过观察滤波器的输入、输出信号的时域波形及其频谱，建立数字滤波的概念。【二】设计原理 抑制载波单频调幅信号的数学表达式为 （2.1）其中，称为载波，为载波频率，称为单频调制信号，为调制正弦波信号频率，且满足>。由（2.1）式可见，所谓抑制载波单频调制信号，就是两

8、个正弦信号相乘，它有2个频率成分：和频+，差频-，这两个频率成分关于载波频率对称。所以，1路抑制载波单频调幅信号的频谱图是关于载波频率对称的两根谱线。显然，当调制频率和（或）载波频率不同时，可以得到包含不同频率成分的单频调幅信号，将几路不同频率成分的单频调幅信号相加后形成混合信号，产生复合信号的函数mstg程序范例如下：function st=mstg%产生信号序列st，并显示st的时域波形和频谱%st=mstg返回三路调幅信号相加形成的混合信号，长度N=800N=800; %信号长度N为800Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz，Tp为采样时间t=0:T

9、:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; %第1路调幅信号载波频率fc1=1000Hzfm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hzfc2=Fs/20; %第2路调幅信号载波频率fc2=500Hzfm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hzfc3=Fs/40; %第3路调幅信号载波频率fc3=250Hzfm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hzxt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第1路调幅信号xt2=cos(2*pi*fm2*t).*cos

10、(2*pi*fc2*t); %产生第2路调幅信号xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号st=xt1+xt2+xt3; %三路信号相加，得到复合信号fxt=fft(st,N); %计算信号st的频谱%以下为绘图命令subplot(2,1,1);plot(t,st);grid;xlabel('t/s');ylabel('s(t)');axis(0,Tp,min(st),max(st);title('(a)s(t)的波形')subplot(2,1,2);stem(f,abs(fxt)/max(a

11、bs(fxt),'.');grid;title('(b)s(t)的频谱')axis(0,Fs/8,0,1.2);xlabel('f/Hz');ylabel('幅度');该函数产生由三路已知载波调幅信号相加构成的复合信号st，并绘图显示st的时域波形和幅频特性曲线如图2.1所示。(a) s(t)的波形(b) s(t)的频谱图2.1 三路调幅信号st的时域波形和幅频特性曲线 由图2.1可见，三路信号时域混叠无法在时域进行分离，但频域是分离的。容易看出，这三路调幅信号的载波频率分别为250Hz、500Hz和1000Hz，因此可以通过设计

12、合适的滤波器的方法在频域分离，这就是本设计的目的。【三】设计内容1. 调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，观察st的时域波形和幅频特性曲线；2. 通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率；假定要求滤波器的通带最大衰减为0.1dB，阻带最小衰减为60dB，编程调用MATLAB滤波器设计函数分别设计这三个数字滤波器，并绘图显示其幅频特性曲线。3. 用所设计的三个滤波器分别对复合信号st进行滤波，分离出st中的三路不同载波频率的调幅信号，并绘图显示滤波

13、后信号的时域波形和频谱，观察分离效果。【四】课程设计报告要求完成课程设计任务后，应按要求提交课程设计报告。设计报告应包含如下几个方面的内容。 1. 课程设计目的2. 课程设计要求。3. 详细设计过程。4. 调试分析。5. 结果分析与体会。6. 附录或参考资料。设计三 语音信号滤波处理【一】设计目的1. 了解语音信号的产生、采集，能绘制语音信号的频率响应曲线及频谱图；2. 学会用MATLAB对语音信号进行分析和处理；3. 掌握用滤波器去除语音信号噪声的方法，观察去噪前后的语音信号。【二】设计原理1. 语音信号的采集在MATLAB软件平台下，利用函数wavread（ ）对语音信号采集，并记录采样频

14、率和采样点数。将语音信号转换成计算机能够运算的有限长序列。wavread函数的调用格式如下：y=wavread(file)读取file所规定的wav文件，返回采样值放在向量y中。y,fs=wavread(file)采样值放在向量y中，fs表示采样频率（Hz）。 2. 用FFT作谱分析FFT即快速傅立叶变换，它是从DFT运算中发展起来的，利用系数 的对称性和周期性减少运算量。长度为N的序列直接计算DFT需要次复乘和次复加，而用FFT进行运算一般需要次复乘和次复加，从而使DFT的运算大大简化。用FFT对连续信号进行谱分析的步骤如下：Sc(t)LPFA/DDFT 图3.1 连续信号谱分析过程引入前置

15、低通滤波器LPF是为了消除或减少时域连续信号转换成序列时可能出现的频谱混叠现象。表示时域有限的窗函数。3.设计滤波器去除语音信号的噪声 通过wavread（）函数将语音信号读入，通过频率采样及fft ( )产生信号，并对之加噪，通过窗函数法设计滤波器滤掉该语音信号的噪声，对比滤波前后的语音波形和频谱。【三】设计内容1. 利用Windows下的录音机录制一段自己的话音，时间在1s内。然后在Matlab软件平台下，利用wavread函数对语音信号进行采样，记住采样频率和采样点数；2. 画出语音信号的时域波形，对采样后的语音信号进行快速傅立叶变换，得到语音信号的频谱特性；对语音信号分别加入正弦噪声和

16、白噪声，画出加噪信号的时域波形和频谱图；3. 根据对加噪语音信号谱分析的结果，确定滤除噪声所需滤波器的技术指标，设计合适的数字滤波器，并画出滤波器的频域响应；4. 用所设计的滤波器对采集的信号进行滤波，在同一个窗口画出滤波前后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；利用sound(x)回放语音信号，验证设计效果；5. *为使编制的程序操作方便，设计一个信号处理用户界面。在所设计的系统界面上可以选择滤波器的类型，输入滤波器的参数，显示滤波器的频率响应，选择信号等。【四】课程设计报告要求完成课程设计任务后，应按要求提交课程设计报告。设计报告应包含如下几个方面的内容。 1.

17、 课程设计目的2. 课程设计要求。3. 详细设计过程。4. 调试分析。5. 结果分析与体会。6. 附录或参考资料。设计四 调制解调系统的设计及实现【一】设计目的1. 通过该设计，掌握调制解调系统的原理及编程实现的方法。2. 在掌握相关知识的基础上，学会自己设计实验，分析验证不同振幅调制方式的优缺点，提高进行信号分析和处理的能力。【二】设计原理两个信号在时域的乘法运算通常用来实现信号的调制，即由一个信号去控制另一个信号的某一个参量。例如用一个低频的正弦波信号去控制另一个频率较高的正弦波信号的幅值，则产生一个振幅调制信号，MATLAB程序范例如下：%用频率为Fc,采样频率为Fs的载波s(t)调制原

18、信号x(t)%注意Fs>2*Fc+BW,其中BW为原信号带宽,本例中BW=Fm%观测信号的持续时间=1sclear;clf;Fm=10;Fc=100;Fs=1000;N=1000;k=0:N-1;t=k/Fs;x=sin(2.0*pi*Fm*t);xf=abs(fft(x,N);s=sin(2.0*pi*Fc*t);sf=abs(fft(s,N);y=x.*s;yf=abs(fft(y,N);subplot(2,1,1);plot(t,y);xlabel('时间(s)');title('已调信号y(t)的波形和频谱');subplot(2,1,2);ste

19、m(yf);xlabel('频率(H)');信号的调制在通信领域应用非常广泛，下面介绍几种常用的振幅调制方式。一常规的双边带调幅（AM）已调信号的时域表达式为：式中外加的直流分量；调制信号（也称未调信号）；调制系数，它的范围在（0，1）之间；载波的频率；载波的初始相位。从频谱来看，调制过程就是将调制信号的频谱由低频搬移到高频，令为调制信号的频谱，则已调信号的频谱为由频域表达式可见，已调信号的频带宽度是调制信号的频带宽度的两倍，同时，已调信号的频谱中包含有正弦载波分量，即有部分功率消耗在载波上，没有用于信息的传递，因此常规双边带调幅方式效率较低。二抑制载波双边带调幅（DSB）DS

20、B是另一种常用的振幅调制方式，其时域和频域的表达式为：可见，已调信号的频带宽度仍是调制信号的频带宽度的两倍，但已调信号中不包含载波分量，因此提高了调制效率。三单边带调幅（SSB）由于DSB信号的上、下两个边带是完全对称的，皆携带了调制信号的全部信息，因此，从信息传输的角度来考虑，仅传输其中一个边带就够了。这就又演变出另一种新的调制方式单边带调制（SSB）。产生SSB信号的方法很多，其中最基本的方法是滤波法，系统框图如图4.1所示。图中的为单边带滤波器，将DSB信号通过单边带滤波器得到单边带已调信号。单边带调幅有上边带(USB)和下边带(LSB)两种方式，对应的频域表达式为： 可见，已调信号的频

21、带宽度就是调制信号的频带宽度，因此占用频带资源相对较小。图4.1产生单边带信号的系统框图用滤波法形成SSB信号，原理框图简洁、直观，但存在的一个重要问题是单边带滤波器不易制作。这是因为，理想特性的滤波器是不可能做到的，实际滤波器从通带到阻带总有一个过渡带。滤波器的实现难度与过渡带相对于载频的归一化值有关，过渡带的归一化值愈小，分割上、下边带就愈难实现。而一般调制信号都具有丰富的低频成分，经过调制后得到的DSB信号的上、下边带之间的间隔很窄，要想通过一个边带而滤除另一个，要求单边带滤波器在附近具有陡峭的截止特性即很小的过渡带，这就使得滤波器的设计与制作很困难，有时甚至难以实现。为此，实际中往往采

22、用多级调制的办法，目的在于降低每一级的过渡带归一化值，减小实现难度。限于篇幅，本节不作详细介绍。四残留边带调制（VSB）残留边带调制是介于单边带调制与双边带调制之间的一种调制方式，它既克服了DSB信号占用频带宽的问题，又解决了单边带滤波器不易实现的难题。用滤波法实现残留边带调制的原理图如图4.2所示。图4.2产生残留边带信号的系统框图图中的为残留边带滤波器，为了保证相干解调时无失真地得到调制信号，残留边带滤波器的传输函数必须满足它的几何含义是，残留边带滤波器的传输函数在载频附近必须具有互补对称性。图4.3示出的是满足该条件的典型实例：残留部分上边带时滤波器的传递函数如图4.3（a）所示，残留部

23、分下边带时滤波器的传递函数如图4.3（b）所示。图4.3残留边带滤波器特性图4.3所示的滤波器，可以看作是对截止频率为的理想滤波器的进行“平滑”的结果，习惯上，称这种“平滑”为“滚降”。显然，由于“滚降”，滤波器截止频率特性的“陡度”变缓，实现难度降低，但滤波器的带宽变宽。在残留边带调制中，除了传送一个边带外，还保留了另外一个边带的一部分。对于具有低频及直流分量的调制信号，用滤波法实现单边带调制时所需要的过渡带无限陡的理想滤波器，在残留边带调制中已不再需要，这就避免了实现上的困难。由于VSB基本性能接近SSB，而VSB调制中的边带滤波器比SSB中的边带滤波器容易实现，所以VSB调制在广播电视、

24、通信等系统中得到广泛应用。振幅解调的原理振幅解调是振幅调制的逆过程，通常称为检波。它的作用是从已调制的高频振荡中恢复出原来的调制信号。检波过程与调制过程正好相反。从频谱来看，检波就是将已调信号的频谱由高频搬移到低频，然后通过滤波器，滤除无用频率分量，取出所需要的原调制信号。检波器可分为包络检波和同步检波两大类。当AM振幅调制信号的包络直接反映了调制信号的变化规律时，可以用二极管包络检波的方法进行解调。而抑制载波的双边带或单边带振幅调制信号的包络不能直接反映调制信号的变化规律，无法用包络检波进行解调，所以采用同步检波方法。即将已调信号乘上与发射端同频同相载波，进行频率变换，经低通滤波恢复出调制信

25、号。【三】设计内容1. 根据范例中给出的调制信号和载波，对已调信号进行同步检波（相干解调），设计低通滤波器提取原信号，画出解调信号的波形和频谱图，验证设计效果。2. 一有限长度的信号，其时域表达式为用常规双边带调幅将其调制在载波上，假设，调制系数，分别画出原信号和已调信号的时域波形图和频谱图。3. 信号，令，载波，用抑制载波双边带调幅来调制信号；用单边带调幅来调制信号（分别采用上边带和下边带两种方式）；对已调信号进行同步检波（解调）。分别画出原信号已调信号解调信号的时域波形和频谱，验证设计效果。思考题：观察123题中各信号的包络，解释说明如何实现包络检波。根据实验结果分析说明各种振幅调制方式的优缺点。*设计滤波器实现VSB调制系统并验证设计效果。【四】课程设计报告要求完成课程设计任务后，应按要求提交课程设计报告。设计报告应包含如下几个方面的内容。1. 课程设计目的2. 课程设计要求。3. 详细设计过程。4. 调试分析。5. 结果分析与体会。6. 附录或参考资料。设计五 倒频系统的软件实现及