CHANGSHA UNIVERSITY OF SCIENCE & TECHNOLOGY
DSP课程设计
题目： DSP
——音频信号频谱分析
学生姓名：翦杰
学号： 2
班级: 07-01
专业：电子信息工程
指导教师：黄亚飞、肖鸿
实习起止时间: 2010年12月27日至2011年1月7日
音频信号频谱分析
摘要
随着信息技术革命的深入和计算机技术的飞速发展，数字信号处理技术已经逐渐发展成为一门关键的技术学科。而DSP芯片的出现则为数字信号处理算法的实现提供了可能。这一方面极大地促进了数字信号处理技术的进一步发展；另一方面，它也使数字信号处理的应用领域得到了极大的拓展。在国外DSP芯片已经被广泛地应用于当今技术革命的各个领域；在我国，DSP技术也正以极快的速度被应用到科技和国民经济的各个领域。本次课程设计介绍了音频信号频谱分析的原理以及其所涉及的硬件结构和软件设计，该设计是基于快速傅立叶变换（FFT）的方法对采集的音频信号进行频谱分析，得到音频信号的频率及功率，FFT算法采用TLC320AD50编写DSP程序实现。现可以完成256点的FFT运算，频率分辨率达到100Hz，输入信号电压（峰峰值）可以达到100mV到4V。
关键词：音频信号；快速傅立叶变换；频谱分析；分辨率
目录
1 绪论 (1)
1.1 课题背景 (1)
1.2 课题目的 (1)
2 硬件电路............................................................................................ ..... ..... . (3)
2.1 系统框图......................................................................................... (3)
2.2 信号处理部分..................................................................................................... . (3)
2.2.1 信号叠加电路.................................................................................. ......... .. (3)
2.2.2 低通滤波器......................................................................................... (3)
2.2.3 50Hz陷波电路......................................................................... (4)
2.2.4 电平转换电路... ......... ............................................... .......................... (4)
2.2.5 信号采集部分............ ........................................................................ (4)
2.3 时钟信号产生电路.. ... .......................................... ................................. (13)
3 软件设计.............................................................. ..... .... ...................................... (14)
3.1 DSP初始化................................................................... ... (14)
3.2 AD50初始化...................................................................................... . .. (14)
3.3 设置DSP中断..... ... ................................................................................... (14)
3.4 数据存储... .... ... ................................................ ....................................... ...... (14)
3.5 数据输出.... .... ... ..................................................... ......................................... .. (14)
3.6 FFT算法... .... ... ............................................................ ................................... (14)
4 总结................................................................................................. .. (17)
参考文献................................................................................... .. (19)
附录................................................................................................. (20)
1 绪论
1.1 课题背景
目前，在微电子技术发展的带动下，DSP芯片的发展日新月异，DSP的功能日益强大，性能价格比不断上升，开发手段不断改进。DSP芯片已经完全走下了“贵族”的圣坛。DSP芯片已经在通信与电子系统、信号处理系统、自动控制、雷达、军事、航空航天、医疗、家用电器、电力系统等许多领域中得到了广泛的应用，而且新的应用领域在不断地被发掘。
TI、AD、AT&T、Motorola和Lucent等公司是DSP芯片的主要生产商。其中TI
公司的TMS320系列的DSP占据了全球DSP市场的50%左右。该系列产品在我国同样被用户广泛使用，市场份额更高，超过90％。
鉴于多数DSP芯片和高速A/D、D/A芯片工艺为贴片封装，对一般用户来说工厂制版成本较高、手工工艺难制版、效果差等的困难，本系统设计了一套基于DSP芯片的最小系统板，并扩展了A/D、D/A实现语音信号的采集和回放，制作语音处理平台。设计的核心芯片采用TI公司的TMS320VC5402PGE100进行设计，其最高处理速度能达到100MIPS(每秒执行100百万条指令)，性能优越、性价比高，适合大多数用户和教学科研。
基于DSP的运用领域和前景，结合我的专业跟个人爱好，本次毕业设计所选课题为DSP语音采集回放处理平台。以TMS320C5402DSP为核心，对外部语音信号进行采集，并对所采集信号进行语音处理，最后通过外部设备回放。该系统适合对单语音信号进行处理。由于设计过程中采用的A/D、D/A芯片是TI公司的TLC320AD50，所以最高采样速率为22.05KHZ。为了验证本次设计的正确性和可用性，对采集的语音信号进行FIR滤波，滤除50HZ交流信号，并抑制频率在3600HZ以上的语音信号。所设计的滤波器是带通滤波器，通带为200HZ—3400HZ，经过实验验证，得到了预期的滤波效果。证明本次单通道语音信号处理平台设计的正确性，可用性。
1.2 课题目的
DSP课程设计是对《数字信号处理》、《DSP原理及应用》等课程的较全面
练习和训练，是实践教学中的一个重要环节。
通过本次课程设计，综合运用数字信号处理、DSP技术课程以及其他有关先
修课程的理论和生产实际知识去分析和解决具体问题，并使所学知识得到进一步
巩固、深化和发展。初步培养学生对工程设计的独立工作能力，掌握电子系统设
计的一般方法。通过课程设计完成基本技能的训练，如查阅设计资料和手册、程
序的设计、调试等，提高学生分析问题、解决问题的能力。
本题目通过TLC320AD50采集音频信号(f.max<10kHz)，编写DSP的FFT处理
程序（自定频谱分辨力），获得幅频特性之后，在点阵液晶(128*64)中大致显示出
幅频图。并在液晶中用文字显示频率幅值前三的频率值。
1、DSP与TLC320AD50接口电路的原理图绘制；
2、DSP控制TLC320AD50的程序编写与调试；
3、TLC320AD50进行语音模拟量到数字信号的转换，实现声音的采集，在CCS 软件中分析信号的幅频特性；
4、编写DSP的FFT处理程序；
5、控制点阵液晶，实现绘图功能，将幅频图显示出来；
6、按要求编写课程设计报告书，正确、完整的阐述设计和实验结果；
7、在报告中绘制程序的流程图，并文字说明。
2 硬件电路
2.1 系统框图
本系统的硬件电路主要由信号预处理模块、信号采集模块、信号频谱分析模块、时钟触发信号、系统控制模块和键盘显示模块6部分组成。系统组成框图如图2.1所示。
图2.1 系统组成框图
2．2 信号处理部分
信号预处理具体电路主要由信号合成、阻抗匹配、0.5f的低通滤波、50H的陷波器和信号抬高等部分组成，如图2.2所示。
图2.2 信号处理框图
2.2.1 信号叠加电路
信号叠加电路采用反向求和电路，将3信号求和的输出，运放采用OP37。
2.2.2 低通滤波器
由于要求音频信号为10KHz以下，为了使信号得到更高的纯度，则需对不用的高频信号进行滤出，从而提高系统测量精度。通过滤波器的快速设计方法设计出12KHz 的低通滤波器。
2.2.3 50Hz陷波电路
在音频信号中，工频50HZ信号对系统的影响最大，因此，在信号处理中加入
50HZ的陷波电路。计算工公式为：f=1/2πRC。
2.2.4 电平转换电路
由于ADC只能采集正电压，则需对信号进行抬高，分析计算得电路需加3V的直流电平。
2.2.5 信号采集部分
语音信号模数／数模转换选择TLC320AD50（以下简称AD50）芯片，AD50使用过采样（over sampling）∑—Δ技术提供从数字信号到模拟信号（DA）和模拟信号到数字信号（AD）的高分辨率低速信号转换。该器件包括两个串行的同步转换通道（用于各自的数据传输），在DA之前有一个插人滤波器（interpolation filter），在AD之后有一个抽取滤波器（decimdtionfilter），由此可降低AD50自身的噪声。此外，AD50还具有片内时序和控制功能。
AD50特点如下。
·输入信号。单端信号输人，幅度在1～4V之间。
·输出信号。单端信号输出，幅度在1～4V之间。
·单一 5V电源供电，也可以使用5V模拟电源和3V数字电源同时供电。
·最大工作功耗为100mW。
·通用16位数据格式，也可以采用2的补码数据格式。
·内部基准电压。
· AD为64倍采样，DA为256倍采样。
·支持各种V．34协议的采样速率。
·具有多种可选的采样频率。
·支持商业级音响应用。
·工作温度范围从－40～850。
AD50的引脚分布如图2.3所示（（1）为DW封装的引脚分布，（2）为PT封装的引脚分布）。
图2.3 AD50的引脚分布AD50的各个引脚的说明如表2-1所示。
表2-1 AD50的引脚说明
续表
在AD50正常工作前，必须对它进行正确的初始化。初始化操作的主要工作是配置AD50的4个控制寄存器CR1、CR2、CR3和CR4。控制寄存器的读写是通过二次通信来实现。在二次通信中，D0～D7为写人控制寄存器的数据或者从控制寄存器读出的数据，D8～D12的内容决定选择哪个控制寄存器，D13位决定是读操作还是写操作。由D8～D13位确定的所有情况如表2-2所示。
表2-2 AD50的D8～D13位确定的所有操作情况
控制寄存器1中各个控制位的功能说明如表2-3所示。
表2-3 控制寄存器1的各个控制位的功能说明
控制寄存器2中各个控制位的功能说明如表2-4所示。
表2-4 控制寄存器2的各个控制位的功能说明
控制寄存器3中各个控制位的功能说明如表2-5所示。
表2-5 控制寄存器3的各个控制位的功能说明
控制寄存器4中各个控制位的功能说明如表2-6、2-7所示。
表2-6 控制寄存器4的各个控制位的功能说明
表2-7 控制寄存器4的第4～第6位对应的N值
AD50控制寄存器的设置必须在二次通信中完成。AD50有硬件和软件两种方式启动二次通信，下面介绍硬仵方式启动AD50的二次通信过程。硬件启动过程如下：C5409通过其内部寄存器将XF引脚变为高电平，从而控制AD50的FC引脚到高电平，然后向缓冲串口写一个16位的控制字，低8位是AD50的控制寄存器初始化值，高8位选择所要初始化的控制寄存器及操作。整个过程的时序如图2.4所示。
图2.4 硬件启动AD50二次通信时序
AD50的工作过程可分为AD通道工作过程和DA通道工作过程。AD通道把模拟信号转换成数字信号，并以2进制补码形式表示。当帧同步信号有效时（FS为低电平），16位（或者15位）数字信号在SCLK的上升沿输出到DOUT引脚，一位数据对应一个SCLK周期。传输时序如图2.5所示。
图2.5 AD50的AD通道数据传输时序
AD50的DA通道把送人的数字信号转换成模拟信号。在SCLK的作用下，数字信号通过DIN引脚进人DA通道，每个SCLK的下降沿输人一位数字信号。DA将输入的数字信号转换成模拟信号输出，DA通道的传输时序如图2.6所示。
图2.6 AD50的DA通道数据传输时序
AD50的前后端信号处理包括两个处理电路：输人模拟信号的处理电路和依据输出模拟信号的处理电路。这两个处理电路的主要作用是将信号进行处理，使之更加适合AD和DA的要求。
为了达到更好的信号处理效果，AD50的模拟信号输人一般采用差分输人方式，即使用两个运算放大器，将单端输入信号转换成差分输人信号，电路连接如图2.7所示。
图2.7 AD50单端信号转换成差分信号
由图可知，单端输人信号经过两个22μF的隔直电容，送人运算放大器的反相端，输出反相信号IMP；IMP再输人到另一个信号的反相端，输出同相信号INP，从而形成差分输人信号INP和IMP。图中的运算放大器选择的是TI公司的TLC4502，也可以选用其他可替代的运算放大器芯片。
除了使用上述的运算放大器外，也可以使用变压器将单端信号转换成差分信号输人到AD50，如图2.8所示。
图2.8 使用变压器将单端信号转换成差分信号
图中推荐的值为50Ω，Rin为22Ω，Cin为10μF，这些元件的值也可以根据具体的信号进行调整，一般清况下电阻值在10～100Ω之间，电容值在10～200pF之间。
AD50的DA输出为差分信号，可以直接驱动600Ω的负载。DA输出处理电路如图2.9所示。
图2.9 AD50的DA后端输出信号处理
AD50的去耦合电容的连接如图2.10所示。
图2.10 AD50的去耦合电容的连接
AD50的模拟电源和数字电源可以使用一个电感连接在一起，连接方法如图 2.11所示。
图2.11 AD50模拟和数字电源的连接
2.3 时钟信号产生电路
通过计算，对模拟信号的采样速率为25.6KHz,则需产生25.6KHZR的方波.作为
A/D转换数据的速率。可采用单片集成电路函数发生器ICL8038。ICL8038工作频率范围在几Hz到几百KHZ之间，可以同时输出三角波、方波和正统波，满足设计要求。
3 软件设计
3.1 DSP的初始化
包括定时器0的初始化，以确保TOUT0引脚能输出正确的时钟信号到AD50，使AD50能根据该时钟信号进行采样。初始化缓冲串口时，应使DSP工作于主设备方式，CLKR0和FSX0引脚设为输出。DSP初始化完成后，在DSP的TOUT0，
CLKR0，FSX0这3个引脚可以看到连续输出的脉冲。
3.2 AD50的初始化
此时主要是设置DSP的XF引脚输出为高，以启动AD50的二次通信，在正确设置AD50的4个控制寄存器以后，AD50开始采样数据。此时检测AD50的DOUT引脚能，能发现引脚有连续的信号输出。
3.3 设置DSP中断
DSP的中断负责从缓冲串口读取数据，可以在CCS中查看读取的数据是否正确，如果每次读取的都是0或者同一个数据，说明AD50没有正常工作，此时需检查DSP 或者AD50是否初始化成功。
3.4 数据存储
DSP从缓冲串口读取数据存放到RAM，连续存放500或者1000个数据，可以通过CCS的画图功能看到采集到的数据。这时主要判断采样频率是否正确、有没有数据丢失等。
3.5 数据输出
将RAM中的数据直接送到AD50的D/A单元，转换成模拟信号后输出。这时从AD50的模拟信号输出端可以看到和输入端一样的信号。
3.6 FFT算法
本系统软件设计主要实现FFT算法和各频率和对应功率的分析和显示，由于DSP 语言表达和运算能力比较强且具有很好的可移植性和硬件控制能力，则采用DSP语言实现。总体程序流程图如图3.1所示。
图3.1 总体程序流程图
本设计中FFT根据DFT的基二分解方法，在第L（L表示从左到右的运算级数，L=1，2，3…M）级中，每个蝶形的两个输入数据相距B（L-1）个点，同一旋转因子对应着间隔为2L点的2（M-L）个蝶形。从输入端开始，逐级进行，在进行L级运算时，依次求出2（L-1）不同的旋转因子，每求出一个旋转因子，就计算完它对应的所有的2（M-L）个蝶形。因此用三重循环程序实现FFT变换。同时采用原址计算，可节省大量的存储单元。FFT算法的程序框图如图3.2所示。
图3.2 FFT算法的程序框图
系统用于实现语音信号的采集和播放，收到数据后，AD50对数据进行模数转换，变成数字语音信号后存储于内部RAM，然后输出到AD50，AD50进行数模转换后输出到相应的后端处理电路，最后由后端电路将输入的语音信号显示。
抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于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)
2007年A题音频信号分析仪 本系统基于Altera Cyclone II 系列FPGA嵌入高性能的嵌入式IP核(Nios)处理器软核，代替传统DSP芯片或高性能单片机，实现了基于FFT的音频信号分析。 音频信号分析仪 山东大学王鹏陈长林秦亦安 摘要：本系统基于Altera Cyclone II 系列FPGA嵌入高性能的嵌入式IP核(Nios)处理器软核，代替传统DSP芯片或高性能单片机，实现了基于FFT的音频信号分析。并在频域对信号的总功率，各频率分量功率，信号周期性以及失真度进行了计算。并在FPGA中嵌入了8阶IIR切比雪夫（Chebyshev）II型数字低通滤波器，代替传统有源模拟滤波器实现了性能优异的音频滤波。配合12位A/D转换芯片AD1674，和前端自动增益放大电路，使在50mV到5V的测量范围下，单一频率功率及总功率测量误差均控制在1%以内。 关键词：FPGA；IP核；FFT；IIR；可控增益放大 Abstract: This system is based on IP core(Nios)soft-core processors embedded in the FPGA of Altera Cyclone II family. Instead of using DSP or microcontroller, we use Nios II to perform a low-cost FFT-based analysis of the audio signal.And we caculated the power of the whole signal,the power of each frequence point that componented the signal.By the way,we anlysised its periodicity and distortion.We also embedded an 8-order Chebyshev II IIR digital low-pass filter to replace the traditional analog Active Filter to perform an excellent audio filter. With 12bit A / D converter chip AD1674, and the front-end automatic gain amplifier, this system’s single-frequency power and total power measurement error is below 1% in 50mV to 5V measurement range. Keyword: FPGA；IP core; FFT;IIR; automatic gain amplifier 一、方案选择与论证 1、整体方案选择 音频分析仪可分为模拟式与数字式两大类。 方案一：以模拟滤波器为基础的模拟式频谱分析仪。有并行滤波法、扫描滤波法、小外差法等。因为受到模拟滤波器滤性能的限制，此种方法对我们来说实现起来非常困难。 方案二：以FFT为基础的的数字式频谱分析仪。通过信号的频谱图可以很方便的得到输入信号的各种信息，如功率谱、频率分量以及周期性等。外围电路少，实现方便，精度高。 所以我们选用方案二作为本音频分析仪的实现方式。
交通信号 控制系统（ATC）设计方案 x x x x有限责任公司
目录 1.概述 (1) 1.1系统简介 (1) 1.2设计原则 (2) 1.3系统设计依据及执行标准 (4) 2．总体设计方案 (6) 2.1控制系统总体功能 (6) 2.2通信系统总体结构 (6) 2.3通信系统主要优势 (8) 3.详细设计方案 (9) 3.1监测点设备 (9) 3.1.1设备功能描述 (9) 3.1.2监测点设备组成、结构及特点 (9) 3.2防雷保护及安全设计 (14) 3.3详细设备说明 (15) 3.3.1高清晰摄像机 (15) 3.3.2标清视频检测 (15) 3.3.3补光设备 (15) 3.3.4嵌入式存储 (15) 3.3.5 GOE210千兆工业以太网交换机 (15) 3.3.6 POE工业以太网光纤收发器 (17) 3.4系统典型配置清单 (18)
1.概述 城市发展交通智能信号灯，减少道路拥堵，最终达到智能化区域交通信号控制系统。智能交通信号灯迎合实现绿色经济的时代潮流，为了解决这个问题，提出智能交通信号灯及网络技术，会根据路口车辆多少，自动调节时间，可减少等候时间在75%以上，从而大大节省了人们的出行时间，减少了路口的无效等候，使出行更快捷。 在智能交通系统中，以往的常规摄像机是对所有通过该地点的机动车辆的车牌进行拍摄、记录与处理。由于受到图像采集设备分辨率的制约，图片仅能反映出车型、车身颜色、车牌号码等简单信息。公安执法部门对部分治安案件、交通肇事案件的取证要求上，希望能掌握更详细更清楚的资料，如驾驶员的面貌特征、车内驾驶室的情况、清晰的车辆信息、货车的装载情况。采用高清晰摄像机做前端采集，可以实现所抓拍的图像中用肉眼清楚地分辨：车辆的颜色、特征、车牌的号码、车牌颜色、司乘人员的面部特征。 如此一来智能化同时也带来了网络数据流量的剧增，对网络通信的可靠传输提出了更高的要求。工业以太网交换机在区域交通信号控制系统网络中稳定性、高可靠性、高安全性成为关键中的关键。 1.1系统简介 区域交通信号控制系统(ATC) 智能化区域交通信号控制系统采用百万像素的数字化网络摄像机（1600×1200 CCD传感器），一台摄像机覆盖两条车道，准确抓拍正常行驶、压线行驶、并行通过的车辆，并自动识别车牌号码，抓拍的车辆图片可清晰地显示车辆特征及前排司乘人员的面部特征。摄像机工作于外触发方式，通过视频分析、环形线圈或者窄波雷达检测通过车辆，在抓拍车辆的同时可获取车辆的行驶速度。两条车道共用一台高清数字摄像机的方式在保障系统性能的前提下，大大降低了系统成本。
哈尔滨工程大学 实验报告 实验名称：离散时间滤波器设计 班级：电子信息工程4班 学号： 姓名： 实验时间：2016年10月31日18：30 成绩：________________________________ 指导教师：栾晓明 实验室名称：数字信号处理实验室哈尔滨工程大学实验室与资产管理处制
实验七音频频谱分析仪设计与实现 一、 实验原理 MATLAB 是一个数据分析和处理功能十分强大的工程实用软件，其数据采集工具箱为实现数据的输入和输出提供了十分方便的函数命令。本实验要求基于声卡和MTLAB 实现音频信号频谱分析仪的设计原理与实现，功能包括： (1)音频信号输入，从声卡输入、从WAV 文件输入、从标准信号发生器输入； (2)信号波形分析，包括幅值、频率、周期、相位的估计、以及统计量峰值、均值、均方值和方差的计算。 (3)信号频谱分析，频率、周期的统计，同行显示幅值谱、相位谱、实频谱、虚频谱和功率谱的曲线。 1、频率(周期)检测 对周期信号来说，可以用时域波形分析来确定信号的周期，也就是计算相邻的两个信号波峰的时间差、或过零点的时间差。这里采用过零点(ti)的时间差T(周期)。频率即为f = 1/T ，由于能够求得多个T 值(ti 有多个)，故采用它们的平均值作为周期的估计值。 2、幅值检测 在一个周期内，求出信号最大值ymax 与最小值ymin 的差的一半，即A = (ymax - ymin)/2，同样，也会求出多个A 值，但第1个A 值对应的ymax 和ymin 不是在一个周期内搜索得到的，故以除第1个以外的A 值的平均作为幅值的估计值。 3、相位检测 采用过零法，即通过判断与同频零相位信号过零点时刻，计算其时间差，然后换成相应的相位差。φ=2π(1-ti/T)，{x}表示x 的小数部分，同样，以φ的平均值作为相位的估计值。 频率、幅值和相位估计的流程如图1所示。 4、数字信号统计量估计 (1) 峰值P 的估计 在样本数据x 中找出最大值与最小值，其差值为双峰值，双峰值的一半即为峰值。 P=0.5[max(yi)-min(yi)] (2)均值估计 i N i y N y E ∑== 1 )( 式中，N 为样本容量，下同。 (3) 均方值估计 () 20 2 1 ∑== N i i y N y E (4) 方差估计 ∑=-=N i i Y E y N y D 0 2))((1)(
应用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');
语音信号的数字滤波 ——利用双线性变换法实现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段上，可以通过以下的正切变换实现
信息采集系统解决方案
信息采集系统解决方案 1系统概述 信息采集是信息服务的基础，为信息处理和发布工作提供数据来源支持。信息数据来源的丰富性、准确性、实时性、覆盖度等指标是信息服务的关键一环，对信息服务质量的影响至关重要。针对交通流信息数据，包括流量、速度、密度等，目前主要是基于微波、视频、地磁等固定车辆检测器以及浮动车等移动式车辆检测器进行采集，各种采集方式都存在响应的利弊。针对车驾管以及出入境数据，包括车辆信息、驾驶人信息、出入境办证进度信息等，主要是通过和公安相关的数据库进行对接，此类信息将在信息分析处理系统进行详细介绍。 针对目前交通信息来源的多样性以及今后服务质量水平发展对信息来源种类扩展要求，需要建设一套统一的，具备良好兼容性和前瞻性的交通信息统一接入接口。一方面，本期项目的各种交通信息来源可以使用该接口进行数据接入，另一方面，当新的或第三方的交通信息来源需要加入到本系统中来时，可以使用该接口进行数据接入，不需要再次投入资源进行额外开发。 统一接入接口建成后，根据各种数据来源系统的网络环境、系统技术特性和交通流信息数据特点，开发相应的交通信息数据对接程序，逐一完成微波采集系统、浮动车分析系统、人工采集等来源的交通信息数据采集接入。 2系统架构及功能介绍 2.1统一接入接口 统一接入接口的建设的关键任务包括接口技术规范制定、路网路段编码规则约定及交通信息数据结构约定等多个方面。
2.1.1接口技术规范 一方面由于本系统接入的交通信息数据来源多样，开发语言和系统运行的环境均存在差异，不具备统一的技术特性；另一方面，考虑到以后可能需要接入更多新的或第三方的信息系统作为数据来源，应当选择较成熟和通用的接口实现技术作为本项目的交通流信息采集统一接入接口实现技术。 根据目前信息系统建设的行业现状，选择Web Service和TCP/UDP Socket 作为数据传输接口的实现技术是较优的选择。Web Service和TCP/UDP Socket 具有实时性强、通用性强、应用广泛、技术支持资源丰富等优势，可以实现跨硬件平台、跨操作系统、跨开发语言的数据传输和信息交换。 项目实施时需要根据现有的信息采集系统的技术特点来具体分析，以选定采用Web Service或TCP/UDP Socket作为接口实现技术，必要时可以两种方式并举，提供高兼容度的接口形式。 为了保护接入接口及其数据传输的安全性，避免恶意攻击访问，避免恶意数据窃取，可以使用身份认证、加密传输等技术来加以保证。 统一数据采集接口的工作流程可以如下进行：
频谱分析仪 摘要频谱分析仪是研究电信号频谱结构的仪器，用于信号失真度、调制度、谱纯度、频率稳定度和交调失真等信号参数的测量，是一种多用途的电子测量仪器。随着软硬件技术的发展，仪器的智能化与虚拟化已成为未来实验室及研究机构的发展方向。虚拟仪器技术的优势在于可由用户定义自己的专用仪器系统，且功能灵活，很容易构建，所以应用面极为广泛。本文介绍了一种使用GUI工具箱用matlab实现的简易虚拟频谱分析仪的设计方法。 关键词matlab，频谱分析仪，时域分析，频域分析
目录 1概述 (3) 2技术路线 (4) 3实现方法 (5) 3.1搭建GUI界面 (5) 3.2信号输入 (6) 3.2.1选择信号输入 (6) 3.2.2声卡输入 (7) 3.2.3读取wav文件 (7) 3.2.4信号发生器输入 (7) 3.3时域分析 (8) 3.4频域分析 (9) 3.5仿真 (10) 3.5.1声卡输入 (10) 3.5.2读取wav文件 (10) 3.5.3信号发生器 (11) 4存在的问题 (15) 5致谢...................................................................................................... 错误！未定义书签。参考文献 (15)
1概述 MATLAB是矩阵实验室(Matrix Laboratory)的简称，是美国MathWorks公司出品的商业数学软件。可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。图形用户界面(Graphical User Interface，简称GUI，又称图形用户接口)是指采用图形方式显示的计算机操作用户界面。与早期计算机使用的命令行界面相比，图形界面对于用户来说在视觉上更易于接受。MATLAB自带了强大的GUl工具[1]。在本文中，将利用MATLAB的GUI工具，设计出数字频谱分析仪。 频谱分析仪是研究电信号频谱结构的仪器，用于信号失真度、调制度、谱纯度、频率稳定度和交调失真等信号参数的测量，可用以测量放大器和滤波器等电路系统的某些参数，是一种多用途的电子测量仪器。它又可称为频域示波器、跟踪示波器、分析示波器、谐波分析器、频率特性分析仪或傅里叶分析仪等。现代频谱分析仪能以模拟方式或数字方式显示分析结果，能分析1赫兹以下的甚低频到亚毫米波段的全部无线电频段的电信号[2]。目前已经有许多较成熟的频谱分析软件，如SpectraLAB、RSAVu、dBFA等[3]。本文将给出的则是通过MATLAB软件实现的基于FFT的数字频谱分析仪。 FFT（Fast Fourier Transformation），即为快速傅氏变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。它对傅氏变换的理论并没有新的发现，但是对于在计算机系统或者说数字系统中应用离散傅立叶变换，可以说是进了一大步[4]。 通过此次设计，能进一步掌握MATLAB软件开发过程的基本理论、基本知识和基本技能，熟悉基于MATLAB平台的若干信号处理系统开发及调试方法，且成本低，易于实现，容易修改，并可以进行仿真。该设计的进行可以为我们以后的学习工作奠定一定的基础。
数字信号处理 —综合实验报告 综合实验名称：应用MatLab对语音信号进行 频谱分析及滤波 系： 学生姓名： 班级： 学号： 成绩： 指导教师： 开课时间学年学期
目录 一．综合实验题目 (1) 二、综合实验目的和意义 (1) 2.1 综合实验目的 (1) 2.2 综合实验的意义 (1) 三．综合实验的主要内容和要求 (1) 3.2 综合实验的要求： (2) 四．实验的原理 (2) 4.1 数字滤波器的概念 (2) 4.2 数字滤波器的分类 (2) （1）根据单位冲激响应h(n)的时间特性分类 (2) 五．实验的步骤 (3) 下面对各步骤加以具体说明。 5.1语音信号的采集 (3) 5.2 语音信号的频谱分析; (3) 5.3 设计数字滤波器和画出其频率响应 (5) 5.3.1设计数字滤波器的性能指标： (5) 5.3.2 用Matlab设计数字滤波器 (6) 5.6 设计系统界面 (19) 六、心得体会 (20) 参考文献： (21)
一．综合实验题目 应用MatLab对语音信号进行频谱分析及滤波 二、综合实验目的和意义 2.1 综合实验目的 为了巩固所学的数字信号处理理论知识，使学生对信号的采集、处理、传输、显示和存储等有一个系统的掌握和理解，再者，加强学生对Matlab软件在信号分析和处理的运用 综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。 2.2 综合实验的意义 语言是我们人类所特有的功能，它是传承和记载人类几千年文明史，没有语言就没有我们今天人类的文明。语音是语言最基本的表现形式，是相互传递信息最重要的手段，是人类最重要、最有效、最常用和最方便的交换信息的形式。 语音信号处理属于信息科学的一个重要分支，大规模集成技术的高度发展和计算机技术的飞速前进，推动了这一技术的发展；它是研究用数字信号处理技术对语音信号进行处理的一门新兴学科，同时又是综合性的多学科领域和涉及面很广的交叉学科，因此我们进行语言信号处理具有时代的意义。 三．综合实验的主要内容和要求 3.1综合实验的主要内容： 录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；给定滤波器的性能指标，采用窗函数法和双线性变换设计滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；回放语音信号；综合实验应完成的工作： (1）语音信号的采集; (2）语音信号的频谱分析;
数字信号处理课程设计 姓名：刘倩 学号：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 算法的应用
Aｄobe Auditｉon系列教程（二)：频谱分析仪 频谱分析仪是研究信号频谱特征的仪器，在电子技术一日千里的今天,是研究、开发、调试维修中的有力武器。现代频谱分析仪都趋向于智能化，虚拟仪器技术广泛应用，有些就是以专用的计算机系统为核心设计的。其结果是结构大大简化、性能飞速提高。当然专业的频谱分析仪就比示波器更加昂贵了,业余爱好者更难用上。不过不必灰心,我们可以充分利用ＡdｏbｅAudition的频谱分析功能,让你拥有精确频谱分析仪的美梦成真!? 1. 频谱显示模式? Ａdoｂe Audition本身有一种“频谱显示”模式。先打开一段波形,或用《妙用Adobe Auｄｉtioｎ:数字存储示波器》一文介绍的方法录制一段波形，即可进行频谱分析。这里我们新建一段２0秒的对数扫频信号(本文大多选用直接建立的波形,以便了解信号原始波形的标准频谱特征），然后选择“Vｉeｗ=>Sｐectral Vieｗ”(视图=＞频谱)，如图1，或点击快捷工具栏的“Togｇle bｅtwｅen Ｓｐectral and Wavｅｆorm vｉews”(切换频谱视图／波形视图）按扭,即可将波形以频谱显示的方式显示出来，如图2。扫频的频谱显示见图３。 图1
图2 图3 可以看到，横轴为时间,纵轴为频率指示。每个时刻对应的波形频谱都被显示出来了，可以看到扫描速度是指数增加的，即将频率轴取对数时扫描速度是线性的。如图中光标处18秒处频谱指示约1１KＨz。实际上频谱指示的颜色是代表频谱能量的高低的,颜色从深蓝到红再到黄，指示谱线电平由低到高的变化。这实际上跟地图的地形鸟瞰显示是比较相似的，看图4频谱复杂变化的声音频谱就更容易理解这点了。 图4
课程设计任务书 课程设计学生日志
课程设计考勤表 课程设计评语表
音频信号的谱分析及去噪 一、研究背景：……………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………… 二、设计方案…………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………
三、设计目的和意义 通过MATLAB编程，用FFT函数绘制出音频信号的频谱。用噪音去干扰音频信号，画出干扰后信号的频谱。这样观察对比，便可分析出噪声对音频信号的干扰。主要目的是通过设计FIR 数字滤波器滤除噪音信号，体会滤波器可提取有用信号消除干扰的作用。去噪，可以减少或消除信号传输过程中的干扰，从而达到有效传输。…………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………………… 四、设计原理 1、FFT原理：运用快速傅里叶变换得信号的频谱，快速傅里叶变换是可以用计
算机编程实现的一种快速求得信号频谱的方法。在MATLAB 编程中可以直接调用FFT 函数。 2、FIR 数字滤波器的设计步骤： 1）将给定的数字滤波器性能指标转化成相应的模拟滤波器性能指标。 2）将模拟滤波器的性能指标变换成模拟低通滤波器的性能指标。 3）用所得的模拟低通滤波器的性能指标，利用某种模拟滤波器逼近方法，设计得出该滤波器。 低通滤波器的频率响应： 设计滤波器时主要考虑四个指标：通带截止频率、阻带起始频率、通带衰减、阻带衰减。 五、 详细设计步骤 1. 用Windows 附件中的录音机录制了一段歌曲，并且由MATLAB 中的FFT 函数绘制出它的频谱图和时域波形。 N=1024;%采样点数 fs=20000; %语音信号采样频率为20000 [x0,fs,NBTS]=wavread('d:\输入噪声\ly.wav'); % ly 为原始语音信号 [m0,n0]=size(x0); x0=wavread('d:\ 输 入 噪 声 \ly.wav'); %读取音频信号 sound(x0,fs); %播放音频信号 figure(1) 阻带 过渡带 通带
数字信号处理课程设计报告《应用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('频域信号'); 运行程序得到的图形：
AdoｂｅＡｕditｉon系列教程(二)：频谱分析仪 频谱分析仪是研究信号频谱特征的仪器,在电子技术一日千里的今天,是研究、开发、调试维修中的有力武器。现代频谱分析仪都趋向于智能化,虚拟仪器技术广泛应用，有些就是以专用的计算机系统为核心设计的。其结果是结构大大简化、性能飞速提高。当然专业的频谱分析仪就比示波器更加昂贵了,业余爱好者更难用上。不过不必灰心,我们可以充分利用AdobｅＡｕdiｔiｏn的频谱分析功能,让你拥有精确频谱分析仪的美梦成真! 1. 频谱显示模式 AｄobｅAudition本身有一种“频谱显示”模式。先打开一段波形,或用《妙用Adobe Ａuｄitｉｏn：数字存储示波器》一文介绍的方法录制一段波形，即可进行频谱分析。这里我们新建一段2０秒的对数扫频信号（本文大多选用直接建立的波形，以便了解信号原始波形的标准频谱特征）,然后选择“Viｅｗ=＞Spe ｃtrａl Ｖiew”（视图=＞频谱）,如图1，或点击快捷工具栏的“Ｔｏｇgle ｂetween Sｐeｃｔraｌ and Ｗaveform views”(切换频谱视图/波形视图)按扭，即可将波形以频谱显示的方式显示出来，如图2。扫频的频谱显示见图3。 图1
图２ 图3 可以看到，横轴为时间,纵轴为频率指示。每个时刻对应的波形频谱都被显示出来了,可以看到扫描速度是指数增加的，即将频率轴取对数时扫描速度是线性的。如图中光标处18秒处频谱指示约11KＨz。实际上频谱指示的颜色是代表频谱能量的高低的,颜色从深蓝到红再到黄,指示谱线电平由低到高的变化。这实际上跟地图的地形鸟瞰显示是比较相似的，看图4频谱复杂变化的声音频谱就更容易理解这点了。
交通信号控制系统解决方案 1概述 交通信号控制系统，是智能交通系统（ITS）在交通管理工作中的基本应用，也是城市智能交通管控系统中最直接、最基础的应用系统。通过建设信号控制系统，实现信号路口联网远程控制、交通流量的采集、路口自适应控制、绿波协调控制以及区域的自适应控制，有效减少车辆的停车次数，节省旅行时间；后台实时调整信号配时，采取多时段控制方式，必要时，可通过智能交通管理中心人工干预，直接控制路口交通信号机执行指定相位，有效的疏导交通，减少行车延误，提高通行能力，缓解日益严峻的城区道路交通拥堵压力，提高城区交通综合管理能力，减少汽车尾气排放，美化环境，提升城区形象。 2系统结构设计 系统结构划分为3级：分别为中心控制级设备、区域控制级设备以及路口控制级设备。交通信号控制系统设备主要包括中心设备、前段设备和通信设备。
（1）中心控制级设备 中心控制级设备作用主要是： ?监控整个系统的运行。 ?协调区域控制级的运行。 ?具备区域控制级的所有功能。（2）区域控制级设备 区域控制级设备作用主要是： ?监控受控区域的运行。
?对路口交通信号进行协调控制。 ?对路口交通信号机的工作状态和故障情况进行监视。 ?通过人机回话对路口交通信号机进行人工干预。 ?监视和控制区域级外部设备的运行。 ?进行交通流量统计处理。 （3）路口控制级设备 路口控制级设备即信号机，其作用主要是： ?控制路口交通信号灯。 ?接收处理来自车辆检测器的交通流信息，并定时向区域计算机发送。 ?接收处理来自区域计算机的命令，并向区域计算机反馈工作状态和故障信息。 ?具有单点优化能力。 3系统功能设计 3.1基础功能 （1）区域自适应控制 系统以控制子区作为基本控制单元，综合考虑子区内的交通运行状态（如交通阻塞、交通拥挤、交通顺畅）、交叉口的关联性大小、交叉口的实际交通量，确定公共信号周期与相位差的决策模型，并运用智能优化算法实时优化子区协调控制配时参数，实现控制子区交叉口的协调控制功能。 系统的区域交叉口协调控制能够确保控制区域内的交通流时刻处于最佳运行状态，相邻交叉口之间协调方向的行驶车流可以获得尽可能不停顿的通行权，大大降低车辆在交叉口频繁加减速所产生的交通污染，减少区域交通总的车辆燃油
原语音信号 [y,fs,bits]=wavread('C:\Users\Administrator\Desktop\111.wav'); >> sound(y,fs,bits); >> n=length(y) n = 92611 >> Y=fft(y,n); >> subplot(2,1,1);plot(y); >> subplot(2,1,2);plot(abs(Y));
加噪声 >> [y,fs,bits]=wavread('C:\Users\Administrator\Desktop\111.wav'); >> sound(y,fs,bits); >> n=length(y) n = 92611 >> Noise=0.2*randn(n,2); >> s=y+Noise; >> sound(s) >> subplot(2,1,1); >> plot(s) >> S=fft(s); >> subplot(2,1,2); >> plot(abs(S)) >> title('加噪语音信号的频谱波形')
FIR 低通滤波器 fp=1000;fc=1200;As=100;Ap=1;fs=30000; >> wc=2*fc/fs;wp=2*fp/fs; >> N=ceil((As-7.95)/(14.36*(wc-wp)/2))+1; >> beta=0.1102*(As-8.7); >> Win=Kaiser(N+1,beta); b=fir1(N,wc,Win); >> freqz(b,1,512,fs); >> s_low=filter(b,1,s); >> plot(s_low);title('信号经过低通滤波器的时域图') >> S_low=fft(s_low,n); >> plot(abs(S_low));title('信号经过低通滤波的频谱') >> sound(s_low,fs,bits)
数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名： 学号： 组员：
摘要 滤波器设计在数字信号处理中占有极其重要的地位，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)
信号处理实验 实验八：音频频谱分析仪设计与实现
一、实验名称：音频频谱分析仪设计与实现 二、实验原理： MATLAB是一个数据信息和处理功能十分强大的工程实用软件，其数据采集工具箱为实现数据的输入和输出提供了十分方便的函数和命令。本实验可以用MATLAB进行音频信号频谱分析仪的设计与实现。 1、信号频率、幅值和相位估计 (1)频率(周期)检测 对周期信号来说，可以用时域波形分析来确定信号的周期，也就是计算相邻的两个信号波峰的时间差、或过零点的时间差。这里采用过零点(ti)的时间差T(周期)。频率即为f = 1/T，由于能够求得多个T值(ti有多个)，故采用它们的平均值作为周期的估计值。 (2)幅值检测 在一个周期内，求出信号最大值ymax与最小值ymin的差的一半，即A = (ymax - ymin)/2，同样，也会求出多个A值，但第1个A值对应的ymax和ymin不是在一个周期内搜索得到的，故以除第1个以外的A值的平均作为幅值的估计值。 (3)相位检测 采用过零法，即通过判断与同频零相位信号过零点时刻，计算其时间差，然后换成相应的相位差。φ=2π(1-ti/T)，{x}表示x的小数部分，同样，以φ的平均值作为相位的估计值。 频率、幅值和相位估计的流程如图所示。
其中tin表示第n个过零点，yi为第i个采样点的值，Fs为采样频率。 2、数字信号统计量估计 (1) 峰值P的估计 在样本数据x中找出最大值与最小值，其差值为双峰值，双峰值的一半即为峰值。P=0.5[max(yi)-min(yi)] (2)均值估计 式中，N为样本容量，下同。 (3) 均方值估计