对语音信号进行分析及处理

一、设计目的

1.进一步巩固数字信号处理的基本概念、理论、分析方法和实现方法；使自身对信号的采集、处理、传输、显示和存储等有一个系统的掌握和理解；

2.增强应用Matlab语言编写数字信号处理的应用程序及分析、解决实际问题的能力；

3．培养自我学习的能力和对相关课程的兴趣；

二、设计过程

1、语音信号的采集

采样频率，也称为采样速度或者采样率，定义了每秒从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。

采样位数可以理解为声卡处理声音的解析度。这个数值越大，解析度就越高，录制和回放的声音就越真实

采样定理又称奈奎斯特定理，在进行模拟/数字信号的转换过程中，当采样频率fs不小于信号中最高频率fm的2倍时，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍。

利用Windows下的录音机，录制了一段发出的声音，内容是“数字信号”，时间在3 s内。接着在D盘保存为WAV格式，然后在Matlab软件平台下．利用函数wavread对语音信号进行采样，并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。

[x1,fs,bits]=wavread('E:\数字信号.wav'); %读取语音信号的数据，赋给变量x1，返回频率fs 44100Hz，比特率为16 。

2 、语音信号的频谱分析

（1）首先画出语音信号的时域波形；

程序段：

x=x1(60001:1:120000); %截取原始信号60000个采样点

1

plot(x) %做截取原始信号的时域图形

title('原始语音采样后时域信号');

xlabel('时间轴 n');

ylabel('幅值 A');

（2）然后用函数fft对语音号进行快速傅里叶变换，得到信号的频谱特性；

y1=fft(x,6000); %对信号做N=6000点FFT变换

figure(2)

subplot(2,1,1),plot(k,abs(y1));

title('|X(k)|');

ylabel('幅度谱');

subplot(2,1,2),plot(k,angle(y1));

title('arg|X(k)|');

ylabel('相位谱');

（3）产生高斯白噪声，并且对噪声进行一定的衰减，然后把噪声加到信号中，再次对信号进行频谱特性分析，从而加深对频谱特性的理解；

d=randn(1,60000); %产生高斯白噪声

d=d/100; %对噪声进行衰减

x2=x+d; %加入高斯白噪声

3、设计数字滤波器

(1)IIR低通滤波器性能指标通带截止频，阻带截止频率Hz?1000f c?Hz1200?f。，阻带最小衰减，通带最大衰减dB100??dB?1st21(2)FIR低通滤波器性能指标通带截

止频率，阻带截止频率Hz?1000f c????Hz?1200f100dB，阻带衰减，通带衰减。

1dB12st（3）IIR高通滤波器的设计指标，，，阻带最小衰Hz?2000f Hz?1000f zp A?30dB，通带最大衰减减。dB1A?sP f?2000Hz，阻带，4（4）（）FIR高通滤波器的设计指标，Hz1000f?zp A?50dB，通带最大衰减。最小衰减dB?A1sP（5）用自己设计的各滤波器分别对采集的信号进行滤波，在Matlab中，FIR滤波器利用函数fftfilt 对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。比较滤波前后语音信号的波形及频谱，在一个窗口同时画出滤波前后2

的波形及频谱。

在Matlab中，函数sound可以对声音进行回放。其调用格式：sound(x，fs，bits)；可以感觉滤波前后的声音有变化。

三、结果分析

1、原始语音信号采样后的时域波形及FFT变换后频谱

图 1.原始信号时域图

图 2.原始信号的频谱图

3

2、加入噪声后信号的时域波形及FFT变换后频谱

加入噪声后的时域信号图3.分析：由图可以看出加入噪声后有明显的不一样，运行sound(x2 ,fs);播放加入高斯噪声后的语音，可以明显听出兹兹的噪声

加入噪声后的频谱图4.图

4

3IIR后的时域波形频谱变化滤波器及原始信号通过、IIR

5. IIR图低通滤波器

低通滤波器滤波前后时域波形图6. IIR分析：经过滤波器后的信号和原始信号很近似 5

图7. IIR低通滤波器滤波前后的频谱

分析：从图7可以看出，经过IIR低通滤波器滤波后将高频部分滤除了。

4、FIR滤波器及原始信号通过FIR后的时域波形频谱变化

《测试信号分析与处理》实验报告

测控1005班齐伟0121004931725 （18号）实验一差分方程、卷积、z变换 一、实验目的 通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。 二、实验设备 1、微型计算机1台； 2、matlab软件1套 三、实验原理 Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。 差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。用x表示滤波器的输入，用y表示滤波器的输出。 a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1) ak,bk 为权系数，称为滤波器系数。 N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。 y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2) 等式定义了数字卷积，*是卷积运算符。输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。 传输函数H(z)是滤波器的第三种实现方法。 H(z)=输出/输入= Y(z)/X(z) （3）即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。 序列x[n]的z变换定义为 X (z)=∑x[n]z-n (4) 把序列x[n] 的z 变换记为Z{x[n]} = X(z)。

应用Matlab对含噪声语音信号进行频谱分析及滤波

应用Matlab对含噪声的语音信号进行频谱分析及滤波 一、实验内容 录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；在语音信号中增加正弦噪声信号(自己设置几个频率的正弦信号)，对加入噪声信号后的语音信号进行频谱分析；给定滤波器的性能指标，采用窗函数法和双线性变换设计数字滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比试听，分析信号的变化。 二、实现步骤 1．语音信号的采集 利用Windows下的录音机，录制一段自己的话音，时间在1 s内。然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，（可用默认的采样频率或者自己设定采样频率）。 2．语音信号的频谱分析 要求首先画出语音信号的时域波形；然后对语音号进行快速傅里叶变换，得到信号的频谱特性。 在采集得到的语音信号中加入正弦噪声信号，然后对加入噪声信号后的语音号进行快速傅里叶变换，得到信号的频谱特性。并利用sound试听前后语音信号的不同。

