对语音信号进行分析及处理资料
- 格式:doc
- 大小:386.98 KB
- 文档页数:17
实验二语音信号的频域特性一、实验目的(1)结合汉语语音信号的各类音素和复元音的特点分析其频域性质;(2)熟悉语音信号的各类音素和复元音的频域参数;(3)熟悉声音编辑软件PRAAT的简单使用和操作。
二、实验记录与思考题1. 观察语音信号的频域特点,总结其规律。
浊音段:其谱线结构是与浊音信号中的周期信号密切相关。
具有与基音及其谐波对应的谱线。
频谱包络中有几个凸起点,与声道的谐振频率相对应。
这些凸起点为共振峰。
清音段:清音的频谱无明显的规律,比较平坦。
2.总结清音/b/p/m/f/d/t/n/l/g/k/h/j/q/x/z/c/s/zh/ch/sh/r/共21个的语谱图的规律,给出辅音的能量集中区;语谱图中的花纹有横杠、乱纹和竖直条。
横杠是与时间轴平行的几条深黑色带纹,它们是共振峰。
从横杠对应的频率和宽度可以确定相应的共振峰频率和带宽。
在一个语音段的语谱图中,有没有横杠出现是判断它是否是浊音的重要标志。
竖直条是语谱图中出现于时间轴垂直的一条窄黑条。
每个竖直条相当于一个基音,条纹的起点相当于声门脉冲的起点,条纹之间的距离表示基音周期,条纹越密表示基音频率越高。
b,p……清音的语谱图为乱纹。
辅音的能量集中区为:高频区3. 总结浊音/a/o/e/i/u/ü/ao/ai/ei/ou/ie /an/en/in/ang/eng/ong/ing/共18个的语谱图的规律,提取这18个浊音的基频、前三个共振峰频率4./r/、/m/、/n/、/l/ 从这几个音素的的基频、共振峰频率5.分析宽带语谱图和窄带语谱图的不同之处,请解释原因;语谱图中的花纹有横杠、乱纹和竖直条等。
横杠是与时间轴平行的几条深黑色带纹,它们是共振峰。
从横杠对应的频率和宽度可以确定相应的共振峰频率和带宽。
在一个语音段的语谱图中,有没有横杠出现是判断它是否是浊音的重要标志。
竖直条(又叫冲直条)是语谱图中出现与时间轴垂直的一条窄黑条。
每个竖直条相当于一个基音,条纹的起点相当于声门脉冲的起点,条纹之间的距离表示基音周期。
语音信号分析与语音合成技术研究近年来,语音信号分析和语音合成技术的研究得到了越来越多的关注。
语音作为一种个性化的沟通方式,在现代社会扮演着重要的角色。
语音识别、语音合成、情感分析等多个领域的技术也在不断发展。
本文将从以下几个方面对语音信号分析和语音合成技术进行探讨。
一、语音信号分析技术语音信号作为一种复杂的信号,其分析和处理涉及到很多方面的知识。
对于语音信号的分析技术,主要可以分为以下几类:1.语音信号的预处理语音信号的预处理其实就是对语音信号的一些附加处理。
比如,对于语音信号的容量比较大,需要将其进行降采样,以减少数据的冗余,同时避免由于数据过多导致的计算误差等。
还可以进行去噪处理,去除环境噪声的影响,使语音信号更加干净、清晰。
2.语音信号的特征提取语音信号是很复杂的信号,其包含了很多信息,如音调、音频、元音等。
因此,在语音信号分析中,在预处理之后,还需要对信号进行特征提取,以便于对语音信号进行分类、识别等处理。
因此,常用的特征提取算法包括:线性预测系数(LPC)、短时傅里叶转换(STFT)、梅尔频率倒谱系数(MFCC)等。
3.语音信号的分类和识别语音信号的特征提取之后,就可以使用分类和识别算法,对语音信号进行分类和识别了。
常用的语音信号分类和识别算法包括隐马尔可夫模型(HMM)、高斯混合模型(GMM)等。
其中,HMM 是一种基于概率的模型,常被用于语音识别中,其主要依靠模型的转移概率、状态概率和观测概率等进行识别分析。
4.语音信号的情感分析由于语音信号不仅包含声音信号,还包含了丰富的文化和情感信息,因此语音信号的情感分析也成为了研究的热点。
情感分析主要是通过语音表现的语言特征,比如声调、语调、语速、音节长度等,对说话人的情感状态进行分类和分析。
二、语音合成技术对于语音合成技术,主要包括如下几个方面的研究:1.语音合成技术的原理语音合成技术主要是将文本转化为语音的过程。
其工作原理是通过提取文本的语音特征,并利用这些特征合成出逼真的人类语音。
语音信号处理语音信号处理是对语音信号进行分析、处理和合成的一种技术。
随着和语音识别技术的快速发展,语音信号处理变得越来越重要。
本文将详细介绍语音信号处理的基本概念、常用技术和应用领域。
基本概念语音信号是指人类通过声音来交流的方式。
语音信号通常采用模拟信号的形式,通过麦克风传感器转换为数字信号,然后使用数字信号处理技术进行分析和处理。
语音信号的特点包括频率、幅度和时域特性。
常用技术预处理语音信号预处理是指在进行语音信号分析和处理之前,对原始语音信号进行预处理以提取和增强感兴趣的特征。
常用的预处理技术包括去噪、滤波、降低共振、归一化等。
特征提取特征提取是从语音信号中提取有用信息的过程,目的是将语音信号转化为可以被机器学习算法处理的形式。
常用的特征包括声谱图、梅尔倒谱系数(MFCC)、线性预测编码(LPC)等。
语音识别语音识别是将语音信号转化为文字或命令的过程。
常用的语音识别技术包括基于模板的方法、隐马尔可夫模型(HMM)、深度学习等。
语音合成语音合成是将文字转化为语音信号的过程。
常用的语音合成技术包括基于拼接的方法、隐马尔可夫模型(HMM)、深度学习等。
应用领域语音信号处理在许多领域中起着重要作用,以下是几个主要应用领域的例子:语音识别系统语音识别系统可以用于实现语音自动接听、语音搜索等应用。
这些系统通过对输入语音信号进行处理和分析,将其转化为文字或命令。
语音合成系统语音合成系统可以将文字转化为语音,实现自动语音播报、电子书朗读等功能。
这些系统通过将输入文本处理和合成为语音信号。
声纹识别系统声纹识别系统通过对语音信号进行处理和分析,将其转化为声纹特征,用于实现语音身份认证等应用。
噪声抑制噪声抑制是指对含噪声的语音信号进行处理,去除噪声以提高语音信号的质量。
语音压缩语音压缩是将语音信号进行压缩,以减小存储空间和传输带宽的需求。
语音压缩技术可以用于语音通信、语音存档等应用。
语音信号处理是一项重要的技术,它在语音识别、语音合成和其他领域中有广泛的应用。
语音信号分析与处理方法综述语音信号是人类交流中重要的组成部分,语音信号分析与处理是一项非常重要的研究方向。
语音信号分析与处理技术可以被应用于语音识别、语音合成、语音增强、语音鉴别等领域,为人类提供更高效、高质量的语音交流体验。
本文将综述语音信号分析与处理的方法,包括基本概念和基础知识,语音信号的预处理、特征提取与分类等技术,以及目前应用较广泛的语音识别和语音合成技术。
本文还将介绍这些技术的理论基础、算法和应用案例。
一. 语音信号的基本概念和基础知识语音信号是人类通过声音来进行交流的信号,包含语音的各种音节和声调信息。
在语音信号分析与处理中,需要掌握语音信号的基本概念和基础知识。
在语音信号中,有许多不同的声音成分和噪声成分,语音信号预处理技术采取一系列措施来消除这些噪声,提高语音信号的质量。
例如,预处理技术包括消噪、滤波、去混响等技术。
二. 语音信号的预处理技术语音信号预处理技术包括噪声去除、滤波、去混响等技术。
这些技术可以帮助清除语音信号中的杂音和噪声,提高语音信号的可读性和质量。
噪声去除技术是预处理语音信号中最常用和最简单的技术。
这种技术采用滤波器、降噪器等设备或软件来消除语音信号中的杂音和噪声。
滤波技术是在语音信号中过滤特定频率成分的一种技术,这可以将一些不必要的频率成分去除,使语音信号更加清晰。
滤波依据的原理是,去除不相关信号的频率成分,使有用信号能够更好地被保留。
去混响技术是将混响(残留回响)去除的技术。
在语音信号分析与处理中,经常需要去除混响,因为混响会导致语音信号的低频部分(主要是语音信号中的音节)发生失真,降低语音信号的质量。
三. 语音信号特征提取语音信号的特征提取是将语音信号转化为一种易于分析和识别的数据结构或特征向量的过程。
一些常用的语音信号特征提取方法包括短时能量、过零率、梅尔倒谱系数等。
短时能量是一种用于计算语音信号能量的技术。
通过计算语音信号每个短时间内的总能量,可以得出一个很好的特征向量。
一、设计目的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个采样点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 f c 1000=,阻带截止频率Hz f st 1200=,通带最大衰减dB 11=δ,阻带最小衰减dB 1002=δ。
语音信号处理实验报告——语音信号分析实验一.实验目的及原理语音信号分析是语音信号处理的前提和基础,只有分析出可表示语音信号本质特征的参数,才有可能利用这些参数进行高效的语音通信、语音合成和语音识别等处理,并且语音合成的音质好坏和语音识别率的高低,都取决于对语音信号分析的准确性和精确性;贯穿语音分析全过程的是“短时分析技术”;因为从整体来看,语音信号的特性及表征其本质特征的参数均是随时间变化的,所以它是一个非平稳态过程,但是在一个短时间范围内一般认为在10~30ms的时间内,其特性基本保持不变,即相对稳定,可将其看做一个准稳态过程,即语音信号具有短时平稳性;所以要将语音信号分帧来分析其特征参数,帧长一般取为10ms~30ms;二.实验过程男声及女声蓝色为时域信号,红色为每一帧的能量,绿色为每一帧的过零率某一帧的自相关函数3.频域分析①一帧信号的倒谱分析和FFT及LPC分析②男声和女声的倒谱分析对应的倒谱系数:,,……对应的LPC预测系数:1,,,,,……原语音波形一帧语音波形一帧语音的倒谱③浊音和清音的倒谱分析④浊音和清音的FFT分析和LPC分析红色为FFT图像,绿色为LPC图像三.实验结果分析1.时域分析实验中采用的是汉明窗,窗的长度对能否由短时能量反应语音信号的变化起着决定性影响;这里窗长合适,En能够反应语音信号幅度变化;同时,从图像可以看出,En可以作为区分浊音和清音的特征参数;短时过零率表示一帧语音中语音信号波形穿过横轴零电平的次数;从图中可以看出,短时能量和过零率可以近似为互补的情况,短时能量大的地方过零率小,短时能量小的地方过零率较大;从浊音和清音的时域分析可以看出,清音过零率高,浊音过零率低;从男声女声的时域信号对比图中可以看出,女音信号在高频率分布得更多,女声信号在高频段的能量分布更多,并且女声有较高的过零率,这是因为语音信号中的高频段有较高的过零率;2.频域分析这里对信号进行快速傅里叶变换FFT,可以发现,当窗口函数不同,傅里叶变换的结果也不相同;根据信号的时宽带宽之积为一常数这一性质,可以知道窗口宽度与主瓣宽度成反比,N越大,主瓣越窄;汉明窗在频谱范围中的分辨率较高,而且旁瓣的衰减大,具有频谱泄露少的有点,所以在实验中采用的是具有较小上下冲的汉明窗;为了使频域信号的频率分辨率较高,所取的DFT及相应的FFT点数应该足够多,但时域信号的长度受到采样率和和短时性的限制,这里可以采用补零的办法,对补零后的序列进行FFT变换;从实验仿真图可以看出浊音的频率分布比清音高;3.倒谱分析通过实验可以发现,倒谱的基音检测与语音加窗的选择也是有关系的;如果窗函数选择矩形窗,在许多情况下倒谱中的基音峰将变得不清晰,窗函数选择汉明窗较为合理,可以发现,加汉明窗的倒谱基音峰较为突出;在典型的浊音清音倒谱对比中,理论上浊音倒谱基音峰应比较突出,而清音不出现这种尖峰,只是在倒谱的低时域部分包含关于声道冲激响应的信息;实验仿真的图形不是很理想;4.线性预测分析从实验中可以发现,LPC谱估计具有一个特点,在信号能量较大的区域即接近谱的峰值处,LPC谱和信号谱很接近;而在信号能量较低的区域即接近谱的谷底处,则相差比较大;在浊音清音对比中,可以发现,对呈现谐波特征的浊音语音谱来说这个特点很明显,就是在谐波成分处LPC谱匹配信号谱的效果要远比谐波之间好得多;在实验中,当P值增加到一定程度,预测平方误差的改善就不很明显了,而且会增加计算量,一般取为8~14,这里P取为10;5.基音周期估计①自互相关函数法②短时平均幅度差法③倒谱分析法共偏移92+32=124个偏移点16000/124=可以发现,上面三种方法计算得到的基音周期基本相同;。
语音信号分析与处理技术语音信号是人类日常交流的重要方式之一,是信息交流和社交活动不可或缺的一部分。
因此,对于语音信号的分析和处理技术的研究具有非常重要的意义。
语音信号处理技术广泛应用于语音识别、语音合成、语音增强、音频编码等领域。
一、语音信号的特征分析语音信号是一种时间和频率的交叉信号,因此,对于语音信号的分析首先要对其进行频域和时域分析。
1.1 时域分析时域分析是将信号分解为相邻时刻的振幅和相位的变化。
时域分析的主要方法是进行窗口函数分析。
常用的窗口函数有矩形窗、汉宁窗、海宁窗、Hamming窗等,窗口函数的选择取决于分析信号的性质。
例如,当信号的重点在时域中的相关性时,计算时输入的窗口函数通常是矩形窗;当需要打包使用的声音时,通常选择海宁窗。
1.2 频域分析频域分析是将语音信号分解为频率块和其振幅的变化。
频域分析的主要方法是通过傅里叶变换将信号转化为频率表示。
在频域分析中,常用的指标包括频率谱、频率曲线和傅里叶系数等。
傅里叶系数的使用能很好地描述信号的频率成分,而且可以用于合成声音的频率表示。
而频率曲线则主要适用于掌握声音的过渡特性。
二、语音信号的处理技术2.1 声音增强声音增强是指通过处理技术,对有噪声的语音信号进行低噪声处理和消除噪声信号。
声音增强的原理是把信号分成干净的声音和噪声两个部分,然后对噪声进行滤波或抑制。
在语音信号增强技术中,常用的技术包括算法去噪、谱减法等。
算法去噪是根据噪声特殊的相关关系,来推测信号的特殊性质。
而谱减法则是基于最小化信号频谱变化来减少噪声信号。
2.2 语音识别语音识别是将说话人的语音转换为文字或机器指令。
语音识别技术在人机交互中发挥着极为重要的作用。
语音识别的基本原理是语音信号的特征提取和建模。
通过将语音数据转换为数字信号,并从数字数据中提取特征,进而应用分类器来实现识别。
常用的语音识别技术包括基于高斯混合模型(GMM)的识别技术、基于隐马尔可夫模型(HMM)的识别技术等。
语音信号处理简介语音信号处理是一种通过对语音信号进行分析、处理和合成的技术,以提取语音中的有用信息并改善语音质量。
它在语音识别、语音合成、语音增强等领域中有着广泛的应用。
本文将介绍语音信号处理的基本概念、常见的处理方法以及应用场景。
基本概念语音信号语音信号是由人类语音产生的声波信号,它是一种时间变化的波形信号。
语音信号包含了说话人的身份特征、语义信息以及情感特征等。
在语音信号处理中,通常使用数字信号来表示和处理语音信号。
语音信号的特性语音信号具有多种特性,包括频域特性和时域特性。
频域特性频域特性描述了语音信号在频率上的分布情况。
常见的频域特性包括频谱、功率谱和频带能量等。
频域特性能够反映语音信号中存在的不同频率成分。
时域特性时域特性描述了语音信号在时间上的变化情况。
常见的时域特性包括时域波形、自相关函数和短时能量等。
时域特性能够反映语音信号的时序关系。
常见的语音信号处理方法语音信号处理涉及到多种方法和技术,下面介绍几种常见的处理方法。
预处理预处理是语音信号处理的第一步,它主要用于降噪和增强语音信号的质量。
常见的预处理方法包括滤波、降噪和增益控制等。
滤波滤波是一种通过选择性地传递或阻止不同频率成分的方法。
在语音信号处理中,常用的滤波器包括低通滤波器和高通滤波器等。
降噪降噪是一种通过消除语音信号中的噪声成分来提高语音质量的方法。
常见的降噪方法包括谱减法、小波降噪和自适应滤波等。
增益控制增益控制是一种通过调整语音信号的幅度来平衡不同部分的能量的方法。
常见的增益控制方法包括自动增益控制(AGC)和压缩扩展(Compand)等。
特征提取特征提取是语音信号处理中最重要的环节之一,它用于从语音信号中提取有用的特征信息。
常见的特征提取方法包括短时能量、短时过零率和线性预测系数等。
语音识别是一种将语音信号转换为文本或命令的技术。
它在语音助手、语音控制和语音翻译等领域中有着广泛的应用。
常见的语音识别方法包括基于模型的方法和基于深度学习的方法等。
一、设计目的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个采样点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 f c 1000=,阻带截止频率Hz f st 1200=,通带最大衰减dB 11=δ,阻带最小衰减dB 1002=δ。
语音信号处理技术及应用
语音信号处理技术是指通过对语音信号进行分析、提取和处理,以达到对语音信号的识别、压缩、增强、转换等各种应用需求。
语音信号处理技术的一些常见方法和算法包括:
1. 语音信号的数字化:将模拟语音信号转换为数字形式,通常使用采样和量化技术。
2. 语音信号的预处理:对于中断、噪声等干扰,可以利用滤波、去噪、增强等方法进行预处理。
3. 语音信号的特征提取:通过对语音信号进行分析,提取出特定的特征参数,如短时能量、频率轮廓、基频、共振峰等。
4. 语音信号的模型建立:通过统计模型、混合高斯模型等方法,对语音信号进行建模,提取语音的概率模型。
5. 语音信号的识别:利用概率模型,将输入的语音信号与预先训练好的模型进行匹配,以实现语音信号的识别。
语音信号处理技术在很多领域都有应用,包括但不限于以下几个方面:
1. 语音识别:利用语音信号处理技术,将输入的语音信号转换为文本。
2. 语音合成:根据文本信息,利用语音信号处理技术生成对应的语音信号。
3. 语音增强:通过去除噪声、增强语音信号,提高语音信号的质量。
4. 语音压缩:将语音信号进行压缩以减少存储空间或传输带宽。
5. 语音转换:将语音信号转换为不同的声音特征,例如男性声转女性声。
语音信号处理技术在语音识别、语音合成、语音增强、语音压缩等领域都发挥着重要的作用,并且在实际应用中已经取得了很大的成果。
语音信号处理与分析语音信号处理与分析是数字信号处理领域的一个重要分支。
它涉及了对语音信号的各种处理技术和分析方法。
语音信号处理与分析的主要目标是提取和控制语音信号中的有用信息,以实现语音识别、语音合成、语音增强、语音编码等一系列语音相关应用。
一、语音信号特点语音信号是人类沟通中最基本的形式之一。
它具有以下几个基本特点:1. 声音频率范围广泛:人类能够听到的声音频率范围约为20Hz到20kHz。
而语音信号一般集中在300Hz到4kHz之间,这个频率范围包含了语音的大部分信息。
2. 时域相关性强:语音信号在时域上呈现出一定的连续性,即相邻时间点的样本值之间存在一定的相关性。
3. 信息量大:语音信号中包含了大量的语义、语法和语音音素信息,涵盖了人类语言交流的各个层面。
二、语音信号处理语音信号处理旨在提取和改善语音信号中的信息,使其更易于分析和理解。
常见的语音信号处理技术包括:1. 语音预处理:对原始语音信号进行降噪、去除回声、均衡化等处理,以增强语音的清晰度和可听性。
2. 特征提取:通过对语音信号进行时频分析,提取出与语音内容相关的特征参数,如短时能量、过零率、共振峰频率等。
3. 语音编码:将语音信号以压缩形式存储或传输,以减少存储空间和传输带宽。
常用的语音编码算法有PCM、ADPCM、MP3等。
4. 语音识别:通过计算机对语音信号进行自动识别,将语音转化为文字。
语音识别广泛应用于语音助手、语音搜索等领域。
5. 语音合成:根据输入的文字信息,生成与人类声音相似的合成语音。
语音合成的应用包括语音助手、有声阅读、机器人交互等。
三、语音信号分析语音信号分析旨在从语音信号中提取有关语音的信息,以揭示语音产生机制和语音特征。
常见的语音信号分析方法包括:1. 短时傅里叶变换(STFT):将语音信号按时间窗进行分段,对每个时间窗进行傅里叶变换,得到时间频率分布谱。
2. 线性预测编码(LPC):通过建立线性预测模型,提取出语音信号中的共振峰频率和预测残差。
语音信号处理实验报告 The Standardization Office was revised on the afternoon of December 13, 2020语音信号处理实验报告——语音信号分析实验一.实验目的及原理语音信号分析是语音信号处理的前提和基础,只有分析出可表示语音信号本质特征的参数,才有可能利用这些参数进行高效的语音通信、语音合成和语音识别等处理,并且语音合成的音质好坏和语音识别率的高低,都取决于对语音信号分析的准确性和精确性。
贯穿语音分析全过程的是“短时分析技术”。
因为从整体来看,语音信号的特性及表征其本质特征的参数均是随时间变化的,所以它是一个非平稳态过程,但是在一个短时间范围内(一般认为在10~30ms的时间内),其特性基本保持不变,即相对稳定,可将其看做一个准稳态过程,即语音信号具有短时平稳性。
所以要将语音信号分帧来分析其特征参数,帧长一般取为10ms~30ms。
二.实验过程男声及女声(蓝色为时域信号,红色为每一帧的能量,绿色为每一帧的过零率)某一帧的自相关函数3.频域分析①一帧信号的倒谱分析和FFT及LPC分析②男声和女声的倒谱分析对应的倒谱系数:,,……对应的LPC预测系数:1,,,,,……原语音波形一帧语音波形一帧语音的倒谱③浊音和清音的倒谱分析④浊音和清音的FFT分析和LPC分析(红色为FFT图像,绿色为LPC图像)三.实验结果分析1.时域分析实验中采用的是汉明窗,窗的长度对能否由短时能量反应语音信号的变化起着决定性影响。
这里窗长合适,En能够反应语音信号幅度变化。
同时,从图像可以看出,En可以作为区分浊音和清音的特征参数。
短时过零率表示一帧语音中语音信号波形穿过横轴(零电平)的次数。
从图中可以看出,短时能量和过零率可以近似为互补的情况,短时能量大的地方过零率小,短时能量小的地方过零率较大。
从浊音和清音的时域分析可以看出,清音过零率高,浊音过零率低。
从男声女声的时域信号对比图中可以看出,女音信号在高频率分布得更多,女声信号在高频段的能量分布更多,并且女声有较高的过零率,这是因为语音信号中的高频段有较高的过零率。
实验三:语音信号谱分析及去噪处理1、实验目的(1)通过对实际采集的语音信号进行分析和处理,获得数字信号处理实际应用的认识。
(2)掌握数字信号谱分析的知识。
(3)掌握数字滤波器设计的知识,并通过对语音信号的去噪处理,获得数字滤波器实际应用的知识。
2、实验内容(1)用麦克风自行采集两段语音信号[高频噪声、人声+高频噪声](.wav格式)。
(2)通过Matlab读入采集信号,观察其采样频率,并绘图采样信号。
(3)通过Matlab对语音信号进行谱分析,分析出噪声的频带。
(4)设计一滤波器,对叠加入噪声的语音信号进行去噪处理。
绘图并发声去噪后的信号。
3、实验步骤(1)利用麦克风采集一段5s以内的语音信号。
利用格式工厂软件对语音信号进行预处理。
通常语音信号为单声道,采样频率为8000Hz,语音信号为.wav格式。
(2)通过Matlab读入语音信号及其采样频率(使用Matlab库函数wavread),在Matlab软件的workspace工作平台上观察读入的语音信号,在Matlab中,对入的语音信号为一维矩阵。
应注意,库函数wavread自动将语音信号幅度归一化[-1,1]区间范围。
使用Matlab库函数plot 绘图语音信号,并使用库函数sound发音语音信号。
(3)分析噪声的频谱。
在这里进行谱分析的目的,是了解噪声信号的频谱特性,为去噪滤波器的技术指标提供依据。
(4)通过Matlab对语音信号进行谱分析。
应注意,对信号进行谱分析,在实验一中已经详细介绍过。
在这里进行谱分析的目的,是了解本段语音信号的频谱特性,为去噪滤波器的技术指标提供依据。
(5)根据语音信号及噪声信号的频谱特性,自行设计一滤波器,对叠加入噪声的语音信号进行去噪处理。
最后绘图并发声去噪后的信号。
应注意,数字滤波器的实际应考虑实际需求,合理制定滤波器的技术指标。
4、实验原理用麦克风采集一段语音信号,绘制波形并观察其频谱,添加一段随机信号,给定相应的滤波器指标,用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器,对该语音信号进行滤波去噪处理,比较滤波前后的波形和频谱并进行分析。
语音信号的分析与处理在日常生活中,我们常常与语音信号打交道,如电话通讯、语音识别、音乐播放、影片配音等。
语音信号是人类声音的一种电信号表示形式,它的特性是非常复杂的,包括语音的声音频率、幅度大小、声音的拐角特征、信号的频率变化以及背景噪声等多方面因素。
因此,对语音信号正确分析和处理是现代通讯研究、智能语音识别和场景识别等领域的重要问题,也是一个迫切需要解决的问题。
一、语音信号的基本特征语音信号具有很多特征,如频率、能量、音色和韵律等。
频率是语音信号的基本特征之一,它是指语音信号中声音的高低频率。
音频信号的波形形状与频率息息相关。
一般组成语音的基元元音频率范围在250 Hz ~ 1000 Hz之间,辅音频率的范围在100 Hz ~ 4 kHz之间。
能量与音量相关,是指语音信号所含有的总能量。
音色是语音信号的另一个特征,它能够指示语音信号的来源。
最后,韵律则是指语音信号的节奏,其包含语音中音节、词语、句子和语气的信息。
二、语音信号的分析方法语音信号的分析方法可以分为时域分析和频域分析两种方式。
其中,时域分析是一种基于时间的分析,它通过观察信号的实时波形来分析语音信号的特征。
频域分析则是一种基于频率的分析,它通过观察信号的频谱特性来分析语音信号的频率、音量和音色。
1. 时域分析时域分析是一种非常基础的语音信号分析方法。
通常,时域分析方法通过分析语音信号的波形特征来判断语音信号的特点。
它能够检查信号在时间上的变化,比如分析语音信号中频率与振幅的变化。
这种方法主要通过时间和采样频率来确定语音信号的基本特征。
2. 频域分析频域分析则是一种付于注意的语音信号分析方法。
它主要通过傅里叶变换(Fourier transform)或小波变换(Wavelet transform)等频率分析方法来研究信号在不同频段上的特征。
通过频域分析可以获得信号在较高频段上的信息,该信息往往无法通过时域分析方法获得。
频域分析方法可以用于语音信号的分析、信号噪声抑制和语音信号质量改进等方面。
一、实验目的1. 理解语音信号的基本特性和处理方法。
2. 掌握语音信号的采样、量化、编码等基本过程。
3. 学习使用相关软件对语音信号进行时域和频域分析。
4. 了解语音信号的降噪、增强和合成技术。
二、实验原理语音信号是一种非平稳的、时变的信号,其频谱特性随时间变化。
语音信号处理的基本过程包括:信号采集、信号处理、信号分析和信号输出。
三、实验仪器与软件1. 仪器:计算机、麦克风、耳机。
2. 软件:Matlab、Audacity、Python。
四、实验步骤1. 信号采集使用麦克风采集一段语音信号,并将其存储为.wav格式。
2. 信号处理(1)使用Matlab读取.wav文件,提取语音信号的采样频率、采样长度和采样数据。
(2)将语音信号进行时域分析,包括绘制时域波形图、计算信号的能量和过零率等。
(3)将语音信号进行频域分析,包括绘制频谱图、计算信号的功率谱密度等。
3. 信号分析(1)观察时域波形图,分析语音信号的幅度、频率和相位特性。
(2)观察频谱图,分析语音信号的频谱分布和能量分布。
(3)计算语音信号的能量和过零率,分析语音信号的语音强度和语音质量。
4. 信号输出(1)使用Audacity软件对语音信号进行降噪处理,比较降噪前后的效果。
(2)使用Python软件对语音信号进行增强处理,比较增强前后的效果。
(3)使用Matlab软件对语音信号进行合成处理,比较合成前后的效果。
五、实验结果与分析1. 时域分析从时域波形图可以看出,语音信号的幅度、频率和相位特性随时间变化。
语音信号的幅度较大,频率范围一般在300Hz~3400Hz之间,相位变化较为复杂。
2. 频域分析从频谱图可以看出,语音信号的能量主要集中在300Hz~3400Hz范围内,频率成分较为丰富。
3. 信号处理(1)降噪处理:通过对比降噪前后的时域波形图和频谱图,可以看出降噪处理可以显著降低语音信号的噪声,提高语音质量。
(2)增强处理:通过对比增强前后的时域波形图和频谱图,可以看出增强处理可以显著提高语音信号的幅度和频率,改善语音清晰度。
语音信号处理实验报告语音信号处理实验报告一、引言语音信号处理是一门研究如何对语音信号进行分析、合成和改善的学科。
在现代通信领域中,语音信号处理起着重要的作用。
本实验旨在探究语音信号处理的基本原理和方法,并通过实验验证其有效性。
二、实验目的1. 了解语音信号处理的基本概念和原理。
2. 学习使用MATLAB软件进行语音信号处理实验。
3. 掌握语音信号的分析、合成和改善方法。
三、实验设备和方法1. 设备:计算机、MATLAB软件。
2. 方法:通过MATLAB软件进行语音信号处理实验。
四、实验过程1. 语音信号的采集在实验开始前,我们首先需要采集一段语音信号作为实验的输入。
通过麦克风将语音信号输入计算机,并保存为.wav格式的文件。
2. 语音信号的预处理在进行语音信号处理之前,我们需要对采集到的语音信号进行预处理。
预处理包括去除噪声、归一化、去除静音等步骤,以提高后续处理的效果。
3. 语音信号的分析语音信号的分析是指对语音信号进行频谱分析、共振峰提取等操作。
通过分析语音信号的频谱特征,可以了解语音信号的频率分布情况,进而对语音信号进行进一步处理。
4. 语音信号的合成语音信号的合成是指根据分析得到的语音信号特征,通过合成算法生成新的语音信号。
合成算法可以基于传统的线性预测编码算法,也可以采用更先进的基于深度学习的合成方法。
5. 语音信号的改善语音信号的改善是指对语音信号进行降噪、增强等处理,以提高语音信号的质量和清晰度。
常用的语音信号改善方法包括时域滤波、频域滤波等。
六、实验结果与分析通过实验,我们得到了经过语音信号处理后的结果。
对于语音信号的分析,我们可以通过频谱图观察到不同频率成分的分布情况,从而了解语音信号的特点。
对于语音信号的合成,我们可以听到合成后的语音信号,并与原始语音信号进行对比。
对于语音信号的改善,我们可以通过降噪效果的评估来判断处理的效果。
七、实验总结通过本次实验,我们深入了解了语音信号处理的基本原理和方法,并通过实验验证了其有效性。
一、设计目的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个采样点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 f c 1000=,阻带截止频率Hz f st 1200=,通带最大衰减dB 11=δ,阻带最小衰减dB 1002=δ。
(2)FIR 低通滤波器性能指标通带截止频率Hz f c 1000=,阻带截止频率Hz f st 1200=, 通带衰减1δ≤1dB ,阻带衰减2δ≥100dB 。
(3)IIR 高通滤波器的设计指标,Hz f z 1000=,Hz f p 2000=,阻带最小衰减dB A s 30=,通带最大衰减dB A P 1=。
(4)(4)FIR 高通滤波器的设计指标,Hz f z 1000=,Hz f p 2000=,阻带最小衰减dB A s 50=,通带最大衰减dB A P 1=。
(5)用自己设计的各滤波器分别对采集的信号进行滤波,在Matlab 中,FIR 滤波器利用函数fftfilt 对信号进行滤波,IIR 滤波器利用函数filter 对信号进行滤波。
比较滤波前后语音信号的波形及频谱,在一个窗口同时画出滤波前后的波形及频谱。
在Matlab中,函数sound可以对声音进行回放。
其调用格式:sound(x,fs,bits);可以感觉滤波前后的声音有变化。
三、结果分析1、原始语音信号采样后的时域波形及FFT变换后频谱图 1.原始信号时域图图 2.原始信号的频谱图2、加入噪声后信号的时域波形及FFT变换后频谱图 3.加入噪声后的时域信号分析:由图可以看出加入噪声后有明显的不一样,运行sound(x2 ,fs);播放加入高斯噪声后的语音,可以明显听出兹兹的噪声图 4.加入噪声后的频谱图3、IIR滤波器及原始信号通过IIR后的时域波形频谱变化图 5. IIR低通滤波器图 6. IIR低通滤波器滤波前后时域波形分析:经过滤波器后的信号和原始信号很近似图7. IIR低通滤波器滤波前后的频谱分析:从图7可以看出,经过IIR低通滤波器滤波后将高频部分滤除了。
4、FIR滤波器及原始信号通过FIR后的时域波形频谱变化图8. FIR低通滤波器图9.FIR低通滤波器滤波前后时域波形图10.FIR低通滤波器滤波前后频谱分析:从图10可以看出,经过FIR低通滤波器滤波后将高频部分滤除了。
5、IIR高通滤波器的设计图11. IIR高通滤波器图12. IIR高通滤波器滤波前后时域波形图13. IIR高通滤波器滤波前后频谱6、FIR高通滤波器的设计图14. FIR高通滤波器图15. FIR高通滤波器滤波前后时域波形图16.FIR高通滤波器滤波前后频谱四、结束语这次的数字信号处理大作业的题目是应用Matlab对语音信号进行频谱分析及滤波,首先通过网络和书籍查找有关本次作业所需的资料,编写相关程序,并通过Matlab软件运行得到相关波形频谱图。
在做作业的过程中,我将上课所学的理论知识运用到实践中。
通过这次应用Matlab对语音信号进行频谱分析及滤波,让我对Matlab的应用以及数字滤波器的设计有了更深层次的理解,每个程序中的语句表示什么意思也有了很清楚的了解。
在实践中增强了我的动手能力,并提高了我的综合能力,使自身得到了很大的锻炼。
另外,在设计滤波器的过程中由于个人知识学得不到位,后面的仿真结果不是很理想,我希望以后多查阅资料,多积累,多思考,只有这样,才能取得更大的进步,才能学有所用,学有所长。
五、程序附录(1)采样+噪声clearclc[x1,fs,bits]=wavread('E:\数字信号.wav'); %读取语音信号的数据,赋给变量x1 sound(x1,fs); %播放原始语音信号N=length(x1)fs %采样频率为44100Hzbits %比特率为16x=x1(60001:1:120000); %截取原始信号60000个采样点N1=length(x)figure(1)plot(x) %做截取原始信号的时域图形title('原始语音采样后时域信号');xlabel('时间轴 n');ylabel('幅值 A');sound(x,fs); %播放截取后语音信号,仍能清晰地听到“数字信号”k=0:5999;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('相位谱');mean(x) %求得语音信号平均幅值-2.1368e-04d=randn(1,60000); %产生高斯白噪声mean(d)d=d/100; %对噪声进行衰减x2=x+d; %加入高斯白噪声sound(x2 ,fs); %播放加入高斯噪声后的语音,可以明显听出噪声figure(3)plot(x2)title('加入噪声后时域信号');xlabel('时间轴 n');ylabel('幅值 A');k=0:5999;y2=fft(x2,6000); %对信号做N点FFT变换figure(4)subplot(2,1,1),plot(k,abs(y2));title('加入噪声后|X(k)|');ylabel('幅度谱');subplot(2,1,2),plot(k,angle(y2));title('加入噪声后arg|X(k)|');ylabel('相位谱');(2)IIR低通滤波器%IIR低通滤波器clearclc[x1,fs,bits]=wavread('E:\数字信号.wav'); %读取语音信号的数据,赋给变量x1 x=x1(60001:1:120000); %截取原始信号60000个采样点d=randn(1,60000); %产生高斯白噪声mean(d);d=d/100; %对噪声进行衰减x2=x+d; %加入高斯白噪声sound(x2 ,fs); %播放加入高斯噪声后的语音,可以明显听出噪声fs=44100;Ts=1/fs;wp=2*pi*45000/fs; %通带截止频率ws=2*pi*50000/fs; %阻带截止频率Rp=1; %通带衰减Rs=100; %阻带衰减wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标ws1=2/Ts*tan(ws/2);[N,Wn]=buttord(wp1,ws1,Rp,Rs,'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))gridxlabel('频率/Hz')ylabel('频率响应幅度')title('IIR低通滤波器')f1=filter(bz,az,x2);figure(2)subplot(2,1,1)plot(x2) %画出滤波前的时域图title('IIR低通滤波器滤波前的时域波形');subplot(2,1,2)plot(f1); %画出滤波后的时域图title('IIR低通滤波器滤波后的时域波形');sound(f1,44100); %播放滤波后的信号F0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x2,1024);subplot(2,1,1);plot(f,abs(y2(1:512))); %画出滤波前的频谱图title('IIR低通滤波器滤波前的频谱')xlabel('频率/Hz');ylabel('幅值');subplot(2,1,2)F1=plot(f,abs(F0(1:512))); %画出滤波后的频谱图title('IIR低通滤波器滤波后的频谱')xlabel('频率/Hz');ylabel('幅值');(3)FIR低通滤波%FIR低通滤波clearclcfs=44100;x1=wavread('E:\数字信号.wav');x=x1(60001:1:120000); %截取原始信号60000个采样点d=randn(1,60000); %产生高斯白噪声d=d/100; %对噪声进行衰减x2=x+d; %加入高斯白噪声sound(x2 ,fs); %播放加入高斯噪声后的语音,可以明显听出噪声wp=2*pi*1000/fs;ws=2*pi*1200/fs;Rp=1;Rs=100;wdelta=ws-wp;N=ceil(8*pi/wdelta); %取整wn=(wp+ws)/2[b,a]=fir1(N,wn/pi,hamming(N+1)); %选择窗函数,并归一化截止频率figure(1)freqz(b,a,512);title('FIR低通滤波器');f2=filter(b,a,x2);figure(2)subplot(2,1,1)plot(x2)title('FIR低通滤波器滤波前的时域波形');subplot(2,1,2)plot(f2);title('FIR低通滤波器滤波后的时域波形');sound(f2,44100); %播放滤波后的语音信号F0=fft(f2,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x2,1024);subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR低通滤波器滤波前的频谱')xlabel('频率/Hz');ylabel('幅值');subplot(2,1,2)F2=plot(f,abs(F0(1:512)));title('FIR低通滤波器滤波后的频谱')xlabel('频率/Hz');ylabel('幅值');(4)IIR高通滤波器wp=2*2000/fs;ws=2*1000/fs;ap=1;as=30;[N,Wc]=buttord(wp,ws,ap,as);[B,A]=butter(N,Wc,'high');[H,W]=freqz(B,A);figure(1)plot(W,abs(H))title('高通滤波的幅值响应')f1=filter(B,A,x2);sound(f1,fs)figure(2)subplot(2,1,1)plot(x2)title('IIR高通滤波器滤波前的时域波形'); subplot(2,1,2)plot(f1)title('IIR高通滤波器滤波后的时域波形');f0=fft(f1,2048);f=fs*(0:1023)/2048;y=fft(x2,2048);y0=fft(f1,2048);figure(3)subplot(2,1,1)plot(f,abs(y(1:1024)))title('IIR高通滤波前的幅频图');subplot(2,1,2)plot(f,abs(y0(1:1024)));title('IIR高通滤波厚的幅频图');(5)FIR高通滤波器clc;[x1,fs,bits]=wavread('E:\数字信号.wav'); t=(0:length(x1)-1)/fs;Au=0.05;d=[Au*cos(2*pi*11025*t)];x2=x1'+d;T=1/fs;wp=2*pi*2000/fs;ws=2*pi*1000/fs;ap=1;as=50;B=wp-ws;N0=ceil(6.6*pi/B);N=N0+mod(N0+1,2);wc=(wp+ws)/2/pi;h=fir1(N-1,wc, 'high',hanning(N));t=1:N;figure(1)plot(t,h)y=fftfilt(h,x2);sound(y,fs)figure(2)subplot(2,1,1)plot(x2)title('FIR高通滤波器滤波前的时域波形'); subplot(2,1,2)plot(y)title('FIR高通滤波器滤波后的时域波形'); f0=fft(y,2048);f=fs*(0:1023)/2048;y0=fft(x2,2048);figure(3)subplot(2,1,1)plot(f,abs(y0(1:1024)))title('FIR高通滤波前的幅频图');subplot(2,1,2)plot(f,abs(f0(1:1024)));title('FIR高通滤波厚的幅频图');。