数字信号处理课程设计三

课题一心电信号分析系统的设计

一、本课题的目的

本设计课题主要研究数字心电信号的初步分析方法及滤波器的应用。通过完成本课题的设计，拟主要达到以下几个目的：

（1）了解基于LabVIEW的虚拟仪器的特点和使用方法，熟悉采用LabVIEW进行仿真的方法。

（2）了解人体心电信号的时域特征和频谱特征。

（3）进一步了解数字信号的分析方法；

（4）通过应用具体的滤波器进一步加深对滤波器的理解。

（5）通过本课题的设计，培养学生运用所学知识分析和解决实际问题的能力。

二、课题任务

利用labVIEW设计一个基于虚拟仪器的简单的心电信号分析系统。对输入的原始心电信号，进行一定的数字信号处理，进行频谱分析。根据具体设计要求完成系统的程序编写、调试及功能测试。

（1）对原始数字心电信号进行读取，由数字信号数据绘制出其时域波形。

（2）对数字信号数据做一次线性插值，使其成为均匀数字信号，以便后面的信号分析。

（3）根据心电信号的频域特征（自己查阅相关资料），设计相应的低通和带通滤波器。

（4）编程绘制实现信号处理前后的频谱，做频谱分析，得出相关结论。

（5）对系统进行综合测试，整理数据，撰写设计报告。

三、主要设备和软件

（1）PC机一台。

（2）LabVIEW软件一套，要求最低版本8.20。

四、设计内容、步骤和要求

必做部分：

1. 利用labVIEW读取MIT-BIH数据库提供的数字心电信号，并还原实际波形

美国麻省理工学院提供的MIT-BIH数据库是一个权威性的国际心电图检测标准库，近年来应用广泛，为我国的医学工程界所重视。MIT-BIH数据库共有48个病例，每个病例数据长30min，总计约有116000多个心拍，包含有正常心拍和各种异常心拍，内容丰富完整。

为了读取简单方便，采用其txt格式的数据文件作为我们的原心电信号数据。利用labVIEW提供的文件I/O函数，读取txt数据文件中的信号，并且还原实际波形。

2.对原始心电信号做线性插值处理

由于原始心电信号数据不是通过等间隔采样得到的，也就是说原始的心电数据并不是均匀的，而用Matlab 中提供的数字滤波器处理数据时，要求数据是等间隔的。因此设计的系统首先应对原始心电信号做线性插值处理，使其变为等间隔的数字信号，否则直接处理后会出现偏差，根据心电信号的特点, 把时间分隔成0.001s 。添加的幅值点采用一次线性插值。对二维数据进行插值,相连幅值间数据的插值根据时间进行,运算公式如下:

1--=?i i t t t ，001.0/t N ?=，1--=?i i A A A ，001.01+=-j j t t ，

N A A A j j /1?+=- 其中i t 是第i 个数据时间点，A i 是与之对应的数据，N 是两数据之间需要的插值数，A ?是

需要插值的两点数据差，，

，，，，，，，1321321-==N j arraysize i 1111----==i j i j A A t t ，， 1=j i ，时数组j j A t ，依次排列，即得到了插值后等间隔的新数据。

3.根据心电信号的频域特征，设计相应的低通和带通滤波器

一般正常人的心电信号频率在0.7～100HZ 范围内，幅度为V 10μ(胎儿)～5mV (成人)。人体心电信号微弱，信噪比小，因此，在采集心电信号时，易受到仪器、人体活动等因素的影响，而且所采集的心电信号常伴有干扰。采集心电数据时，由于人的说话呼吸，常常会混有约为0.1Hz 到0.25Hz 频段的干扰，对于这些低频干扰，可以让信号通过一个高频滤波器，低截止频率设置为0.25,来滤除低频信号，对于高频信号干扰，可以让信号再通过一个低频滤波器，其中截止频率设置为99Hz 。也可以直接应用带通滤波器设计。

根据以上说明，利用labVIEW 中的信号处理函数设计相应滤波器，滤除数字信号中的干扰信号。

4.对处理前后的心电信号分别做频谱分析，分析结果

利用labVIEW 对处理前后的心电信号编程显示其频谱，分析比对滤波处理前后的频谱，得出结论。

如果分析频谱，滤波效果不明显，则需变动滤波器参数指标，重新设计滤波器。通过频谱分析，多次试验确定最合适的滤波器。

5.系统界面设计

综合前面几步，设计出一个完整的系统，并且本着简洁的原则，设计友好的人机交互界面。

选作部分：

1.三种滤波器设计