分别设计IIR和FIR滤波器，对加入噪声信号的语音信号进行去噪，画出并分析去噪后的语音信号的频谱，并进行前后试听对比。 3.数字滤波器设计 给出数字低通滤波器性能指标:如，通带截止频率fp＝10000 Hz，阻带截止频率fs＝12000 Hz（可根据自己所加入噪声信号的频率进行阻带截止频率设置），阻带最小衰减Rs＝50 dB，通带最大衰减Rp＝3 dB（也可自己设置），采样频率根据自己语音信号采样频率设定。

报告内容 一、实验原理 含噪声语音信号通过低通滤波器，高频的噪声信号会被过滤掉，得到清晰的无噪声语音信号。 二、实验内容 录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；在语音信号中增加正弦噪声信号(自己设置几个频率的正弦信号)，对加入噪声信号后的语音信号进行频谱分析；给定滤波器的性能指标，采用窗函数法和双线性变换设计数字滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比试听，分析信号的变化。给出数字低通滤波器性能指标:如，通带截止频率fp＝10000 Hz，阻带截止频率fs＝12000 Hz （可根据自己所加入噪声信号的频率进行阻带截止频率设置），阻带最小衰减Rs＝50 dB，通带最大衰减Rp＝3 dB（也可自己设置），采样频率根据自己语音信号采样频率设定。 三、实验程序 1、原始信号采集和分析 clc;clear;close all; fs=10000; %语音信号采样频率为10000 x1=wavread('C:\Users\acer\Desktop\voice.wav'); %读取语音信号的数据，赋给x1 sound(x1,40000); %播放语音信号 y1=fft(x1,10240); %对信号做1024点FFT变换 f=fs*(0:1999)/1024; figure(1); plot(x1) %做原始语音信号的时域图形 title('原始语音信号'); xlabel('time n'); ylabel('fuzhi n'); figure(2); plot(f,abs(y1(1:2000))); %做原始语音信号的频谱图形 title('原始语音信号频谱') xlabel('Hz'); ylabel('fuzhi');

语音信号分析与处理2011

