模拟信号数字化
一、课程设计的目的
通过本次课程设计的学习我们不仅要能够加深对在课堂中的理论知识的理解并且巩固理论课上所学的有关PCM编码和译码的基本的概念、理论和方法。利用Matlab软件平台运行，对模块进行调试、仿真和分析。对结果记性分析，输出结果应和结论一致，且符合理论。通过这些过程不但可以能锻炼我们发现问题、分析问题和解决问题的能力，而且还能够通过实验学习和掌握到书本中不能学习到的实践知识和经验。并且能够通过实验验证一些基本的简单的理论，为以后自己学习工作积累充分的实践知识和经验。
二、课程设计的要求
在本次课程设计中利用Matlab软件集成环境下的Simulink仿真平台通过理论分析，设计一个PCM编译码系统用示波器观察编码与译码前后的结果，即在示波器上信号波形，根据运行的数据和波形来分析该系统性能。结合理论分析验证结果的正确性。
（1）利用Simulink对系统建模设计一个PCM的编译解码器。
（2）输入模拟的话音信号并且观察其输出波形。
（3）对系统设计进行仿真演示。
三、课程设计涉及原理
（一）PCM原理
数字信号是对连续变化的模拟信号进行抽样、量化和编码产生的，称为PCM（Pulse-code modulation），即脉冲编码调制。PCM是实现语音信号数字化的一种方法，就是把一个时间连续，取值连续的模拟信号变换成时间离散，取值离散的数字信号后在信道中传输。脉冲编码调制就是对模拟信号先抽样，再对样值幅度量化，编码的过程。PCM系统的原理框图简要介绍。原理框图如下。
（二）PCM的工作包括：抽样，量化，编码。
1、抽样：把连续信号变为时间轴上离散的信号过程称为抽样抽样必须遵守奈奎斯特抽样定理。离散信号才可以完全代替连续信号。
2、量化：量化是把幅度连续变化的模拟量变成用有限为二进制数字量的过程。量化分为均匀量化和非均匀量化由于均匀量化的缺点所以实际中往往用非均匀量化。在实际中使用的是两种对数形式的压缩特性：A律和U律，U律PCM用于北美和日本，A律PCM用于欧洲和我国，由于A律压缩实现复杂，常使用13折线法编码。
A律13折线图
3、编码：编码就是把两户后的信号变换成代码，其反过程就成为译码。在现有的编码方法中，若安编码的速度来分，分为：低速编码和高速编码。在实际的PCM设备中，量化和编码是一起进行的。通信中采用高速编码方式。编码器的种类在大体上可归结为三类：逐次比较型、折叠级联型、混合型。
四、实验内容
任务一：PCM编码
A律PCM数字电话系统中规定：传输语音的信号频段为300~3400Hz。采样频率为8000次/s，对样值进行13折线压缩后编码为8位二数字序列。因此，PCM编码输出的数码速率为64kbps。
PCM编码的二进制序列中，每个采样值用8位二进制码表示，其中最高位比特位表示样值的正负性，规定负值用0表示，正值用1表示。接下来的3位比特表示样值的绝对值所在的8段折线的段落号，最后4位是样值处于段落内16个均匀间隔上的间隔序号。
在实验中用一个常数表示信号源，用Saturation将信号幅度限制在（-1,1）范围内，Relay模块的门限设置为0，其输出作为PCM编码输出的最高位，确定极性码。样值取绝对值后，以Look-Up Table模块进行13折线压缩，并用增益模块将样值范围放大到0~127，然后用间距为1的Quantier模块进行四舍五入取整量化，并用Inter to Bit Converter将整数转换成长度为8个比特的二进制数据，最后用Display模块显示编码结果。
13折线近似的PCM编码器测试模型和仿真结果
具体参数设置如下：
任务二：PCM解码
PCM解码子系统就是PCM编码部分。PCM解码器中首先分离并行数据中的最高为和7为数据位，然后将7位数据转换成整数值，再进行归一化、扩张后与双极性的极性码相乘得出解码值。
13折线A律解码器
具体参数设置如下：
任务三：PCM串行传输模型
在编码与解码的基础上，建立PCM串行传输模型，并在传输信道中加入指定错误概率的随机误码。仿真模型如下图，其中PCM编码和解码子系统见以上2个任务。PCM编码输出经过并串转换后得到二进制码流送入二进制对称信道。在解码端信道输出的码流经过串并转换后送入PCM编码，之后输出模型并显示波形。实验模型中信道传输速率最高为64kbps，仿真步进为1/64000s。信道错误比特率设为0.01，观察信道误码的位置，信号源采用200Hz 的正弦波。
PCM串行传输模型
具体参数设置如：
经过编码与解码之后，然后通过低通滤波器，最后在示波器Scope上得出输出波形，此外，还需注意在通过示波器前需经过To Frame和Buffer，其中To Frame是用来形成帧信号，Buffer 是缓冲器。
仿真波形
五、心得体会与总结
通过本次课程设计，我教系统地掌握了有关PCM脉冲编码调制的设计思路和设计方法，主要对MATLAB的仿真方法，开发环境等有了一定的了解并对其进行测试和应用。掌握了利用MATLAB集成环境下的Simulink仿真平台,完成对脉冲编码调制系统的建模与仿真分析的基本方法。以前对PCM编码的了解只是理论层面上的，通过这次课程设计，加深了对理论知识的理解，并对其在实际中的应用有了一定的了解。在课程设计刚开始的时候，由于对MATLAB的使用不太熟悉，觉得无从下手，课程设计进展困难。为此，我借助网络、图书等资源，一步一步熟悉MATLAB，并在熟练的同学的帮助下，渐渐摸索找到了方法，我认识到，在以后的学习过程中，要时刻保持着刻苦钻研的精神和坚持不懈的毅力。本此课程设计的成功离不开自己的努力，离不开同学的帮助，更离不开知道老师的教育和指导！xxx老师一方面在理论课上很详细很专业的为我们讲解了本次课程设计的理论知识，让我对此有了宏观上了解，并能够掌握这些理论知识，为以后的实际操作提供了坚实的基础。另一方面在实际操作时，也给我们很多技术上的指导，让我们能在此过程中，学到更多的操作技能。总的来说，这次的课程设计让我受益匪浅，不仅是我体验到了动手操作的乐趣，而且培养了我
的设计思维，增加了实践能力，让我体会到艰辛的同时，更让我收获了成功的喜悦。
抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于2倍信号最高频率。抽样频率小于2倍频谱最高频率时，信号的频谱有混叠。抽样频率大于2倍频谱最高频率时，信号的频谱无混叠。 语音信号处理是研究用数字信号处理技术和语音学知识对语音 信号进行处理的新兴学科，是目前发展最为迅速的学科之一，通过语音传递信息是人类最重要，最有效，最常用和最方便的交换信息手段，所以对其的研究更显得尤为重要。 Matlab语言是一种数据分析和处理功能十分强大的计算机应用 软件，它可以将声音文件变换成离散的数据文件，然后用起强大的矩阵运算能力处理数据。这为我们的本次设计提供了强大并良好的环境！ 本设计要求通过利用matlab对模拟信号和语音信号进行抽样，通过傅里叶变换转换到频域，观察波形并进行分析。 关键词：抽样Matlab
目录 一、设计目的： (2) 二、设计原理： (2) 1、抽样定理 (2) 2、MATLAB简介 (2) 3、语音信号 (3) 4、Stem函数绘图 (3) 三、设计内容： (4) 1、已知g1(t)=cos(6πt),g2(t)=cos(14πt),g3(t)=cos(26πt),以抽样频率 fsam=10Hz对上述三个信号进行抽样。在同一张图上画出g1(t)，g2(t)，g3(t)及其抽样点,对所得结果进行讨论。 (4) 2、选取三段不同的语音信号，并选取适合的同一抽样频率对其进 行抽样，画出抽样前后的图形，并进行比较，播放抽样前后的语音。 (6) 3、选取合适的点数，对抽样后的三段语音信号分别做DFT，画图 并比较。 (10) 四、总结 (12) 五、参考文献 (13)
《数字信号处理》课程设计报告 设计题目： IIR滤波器的设计 专业： 班级： 姓名： 学号： 指导教师： 2010年月日
1、设计目的 1、掌握IIR 滤波器的参数选择及设计方法； 2、掌握IIR 滤波器的应用方法及应用效果； 3、提高Matlab 下的程序设计能力及综合应用能力。 4、了解语音信号的特点。 2、设计任务 1、学习并掌握课程设计实验平台的使用，了解实验平台的程序设计方法； 2、录制并观察一段语音信号的波形及频谱，确定滤波器的技术指标； 3、根据指标设计一个IIR 滤波器，得到该滤波器的系统响应和差分方程，并根据差分方程将所设计的滤波器应用于实验平台，编写相关的Matlab 程序； 4、使用实验平台处理语音信号，记录结果并进行分析。 3、设计内容 3.1设计步骤 1、学习使用实验平台，参见附录1。 2、使用录音机录制一段语音，保存为wav 格式，录音参数为：采样频率8000Hz、16bit、单声道、PCM 编码，如图1 所示。 图1 录音格式设置 在实验平台上打开此录音文件，观察并记录其波形及频谱（可以选择一段较为稳定的语音波形进行记录）。 3、根据信号的频谱确定滤波器的参数：通带截止频率Fp、通带衰减Rp、阻带截止频率Fs、阻带衰减Rs。 4、根据技术指标使用matlab 设计IIR 滤波器，得到系统函数及差分方程，并记录得到系统函数及差分方程，并记录其幅频响应图形和相频响应图形。要求设计 第 1页出的滤波器的阶数小于7，如果不能达到要求，需要调整技术指标。 5、记录滤波器的幅频响应和系统函数。在matlab 中，系统函数的表示公式为：
因此，必须记录系数向量a 和b。系数向量a 和b 的可以在Matlab 的工作空间（WorkSpace）中查看。 6、根据滤波器的系统函数推导出滤波器的差分方程。 7、将设计的滤波器应用到实验平台上。根据设计的滤波器的差分方程在实验平台下编写信号处理程序。根据运行结果记录处理前后的幅频响应的变化情况，并试听处理前后声音的变化，将结果记录，写入设计报告。 3.2实验程序 （1）Rs=40; Fs=1400; Rp=0.7; Fp=450; fs=8000; Wp=2*pi*Fp;Ws=2*pi*Fs; [N,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); [b1,a1]=butter(N,Wn,'s'); [b,a]=bilinear(b1,a1,fs); [H,W]=freqz(b,a); figure; subplot(2,1,1);plot(W*fs/(2*pi),abs(H));grid on;title('频率响应'); xlabel('频率');ylabel('幅值');、 subplot(2,1,2); plot(W,angle(H));grid on;title('频率响应'); xlabel('相位(rad)');ylabel('相频特性'); 3.3实验结果（如图）： N =5 Wn=6.2987e+003 第 2页
语音信号的数字滤波 ——利用双线性变换法实现IIR数字滤波器的设计一．课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法，掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。 二．设计方案论证 1.IIR数字滤波器设计方法 IIR数字滤波器是一种离散时间系统，其系统函数为 假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数和，它 是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。 2.用双线性变换法设计IIR数字滤波器 脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用z=e sT转换 平面的-π/T～π到Z平面上。也就是说，第一步先将整个S平面压缩映射到S 1 /T一条横带里；第二步再通过标准变换关系z=e s1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图1所示。 图1双线性变换的映射关系 为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上，可以通过以下的正切变换实现
唐山学院 通信原理课程设计 题目模拟信号数字化传输系统的设计与仿真分析系 (部) 班级 姓名 学号 指导教师 2017 年 6 月 26 日至2017 年7月 8 日共 2 周
通信原理课程设计任务书 一、设计题目、内容及要求 设计题目：模拟信号数字化传输系统的设计与仿真分析 内容及要求： 1.了解Matlab/Simulink的运行环境及应用领域； 2.逐步熟悉模拟信号数字化传输系统的仿真过程，由简到难； 3.系统仿真及波形分析 (1) 模拟信号抽样过程原理与仿真分析； (2) 模拟信号量化过程原理与仿真分析； (3) PCM编译码系统设计与仿真分析； (4) DPCM编译码系统设计与仿真分析。 (5) 在高斯信道下对PCM系统的性能进行仿真分析。（可选） 二、设计原始资料 通信原理；软件Matlab；计算机一台 三、要求的设计成果（课程设计说明书、设计实物、图纸等） 设计说明书1份，不少于2000字，应包含模拟信号数字化传输系统原理、相关系统设计、相关软件Matlab/Simulink介绍、系统仿真及波形分析。 四、进程安排 第1-2天课设理论讲解及仿真软件介绍、学生练习使用软件 第3-4天相关系统设计 第5-6天系统仿真及波形分析 第7-8天整理、撰写说明书 第9-10天进行测试或答辩 五、主要参考资料 [1]樊昌信、曹丽娜.通信原理.北京:国防工业出版社，2006 [2]刘学勇.详解MATLAB/Simulink通信系统建模与仿真.北京:电子工业出版社，2011 [3]邵玉斌.MATLAB/Simulik通信系统建模与仿真实例分析.北京:清华大学出版社，2008 [4]张水英，徐伟强.通信原理及MATLAB/Simulink仿真.北京:人民邮电出版社，2012 [5]邵佳，董辰辉. MATLAB/Simulink通信系统建模与仿真实例精讲.北京:电子工业出版社，2009 指导教师（签名）：教研室主任（签名）：
信号与线性系统课程设计 报告 课题五基于FIR滤波的语音信号处理系统设计 班级: 姓名： 学号： 组号及同组人： 成绩： 指导教师： 日期：
课题五基于FIR滤波的语音信号处理系统设计 摘要：MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB 可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB特点：1) 高效的数值计算及符号计算功能,能使用户从繁杂的数学运算分析中解脱出来;2) 具有完备的图形处理功能,实现计算结果和编程的可视化;3)友好的用户界面及接近数学表达式的自然化语言,使学者易于学习和掌握;4)功能丰富的应用工具箱,为用户提供了大量方便实用的处理工具。 关键词：GUI界面，信号采集，内插恢复，重采样，滤波器 一、课程设计目的及意义 本设计课题主要研究数字语音信号的初步分析方法、FIR数字滤波器的设计及应用。通过完成本课题的设计，拟主要达到以下几个目的： （1）熟悉Matlab软件的特点和使用方法。 （2）熟悉LabVIEW虚拟仪器的特点以及采用LabVIEW进行仿真的方法。 （3）掌握信号和系统时域、频域特性分析方法。 （4）掌握FIR数字滤波器的设计方法（窗函数设计法、频率采样设计法）及应用。 （5）了解语音信号的特性及分析方法。 （6）通过本课题的设计，培养学生运用所学知识分析和解决实际问题的能力。 二、课题任务 （一）简单数字语音信号处理系统的Matlab设计。 使用GUI进行系统的图形用户界面设计，在该界面中包括对语音信号的读取，对信号的时域，频域分析，添加噪声，设计FIR数字滤波器(利用窗函数设计法、频率采样设计法任选)实现噪声滤除。具体任务如下： （1）对语音信号进行采集（读取），对数字语音信号加入干扰噪声，画出原始信号及带噪信号的时域波形，利用FFT进行频域分析，画出相应波形，并对语音进行播放。 （3）根据对语音信号及噪声的实际情况分析，选择适当的FIR数字滤波器进行设计，并对噪声进行滤除。
基于matlab的模拟信号数字化仿真 作者：李亚琼 学号：1305160425
摘要 本文研究的主要内容模拟信号数字化Matlab软件仿真。若信源输出的是模拟信号，如电话传送的话音信号，模拟摄像机输出的图像信号等，若使其在数字信道中传输，必须在发送端将模拟信号转换成数字信号，即进行A/D变换，在接收端则要进行D/A变换。模拟信号数字化由抽样、量化、编码三部分组成。由于数字信号的传送具有稳定性好，可靠性高，方便传送和传送等诸多优点，使得被广泛应用到各种技术中。不仅如此，Matlab仿真软件是常用的工具之一，可用于通信系统的设计和仿真。在科研教学方面发挥着重要的作用。Matlab有诸多优点，编程简单，操作容易、处理数据迅速等。 本文主要阐述的是模拟信号数字化的理论基础和实现方法。利用Matlab提供的可视化工具建立了数字化系统的仿真模型，详细讲述了抽样、量化、编码的设计，并指出了在仿真建模中要注意的问题。在给定的仿真条件下，运行了仿真程序，得到了预期的仿真结果。 关键词：Matlab、模拟信号数字化、仿真 1.1基本原理 模拟信号的数字传输是指把模拟信号先变换为数字信号后，再进行传输。由于与模拟传输相比，数字传输有着众多优点，因而此技术越来越受到重视。此变化成为A/D变换。A/D变换是把模拟基带信号变换喂数字基带信号，尽管后者的带宽会比前者大得很多，但本质上仍属于基带信号。这种传输可直接采用基带传输，或经过熟悉调制后再做频带传输。A/D变化包括抽样、量化、编码三个步骤，如图。 图1.模拟信号数字化 1.1.1抽样定理 抽样就是把模拟信号在时间上的连续变成离散的抽样值。而能不能用这一系列抽样值重新恢复原信号，就需要抽样定理来解决了。所以说，如果我们要传输模拟信号，可以通过传输抽样定理的抽样值来实现而不是非要传输原本的模拟信号。模拟信号数字化的理论基础就是抽样定理，抽样定理的作用不言而喻。 抽样定理：设时间连续信号) f，其最高截止频率为m f，如果用时间间 (t
数字信号处理课程设计 姓名：刘倩 学号：201014407 专业：信息与计算科学 实验一：常见离散信号产生和实现 一、实验目的： 1、加深对常用离散信号的理解； 2、掌握matlab 中一些基本函数的建立方法。 二、实验原理： 1.单位抽样序列 在MATLAB 中可以利用zeros()函数实现。 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： 2．单位阶越序列 在MATLAB 中可以利用ones()函数实现。 3．正弦序列 在MATLAB 中 4．复指数序列 在MATLAB 中 5．指数序列 在MATLAB 中
实验内容：由周期为10的正弦函数生成周期为20的余弦函数。 实验代码： n=0:30; y=sin(0.2*pi*n+pi/2); y1=sin(0.1*pi*n+pi/2); subplot(121) stem(n,y); xlabel ('时间序列n');ylabel('振幅');title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); subplot(122) stem(n,y1); xlabel ('时间序列n');ylabel('振幅'); title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); 实验结果： 实验二：离散系统的时域分析 实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。实验原理：离散系统 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号， 记系统单位冲激响应 则系统响应为如下的卷积计算式：
当N k d k ,...2,1,0==时，h[n]是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。 在MATLAB 中，可以用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应。 实验内容：用MATLAB 计算全解 当n>=0时，求用系数差分方程y[n]+y[n-1]-6y[n-2]=x[n]描述的一个离散时间系统对阶跃输入x[n]=8μ[n]的全解。 实验代码： n=0:7; >> [y,sf]=filter(1,[1 1 -6],8*ones(1,8),[-7 6]); >> y1(n+1)=-1.8*(-3).^n+4.8*(2).^n-2; >> subplot(121) >> stem(n,y); >> title('由fliter 函数计算结果'); >> subplot(122) >> stem(n,y1); >> title('准确结果'); 实验结果： 结果分析：有图可得由fliter 函数得出的结果与计算出的准确结果完全一致。 实验三FFT 算法的应用
数字信号处理课程设计报告《应用Matlab对信号进行频谱分析及滤波》 专业： 班级： 姓名： 指导老师： 二0 0五年一月一日
目录 设计过程步骤（） 2.1 语音信号的采集（） 2.2 语音信号的频谱分析（） 2.3 设计数字滤波器和画出其频谱响应（） 2.4 用滤波器对信号进行滤波（） 2.5滤波器分析后的语音信号的波形及频谱（） ●心得和经验（）
设计过程步骤 2.1 语音信号的采集 我们利用Windows下的录音机，录制了一段开枪发出的声音，时间在1 s内。接着在C盘保存为WAV格式，然后在Matlab软件平台下．利用函数wavread对语音信号进行采样，并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。通过wavread函数和sound的使用，我们完成了本次课程设计的第一步。其程序如下： [x,fs,bite]=wavread('c:\alsndmgr.wav',[1000 20000]); sound(x,fs,bite); 2.2 语音信号的频谱分析 首先我们画出语音信号的时域波形；然后对语音信号进行频谱分析，在Matlab中，我们利用函数fft对信号进行快速傅里叶变换，得到信号的频谱特性性。到此，我们完成了课程实际的第二部。 其程序如下： n=1024; subplot(2,1,1); y=plot(x(50:n/4)); grid on ; title('时域信号') X=fft(x,256); subplot(2,1,2); plot(abs(fft(X))); grid on ; title('频域信号'); 运行程序得到的图形：
一、语音信号去噪处理 1.设计要求： （1）在windows系统下的录音机录制一段1s左右的语音信号作为原声信号，在MATLAB软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数； （2）画出语音信号的时域波形，对采样后的语音进行fft变换，得到信号的频谱特性；对语音信号分别加入正弦噪声和白噪声，画出加噪信号的时域波形和频谱图； （3）根据对加噪语音信号谱分析结果，确定滤除噪声滤波器的技术指标，设计合适的数字滤波器，并画出滤波器的频域响应； （4）用所设计的滤波器对加噪的信号进行滤波，在同一个窗口画出滤波前后信号的时域图和频谱图，对滤波前后的信号进行对比，分析信号变化； （5）利用sound(x)回放语音信号，验证设计效果。 2.设计步骤： （1）找到7s的语音信号，利用函数wavread对语音信号进行信号读取；（2）计算样本时刻和频谱图的频率，并进行N+1点FFT变换； （3）加噪声为5000Hz的正弦信号正弦噪声，采用awgn函数加信噪比为10的高斯白噪声； （4）设计滤波器； （5）绘出相应的时域、频域图； （6）利用sound函数进行原始信号的语音播放，加噪声音播放，以及滤波之后的语言播放。 3.设计实现： （1）时域图与频谱图（加正弦） 录入原始信号的时域图： 加入正弦信号后的时域图：
滤波后的时域图： 录入原始信号的频域图： 加入正弦信号后的频率图： 滤波后的频域图： 采用巴斯低通滤波器滤除正弦波：
（2）具体代码实现： [x,fs,bits]=wavread('E:\mcpass.wav');%原信号 n=size(x,1); %提取采样信号的长度 t=(0:length(x)-1)/fs; %计算样本时刻 f=fs*(0:(n+1)/2-1)/n+1; %计算频域图的频率 X=fft(x,n+1); %进行N+1点FFT变换 ts=0:1/fs:(size(x)-1)/fs; %将所加噪声信号的点数调整到与原始信号相同 s=x+0.05*sin(2*pi*5000*ts)'; %加噪声为5000Hz的正弦信号正弦噪声 S=fft(s,n+1); %加正弦噪声后的频域 %正弦滤波 wp=2000/fs*2*pi; %2000为通带截止频率 ws=3000/fs*2*pi; %3000为阻带下限截止频率 Rp=4; %通带波纹 Rs=25; %阻带波纹 T=1/fs;Fs=1/T; %定义采样间隔 Wp=2/T*tan(wp/2); %计算对应的数字频率 Ws=2/T*tan(ws/2); [N,wn]=buttord(Wp,Ws,Rp,Rs,'s'); %计算滤波器介数和截止频率 [c,d]=butter(N,wn,'s'); %计算滤波器系统函数分子分母系数 [B,A]=bilinear(c,d,Fs); %双线性变换得到数字滤波器系统函数分子分母系数[Hb,Wc]=freqz(B,A); sf=filter(B,A,s); %对加噪信号进行滤波 Sf=fft(sf,n+1); %对滤波后进行N+1点FFT变换 %绘图部分 figure(3); plot(fs*Wc/(2*pi),20*log10(abs(Hb)));title('巴斯低通滤波器频域响应图'); xlabel('频率(Hz)'); ylabel('幅度'); figure(1); subplot(3,1,1); plot(t,x);title('原信号时域') xlabel('时间(s)'); ylabel('幅度'); figure(2); subplot(3,1,1); plot(f,abs(X(1:(n+1)/2)));title('原信号频域') xlabel('频率(Hz)'); ylabel('幅度'); figure(1); subplot(3,1,2); plot(t,s);title('加正弦信号后的时域') xlabel('时间(s)'); ylabel('幅度');
课程设计报告 课程设计名称：通信系统原理系：三系 学生姓名：张梦瑶 班级：11通信工程1班学号：20110306111 成绩： 指导教师：吴琼 开课时间：2013-2014学年一学期
一、课程设计目的 本课程是为通信工程专业本科生开设的专业必修课，结合学生的专业方向的理论课程，充分发挥学生的主动性，使学生掌握应用MATLAB或者SYSTEMVIEW等仿真软件建立通信系统，巩固理论课程内容，规范文档的建立，培养学生的创新能力，并能够运用其所学知识进行综合的设计。 通信系统原理的课程设计是对通信系统仿真软件、课程学习的综合检验，配合理论课的教学，让学生亲自参加设计、仿真、验证通信系统的一般原理、调制解调原理、信号传输及受噪声影响等方面的知识点。 二、设计选题 模拟信号的数字化处理 三、具体要求 a.熟悉模拟信号数字化的处理步骤：抽样、量化、编码； b.模拟信号的抽样过程，理解抽样频率的变化对抽样信号的影响； c.用MATLAB或其它EDA工具软件对PCM编码进行使用A律和μ律的压缩和扩张进行软件仿真； d.PCM的8位编码C1C2C3C4C5C6C7C8 e．仿真实现增量调制的过程和并理解噪声产生的原理。 四、进度安排 1、星期一查阅资料，确定选题和软件，思考总体设计方案； 2、星期二熟悉软件的编程环境； 3、星期三总体设计方案的确定与设计； 4、星期四各部分的具体实现（程序调试并程序注释）； 5、星期五整理完成设计报告的电子版，并答辩。 五、成绩评定办法 总成绩由平时成绩（考勤与课堂表现）、程序设计成绩和报告成绩三部分组成，各部分比例为30%,50%,20%。 （1）平时成绩：无故旷课一次，平时成绩减半；无故旷课两次平时成绩为0分，无故旷课三次总成绩为0分。迟到15分钟按旷课处理。
数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名： 学号： 组员：
摘要 滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR 滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。 关键词数字滤波器 MATLAB 窗函数法巴特沃斯
目录 摘要 (1) 1 引言 (1) 1.1课程设计目的 (1) 1.2 课程设计内容及要求 (1) 1.3课程设计设备及平台 (1) 1.3.1 数字滤波器的简介及发展 (1) 1.3.2 MATLAB软件简介 (2) 2 课程设计原理及流程 (4) 3.课程设计原理过程 (4) 3.1 语音信号的采集 (4) 3.2 语音信号的时频分析 (5) 3.3合成后语音加噪声处理 (7) 3.3.1 噪声信号的时频分析 (7) 3.3.2 混合信号的时频分析 (8) 3.4滤波器设计及消噪处理 (10) 3.4.1 设计IIR和FIR数字滤波器 (10) 3.4.2 合成后语音信号的消噪处理 (13) 3.4.3 比较滤波前后语音信号的波形及频谱 (13) 3.4.4回放语音信号 (15) 3.5结果分析 (15) 4 结束语 (15) 5 参考文献 (16)
￥ 课程设计报告? < 课程名称通信原理 设计题目 DSB与2ASK调制与解调 专业通信工程 班级 学号 姓名 完成日期 …
课程设计任务书 设计题目：DSB与2ASK调制与解调 设计内容与要求： 设计内容： 1.根据DSB的调制原理设计线路，进行仿真模拟调制DSB的调制和解调过程，并通过仿真软件观察信号以及的调制过程中信号波形和频谱的变化。 2. 根据ASK的调制原理设计线路，进行仿真模拟调制DSB的调制和解调过程，并通过仿真软件观察信号以及的调制过程中信号波形和频谱的变化。 3.在设计过程中分析信号变化的过程和思考仿真过程的设计原理。 ; 设计要求： 1.独立完成DSB与ASK的调制与解调； 2.运用仿真软件设计出DSB与ASK的调制线路 3.分析信号波形和频谱 指导教师：范文 2012年12月16日 课程设计评语 （ 成绩： 指导教师：_______________
年月日
一．调制原理： 调制: 将各种数字基带信号转换成适于信道传输的数字调制信号(已调信号或频带信号)； 时域定义：调制就是用基带信号去控制载波信号的某个或几个参量的变化，将信息荷载在其上形成已调信号传输，而解调是调制的反过程，通过具体的方法从已调信号的参量变化中将恢复原始的基带信号。 频域定义：调制就是将基带信号的频谱搬移到信道通带中或者其中的某个频段上的过程，而解调是将信道中来的频带信号恢复为基带信号的反过程. 根据所控制的信号参量的不同，调制可分为： 调幅，使载波的幅度随着调制信号的大小变化而变化的调制方式。 调频，使载波的瞬时频率随着调制信号的大小而变，而幅度保持不变的调制方式。 调相，利用原始信号控制载波信号的相位。 调制的目的是把要传输的模拟信号或数字信号变换成适合信道传输的信号，这就意味着把基带信号（信源）转变为一个相对基带频率而言频率非常高的代通信号。该信号称为已调信号，而基带信号称为调制信号。调制可以通过使高频载波随信号幅度的变化而改变载波的幅度、相位或者频率来实现。调制过程用于通信系统的发端。在接收端需将已调信号还原成要传输的原始信号，也就是将基带信号从载波中提取出来以便预定的接受者（信宿）处理和理解的过程。该过程称为解调。
（一）用窗函数法设计FIR数字滤波器 一、设计题目用窗函数法设计FIR数字低通滤波器 二、设计目的 1. 熟悉设计线性相位数字滤波器的一般步骤。 2. 掌握用窗函数法设计FIR数字滤波器的原理和方法。 3. 熟悉各种窗函数的作用以及各种窗函数对滤波特性的影响。 4. 学会根据指标要求选取合适的窗函数。 三、设计原理 窗函数法又称为傅里叶级数法，FIR数字滤波器的设计问题就是要所设计的FIR数字滤波器的响应H(ejw)去逼近所要求的理想滤波器的响应Hd(ejw)。从单位取样响应序列来看，就是使所设计的滤波器的h(n)逼近理想单位取样响应序列hd(n)。而且Hd(ejw)=逐段恒定的，且在频带边界处有不连续点，因此序列hd(n)是无限长的，通过直接截取无限长序列以得到有限长序列的办法，可以形象的比喻为h(n)通过一个窗口所看到的一段hd(n)。因此，h(n)也可以表达为hd(n)和一个窗函数w(n)的乘积，h(n)=w(n)hd(n)。这里的窗函数就是矩形序列RN(n)。 四、实现方法 用MATLAB编程实现给定指标要求的滤波器设计 五、设计内容及要求 1、各窗函数图（假设N=67;） N=67;
n=0:N-1; wn1=ones(1,N); stem(n, wn1);矩形窗 figure; wn2=hamming(N); stem(n, wn2);海明窗 figure; wn3=BARTLETT(N); stem(n, wn3);巴特列特 figure; wn4= Hanning(N); stem(n, wn4);汉宁窗 将窗函数分别画出来 2、计算理想低通滤波器单位冲激响应的源程序function[hd]=ideal(wc,N) q=(N-1)/2; n=0:N-1; m=n-q+eps; hd=sin(wc*m)./(pi*m); 3、计算频率响应的源程序 function[H]=fr(b,a,w); m=0:length(b)-1; l=0:length(a)-1; num=b*exp(-j*m'*w); den=a*exp(-j*l'*w); H=num./den;
实验一：基于DFT的数字谱分析以及可能出现的问题 一、实验目的： 1.进一步加深对DFT的基本性质的理解。 2.掌握在MATLAB环境下采用FFT函数编程实现DFT的语句用法。 3.学习用DFT进行谱分析的方法，了解DFT谱分析中出现的频谱泄露和栅栏效应现 象，以便在实际中正确应用DFT。 二、实验步骤: 1.复习DFT的定义、物理含义以及主要性质。 2.复习采用DFT进行谱分析可能出现的三个主要问题以及改善方案。 3.按实验内容要求，上机实验，编写程序。 4.通过观察分析实验结果，回答思考题，加深对DFT相关知识的理解。 三、上机实验内容: 1.编写程序产生下列信号供谱分析用： 离散信号： x1=R10(n) x2={1,2,3,4,4,3,2,1}，n=0,1,2,3,4,5,6,7 x3={4,3,2,1, 1,2,3,4}，n=0,1,2,3,4,5,6,7 连续信号： x4=sin(2πf1t)+sin(2πf2t) f1=100Hz, f2=120Hz，采样率fs=800Hz 2.对10点矩形信号x1分别进行10点、16点、64点和256点谱分析，要求256点 频谱画出连续幅度谱，10点、16点和64点频谱画出离散幅度谱，观察栅栏效应。 3.产生信号x2和x3分别进行8点、16点谱分析，画出离散幅度谱，观察两个信 号的时域关系和幅度谱的关系。 4.对双正弦信号x4以采样率fs=800Hz抽样，生成离散双正弦信号并画出连续波形； 对离散双正弦信号进行时域截断，截取样本数分别为1000、250、50。对不同样本的双正弦信号分别进行1024点谱分析，画出连续幅度谱，观察频谱泄露现象。
实验1 1、音乐信号的音谱和频谱观察 ○1使用wavread语句读取音乐信号，获取抽样率； ○2输出音乐信号的波形和频谱，观察现象； ○3使用sound语句播放音乐信号，注意不同抽样率下的音调变化，解释现象。 clear all;close all;clc [a,fs,bit]=wavread('c:\MATLAB6p5\work\陪你一起看草原.wav'); size(a); y1=a(:,1); a1=y1(10000:60000) figure; subplot(2,1,1),plot(a); subplot(2,1,2),plot(a1); x1=resample(a1,2,1); %y=resample(x,p,q)返回量的长度是向量x的p/q倍sound(x1,fs); %sound(a,fs); N1=length(a1); F1=fft(a1,N1); w=2/N1*[0:N1-1]; %频谱图横坐标设置 figure; plot(w,abs(F1)); N2=length(a1); t=0:1/N2:1/N2*(N2-1); title('傅利叶变换'); %傅利叶变换; figure; plot(a1); title('时域波形'); %时域波形;
1，以二倍的抽样率听声音信号时，音乐播放的特别快，像被压缩了，播放的时间比原信号短。 2，以二分之一的抽样率听声音信号时，音乐播放的特别慢，像被拉长了，播放的时间比原信号长。 3，原信号频谱截止频率为0.5*pi 实验2 2、音乐信号的抽取（减抽样） ○1观察音乐信号频率上限，选择适当的抽取间隔对信号进行减抽样（给出两种抽取间隔， 代表混叠和非混叠）； ○2输出减抽样音乐信号的波形和频谱，观察现象，给出理论解释； ○3播放减抽样音乐信号，注意抽样率的改变，比较不同抽取间隔下的声音，解释现象。 clear all;close all;clc [a,fs,bit]=wavread('c:\MATLAB6p5\work\陪你一起看草原.wav'); size(a); y1=a(:,1); a1=y1(10000:60000)
课程设计报告书 数字信号处理课程设计报告书基于Matlab对语音信号进行频谱分析及滤波 学院： 专业： 设计者： 学号： 设计周数： 完成日期：
目录 目录 (2) 一、课程设计目的 (2) 二、课程设计题目 (3) 三、课程设计题目描述和要求 (3) 五、课程设计流程 (5) 六、基于MATLAB的滤波器设计结果分析 (18) 八、设计经验总结 ...................................................................... 错误！未定义书签。 七、参考资料及网址 .................................................................. 错误！未定义书签。 八、附录 ...................................................................................... 错误！未定义书签。
[作者] 钟伟雄 [机构] 广东工业大学在校本科生 摘要： 数字信号处理是通信工程专业的一门相当重要的学科，对日后就业和科研有重大意义，通过MATLAB，我们可以清晰地理解数字信号处理中难以理解的一面，对理论的知识加以深化。 关键字：MATLAB 数字信号处理 GUI 频谱相位滤波器 一、课程设计题目 应用Matlab对语音信号进行频谱分析及滤波 二、课程设计目的 数字信号处理是一门以算法为核心，理论和实践性较强的学科。是电子信息工程、通信工程专业、电子信息科学与技术专业的一门重要的专业技术基础课。数字信号处理课程是在学习完数字信号处理的相关理论后，进行的综合性训练课程，其目的是： 1.使学生进一步巩固数字信号处理的基本概念、理论、分析方法和实现方法； 2.增强学生应用Matlab语言编写数字信号处理的应用程序及分析、解决实际问题的能 力； 三、课程设计内容描述和要求 为了巩固所学的数字信号处理理论知识，使学生对信号的采集、处理、传输、显示和存储等有一个系统的掌握和理解，安排了以下的课程设计的内容：录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；给定滤波器的性能指标，采用窗函数法和双线性变换设计滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集
计算机与信息工程学院验证性实验报告 一、实验目的 1．掌握信号的采样的方法。 2．分析信号的采样频率对频率特性的影响。 二．实验原理及方法 在现实世界里，声音、图像等各种信号多为模拟信号，要对它进行数字化处理，首先要将模拟信号经过采样、量化、编码，变成数字信号，即进行 A/D 转换，然后用数字技术进行数字信号处理，最后经过 D/A 转换成为模拟信号，这一处理过程称为模拟信号的数字信号处理． 在这一过程中最主要的是采样定理．采样定理是指对于一个c Ω≤Ω的带限信号，只要采样频率高于带限信号最高频率的两倍，即2c Ω>Ω时，则可以由其采样信号惟一正确地重建原始信号． 严格地说，在 MATLAB 中不能分析模拟信号，但当采样时间间隔充分小的时候，可以产生平滑的曲线，当时间足够长，可显示所有的模型，即近似的分析． 三．实验内容及步骤 已知1000()t a x t e -=，分别取采样频率为5000s f Hz = 和1000s f Hz =，绘出() X j ω 曲线,并对比两次结果说明采样频率对曲线的影响。
参考流程图： 四．实验报告要求 1．简述实验原理及目的． 2．写出程序并绘制图形，分析图形并指出频谱混迭的原因 五、实验程序 n=0:44;T=1; w=0.1*pi; x=sin(w*T*n); x1=fft(x) subplot(2,1,1),stem(n,x,'.'),hold on xlabel('nT'),ylabel('x') subplot(2,1,2),stem(n,x1,'.') title('离散时间变换') figure w=0.4*pi; x=sin(w*T*n); x1=fft(x) subplot(2,1,1),stem(n,x,'.'),hold on title('抽样后的图') xlabel('nT'),ylabel('x') subplot(2,1,2),stem(n,x1,'.') title('离散时间变换')
课程设计 课程设计名称：数字信号处理课程设计 专业班级：电信1203 学生姓名：刘海峰 学号： 201216020307 指导教师：乔丽红 课程设计时间：2015/07/01-2015/07/06 电子信息工程专业课程设计任务书
说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页
一. 技术要求 ?双线性变换法设计切比雪夫II型数字IIR低通滤波器， ?要求通带边界频率为400Hz， ?阻带边界频率分别为500Hz， ?通带最大衰减1dB， ?阻带最小衰减40dB， ?抽样频率为2000Hz， 二. 设计原理 IIR滤波器的设计包括三个步骤：①给出所需要的滤波器的技术指标； ②设计一个H(z)使其逼近所需要的技术指标：③实现所设计的H(z)，IIR数字滤波器设计的最通用的方法是借助于模拟滤波器的设计方法。所以IIR数字低通滤波器的设计步骤是：①按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标；②根据转换后的技术指标设计模拟低通滤波器G(s)：③再按一定规则将G(s)转换成H(z)。 在此过程中，我们用到了很多MATLAB中的函数，如设计切比雪夫低通滤波器的函数afd_chebl、由直接型转换为级联型的函数dir2cas、双线性变换的函数bilinear等。其中afd _chebl用于实现用模拟指标设计一个低通模拟滤波器，bilinear用于利用双线性变换法将模拟低通滤波器转换为数字低通滤波器。
三.程序流程图
四：源代码（完美版） %归一化低通滤波器技术指标 clc; clear all; Ap=1; %最大通带衰减 As=40; %最小阻带衰减 W=2000; %抽样周期 Wp=400; %通带边界频率 Ws=500; %阻带边界频率 wp=2*pi*Wp/W; %归一化通带边界频率 ws=2*pi*Ws/W; %归一化阻带边界频率 Wp1=tan(wp/2); %模拟低通滤波器通带边界频率 Ws1=tan(ws/2); %模拟低通滤波器阻带边界频率 %归一化切比雪夫II型低通模拟滤波器 [N,Wn]=cheb2ord(Wp1,Ws1,Ap,As,'s'); %确定滤波器阶数和频率尺度缩放因子 [BT,AT]=cheby2(N,As,Wn,'s');%传输函数的系数 [Z,P,K]=cheb2ap(N,As);%最小阻带衰减为As（DB）的N阶归一化模拟切比雪夫2型低通滤波器的零点、极点和增益因子 [H,W]=zp2tf(Z,P,K);%传输函数有理化形式 figure;
数字信号处理 课程设计报告 题目：数字调音台 学院：信息工程学院 专业：通信工程 指导教师： 指导时间:
目录 1背景与内容 (1) 2设计目的 (1) 3设计原理 (1) 3.1研究语音信号的特点； (1) 3.2IIR数字滤波器 (2) 3.3滤波器设计所涉及的matlab函数 (2) 3.4总体方案 (3) 4设计过程 (4) 4.1系统流程图 (4) 4.2具体设计 (4) 5实验代码及结果 (7) 5.1MATLAB源程序 (7) 5.2系统运行与测试 (19) 6实验体会及总结 (24) 7参考文献 (24)
1.背景与内容 根据数字信号处理课程设计教学大纲的要求，利用Matlab（工程设计软件）设计数字调音台，并且给出了具体的软件实施方案。Matlab功能强大、简单易学、编程效率高,深受广大科技工作者的欢迎。特别是Matlab还具有信号分析工具箱,不需具备很强的编程能力,就可以很方便地进行信号分析、处理和设计。此次设计利用Matlab对音频信号进行采样后；并画出音频信号的时域波形、频谱图；然后利用IIR设计高通、低通、带通以及帯阻滤波器对现有的音频信号进行滤波，实现音频音色的改变实现调音台的功能，并画出滤波后信号的时域波形和频谱，对滤波前后的信号进行对比，分析信号的变化；最后利用GUI设计一个数字信号系统界面。 1、选一首歌曲，采用已经学过的信号处理知识，设计高通、带通、低通、带阻几种滤波器，对其进行数字信号滤波，观察信号的变化； 2、实时显示处理前的信号频谱和处理后的信号频谱； 3、设计系统界面：为了使编制的程序操作方便，设计处理系统的用户界面，在所设计的系统界面上可以选择滤波器的类型，输入滤波器的参数、显示滤波器的频率响应等。 2.设计目的 1、通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉用IIR 数字滤波器和FIR滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法 2、掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。 3、通过对Matlab的使用，掌握Matlab的程序设计方法，掌握Matlab设计IIR数字滤波器和FIR滤波器的方法，并学会用Matlab对信号分析和处理。 4、学习使用MATLAB GUI设计工具应用平面。 3.设计原理 3.1研究语音信号的特点； 语音信号的带宽约为5KHz，主要能量集中在低频段。语音信号是典型的随机信号，人的每次发音过程都是一个随机过程。很难得到两次完全相同的发音样