分别设计Butterworth 、Chebyshev 、Inverse Chebyshev 三种滤波器，并对滤波后的信号

分别做频谱分析，比较几种滤波器的差别。

2.设计50HZ工频陷波器

由于电子设备采集到的信号经常会混有电源线干扰。电源线干扰是以50 Hz为中心的窄带噪声，带宽小于1Hz。设计相应滤波器滤除电源线干扰，并对处理后的信号做频谱分析。

五、课程设计报告要求

（1）设计报告书包括内容：课程设计题目，课程设计目的和意义，设计方案，详细设计步骤，设计结果（原理图等），测试和仿真结果（图形或数据）及其分析，其它有明确要求的设计内容，结论，参考文献等。

（2）提交课程设计报告时应同时提交相关设计和仿真分析材料（框图、程序、结果等）的电子版。

六、参考文献

[1] 陈锡辉，张银鸿编著.LabVIEW 8.20程序设计从入门到精通[M].北京：清华大学出版社，2007.

[2] 丁玉美.数字信号处理（第二版）.西安电子科技大学出版社，2001

[3] 吴大正. 信号与线性系统分析（第四版）. 高等教育出版社，2005,8

[4] 谢嘉奎. 电子线路--线性部分(第四版). 高等教育出版社，2003,2

[5] 陈后金. 信号分析与处理实验. 高等教育出版社，2006，8

七、附录——设计原理

附录：设计原理

1.心电信号的读取

txt格式的数据文件内容及格式如图1-1所示（以100.txt为例）：

图1.1 txt格式心电数据文件

其中文件的第一列为采样时间，第二列是在以MLII这种导联方式所得到的采样数据，第三列式以V5这种导联方式所得到的采样数据，全文件记录了约为10s的心电数据，3600个采样数据，每一行数据之间用Tab符分隔。

由于数据文件中后两列数据是对同一种心电信号进行不同的导联方式所得到的采样数据，所以可以采用任意其中的一种采样数据（比如选择MLII），摒弃另外一种，即可完成对此心电信号的分析。全部的心电文件记录时间约为10s，共计12个左右周期的心电信号。

根据txt格式的数据文件的特点，利用labvIEW提供的I/O文件函数，在本课题中，主要是围绕LabVIEW中的read from spreadsheet file读表单文件函数来设计心电信号的读取部分的VI，并利用XY Graph来对数据做图形化显示。让心电数据文件中的第一列时间数据作为x轴，对应的MLII方式的幅值作为y轴，以此得到绘制的原心电波形。

图1.2 读表单文件函数VI

图1.3读表单文件函数使用举例

实际设计心电信号数据文件时需要注意：

（1）数据文件的前两行为解释说明文字，不是真正的信号数据，读取信号程序要能够自动忽略前两行文字，只读取真正的数字信号数据（严禁自己手动删除原心电数据文件中的前两行数据，必须通过程序来实现忽略前两行文字的目的）。

（2）利用数组函数分别将文件的前两列分别读入一个一维数组。labvIEW默认的从文本文件中读取的数据都是字符串，因此在使用心电信号数据前需要将其转换为数值才可以。注意：第一列时间数据均为0:00.007这种格式，因此需要将字符串0:00.007先转化为字符串0.007，即去除字符串中冒号（：）以前的部分，然后再将其转为数值。

（3）最后利用已经转为数值的分别代表心电信号时间和幅值的两个一维数组，图形化还原原始心电信号波形，在此推荐利用labvIEW中XY Graph。

2.心电信号的线性插值处理

根据上文中提到的插值公式，以此为原理，设计labvIEW程序，对心电信号数据做线性插值处理。插值完以后的数据应该是时间均匀的、以0.001秒为间隔的。

此步骤主要是基于labvIEW中的数组操作函数来实现，建议一定首先熟悉并掌握labvIEW中的所有数组操作函数的作用和操作方法（比如array size函数、index array函数、insert into array函数等）。

其中一种插值方法的思路是：第一步中读取的心电信号数据的时间数据和幅值数据分别存放在一个一维数组中。然后利用for循环结构把所有数据依次读取进来。判断时间数据数组中前后两个相邻的数据间隔是否为0.001s，如果是则判断下一对相邻两个数据；如果间隔大于0.001s则在一个CASE结构里面做插值处理。

注意对时间数据做插值的同时一定不要忘记对幅值数据同样做插值处理，时间数据和幅值数据一定是相互对应的。

3.设计相应的数字滤波器

原心电信号里面是包含有噪声的，因此需要对数字心电信号做一定滤波处理。

