语言信号处理.

1、 齐鲁理工学院语音信号处理(论文) 课 程 名 称 语音信号处理课程设计 二 级 学 院 机电工程学院 专 业 通信工程专业 班 级 2013级通信工程一班 学 生 姓 名 刘永婷 学 号 201310537031 设 计 地 点 电气信息工程训练中心 指 导 教 师 赵 安 兴 目录 一、设计目的3二、设计要求3三、程序及其仿真3四、课程设计心得13五、参考文献131、 设计目的通过本课程设计, 主要训练和培养学生综合应用所学过的信号及信息处理等课程的相关知识，独立完成信号仿真以及信号处理的能力。包括：查阅资料、合理性的设计、分析和解决实际问题的能力，数学仿真软件Matlab和C语言程序设计

2、的学习与应用，培养规范化书写说明书的能力。2、 设计要求（1） 熟悉有关频谱分析的有关指示，对语音信号做频谱分析（2） 熟悉有关滤波器设计理论知识，选择合适的滤波器技术指标和类型，设计滤波器，得到滤波器参数（3） 实现信号频谱分析和滤波等有关Matlab函数（4） 写好总结报告，写出基本原理，有关程序，得到的图表，结果分析，总结。三、程序及其仿真 %用录音工具录一段自己的声音，转换成wav格式的文件，然后保存到MATLAB的文件夹下的work文件夹中。在MATLAB软件平台下，利用wavread函数对语音信号进行采样，记住采样频率和采样点数。1、FIR低通滤波 y,fs,nbits=wavre

3、ad(1.wav); N=length(y); noise=0.2*randn(N,1); Si=y+noise; S=fft(Si); z21=fftfilt(b21,Si); sound(z21); m21=fft(z21); subplot(2,2,1); plot(abs(S),g); title(滤波前信号频谱)； xlabel(时间); ylabel(幅度)； grid; subplot(2,2,2); plot(abs(m21),r); title(滤波后信号频谱)； xlabel(时间); ylabel(幅度)； grid; subplot(2,2,3); plot(S); t

4、itle(滤波前信号波形)； xlabel(时间); ylabel(幅度)； grid; subplot(2,2,4); plot(z21); title(滤波后信号波形)； xlabel(时间); ylabel(幅度)； grid; 2、%用录音工具录一段自己的声音，转换成wav格式的文件，然后保存到MATLAB的文件夹下的work文件夹中。在MATLAB软件平台下，利用wavread函数对语音信号进行采样，记住采样频率和采样点数。 y,fs,nbits=wavread(1.wav); sound(y,fs,nbits); N=length(y); Y=fft(y,N); subplot(2

5、,1,1); plot(y); title(原始信号波形); xlabel(时间)； ylabel(幅度)； subplot(2,1,2); plot(abs(Y); title(原始信号频谱)； xlabel(时间); ylabel(幅度)；3、构造受干扰信号并对其进行频谱分析y,fs,nbits=wavread(1.wav); N=length(y); noise=0.2*randn(N,1) Si=y+noise; sound(Si) subplot(2,1,1); plot(Si); title(加噪语音信号的时域波形); xlabel(时间); ylabel(幅度)； S=fft(S

6、i); subplot(2,1,2); plot(abs(S); title(加噪语音信号的频域波形); xlabel(时间); ylabel(幅度)； 4、 global x; global fs; global nbits; N=length(x); x1=resample(x,1,2); sound(x1,fs/2,nbits); x2=resample(x,1,4); sound(x2,fs/4,nbits); n1=0:N/2-1; n2=0:N/4-1; w1=4*n1*pi/N; w2=8*n2*pi/N; y11=fft(x1); y12=fft(x2); subplot(2,

7、2,1); plot(n1,x1); title(8KH重采样语音信号时域图)； xlabel(时间t)； ylabel(幅值)； subplot(2,2,2); plot(w1/pi,abs(y11); title(8KH重采样语音信号频域图)； xlabel(频域Hz)； ylabel(幅度)； subplot(2,2,3); plot(n2,x2); title(4KH重采样语音信号时域图)； xlabel(时间t)； ylabel(幅值)； subplot(2,2,4); plot(w2/pi,abs(y12); title(4KH重采样语音信号频域图)； xlabel(频率Hz)；

8、ylabel(幅度)；5、 global a; global fs; global b; global y; global w; global nbits; global fs; global N; n=0:N-1; y1=filter(b,a,y); y2=fft(y1); sound(y1,fs,nbits); subplot(2,1,1); plot(n,y1) title(滤波后信号频谱); xlabel(频率Hz); ylabel(幅度);四、课程设计心得这次的课程设计主要是跟随了课堂学习的内容。虽然原来已经学习了信号与系统，对于这方面有了一些基础的知识，也学了数字信号处理，但是刚开

9、始面对这个题目却也不知从何下手了。首先一个录音就为难我们半天，在老师的指导下终于迈出了这艰难的第一步，并且得到了原信号的波形图和频谱图。于是我们斗志高昂，继续进行下去。通过大致对这本书的翻阅和有目的性的查找，我们对于这次实验中要求的很多函数有了进一步的了解。由于实验结果主要要求就是图形了，因而对于图形函数的运用很多，最熟悉的就是figure函数了，然后还有画图函数plot，在IIR滤波器的设计中还使用了grid网格线命令，同时还有xlabel和ylabel以及title等命令，应该说是收获颇丰吧。而对于滤波器的设计也认识不少函数，像buttap，cheb1ap等等。很多程序里的东西，特别是程序语法，总是有错误，但是总是不知道错在哪里。在细心的检查下，终于找出了错误和警告，感到很满足。总的来说，通过这次的课程设计让我对语音信号有了全面的认识，提高了动手实践能力。同时我相信，进一步加强对MATLAB的学习和研究对我今后的学习将会起到很大的帮助。五、参考文献(1).程佩青。数字信号处