数字信号处理实验二：语音信号分析与处理 学号 姓名 注：1）此次实验作为《数字信号处理》课程实验成绩的重要依据，请同学们认真、独立完成，不得抄袭。 2）请在授课教师规定的时间内完成； 3）完成作业后，请以word 格式保存，文件名为：学号+姓名 4）请通读全文，依据第2及第3 两部分内容，认真填写第4部分所需的实验数据，并给出程序内容。 1. 实验目的 (1) 学会MATLAB 的使用，掌握MATLAB 的程序设计方法 (2) 掌握在windows 环境下语音信号采集的方法 (3) 掌握MATLAB 设计FIR 和IIR 滤波器的方法及应用 (4) 学会用MATLAB 对语音信号的分析与处理方法 2. 实验内容 录制一段自己的语音信号，对录制的语音信号进行采样，画出采样后语音信号的时域波形和频谱图，确定语音信号的频带范围；使用MATLAB 产生白噪声信号模拟语音信号在处理过程中的加性噪声并与语音信号进行叠加，画出受污染语音信号的时域波形和频谱图；采用双线性法设计出IIR 滤波器和窗函数法设计出FIR 滤波器，画出滤波器的频响特性图；用自己设计的这两种滤波器分别对受污染的语音信号进行滤波，画出滤波后语音信号的时域波形和频谱图；对滤波前后的语音信号进行时域波形和频谱图的对比，分析信号的变化；回放语音信号，感觉与原始语音的不同。 3. 实验步骤 1）语音信号的采集与回放 利用windows 下的录音机或其他软件录制一段自己的语音（规定：语音内容为自己的名字，以wav 格式保存，如wql.wav ），时间控制在2秒之内，利用MATLAB 提供的函数wavread 对语音信号进行采样，提供sound 函数对语音信号进行回放。 [y,fs,nbits]=wavread(file), 采样值放在向量y 中，fs 表示采样频率nbits 表示采样位数。Wavread 的更多用法请使用help 命令自行查询。 2）语音信号的频谱分析 利用fft 函数对信号进行频谱分析 3）受白噪声干扰的语音信号的产生与频谱分析 ①白噪声的产生： N1=sqrt （方差值）×randn(语音数据长度，2)（其中2表示2列，是由于双声道的原因） 然后根据语音信号的频谱范围让白噪声信号通过一个带通滤波器得到一个带限的白噪声信号 N2； 带通滤波器的冲激响应为： h B （n ）= ))((sin ))((sin 1122απ ωπωαπωπω---n c n c c c c c

广东工业大学《测试技术与信号分析》测试实验报告

测试技术与信号处理实验报告 机械转子底座的振动测量和分析 一、实验目的 1．掌握磁电式速度传感器的工作原理、特点和应用。

2．掌握振动的测量和数据分析。 二、实验内容和要求 先利用光电式转速传感器测量出电机的转速；然后利用磁电式速度传感器测量机械转子底座在该电机转速下的振动速度；对测量出的振动速度信号进行频谱分析；找出振动信号的主频与电机转速之间的关系。 三、实验步骤 1．启动实验程序“机械转子系统的振动测量.exe”; 输入个人信息，也可以启动之后通过单击“修改”按钮修改个人信息。 2．单击“采样设置”按钮，输入采集卡连接磁电速度传感器的采样通道号，批量采样频率（建议设为10KHz）、批量采样点数（建议设为10000）。 3．打开转子电机的电源，单击“单点采样”。 4．旋转调节旋钮改变转子的转速,观察图形区显示的磁电速度传感器采集到的转子底座振动信号；如果振动信号比较小，可适当提高转子的转速。 5．转子转速的测量： (1) 单击“采样设置”按钮，输入采集卡连接光电转速传感器的 采样通道号、批量采样频率（建议值为10KHz）、批量采样点 数（建议值为10000）。 (2) 单击“批量采样”按钮，开始采样；采样完成之后，采集到 的波形信号会显示在图形窗口，系统会自动计算出转子的速度

并显示出来。记录下此时的转子的转速（单位：r/s）。 (3) 再重复步骤(2)测量2次。以三次测量的平均值作为此时转子 的转速。 转速的测量结果 单点采样采集通道6,测量3组数据 6．振动信号的测量和频谱分析： (1) 单击“采样设置”按钮，输入采集卡连接磁电速度传感器的 采样通道号、批量采样频率（建议设为10KHz）、批量采样点 数（建议设为10000）。 (2) 单击“批量采样”按钮，开始采样；采样完成之后，采集到 的波形信号会显示在图形窗口。如果信号不正常，重复点击“批 量采样”按钮 (3) 单击“保存”按钮，将采集到的磁电传感器的信号数据保存 为文本文件。文件必须保存到“C:\ExperiData\”目录下。可单 击“保存设置”更改文件名。 (4) 打开刚保存的文本文件，文件前面几行保存了个人信息、采 样频率、采样通道、保存的数据个数等信息。文件中共有四列 数据，第一列为数据的序号，第二列为磁电传感器检测到的数 据。

测试技术复习资料第七章测试信号的处理与分析考试重点

测试技术复习资料 第七章 测试信号的处理与分析 考试重点 一、选择题 1. 两个正弦信号间存在下列关系：（ B ） A. 同频相关，不同频也相关 B. 同频相关，不同频不相关 C. 同频不相关，不同频相关 D. 同频不相关，不同频也不相关 2. 自相关函数是一个（ B ）函数。 A. 奇 B. 偶 C. 非奇非偶 D. 三角 3. 如果一信号的自相关函数)(τx R 呈现一定周期的不衰减，则说明该信号（ B ）。 A. 均值不为0 B. 含有周期分量 C. 是各态历经的 D. 不含有周期分量 4. 正弦信号的自相关函数是（ A ），余弦函数的自相关函数是（C ）。 A. 同频余弦信号 B. 脉冲信号 C. 偶函数 D. 正弦信号 5.经测得某信号的相关函数为一余弦曲线，则其（ C ）是正弦信号的（ D ）。 A. 可能 B. 不可能 C. 必定 D. 自相关函数 6. 对连续信号进行采样时，采样频率越高，当保持信号的记录的时间不变时，则（ C ）。 A. 泄漏误差就越大 B. 量化误差就越小 C. 采样点数就越多 D. 频域上的分辨率就越低 7. 把连续时间信号进行离散化时产生混叠的主要原因是（ B ）。 A. 记录时间太长 B. 采样间隔太宽 C. 记录时间太短 D. 采样间隔太窄 8. 若有用信号的强度、信噪比越大，则噪声的强度（C ）。 A. 不变 B. 越大 C. 越小 D. 不确定 9. A/D 转换器是将（ B ）信号转换成（ D ）信号的装置。 A. 随机信号 B. 模拟信号 C. 周期信号 D. 数字信号 10. 两个同频方波的互相关函数曲线是（ C ）。 A. 余弦波 B. 方波 C. 三角波 D. 正弦波 11. 已知x （t ）和y （t ）为两个周期信号，T 为共同的周期，其互相关函数的表达式为（ C ）。 A. dt t y t x T T )()(210?+τ B. dt t y t x T T )()(210?+τ C. dt t y t x T T )()(10?+τ D. dt t y t x T T )()(210?-τ 12. 两个不同频率的简谐信号，其互相关函数为（ C ）。 A. 周期信号 B. 常数 C. 零 D. 非周期信号 13. 数字信号处理中，采样频率s f 与限带信号最高频率h f 间的关系应为（ B ）。 A. s h f f = B. 2s h f f > C. s h f f < D. h s f f 7.0= 14. 正弦信号)sin()(0?ω+=t x t x 的自相关函数为（ B ）。 A. ωτsin 2 x B. ωτcos 220x C . ωτsin 220x D. ωτcos 20x

基于MATLAB的语音信号的分析与处理

基于MATLAB的语音信号分析与处理 [摘要]语音信号的处理是一门非常重要的学科，如今普遍应用在电话通信、助听器等方面。本次设计主要是为了在MATLAB软件的帮助下处理一段加噪的声音信号，该过程会涉及到采样定理，傅立叶变换等理论和算法在设计过程中的实际应用。在本次设计中，我们关注的是在驱除噪声污染是所需要的滤波器的选择，充分比较各种优缺点后，再利用滤波器来驱除杂音。通过滤波前后的声音的频谱图的比较，来了解滤波器的特性和作用，并得到本次设计的结果。 [关键词]:语音信号；MATLAB；傅立叶变换；滤波器

目录 摘要 (Ⅰ) 1 绪论 (1) 1.1 课题的研究背景及意义 (1) 1.1 设计任务 (1) 2 语音信号处理的基本理论知识 (1) 2.1 采样频率和采样位数 (1) 2.2 采样定理 (1) 2.3 IIR数字滤波器 (2) 2.4 FIR数字滤波器 (2) 2.5 IIR数字滤波器和FIR数字滤波器的比较 (3) 2.6 倒谱 (3) 3 语音信号处理和理论方案 (3) 3.1 语音信号的采集 (3) 3.2 语音信号的处理 (4) 3.3 系统框图 (4) 4 语音信号处理的实例分析. (5) 4.1语音文件在MATLAB平台上的录入与打开 (5) 4.2原始语音信号频谱分析及仿真 (5) 4.3加噪语音信号频谱分析及仿真 (8) 4.4去噪及仿真 (13) 4.5 语音信号的回放 (18) 4.6结合去噪后的频谱图对比两种方式滤波的优缺点 (18) 5 总结. (19) 6 致谢. (19) 参考文献 (20)

1绪论 1.1课题的研究背景及意义 语言是我们人类所特有的功能，它是传承和记载人类几千年文明史，没有语言就没有我们今天人类的文明。语音是语言最基本的表现形式，是相互传递信息最重要的手段，是人类最重要、最有效、最常用和最方便的交换信息的形式。 语音信号处理属于信息科学的一个重要分支，大规模集成技术的高度发展和计算机技术的飞速前进，推动了这一技术的发展；它是研究用数字信号处理技术对语音信号进行处理的一门新兴学科，同时又是综合性的多学科领域和涉及面很广的交叉学科，因此我们进行语言信号处理具有时代的意义。 1.2设计任务 本设计先完成语音信号的采集，然后设计低通，带通等滤波器对采集到的语音信号进行滤波处理，分析语音信号各频率段的特性。并对所采集的语音信号加入不同的干扰噪声，对加入噪声的信号进行频谱分析，针对受干扰语音信号的特点设计不同的滤波器，对加噪信号进行滤波，恢复原信号。把原始语音信号、加噪语音信号和滤波后的信号进行时域变换和频域变换，画出它们的时域波形和频域波形图，从视觉角度比较分析滤波的效果。 2语音信号处理的基本理论知识 2.1采样频率和采样位数 采样频率是指计算机每秒钟采集多少个声音样本，是描述声音文件的音质、音调，衡量声卡、声音文件的质量标准。采样频率越高，即采样的间隔时间越短，则在单位时间内计算机得到的声音样本数据就越多，对声音波形的表示也越精确。采样频率与声音频率之间有一定的关系，根据奎斯特理论，只有采样频率高于声音信号最高频率的两倍时，才能把数字信号表示的声音还原成为原来的声音。这就是说采样频率是衡量声卡采集、记录和还原声音文件的质量标准。 采样位数即采样值或取样值，用来衡量声音波动变化的参数，是指声卡在采集和播放声音文件时所使用数字声音信号的二进制位数。采样频率是指录音设备在一秒钟内对声音信号的采样次数，采样频率越高声音的还原就越真实越自然。 采样位数和采样率对于音频接口来说是最为重要的两个指标，也是选择音频接口的两个重要标准。无论采样频率如何，理论上来说采样的位数决定了音频数据最大的力度范围。每增加一个采样位数相当于力度范围增加了6dB。采样位数越多则捕捉到的信号越精确。 2.2 采样定理 在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中，最高频率fmax

现代测试技术习题解答--第二章--信号的描述与分析---副本

第二章 信号的描述与分析 补充题2-1-1 求正弦信号0()sin()x t x ωt φ=+的均值x μ、均方值2 x ψ和概率密度函数 p (x )。 解答： (1)0 00 11lim ()d sin()d 0T T x T μx t t x ωt φt T T →∞== +=? ? ，式中02π T ω = —正弦信号周期 (2) 2 222 2 2 0000 1 1 1cos 2() lim ()d sin ()d d 22 T T T x T x x ωt φψx t t x ωt φt t T T T →∞-+== += = ? ? ? (3)在一个周期内 012ΔΔ2Δx T t t t =+= 000 2Δ[()Δ]lim x x T T T t P x x t x x T T T →∞<≤+=== Δ0Δ000 [()Δ]2Δ2d ()lim lim ΔΔd x x P x x t x x t t p x x T x T x →→<≤+==== 正弦信号 x

2-8 求余弦信号0()sin x t x ωt 的绝对均值x μ和均方根值rms x 。 2-1 求图示2.36所示锯齿波信号的傅里叶级数展开。

2-4周期性三角波信号如图2.37所示，求信号的直流分量、基波有效值、信号有效值及信号的平均功率。

2-1 求图示2.36所示锯齿波信号的傅里叶级数展开。 补充题2-1-2 求周期方波（见图1-4）的傅里叶级数（复指数函数形式），划出|c n|–ω和φn–ω

图，并与表1-1对比。 解答：在一个周期的表达式为 00 (0)2 () (0) 2 T A t x t T A t ? --≤

《测试信号分析与处理》实验报告

《测试信号分析与处理》 实验一差分方程、卷积、z变换 一、实验目的 通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。 二、实验设备 1、微型计算机1台； 2、matlab软件1套 三、实验原理 Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。 差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。用x表示滤波器的输入，用y表示滤波器的输出。 a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1) ak,bk 为权系数，称为滤波器系数。

N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。 y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2)等式定义了数字卷积，*是卷积运算符。输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。 传输函数H(z)是滤波器的第三种实现方法。 H(z)=输出/输入= Y(z)/X(z) （3） 即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。 序列x[n]的z变换定义为 X (z)=∑x[n]z-n (4)把序列x[n] 的z 变换记为Z{x[n]} = X(z)。 由X(z) 计算x[n] 进行z 的逆变换x[n] = Z-1{X(z)}。 Z 变换是Z-1的幂级数，只有当此级数收敛，Z 变换才有意义，而且同一个Z 变换等式，收敛域不同，可以代表不同序列的Z 变换函数。 这三种数字滤波器的表示方法之间可以进行相互转换。 四、实验步骤 1、熟悉matlab软件基本操作指令。读懂下列matlab程序指令，键入程序并 运行，观察运行结果。 Conv.m% 计算两个序列的线性卷积； %-----------------------------------------------------------------

