信号峰值检测仪课程设计报告.
东北石油大学课程设计
2017年7月10日
目录
第1章概述 (1)
1.1 引言 (1)
第2章总体设计方案 (3)
2.1 系统设计方案................................................. 错误！未定义书签。
2.2 系统设计框图 (4)
第3章单元电路的设计 (5)
3.1 总电路的设计 (5)
3.2 待测电信号放大电路模块的设计 (6)
3.3 信号峰值检测电路模块的设计 (8)
3.4 NE555芯片25khz时钟信号产生电路模块的设计 (9)
3.5模数转换模块的设计 (11)
3.6 锁存器锁存数字信号连接数码管显示模块的设计 (12)
第4章电路的仿真 (14)
结论 (17)
参考文献 (18)
附录元器件清单 (19)
第1章概述
1.1 引言
峰值检测是电子测量自动化仪表以及其它相关技术领域常会遇到的问题。峰值反映了信号极为重要的方面，尤其是小信号。设计完善的峰值检测系统，不仅可以用于对微弱信号进行检测，还可以通过传感器对其他非电信号如微弱的机械振动实现自动检测和控制，从而构成完整的检测系统，因此峰值检测具有广泛的使用价值。
峰值检测通过对输入信号的峰值进行提取，产生输出，从而实现对信号的监控，保证系统中其他结构工作的稳定性，避免了过大输入对系统造成的损伤，延长了设备的使用寿命。峰值检测是电子测量、自动化仪表以及其它相关技术领域常会遇到的问题。峰值反映了信号极为重要的方面。尤其是小信号。设计完善的峰值检测系统，不仅可以用于对微弱信号进行检测，还可以通过传感器对其它非电信号如微弱的机械振动实现自动检测和控制，从而构成完整的测控系统，因此峰值检测具有广泛的实用价值。
峰值检测技术是数字存储示波器及数据采集卡中的重要技术之一，用来实现波形的峰值捕捉。在科研、生产的许多领域都需要用到峰值检测设备，比如检测某建筑物中梁的最大承受力，检测一根钢丝绳的最大允许拉力等，这就需要用到相应的检测设备。目前常用的方法是先求得检测信号的平均值，但使用平均值掩盖了被检测信号的突然脉冲，从而可能引起系统的失灵及不稳定。若用由二极管和电阻电容构成的普通峰值检波电路来检波，效果会很差，主要表现在两个方面：第一，若选择RC电路时间常数大一些，则输出信号的波形会好一些。但检波输出之后的信号幅值和检波之前的信号幅值有明显的差距,输出信号幅值明显降低，峰值检波效率变差,同时,信号快变部分的丢失变得严重。第二，若选择RC电路时间常数小一些。则会发现检波前后的信号幅值的差异变小。信号之中的快变分量明显变好，但输出信号的波形明显变差。不利于对信号的A/D变换。
峰值检测电路是将某一段时间内信号的最值反映出来，遇到信号峰值就跟着，若没有更大的就保持，其工作状态主要包括跟踪、保持和复位，可以广泛应用于信号采集和处理、仪器仪表、自动控制等众多领域。如用来实现波形的毛刺捕捉、冲刺信号峰值检测等。因此峰值检测器有广泛的市场应用，且随着电子产
业和信息产业的不断发展，对信号的采集处理方面的要求愈加突出，它将发挥愈来愈大的作用。
第2章总体设计方案
2.1 系统设计方案
1.基本原理分析
峰值检测器是用来检测交流电压峰值的电路。最简单的峰值检测器依据半波整流原理构成电路，如下图2-1。交流电源在正半周的一段时间内，通过二极管对电容充电，使电容上的电压逐渐趋近于峰值电压。只要RC够大，可以认为其输出的直流电压数值上十分接近于交流电压的峰值。交流电源在正半周的一段时间内，通过二极管对电容充电，使电容上的电压逐渐趋近于峰值电压。这种简单电路的工作过程是，在交流电压的每一周期中，可分为电容充电和放电两个过程。在交流电压的作用下，在正半周的峰值附近一段时间内，通过二极管对电容C 充电，而在其它时段电容C上的电压将对电阻R放电。
图2-1 简单峰值检测电路
基于该简单的峰值检测器，如采用集成电路对其功能实现，性能进行优化，则更具有实用性。峰值检测器之原理可概括为：当输入信号为峰值更大值时，电容刷新进行充电，相关原件记录新的峰值，如若未有新峰值，相关元件则进行保持功能。
2.设计思路
本课设的关键任务是检测出峰值并使之保持稳定和数字显示。将输入量的峰值检测出来并保持，在通过数据锁存控制电路所存峰值的数字量。通过A/D数模转换，将模拟量转换为数字量，然后需要将数字量显示在数码管上。该方案由待测电信号放大电路、信号峰值检测电路、NE555芯片25KHZ时钟信号产生电路、模数转换电路、锁存器锁存数字信号连接数码管显示电路五部分组成。
2.2 系统设计框图
图2-2信号峰值检测系统框图
第3章单元电路的设计
3.1 总电路的设计
将待测信号源经LM358构成的运算放大电路放大输出后，接一单刀双掷开关，开关一端接ADC16模数转换芯片转化为16位二进制数字信号，开关另一端接峰值信号检测电路并且输出信号峰值再接入ADC16模数转换芯片转换成16位数字信号，拨动开关可选择测量实时信号和峰值信号，NE555构成的多协振荡器输出25KHZ方波为ADC16提供时钟信号以保证ADC16能正常工作输出相应的数字信号，由于模数转换的过程是需要一定时间的，转换过程中会出现乱码所以在16输出端接入两个74LS273锁存器当ADC16转换完成后ADC16的SOC 引脚会输出低电平再接入7406非门就会产生一个上升沿驱动74LS273锁存并输出转换完成的16位数字量，最后将16位数据通过16跟数据线接在4个16进制数码管上显示出相应的值，完成信号实时检测和峰值检测。仿真电路有改进的地方也有不完善的地方，仿真电路用到了16位模数转换芯片其分辨率高达0.15MV，采用了运放电路将待测信号放大5倍后进行测量也提高了测量的精度，再有就是单刀双掷开关的设计可以选择性的测量信号峰值和实时信号，仿真也解决了采样保持的问题NE555构成的多谐振荡器发出25KHZ的时钟信号对于测量50HZ的电压信号其精度也非常高。仿真也有部分功能不够完善的部分，数码管采用的是4位16进制输出数字没有十进制看起来直观。仿真使用Multism13.0仿真软件进行，ADC模数转换芯片的转换范围是-5V~+5V当信号为5V时数码管显示FFFF为-5V时显示0000。测量结果为[0.00015*(数码管显示数字)-5]/5V。图3-1是由仿真软件Multism13.0绘制的应用ADC16设计的峰值信号检测仪，ADC16的引脚图如图3-2。
图3-1 由ADC16构成的信号峰值检测仪
图3-2 ADC16的芯片引脚示意图
3.2 待测电信号放大电路模块的设计
该模块如图3-3所示用LM358组成的非反向放大器以实现对待测信号的放大和输出，LM358 内部包括有两个独立的、高增益、内部频率补偿的双运算放大器，适合于电源电压范围很宽的单电源使用，也适用于双电源工作模式，在推荐的工作条件下，电源电流与电源电压无关。它的使用范围包括传感放大器、直流增益模块和其他所有可用单电源供电的使用运算放大器的场合。LM358正负极分别用+12V和-12V供能，放大倍数为5*（R5/R4+1），当输入端接收到待测
信号时能够稳定的进行放大并输出，当待测电信号在0~1v内变化时，输出电压的范围为0-5V，用以供给后续电路部分模数转换器正常工作，符合ADC输入电压范围。如图3-4所示用LM358组成的非反向放大器放大倍数的验证。图3-5所示为验证结果。
V1
图3-3 LM358组成的非反向放大器
图3-4所示用LM358组成的非反向放大器放大倍数的验证
3.3信号峰值检测电路模块的设计
该模块用于检测并稳定输出信号的峰值，如图3-6所示峰值信号检测器是由两级运算放大器组成的。第一级运放U4将输入信号的峰值传递到电容C 上，并保持下来。第二级运放U7组成缓冲放大器，将输出与电容隔离开来。为了获得优良的保持性能和传输性能，运放应具有输入阻抗高、响应速度快和跟随精度好等优点。
当输入电压V1上升时，V01跟随上升使二极管D2、D1导通，电容C 充电，VC 上升。当输入电压V1下降时，V01跟随下降，D2、D1截止，VC 值保持不变。当V1再次上升使V01上升并使D2、D1导通，再次对电容C 充电（VC 高于前次充电时电压）；当V1再次下降时，D2、D1又截止，VC 将峰值再次保持。输出VO 反映VC 的大小，并可通过AD 转换电路和译码显示峰值。峰值检测电路的Multism13.0仿真电路如图3-7所示，仿真波形如图3-8所示。
R7
VO
图3-6峰值信号检测电路
V2
1Vpk
1kHz
0° XSC1
图3-7 峰值检测电路的Multism13.0仿真电路
图3-8 峰值检测波形图
示波器A端接入一个正弦波，B端接的是输出波形，从仿真图中能直观看出此电路能实现峰值检测功能。
3.4 NE555芯片25KHZ时钟信号产生电路模块的设计
该模块电路是由NE555为核心构成的多谐振荡器，输出25KHZ的方波为AD 模块提供时钟信号，驱动ADC16正常转换。NE555是属于555系列的计时IC的其中的一种型号，555系列IC的接脚功能及运用都是相容的，只是型号不同的因其价格不同其稳定度、省电、可产生的振荡频率也不大相同；而555是一个用途很广且相当普遍的计时IC，只需少数的电阻和电容，便可产生数位电路所需的各种不同频率的脉波讯号。只需简单的电阻器、电容器，即可完成特定的振荡延时作用。其延时范围极广，可由几微秒至几小时之久。如图3-9所示为NE555构成的多谐振荡器，当R1=1.154KΩ，R2=2.309 KΩ时该模块电路便能输出25KHZ的方波。图3-10所示为方波波形图。
图3-9 NE555构成的多谐振荡器
图3-10 多谐振荡器输出的25KHZ 方波波形
3.5模数转换模块的设计
该模块是整个电路的核心部分，用到了ADC16模数转换芯片，将采集到的模拟信号转换成数字信号并输出。ADC16是8位A/D 转换器，所以其分辨率高达
图3-11 信号为5V时的检测仿真图
信号与系统课程设计报告材料
课程设计报告 课程名称信号与系统课程设计指导教师 设计起止日期 学院信息与通信工程 专业电子信息工程 学生 班级/学号 成绩 指导老师签字
目录 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);
信号分析课程设计报告书
信号分析课程设计 信号系统的时域分析 编程实现的卷积积分或卷积和 一、课程设计题目： 基于 MATLAB 的连续时间LTI 系统的时域分析 二、基本要求： ① 掌握连续时不变信号处理的基本概念、基本理论和基本方法； ② 学会 MATLAB 的使用，掌握 MATLAB 的程序设计方法； ③ 学会用 MATLAB 对信号进行分析和处理； ④ 编程实现卷积积分或卷积和，零输入响应，零状态响应； ⑤ 撰写课程设计论文，用信号处理基本理论分析结果。 三、设计方法与步骤： 一般的连续时间系统分析有以下几个步骤: ①求解系统的零输入响应; ②求解系统的零状态响应; ③求解系统的全响应; ④分析系统的卷积；⑤画出它们的图形. 下面以具体的微分方程为例说明利用MATLAB 软件分析系统的具体方法. 1．连续时间系统的零输入响应 描述n 阶线性时不变（LTI ）连续系统的微分方程为： 已知y 及各阶导数的初始值为y(0),y (1)(0)，… y (n-1)(0)， 求系统的零输入响应。 建模 当LIT 系统的输入为零时，其零输入响应为微分方程的其次解（即令微分方程的等号右端为零），其形式为（设特征根均为单根） 其中p 1,p 2,…,p n 是特征方程a 1λ n +a 2λn-1+…+a n λ+a n =0的根，它们可以 用root(a)语句求得。各系数 由y 及其各阶导数的初始值来确定。对此有 1121111n n m n n m m n n m d y d y dy d u du a a a a y b b b u dt dt dt dt dt -++-++?????++=+????++1212()n p t p t p t n y t C e C e C e =++????+120n C C C y ++????+=11220 n n p C p C p C Dy ++????+=
信号与系统课程设计报告
信号与系统课程设计报告 实验题目：信号的运算与处理 内容简介： 设计一个信号，对其进行信号运算和处理，利用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软件！（附上调试过程图片） 左边为文件、历史窗口，底下是命令窗口，最右下角为实验仿真波形，中间为运算程序，绘图画图程序。
数字信号处理课程设计报告 杨俊
课程设计报告 课程名称数字信号处理 课题名称数字滤波器设计及在语音信号分析中的应用 专业通信工程 班级1281 学号201213120101 姓名杨俊 指导教师彭祯韩宁 2014年12月5日
湖南工程学院 课程设计任务书 课程名称数字信号处理 课题数字滤波器设计 及在语音信号分析中的应用专业班级通信工程1281班 学生姓名杨俊 学号201213120101 指导老师彭祯韩宁 审批 任务书下达日期2014 年12月5日 任务完成日期2014 年12月13日
《数字信号处理》课程设计任务书 一、课程设计的性质与目的 《数字信号处理》课程是通信专业的一门重要专业基础课，是信息的数字化处理、存储和应用的基础。通过该课程的课程设计实践，使学生对信号与信息的采集、处理、传输、显示、存储、分析和应用等有一个系统的掌握和理解；巩固和运用在《数字信号处理》课程中所学的理论知识和实验技能，掌握数字信号处理的基础理论和处理方法，提高分析和解决信号与信息处理相关问题的能力，为以后的工作和学习打下基础。 数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应（IIR）滤波器和有限冲激响应（FIR）滤波器。 二、课程设计题目 题目1：数字滤波器设计及在语音信号分析中的应用。 1、设计步骤： （1）语音信号采集 录制一段课程设计学生的语音信号并保存为文件，要求长度不小于10秒，并对录制的信号进行采样；录制时可以使用Windows自带的录音机，或者使用其它专业的录音软件，录制时需要配备录音硬件（如麦克风），为便于比较，需要在安静、干扰小的环境下录音。 然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。 （2）语音信号分析 使用MATLAB绘出采样后的语音信号的时域波形和频谱图。根据频谱图求出其带宽，并说明语音信号的采样频率不能低于多少赫兹。 （3）含噪语音信号合成 在MATLAB软件平台下，给原始的语音信号叠加上噪声，噪声类型分为如下几种：①白
数字信号处理课程设计报告
《数字信号处理》课程设计报告 设计题目： 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页
信号发生器课程设计报告
目录 一、课题名称 (2) 二、内容摘要 (2) 三、设计目的 (2) 四、设计内容及要求 (2) 五、系统方案设计 (3) 六、电路设计及原理分析 (4) 七、电路仿真结果 (7) 八、硬件设计及焊接测试 (8) 九、故障的原因分析及解决方案 (11) 十、课程设计总结及心得体会 (12)
一、课题名称：函数信号发生器的设计 二、内容摘要： 函数信号发生器作为一种常用的信号源,是现代测试领域内应用最为广泛的通用仪器之一。在研制、生产、测试和维修各种电子元件、部件以及整机设备时，都要有信号源，由它产生不同频率不同波形的电压、电流信号并加到被测器件或设备上，用其他仪器观察、测量被测仪器的输出响应，以分析确定它们的性能参数。信号发生器是电子测量领域中最基本、应用最广泛的一类电子仪器。它可以产生多种波形信号,如正弦波,三角波,方波等,因而此次课程设计旨在运用模拟电子技术知识来制作一个能同时输出正弦波、方波、三角波的信号发生器。 三、设计目的： 1、进一步掌握模拟电子技术知识的理论知识，培养工程设计能力和综合分析能力、解决问题的能力。 2、基本掌握常用电子电路的一般设计方法，提高电子电路的设计和实验能力。 3、学会运用Multisim仿真软件对所做出来的理论设计进行仿真测试，并能进一步解决出现的基本问题，不断完善设计。 4、掌握常用元器件的识别和测试，熟悉万用表等常用仪表，了解电路调试的基本方法，提高实际电路的分析操作能力。 5、在仿真结果的基础上，实现实际电路。 四、设计内容及要求： 1、要求完成原理设计并通过Multisim软件仿真部分 （1）RC桥式正弦波产生电路，频率分别为300Hz、1KHz、10KHz、500KHz，输出幅值300mV~5V可调、负载1KΩ。 （2）占空比可调的矩形波电路，频率3KHz，占空比可调范围10%~90%，输出幅值3V、负载1KΩ。 （3）占空比可调的三角波电路，频率1KHz，占空比可调范围10%~90%，输出幅值3V、负载1KΩ。 （4）多用信号源产生电路，分别产生正弦波、方波、三角波，频率范围
数字信号处理课程设计报告
抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于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)
随机信号分析课程设计完整版
随机信号分析课程设计 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】
随机信号分析课程设计 一、题目： 设计一个抑制载波的复AM 信号，载波为40MHz ，接收带宽10MHz ，调制信号带宽50KHz ，加入高斯白噪声，带内信噪比10dB ： 1．1画出加噪后信号时域波形； 画出功率谱密度； 画出其同相、正交分量的功率谱； 统计方法画出包络概率密度。 二、问题分析： 画出加噪后信号时域波形： 首先，由已知条件先采样产生抑制载波的实AM 离散信号sr ，经过Hilbert 变换求得其解析信号s0，并经过低通滤波器，截止频率fs 限制在接收带宽下，加入噪声v ，得到此时的复AM 信号s ，再画出此时得到的复AM 信号时域波形。 画出功率谱密度： 将信号s 进行fft 变换后求得其功率谱密度，画出图形。 画出其同相、正交分量的功率谱： 将信号s 分解为正交分量)(cos )()(t t A t A s φ=和同相分量 )(sin )(t t A A c φ=，进行fft 变换得到功率谱密度，画出图形。 三、程序代码： f0=4*10^7;%载波信号频率40MHz f1=5*10^4;%调制信号频率50kHz fs=1*10^7;%接收机带宽采样频率10Mhz N=40001;%采样点数 %t=(0:1:N-1)/fs; f =10*f0; %画图范围设置 t0 = 5/f1;
t = 0:1/f:t0; k=1; sr=k*cos(2*pi*f1*t).*cos(2*pi*f0*t);%实am信号 % figure(2) % plot(t,sr) s0=hilbert(sr);%复am信号 h=sin(fs*t)/(pi*t); s0=conv(s0,h); am=max(abs(s0)); % % -------加噪方案（由加噪后信噪比确定高斯白噪声）----- snr=10; %设定加入白噪声后的信噪比为10db(均值为0) Pv=(am/(10^(snr/20)))^2;%噪声方差 % % -------------------------------------------------- % % % ---------加噪声------------- v=rand(1,N); v=v*sqrt(Pv);%白噪声 s=s0+v;%信号加噪声 % % ----------信号画图------------- figure(1) subplot(2,1,1),plot(t,s0); axis([*10^(-4) *10^-4 -10 10]) title('原始信号') subplot(2,1,2),plot(t,s); title(['加噪信号信噪比= ',num2str(snr),' dB. 噪声方差= ',num2str(Pv)]) axis([*10^(-4) *10^-4 -10 10]) %%----------画功率谱--------------- s1=detrend(s);%去趋势 ffs=abs(fft(s1)); theta=angle(s1)-2*pi*f0*t; a=abs(s1); ffs=ffs.*conj(ffs)*2/N;%频谱 %ffs=ffs.^2;%功率谱 figure(2) plot(ffs(1:N/2)); title('加噪信号功率谱') axis([3500 4500 0 4*10^4]) xlabel('*10^4') %%------------画正交同相分量功率谱---------- ac=s.*cos(2*pi*f0*t)-j*(hilbert(s)-s).*sin(2*pi*f0*t); as=-s.*sin(2*pi*f0*t)-j*(hilbert(s)-s).*cos(2*pi*f0*t); as1=detrend(as); ffas=abs(fft(as1)); ffas=(abs(ffas)).^2*2/N; ac1=detrend(ac); ffac=abs(fft(ac1)); ffac=(abs(ffac)).^2*2/N;
数字信号课程设计报告
一、语音信号去噪处理 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('幅度');
信号与系统课程设计报告 信号与系统课程设计题目
信号与系统课程设计报告信号与系统课程设计题目 信号与系统课程设计报告 ——频分复用通信系统的仿真设计 指导老师：XXX 小组成员： 摘要： 通过对信号与系统这门课程第八章通信系统学习，我们对频分复用（FDMA ）技术产生了浓厚的兴趣，于是决定自己利用MATLAB 强大的仿真功能来对频分复用系统进行仿真。本文首先录制三段不同的语音信号。然后通过推导，确定合适的载波信号的频率，对信号进行调制，调制后整合到一个复用信号上。再在复用信号上加一个随机的高斯白噪声得到在信道中传输的信号。之后根据通过对复用信号的频谱分析，得出切比雪夫滤波器的各项参数，通过设计好的滤波器进行信号分离后分别根据载波信号进行解调，再通过一个低通滤波器，得到原始信号。通过此次对FDMA 的仿真，我们更清楚了解了频分复用的工作原理，以及AM 调制解调方法，和滤波器的设计方法。频分复用技术对与通信系统节省资源有着重要的意义。
关键词： 频分复用 MATLAB 高斯白噪声 引言： 在电话通信系统中，语音信号频谱在300—3400Hz 内，而一条干线的通信资源往往远大于传送一路语音信号所需的带宽。这时，如果用一条干线只传一路语音信号会使资源大大的浪费，所以常用的方法是“复用”，使一条干线上同时传输几路电话信号，提高资源利用率。 本文是基于MATLAB 的简单应用，首先录制三段不同的语音信号。然后选择合适的高频载波，对信号进行调制，调制后整合到一个复用信号上。确定合适的信噪比，在复用信号上加一个随机的高斯白噪声得到在信道中传输的信号。之后根据载波信号设计合适的带通滤波器将三种信号进行分离，信号分离后分别进行同步解调，再通过一个低通滤波器，得到通过频分复用系统传输后得到的各个信号，将得到的信号与原信号对比，要保证信号与原信号吻合较好。 正文：
信号课程设计
实验一 时域采样与频域采样定理的验证实验 1. (1) 时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息； (2) 要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。 2. 时域采样定理的要点是： ① 对模拟信号()a x t 以T 进行时域等间隔理想采样，形成的采样信号的频谱 会 以采样角频率Ωs （Ωs=2π/T ）为周期进行周期延拓。公式为 ② 采样频率Ωs 必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。 利用计算机计算上式并不方便，下面我们导出另外一个公式，以便在计算机上进行实验。 理想采样信号 和模拟信号()a x t 之间的关系为： 对上式进行傅里叶变换，得到： 上式中，在数值上x a (nT)＝x(n)，再将ω=ΩT 代入，得到： 上式的右边就是序列的傅里叶变换，即 上式说明采样信号的傅里叶变换可用相应序列的傅里叶变换得到，只要将自变量ω用ΩT 代替即可。 频域采样定理的要点是： ① 对信号x(n)的频谱函数在［0，2π］上等间隔采样N 点，得到: ?(j )a X Ωa a a s 1??(j )FT[()](j j ) k X x t X k T ΩΩΩ∞ =-∞ ==-∑a ?()x t a a ?()()()n x t x t t nT δ∞ =-∞ =-∑j a a ?(j )[()()]e d t n X x t t nT t ΩΩδ∞ ∞--∞=-∞ =-∑?j a ()()e d t n x t t nT t Ωδ∞ ∞ --∞ =-∞ -∑?＝j a a ?(j )()e nT n X x nT ΩΩ∞ -=-∞ =∑j a ?(j )(e )T X X ωωΩΩ==j 2π()(e ) , 0,1,2, ,1 N k N X k X k N ωω= ==-
数字信号处理课程规划报告
数字信号处理课程设计报告《应用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('频域信号'); 运行程序得到的图形：
信号与系统课程设计报告分析
成绩评定表 课程设计任务书
摘要 本文研究的是傅里叶变换的对称性和时移特性，傅里叶变换的性质有：对称性、线性（叠加性）、奇偶虚实性、尺度变换特性、时移特性、频移特性、微分特性、积分特性、卷积特性（时域和频域）；从信号与系统的角度出发,给出了激励信号的具体模型；应用Matlab软件进行仿真，将研究的信号转化成具体的函数形式，在Matlab得到最终变换结果。使用傅里叶变换的方法、卷积的求解方法以及函数的微分等方法研究题目。 关键词: 傅里叶变换;对称性;时移特性;Matlab 目录 1、Matlab介绍................................................................................... 错误!未定义书签。
2.利用Matlab实现信号的频域分析—傅里叶变换的对称性与时移特性设计 (5) 2.1.傅里叶变换的定义及其相关性质 (5) 2.2.傅里叶变换的对称性验证编程设计及实现 (7) 2.3.傅里叶变换的时移特性验证编程设计及实现 (10) 3.总结 (13) 4.参考文献 (13) 1 、Matlab介绍 MATLAB作为一种功能强大的工程软件，其重要功能包括数值处理、程序设计、可视化显示、图形用户界面和与外部软件的融合应用等方面。 MATLAB软件由美国Math Works公司于1984年推出，经过不断的发展和完善，如今己成为覆盖多个学科的国际公认的最优秀的数值计算仿真软件。MATLAB具备强大的数值计算能力，许多复杂的计算问题只需短短几行代码就可在MATLAB中实现。作为一个跨平台的软件，MATLAB已推出Unix、Windows、Linux和Mac等十多种操作系统下的版本，大大方便了在不同操作系统平台下的研究工作。 MATLAB软件具有很强的开放性和适应性。在保持内核不变的情况下，MATLAB可以针对不同的应用学科推出相应的工具箱(toolbox)，目前己经推出了
课题信号与系统课程设计报告书
信号与系统课程设计 课程名称：信号与系统 题目名称：滤波器的设计与实现 学院：电气与电子工程学院 专业班级：电气工程及其自动化 学号：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、各种滤波器的幅频特性：
信号课程设计(郑巧健)
<<信号与线性系统>> 课程设计 姓名:郑巧健 学号:120901228 班级：自动化1204
目录 实验一连续信号的时域分析 3 实验二连续时间系统的时域分析 8 实验三连续信号的频域分析 15 实验四连续系统的频域分析 21 实验五信号采样与重建 28 实验六离散时间信号和系统分析 32
实验一 一、实验目的 1、熟悉MATLAB软件。 2、掌握常用连续信号与离散信号的MA TLAB表示方法。 二、实验设备 安装有matlab6.5以上版本的PC机一台。 三、实验原理
四、实验内容 1、用MATLAB表示连续信号：Aeαt，0 A cos(ωt +?)，0 A sin(ωt +?)。 clc clear all close all
A=10;a=2; f=A*exp(a*t) ezplot(f,[-10,10]); xlabel('t'); title('f=10*exp(2*t)'); grid on >> clc clear all close all syms t A=10;w=2;q=5; f=A*sin(w*t+q) ezplot(f,[-10,10]); xlabel('t'); title('f=10*sin(2*t+5)'); grid on >> clc clear all close all syms t A=10;w=2;q=5; f=A*cos(w*t+q) ezplot(f,[-10,10]); xlabel('t'); title('f=10*cos(2*t+5)'); grid on 2、用MATLAB表示抽样信号(sinc(t))、矩形脉冲信号(rectpuls(t, width)) 及三角脉冲信号(tripuls(t, width, skew))。 >> clc clear all close all t=-10:0.01:10 f=sinc(t)
郑州大学数字信号处理课程设计报告
实验一：基于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点谱分析，画出连续幅度谱，观察频谱泄露现象。
数字信号处理实验报告92885
目录 实验1 离散时间信号的频域分析-----------------------2 实验2 FFT算法与应用-------------------------------7 实验3 IIR数字滤波器的设计------------------------12 实验4 FIR数字滤波器的设计------------------------17
实验1 离散时间信号的频域分析 一．实验目的 信号的频域分析是信号处理中一种有效的工具。在离散信号的时域分析中，通常将信号表示成单位采样序列δ（n ）的线性组合，而在频域中，将信号表示成复变量e n j ω-或 e n N j π2-的线性组合。通过这样的表示，可以将时域的离散序 列映射到频域以便于进一步的处理。 在本实验中，将学习利用MATLAB 计算离散时间信号的DTFT 和DFT,并加深对其相互关系的理解。 二、实验原理 (1)DTFT 和DFT 的定义及其相互关系。序列x(n)DTFT 定义为()jw X e = ()n x n e ∞ =∞ ∑ω jn -它是关于自变量ω的复函数，且是以2π为周期的连续函数。 ()jw X e 可以表示为()()()jw jw jw re im X e X e jX e =+,其中，()jw re X e 和()jw im X e 分别是 ()jw X e 实部和虚部；还可以表示为 ()jw X e =()|()|jw j w X e e θ，其中， |()|jw X e 和{} ()arg ()j w X e ωθ=分别是()jw X e 的幅度函数和相位函数；它们都是ω的实函数，也是以2π为周期的周期函数。 序列()x n 的N 点DFT 定义为2211 ()()()()N N j k j kn kn N N N N n X k X e x n e x n W π π ---==== ∑∑，()X k 是周期为N 的序列。()j X e ω与()X k 的关系：()X k 是对()j X e ω）在一个周期 中的谱的等间隔N 点采样，即 2k |()()|jw w N X k X e π = = ，而()j X e ω 可以通过对()X k 内插获得，即
数字信号课程设计
《数字信号》课程设计报 告 学院：信息科学与工程 专业班级：通信1201
一、 目的与要求 是使学生通过上机使用Matlab 工具进行数字信号处理技术的仿真练习，加深对《信号分析与处理（自）》课程所学基本理论和概念的理解，培养学生应用Matlab 等工具进行数字信号处理的基本技能和实践能力，为工程应用打下良好基础。 二、 主要内容 1．了解Matlab 基本使用方法，掌握Matlab 数字信号处理的基本编程技术。掌握数字信号的基本概念。 2.用Matlab 生成几种典型数字信号（正弦信号、矩形信号、三角波信号等），并做幅频特性分析 2．Matlab 编程实现典型离散信号（正弦信号、矩形信号、三角信号）的离散傅立叶变换，显示时域信号和频谱图形（幅值谱和相位谱）；以正弦周期信号为例，观察讨论基本概念（混叠、泄漏、整周期截取、频率分辨率等）。 3．设计任意数字滤波器，并对某类型信号进行滤波，并对结果进行显示和分析。 4．利用matlab 求解差分方程，并做时域和频域分析。用matlab 函数求解单位脉冲响应，并利用窗函数分离信号。 5.用matlab 产生窗函数，并做世玉和频域分析。 6.显示图像，理解图像的模型，将图像进行三原色分解和边缘分析。 三．课程设计题目 一、 1） 生成信号发生器：能产生频率（或基频）为10Hz 的周期性正弦波、三角波和方波信号。绘出它们的时域波形 2） 为避免频谱混叠，试确定各信号的采样频率。说明选择理由。 3）对周期信号进行离散傅立叶变换，为了克服频谱泄露现象，试确定截取数据的长度，即信号长度。分析说明选择理由。 4）绘出各信号频域的幅频特性和相频特性 5）以正弦周期信号为例，观察讨论基本概念（频谱混叠、频谱泄漏、整周期截取等）。 二、已知三个信号()i a p n ，经调制产生信号3 1 ()()cos(/4)i i s n a p n i n π==∑，其中i a 为常 数，()p n 为具有窄带特性的Hanning 信号。将此已调信号通过信道传输，描述该信道的差分方程为 得到接收信号()()*()y n s n h n = 1）分析Hanning 信号()p n 的时域与频域特性 2）分析已调信号()s n 的时域与频域特性 () 1.1172(1)0.9841(2)0.4022(3)0.2247(4) 0.2247()0.4022(1)0.9841(2) 1.1172(3)(4)y n y n y n y n y n x n x n x n x n x n --+---+-= --+---+-
《数字信号处理》课程设计,基于MATLAB的音乐信号处理和分析解析
《数字信号处理》课程设计设计题目：基于MATLAB的音乐信号处理和分析 院系：物理工程学院 专业：电子信息科学与技术 学号： 姓名：
一、课程设计的目的 本课程设计通过对音乐信号的采样、抽取、调制解调、滤波、去噪等多种处理过程的理论分析和MATLAB实现，使学生进一步巩固数字信号处理的基本概念、理论以及频谱分析方法和数字滤波器设计方法；使学生掌握的基本理论和分析方法只是得到进一步扩展；使学生能有效地将理论和实际紧密结合；增强学生软件编程实现能力和解决实际问题的能力。 二、课程设计的基本要求 1 学会MATLAB的使用，掌握MATLAB的基本编程语句。 2 掌握在Windows环境下音乐信号采集的方法。 3 掌握数字信号处理的基本概念、基本理论和基本方法。 4 掌握MATLAB设计FIR和IIR数字滤波器的方法。 5 掌握使用MATLAB处理数字信号、进行频谱分析、涉及数字滤波器的编程方法。 三、课程设计内容 实验1音乐信号的音谱和频谱观察 使用windows下的录音机录制一段音乐信号或采用其它软件截取一段音乐信号（要求：时间不超过5s、文件格式为wav文件） ①使用wavread语句读取音乐信号，获取抽样率；（注意：读取的信号时双声道信号，即为双列向量，需要分列处理）； ②输出音乐信号的波形和频谱，观察现象； 使用sound语句播放音乐信号，注意不同抽样率下的音调变化，解释现象。 程序如下： [Y,FS,NBITS]=WAVREAD('怒放的生命 - 汪峰5s'); %读取音乐信号 plot(Y); %显示音乐信号的波形和频谱 sound(Y,FS); %听音乐（按照原来的抽样率） Y1=Y(:,1); %由双声道信号变为单声道信号 size(Y1) figure subplot(2,1,1);