LabVIEW提供的IIR滤波器类型有Butterworth、Chebyshev、Inverse Chebyshev、Elliptic 和Bessel滤波器。它们都有各自的特点，用途也不尽相同。

LabVIEW还提供了高级IIR和FIR滤波器子面板。在高级面板中，滤波器的设计部分和执行部分是分开的。由于滤波器的设计很费时间，而滤波过程则很快。在含有循环结构的程序中，可以将滤波器的设计放在循环外，将设计好的滤波器参数传递到循环内，在循环内进行滤波，从而提高程序的运行效率。

Labview提供的滤波器函数面板面板位于Functions Palette的Signal Processing| Filters面板下，如图3.1所示。

图3.1 滤波器函数面板

选择合适的滤波器为心电信号设计一个低通和高通滤波器，或者带通滤波器。

4. 频谱分析

应该对线性插值后的心电信号和滤波处理后的心电信号做傅里叶变换，画出其频谱，比对前后差异，分析滤波器性能。

labvIEW中频域分析函数被划分为两个面板：Transforms面板实现的函数功能主要有傅立叶变换、Hilbert变换、小波变换、拉普拉斯变换等；Spectral Analysis面板包含的函数主要包括功率谱分析、联合时频分析等。

图4.1 transforms面板

图4.2 spectral analysis面板

5.低通滤波器和FFT举例

信号源由一个正弦信号与一个经过高通滤波的高频信号叠加而成。高通滤波器的截止频率为100Hz，即滤掉频率小于100 Hz的低频噪声分量。信号滤波器为Butterworth滤波器，截止频率设为30Hz，即滤掉频率大于30Hz的噪声分量。从图中可以清楚地看到滤波后的信号基本还原了正弦信号。

图4.3 低通滤波VI程序面板

图4.4 低通滤波VI前面板

数字信号处理课程设计报告

抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于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)

dsp课程设计实验报告

