- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
[[数字信号处理》课程研究性学习报告
《数字信号处理》课程研究性学习报告
指导教师 薛健
时间 2014.6
【目的】
(1) 掌握IIR和FIR数字滤波器的设计和应用;
(2) 掌握多速率信号处理中的基本概念和方法 ;
(3) 学会用Matlab计算小波分解和重建。
(4)了解小波压缩和去噪的基本原理和方法。
【研讨题目】
一、
(1)播放音频信号 yourn.wav,确定信号的抽样频率,计算信号的频谱,确定噪声信号的频率范围;
(2)设计IIR数字滤波器,滤除音频信号中的噪声。通过实验研究,的选择对滤波效果及滤波器阶数的影响,给出滤波器指标选择的基本原则,确定你认为最合适的滤波器指标。
(3)设计FIR数字滤波器,滤除音频信号中的噪声。与(2)中的IIR数字滤波器,从滤波效果、幅度响应、相位响应、滤波器阶数等方面进行比较。
【设计步骤】
【仿真结果】
【结果分析】
由频谱知噪声频率大于3800Hz。FIR和IIR都可以实现滤波,但从听觉上讲,人对于听觉不如对图像(视觉)明感,没必要要求线性相位,因此,综合来看选IIR滤波器好一点,因为在同等要求下,IIR滤波器阶数可以做的很低而FIR滤波器阶数太高,自身线性相位的良好特性在此处用处不大。
【自主学习内容】
MATLAB滤波器设计
【阅读文献】
老师课件,教材
【发现问题】 (专题研讨或相关知识点学习中发现的问题):
过渡带的宽度会影响滤波器阶数N
【问题探究】
通过实验,但过渡带越宽时,N越小,滤波器阶数越低,过渡带越窄反之。这与理论相符合。
【仿真程序】
信号初步处理部分:
[x1,Fs,bits] = wavread(yourn.wav);
sound(x1,Fs);
y1=fft(x1,1024);
f=Fs*(0:511)/1024;
figure(1)
plot(x1)
title(原始语音信号时域图谱);
xlabel(time n);
ylabel(magnitude 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(magnitude);
IIR:
fp=2500;fs=3500;
wp = 2*pi*fp/FS;
ws = 2*pi*fs/FS;
Rp=1;
Rs=15;
Ts=1/Fs;
wp = 2*pi*fp/FS;
ws = 2*pi*fs/FS;
wp1=2/Ts*tan(wp/2);
ws1=2/Ts*tan(ws/2);
t=0:1/11000:(size(x1)-1)/11000;Au=0.03;d=[Au*cos(2*pi*5000*t)];x2=x1+d;
[N,Wn]=buttord(wp1,ws1,Rp,Rs,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(4)
plot(W*Fs/(2*pi),abs(H))
grid
xlabel(频率/Hz)
ylabel(频率响应幅度)
title(Butterworth)
f1=filter(bz,az,x2);
figure(5)
subplot(2,1,1)
plot(t,x2)
title(滤波前时域波形);
subplot(2,1,2)
plot(t,f1);
title(滤波后时域波形);
sound(f1,FS);
FIR
[x1,Fs,bits] = wavread(I:/dsp_2014_project3/yourn);
fp=2500;fs=3500;
wp = 2*pi*fp/Fs;
w
文档评论(0)