实验三声音信号的采集
1.用Adobe Audition 3.0录制一段数字音频（英语朗诵）。
操作步骤：
（一）配置录音选项与麦克风测试：（如果已经能正常录音，可跳过不看这部分。）
首先，要明白“播放控制台（Volume Control）”与“录音控制台(Recording Control)”这两个概念，只有弄清楚它们的区别，才能正确的录音。
1.“播放控制台”：
使用“音量控制”可以调节计算机或其他多媒体应用程序（如CD 唱机、DVD 播放器和录音机）所播放声音的音量、平衡、低音、高音设置（高音与低音控制，通常是高级声卡才有此控制。）。也可以使用“音量控制”对话框调节系统声音、麦克风、CD 音频、线路输入、合成器和波形输出的级别。
打开“音量控制”对话框的方法：单击“开始”，指向“程序”，指向“附件”，指向“娱乐”，然后单击“音量控制”图标，就可以打开“音量控制”对话框，也就是我说的“播放控制台”。如图3.1所示：（上面操作也可鼠标指向任务栏右边的“喇叭”图标，单击鼠标右键，在打开的快捷菜单中选“打开音量控制”菜单）
图3.1 “音量控制”对话框
在上面这个图里，我们可以看到主音量、波形、软件合成器、麦克风音量、CD音量、线路音量等控制，它们的含义分别是：
（1）“主音量”：它能控制系统所有设备的播放音量，如果把它静音，整个系统将没有声音被放出来（ASIO、E-WDM 等专业驱动除外)。
（2）“波形”：它控制的是系统内部数字音频流的声音，如：WAV、MP3、WMA等格式的音量，而MIDI 播放的音量是不受它控制的。
（3）“软件合成器”，它专门控制MIDI 合成器的音量，如果静音，系统MIDI播放就无声，而其他WAV、MP3、WMA播放就不受此影响。
（4）“麦克风音量”：它是控制麦克风捕获声音并被送到声卡输出的音量，你的声卡MIC孔一旦插上麦克风并对麦克风说话，麦克风捕获的声音将被输出到声卡的输出系统，你就可以
直接在音箱或耳机里听到自己说话的声音，在录音过程中，这个通常是控制麦克风实时监听音量的，方便在录音的时候能从耳机里听到自己的声音，而不是单控制麦克风录音音量，无论你如何调节它，都不能改变麦克风录音音量大小（除非录音控制台里选择了“混音器”内录方式）。如果要控制麦克风录音音量，请参阅“录音控制台”相关说明，后面会讲到。
（5）“CD音量”：通常情况下，CD-ROM会有一条四芯模拟音频线（L、G、G、R）连接到声卡，这个项目就是控制CD-ROM驱动器模拟线路的音量的。现在，有很多播放器（Windows Media Player 9.0、Foobar2000等）都可以通过操作系统内部的DirectX 通道来播放CD音频，CD音频数据通道来自CD-ROM的IDE线路，从而实现“软数字播放”，如果是这样，CD Audio模拟音量将不控制CD播放音量。有些中高档声卡，还有一个专门的SPDIF CD 同轴接口，它是连接CD-ROM的两芯数字接口的，从而实现“硬数字播放”，硬数字播放的音质，和声卡的DAC（数模转换器）的品质有关。
（6）“线路音量”：所有声卡都有Line In接口，我们可以通过一根立体声音频线来把外部模拟设备（如：收音机、录音机、CD随身听、调音台等）连接到声卡，这样，我们就可以通过电脑来听这些设备播放的声音了。同Microphone控制一样，它仅仅是控制外部线路设备的播放音量，而不是控制录音音量。
2.“录音控制台”：
使用“录音控制”可以调节系统声音、麦克风、CD 音频、线路输入、合成器和波形输出的录音级别。我们通常使用它来控制这些设备的录音音量。
打开“录音控制台”的方法：先根据上面的方法打开“播放控制台”，在“播放控制台”控制条里的“选项(P)”菜单上，单击“属性”，在“混音器”下拉列表框中选“Realtek HD Audio Input”（提示：播放控制台的该项值为“Realtek HD Audio output”），单击“确定”按钮，在“显示下列音量控制”框里选择要调节其输入音量的设备（见图 3.2），然后单击“确定”。或者在“Adobe Audition 1.0”录音软件里，我们点“选项”菜单下的“录制调音台”，同样可以打开它。
图3.2 录音设备选择对话框
我们来看看“录音控制台”里面各个设备的含义：
（1）立体声混音：叫“立体声混音器录音音量控制”，一旦选择它，录音软件将录制系统发出的所有声音，即“播放控制台”里所有没有被静音的设备。我们通常叫系统内录。比如说，你使用播放器听MP3，选择“录音控制台”里的“Stereo Mix”，然后打开录音软件来录音，你就会发现录音软件会录下播放器所播放的MP3音乐，这个过程是内录的。通过这个选项，我们可以把那些网络限制下载的音乐直接录下来。
[注意：“播放控制台”里的麦克风如果没有静音，而“录音控制台”又刚好选择了“立体声混音”，这样录音软件同样会记录麦克风的声音，这样做，会对Adobe Audition（CoolEditPro）录音造成麻烦，因为人声音轨不但录下了你的人声，还同时录下了伴奏音轨的声音。对于多轨录音而言，这样做是不可取的，因为在后期处理前，你已经通过立体声混音器混合了人声与伴奏，就不能单独对人声或伴奏进行后期处理了，这个也是很多新手比较容易犯的错误。多轨录音是通过后期的“混缩”功能，来混合人声与伴奏的。]
（2）麦克风音量：对于一般的声音录音，这个才是真正的麦克风录音音量控制，我们必须选择它才对。只有这样，录音软件才会单独记录麦克风的声音，而不会记录其他声音。多轨录音中，人声音轨就不会录下伴奏音轨的声音，而是纯纯的人声了，这个才是我们想要的。
（3）CD 音量：选择它，你可以用录音软件来单独录制CD-ROM播放的音乐，不过，现在有强大的CD抓取软件，我们也可以不需要它了。
（4）线路音量：选择它，你可以使用录音软件来记录外部模拟设备的音乐，平时，把收音机连接到声卡的Line In 口上，这样就可以使用录音软件来录制电台的广播节目了。在专业的录音领域，录音师都使用调音台或麦克风放大器来连接专业麦克风的，而调音台或麦克风放大器就通过Line In 连接到声卡上。当然，高档的声卡与调音台，还可以通过数字光纤、数字同轴来相互连接。
（5）Front Mic：Front Mic是前置麦克风插孔，如果将这个选项打开，有时后会有杂音。（二）开始录音
（1）在“录音属性设置”对话框中选择“麦克风”，如图3.2所示，单击“确定”按钮，并调节好音量。
（2）连接好麦克风。
（3）启动Adobe Audition3.0:双击电脑桌面上的“”图标，进入audition 3.0工作窗口，选择“文件→新建”命令，出现新建波形对话框，选择默认值（采样率44100;通道立体声；分辨率 16位），单击“确定”按钮，如图3.3所示。
图3.3 新建波形对话框
（4）单击窗口左下角的“录音”按钮，开始“录音”，（注意：开始半分钟不要说话，录一段环境噪音约半分钟，以备后面降噪音时提供噪音样本），接着说话，录下你想录的语音内容，单击“停止”按钮，结束录音。单击窗口左下角的“播放”按钮，试听录音效果（如果操作不满意，执行“文件”菜单中的“关闭”命令，在出现“是否保存文件？”的对话框中，见图3.4，单击“否”按钮，不保存文件退出录音）
图3.4保存文件对话框
（5）降噪处理：
a.首先选择录制的噪音区，如图3-5所示
图3-5噪音选择区
b.选择“效果—修复—采集降噪预置噪声”菜单命令，单击“确定”按钮。
C.选择要去噪音的声波区，选择“效果—修复—降噪器”菜单命令，在“降噪器”对话框中单击“确定”按钮，如图3-6所示。
图3-6 降噪器对话框
D．剪去无用的波段：选中波段，按“Delete”键。
（6）保存录音文件：执行“文件”菜单中的“另存为(E)…”命令,在出现的“另存为”对话框中选择好文件保存路径、文件格式及文件名，单击“保存”按钮。如图3.7所示。执行“文件中“关闭”命令，结束录音操作。
图3.7保存文件对话框
2.用Adobe Audition
3.0从视频文件中提取声音。
（1）启动Adobe Audition 3.0,执行“选项”菜单中的“Windows 录音控制台（X）…”命令，在出现的“录音控制”属性对话窗口中，执行“选项”→“属性”命令，在“显示下列音量控制”列表框中勾选“立体声混音”，如图3.8所示。单击“确定”退出。
图3.8录音输入设备选择窗口
（2）执行“文件”→“打开视频中的音频文件（V）…”命令，选择一个视频文件，如cz.mpg，执行“文件”菜单中的“另存为(E)…”命令,在出现的“另存为”对话框中选择好文件保存路径、文件格式及文件名，单击“保存”按钮。如图3.9所示。执行“文件中“关闭”命令，结束操作
图3.9另存文件对话框
3.录制网上音乐（包括电脑中播放的音乐）
1）录音设备选择：
执行“选项/Windows录音控制台(X)”命令,选择“选项/属性”,在属性窗口中选择录音控制平台,勾选“麦克风音量”,见图3-10,单击“确定”返回。
图3-10录音设备选择对话框
2）执行“文件—新建…”,设置采样率、通道、分辨率
见图3-11,单击“确定”返回单轨视图窗口。
图3-11 新建波形对话
3）单击“传送器”面板中的录音按钮，同时播放电脑中的音乐文件或点播网上的音乐，这里开始录制音乐，按“传送器”面板中的“停止”按钮结束录音。
（4）选择“文件->另存为”命令（图3-9），选择好保存路径和输入文件名（sy3-3.wma），单击“保存”按钮保存录制的音乐文件。
基于MATLAB 的声音信号采集系统 野龙平 （陕西师范大学电子信息科学与技术，陕西） 摘要: 声音是各种信号传递与交流最直接的体现，因此对声音信号的研究有十分重要的意义。本文主要针对Matlab指令系统对声音信号的采集，作者利用Matlab 提供的数据采集工具箱, 介绍了倆种采集方法，简单分析并比较其优缺点。基于matlab的数据采集系统, 具有实现简单、性价比和灵活度高的优点。 关键词: Matlab; 数据采集 0 引言 随着科技的发展，对于语音信号的采集已经有很多种方法，如基于单片机技术、VC,C++等编程、纯硬件电路，本文介绍的方法主要通过一款软件MATLAB。它是MathWorks 公司推出的一种面向工程和科学运算的交互式计算软件, 其中包含了一套非常实用的工具-- 数据采集工具箱。使用此工具箱更容易将实验测量、数据分析和可视化的应用集合在一起。数据采集工具箱提供了一整套的命令和函数, 通过调用这些命令和函数, 可以直接控制数据采集设备的数据采集。 作者简单介绍了一种用声卡进行语音信号采集和MATLAB 的数据采集工具箱进行分析处理的语音信号采集系统。经实验证明, 该系统可实现在线连续采集语音信号并进行分析和处理, 具有实现简单、性价比和灵活度高的特点。 1 语音数据采集系统设计 MATLAB 中提供了强大的数据采集工具箱( DAQ- Data Acquisition Toolbox) , 可满足控制声卡进行数据采集的要求。用户通过调用MATLAB 命令, 可对采集的数据进行分析和处理, 为用户带来了极大的方便。 语音数据采集过程如图1 所示。
图1 声卡采集声音信号有两种方式: 传输线输入方式(LineIn) 和麦克风输入( MicIn) 方式。LineIn 方式是通过传输线把其他声音设备, 如录音机等设备的音频输出信号连接到声卡, 通过声卡记录数据存入计算机。 本系统采用MicIn 方式, 即用麦克风接收语音通过声卡将音频信号存入计算机。利用MicIn 方式通过声卡采集数据有两种方法: 方法一是采用对声卡产生一个模拟输入对象进行采集, 方法二是直接利用MATLAB 数据采集箱中提供的的函数命令进行采集。 1. 1 方法1 本系统是以声卡为对象利用MATLAB 数据采集工具箱提供的环境完成数据采集过程, 麦克风成为数据采集系统中的传感器。数据采集过程与其他硬件设备无关, 只与声卡有关, 因此应对声卡产生一个模拟输入对象(AI) 。 数据采集过程的具体实现: 1) 初始化: 创建设备对象。 ai= analoginput(‘ winsound ’ ) 2) 配置: 根据数据采集硬件设备的特性, 增加通道和控制数据采集的行为。为AI 添加1 个通道, 设置采样频率和采样时间。 addchannel( ai, 1) freq= 8000; % 采样频率fs8000Hz set( AI, sampleRate, freq) %为模拟输入设备设置采样频率 duration= 2; %采样时间为2 秒 set (AI, SamplesPerTriffer, duration* freq) ; % 为模拟输入设备设置触发时间 3) 执行: 启动设备对象, 采集数据。 start( ai) ; %启动设备对象 data= getdata( ai) ; % 获得采样数据 4) 终止: 删除设备对象。 stop( ai) ; % 停止设备对象 语音信号输入 声卡 Matlab 数据采集箱 计算机
语音信号处理实验指导书 实验一 语音信号采集与简单处理 一、 实验目的、要求 （1）掌握语音信号采集的方法 （2）掌握一种语音信号基音周期提取方法 （3）掌握短时过零率计算方法 （4）了解Matlab 的编程方法 二、 实验原理 基本概念： （a ）短时过零率： 短时内，信号跨越横轴的情况，对于连续信号，观察语音时域波形通过横轴的情况；对于离散信号，相邻的采样值具有不同的代数符号，也就是样点改变符号的次数。 对于语音信号，是宽带非平稳信号，应考察其短时平均过零率。 其中sgn[.]为符号函数 ?? ?? ?<=>=0 x(n)-1sgn(x(n))0 x(n)1sgn(x(n)) 短时平均过零的作用 1.区分清/浊音： 浊音平均过零率低，集中在低频端； 清音平均过零率高，集中在高频端。 2.从背景噪声中找出是否有语音，以及语音的起点。 （b ）基音周期 基音是发浊音时声带震动所引起的周期性，而基音周期是指声带震动频率的倒数。基音周期是语音信号的重要的参数之一，它描述语音激励源的一个重要特征，基音周期信息在多个领域有着广泛的应用，如语音识别、说话人识别、语音分析与综合以及低码率语音编码，发音系统疾病诊断、听觉残障者的语音指导等。因为汉语是一种有调语言，基音的变化模式称为声调，它携带着非常重要的具有辨意作用的信息，有区别意义的功能，所以，基音的提取和估计对汉语更是一个十分重要的问题。 ∑--= -=1 )]1(sgn[)](sgn[21N m n n n m x m x Z
由于人的声道的易变性及其声道持征的因人而异，而基音周期的范围又很宽，而同—个人在不同情态下发音的基音周期也不同，加之基音周期还受到单词发音音调的影响，因而基音周期的精确检测实际上是一件比较困难的事情。基音提取的主要困难反映在：①声门激励信号并不是一个完全周期的序列，在语音的头、尾部并不具有声带振动那样的周期性，有些清音和浊音的过渡帧是很难准确地判断是周期性还是非周期性的。②声道共振峰有时会严重影响激励信号的谐波结构，所以，从语音信号中直接取出仅和声带振动有关的激励信号的信息并不容 易。③语音信号本身是准周期性的(即音调是有变化的)，而且其波形的峰值点或过零点受共振峰的结构、噪声等的影响。④基音周期变化范围大，从老年男性的50Hz 到儿童和女性的450Hz ，接近三个倍频程，给基音检测带来了一定的困难。由于这些困难，所以迄今为止尚未找到一个完善的方法可以对于各类人群(包括男、女、儿童及不向语种)、各类应用领域和各种环境条件情况下都能获得满意的检测结果。 尽管基音检测有许多困难，但因为它的重要性，基音的检测提取一直是一个研究的课题，为此提出了各种各样的基音检测算法，如自相关函数(ACF)法、峰值提取算法(PPA)、平均幅度差函数(AMDF)法、并行处理技术、倒谱法、SIFT 、谱图法、小波法等等。 三、使用仪器、材料 微机（带声卡）、耳机，话筒。 四、 实验步骤 （1）语音信号的采集 利用Windows 语音采集工具采集语音信号，将数据保存wav 格式。 采集一组浊音信号和一组清音信号，信号的长度大于3s 。 （2）采用短时相关函数计算语音信号浊音基音周期，考虑窗长度对基音周期计算的影响。采用倒谱法求语音信号基音周期。 （3）计算短时过零率，清音和浊音的短时过零率有何区别。 五、实验过程原始记录（数据，图表，计算） 短时过零率 短时相关函数 P j j n s n s j R N j n n n n ,,1) ()()(1 =-=∑-= ∑--=-=10 )]1(sgn[)](sgn[21N m n n n m x m x Z
《matlab与信号系统》实验报告 学院： 学号： 姓名： 考核实验——语音信号采集与处理初步 一、课题要求 1．语音信号的采集 2．语音信号的频谱分析 3．设计数字滤波器和画出频率响应 4．用滤波器对信号进行滤波 5．比较滤波前后语音信号的波形及频谱 6．回放和存储语音信号 （第5、第6步我放到一起做了） 二、语音信号的采集 本段音频文件为胡夏演唱的“那些年”的前奏（采用Audition音频软件进行剪切，时长17秒）。运行matlab软件，在当前目录中打开原音频文件所在的位置，采用wavread函数对其进行采样，并用sound函数可进行试听，程序运行之后记下采样频率和采样点。 利用函数wavread对语音信号的采集的程序如下： clear; [y,fs,bits]=wavread('music.wav'); %x：语音数据；fs：采样频率；bits：采样点数sound(y,fs,bits); %话音回放 程序运行之后，在工作区间中可以看到采样频率fs=44100Hz，采样点bits=16
三、语音信号的频谱分析 先画出语音信号的时域波形，然后对语音号进行快速傅里叶变换，得到信号的频谱特性。语音信号的FFT频谱分析的完整程序如下： clear; [y,fs,bits]=wavread('music.wav'); %x：语音数据；fs：采样频率；bits：采样点数sound(y,fs,bits); %话音回放 n = length (y) ; %求出语音信号的长度 Y=fft(y,n); %傅里叶变换 subplot(2,1,1); plot(y); title('原始信号波形'); subplot(2,1,2); plot(abs(Y)); title('原始信号频谱'); 程序结果如下图： 四、设计数字滤波器和画出频率响应 根据语音信号的特点给出有关滤波器的性能指标： 1）低通滤波器性能指标，fp=1000Hz，fc=1200 Hz，As=100dB，Ap=1dB； 2）高通滤波器性能指标，fc=4800 Hz，fp=5000 Hz As=100dB，Ap=1dB。
虚拟仪器技术 姓名：史昌波 学号：2131391 指导教师：孙来军 院系(部所)：电子工程学院专业：控制工程
目录 1、前言 (3) 2、声卡的硬件结构和特性 (3) 2．1声卡的作用和特点 (3) 2．2声卡的构造 (5) 3、LABVIEW中与声卡相关的函数节点 (5) 4、LABVIEW程序设计 (6) 4．1程序原理 (6) 4．2程序结构 (7) 4．3结果分析 (9) 5、结束语 (9) 6、参考文献 (10)
基于声卡的数据采集与分析 1、前言 虚拟仪器技术是利用高性能的模块化硬件，结合高效灵活的软件来完成各种测试、测量和自动化的应用。在虚拟仪器系统中，硬件解决信号的输入和输出，软件可以方便地修改仪器系统的功能，以适应不同使用者的需要。其中硬件的核心是数据采集卡。目前市售的数据采集卡价格与性能基本成正比，一般比较昂贵1。 随着DSP(数字信号处理)技术走向成熟，计算机声卡可以成为一个优秀的数据采集系统，它同时具有A／D和D／A转换功能，不仅价格低廉，而且兼容性好、性能稳定、灵活通用，驱动程序升级方便，在实验室中，如果测量对象的频率在音频范围，而且对指标没有太高的要求，就可以考虑使用声卡取代常规的DAQ设备。而且LABVIEW中提供了专门用于声卡操作的函数节点，所以用声卡搭建数据采集系统是非常方便的2。 2、声卡的硬件结构和特性 2．1声卡的作用和特点 声卡的主要功能就是经过DSP（数字信号处理）音效芯片的处理，进行模拟音频信号的与数字信号的转换，在实际中，除了音频信号以外，很多信号都在音频范围内，比如机械量信号，某些载波信号等，当我们对这些信号进行采集时，使用声卡作为采集卡是一种很好的解决方案。 声卡的功能主要是录制与播放，编辑与合成处理，MIDI接口三个部分3。 （1）录制与播放
语音信号处理实验 班级： 学号： 姓名： 实验一基于MATLAB的语音信号时域特征分析（2学时）
1）短时能量 （1）加矩形窗 a=wavread('mike.wav'); a=a(:,1); subplot(6,1,1),plot(a); N=32; for i=2:6 h=linspace(1,1,2.^(i-2)*N);%形成一个矩形窗，长度为2.^(i-2)*N En=conv(h,a.*a);% 求短时能量函数En subplot(6,1,i),plot(En); if(i==2) ,legend('N=32'); elseif(i==3), legend('N=64'); elseif(i==4) ,legend('N=128'); elseif(i==5) ,legend('N=256'); elseif(i==6) ,legend('N=512'); end end
00.51 1.52 2.5 3 x 10 4 -1 1 x 10 4 024 x 10 4 05 x 10 4 0510 x 10 4 01020 x 10 4 02040 （2）加汉明窗 a=wavread('mike.wav'); a=a(:,1); subplot(6,1,1),plot(a); N=32; for i=2:6 h=hanning(2.^(i-2)*N);%形成一个汉明窗，长度为2.^(i-2)*N En=conv(h,a.*a);% 求短时能量函数En subplot(6,1,i),plot(En); if(i==2), legend('N=32'); elseif(i==3), legend('N=64'); elseif(i==4) ,legend('N=128');
河南科技大学 课程设计说明书课程名称微机应用技术课程设计 题目语音信号采集与处理课程设计 学院医学技术与工程学院 班级生物医学工程1201班 学生姓名 指导教师杨晓利 日期 2014年3月29日
课程设计任务书 （指导教师填写） 课程设计名称微机应用技术课程设计 学生姓名 专业班级生物医学工程1201班 设计题目语音信号采集与处理课程设计 课程设计目的 1．学会MATLAB的使用，掌握MATLAB的程序设计方法； 2．掌握在Windows环境下语音信号采集的方法； 3．掌握信号处理的基本概念、基本理论和基本方法； 4．学会用MATLAB对信号进行分析和处理。 设计内容、技术条件和要求 1．语音信号的采集 用windows自带的录音机（开始—程序—附件—娱乐—录音机，文件—属性—立即转换—8000KHz，8位，单声道）或其他软件，录制一段语音信号，时间控制在2秒左右。然后在MATLAB软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。通过wavread函数的使用，理解采样频率、采样位数等概念。 wavread函数调用格式： y=wavread(file)，读取file所规定的wav文件，返回采样值放在向量y中。 [y,fs,nbits]=wavread(file)，采样值放在向量y中，fs表示采样频率（Hz），nbits表示采样位数。 y=wavread(file,N)，读取前N点的采样值放在向量y中。 y=wavread(file,[N1,N2])，读取从N1点到N2点的采样值放在向量y中。2．语音信号的频谱分析 首先画出语音信号的时域波形；然后对语音信号进行频谱分析，在MATLAB中，可以利用函数fft对信号进行快速傅里叶变换，得到信号的频谱特性。 时间进度安排 第1周：查阅资料； 第2周：实现设计内容 第3周：整理资料，撰写课程设计任务书 1
通信工程学院12级1班罗恒2012101032 实验一语音信号的低通滤波和短时分析综合实验 一、实验要求 1、根据已有语音信号，设计一个低通滤波器，带宽为采样频率的四分之一，求输出信号； 2、辨别原始语音信号与滤波器输出信号有何区别，说明原因； 3、改变滤波器带宽，重复滤波实验，辨别语音信号的变化，说明原因； 4、利用矩形窗和汉明窗对语音信号进行短时傅立叶分析，绘制语谱图并估计基音周期，分析两种窗函数对基音估计的影响； 5、改变窗口长度，重复上一步，说明窗口长度对基音估计的影响。 二、实验目的 1.在理论学习的基础上，进一步地理解和掌握语音信号低通滤波的意义，低通滤波分析的基本方法。 2.进一步理解和掌握语音信号不同的窗函数傅里叶变化对基音估计的影响。 三、实验设备 1.PC机； 2.MATLAB软件环境； 四、实验内容 1.上机前用Matlab语言完成程序编写工作。 2.程序应具有加窗（分帧）、绘制曲线等功能。 3.上机实验时先调试程序，通过后进行信号处理。 4.对录入的语音数据进行处理，并显示运行结果。 5. 改变滤波带宽，辨别与原始信号的区别。 6.依据曲线对该语音段进行所需要的分析，并且作出结论。 7.改变窗的宽度（帧长），重复上面的分析内容。 五、实验原理及方法 利用双线性变换设计IIR滤波器（巴特沃斯数字低通滤波器的设计），首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得所要设计的IIR滤波器的系统函数H(z)。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里主要是边界频率Wp和Ws 的转换，对ap和as指标不作变化。边界频率的转换关系为∩=2/T tan(w/2)。接着，按照模拟低通滤波器的技术指标根据相应设计公式求出滤波器的阶数N和3dB截止频率∩c ；根据阶数N查巴特沃斯归一化低通滤波器参数表，得到归一化传输函数Ha(p)；最后，将p=s/ ∩c 代入Ha(p)去归一，得到实际的模拟滤波器传输函数Ha(s)。之后，通过双线性变换法转换公式s=2/T((1-1/z)/(1+1/z))得到所要设计的IIR滤波器的系统函数H(z)。
通信与信息工程学院 信息处理综合实验报告 班级：电子信息工程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的工作空间。
数字信号处理 课程设计 题目：基于MATLAB的语音信号的采集与处理学院：皖西学院 专业：通信工程 班级：通信1001班 学号：2010013461 2010013494 姓名：刘敏纵大庆指导教师：何富贵
摘要: 本次课程设计题目为<<基于MATLAB的语音信号的采集与处理>>。首先我们利用计算机上的录音软件获得语音信号，然后利用MATLAB对语音信号进行分析和处理，采集语音信号后，利用MATLAB软件平台进行频谱分析；并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号！
1.背景 2. 设计目的 (2) 3. 设计原理 (2) 4. 设计过程 .......................................... ,,, 3 5. 实验代码及结果 (4) 5.1 语音信号的采集 (4) 5.2 语音信号加噪与频谱分析 ..................................... ,,,, 7 5.3 巴特沃斯滤波器的设计 .. (9) 5.4 比较滤波前后语音信号波形及频谱 (10) 6. 收获与体会 (12) 参考文献 (13)
1. 引言 数字信号处理是利用计算机或专用处理设备，以数值计算的方法对信号进行采集、抽样、 变换、综合、估值与识别等加工处理，借以达到提取信息和便于应用的目的。它在语音、雷达、图像、系统控制、通信、航空航天、生物医学等众多领域都获得了极其广泛的应用。具有灵活、精确、抗干扰强、度快等优点。 数字滤波器，是数字信号处理中及其重要的一部分。随着信息时代和数字技术的发展，受到人们越来越多的重视。数字滤波器可以通过数值运算实现滤波，所以数字滤波器处理精 度高、稳定、体积小、重量轻、灵活不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。数字滤波器种类很多，根据其实现的网络结构或者其冲激响应函数的时域特性，可分为两种，即有限冲激响应(FIR, Finite Impulse Response)滤波器和无限冲激响应(IIR，Infin ite Impulse Resp on se) 滤波器。 FIR滤波器结构上主要是非递归结构，没有输出到输入的反馈，系统函数H (z)在 z= R处收敛，极点全部在z = 0处(因果系统)，因而只能用较高的阶数达到高的选 择性。FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低，但是线性相位，就是 不同频率分量的信号经过FIR滤波器后他们的时间差不变，这是很好的性质。FIR数字滤波器是有限单位脉冲响应有利于对数字信号的处理，便于编程，用于计算的时延也小，这对实时的信号处理很重要。FIR滤波器因具有系统稳定，易实现相位控制，允许 设计多通带(或多阻带)滤波器等优点收到人们的青睐[1]。 IIR滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由 延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。同时，IIR数字滤波器在设计上可以借助成熟的 模拟滤波器的成果，如巴特沃斯滤波器等。
摘要：科技的发展带来许多电磁干扰或射频干扰的恶劣环境，要想解决电磁干扰问题的，必须从本质上改变麦克风的工作模式。文章提出了利用激光的传输频带宽、通信容量大和抗电磁干扰能力强等优点，研制一种基于光相位干涉的高灵敏度声音采集系统。光纤声音采集能够传送非常高的声音质量，适用于多种恶劣环境下的声音采集。 关键词：光纤声音采集、干涉型光纤传感器 引言： 麦克风在声场和电场中起着重要的沟通界面，它可将声音信号传至任何地方或者记忆装置。传统型的使用电磁场或静电场来产生动作，外部的强电磁场影响会阻绝这些装置的功能。本项目研制的光纤声音采集系统是一种新颖的声音信号传感器，在反射式强度型光纤传感器的原理基础上，利用激光来采集声音信号，由于它与传统的麦克风有着本质的区别，所以在使用方面具有很大的优越性。系统由非导磁材料制成，其主要工作本体是光，即使在强电磁场或高射频环境中也能正常工作。把光纤应用于麦克风，充分利用了光纤传感器体积小、结构简单、灵敏度高、抗电磁干扰且光纤本身低损耗、耐腐蚀、安全可靠等优良特性。 1、系统结构 本系统利用干涉型光纤传感器的原理，开发基于光相位干涉的高灵敏度声音采集系统，由光纤传感探头、光路系统、光信号调制解调器等部分组成。 干涉型光纤传感器通常将被测量转化为光信号的相位，因此，相位测量是该类型传感器信号处理的基本要求。若直接对相位进行测量，那么有两个问题将限制系统的性能：一是系统受到环境的干扰时被测相位会产生随机漂移，从而引入测量误差，此外，相位漂移还会导致信号衰减；二是直接测相意味着直流检测，信号处理易受电路直流漂移的影响。针对这两个问题引出了相位生成载波技术。相位生成载波调制是在被测信号带宽以外的某一频带之外引入大幅度的相位调制，被测信号则位于调制信号的边带上，这样就把外界干扰的影响转化为对调制信号的影响，且把被测信号频带与低频干扰频带分开，以利于后续的噪声分离。 项目研制的光纤声音采集系统，在对传统michelson干涉仪加以改进的基础上，通过构造由光纤耦合器和振动膜组成的动态michelson干涉光路，能够将外界声压对振膜的作用转化为对光路相位的调制，得到的干涉光信号直接光电转换后即可解调还原声音信号。在多种干涉型光纤传感器的解调方法中，相位生成载波解调技术（pgc）由于是一种无源解调技术，并具有高灵敏度、大动态范围和好的线性度而得到广泛的应用。 2、系统原理 2.1光纤传感探头原理： 激光器发出的激光经耦合器到达传输光纤，由光纤出射的光束照射到振动膜上，传输光纤出射端面m1与振动膜构成一个干涉腔，从两表面反射回的光进行干涉，干涉光再经耦合器由光电探测器接收，外界声音信号通过改变干涉腔的光纤出射端面m1和振动膜之间的距离对光相位进行调制。系统中半导体激光器发出的光源光频随输入的调制电流线性变化，振动膜采用硅微技术进行研制。 2.2解调原理： 光纤声音采集系统中的调制解调器是由光源，光电转换器，高增益微弱信号放大电路，背景噪声消除器等组成。 光源向光纤传感头发射一稳定的激光，传感头内的振动薄膜被周围声音振动信号带动，从而对发射到振动薄膜上的激光进行相位调制后再反射回去，被调制的激光在光路系统里发生干涉，形成携带微弱声音强度的激光信号，光电转换系统的探测器将此激光信号转换成电信号，再经高增益微弱信号放大，pgc解调，噪声滤除，后将解调后的电信号还原成声音信号输出。
实验一语音信号的采集及预处理 一、实验目的 在理论学习的基础上，进一步地理解和掌握语音信号预处理及短时加窗的意义及基于matlab的实现方法。 二、实验原理 1.语音信号的录音、读入、放音等：练习matlab中几个音频处理函数，利用函数wavread 对语音信号进行采样,记住采样频率和采样点数，给出以下语音的波形图（2.wav）。利用wavplay或soundview放音。也可以利用wavrecord自己录制一段语音，并进行以上操作(需要话筒)。 2.语音信号的分帧：对语音信号进行分帧，可以利用voicebox工具箱中的函数enframe。 voicebox工具箱是基于GNU协议的自由软件，其中包含了很多语音信号相关的函数。3．语音信号的加窗：本步要求利用window函数设计窗口长度为256(N=256)的矩形窗(rectwin)、汉明窗(hamming)及汉宁窗(hann))，利用wvtool函数观察其时域波形图及频谱特性，比较得出结论。观察整个信号加矩形窗及汉明窗后的波形，利用subplot与reshape函数将分帧后波形、加矩形窗波形及加汉明窗波形画在一张图上比较。取出其中一帧，利用subplot与reshape函数将一帧语音的波形、加矩形窗波形及加汉明窗波形画在一张图上比较将得出结论。 4.预加重：即语音信号通过一个一阶高通滤波器1 9375 1- -z。 .0 三、实验步骤、实验程序、图形及结论 1.语音信号的录音、读入、放音等 程序： [x,fs,nbit]=wavread('D:\2.wav'); %fs=10000,nbit=16 y=soundview('D:\2.wav') 2.语音信号的分帧 程序： [x,fs,nbit]=wavread('D:\2.wav'); len=256; inc=128; y=enframe(x,len,inc); figure; subplot(2,1,1),plot(x) subplot(2,1,2),plot(y)
语音信号处理实验讲义 时间：2011-12
目录 实验一语音信号生成模型分析 (3) 实验二语音信号时域特征分析 (7) 实验三语音信号频域特征分析 (12) 实验四语音信号的同态处理和倒谱分析 (16)
实验一 语音信号生成模型分析 一、实验目的 1、了解语音信号的生成机理，了解由声门产生的激励函数、由声道产生的调制函数和由嘴唇产生的辐射函数。 2、编程实现声门激励波函数波形及频谱，与理论值进行比较。 3、编程实现已知语音信号的语谱图，区分浊音信号和清音信号在语谱图上的差别。 二、实验原理 语音生成系统包含三部分：由声门产生的激励函数()G z 、由声道产生的调制函数()V z 和由嘴唇产生的辐射函数()R z 。语音生成系统的传递函数由这三个函数级联而成，即 ()()()()H z G z V z R z = 1、激励模型 发浊音时，由于声门不断开启和关闭，产生间隙的脉冲。经仪器测试它类似于斜三角波的脉冲。也就是说，这时的激励波是一个以基音周期为周期的斜三角脉冲串。单个斜三角波的频谱表现出一个低通滤波器的特性。可以把它表示成z 变换的全极点形式 12 1()(1) cT G z e z --= -? 这里c 是一个常数，T 是脉冲持续时间。周期的三角波脉冲还得跟单位脉冲串的z 变换相乘： 112 1 ()()()1(1)v cT A U z E z G z z e z ---=?= ?--? 这就是整个激励模型，v A 是一个幅值因子。 2、声道模型 当声波通过声道时，受到声腔共振的影响，在某些频率附近形成谐振。反映在信号频谱图上，在谐振频率处其谱线包络产生峰值，把它称为共振峰。 一个二阶谐振器的传输函数可以写成 12 ()1i i i i A V z B z C z --= -- 实践表明，用前3个共振峰代表一个元音足够了。对于较复杂的辅音或鼻音共振峰要到5个以上。多个()i V z 叠加可以得到声道的共振峰模型 12 1 11 ()()11R r r M M i r i N k i i i i k k b z A V z V z B z C z a z -=---======---∑∑∑ ∑ 3、辐射模型 从声道模型输出的是速度波，而语音信号是声压波。二者倒比称为辐射阻抗，它表征了
工程设计论文 题目：基于MATLAB的语音信号采集与处理 姓名： 班级： 学号：
指导老师： 一.选题背景 1、实践意义： 语音信号是一种非平稳的时变信号，它携带着各种信息。在语音编码、语音合成、语音识别和语音增强等语音处理中无一例外需要提取语音中包含的各种信息。语音信号分析的目的就在于方便有效地提取并表示语音信号所携带的信息。所以理解并掌握语音信号的时域和频域特性是非常重要的。 通过语音相互传递信息是人类最重要的基本功能之一.语言是人类特有的功能.声音是人类常用工具,是相互传递信息的最重要的手段.虽然,人可以通过多种手段获得外界信息,但最重要,最精细的信息源只有语言,图像和文字三种.与用声音传递信息相比,显然用视觉和文字相互传递信息,其效果要差得多.这是因为语音中除包含实际发音容的话言信息外,还包括发音者是谁及喜怒哀乐等各种信息.所以,语音是人类最重要,最有效,最常用和最方便的交换信息的形式.另一方面,
语言和语音与人的智力活动密切相关,与文化和社会的进步紧密相连,它具有最大的信息容量和最高的智能水平。 语音信号处理是研究用数字信号处理技术对语音信号进行处理的一门学科,处理的目的是用于得到某些参数以便高效传输或存储;或者是用于某种应用,如人工合成出语音,辨识出讲话者,识别出讲话容,进行语音增强等. 语音信号处理是一门新兴的学科,同时又是综合性的多学科领域,是一门涉及面很广的交叉学科.虽然从事达一领域研究的人员主要来自信息处理及计算机等学科.但是它与语音学,语言学,声学,认知科学,生理学,心理学及数理统计等许多学科也有非常密切的联系. 语音信号处理是许多信息领域应用的核心技术之一,是目前发展最为迅速的信息科学研究领域中的一个.语音处理是目前极为活跃和热门的研究领域,其研究涉及一系列前沿科研课题,巳处于迅速发展之中;其研究成果具有重要的学术及应用价值. 数字信号处理是利用计算机或专用处理设备，以数值计算的方法对信号进行采集、抽样、变换、综合、估值与识别等加工处理，借以达到提取信息和便于应用的目的。它在语音、雷达、图像、系统控制、通信、航空航天、生物医学等众多领域都获得了极其广泛的应用。具有灵活、精确、抗干扰强、度快等优点。 数字滤波器, 是数字信号处理中及其重要的一部分。随着信息时代和数字技术的发展，受到人们越来越多的重视。数字滤波器可以通过数
基于LabView的双声道声卡数据采集系统班级：热动1007 姓名：刘堂俊学号：U201011568 在虚拟仪器系统中，信号的输入环节一般采用数据采集卡实现。商用的数据采集卡具有完整的数据采集电路和计算机借口电路，但一般比较昂贵，计算机自带声卡是一个优秀的数据采集系统，它具有A/D和D/A转换功能，不仅价格低廉，而且兼容性好、性能稳定、通用性强，软件特别是驱动程序升级方便。如被测对象的频率在音频范围内，同时对采样频率要求不是太高，则可考虑利用声卡构建一个数据采集系统。 1.从数据采集的角度看声卡 1.1声卡的作用 从数据采集的角度来看，声卡是一种音频范围内的数据采集卡，是计算机与外部的模拟量环境联系的重要途径。声卡的主要功能包括录制与播放、编辑和处理、MIDI接口三个部分。 1.2声卡的硬件结构 图1是一个声卡的硬件结构示意图。一般声卡有4~5个对外接口。 图1 声卡的硬件结构示意图 声卡一般有Line In 和Mic In 两个信号输入，其中Line In为双通道输入，Mic In仅作为单通道输入。后者可以接入较弱信号，幅值大约为0.02~0.2V。声音传感器（采用通用的麦克风）信号可通过这个插孔连接到声卡。若由Mic In 输入，由于有前置放大器，容易引入噪声且会导致信号过负荷，故推荐使用Line In ，其噪声干扰小且动态特性良好，可接入幅值约不超过1.5V的信号。另外，输出接口有2个，分别是Wave Out和SPK Out。Wave Out（或LineOut）给出的信号没有经过放大，需要外接功率放大器，例如可以接到有源音箱；SPK Out给出的信号是通过功率放大的信号，可以直接接到喇叭上。这些接口可以用来作为双通道信号发生器的输出。 1.3声卡的工作原理 声音的本质是一种波，表现为振幅、频率、相位等物理量的连续性变化。声卡作为语音信号与计算机的通用接口，其主要功能就是将所获取的模拟音频信号转换为数字信号，经过DSP音效芯片的处理，将该数字信号转换为模拟信号输出。输入时，麦克风或线路输入（Line In）获取的音频信号通过A/D转换器转换成数字信号，送到计算机进行播放、录音等各种处理；输出时，计算机通过总线将数字化的声音信号以PCM（脉冲编码调制）方式送到D/A转换器，变成模拟的音频信号，进而通过功率放大器或线路输出（Line Out）送到音箱等设备转换为声波。
通信工程学院12级1班 罗恒 2012101032 实验二 基于MATLAB 的语音信号频域特征分析 一、 实验要求 要求根据已有语音信号，自己设计程序，给出其倒谱、语谱图的分析结果，并根据频域分析方法检测所分析语音信号的基音周期或共振峰。 二、 实验目的 信号的傅立叶表示在信号的分析与处理中起着重要的作用。因为对于线性系统来说，可以很方便地确定其对正弦或复指数和的响应，所以傅立叶分析方法能完善地解决许多信号分析和处理问题。另外，傅立叶表示使信号的某些特性变得更明显，因此，它能更深入地说明信号的各项红物理现象。 由于语音信号是随着时间变化的，通常认为，语音是一个受准周期脉冲或随机噪声源激励的线性系统的输出。输出频谱是声道系统频率响应与激励源频谱的乘积。声道系统的频率响应及激励源都是随时间变化的，因此一般标准的傅立叶表示虽然适用于周期及平稳随机信号的表示，但不能直接用于语音信号。由于语音信号可以认为在短时间内，近似不变，因而可以采用短时分析法。 三、 实验设备 1.PC 机； 2.MATLAB 软件环境； 四、 实验内容 1.上机前用Matlab 语言完成程序编写工作。 2.程序应具有加窗（分帧）、绘制曲线等功能。 3.上机实验时先调试程序，通过后进行信号处理。 4.对录入的语音数据进行处理，并显示运行结果。 5.依次给出其倒谱、语谱图的分析结果。 6. 根据频域分析方法检测所分析语音信号的基音周期或共振峰。 五、 实验原理及方法 1、短时傅立叶变换 由于语音信号是短时平稳的随机信号，某一语音信号帧的短时傅立叶变换的定义为： 其中w(n -m)是实窗口函数序列，n 表示某一语音信号帧。令n -m=k'，则得到 ()()()jw jwm n m X e x m w n m e ∞-=-∞= -∑
数字信号处理 课程设计报告 课设题目：语音信号的采集与处理 学院： 专业： 班级： 姓名： 学号： 指导教师： 2011 年7月1日
课程设计报告撰写要求 1、页面设置 纸张大小设置为纵向A4，页边距设置为：上3.8厘米，下 3.5厘米，左3厘米，右3厘米，页眉设置为3厘米，页脚设置为2.7厘米，文档网络设置为指定行和字符网格，每行34字，每页34行。 2、段落及字体设置 除各级标题外，首行缩进2字符；图、表及图题、表题首行不缩进，居中放置；图表不应超出版心范围；行距采用单倍行距。 正文中文采用小四号宋体，英文采用新罗马字体（Times New Roman），段前0磅，断后0磅； 一级标题采用小二号黑体，段前12磅，段后12磅 二级标题采用小三号黑体，段前6磅，段后6磅 三级标题采用四号黑体，段前6磅，段后0磅 3、装订要求 采用左侧装订，订两钉。
目录 一. 课程设计任务 (1) 二. 课程设计原理及设计方案 (2) 三. 课程设计的步骤和结果 (6) 四. 课程设计总结 (39) 五. 设计体会 (40) 六. 参考文献 (41)
一. 课程设计任务 1、语音信号的采集 利用Windows下的录音机，录制一段自己的话音，时间在1s内，然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。 2、语音信号的频谱分析 在Matlab中，可以利用函数fft对信号进行快速傅立叶变换，得到信号的频谱特性，要求学生首先画出语音信号的时域波形，然后对语音信号进行频谱分析。 3、设计数字滤波器和画出其频率响应给出各滤波器的性能指标； 给定滤波器的性能指标如下： (1)低通滤波器的性能指标：fb=1000Hz,fc=1200Hz,As=100dB,Ap=1dB. (2)高通滤波器的性能指标：fc=4800Hz,fb=5000Hz,As=100dB,Ap=1dB. (3)带通滤波器的性能指标：fb1=1200Hz, fb2=3000Hz,fc1=1000Hz, fc2=3200Hz,As=100dB,Ap=1dB. 采用窗函数法和双线性变换法设计上面要求的3种滤波器，并画出滤波器的频率响应； 4、用滤波器对信号进行滤波 然后用自己设计的滤波器对采集到的信号进行滤波，画出滤波后信号的时域波形及频谱，并对滤波前后的信号进行对比，分析信号的变化； 5、回放语音信号，分析滤波前后的语音变化； 6、设计系统界面 为了使编制的程序操作方便，设计处理系统的用户界面，在所设计的系统界面上可以选择滤波器的类型，输入滤波器的参数、显示滤波器的频率响应，选择信号等。 - 1 -
音频信号采集及处理程序代码及实验结果图： [voice,fs]=audioread('notify.wav');%声音读取 sound(voice,fs); %声音回放 n=length(voice);%计算长度 voice1=fft(voice,n); %快速傅里叶变换 figure(1);subplot(2,1,1);plot(voice); %绘出时域波 xlabel('t');ylabel('amp');%坐标名称 title('初始音频信号时域波形');grid on; subplot(2,1,2);plot(abs(fftshift(voice1))); %绘出原始音频信号频谱 title('初始音频信号频域波形'); xlabel('f');ylabel('amp');grid on; t=0:1/fs:(n-1)/fs; noise=0.05*sin(2*pi*100000*t');%100kHz正弦波噪声 s=voice+noise;%加噪后的音频信号 pause;sound(s,fs); %播放加噪的语音 n=length(s); S=fft(s,n);%计算频谱 figure(2);subplot(2,1,1)；plot(s);%画出加噪之后的音频信号时域波 形 title('加噪声后的音频信号时域波形')； xlabel('t');ylabel('amp');grid on; subplot(2,1,2)；plot(abs(fftshift(S)));%零频移到频谱中心后，绘制加噪 之后的音频信号频谱 xlabel('f');ylabel('amp'); title('加噪声后的音频信号频域波形');grid on; pause; rp=2; rs=80; Ft=8000;Fp=1000;Fs=1300; wp=2*pi*Fp/Ft; ws=2*pi*Fs/Ft; %求出待设计的模拟滤波器的边界频率 [n,wn]=buttord(wp,ws,rp,rs,'s'); %低通滤波器的阶数和截止频率 [b,a]=butter(n,wn,'s'); %S域频率响应的参数即：滤波器的传输函数 [bz,az]=bilinear(b,a,0.5); %利用双线性变换实现频率响应S域到Z域的变换 [h,w]=freqz(bz,az); figure(3);plot(w*fs/(2*pi),abs(h));%绘制IIR低通滤波器特性曲线 title('IIR低通滤波器特性曲线');grid on; z=filter(bz,az,s); %滤波 pause;sound(z,fs); %回放滤波后的信号 Z=fft(z); %滤波后的信号频谱 figure(4);subplot(2,2,2);plot(z);%绘制低通滤波后的音频信号时域