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

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

指导书

适用专业： 通信工程、电子信息工程

教师姓名： 余尤好

莆田学院电信系电信教研室

2011.02

目 录

设计一 正余弦信号的谱分析……………………………………………………………………1 设计二 数字滤波器的设计及实现………………………………………………………………3 设计三 语音信号滤波处理………………………………………………………………………6 设计四 调制解调系统的设计及实现……………………………………………………………7 设计五 倒频系统的软件实现及使用…………………………………………………………11

设计一 正余弦信号的谱分析

【一】 设计目的

1. 用DFT实现对正余弦信号的谱分析； 2. 观察DFT长度和窗函数长度对频谱的影响； 3. 对DFT进行谱分析中的误差现象获得感性认识。

【二】 设计原理

一、谱分析原理

数字信号处理方法的一个重要用途是在离散时间域中确定一个连续时间信号的频谱，通常称为频谱分析，更具体地说，它也包括确定能量谱和功率谱。数字频谱分析可以使用在很广阔领域，频谱分析方法是基于以下的观测：如果连续时间信号ga(t)是频带有限的，那么对其离散时间等效信号g(n)的DFT进行谱分析。然而，在大多数情况下，ga(t)是在

???t??范围内定义的，因此g(n)也就定义在???n??的无线范围内，要估计一个

无限长信号的频谱是不可能的。实用的方法是：先让模拟连续信号ga(t)通过一个抗混叠的模拟滤波器，然后把它采样成一个离散序列g(n)。假定反混叠滤波器的设计是正确的，则混叠效应可以忽略，又假设A/D变换器的字长足够长，则A/D变换中的量化噪声也可忽略。

假定表征正余弦信号的基本参数，如振幅、频率和相位不随时间改变，则此信号的傅立叶变换G(ej?)可以用计算它的DTFT得到

j? G(e)?n????g(n)e??j?n （1.1）

实际上无限长序列g(n)首先乘以一个长度为M的窗函数w(n)，使它变成一个长为M的有限长序列，g1(n)?g(n)w(n)，对g1(n)求出的DTFT G1(e模拟信号ga(t)的频谱估计，然后求出G1(ej?j?)应该可以作为原连续

)在0???2?区间等分为N点的离散傅立

叶变换DFT。为保证足够的分辨率，DFT的长度N选的比窗长度M大，其方法是在截断了的序列后面补上N－M个零。计算采用FFT算法。 二、MATLAB函数介绍

1. 输入函数input( ) 格式：R=input(string)

功能：在屏幕上显示input括号后的’string’内容，提示用户从键盘输入某值，并将输入的值赋给R。

例如，在命令窗口输入R=input(‘How many apples’) 会显示How many apples 从键盘输入 3

会显示 R=3

2. 一维快速傅里叶变换函数fft( ) 格式： y=fft(x) y=fft(x,n)

说明：fft函数用于计算矢量或矩阵的傅里叶变换。 ３．幅度函数abs( ) 格式：abs（X）

功能：对X取绝对值，当X是复数时，得到X的复模值。

例1-1 对余弦序列x(n)?cos(?n/4)进行8点谱分析，绘出序列和幅频特性曲线。

MATLAB命令如下：

N=input('输入谱分析的长度') n=0:N-1; x=cos(pi*n/4);

subplot(1,2,1),stem(n,x) %绘制序列的波形 xlabel('n');ylabel('x1(n)'); title('余弦序列');

X=abs(fft(x,N)); %求余弦序列的N点FFT subplot(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，其他参数不变，重新画出该序列的频谱图，观察频率泄漏现象，分析原因。

2. 考察DFT的长度对双频率信号频谱分析的影响。设待分析的信号为 x(n)?0.5sin(2?f1n)?sin(2?f2n)0?n?15 （1.2）

令两个长度为16的正余弦序列的数字频率为f1?0.22及f2?0.34。取N为四个不同值16，32，64，128。画出四个DFT幅频图，分析DFT长度对频谱分辨率的影响。

3. 在上题中若把两个正弦波的频率取得较近，令f1?0.22，f2?0.25，试问怎样选

择FFT参数才能在频谱分析中分辨出这两个分量？

【四】课程设计报告要求

完成课程设计任务后，应按要求提交课程设计报告。设计报告应包含如下几个方面的内容。

1. 课程设计目的 2. 课程设计要求。 3. 详细设计过程。 4. 调试分析。 5. 结果分析和体会。 6. 附录或参考资料。

设计二 数字滤波器的设计及实现

【一】设计目的

1. 熟悉IIR数字滤波器和FIR数字滤波器的设计原理和方法；

2. 学会调用MATLAB信号处理工具箱中的滤波器设计函数设计各种IIR和FIR数字

滤波器，学会根据滤波要求确定滤波器指标参数；

3. 掌握用IIR和FIR数字滤波器的MATLAB实现方法，并能绘制滤波器的幅频特性、

相频特性；

4. 通过观察滤波器的输入、输出信号的时域波形及其频谱，建立数字滤波的概念。

【二】设计原理

抑制载波单频调幅信号的数学表达式为 s(t)?cos(2?f0t)cos(2?fct)?1?cos(2?(fc?f0)t)?cos(2?(fc?f0)t)? （2.1） 2

本文来源：https://www.bwwdw.com/article/ja3x.html