没有合适的资源?快使用搜索试试~ 我知道了~
首页数字信号处理课程设计.doc
资源详情
资源评论
资源推荐

数字信号处理实践与设计
封面自己打!
目 录
课题一 语音信号中的抽取、混叠和数字限带滤波......................................................1
1 需求分析与方案设计.............................................................................................1
1.1 需求分析.......................................................................................................1
1.2 方案设计.......................................................................................................1
1.2.1 设计原理.............................................................................................1
1.2.2 方案构思.............................................................................................3
2 模块设计及原理分析.............................................................................................5
2.1 各个函数模块设计.......................................................................................5
2.1.1 作频谱模块.........................................................................................5
2.1.2 音频数据的读取模块.........................................................................5
2.1.3.resample 函数抽取模块......................................................................5
2.1.4 直接抽取模块.....................................................................................5
2.1.5 抽取滤波器设计模块.........................................................................5
2.1.6 卷积滤波抽取模块.............................................................................6
2.1.7 多项滤波模块.....................................................................................6
2.1.8 作出四个时域波形模块.....................................................................6
2.1.9 计算对应指标模块.............................................................................7
2.1.10.卷积滤波插值模块...........................................................................7
2.1.11 多相滤波插值模块............................................................................7
2.2 主函数模块...................................................................................................7
2.2.1 主函数调用子函数实现基本要求.....................................................7
2.2.2 主函数调用子函数实现提高要求.....................................................7
3 联合调试与总体测试.............................................................................................9
3.1 原音频信号的频谱.......................................................................................9
3.2 四种抽取方式的时频域的分析...................................................................9
3.3 有限长滤波器设计.......................................................................................9
3.4 技术指标分析...............................................................................................9
3.5 插值.............................................................................................................10
3.6 直接卷积与多项滤波速率转换..................................................................11
心 得 体 会......................................................................................................................12
参 考 文 献......................................................................................................................13
.........................................................................................................................................14
Ⅰ

1
数字信号处理实践与设计
课题一 语音信号中的抽取、混叠和数字限带滤波
1 需求分析与方案设计
1.1 需求分析
在数字音频系统中,广播,数字压缩光盘和数字音频磁带就使用了不同的采
样率。音频系统经常需要在几种不同的采样率之间转换,多速率处理系统中含有
两类操作,即抽取和插值,这就需要不同的插值与抽取以适应多速率转换。抽取
是在原序列的基础上抽取点数,降低采样率;抽取是在原序列的基础上插入点数,
提高采样率。而抽取和插值结合就能实现速率的分数倍转换。多速率信号处理系
统已广泛应用于音频信号处理、多媒体通信等许多领域,在实际中会采取大量的
采样率转换。因此,关于采样率转换的理论与技术研究也显得越来越重要,建立
在采样率转换基础上的多速率信号处理已成为数字信号处理学科中的重要内容之
一。
基于上述实际应用以及设计指标要求,我们的设计需要满足如下需求:
1.能够实现对音频信号的抽取,并能控制抽取的起止时间,而且能够不失真还原
原始音频信号。利用 MATLAB 的读取音频文件的函数能够从音频文件读取
44100Hz 采样频率的 x(n),可以通过参数设置读取的起始时间和持续时间,即能
够获取音频文件的时域序列。
2.能够实现对音频信息的不同抽取方式,本次设计需要设计四种抽取方式:调用
MATLAB resample 函数对 x(n)进行抽取得到 y1(m) 、直接对 x(n)进行抽取得到
y2(m)、先卷积滤波然后抽取得到 y3(m)、先多相滤波然后抽取得到 y4(m)。滤波
器需要自主设计,能够满足抽样抗混叠滤波。
3.能够做出抽取信号的时域波形与频谱图,分析和对比不同帧片段的 y1(m)、
y2(m)、 y3(m)、y4(m)的时域和频域信息,即能够观察当抽取因子 D 在任意参数
下,四种方式抽取信号的混叠情况;并能够对四种方式抽取的技术指标进行分析,
即对比计算信噪比、运算速率(MATLAB 中用计算时间来呈现)技术参数,并满
足 D=30 下信噪比 SNR3-SNR2>30dB,|SNR3-SNR4|<1dB 的指标要求。
4.能够实现信号的任意 I 插值,并实现多相滤波结构完成速率变换和分数倍采样
频率变换,即能够实现抽样插值定理,自主设计插值-抽样中的滤波器,完成分数
倍转换;完成多项滤波结构,实现高采样率下的抽样与插值。
1.2 方案设计
1.2.1 设计原理
1.整数因子的抽取:满足关系式 x
D
=(m)=x(mD),如果采样速率为 fs,那得
到的序列
x
D
(n)的采样率为 fs/D,抽取后信号不产生混叠的频率范围变为 fm<fs/2D。而
抽样序列的频谱为

