《数字信号处理大作业.pdf》由会员分享，可在线阅读，更多相关《数字信号处理大作业.pdf（6页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、数字信号处理课题设计基于 Matlab 的音频采样姓姓名：名：谌谌 海海 龙龙学学号：号：指导教师：指导教师：班班级：级：1313 电子电子 4 4 班班成成绩：绩：重庆大学通信工程学院2015 年 11 月基于基于 MATLABMATLAB的音乐采样实验的音乐采样实验一、实验内容及原理一、实验内容及原理内容内容：1、用 fs=44,100Hz 采集一段音乐；2、改变采样率，用 fs=5,512Hz 采集一段音乐，体会混叠现象；3、录制一段自己的声音，试验当 fs=？时，发生混叠？4、在噪声环境中录制一段自己的声音，试采用一种方法将噪声尽可能地消除。原理原理：根据奈奎斯特采样定律，如果连续时间
2、信号xat是最高截止频率为m的带限at通过一个增益为 T，截止频信号，采样频率s 2m，那么让采样信号x率为s/2的理想低通滤波器，可以无失真地恢复出原连续时间信号xat。否则，s 2m会造成采样信号中频谱混叠现象，不能无失真地恢复连续时间信号。a a t t x xa a t t x x采样脉冲p p t t 采样过程X(j)10j)P(amm连续信号采样信号 2ss0s s 2ss0(j)Xas2s1/T0smXa(j)1/Tss2s3s采样信号的频谱混叠x xa a t t 1X Xa a j j op(t)(1)t mmo mm P P j j Et相乘 s s oo a a t t
3、x xT s s s s X Xa a1 1T T 卷积oTt s so mm s s 采样信号的频谱人发出的声音频率在 3003400Hz，可以通过加低通滤波器消除高频噪声。二、实验过程二、实验过程脚本文件 ffts.m 的功能是画 y 的频谱图function ffts(y,fs)NFFT1=2nextpow2(length(y);fy=fft(y,NFFT1)/length(y);p=2*abs(fy(1:NFFT1);f=fs/2*linspace(0,1,NFFT1);plot(f,p);1、用 f s=44,100Hz 采集一段音乐；脚本文件：work1.mclear;clc;cl
4、ose all;y,fs,nbits=wavread(work1.wav);y1=y(1:600000);ffts(y1,fs);axis(0 23000 0 0.02);title(用44100Hz抽样work1.wav的频谱);sound(y1,44100);wavwrite(y1,44100,work1_1.wav);分析：根据采样定理，当采样频率fs2fm时，会产生频率混叠现象，听起来的声音会有种很多时刻的声音重叠的效果。由于人耳所能听到的声音在20Hz到20000Hz之间，因此当采样频率为44100Hz时刚好满足采样定理所要求的条件fs2fm，所以这也验证了抽样定理中对于采样频率的要
5、求fs2fm，其中fm为信号的频谱中最高频率分量2、改变采样率，用 fs=5,512Hz 采集一段音乐，体会混叠现象；脚本文件：work2.mclear;clc;close all;y,fs,nbits=wavread(work1.wav);y1=y(1:8:600000);ffts(y1,5512);axis(0 3000 0 0.02);title(用 5512Hz 抽样 work1 的频谱);sound(y1,5512);wavwrite(y1,5512,work2_2.wav);分析：由 1 可知，当用 5510Hz 对音频进行采样时，不满足采样定理约束的条件，因此频谱会出现混叠，采样
6、后的声音听起来也会出现混叠的效果，音质远远不如当采样频率为 44100Hz 时。3、录制一段自己的声音，试验当 fs=？时，发生混叠？脚本文件：work3.mclear;clc;close all;y,fs=audioread(voice.m4a);Fs=16000 8000 5333 4000 3200 2666 2285 2000;de=1 2 3 4 5 6 7 8;for m=1:8y1=y(1:de(m):200000);subplot(4,2,m)ffts(y,Fs(m);title(用,num2str(Fs(m),Hz抽样voice的频谱);sound(y1,Fs(m);paus
7、e(15);endy1=y(1:de(3):200000);wavwrite(y1,Fs(3),voice.wav);分析：分别用 16000Hz、8000 Hz、5333 Hz、4000 Hz、3200 Hz、2666Hz、2285 Hz、2000 Hz8种采样频率进行采样，播放采样后的音频聆听效果，观察频谱图，得出在采样频率为 5300Hz 时恰好发生频谱混叠。采样频率低于 5333Hz时，采样频率越低，混叠越严重。4、在噪声环境中录制一段自己的声音，试采用一种方法将噪声尽可能地消除。脚本文件：work4.mclear;clc;close all;y,fs=audioread(work4.
8、m4a);b=fir1(62,2000*2*pi/fs,hann(62+1);yy=filter(b,1,y);sound(yy,fs);subplot(211)ffts(y,fs);title(录音频谱);subplot(212)ffts(yy,fs);title(消除噪音后的频谱);wavwrite(yy,fs,work4_4.wav);三、总结和分析三、总结和分析通过本次实验，进一步加深DFT算法原理和基本性质的理解，熟悉了FFT算法原理，理解和掌握了音频信号各参数的意义，通过对音频信号加入噪声而去掉噪声的过程学会了设计低通滤波器的方法，巩固了课程理论知识。理论的学习加上实验的理解，使我对奈奎斯特采样定律有了更加深刻的理解。通过查看帮助学会了如何运用FFTl对一个信号求频谱，进一步加深了我对matlab帮助的使用。在本次试验中自己编写的一个函数，使我对于写函数、用函数有了更深的理解。在实验过程中也遇到了许多问题，通过查阅资料、与同学讨论，最终解决了这些问题。本次试验也让我更加深刻的明白了Matlab的重要性，特别是对于专业上的一些问题，可以更好的让我们理解理论上的一些深奥的问题。