数字信号处理在语音信号分析中的应用.docxVIP

PAGE4

《数字信号处理》

课程设计报告

数字信号处理在语音信号分析中的应用

专业班级：

姓名：

学号：

目录

TOC\o1-3\h\z\u摘要 3

1、绪论 3

2、课程设计的具体内容 4

2.1.1、读取语音信号的任务 4

2.1.2、任务分析和解决方案 5

2.1.4、运行结果和相应的分析 5

2.2、IIR滤波器设计和滤波处理 6

2.2.1、设计任务 6

2.2.2、任务分析和解决方案 7

2.2.3、编程得到的MATLAB代码 7

2.2.4、运行结果和相应的分析 7

2.3、FIR滤波器设计和滤波处理 9

2.3.1、设计任务 9

2.3.2、任务分析和解决方案 9

2.3.3、编程得到的MATLAB代码 9

2.3.4、运行结果和相应的分析 11

3、总结 13

4、存在的不足及建议 13

5、参考文献 13

title(时域图);

xlabel(时间（n）);

ylabel(幅值（n）);

频域程序：

[z2,Fs]=audioread(F:\cjl\chenjunlin.wma);

y1=fft(z2);

y1=fftshift(y1);

sound(z2,Fs);%播放语音信号

figure(2);

derta_Fs=Fs/length(z2);

plot([-Fs/2:derta_Fs:Fs/2-derta_Fs],abs(y1));

title(频谱图);

xlabel(时间（n）);

ylabel(幅值（n）);

2.1.4、运行结果和相应的分析

由图可知人的声音分布在低频。而声音比较底和深沉。

2.2、IIR滤波器设计和滤波处理

2.2.1、设计任务

要求用自己设计的各种滤波器(高通、低通)分别对采集的信号进行滤波。

2.2.2、任务分析和解决方案

在MATLAB中，利用函数filter对信号进行滤波。利用MATLAB中的函数freqz画出各滤波器的频率响应。

2.2.3、编程得到的MATLAB代码

IIR数字低通滤波器

fs=22050;

[z3,Fs]=audioread(F:\cjl\chenjunlin.wma);

Ts=1/fs;R1=10;

wp=2*pi*1000/fs;%通带截止频率?

ws=2*pi*1200/fs;%阻带截止频率

Rp=1;%通带衰减?

Rl=100;%阻带衰减?

wp1=2/Ts*tan(wp/2);%将模拟指标转换成数字指标?

ws1=2/Ts*tan(ws/2);

[N,Wn]=buttord(wp1,ws1,Rp,R1,s);%选择滤波器的最小阶数

[Z,P,K]=buttap(N);%创建butterworth模拟低通滤波器?

[Bap,Aap]=zp2tf(Z,P,K);

[b,a]=lp2lp(Bap,Aap,Wn);%将模拟原型低通滤波器转换为低通滤波器

[bz,az]=bilinear(b,a,fs);%用双线性变换法实现模拟滤波器到数字滤波器的转换?

[H,W]=freqz(bz,az);%绘制频率响应曲线?

figure(1)

plot(W*fs/(2*pi),abs(H))

grid

xlabel(频率／Hz)

ylabel(频率响应幅度)

title(IIR低通滤波器)

f1=filter(bz,az,z3);

figure(2)

subplot(2,1,1)

plot(z3)%画出滤波前的时域图?

title(IIR低通滤波器滤波前的时域波形);

subplot(2,1,2)

plot(f1);%画出滤波后的时域图?

title(IIR低通滤波器滤波后的时域波形);

sound(f1,44100);%播放滤波后的信号?

F0=fft(f1);

f=fs*(0:511)/1024;

figure(3)

y2=fft(z3);

subplot(2,1,1);

derta_Fs=Fs/length(z3);

plot([-Fs/2:derta_Fs:Fs/2-derta_Fs],abs(y2));%画出滤波前的频谱图?

title(IIR低通滤波器滤波前的频谱)

xlabel(频率/Hz);

ylabel(幅值);

subplot(2,1,2)

derta_Fs=Fs/length(z2);

plot([-Fs/2:derta_Fs:Fs/2-derta_Fs],abs(F0));%画出滤波后的频谱图?

title(IIR低通滤波器滤波后的频谱)

xlabel(频率/Hz);

ylabel(幅值);

I