- 1、本文档共10页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理
课程设计报告
姓 名:韩 斌
学 号:200712130
专 业:通信工程
院 系:信息工程学院
指导教师:陈黎霞 王玲
一、课题名称 基于matlab的IIR滤波器设计
二、课题意义
数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。数字滤波器和模拟滤波器有着相同的滤波概念,根据其频率响应特性可分为低通、高通、带通、带阻等类型,与模拟滤波器相比,数字滤波器除了具有数字信号处理的固有优点外,还有滤波精度高(与系统字长有关)、稳定性好(仅运行在0与l两个电平状态)、灵活性强等优点。
三、数字滤波器原理
在数字滤波中,我们主要讨论离散时间序列。如图1所示。设输入序列为,离散或数字滤波器对单位抽样序列的响应为。因在时域离散信号和系统中所起的作用相当于单位冲激函数在时域连续信号和系统中所起的作用。
图1 数字滤波器原理
数字滤波器的序列将是这两个序列的离散卷积,即
(1)
同样,两个序列卷积的z变换等于个自z变换的乘积,即
(2)
用代入上式,其中T为抽样周期,则得到
(3)
式中和 分别为数字滤波器输入序列和输出序列的频谱,而为单位抽样序列响应的频谱。由此可见,输入序列的频谱经过滤波后,变为 ,按照的特点和我们处理信号的目的,选取适当的使的滤波后的符合我们的要求。
四、软件设计过程
采集声音信号
语音信号的采样及频谱分析
添加噪声
设计滤波器
语音信号滤波
对滤波前后信号的波形和频谱进行比较
GUI界面(选做)
五、各部分源程序及相应结果
原始信号的采集及分析
源程序:
[y1,fs,bits]=wavread(d:\q.wav);
sound(y1,44000); %播放语音信号
y=fft(y1); %对信号做FFT变换
f=fs*(0:511)/1024;
figure(1)
subplot(2,1,1);
plot(abs(y(1:512))) %做原始语音信号的FFT频谱图
title(原始语音信号FFT频谱)
subplot(2,1,2);%做原始语音信号的时域图形
plot(y1)
title(原始语音信号波形);
xlabel(时间 n);
label(幅值 n);
信号加燥及相关分析
源程序:
[y,Fs,bits]=wavread(D:\q.wav);%读出信号,采样率和采样位数。
y=y(:,1);%取单声道作分析
yl=length(y)%求语音信号长度
yy=fft(y,yl);%傅里叶变换
t=[0:1/8000:4 zeros(1,yl-32001)];
m=0.07*sin(10000*pi*t);%产生噪声
n=y+m;%加入噪声
nl=length(n)%求语音信号长度
nn=fft(n,nl);%傅里叶变换
figure(1);
subplot(2,1,1);plot(n);title(噪声信号波形)
subplot(2,1,2);plot(y);title(原信号波形)
figure(2);
subplot(2,1,1);plot(abs(nn));title(噪声信号频谱);
subplot(2,1,2);plot(abs(yy));title(原信号频谱);
sound(n,fs)
滤波器的设计及相关分析
源程序:
[y,Fs,bits]=wavread(D:\q.wav);%读出信号,采样率和采样位数。
y=y(:,1);%取单声道作分析
yl=length(y)%求语音信号长度
yy=fft(y,yl);%傅里叶变换
t=[0:1/8000:4 zeros(1,yl-32001)];
m=0.07*sin(10000*pi*t);%产生噪声
k=y+m;%加入噪声
kl=length(k)%求语音信号长度
K=fft(k,kl);%傅里叶变换
fb=00;fc=1000;As=100;Ap=1;fs=44000;
wc=2*fc/fs; wb=2*fb/fs;
[n,wn]=ellipord(wc,wb,Ap,As);
[b,a]=ellip(n,Ap,As
您可能关注的文档
最近下载
- XGZ系列刮板输送机图册.pdf
- 第10课《 再塑生命的人》课件(共22张ppt)部编版语文七年级上册.pptx VIP
- 近5年上海市高考语文真题之书序、游记类文言文阅读汇编.docx VIP
- 00708装饰材料与构造-全国自考-2021年4月高等教育自学考试真题(试卷).pdf VIP
- XGZ系列铸石刮板输送机.doc
- 小学篮球原地运球、行进间运球教案.doc
- 中国林业研究院806有机化学考研真题试题2020年.pdf VIP
- 中国林业研究院806有机化学2019考研真题试题.pdf VIP
- 2024年初级咖啡师技能鉴定考试题库大全-上(单选题汇总).docx
- 2024-2030年中国氦(液氦和气态氦)行业市场发展趋势与前景展望战略分析报告.docx
文档评论(0)