2
数字信号处理实践与设计
图 1.1 抽取后混叠示意图
可以看到抽取序列的频谱是原序列频谱 D 倍展宽后,按 2π 的整数倍移位并
叠加而成,所以需求中当 D=30,采样频率变为 1470Hz,小于音频信号的最高
频率,必然会发生失真。为保证抽取后不丢失信息或减少丢失信息,应在抽取
前先用一个数字低通滤波器进行滤波,对应理想数字滤波器的幅度响应应满足
|H
D
(e
jw
)| = D, (|w|<π/D); 0,(π/D≤|w|≤π)
而在实际应用中,我们无法实现理想数字低通滤波器,所以我们需要利用窗函数
法,自主设计 FIR 数字低通滤波器,使其满足对应的幅度响应,达到理想低通滤
波器的效果。
2. 整数因子插值:满足 x
I
(m)=x(m/I),(m=0,±I,±2I...);0,(其他),
得到频谱 X(e
jwI
)
内插后信号频谱被压缩了 I 倍 ,需要通过低通滤波器消除内插带来的镜像,滤波
器满足
|H
I
(e
jw
)| = 1, (|w|<π/I); 0,(π/I≤|w|≤π)
3.采样率的分数倍转换:即通过 D 抽取和 I 插值结合来实现,为了保证不丢失信
息,先插值再抽取,考虑到抽样与插值的处理过程均需要使用低通滤波器,所以
如图分析可以将两个低通滤波器合二为一,低通截止频率取二者之间的较小值。
图 1.2 采样率分数倍转换原理图
4.多项滤波结构
抽取和插值是多速率信号处理系统中的两个最基本的运算,当对运算速度要

3
数字信号处理实践与设计
求相当高。但在抽取模型中,低通滤波器在减速之前实现的,而插值模型,低通
滤波器又位于提速之后,对实时处理不利。所以,为了提高运算速率,抽取利用
如图 1.3 的多项滤波器结构,数字滤波器位于抽取器之后,即滤波是降速后进行
的,降低了处理速度的要求,提升实时处理能力。
图 1.3 抽取的多项滤波器结构
同理,对于插值的多项滤波器结构使滤波器位于内插器前,使滤波在提速前,
降低了对处理速度的要求的同时,也提高了运算精度和降低对字长的要求。
1.2.2 方案构思
1.基础要求设计构思
根据设计需求,首先完成 44100Hz 可调节起止区间的音频信号采样,为了观察
到频谱,可以利用 audioread 函数读取数据,并利用 fft 函数得到频谱,fftshift 进
行频谱搬移,并用进行 plot 绘图。
其次完成四种抽取方式:
(1)调用 MATLAB 自带的 resample()函数进行抽取,并绘出时域波形,用
fft、fftshift、plot 函数绘出频谱。
(2)直接按步长抽取,利用 a:b:c 的形式抽取,以 b 为间隔,从 a-c 中抽
取,并绘出时域波形,用 fft、fftshift、plot 函数绘出频谱。
(3)卷积滤波后抽取,利用卷积函数 conv,时域音频信号与滤波器单位冲
激函数 h(n)卷积后抽取,并绘出时域波形,用 fft、fftshift、plot 函数绘出频谱。
滤波器设计要达到对应设计原理中的要求,本次设计中用窗函数法设计有限长数
字滤波器。
(4)多相滤波后抽取,滤波器设计同上,利用 floor 函数以及 reshape 函数实
现多项滤波操作,并绘出时域波形,用 fft、fftshift、plot 函数绘出频谱。
这样可以对比抽取绘制其时域波形以及频谱得到对应的时域波形以及频谱,
通过 tic,toc 函数可以得到对应的计算时间,来分析抽取效率。
计算信噪比则根据公式:
剩余15页未读,继续阅读










Green_endless
- 粉丝: 11
- 资源: 5
上传资源 快速赚钱
我的内容管理 收起
我的资源 快来上传第一个资源
我的收益
登录查看自己的收益我的积分 登录查看自己的积分
我的C币 登录后查看C币余额
我的收藏
我的下载
下载帮助

会员权益专享
最新资源
- ARM Cortex-A(armV7)编程手册V4.0.pdf
- ABB机器人保养总结解析.ppt
- 【超详细图解】菜鸡如何理解双向链表的python代码实现
- 常用网络命令的使用 ipconfig ping ARP FTP Netstat Route Tftp Tracert Telnet nslookup
- 基于单片机控制的DC-DC变换电路
- RS-232接口电路的ESD保护.pdf
- linux下用time(NULL)函数和localtime()获取当前时间的方法
- Openstack用户使用手册.docx
- KUKA KR 30 hA,KR 60 hA机器人产品手册.pdf
- Java programming with JNI
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈



安全验证
文档复制为VIP权益,开通VIP直接复制

评论0