DSP 课程设计实验 一、语音信号的频谱分析： 要求首先画出语音信号的时域波形，然后对语音信号进行频谱分析。在MATLAB 中，可以利用函数fft 对信号进行快速傅立叶变换，得到信号的频谱特性，从而加深对频谱特性的理解。 其程序为： >> [y,fs,bits]=wavread('I:\',[1024 5120]); >> sound(y,fs,bits); >> Y=fft(y,4096); >> subplot(221);plot(y);title('原始信号波形'); | >> subplot(212);plot(abs(Y));title('原始信号频谱'); 程序运行结果为： 二、设计数字滤波器和画出频率响应： 根据语音信号的特点给出有关滤波器的性能指标： 低通滤波器性能指标，p f =1000Hz ，c f =1200Hz ，s A =100dB ，p A =1dB ； 高通滤波器性能指标，c f =4800Hz ，p f =5000Hz ，s A =100dB ，p A =1dB ； 带通滤波器性能指标，1p f =1200Hz ，2p f =3000Hz ，1c f =1000Hz ，2c f =3200Hz ，s A =100dB ， p A =1dB ；

】 要求学生首先用窗函数法设计上面要求的三种滤波器，在MATLAB中，可以利用函数firl 设计FIR滤波器；然后再用双线性变换法设计上面要求的三种滤波器，在MATLAB中，可以利用函数butte、cheby1和ellip设计IIR滤波器；最后，利用MATLAB中的函数freqz画出各种滤波器的频率响应，这里以低通滤波器为例来说明设计过程。 低通： 用窗函数法设计的低通滤波器的程序如下： >> fp=1000;fc=1200;As=100;Ap=1;fs=22050; >> wc=2*fc/fs;wp=2*fp/fs; >> N=ceil(/*(wc-wp)/2))+1; >> beta=*; >> Win=Kaiser(N+1,beta); 、 >>b=firl(N,wc,Win); >>freqz(b,1,512,fs); 程序运行结果： 这里选用凯泽窗设计，滤波器的幅度和相位响应满足设计指标，但滤波器长度（N=708）太长，实现起来很困难，主要原因是滤波器指标太苛刻，因此，一般不用窗函数法设计这种类型的滤波器。 用双线性变换法设计的低通滤波器的程序如下： >> fp=1000;fc=1200;As=100;Ap=1;fs=22050; >> wc=2*fc/fs;wp=2*fp/fs; 》 >> [n,wn]=ellipord(wp,wc,Ap,As); >> [b,a]=ellip(n,Ap,As,wn); >> freqz(b,a,512,fs); ^

DSP课程设计总结报告

课程设计总结报告课程名称DSP控制器及其应用 设计题目万年历设计 业专电子信息工程 班级 姓名 学号

指导教师 报告成绩 信息工程学院 年六月十三日二〇一四 录目 言前 (3) 设计要求第一章4.....................................................................................基本要求1.14.....................................................................................

系统的组成和工作原理第二章5............................................................. 芯片的工作原理VC5509APGE2.1DSPTMS3205.............................. 液晶显示器的工作原理2.2LCD16026..............................................主电路图及程序流程图第三章.. (7) 主电路图3.17...................................................................................... 程序总流程图3.27.............................................................................. 程序分块流程图3.38..........................................................................软件程序设计第四章9.............................................................................

数字信号处理课设+语音信号的数字滤波

语音信号的数字滤波 ——利用双线性变换法实现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段上，可以通过以下的正切变换实现

dsp课程设计实验报告总结

DSP课程设计总结（2013-2014学年第2学期） 题目： 专业班级：电子1103 学生姓名：万蒙 学号：11052304 指导教师： 设计成绩： 2014 年6 月

目录 一设计目的----------------------------------------------------------------------3 二系统分析----------------------------------------------------------------------3 三硬件设计 3.1 硬件总体结构-----------------------------------------------------------3 3.2 DSP模块设计-----------------------------------------------------------4 3.3 电源模块设计----------------------------------------------------------4 3.4 时钟模块设计----------------------------------------------------------5 3.5 存储器模块设计--------------------------------------------------------6 3.6 复位模块设计----------------------------------------------------------6 3.7 JTAG模块设计--------------------------------------------------------7 四软件设计 4.1 软件总体流程-----------------------------------------------------7 4.2 核心模块及实现代码---------------------------------------8 五课程设计总结-----------------------------------------------------14

DSP课程设计 (3)

深圳大学考试答题纸 (以论文、报告等形式考核专用) 二○一四～二○一五学年度第1 学期 课程编号 课程 名称 单片机/ARM /DSP技术实践 主讲 教师 评 分 学号姓名 专业年级 题 目： 基于DSP2812的课程设计

一、实验要求 由外接的信号发生器产生一正弦信号（电压范围：0~3V），通过DSP的AD功能对此正弦信号进行采集，通过DSP的SCI功能与PC机之间进行通信，把所采集的AD信号发送至PC机端，在超级终端上进行实时显示。 二、实验原理 2.1 ADC概述 ADC，即模/数转换器，将模拟量转换成数字量，提供给控制器使用。TMS320F2812片上有一个12位分辨率、具有流水线结构的模/数转换器，其机构框图如图1所示。其前端为2个8选1多路切换器和2路同时采样/保持器，构成16个模拟输入通道，模拟通道的切换由硬件自动控制，并将各模拟通道的转换结果顺序存入16个结果寄存器中。 图1 ADC机构框图 2.2 ADC模块特点 （1）带2个8选1多路切换器和双采样/保持器的12位的ADC，共有16个模拟输入通道； （2）模拟量输入范围：0.0V-3.0V；

（3）转换率：在25MHZ的ADC时钟下为80ns; (4)转换结果存储在16个结果存储器中； （5）转换结果=4095*（输入的模拟信号-ADCLO）/3； （6）多种A/D触发方式：软件启动、EVA和EVB； （7）灵活中断方式：可以在每次转换结束或每隔一次转换结束触发中断； 3．AD C转换步骤 （1）初始化DSP系统； （2）设置PIE中断矢量表； （3）初始化ADC模块； （4）将ADC中断的入口地址装入PIE中断矢量表中，开中断； （5）软件启动ADC转换； （6）等待ADC中断； （7）在ADC中断中读取ADC转换结果，软件启动下一次ADC中断。 三、实验实现 3.1硬件方案设计 本实验以TMS320F2812为核心控制部件，利用软件编程，通过ADC模块对试验箱上的信号发生器发出的正弦信号进行采集，由于试验箱上的信号发生器只能调节到2V，所以此次实验只针对2V的正弦信号，再通过串口线与PC机连接，将采集转换的数字信号传送到PC机端的串口助手，并还原成采集时的电压值。硬件框架图如图2所示。本次ADC采用SEED-DEC2812的AD接口的ADCINA6通道。 图2 硬件框架图

数字信号处理课设共18页文档

数字信号处理课程设计 姓名：刘倩 学号：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 算法的应用

DSP实验报告

电气信息工程学院 D S P技术与综合训练 实验报告 班级 08通信1W 姓名丁安华 学号 08313115 指导老师倪福银刘舒淇 2011年09 月

目录 实验一 LED演示 1.1.实验目的 -------------------------------------------------P2 1. 2.实验设备-------------------------------------------------P2 1. 3.实验原理-------------------------------------------------P2 1. 4.实验程序设计流程------------------------------------------P3 1. 5.实验程序编写----------------------------------------------P4 1. 6.实验步骤-------------------------------------------------P7 1. 7.实验结果与分析--------------------------------------------P7实验二键盘输入 2.1.实验目的 -------------------------------------------------P8 2.2.实验设备-------------------------------------------------P8 2. 3.实验原理-------------------------------------------------P8 2. 4.实验程序设计流程------------------------------------------P9 2. 5.实验程序编写----------------------------------------------P10 2. 6.实验步骤-------------------------------------------------P14 2. 7.实验结果与分析--------------------------------------------P14实验三液晶显示器控制显示 3.1.实验目的 -------------------------------------------------P15 3.2.实验设备-------------------------------------------------P15 3.3.实验原理-------------------------------------------------P15 3. 4.实验程序设计流程------------------------------------------P17 3. 5.实验程序编写----------------------------------------------P18 3. 6.实验步骤-------------------------------------------------P22 3. 7.实验结果与分析--------------------------------------------P23实验四有限冲激响应滤波器（FIR）算法 4.1.实验目的 -------------------------------------------------P23 4.2.实验设备-------------------------------------------------P23 4.3.实验原理-------------------------------------------------P24 4.4.实验程序设计流程------------------------------------------P25 4. 5.实验程序编写----------------------------------------------P25 4. 6.实验步骤-------------------------------------------------P27 4. 7.实验结果与分析--------------------------------------------P28

数字信号处理课程规划报告

数字信号处理课程设计报告《应用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('频域信号'); 运行程序得到的图形：

DSP实验报告

实验一 程序的控制与转移 一、实验目的 1、掌握条件算符的使用。 2、掌握循环操作指令（BNAZ ）和比较操作指令（CMPR ） 二、实验设备 计算机、ZY13DSP12BD 实验箱、5402EVM 板。 三、实验原理 程序控制指令主要包括分支转移、子程序调用、子程序返回、条件操作及循环操作等。通过传送控制到程序存储器的其他位置，转移会中断连续的指令流。转移会影响在PC 中产生和保护的程序地址。其中转移可以分为两种形式的，一种是有条件的，另一种是无条件的。 四、实验内容 编写程序，实现计算y= ∑=5 1 i i x 的值。 五、实验步骤 1、用仿真机将计算机与ZY13DSP12BD 实验箱连接好，并依次打开实验箱电源、仿真机电源，然后运行CCS 软件。 2、新建一个项目：点击Project －New ，将项目命名为example2，并将项目保存在自己定义的文件夹下。 3、新建一个源文件example2.asm 。将该文件添加到工程example2.pjt 中。 4、在工程管理器中双击example2.asm ，编写源程序： .tiltle ”example2.asm ” .mmregs STACK .usect ”STACK ”,10H ;堆栈的设置 .bss x,5 ;为变量分配6个字的存储空间 .bss y,1 .def start .data table: .word 10,20,3,4,5 ;x1,x2,x3,x4,x5 .text Start: STM #0,SWWWSR ;插入0个等待状态 STM #STACK+10H,sp ;设置堆栈指针 STM #x,AR1 ;AR1指向x RPT #4 ;下一条被重复执行5遍 MVPD table,*AR1+ ;把程序存储器中的数据传送到数据存储器 LD #0,A ;A 清零 CALL SUM ;调用求和函数 end: B end SUM: STM #x,AR3 ;AR3指向x STM #4,AR2 ;AR2=4 loop: ADD *AR3+,A ;*AR3+A-->A,然后AR3+ BANZ loop,*AR2- ;如果AR2的值不为0,则跳到loop 处;否则执行下一条指令 STL A,*(y) ;把A 的低16位赋给变量y

DSP课程设计报告

共享知识分享快乐 盛年不重来，一日难再晨。及时宜自勉，岁月不待人。 数据采集处理和控制系统设计 一课程设计要求 1.基本DSP硬件系统设计要求 ①基本DSP硬件系统以TMS320C54x系列为核心处理器，包括最小系统、存储器扩展、显示器、键盘、AD、DA等电路模块； ②硬件设计画出主要芯片及电路模块之间的连接即可，重点考查电路模块方案设计与系统地址分配； ③设计方案以电路示意图为主，辅以必要的文字说明。 2.基本软件设计要求 ①看懂所给例程，画出例程输出波形示意图； ②修改例程程序，使之输出其它波形，如方波、三角波、锯齿波等均可； ③设计方案以程序实现为主，辅以必要的文字说明。 3.课程设计报告要求 ①硬件系统设计：设计思路、设计系统功能、主要芯片选型及使用方法、设计方案说明、电路示意图 ②软件系统设计：示例程序功能解读及输出波形示意图、设计软件功能、设计思路、实现源码（带程序注释） ③报告总结 二系统分析 利用实验箱的模拟信号产生单元产生不同频率的信号，或者产生两个频率的信号的叠加。在DSP 中采集信号，并且对信号进行频谱分析，滤波等。通过键盘或者串口命令选择算法的功能，将计算的信号频率或者滤波后信号的频率在LCD 上显示。主要功能如下： (1)对外部输入的模拟信号采集到DSP 内存，会用CCS 软件显示采集的数据波形。 (2)对采集的数据进行如下算法分析： ①频谱分析：使用fft 算法计算信号的频率。 ②对信号进行IIR 滤波或FIR 滤波，并且计算滤波前后信号的频率。 ③外部键盘或者从计算机来的串口命令选择算法功能，并且将结果在 LCD 上显示。 绘制出DSP系统的功能框图、使用AD（Altium Designer）绘制出系统的原理图和PCB 版图。 在 DSP 中采集信号，用CCS 软件显示采集的数据波形，以及对采集的数据进行算法分析。 三硬件设计 3.1 硬件总体结构

数字信号处理课程设计

数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名： 学号： 组员：

摘要 滤波器设计在数字信号处理中占有极其重要的地位，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)

DSP技术与课程设计实验报告二(精)

东南大学自动化学院 实验报告 课程名称： D SP 原理及C 程序开发 第二次实验 实验名称：基于DSP 系统的实验——指示灯、拨码开关和定时器院（系）：自动化专业：自动化 姓名：学号： 实验室：实验组别： 同组人员：实验时间：2012 年 4 月 18日 评定成绩：审阅教师： 第一部分实验：基于DSP 系统的实验——指示灯和拨码开关 一．实验目的 1. 了解ICETEK –F28335-A 评估板在TMS320F28335DSP 外部扩展存储空间上的扩展。 2. 了解ICETEK –F28335-A 评估板上指示灯和拨码开关扩展原理。 3. 学习在C 语言中使用扩展的控制寄存器的方法。 二．实验设备 计算机，ICETEK –F28335-A 实验箱（或ICETEK 仿真器+ICETEK–F28335-A 评估板+相关连线及电源）。 三．实验原理

1．TMS320F28335DSP 的存储器扩展接口 存储器扩展接口是DSP 扩展片外资源的主要接口，它提供了一组控制信号和地址、数据线，可以扩展各类存储器和存储器、寄存器映射的外设。 -ICETEK –F28335-A 评估板在扩展接口上除了扩展了片外SRAM 外，还扩展了指示灯、DIP 开关和D/A 设备。具体扩展地址如下： 0x180004- 0x180005：D/A 转换控制寄存器 0x180001：板上DIP 开关控制寄存器 0x180000：板上指示灯控制寄存器 -与ICETEK –F28335-A 评估板连接的ICETEK-CTR 显示控制模块也使用扩展空间控制主要设备： 208000-208004h ：读-键盘扫描值，写-液晶控制寄存器 208002-208002h ：液晶辅助控制寄存器 208003-208004h ：液晶显示数据寄存器 2．指示灯与拨码开关扩展原理

dsp课程设计报告(2)分析

华北水利水电大学North China University of Water Resources and Electric Power DSP课程设计 题目： FIR数字低通滤波器 学院信息工程学院 专业电子信息工程 姓名 学号 指导教师

摘要 (1) 一. 绪论 (1) 1.1设计背景 (1) 1.2设计方法分析 (1) 二．FIR滤波器设计过程 (2) 2.1 FIR滤波器原理 (2) 2.2 FIR滤波器的实现方法 (3) 2.3 FIR滤波器的MATLAB实现 (4) 2.4 设计流程图 (6) 三．MATLAB和 CCS操作步骤及仿真结果 (7) 3.1 matlab中的.M文件的编写 (7) 3.2 工程文件的建立 (12) 3.3 仿真结果及分析 (12) 四．心得与总结 (12)

摘要 当前,数字信号处理技术受到了人们的广泛关注,其理论及算法随着计算机技术和微电子技术的发展得到了飞速地发展,并被广泛应用于语音和图象处理、数字通信、谱分析、模式识别和自动控制等领域。数字滤波器是数字信号处理中最重要的组成部分之一,几乎出现在所有的数字信号处理系统中。设计中通过MATLAB环境中图形化的方式建立数字信号处理的模型进行DSP的设计和仿真验证，将设计的图形文件.mdl直接转换成C语言程序在CCS中运行。利用MATLAB 软件开发产品加速了开发周期，比直接在CCs中编程方便快捷了很多，对于任何复杂功能的DSP系统，只需要进行少量的添加和修改就能完成功能正确的C语言程序设计。 一. 绪论 1.1设计背景 一个实际的应用系统中，由于设备或者是外界环境的原因，总存在各种干扰，使信号中混入噪声，譬如音频信号中高频成分的噪声使得音乐听起来刺耳，失去了原有悦耳的音质。为了提高信号质量，可以对信号进行滤波，从噪声中提取信号，即对一个具有噪声和信号的混合源进行采样，然后经过一个数字滤波器，滤除噪声，提取有用信号。DSP（数字信号处理器）与一般的微处理器相比有很大的区别，它所特有的系统结构、指令集合、数据流程方式为解决复杂的数字信号处理问题提供了便利，本文选用TMS320C54X作为DSP处理芯片，通过对其编程来实现FIR滤波器。对数字滤波器而言,从实现方法上,有FIR滤波器和无限冲激响应(IIR)滤波器之分。由于FIR滤波器只有零点,因此这一类系统不像IIR系统那样易取得比较好的通带与阻带衰减特性。但是FIR系统有自己突出的优点:①系统总是稳定的;②易实现线性相位;③允许设计多通带(阻带)滤波器。其中后两项是IIR系统不易实现的。 1.2设计方法分析 FIR滤波器的设计方法分析 数字滤波器依据冲激响应的宽度划分为有限冲激响应(FIR)滤波器和无限冲激响应滤波(IIR)。FIR 滤波器是有限长单位冲激响应滤波器，在结构上是非递归型的，有限冲激响应滤波器(FIR)，具有以下的优点：（1）可以在幅度特性随意设计的同时，保证精确、严格的线性相位；（2）由于FIR滤波器的单位脉冲响应h（n）是有限长序列，因此F I R 滤波器没有不稳定的问题；（3）由于FIR 滤

dsp实验报告 哈工大实验三 液晶显示器控制显示实验

实验三液晶显示器控制显示实验 一. 实验目的 通过实验学习使用2407ADSP 的扩展I/O 端口控制外围设备的方法，了解液晶显示器的显示控制原理及编程方法。 二. 实验设备 计算机，ICETEK-LF2407-EDU 实验箱。 三.实验原理 ICETEK-LF2407-A 是一块以TMS320LF2407ADSP 为核心的DSP 扩展评估板，它通过扩展接口与实验箱的显示/控制模块连接，可以控制其各种外围设备。 液晶显示模块的访问、控制是由2407ADSP 对扩展I/O 接口的操作完成。 控制I/O 口的寻址：命令控制I/O 接口的地址为0x8001，数据控制I/O 接口的地址为0x8003 和0x8004，辅助控制I/O 接口的地址为0x8002。 显示控制方法： ◆液晶显示模块中有两片显示缓冲存储器，分别对应屏幕显示的象素，向其中写入数 值将改变显示，写入“1”则显示一点，写入“0”则不显示。其地址与象素的对应 方式如下： ◆发送控制命令：向液晶显示模块发送控制命令的方法是通过向命令控制I/O 接口 写入命令控制字，然后再向辅助控制接口写入0。下面给出的是基本命令字、解释 和 C 语言控制语句举例。 ?显示开关：0x3f 打开显示；0x3e 关闭显示； ?设置显示起始行：0x0c0+起始行取值，其中起始行取值为0 至63； ?设置操作页：0x0b8+页号，其中页号取值为0-7； ?设置操作列：0x40+列号，其中列号为取值为0-63； ◆写显示数据：在使用命令控制字选择操作位置(页数、列数)之后，可以将待显示的 数据写入液晶显示模块的缓存。将数据发送到相应数据控制I/O 接口即可。

数字信号处理课程设计 1

（一）用窗函数法设计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;

DSP实验报告

学校代码学号分类号密级 DSP实验报告 院系名称 专业名称 年级 学生姓名 指导老师 年月日

实验一数据存储实验 一、实验目的 1. 掌握 TMS320C54X 程序空间的分配； 2. 掌握 TMS320C54X 数据空间的分配； 3. 能够熟练运用TMS320C54X 数据空间的指令。 二、实验设备 计算机，CCS 3.1版软件，DSP仿真器，E300实验箱，DSP-54XP CPU板。 三、实验系统相关资源介绍 1. 本实验指导书是以TMS320VC5416为例，介绍其相关的内部和外部存储器资源。对于其他类型的CPU请参考查阅相关的数据手册。） 下面给出TMS320VC5416的存储器分配表： 对于数据存储空间而言，映射表相对固定。值得注意的是内部寄存器都映射到数据存储器空间内。因此在编程时这些特定的空间不能作其他用途。 对于程序空间而言，其映射表和CPU 的工作模式有关。当MP/MC 引脚为高电平时，CPU 工作在微处理器模式；当MP/MC引脚为低电平时，CPU工作在微计算机模式。具体的MP和MC模式下的程序和数据映射关系如上图所示。 2. 样例程序实验操作简单说明： 本实验程序将对0x1000 开始的8 个地址空间，填写入0xAAAA 的数据，然后读出，并存储到以0x1008开始的8个地址空间，在CCS中可以观察DATA存储器空间地址0x1000～0x100F 值的变化。 四、实验步骤与内容 1. 在进行 DSP实验之前，需先连接好仿真器、实验箱及计算机，连接方法如下所示： 2. E300 底板的开关SW4 的第1位置ON，其余位置OFF，SW5全部置ON，其余开关不做设置要求。 3. 上电复位 在硬件安装完成后，确认安装正确、各实验部件及电源连接无误后，启动计算机，接通仿真器电源，此时，仿真器上的“红色指示灯”应点亮，否则DSP开发系统与计算机连接存

DSP课程设计报告—刘雅琪

一、课程设计的目的和要求 1.1课程设计目的： 本课程是DSP技术类课程配套的课程设计，要求学生通过高级语言或汇编语言编程实现较复杂的功能。通过课程设计，使学生加深对DSP芯片 TMS320C54的结构、工作原理的理解，获得DSP应用技术的实际训练，掌握设计较复杂DSP系统的基本方法。 1.2课程设计要求 1、认真查阅资料 2、课程设计前认真预习 3、遵守课程设计时间安排 4、认真保质保量完成设计要求 5、认真书写报告 二、系统功能介绍及总体设计方案 2．1 功能介绍 随着信息技术和计算机技术的飞速发展, 数字信号处理技术在众多领域得到广泛应用。数字滤波器由于其精度高、稳定性好、使用灵活等优点, 广泛应用在各种数字信号处理领域。数字滤波器根据冲击响应函数的时域特性, 可以分为FIR （有限长冲激响应滤波器）和IIR（无限长冲激响应滤波器）。FIR滤波器与IIR 滤波器相比,具有严格的线性相位,幅度特性可任意等优点。而且, FIR 滤波器的单位抽样响应是有限长的, 故一定是稳定的, 他又可以用快速傅里叶变换（FFT）算法来实现过滤信号，可大大提高运算效率。 本课程设计的是一个等波纹FIR 低通滤波器，其具体参数为：采样频率 F s=1000Hz，通带频率F pass=150Hz截止频率F sg=250Hz，通带衰减A pass =0.5dB 阻带衰减A stop=80dB。 2．2 总体设计方案： 先进行Matlab 程序设计产生待滤波数据（借助设计工具FDATOOL 产生设计系数），将其导入CCS在CCS!进行仿真调试运行，得到了输入和输出的波形及

其频谱。 三、主要设计内容和步骤 3?1 FIR 数字滤波器的原理分析 3.1.1FIR 数字滤波器 数字滤波器原理一般具有如下差分方程 N A N -1 y(n) f a k X( n- k) 、b k y( n-k) k z 0 k =0 式中x(n)为输入序列，y(n)为输出序列，兔和b k 为滤波器系数，N 是滤波器阶 N J 数。当所有的b k 均为零，则有 y (n )八 时(n-k) ⑵ k=0 (2) 式是FIR 滤波器的差分方程，其一般形式为 N -1 y(n)「h k x( n-k) (3) k ￡ 对(3)式进行z 变换，整理后可得FIR 滤波器的传递函数 FIR 的直接型结构: x(n) ——B —— h(0) n "(2) ' b ——■—— 1 h(N ￡3) ?_m h(N ￡2) I ---- ■—1 ? ■ ” ￡ ■ 1 h(N ￡l) z ￡i z ￡i y(n) H(z) = 迩 2 X o -k h)k)z

数字信号处理课程设计(对音乐信号的各种处理)

实验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)