数字信号处理大作业
更新时间:2024-01-29 16:52:01 阅读量: 教育文库 文档下载
《数字信号处理》课题设计
基于Matlab的音频采样
姓名:谌海龙 学号:20134361 指导教师:仲元红 班级:13电子4班
成绩:
重庆大学通信工程学院
2015年11月
基于MATLAB的音乐采样实验
一、实验内容及原理
内容:
1、用fs=44,100Hz采集一段音乐;
2、改变采样率,用fs=5,512Hz采集一段音乐,体会混叠现象; 3、录制一段自己的声音,试验当fs=?时,发生混叠?
4、在噪声环境中录制一段自己的声音,试采用一种方法将噪声尽可能地消除。
原理:
根据奈奎斯特采样定律,如果连续时间信号xa?t?是最高截止频率为?m的带限信
?a?t?通过一个增益为T,截止频率为号,采样频率?s?2?m,那么让采样信号x?s/2的理想低通滤波器,可以无失真地恢复出原连续时间信号xa?t?。否则,?s?2?m会造成采样信号中频谱混叠现象,不能无失真地恢复连续时间信号。
采样过程
连续信号xa?t??p?t?采样信号?a?t?x采样脉冲
xa?t??2?sXa(j?)1??m0P(j?)?m??s??s0?s2?s?(j?)Xa1/T0???sm?Xa(j?)1/T?s??2?s??s0?s2?s3?s?采样信号的频谱混叠
1Xa?j??op(t)t??mo?m?P?j??
采样信号的频谱
人发出的声音频率在300~3400Hz,可以通过加低通滤波器消除高频噪声。
二、实验过程
脚本文件ffts.m的功能是画y的频谱图 function ffts(y,fs)
NFFT1=2^nextpow2(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、用fs=44,100Hz采集一段音乐;
脚本文件:work1.m clear;clc;close 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');
分析:根据采样定理,当采样频率fs<2fm时,会产生频率混叠现象,听起来的
声音会有种很多时刻的声音重叠的效果。由于人耳所能听到的声音在20Hz到20000Hz之间,因此当采样频率为44100Hz时刚好满足采样定理所要求的条件fs>2fm,所以这也验证了抽样定理中对于采样频率的要求fs≥2fm,其中fm为信号的频谱中最高频率分量
2、改变采样率,用fs=5,512Hz采集一段音乐,体会混叠现象;
脚本文件:work2.m clear;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对音频进行采样时,不满足采样定理约束的条件,因此频谱会出现混叠,采样后的声音听起来也会出现混叠的效果,音质远远不如当采样频率为44100Hz时。
3、录制一段自己的声音,试验当fs=?时,发生混叠?
脚本文件:work3.m clear;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:8
y1=y(1:de(m):200000); subplot(4,2,m) ffts(y,Fs(m));
title(['用',num2str(Fs(m)),'Hz抽样voice的频谱']); sound(y1,Fs(m));
pause(15); end
y1=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.m clear;clc;close all;
[y,fs]=audioread('work4.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的重要性,特别是对于专业上的一些问题,可以更好的让我们理解理论上的一些深奥的问题。
subplot(211) ffts(y,fs); title('录音频谱'); subplot(212) ffts(yy,fs);
title('消除噪音后的频谱'); wavwrite(yy,fs,'work4_4.wav');
三、总结和分析
通过本次实验,进一步加深DFT算法原理和基本性质的理解,熟悉了FFT算法原理,理解和掌握了音频信号各参数的意义,通过对音频信号加入噪声而去掉噪声的过程学会了设计低通滤波器的方法,巩固了课程理论知识。理论的学习加上实验的理解,使我对奈奎斯特采样定律有了更加深刻的理解。通过查看帮助学会了如何运用FFTl对一个信号求频谱,进一步加深了我对matlab帮助的使用。在本次试验中自己编写的一个函数,使我对于写函数、用函数有了更深的理解。在实验过程中也遇到了许多问题,通过查阅资料、与同学讨论,最终解决了这些问题。本次试验也让我更加深刻的明白了Matlab的重要性,特别是对于专业上的一些问题,可以更好的让我们理解理论上的一些深奥的问题。
正在阅读:
数字信号处理大作业01-29
煤矿安全知识培训教案107-17
自然感觉02-14
PHC预应力管桩监理实施细则06-03
07年6月-09年6月英语六级真题听力原文 - Word版学习啊06-27
2017年湘西自治州初中毕业学业考试数学试题卷(无答案)(基础题04-07
小学一年级数学上册期末总复习应用题练习卷12-27
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 信号处理
- 作业
- 数字
- 危险性较大工程施工安全管理制度
- 西北工业大学深入学习实践科学发展观活动整改实施方案
- 酒店服务语言技巧案例培训教案
- 实验五 接口与包
- EDA技术实验书
- 拉西瓦水电站水垫塘和二道坝混凝土施工技术要求(2007.3)
- 编译原理课后习题答案
- 《仪表与自动化》解题参考
- 中国财政收入和经济增长关系的实证分析
- 软件项目管理大作业
- 教师学习传统文化心得体会
- 2018年社区关工委工作总结
- 江南大学2018年上半年设计概论第1阶段练习题题目
- 泄露天机--33位06年高考状元学习方法揭密
- 工程造价实习报告
- 华为ptn题目汇总
- 第四章 JSP的内置对象习题
- 关于我国当前社会公德建设的现状的调查报告
- 卫生部等4部门关于商业保险机构参与新型农村合作医疗经办服务的指导意见
- 新视野大学英语四Quiz1-ReadingUnit1-2题目及答案