信号与系统课程设计语音滤波系统
（课题三）
课题三语音信号处理系统设计
一、本课题的目的
本设计课题主要研究语音信号抽样和恢复的软硬件实现方法、滤波器的设计及应用。通过完成本课题的设计，拟主要达到以下几个目的：
1．通过硬件实验观察连续时间信号抽样及恢复的波形特点。加深理解时域抽样定理的内容。
2．掌握利用MATLAB实现连续时间信号抽样及恢复的基本原理和方法。
3．掌握利用MATLAB分析模拟及数字系统时域、频域特性的方法；
4．了解模拟滤波器系统的设计方法、基于运算电路的模拟系统有源实现方法；通过实验平台掌握模拟系统的频率特性测试方法。
5．熟悉由模拟滤波器转换为数字滤波器的原理。
6．掌握数字滤波器的设计方法。通过设计具体的滤波器掌握滤波器设计方法、步骤。
7．了解数字滤波器的应用，了解语音信号的频率特性。
8．培养学生运用所学知识分析和解决实际问题的能力。
二、课题任务
设计一个语音信号处理系统，实现对语音信号的抽样、滤波、频谱分析以及信号的回复。要求通过硬件实验掌握其电路工作原理、测试方法以及数据处理方法，根据系统的设计技术指标通过程序设计实现系统仿真。
硬件部分：
1．利用信号与系统实验箱实现信号的抽样和恢复。
2．利用信号与系统实验箱熟悉四阶巴特沃思滤波器（或切比雪夫滤波器）的工作原理并观察记录各型滤波器的幅频特性。
软件部分：
1．根据抽样定理及语音信号频谱范围设计一个最小3阶模拟滤波器对语音信号进行预滤波，用直接、级联或并联结构实现所设计系统，对系统的时域、频域特性进行仿真测试，对结果进行分析比较。（要求保留4000Hz以内频率的信号，可采用巴特沃斯或者切比雪夫滤波器）
2.设计一个系统，要求：（1）实现连续信号的抽样。（2）针对语音信号频谱及噪声频率，设计巴特沃思数字滤波器（或切比雪夫滤波器）滤除噪声，进行频谱分析并与原始信号进行比较。（3）由滤波后信号恢复出连续信号，进行谱分析并进行回放。
3. 利用MATLAB软件的系统仿真功能(Simulink)实现系统工作过程的仿真测试，并对其结果加以分析。
注：巴特沃斯滤波器和切比雪夫滤波器任选一种。
三、主要设备和软件
1．信号与系统实验箱，一台
2. 同步信号源模块（DYT3000-57），一块
3. 信号的抽样与恢复模块（DYT3000-68），一块
4．四阶巴特沃斯滤波器模块（或四阶切比雪夫滤波器模块），一块
5．DYT3000-60模拟滤波器设计模块，两~三块
5. PC机，一台
6．基于PC的DSO-2090 USB数字示波器，一台
7．MA TLAB6.5以上版本，一套
8．数字示波器与MATLAB接口软件模块，一套
四、设计内容、步骤和要求
必做部分：
1．信号的抽样和恢复
（1）用本课题使用的硬件系统模块，实现给定信号的抽样和恢复
1）熟悉信号的抽样与恢复的工作原理。接好电源线，将信号的抽样与恢复模块和同步信号源模块插入信号系统实验平台插槽中，打开实验箱电源开关，通电检查模块灯亮，实验箱开始正常工作。
2）将同步信号源模块产生的V PP＝1V、f0＝1KHz的正弦波和f0＝2KHz的方波分别送入待抽样信号输入点S_IN和抽样脉冲信号输入点SQU_IN，用示波器分别观察抽样信号输出点PAM_OUT和恢复后的信号输出点S_OUT的波形并将实验数据记录下来（实验中低通滤波器的截止频率f C＝1kHz）。
3）改变抽样脉冲信号的频率，分别将f0＝2kHz、4kHZ、8kHz、16KHz的方波送入抽样脉冲信号输入点SQU_IN，重复实验步骤2，比较在不同的抽样频率下恢复后的信号波形之间的差别并得出结论。
2KHZ—SQU 抽样输出PAM_OUT 恢复输出S_OUT
8KHZ—SQU 抽样输出恢复输出S_OUT
16KHZ—SQU 抽样输出恢复输出S_OUT
32KHZ—SQU 抽样输出恢复输出S_OUT
4）将同步信号源模块产生的V PP＝1V、f0＝1KHz的三角波作为待抽样信号送入S_IN，重复上述实验步骤。
注：使用2k正弦波作被抽样信号时效果较好，可以自行比较。
三角
2KHZ—SQU 抽样输出恢复输出
8KHZ—SQU 抽样输出恢复输出
16KHZ—SQU 抽样输出恢复输出
32KHZ—SQU 抽样输出恢复输出
分析：
（2）用软件实现指导给定信号的抽样和恢复（给定的信号频率及抽样频率太小）
对连续信号)5.0cos()(t t x π=以抽样频率f s =2Hz 抽样得到)()(nT x n x =，编程完成下列任务：画出信号)(t x ，40≤≤t 及其在相应范围内的抽样序列；利用抽样内插函数
)1
(),()(s
r f T T t Sa t h ==π恢复连续时间信号。画出信号x(t)和重建信号)(t x r 的波形，比较
这两个信号。若信号x(t)与)(t x r 存在较大差异，应如何改善？
% function y=sinc(x) y=ones(size(x)); i=find(x);
y(i)=sin(pi*x(i))./(pi*x(i));
t=0:0.001:4;
x=cos(0.5*pi*t); subplot(3,1,1); plot(t,x); axis([0,4,-1.1,1.1]); title('原信号');
fs=2; %采样频率 ts=1/fs;
tp=4; %重构时间区间为4s n=0:tp/ts;
xn=cos(n*ts*0.5*pi); subplot(3,1,2) stem(n*ts,xn);
title('fs=2Hz 抽样信号');
dt=ts/50; %每个采样间隔上取点数 t1=0:dt:4; %生成序列t1
tp=4; %重构时间区间为4s n=0:tp/ts;
tmn=ones(length(n),1)*t1-n'*ts*ones(1,length(t1)); xr1=sinc(fs*tmn); x2=xn*xr1; subplot(3,1,3)
plot(t1,x2);title('用抽样内插函数恢复的结果');
figure(2)
subplot(3,1,1);
plot(abs(fft(x)));title('原信号的频谱'); subplot(3,1,2);
plot(abs(fft(xn)));title('抽样信号的频谱'); subplot(3,1,3);
plot(abs(fft(x2)));title('恢复后信号频谱');
（3）用软件对给定的数字语音信号（采样频率为16000Hz）进行恢复（可以利用interp1函数），并重新进行采样，重新采样的频率分别为8000Hz和4000Hz，对恢复的语音信号及重采样的语音信号进行回放（利用sound函数），比较语音的变化并记录处理过程中所得各种波形及频谱图。
[y,fs]=wavread('bu14.wav')
figure(1)
subplot(2,2,1);plot(y);title('原语音信号波形')
y0=abs(fft(y))
subplot(2,2,2);plot(y0);title('原语音信号的频谱')
x=0:length(y)-1;
x1= 0:2.5:length(y)-1;
y1 = interp1(x,y,x1);
subplot(2,2,3);plot(x,y,x1,y1);title('恢复后的波形');
z=fft(y1);
y2=abs(z);
subplot(2,2,4);plot(y2);title('恢复后波形的频谱');
sound(y1);
figure(2)
y3= decimate(y,2);
subplot(2,2,1);plot(y3);title('8k采样波形');
y4=abs(fft(y3))
subplot(2,2,2);plot(y4);title('8k采样波形频谱');
x2=0:length(y3)-1;
x3= 0:2.5:length(y3)-1;
y5= interp1(x2,y3,x3);
subplot(2,2,3);plot(x2,y3,x3,y5);title('恢复后的波形');
z=fft(y5);
y6=abs(z);
subplot(2,2,4);plot(y6);title('恢复后波形的频谱');
figure(3)
y7=decimate(y,4);
subplot(2,2,1);plot(y7);title('4k采样波形');
y8=abs(fft(y7));
subplot(2,2,2);plot(y8);title('4k采样波形频谱');
x4=0:length(y7)-1;
x5= 0:2.5:length(y7)-1;
y9= interp1(x4,y7,x5);
subplot(2,2,3);plot(x4,y7,x5,y9);title('恢复后的波形'); z=fft(y9);
y10=abs(z);
subplot(2,2,4);plot(y10);title('恢复后波形的频谱');
2. 滤波器设计
（1）四阶巴特沃思滤波器（或四阶切比雪夫滤波器）参数测量
（a）四阶巴特沃思滤波器参数测量
四阶巴特沃思滤波器使用信号源单元和四阶巴特沃斯滤波器模块
?熟悉四阶巴特沃斯滤波器设计、工作原理。接好电源线，将四阶巴特沃斯滤波器模块插入信号系统实验平台插槽中，打开实验箱电源开关，通电检查模块灯亮，实验箱开始正常工作。
?扫频信号产生
1）将信号源单元的开关k2向下拨，切换至扫频输出。
2）用示波器观察OUT1测试点波形，为一频率可调的锯齿波信号，该信号为扫频压控信号。调节电位计OUT1 Freq使锯齿波信号为50Hz，并将“波形选择”跳线的第一组引脚连接。
3）用示波器观察OUT2正弦波扫频输出信号，电位计GAIN Adj则用于调节输出正弦波的幅度。
4）调节电位计OUT1 Freq，改变送入锯齿波压控信号的频率，重复上述步骤，用示波器观察扫频输出信号的变化。
5）分别连接“频率选择”跳线的四组引脚，改变扫频信号的扫频段，重复上述步骤，用示波器观察不同扫频段的扫频输出信号。
?输入输出参考点说明
OUT1 扫频压控锯齿波信号输出点
OUT2 正弦波扫频信号输出点。
?扫频源法观察滤波器的特性：
1) 将信号源单元“波形选择”跳线的第1组引脚连接，并将开关k2向下拨，切换至扫频输出，按照前述步骤得到扫频正弦波信号，并用示波器观察OUT1点锯齿波频率，将其调为50Hz，作为扫频压控信号。
2) 将OUT2输出的扫频信号送入四阶巴特沃斯低通滤波器信号输入点BLP_IN，用示波器观察输出点BLP_OUT的输出信号波形。
3) 将锯齿波压控信号和低通滤波器输出信号分别接示波器的X轴和Y轴，观察李沙育图形。
4) 将扫频信号分别送入四阶巴特沃斯高通和带通滤波器，重复上述实验步骤，分别观察各种滤波器的输出信号波形。
?描点法观察滤波器的幅频特性曲线：
1）将信号源单元产生的固定频率正弦波送入低通滤波器的信号输入端BLP_IN，用示波器观察BLP_OUT的输出波形，测量波形的电平值（有效值），记录此时的电平值及频率。
2）调节电位计OUT2 Freq，改变输入正弦波信号的频率（保持信号幅度不变），重复步骤1。
3）整理实验数据，以频率为X轴，以幅度（电平）为Y轴，绘出幅频特性图。
4）将频率正弦波信号分别送入四阶巴特沃斯高通和带通滤波器，重复上述实验步骤，绘出各种滤波器的幅频特性曲线。
?研究各滤波器对方波信号或其它非正弦波信号输入的响应（实验步骤自拟）。
输入、输出点参考说明
BLP_IN、BLP_OUT：四阶巴特沃斯低通滤波器信号输入点、信号输出点。
BHP_IN、BHP_OUT：四阶巴特沃斯高通滤波器信号输入点、信号输出点。
BBP_IN、BBP_OUT：四阶巴特沃斯带通滤波器信号输入点、信号输出点。
高通
（2）模拟及数字滤波器设计
1）根据以上测量的指标，设计模拟巴特沃思（或切比雪夫）低通、高通、带通滤波器，画出幅频特性（模拟滤波器幅频特性freqs）。
四阶巴特沃思滤波器参数测量
模拟巴特沃斯低通滤波器
F(kHz） 5.15 10.72 15.31 20.49 25.63 29.73 32.31 33.42 35.27 37.52
Vp-p（v）0.76 0.72 0.65 0.54 0.45 0.28 0.24 0.22 0.19 0.16
1 0.94 0.85 0.71 0.59 0.37 0.3
2 0.28 0.25 0.21
X=[5.15 10.72 15.31 20.49 25.63 29.73 32.31 33.42 35.27 37.52 ]
Y=[0.76 0.72 0.65 0.54 0.45 0.28 0.24 0.22 0.19 0.16 ]
plot(X,Y)
grid on
模拟巴特沃斯高通滤波器
F(kHz) 17.62 18.70 20.53 23.88 26.43 28.34 30.32 32.41 34.29 35.82 38.72 40.33 Vp-p(v) 0.20 0.24 0.28 0.36 0.46 0.54 0.57 0.60 0.62 0.63 0.63 0.64
0.31 0.37 0.44 0.56 0.72 0.84 0.89 0.94 0.96 0.98 0.98 1
X=[17.62 18.70 20.53 23.88 26.43 28.34 30.32 32.41 34.29 35.82 38.72
40.33 ]
Y=[0.20 0.24 0.28 0.36 0.46 0.54 0.57 0.60 0.62 0.63 0.63 0.64] plot(X,Y)
grid on
模拟巴特沃斯带通滤波器
F(kHz) 16.51 19.96 20.23 21.53 22.72 23.56 24.43 25.62 27.19 28.75 30.35 Vp-p(v) 0.12 0.42 0.58 0.66 0.69 0.71 0.72 0.73 0.58 0.50 0.35
0.16 0.58 0.79 0.90 0.95 0.97 0.99 1 0.80 0.69 0.47 32.71 33.21 35.66
0.26 0.22 0.13
0.36 0.30 0.18
X=[16.51 19.96 20.23 21.53 22.72 23.56 24.43 25.62 27.19 28.75 30.35 32.71 33.21 35.66 ]
Y=[0.12 0.42 0.58 0.66 0.69 0.71 0.72 0.73 0.58 0.50 0.35 0.26 0.22 0.13 ]
plot(X,Y)
grid on
（2）模拟及数字滤波器设计
1）根据以上测量的指标，设计模拟巴特沃思低通、高通、带通滤波器，画出幅频特性（模拟滤波器幅频特性freqs）。
模拟巴特沃思低通滤波器
N=4;wc=21700*2*pi
[B,A]=butter(N,wc,'s') % B，A分别为基本二阶节的分子和分母的系数
[H,w]=freqs(B,A);
plot(w/2/pi/1000,abs(H));
xlabel('频率(khz)');ylabel('幅度');title('低通');
wc =1.3635e+005
B =1.0e+020 *
0 0 0 0 3.4559
A = 1.0e+020 *
0.0000 0.0000 0.0000 0.0001 3.4559
模拟巴特沃思高通滤波器
wc=25000*2*pi;
[B,A]=butter(N,wc,'high','s')
[H,w]=freqs(B,A);
plot(w/2/pi/1000,abs(H));
xlabel('频率(khz)');ylabel('幅度');title('高通');
B =
1 0 0 0 0
A =
1.0e+020 *
0.0000 0.0000 0.0000 0.0001 6.0881
模拟巴特沃思带通滤波器
wc=[20.7,28.5]*1000*2*pi
[B,A]=butter(N,wc,'bandpass','s')
[H,w]=freqs(B,A);
plot(w/2/pi/1000,abs(H));
xlabel('频率');ylabel('幅度');title('带通');
wc = 1.0e+005 *
1.3006 1.7907
B = 1.0e+018 *
Columns 1 through 6
0 0 0 0 5.7690 0 Columns 7 through 9
0 0 0
A = 1.0e+041 *
Columns 1 through 6
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Columns 7 through 9
0.0000 0.0000 2.9424
2）将模拟各型滤波器转换为数字滤波器，画出幅频特性（求数字滤波器幅频特性freqz，采样频率fs=200kHz）。
3）利用MATLAB软件实现设计，记录相关波形，加以分析修改。
function [db,mag,pha,w]=freqz_m(b,a);
[H,w]=freqz(b,a,1000,'whole'); %在0-2*pi之间选取N个点计算频率响应
H=(H(1:501))'; %频率响应
w=(w(1:501))'; %频率
mag=abs(H); %响应幅度
db=20*log10((mag+eps)/max(mag)); %增益
pha=angle(H); %相位
function [b,a]=imp_invr(c,d,T)
[R,p,k]=residue(c,d); %部分分式展开
p=exp(p*T); %从模拟到数字极点对应关系 ,部分分式系数相同
[b,a]=residuez(R,p,k); %将部分分式的形式变换成多项式之比的形式
b=real(b'); %求出数字滤波器系数
a=real(a');
数字巴特沃斯低通滤波器
N=4;
wc=21700*2*pi;
fs=2*10^5;
课程设计报告 课程名称信号与系统课程设计指导教师 设计起止日期 学院信息与通信工程 专业电子信息工程 学生 班级/学号 成绩 指导老师签字
目录 1、课程设计目的 (1) 2、课程设计要求 (1) 3、课程设计任务 (1) 4、课程设计容 (1) 5、总结 (11) 参考文献 (12) 附录 (12)
1、课程设计目的 “信号与系统”是一门重要的专业基础课，MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。本课程设计基于MATLAB完成信号与系统综合设计实验，以提高学生的综合应用知识能力为目标，是“信号与系统”课程在实践教学环节上的必要补充。通过课设综合设计实验，激发学生理论课程学习兴趣，提高分析问题和解决问题的能力。 2、课程设计要求 （1）运用MATLAB编程得到简单信号、简单信号运算、复杂信号的频域响应图； （2）通过对线性时不变系统的输入、输出信号的时域和频域的分析，了解线性时不变系统的特性，同时加深对信号频谱的理解。 3、课程设计任务 （1）根据设计题目的要求，熟悉相关容的理论基础，理清程序设计的措施和步骤； （2）根据设计题目的要求，提出各目标的实施思路、方法和步骤； （3）根据相关步骤完成MATLAB程序设计，所编程序应能完整实现设计题目的要求； （4）调试程序，分析相关理论； （5）编写设计报告。 4、课程设计容 （一）基本部分 （1）信号的时频分析 任意给定单频周期信号的振幅、频率和初相，要求准确计算出其幅度谱，并准确画出时域和频域波形，正确显示时间和频率。 设计思路： 首先给出横坐标，即时间，根据设定的信号的振幅、频率和初相，写出时域波形的表达式；然后对时域波形信号进行傅里叶变化，得到频域波形；最后使用plot函数绘制各个响应图。 源程序： clc; clear; close all; Fs =128; % 采样频率 T = 1/Fs; % 采样周期 N = 600; % 采样点数 t = (0:N-1)*T; % 时间，单位：S x=2*cos(5*2*pi*t);
******************* 实践教学 ******************* 兰州理工大学 计算机与通信学院 2013年春季学期 信号处理课程设计 题目：基于MATLAB的FIR滤波器语音信号去噪专业班级：通信工程（1）班 姓名：王兴栋 学号：10250114 指导教师：陈海燕 成绩：
摘要 语音信号在数字信号处理中占有极其重要的地位，因此选择通过对语音信号的研究来巩固和掌握数字信号处理的基本能力十分具有代表性。对数字信号处理离不开滤波器，因此滤波器的设计在信号处理中占有极其重要的地位。而MATLAB 软件工具箱提供了对各种数字滤波器的设计。本论文“在MATLAB平台上实现对语音信号的去噪研究与仿真”综合运用了数字信号处理的各种基本知识，进而对不带噪语音信号进行谱分析以及带噪语音信号进行谱分析和滤波处理。通过理论推导得出相应的结论，再通过利用MATLAB作为编程工具来进行计算机实现比价已验证推导出来的结论。在设计过程中，通过设计FIR数字滤波器和IIR数字滤波器来完成滤波处理。在设计过程中，运用了MATLAB对整个设计中的图形的绘制和一些数据的计算以及仿真。 关键字滤波器；MATLAB；仿真；滤波
前言 语音是语言的声学表现，是人类交流信息最自然、最有效、最方便的手段。随着社会文化的进步和科学技术的发展，人类开始进入了信息化时代，用现代手段研究语音处理技术，使人们能更加有效地产生、传输、存储、和获取语音信息，这对于促进社会的发展具有十分重要的意义，因此，语音信号处理正越来越受到人们的关注和广泛的研究。 语音信号是信息技术处理中最重要的一门科学，是人类社会几步的标志。那么什么是语音？语音是人类特有的功能，也是人类获取外界信息的重要工具，也是人与人交流必不可少的重要手段。那么什么又是信号？那信号是什么呢？信号是传递信息的函数。离散时间信号——序列——可以用图形来表示。 语音信号处理是一门用研究数字信号处理研究信号的科学。它是一新兴的信息科学，同时又是综合多个学科领域的一门交叉科学。语音在我们的日常生活中随时可见，也随处可见，语音很大程度上可以影响我们的生活。所以研究语音信号无论是在科学领域上还是日常生活中都有其广泛而重要的意义。 本论文主要介绍的是的语音信号的简单处理。本论文针对以上问题，运用数字信号学基本原理实现语音信号的处理，在matlab7.0环境下综合运用信号提取，幅频变换以及傅里叶变换、滤波等技术来进行语音信号处理。我所做的工作就是在matlab7.0软件上编写一个处理语音信号的程序，能对语音信号进行采集，并对其进行各种处理，达到简单语音信号处理的目的。 对语音信号的研究，本论文采用了设计两种滤波器的基本研究方法来达到研究语音信号去噪的目的，最终结合图像以及对语音信号的回放，通过对比，得出结论。
基于MATLAB的语音信号滤波处理 题目：基于MATLAB的语音信号滤波处理 课程：数字信号处理 学院：电气工程学院 班级： 学生： 指导教师： 二O一三年十二月
目录CONTENTS 摘要 一、引言 二、正文 1.设计要求 2.设计步骤 3.设计内容 4.简易GUI设计 三、结论 四、收获与心得 五、附录
一、引言 随着Matlab仿真技术的推广，我们可以在计算机上对声音信号进行处理，甚至是模拟。通过计算机作图，采样，我们可以更加直观的了解语音信号的性质，通过matlab编程，调用相关的函数，我们可以非常方便的对信号进行运算和处理。 二、正文 2.1 设计要求 在有噪音的环境中录制语音，并设计滤波器去除噪声。 2.2 设计步骤 1.分析原始信号，画出原始信号频谱图及时频图，确定滤波器类型及相关指标； 2.按照类型及指标要求设计出滤波器，画出滤波器幅度和相位响应，分析该滤波器是否符合要求； 3.用所设计的滤波器对原始信号进行滤波处理，画出滤波后信号的频谱图及时频图； 4.对滤波前的信号进行分析比对，评估所设计滤波器性能。 2.3 设计内容 1.原始信号分析
分析信号的谱图可知，噪音在1650HZ和3300HZ附近的能量较高，而人声的能量基本位于1000HZ以下。因此，可以设计低通滤波器对信号进行去噪处理。 2.IIR滤波器设计 用双线性变换法分别设计了巴特沃斯低通滤波器和椭圆低通滤波器和带阻滤波器： ①巴特沃斯滤波器 fp=800;fs=1300;rs=35;rp=0.5; 程序代码如下： fp=800;fs=1300;rs=35;rp=0.5;Fs=44100; wp=2*Fs*tan(2*pi*fp/(2*Fs));ws=2*Fs*tan(2*pi*fs/(2*Fs)); [n,wn]=buttord(wp,ws,rp,rs,'s'); [b,a]=butter(n,wn,'s'); [num,den]=bilinear(b,a,Fs); [h,w]=freqz(num,den,512,Fs);
数字信号处理课程设计 题目：基于matlab的语音信号滤波处理学院：物理与电子信息工程 专业：电子信息工程 班级： B07073041 学号： 200932000066 姓名：高珊 指导教师：任先平
摘要： 语音信号处理是研究用数字信号处理技术和语音学知识对语音信号进行处理的新兴学科，是目前发展最为迅速的学科之一，通过语音传递信息是人类最重要，最有效，最常用和最方便的交换信息手段，所以对其的研究更显得尤为重要。 Matlab语言是一种数据分析和处理功能十分强大的计算机应用软件，它可以将声音文件变换成离散的数据文件，然后用起强大的矩阵运算能力处理数据。这为我们的本次设计提供了强大并良好的环境！ 本设计要求自己录制一段自己的语音后，在MATLAB软件中采集语音信号、回放语音信号并画出语音信号的时域波形和频谱图。再在Matlab中分别设计不同形式的FIR数字滤波器。之后对采集的语音信号经过不同的滤波器（低通、高通、带通）后，观察不同的波形，并进行时域和频谱的分析。对比处理前后的时域图和频谱图，分析各种滤波器对于语音信号的影响。最后分别收听进行滤波后的语音信号效果，做到了解在怎么样的情况下该用怎么样的滤波器。
目录 1.设计内容 (4) 2.设计原理 (4) 2.1语音信号的时域分析 (4) 2.2语音信号的频域分析 (5) 3.设计过程 (5) 3.1实验程序源代码 (6) 3.1.1原语音信号时域、频域图 (6) 3.1.2低通滤波器的设计 (6) 3.1.3高通滤波器的设计 (7) 3.1.4带通滤波器的设计 (8) 3.1.5语音信号的回放 (9) 3.2调试结果描述 (10) 3.3所遇问题及结果分析 (15) 3.3.1所遇主要问题 (16) 3.3.2结果分析 (16) 4.体会与收获 (17) 5.参考文献 (17)
信号与系统课程设计报告 实验题目：信号的运算与处理 内容简介： 设计一个信号，对其进行信号运算和处理，利用Matlab仿真。 课设方式： 利用电子技术、电路理论和信号与系统的知识学习验证信号的运算和处理，如延时、相加、微分、抽样等。自已设计信号及运算方式，并利用Matlab仿真。 分析计算结果。 课程设计要求： 独立完成； 完成信号设计（任意信号均可）及其某种运算（任意运算均可，也可多做几种，或做组合运算）的验证； 学会利用Matlab仿真；提交课程设计报告。 例如： 设计一个信号为f(t)=3sin2t 对其做微分运算得到f/(t) , 用MATLAB 编程实现计算过程，画出f(t)和f/(t)
本次课程设计本人选的信号运算是： 设计一个信号为y1=y(x)=sin2x,对其作微分运算得到dy1,用MATLAB对其实现运算过程，后画出y1,dy1,y1+dy1的图像 实验步骤（操作过程） 1、 首先打开MATLAB软件，在其命令窗口直接输入以下程序，对y(x)进 行微分运算。得到dy1 clear >> syms x y1; >> y1=sin(2*x); >> dy1=diff(y1,'x') dy1 =2*cos(2*x) 运算过程如下图所示： 2、 接着便是对其进行验证，点击fire，新建一个文件，输入以下程序（绘制出y1=sin2x, dy1=2cos2x, 以及y1+ dy1=sin2x+2cos2x。的波形）
3、保存文件，后缀名为.m,随后按F5执行输出输出图形。实验结果如下图所示 、
结果分析 如图所示绿色波形为y1=sin2x，蓝色为dy1=2cos2x,红色波形为y1+dy1。仿真结果与运算结果一致。 实验心得体会（调试过程） 总的来说，这次课程设计难度并不是太高，而我选取的正玄信号也是较为简单常用的一种函数，对其进行微分运算之后，得到了余弦函数，其仿真结果波形也如上所示，与预期一致。在设计过程中，还是出现了几个小问题的，一个是变量的定义，之前没有定义x，直接取范围结果出错了，还有一个是注意各种函数的调用以及运算格式，还是希望能在之后再接再厉，掌握好matlab软件！（附上调试过程图片） 左边为文件、历史窗口，底下是命令窗口，最右下角为实验仿真波形，中间为运算程序，绘图画图程序。
数字信号处理实验报告 实验名称：基于MATLAB对语音信号 进行分析及滤波处理 院系：物联网工程学院 班级：电子信息工程1101 姓名：
一、实验目的 综合计算运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应的结论，培养发现问题、分析问题和解决问题的能力。并利用MATLAB作为工具进行实现，从而复习巩固课堂所学的理论知识，提高对所学知识的综合应用能力，并从实践上初步实现对数字信号的处理。此外，还系统的学习和实现对语音信号处理的整体过程，从语音信号的采集到分析、处理、频谱分析、显示和储存。 二、实验要求 1.分析原始语音信号的时域特性和频谱特性。 2.设计一个IIR滤波器，用该滤波器对语音信号进行滤波处理，分析滤波后信号的时域特性和频谱特性。 3.设计一个FIR滤波器，用该滤波器对语音信号进行滤波处理，分析滤波后信号的时域特性和频谱特性。 三、实验原理 1.采样定理：在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中最高频率fmax的2倍时(fs.max>=2fmax)，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。 采样频率越高，即采样的时间间隔越短，对声音波形的表示越精确。 2.时域信号的FFT分析：信号的频谱分析就是计算机信号的傅里叶变换。连续信号与系统的傅里叶分析显然不便于用计算机进行计算，使其应用受到限制。而FFT是一种时域和频域均离散化的变换，适合数值运算，成为用计算机分析离散信号和系统的有力工具。对连续信号和系统，可以通过时域采样，应用DFT进行近似谱分析。 3.IIR数字滤波器设计原理利用双线性变换设计IIR滤波器（巴特沃斯数字低通滤波器的设计），首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得所要设计的IIR滤波器的系统函数H(z)。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里主要是边界频率Wp和Ws的转换，对ap和as指标不作变化。边界频率的转换关系为∩=2/T tan(w/2)。接着，按照模拟低通滤波器的技术指标根据相应设计公式
课程设计任务书 学生姓名：苗强强专业班级：电信1204 指导教师：阙大顺沈维聪工作单位：信息工程学院 题目: 程控宽带放大器的设计 初始条件： 程控宽带放大器是电子电路中常用模块，在智能仪器设备及嵌入式系统中有广 泛的应用。因此对于电子信息专业的技术人员来说，熟练掌握该项技术很有必要。 要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体 要求） （1）输入阻抗>1KΩ,单端输入，单端输出，放大器负载电阻为600Ω； （2）3dB通频带10kHz~6MHz，在20kHz~5MHz频带内增益起伏<1dB。 （3）增益调节范围10 dB~40 dB，（通过键盘操作调节）。 （4）发挥部分：当输入频率或输出负载发生变化时，通过微处理器自动调节，保持 放大器增益不变。 （5）电路通过仿真即可。 时间安排： 1. 任务书下达，查阅资料 1天 2. 制图规范、设计说明书讲解 2天 3. 设计计算说明书的书写 5天 4. 绘制图纸 1天 5. 答辩 1天 指导教师签名：年月日 系主任（或责任教师）签名：年月日
滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。 关键词数字滤波器 MATLAB IIR滤波器 FIR滤波器
目录 1 绪论 (2) 2课程设计的具体实现 (4) 2.1 语音信号的录制 (4) 2.2 FIR滤波器的设计方法 (4) 3语音信号的时频分析 (5) 3.1语音信号载入MATLAB (5) 3.2语音信号时域和频域分析 (6) 3.3设计FIR数字滤波器 (7) 3.3.1 窗函数hamming设计带通滤波器 (8) 3.3.2 窗函数hanning设计带通滤波器 (9) 3.3.3 窗函数Blackman设计带通滤波器 (10) 3.3.4 窗函数Boxcar设计带通滤波器 (11) 3.4滤波之后的时域和频域分析 (12) 3.4.1 窗函数hamming设计带通滤波器滤波 (13) 3.4.2 窗函数hanning设计带通滤波器滤波 (13) 3.4.3 窗函数Blackman设计带通滤波器滤波 (14) 3.4.4 窗函数Boxcar设计带通滤波器滤波 (15) 3.5回放语音信号 (16)
结论 (16) 参考文献 (18) 1 绪论 数字信号处理是利用计算机或专用处理设备，以数值计算的方法对信号进行采集、抽样、变换、综合、估值与识别等加工处理，借以达到提取信息和便于应用的目的。它在语音、雷达、图像、系统控制、通信、航空航天、生物医学等众多领域都获得了极其广泛的应用。具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等优点。 数字滤波器, 是数字信号处理中及其重要的一部分。随着信息时代和数字技术的发展，受到人们越来越多的重视。数字滤波器可以通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。数字滤波器种类很多，根据其实现的网络结构或者其冲激响应函数的时域特性，可分为两种，即有限冲激响应( FIR，Finite Impulse Response)滤波器
一、课程设计题目 去除干扰蜂鸣音 1.目的：掌握信号时频域分析方法，正确理解采样定理，准确理解滤波器的概念。 2.内容：提供一个包含某人说话语音片段的声音文件，但该语音信号被一个包含有几个谐波分量的蜂鸣信号干扰了。 用Matlab 的wavread 命令读取该声音文件。注意，该命令可以同时得到声音文件的采样率和采样位宽，请查阅Matlab 的帮助文件。 （1） 用快速傅立叶变换（FFT ）计算并画出声音信号的频谱，列写出蜂鸣信号的谐波频率。 （2） 思考如何将这些蜂鸣音去除？将去除了蜂鸣音的语音片段播放出来，仔细聆听并写下语音片 段中人物所说的话。注意：由于只能播放实信号，因此记得提取信号的实部。 Matlab 命令：wavread, wavplay, fft, fftshift, fir1, filter, plot, figure. 二、设计思路 用waveread()函数读取音频和其采样率和采样位宽，对读取的音频信号使用fft()函数进行快速傅立叶变换并绘出得到的频谱。观察频谱分析噪声（蜂鸣信号）的谐波频率分布，选择合适的滤波模式将噪声信号的谐波滤去，便可以得到去除噪声后的人声。 设计滤波器的频域特性便成了除去噪声并留下原声的关键，我们注意到所学的采样定理以及一维sinc 函数（辛格函数）x x x Sinc ππ) sin()(=，然而汉宁窗可以看作是3个矩形时间窗的频谱之和，或者 说是 3个)(x Sinc 型函数之和，而括号中的两项相对于第一个谱窗向左、右各移动了π/T ，从而使旁瓣互相抵消，消去高频干扰和漏能。它适用于非周期性的连续信号。利用它的选择特性使用fir1()建立一个汉宁窗函数，并用filter()函数进行滤波，去除噪声部分。最后用play()函数播放音频检查效果。 三、设计过程 1.音频的读取和分析 先将原始音频文件读入， [audio0, Fs, nbits] = wavread('C:\Users\Administrator\Desktop\signal\buzz.wav')；%按路径读取音频存入audio0变量，并用Fs 变量记录采样率，nbits 变量记录采样位宽。 其中， Fs=11025 #采样率为11025Hz nbits=32 #采样带为32 p0=audioplayer(audio0,Fs);%将audio0载入音频播放器 play(p0)；%并进行播放 subplot(2,1,1);%分屏绘图 plot(audio0);%绘制原始音频时域图,如下图所示 title('时域');%标注题目
通信与信息工程学院 信息处理综合实验报告 班级：电子信息工程1502班 指导教师： 设计时间：2018/10/22-2018/11/23 评语： 通信与信息工程学院 二〇一八年 实验题目：语音信号分析与处理 一、实验内容 1. 设计内容 利用MATLAB对采集的原始语音信号及加入人为干扰后的信号进行频谱分析，使用窗函数法设计滤波器滤除噪声、并恢复信号。 2．设计任务与要求 1. 基本部分
（1）录制语音信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （2）对所录制的语音信号加入干扰噪声，并对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （3）分别利用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman 窗几种函数设计数字滤波器滤除噪声，并画出各种函数所设计的滤波器的频率响应。 （4）画出使用几种滤波器滤波后信号时域波形和频谱，对滤波前后的信号、几种滤波器滤波后的信号进行对比，分析信号处理前后及使用不同滤波器的变化；回放语音信号。 2. 提高部分 （5）录制一段音乐信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （6）利用MATLAB产生一个不同于以上频段的信号；画出信号频谱图。 （7）将上述两段信号叠加，并加入干扰噪声，尝试多次逐渐加大噪声功率，对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （8）选用一种合适的窗函数设计数字滤波器，画出滤波后音乐信号时域波形和频谱，对滤波前后的信号进行对比，回放音乐信号。 二、实验原理 1.设计原理分析 本设计主要是对语音信号的时频进行分析，并对语音信号加噪后设计滤波器对其进行滤波处理，对语音信号加噪声前后的频谱进行比较分析，对合成语音信号滤波前后进行频谱的分析比较。 首先用PC机WINDOWS下的录音机录制一段语音信号，并保存入MATLAB软件的根目录下，再运行MATLAB仿真软件把录制好的语音信号用audioread函数加载入MATLAB仿真软件的工作环境中，输入命令对语音信号进行时域，频谱变换。 对该段合成的语音信号，分别用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman窗几种函数在MATLAB中设计滤波器对其进行滤波处理，滤波后用命令可以绘制出其频谱图，回放语音信号。对原始语音信号、合成的语音信号和经过滤波器处理的语音信号进行频谱的比较分析。 2.语音信号的时域频域分析 在Matlab软件平台下可以利用函数audioread对语音信号进行采样,得到了声音数据变量y,同时把y的采样频率Fs=44100Hz放进了MATALB的工作空间。
DSP课程设计报告语音信号的FIR滤波器处理 姓名：张伟槟 班级：12电信3 学号：201231190425 指导老师：徐梅宣 日期：2015.06.02～2014.06.19 华南农业大学电子工程学院
摘要 随着信息与数字技术的发展，数字信号处理已经成为当今极其重要而学科与技术领域之一。它在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理的基本方法中，通常会涉及到变换、滤波、频谱分析、调制解调和编码解码等处理。其中滤波是应用非常广泛的一个环节，数字滤波器的理论和相关设计也一直都是人们研究的重点之一。FIR滤波器的是非递归的，稳定性好，精度高；更重要的是，FIR滤波器在满足幅频响应要求的同时，可以获得严格的线性相位特征。因此，它在高保真的信号处理，如数字音频、图像处理、数据传输和生物医学等领域得到广泛应用。 CSS集成开发环境使用CCS内置的软件仿真simulator对程序进行编译，调试和运行，主要用于检测目标程序运行的正确性和连贯性,并能通过仿真器与目标板连接，在目标板上实时观察效果。 在本次设计中，我们选择的课题是基于DSP的语言信号的FIR滤波处理。首先利用MATLAB进行了仿真，得到滤波前后的时域波形和频谱。然后通过调用MATLAB的分析工具FDATOOL，根据仿真结果导出了滤波器的相关参数，将原始信号数据和滤波器参数输入CCS进行DSP编程。最后在DSP中实现了FIR低通滤波，并通过CCS的频谱分析功能查看了最终DSP的滤波效果。 关键词:语音信号处理 FIR滤波 MATLAB CCS
目录 1.前言 (4) 2.设计需求 (4) 3.滤波器设计原理介绍 (4) 3.1数字滤波器设计原理 (4) 3.2FIR滤波器的基本结构 (5) 3.3 滤波器的特点 (5) 3.4窗函数的介绍 (6) 4.FIR滤波器的设计 (7) 4.1FIR滤波器设计 (7) 4.2窗函数设计的基本方法 (7) 4.3滤波器的1 z算法实现 (7) 5.FIR滤波器的MATLAB设计 (8) 6.FIR滤波器的DSP实现 (10) 6.1CCS设计与仿真 (10) 6.2CCS仿真结果 (13) 6.调试问题及解决办法 (16) 7.设计感想 (16) 参考文献 (17) 8.附录 (17)
南通大学电子信息学院信号与系统课程设计报告 班级： 姓名： 学号： 指导老师： 设计时间： 2014—2015学年第一学期
一、连续信号的时域分析 二、 1. 信号的产生 (1)阶跃函数 function [t,y1]=jieyue(t1,t2,t0) dt=0.01; ttt=t1:dt:t0-dt; tt=t0:dt:t2; t=t1:dt:t2; n=length(ttt); nn=length(tt); u=zeros(1,n); uu=ones(1,nn); y1=[u,uu]; return 冲激函数 function [t,y2]=chongji(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); y2(1:n)=0; y2(1,(t0-t1)/dt+1)=1/dt; （2）调用上述函数产生信号）2-t ε（，）（4-t δ，-t e ）（t ε，-6s ≤t ≤6s,并画出波形。 Command Window subplot(3,1,1); [t1,y1]=jieyue(-6,6,2);
stairs(t1,y1); axis([-6 6 0 1.5]); subplot(3,1,2); [t2,y2]=chongji(-6,6,4);plot(t2,y2); subplot(3,1,3); [t3,y3]=jieyue(-6,6,0); y3=exp(-(t3)).*y3;plot(t3,y3); 波形如下图所示： （3）根据f(t)画出f(2t)和f(1-0.5t)的波形 t=-3:0.01:3; y=tripuls(t,4,0.6); subplot(3,1,1); plot(t,y);
《计算机仿真技术》 基于MATLAB的加噪语音信号的滤波学生姓名： 专业：电子信息工程 班级： 学号： 指导教师： 完成时间：2017年12月
一．滤波器的简述 在MATLAB环境下IIR数字滤波器和FIR数字滤波器的设计方法即实现方法，并进行图形用户界面设计，以显示所介绍迷你滤波器的设计特性。 在无线脉冲响应（IIR）数字滤波器设计中，先进行模拟滤波器的设计，然后进行模拟-数字滤波器转换，即采用脉冲响应不变法及双线性Z变化法设计数字滤波器，最后进行滤波器的频带转换。在有限脉冲响应（FIR）数字滤波器设计中，讨论了FIR线性相位滤波的特点和用窗口函数设计FIR数字滤波器两个问题。两类滤波器整个过程都是按照理论分析、编程设计、集体实现的步骤进行的。为方便分析直观者直观、形象、方便的分析滤波器的特性，创新的设计出图形用户界面---滤波器分析系统。整个系统分为两个界面，其内容主要包括四个部分：System（系统）、Analysis（分析）、Tool（工具）、Help（帮助）。 数字滤波在DSP中占有重要地位。数字滤波器按实现的网络结构或者从单位脉冲响应，分为IIR（无限脉冲响应）和FIR（有限脉冲响应）滤波器。如果IRR 滤波器和FIR滤波器具有相同的性能，那么通常IIR滤波器可以用较低的阶数获得高的选择性，执行速度更快，所有的存储单元更少，所以既经济又高效。二．设计要求 1.在matlab平台上录制一段语音信号； 2.完成语音信号的谱分析； 3.对语音信号进行加噪以及加噪后信号的谱分析； 4.选择合适的滤波器进行滤波，确定相关指标； 5.实现滤波过程，显示滤波后的结果，并进行谱分析。 三．实验内容与步骤 1、语音信号的录入
DSP技术及应用课程设计 语音滤波处理 院系：机电工程学院 专业（班级）：电子信息工程2班 姓名：洪育钦 学号: 20134082033 指导教师：韩杨杨 职称：助教 完成日期： 2016年 6月 16日 成绩：
目录 1.引言 (1) 1.2 设计的目的 (1) 1.3设计指标要求 (2) 2滤波器的基础知识 (2) 2.1 滤波器的定义 (2) 2.2 滤波器的功能 (3) 2.3滤波器的特点 (3) 2.4 FIR数字滤波器的窗函数设计 (4) 3.FIR滤波器设计与实现过程 (5) 3.1 FIR滤波器设计指标 (5) 3.2 FIR滤波器设计过程 (6) 3.2.2 FIR滤波器的matlab实现 (6) 3.2.3 FIR滤波器的CCS实现 (8) 3.3 FIR滤波器设计结果 (9) 4.设计过程中遇到的问题及解决的办法 (11) 参考文献 (12) 附录 (12)
语音信号的FIR滤波器处理 1.引言 随着信息时代和数字世界的到来，数字信号处理已成为如今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。 数字信号处理器，也称DSP芯片，是针对数字信号处理需要而设计的一种具有特殊结构的微处理器，它是现代电子技术、相结合的产物。一门主流技术，随着信息处理技术的飞速发展，计算机技术和数字信号处理技术逐渐发展成为使它在电子信息、通信、软件无线电、自动控制、仪表技术、信息家电等高科技领域得到了越来越广泛的应用。 数字滤波是语音处理、图像处理、频谱分析等应用中的基本处理算法。DSP 是一种处理数字信号的专用微处理器, 主要应用于实时快速地实现各种信号的数字处理算法。用 DSP 芯片实现数字滤波具有稳定性好、精确度高、不受环境影响等优点。 数字滤波器分为有限冲激响应滤波器（FIR滤波器）和无限冲激响应滤波器（IIR滤波器）。FIR 滤波器属于经典滤波器，优点就是由于不存在系统极点，FIR 滤波器是绝对稳定的系统，FIR 滤波器还确保了线性相位，在信号处理中占有极其重要的地位。数字滤波器一直以来就是数字信号处理器（DSP）最广为人知的应用，FIR 滤波器的单位冲激响应 b(n)为有限长序列,若 b(n)为实数, 且满足偶对称:b(n)= b(N- 1- n)的条件, 称为系数对称FIR 滤波器。系数对称 FIR 滤波器在数字信号处理中应用十分广泛。 1.2 设计的目的
一、 1.正弦信号 A = input('input A=') ;% 给正弦信号的幅度A赋值 w = input('input w=') ; % 给正弦信号的频率w赋值 theta =input('input theta='); % 给正弦信号的初始相位theta 赋值disp(['这个信号是周期信号']) T=2*pi/w t = 0 : 0.01 : 3*T ; % 定义时间点 ft = A * sin( w * t + theta ) ; % th计算函数值 plot( t ,ft ) ; % 画图 title( '正弦信号' ) ; % 为图像加标题注释 grid on ; % 在图上画方格
2.复指数信号 j00 = sqrt( - 1 ) ; % 定义复数j a = input('input a='); % 复指数信号赋值w = input('input w='); K = input('input K='); if a==0 disp('这是一个周期信号') T=2*pi/w else if a>0 disp('这不是一个周期信号') else disp('这不是一个周期信号') end end t = -1.5*abs(a) : 0.01 : 1.5*abs(a) % 定义时间点 ft = K*exp( ( a + j00 * w ) * t ) ; subplot( 2 , 2 , 1 ) ; plot( t , real( ft ) ) ; title( '实部' ) ; %画图subplot( 2 , 2 , 2 ) ; plot( t , imag( ft ) ) ; title( '虚部' ) ; subplot( 2 , 2 , 3 ) ; plot( t , abs( ft ) ) ; title( '模' ) ; subplot( 2 , 2 , 4 ) ; plot( t , angle( ft ) ) ; title( '相角' ) ;
基于数字滤波器的语音信号处理系统的设计 前言：语音信号的数字滤波就是利用快速傅里叶变换fft对语音信号进行频谱分析，利用分析得到的指标设计数字滤波器，将不需要的部分通过设计好的数字滤波器滤除掉，以达到优化语音的目的。数字滤波器是一种对数字信号进行处理的系统，数字滤波器完成数字信号滤波处理功能，其输入是一组数字量，其输出是经过变换的另一组数字量。它通过一定的运算关系或电路形式来改变输入信号所包含的频率成分的相对比例或滤出某些频率成分，数字滤波器具有稳定性高、精度高、灵活性大等优点。随着数字技术的发展，用数字技术实现滤波器的功能越来越受到人们的注意和广泛的应用。 摘要：本文介绍了两种数字滤波器的基本设计原理及步骤1、利用双线性变换法以及脉冲响应不变法设计IIR数字滤波器2、用窗函数法实现FIR数字滤波器的设计。并且把理论与实践结合，通过一段语音，利用MATLAB软件分析语音信号频谱，确定指标设计滤波器，编写程序设计滤波器，对语音信号进行滤波处理，分析滤波后的时域和频域特征，回放语音信号分析其经过信号处理后的效果。 正文 1、总体设计 1.1 总体设计内容 1）利用脉冲响应不变法和双线性变换法设计IIR数字滤波器，实现数字低通、数字高通、数字带通、数字带阻。 2）利用窗函数实现FIR数字滤波器的设计 3）设计完整的语音信号处理系统，可实现语音信号的读入、语音信号的频域分析、滤波器选择、语音信号的时/频域滤波操作、语音信号的短时频域分析等功能。 1.2 课设流程图设计
2、模块设计 2.1 语音信号的采集 用windows工具中的录音机录一段（语音信号），将语音信号的文件命名为test.wav. 将语音文件保存后，在MATLAB软件平台下，首先调用wavread函数可采集到录制的音乐信号，并得到其采样率fs和比特数bits。具体调用如下：>> [y,fs,bits]=wavread('d:\test.wav')%读取语音信号。 运行后得出fs，bit。其中参数y是每个样本的值，fs是生成该波形文件时的采样率，bits 是波形文件每样的编码位数。 sound(y); 用于对声音的回放。向量y则就代表了一个信号，也即一个复杂的“函数表达式”，也可以说像处理一个信号的表达式一样处理这个声音信号。 2.2 语音信号的时域及频域分析 语音信号分析是语音信号处理的前提和基础，只有分析出可表示语音信号本质特征的参数，才有可能利用这些参数进行高效的语音通信、语音合成和语音识别等处理。而且，语音合成的音质好坏，语音识别率的高低，也都取决于对语音信号分桥的准确性和精确性。因此语音信号分析在语音信号处理应用中具有举足轻重的地位。 2.2.1语音信号的时域分析 语音信号的时域分析就是分析和提取语音信号的时域参数。进行语音分析时，最先接触到并且也是最直观的是它的时域波形。语音信号本身就是时域信号，因而时域分析是最早使用，也是应用最广泛的一种分析方法，这种方法直接利用语音信号的时域波形。 运用matlab软件，调用plot函数，在读取语音信号的基础上画出语音信号的时域波形图。理想状态下，语音信号的时域波形应该是关于时间轴对称的，但是实际情况下，由于存在着一些噪音，使得波形看上去不是那么完美。 2.2.2语音信号的频域分析 对语音信号进行频谱分析。MATLAB提供了快速傅里叶变换算法FFT计算DFT的函数fft,调用格式如下： Xk=fft(xn,N) 参数xn为被变换的时域序列向量，N是DFT变换区间的长度，当N大于xn时，fft函数自动在xn后面补零，当N小于xn的长度时，fft函数计算xn的前N个元素，忽略其后面的元素。在本次设计中，调用fft函数得到语音信号的快速傅里叶变换序列，再用plot函数画出信号的频谱图。 2.3 IIR数字滤波器的设计 IIR滤波器系统函数的极点可以再单位圆内的任何位置，实现IIR滤波器的阶次较低，所用的存储单元较少，效率高，又由于IIR数字滤波器能够保留一些模拟滤波器的优良特性，因 此应用很广。设计IIR数字滤波器的方法主要有基于脉冲响应不变法和双线性Z变换法设计，
DSP实验课程设计实验报告 姓名：学号：班级： 1.课程设计题目： 基于MATLAB的有噪声的语音信号处理的课程设计。 2.课程设计的目的： 综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应的结论，再利用MATLAB做为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。 3.课程设计的要求： (1)熟悉离散信号和系统的时域特性。 (2)掌握序列快速傅里叶变换FFT方法。 (3)学会MATLAB的使用，掌握MATLAB的程序设计方法。 (4)利用MATLAB对语音信号进行频谱分析。 (5)掌握MATLAB设计各种数字滤波器的方法和对信号进行滤波的方法。 4.课程设计的内容： 录制一段语音信号，对语音信号进行频谱分析，利用MATLAB中的随机函数产生噪声加入到语音信号中，使语音信号被污染，然后进行频谱分析，设计FIR和IIR数字滤波器，并对噪声污染的语音信号进行滤波，分析滤波后的信号的时域和频域特征，回放语音信号。 5.课程设计的步骤： （1）语音信号的获取 通过录音软件录制一段语音“数字信号处理”，命名为“OriSound”，时长大约1到2秒，在MATLAB中，通过使用wavread函数，对语音进行采样: [y,fs,nbits]=wavread('OriSound'); %语音信号的采集 采样值放在向量y中，采样频率为fs，采样位数为nbits。
（2）语音信号的频谱分析 画出语音信号的时域波形，然后对语音信号进行频谱分析，在MATLAB中，通过使用fft 函数对信号进行快速傅里叶变换，得到信号的频谱特性。 因此采集语音并绘出波形和频谱的模块程序如下： [y,fs,nbits]=wavread('OriSound'); %语音信号的采集 sound(y,fs,nbits); %语音信号的播放 n=length(y) ; %计算语音信号的长度 Y=fft(y,n); %快速傅里叶变换 figure; subplot(2,1,1); %绘出时域波形 plot(y); title('原始信号波形','fontweight','bold'); axis([ 00000 80000 -1 1]); %通过尝试确定合适的坐标参数 grid; subplot(2,1,2); %绘出频域频谱 plot(abs(Y)); title('原始信号频谱','fontweight','bold'); axis([ 0 150000 0 4000]); %通过尝试确定合适的坐标参数 grid; 结果如下：
信号与系统课程设计 课程名称：信号与系统 题目名称：滤波器的设计与实现 学院：电气与电子工程学院 专业班级：电气工程及其自动化 学号：3 学生：宗喜 指导教师：黄劲 2015年12 月20 日
目录 一、设计要求 (2) 二、设计原理 (2) 三、设计思路 (3) 四、设计容 (3) A、一阶有源滤波电路 (3) B、二阶有源滤波电路 (5) 1、二阶低通滤波电路 (5) 2、二阶高通滤波电路 (6) 3、二阶带通滤波电路 (8) C、用仿真软件设计滤波器 (10) 1、给定性能参数设计滤波器 (10) a、二阶低通滤波器 (10) b、二阶高通滤波器 (11) c、二阶带通滤波器 (12) 2、不同阶数滤波器性能比较 (12) D、滤波器的Matlab设计仿真 (13) 1、二阶低通滤波器 (13) 2、二阶高通滤波器 (14) 五、参考文献 (16)
一、设计要求 自已设计电路系统，构成低通滤波器、高通滤波器和带通滤波器。利用Matlab或其他仿真软件进行仿真。 有源滤波器由是有源元件和无源元件(一般是R和C)共同组成的电滤波器。和无源滤波器相比，它的设计和调整过程较简便，此外还能提供增益。因此，本课程设计中选择了二阶有源滤波器作为主要研究对象。 1、自行设计电路图，确定前置放大电路，有源滤波电路，功率放大电路的方案， 并使用绘图软件（Electronics Worrkbench）画出设计电路，包括低通、高通和带通。 2、所设计的滤波器不仅有滤波功能，而且能起放大作用，负载能力要强。 3、根据给定要求和电路原理图计算和选取单元电路的元件参数。 4、用Matlab或其他仿真软件（FilterLab）对滤波器进行仿真，记录仿真结果。 二、设计原理 1、电容器C具有通高频阻低频的性能。 2、由源滤波器由放大电路部分和滤波电路部分组成。 3、仿真软件可以将滤波器的性能直观的表现出来。 4、各种滤波器的幅频特性：
广西工学院 数字信号处理课程设计设计题目：语音信号的处理与滤波 系别： 学号： 姓名： 班级： 指导教师： 完成日期：
目录 1. 摘要 (3) 2.MATLAB简介 (3) 3.设计目的 (4) 4.设计内容 (4) 5.设计原理 (4) 6.设计步骤 (5) 7.总结与分析 (10) 8.参考资料 (10)
摘要 本课题分析了数字信号处理课程的重要性及特点，可以帮助理解与掌握课程中的基本概念、基本原理、基本分析方法；并利用MATLAB对语音信号进行分析和处理，要求采集语音信号后，在MATLAB软件平台进行频谱分析。用设计的数字滤波器对语音信号进行了滤波 MATLAB简介 MATLAB 是一种对技术计算高性能的语言。它集成了计算，可视化和编程于一个易用的环境中，在此环境下，问题和解答都表达为我们熟悉的数学符号。典型的应用有： ?数学和计算 ?算法开发 ?建模，模拟和原形化 ?数据分析，探索和可视化 ?科学与工程制图 ?应用开发，包括图形用户界面的建立 MATLAB是一个交互式的系统，其基本数据元素是无须定义维数的数组。这让你能解决很多技术计算的问题，尤其是那些要用到矩阵和向量表达式的问题。而要花的时间则只是用一种标量非交互语言(例如C或Fortran)写一个程序的时间的一小部分。 . 名称“MATLAB”代表matrix laboratory（矩阵实验室）。MATLAB最初是编写来提供给对由LINPACK和EINPACK工程开发的矩阵软件简易访问的。今天，MATLA B使用由LAPACK和ARPACK工程开发的软件，这些工程共同表现了矩阵计算的软件中的技术发展。 MATLAB已经与许多用户输入一同发展了多年。在大学环境中，它是很多数学类、工程和科学类的初等和高等课程的标准指导工具。在工业上，MATLAB是高产研究、开发和分析所选择的工具。 MATLAB以一系列称为工具箱的应用指定解答为特征。对多数用户十分重要的是，工具箱使你能学习和应用专门的技术。工具箱是是MATLAB函数（M-文件）的全面的综合，这些文件把MATLAB的环境扩展到解决特殊类型问题上。具有可用工具箱的领域有：信号处理，控制系统神经网络，模糊逻辑，小波分析，模拟等等。