基于matlab的语音信号的基本处理.docVIP

..... word格式.整理版 专题研讨四、信号与系统综合应用 确定题目（根据个人兴趣、结合实际确定题目，可从下面参考题目中选择，也可自由确定）： 基于matlab的语音信号的基本处理 参考题目： 题目1：含噪信号滤波 题目2：双音多频信号的产生与检测 题目3：磁盘驱动系统仿真 题目4：卡尔曼滤波器的应用 题目5：应用反馈扩大放大器的带宽 (以上只是本专题的部分题目) 开题报告 课题题目 基于matlab的语音信号的基本处理 关键词 matlab 语音信号 基本处理 简要背景说明（课题是如何提出来的）： 学习了信号与系统及配套matlab实验,对信号的基本处理有了初步了解,尤其是实际应用的部分.我们对课件和部分matlab实验中的语音信号的处理十分感兴趣,我们想自己动手做一做基于matlab的语音信号的基本处理实践. 课题涉及的基本理论： 傅利叶变换与频域分析 抽样定理 系统响应 数字滤波 课题实施方案： 将课题分为四个部分: eq \o\ac(○,1)语音信号的制作及描述; eq \o\ac(○,2)语音信号抽取及倍插; eq \o\ac(○,3)语音信号的加噪; eq \o\ac(○,4)滤波器.分别实施 课题实施计划及阶段目标： 集中快速完成前三个部分,在此基础上钻研第 eq \o\ac(○,4)部分 本人承担任务： eq \o\ac(○,4)滤波器 ,并参与其他三个部分 参考文献： 教材及三本参考电子书 指导老师的意见和建议： 课题实施过程记录 包括仿真程序、仿真结果、结果分析、方案完善等 eq \o\ac(○,1)语音信号的制作及描述 制作语音文件: 用windows录音机录制一小段语音文件333.wav”,内容为信号与信息系统,由一同学播音. 用matlab播放”333.wav” 仿真程序: [y,Fs,bits]=wavread(333.wav); sound(y,Fs);pause; 绘画出语音文件的时域和频域波形: 仿真程序: [y,Fs,bits]=wavread(333.wav); plot(y); 仿真结果: 结果分析: 随着时间变化,声音能量图形 方案完善: 1.时间轴有问题,与实际的时间不一样. 2.语音信号的频域分析更清楚. 仿真程序: [y,Fs,bits]=wavread(333.wav); y=y(:,1); sLength=length(y); Y = fft(y,sLength); Pyy = Y.* conj(Y) / sLength; halflength=floor(sLength/2); f=Fs*(0:halflength)/sLength; figure; plot(f,Pyy(1:halflength+1));xlabel(Frequency(Hz)); t=(0:sLength-1)/Fs; figure; plot(t,y);xlabel(Time(s)); 仿真结果: 结果分析: 频域分析是从另一个角度观察信号;语音信号的一般频域范围200~2000Hz eq \o\ac(○,2)语音信号抽取及倍插 仿真程序: [y,Fs,bits]=wavread(111.wav); sound(y,Fs/2);pause; [y,Fs,bits]=wavread(111.wav); sound(y,2*Fs);pause; 仿真结果与分析: 以Fs/2及2*Fs播放的语音信号存在失真, 方案完善: 需要做出波形,做更直观的观察. eq \o\ac(○,3)语音信号的加噪 1)语音信号加高频噪音及播放. 仿真程序: [y,Fs,bits]=wavread(333.wav); y=y(:,1); sLength=length(y); t=(0:sLength-1)/Fs; f=50000;A=0.5; d=A*cos(2*pi*f*t); y=y+d; sound(y,Fs); 仿真结果: 播放时伴有尖锐的”吱吱”声. 结果分析: 由于加入高频成分余弦信号,信号叠加后出现了尖锐的噪音. 2)加噪后的语音信号的时域和频域波形. 仿真程序: [y,Fs,bits]=wavread(333.wav); y=y(:,1); sLength=length(y); t=(0:sLength-1)/Fs; f=50000;A=0.05; d=A*cos(2*pi*f*t); y1=y+d; sLength=length(y1); Y =