《测试技术与信号处理》习题答案-华科版

《测试技术与信号处理》习题答案 第二章 信号分析基础 1、请判断下列信号是功率信号还是能量信号： （1）)()(10cos 2 ∞<<-∞=t e t x t π （2）)()(||10∞<<-∞=-t e t x t 【解】（1）该信号为周期信号，其能量无穷大，但一个周期内的平均功率有限，属功率信号。 （2）信号能量：? ∞ ∞ -= =10 1 )(2dt t x E ，属于能量信号。 2、请判断下列序列是否具有周期性，若是周期性的，请求其周期。)8 ()(π-=n j e n x 【解】设周期为N ，则有：8 )8 8()()(N j N n j e n x e N n x ?==+-+π 若满足)()(n x N n x =+，则有1)8/sin()8/cos(8/=-=-N j N e jN 即：k N π28/=，k N π16=，k = 0，1，2，3，… N 不是有理数，故序列不是周期性的。 3、已知矩形单脉冲信号x 0(t)的频谱为X 0(ω)=A τsinc(ωτ/2) ，试求图示三脉冲信号的频谱。 【解】三脉冲信号的时域表达式为：)()()()(000T t x t x T t x t x -+++= 根据Fourier 变换的时移特性和叠加特性，可得其频谱： )]cos(21)[2 ( sin )()()()(000T c A e X X e X X T j T j ωωτ τωωωωωω+=++=- 4、请求周期性三角波（周期为T ，幅值为0—A ）的概率分布函数F(x)与概率密度函数p(x) 。 【解】在一个周期T 内，变量x (t )小于某一特定值x 的时间间隔平均值为：T A x t i = ? 取n 个周期计算平均值，当∞→n 时，可有概率分布函数：A x nT t n x F i n =?=∞→lim )( 概率密度函数：A dx x dF x p 1 )()(== t -τ/2 0 τ/2 -T T

数字信号处理 语音信号分析与处理及其MATLAB实现..

摘要 (2) 1 设计目的与要求 (3) 2 设计步骤 (4) 3 设计原理及内容 (5) 3.1 理论依据 (5) 3.2 信号采集 (6) 3.3 构造受干扰信号并对其FFT频谱分析 (8) 3.4 数字滤波器设计 (9) 3.5 信号处理 (10) 总结 (12) 致谢 (13) 参考文献 (14)

用MATLAB对语音信号进行分析与处理，采集语音信号后，在MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。 数字滤波器是数字信号处理的基础，用来对信号进行过滤、检测和参数估计等处理。IIR数字滤波器最大的优点是给定一组指标时，它的阶数要比相同组的FIR滤波器的低的多。信号处理中和频谱分析最为密切的理论基础是傅立叶变换（FT）。离散傅立叶变换（DFT）和数字滤波是数字信号处理的最基本内容。 关键词：MATLAB;语音信号；加入噪声；滤波器；滤波

1. 设计目的与要求 （1）待处理的语音信号是一个在20Hz~20kHz频段的低频信号 （2）要求MATLAB对语音信号进行分析和处理，采集语音信号后，在MATLAB平台进行频谱分析；并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器进行滤除噪声，恢复原信号。

2. 设计步骤 （1）选择一个语音信号或者自己录制一段语音文件作为分析对象； （2）对语音信号进行采样，并对语音信号进行FFT频谱分析，画出信号的时域波形图和频谱图； （3）利用MATLAB自带的随机函数产生噪声加入到语音信号中，对语音信号进行回放，对其进行FFT频谱分析； （4）设计合适滤波器，对带有噪声的语音信号进行滤波，画出滤波前后的时域波形图和频谱图，比较加噪前后的语音信号，分析发生的变化； （5）对语音信号进行回放，感觉声音变化。

音频信号分析与处理

实验三音频信号的分析与处理1 一、实验目的 1.掌握音频信号的采集以及运用Matlab软件实现音频回放的方 法； 2.掌握运用Matlab实现对音频信号的时域、频谱分析方法； 3.掌握运用Matlab设计RC滤波系统的方法； 4.掌握运用Matlab实现对加干扰后的音频信号的进行滤波处理 的方法； 5.锻炼学生运用所学知识独立分析问题解决问题的能力，培养学 生创新能力。 二、实验性质 设计性实验 三、实验任务 1.音频信号的采集 音频信号的采集可以通过Windows自带的录音机也可以用专用的录制软件录制一段音频信号（尽量保证无噪音、干扰小），也可以直接复制一段音频信号，但必须保证音频信号保存为.wav的文件。 2.音频信号的时域、频域分析 运用Matlab软件实现对音频信号的打开操作、时域分析和频域分析，并画出相应的图形（要求图形有标题），并打印在实验报告中（注意：把打印好的图形剪裁下来，粘贴到实验报告纸上）。 3.引入干扰信号 在原有的音频信号上，叠加一个频率为100KHz的正弦波干扰信号（幅度自定，可根据音频信号的情况而定）。 4.滤波系统的设计 运用Matlab实现RC滤波系统，要求加入干扰的音频信号经过RC滤波系统后，能够滤除100KHz的干扰信号，同时保留原有的音频信号，要求绘制出RC滤波系统的冲激响应波形，并分析其频谱。

% 音频信号分析与处理 %% 打开和读取音频文件 clear all; % 清除工作区缓存 [y, Fs] = audioread('jyly.wav'); % 读取音频文件 VoiceWav = y(300000 : 400000, 1); % 截取音频中的一段波形 clear y; % 清除缓存 hAudio = audioplayer(VoiceWav, Fs); % 将音频文件载入audioplayer SampleRate = get(hAudio, 'SampleRate'); % 获取音频文件的采样率KHz T = 1/SampleRate; % 计算每个点的时间，即采样周期SampLen = size(VoiceWav,1); % 单声道采样长度 %% 绘制时域分析图 hFig1 = figure('Units', 'normalized', 'Position', [0 0.05 0.49 0.85]); t = T: T: (SampLen* T); subplot(2, 1, 1); % 绘制音频波形 plot(t, VoiceWav); % 绘制波形 title('音频时域波形图'); axis([0, 2.3, -0.5, 0.5]); xlabel('时间(s)'); ylabel('幅值(V)'); % 显示标题 %% 傅里叶变换 subplot(2, 1, 2); % 绘制波形 myfft(VoiceWav, SampleRate, 'plot'); % 傅里叶变换 title('单声道频谱振幅'); % 显示标题 xlabel('Frequency (Hz)'); ylabel('|Y(f)|'); play(hAudio); % 播放添加噪声前的声音 pause(3); %% 引入100KHz的噪声干扰 t = (0: SampLen-1)* T; noise = sin(2 * pi * 10000 * t); % 噪声频率100Khz，幅值-1V到+1V hFig2 = figure('Units', 'normalized', 'Position', [0.5 0.05 0.5 0.85]); subplot(2, 1, 1); % 绘制波形 plot(t(1: 1000), noise(1: 1000)); title('100KHz噪声信号'); % 显示标题 noiseVoice = VoiceWav+ noise'; % 将噪声加到声音里面 hAudio = audioplayer(noiseVoice, Fs); % 将音频文件载入audioplayer subplot(2, 1, 2); % 绘制波形 [fftNoiseVoice, f] = myfft(noiseVoice, SampleRate, 'plot'); title('音乐和噪声频谱'); % 显示标题 play(hAudio); % 播放添加噪声后的声音 pause(3);

机械工程测试技术基础(第三版)试题(卷)与答案解析集

机械工程测试技术基础(第三版)试卷集. 一、填空题 １、周期信号的频谱是离散的，而非周期信号的频谱是的。 ２、均方值Ψx2表示的是信号的强度，它与均值μx、方差σx2的关系是。 ３、测试信号调理电路主要有、、。 ４、测试系统的静态特性指标有、、。 ５、灵敏度表示系统输出与输入之间的比值，是定度曲线的。 ６、传感器按信号变换特性可分为、。 ７、当时，可变磁阻式电感传感器的输出和输入成近似线性关系，其灵敏度S趋于。 ８、和差特性的主要内容是相临、相反两臂间阻值的变化量符合、的变化，才能使输出有最大值。 ９、信号分析的过程主要包括：、。 10、系统动态特性在时域可用来描述，在复数域可用来描述，在频域可用来描述。 11、高输入阻抗测量放大电路具有高的共模抑制比，即对共模信号有抑制作用，对信号有放大作用。 12、动态应变仪上同时设有电阻和电容平衡旋钮，原因是导线间存在。 13、压控振荡器的输出电压是方波信号，其与输入的控制电压成线性关系。 14、调频波的解调又称，其解调电路称为。 15、滤波器的通频带宽和响应时间成关系。 16、滤波器的频率分辨力主要由其决定。 17、对于理想滤波器，滤波器因数λ＝。 18、带通滤波器可由低通滤波器（f c2）和高通滤波器（f c1）而成（f c2> f c1）。 19、测试系统的线性度和滞后度是由误差引起的；而重复性误差是 由误差引起的。 二、问答题（共30分） １、什么是测试？说明测试系统的构成及各组成部分的作用。（10分） ２、说明电阻丝应变片和半导体应变片的异同点，各有何优点？（10分） ３、选用传感器的原则是什么？（10分） 三、计算题（共55分） １、已知信号x(t)=e-t (t≥0)， (1) 求x(t)的频谱函数X(f)，并绘制幅频谱、相频谱。 (2) 求x(t)的自相关函数R x (τ) 。（15分） ２、二阶系统的阻尼比ξ=0.2，求ω=ωn时的幅值误差和相位误差，如果使幅值误差不大于10％，应取多大阻尼比？。（10分）３、一电容传感器，其圆形极板r = 4mm，工作初始间隙δ0 =0.3mm， （1）工作时如果传感器的工作间隙变化Δδ=±2μm，求电容的变化量。 （2）如果测量电路灵敏度S1=100mv/pF，读数仪表灵敏度S2=5格/mv，在 Δδ=±2μm时，读数仪表的指示值变化多少格？ （ε0 = 8.85×10-12 F/m）(8分) ４、已知RC低通滤波器的R=1KΩ，C=1MF，当输入信号μx= 100sin1000t时， 求输出信号μy 。（7分） ５、（1）在下图中写出动态应变仪所包含的各个电路环节。 （2）如被测量x(t) = sinωt，载波y(t)=sin6ωt，画出各环节信号的波形图。(15分 一、填空题： 1、连续 2、￠x2=H x2+óx2 3、电桥、放大、调制解调电路 4、非线性度、灵敏度、回程误差 5、斜率 6、组合型、一体化型 7、Δó〈〈ó0定位8、相邻相反相对相同9、信号分析、信号处理 10、传递函数、频率函数、脉冲响应函数11、差模12、分布电容13、频率14、鉴频、鉴频器15、反比16、带宽B 17、1 18、串联19、

信号分析与处理

信号分析与处理 第一章绪论：测试信号分析与处理的主要内容、应用；信号的分类，信号分析与信号处理、测试信号的描述，信号与系统。 测试技术的目的是信息获取、处理和利用。 测试过程是针对被测对象的特点，利用相应传感器，将被测物理量转变为电信号，然后，按一定的目的对信号进行分析和处理，从而探明被测对象内在规律的过程。 信号分析与处理是测试技术的重要研究内容。 信号分析与处理技术可以分成模拟信号分析与处理和数字信号分析与处理技术。 一切物体运动和状态的变化，都是一种信号，传递不同的信息。 信号常常表示为时间的函数，函数表示和图形表示信号。 信号是信息的载体，但信号不是信息，只有对信号进行分析和处理后，才能从信号中提取信息。 信号可以分为确定信号与随机信号；周期信号与非周期信号；连续时间信号与离散时间信号；能量信号与功率信号；奇异信号； 周期信号无穷的含义，连续信号、模拟信号、量化信号，抽样信号、数字信号 在频域里进行信号的频谱分析是信号分析中一种最基本的方法：将频率作为信号的自变量，在频域里进行信号的频谱分析； 信号分析是研究信号本身的特征，信号处理是对信号进行某种运算。 信号处理包括时域处理和频域处理。时域处理中最典型的是波形分析，滤波是信号分析中的重要研究内容； 测试信号是指被测对象的运动或状态信息，表示测试信号可以用数学表达式、图形、图表等进行描述。 常用基本信号（函数）复指数信号、抽样函数、单位阶跃函数单位、冲激函数（抽样特性和偶函数）离散序列用图形、数列表示，常见序列单位抽样序列、单位阶跃序列、斜变序列、正弦序列、复指数序列。 系统是指由一些相互联系、相互制约的事物组成的具有某种功能的整体。被测系统和测试系统统称为系统。输入信号和输出信号统称为测试信号。系统分为连续时间系统和离散时间系统。

(完整版)语音信号分析与处理系统设计

语音信号分析与处理系统设计

语音信号分析与处理系统设计 摘要 语音信号处理是研究用数字信号处理技术和语音学知识对语音信号进行处理的新兴的学科，是目前发展最为迅速的信息科学研究领域的核心技术之一。通过语音传递信息是人类最重要、最有效、最常用和最方便的交换信息形式。 Matlab语言是一种数据分析和处理功能十分强大的计算机应用软件，它可以将声音文件变换为离散的数据文件，然后利用其强大的矩阵运算能力处理数据，如数字滤波、傅里叶变换、时域和频域分析、声音回放以及各种图的呈现等，它的信号处理与分析工具箱为语音信号分析提供了十分丰富的功能函数，利用这些功能函数可以快捷而又方便地完成语音信号的处理和分析以及信号的可视化，使人机交互更加便捷。信号处理是Matlab重要应用的领域之一。 本设计针对现在大部分语音处理软件内容繁多、操作不便等问题，采用MATLAB7.0综合运用GUI界面设计、各种函数调用等来实现语音信号的变频、变幅、傅里叶变换及滤波，程序界面简练，操作简便，具有一定的实际应用意义。 最后，本文对语音信号处理的进一步发展方向提出了自己的看法。 关键字：Matlab；语音信号；傅里叶变换；信号处理；

目录 1 绪论 (1) 1.1课题背景及意义 (1) 1.2国内外研究现状 (1) 1. 3本课题的研究内容和方法 (2) 1.3.1 研究内容 (2) 1.3.2 运行环境 (2) 1.3.3 开发环境 (2) 2 语音信号处理的总体方案 (3) 2.1 系统基本概述 (3) 2.2 系统基本要求 (3) 2.3 系统框架及实现 (3) 2.4系统初步流程图 (4) 3 语音信号处理基本知识 (6) 3.1语音的录入与打开 (6) 3.2采样位数和采样频率 (6) 3.3时域信号的FFT分析 (6) 3.4数字滤波器设计原理 (7) 3.5倒谱的概念 (7) 4 语音信号处理实例分析 (8) 4.1图形用户界面设计 (8) 4.2信号的采集 (8) 4.3语音信号的处理设计 (8) 4.3.1 语音信号的提取 (8) 4.3.2 语音信号的调整 (10)

测试信号处理与分析.

2013—2014学年第二学期 《测试信号处理与分析》 实训报告 学院：机械与汽车工程学院 专业：测控技术与仪器 班级：11级测控二班 姓名： 学号： 指导教师

【摘要】：现代信号分析处理技术发展的非常迅速，各种信号专业处理软件也出现在了人们的视野中，这些软件给人们带来了极大方便。本次实训，我们主要学习了INV1612型软件以及DASP信号分析处理系统，切实的感受到了方便。本次实训为期两周，包括在实验室做简支梁的振动信号测试及柔性转子的共振试验，在创新实验室测铣床的振动信号，和在圆楼三楼的AutoCAD机房对铣床的振动信号进行分析等等。在此过程中要基本掌握简支梁的震动信号的测试方法和数据分析，INV1612型多功能柔性转子测试系统、INV1601型振动与控制教学实验系统及MATLAB软件对信号的采集和处理的方法，同时在试验中遇到的问题及我们一起解决的过程。在本次实训中，充分要求了动手能力，实训中的每一项数据都要求自己动手去采集处理，从中我学会了很多知识与方法。 【关键词】：测试信号实训软件知识方法 一、简支梁 1、简支梁的概念 一种简易的支架，包括两个在一平面上可交叉扣合的条形支架，所述每个条形支架的两端为一端高一端低的结构，所述低的一端为钩状结构，钩状体与条形支架主体之间可伸缩的连接，使得每个支架针对不同大小的支撑物在长度方向上可调，且本支架结构简单，节省材料。简支梁就是承载两端竖向荷载，而不提供扭矩的支撑结构。只有两端支撑在柱子上的梁，主要承受正弯矩，一般为静定结构。体系温变、混凝土收缩徐变、张拉预应力、支座移动等都不会在梁中产生附加内力，受力简单，简支梁为力学简化模型。对于简支梁来说，梁的两端搭在两个支撑物上，两端铰接，现实看是只有两端支撑在柱子上的梁，主要承受弯距的单跨结构.一般为静定结构。 2、用“双踪示波比较法”测量简谐振动的频率 用“双踪示波比较法”测量简谐振动的频率，实验仪器有：INV1601B 型振动教学实验仪、INV1601T 型振动教学实验台、速度传感器、调速电机及调压器。

测试技术与信号分析汇总

1．在系统特性测量中常用白噪声信号作为输入信号，然后测量系统的输出，并将输出信号的频谱作为系统频率特性。请用卷积分定理解释这样做的道理。 答：白噪声是指功率谱密度在整个频域内均匀分布的噪声，所有频率具有相同能量的随机噪声称为白噪声。在其频谱上是一条直线。系统频率特性：传递函数的一种特殊情况，是定义在复平面虚轴上的传递函数。时域卷积分定理：两个时间函数的卷积的频谱等于各个时间函数的乘积，即在时域中两信号的卷积等效于在频域中频谱相乘。频域卷积分定理：两个时间函数的频谱的卷积等效于时域中两个时间函数的乘积。y(t)=h(t)*x(t),对y(t)作付式变换，转到相应的频域下Y(f)=H(f)X(f)，由于x(t)是白噪声，付式变换转到频域下为一定值，假定X(f)=1，则有Y(f)=H(f)，此时就是传递函数。 2．用1000Hz的采样频率对200Hz的正弦信号和周期三角波信号进行采样，请问两个信号采样后是否产生混叠？为什么？ 采样频率ωs（2π／Ts）或fs(1／Ts)必须大于或等于信号x(t)中的最高频率ωm的两倍，即ωs＞2ωm，或fs＞2fm。 为了保证采样后的信号能真实地保留原始模拟信号的信息，采样信号的频率必须至少为原信号中最高频率成分的2倍。这是采样的基本法则，称为采样定理。 但在对信号进行采样时，满足了采样定理，只能保证不发生频率混叠，对信号的频谱作逆傅立叶变换时，可以完全变换为原时域采样信号，而不能保证此时的采样信号能真实地反映原信号。工程实际中采样频率通常大于信号中最高频率成分的3到5倍。 理论上周期三角波的频谱里包含所有奇次谐波分量，也就是说200Hz的周期三角波信号包含600Hz、1kHz、1.4kHz等等谐波，所以用1000Hz采样频率对200Hz周期三角波信号采样，会发生混叠。而对200Hz正弦信号采样不会发生混叠。 3．什么是能量泄露和栅栏效应？能量泄漏与栅栏效应之间有何关系？ 能量泄漏：将截断信号的谱XT(ω)与原始信号的谱X（ω）相比较可知，它已不是原来的两条谱线，而是两段振荡的连续谱．这表明原来的信号被截断以后，其频谱发生了畸变，原来集中在f0处的能量被分散到两个较宽的频带中去了，这种现象称之为频谱能量泄漏（Leakage）。 栅栏效应：对采样信号的频谱,为提高计算效率,通常采用FFT算法进行计算,设数据点数为N = T/dt = T.fs则计算得到的离散频率点为Xs(fi) , fi = i.fs/N , i = 0,1,2,…,N/2。这就相当于透过栅栏观赏风景,只能看到频谱的一部分,而其它频率点看不见,因此很可能使一部分有用的频率成分被漏掉,此种现象被称为栅栏效应。 频谱的离散取样造成了栅栏效应，谱峰越尖锐，产生误差的可能性就越大。例如，余弦信号的频谱为线谱。当信号频率与频谱离散取样点不等时，栅栏效应的误差为无穷大。 实际应用中，由于信号截断的原因，产生了能量泄漏，即使信号频率与频谱离散取样点不相等，也能得到该频率分量的一个近似值。从这个意义上说，能量泄漏误差不完全是有害的。如果没有信号截断产生的能量泄漏，频谱离散取样造成的栅栏效应误差将是不能接受的。 能量泄漏分主瓣泄漏和旁瓣泄漏，主瓣泄漏可以减小因栅栏效应带来的谱峰幅值估计误差，有其好的一面，而旁瓣泄漏则是完全有害的。 4．简述传递函数、频响函数和脉冲响应函数间的联系与区别。 传递函数：零初始条件下线性系统响应（即输出）量的拉普拉斯变化（或z变换）与激励（即输入）量的拉普拉斯变换之比。记作G（s）＝Y（s）／U（s），其中Y（s）、U（s)分别为输出量和输入量的拉普拉斯变换。 频响函数：(1)简谐激励时，稳态输出相量与输入相量之比。(2)瞬态激励时，输出的傅里叶变换与输入的傅里叶变换之比。(3)平稳随机激励时，输出和输入的互谱与输入的自谱之比。