Protues8255实验报告

华北电力大学

实验报告|

|

实验名称8255应用实验

课程名称微机原理及应用

老师

| 专业班级

姓名

学号

一．实验要求

1、当K0为高电平时，LED由高到低循环亮灭；当K0为低电平时，LED由低到高循环亮灭；当K1为高电平时，LED全灭且程序结束

（扩展：当K1为高电平时，LED全灭但程序不结束，等K1再次为低电平时灯又开始循环）。

2、在8259的IR2端输入中断请求信号，该信号由8253的方波信号产生（频率1Hz）。每来一个上升沿，申请中断一次，CPU响应后通过8255 PA口使发光二极管亮，第1次中断，LED0亮，第2次中断，LED1亮，……第8次中断，LED7亮，中断8次后结束。

3、在8259的IR2端输入中断请求信号，该信号由8253的方波信号产生（频率1Hz）。每来一个上升沿，申请中断一次，CPU响应后通过8255 PA口使7段数码管亮，第1次中断，显示0，第2次中断，显示1，…… 第10次中断，显示9，程序结束。

【8255的片选地址为8000H,8259的片选地址为9000H，8253的片选地址为A000H】

二思路

1.硬件

（1）74HC138译码电路如图所示，A15为1，E2、E3接地保证74HC138正常工作，此时Y0、Y1、Y2对应地址分别为8000H、9000H、0A000H. 此芯片实验1.2.3都需用到

（2）8259的片选地址为9000H,所以CS接Y1

8086有16位数据总线，其低8位作为偶存储体来传输数据，8086的A0要一直为0，所以8259的A0要接8086的A1

IR2端输入中断请求信号，该信号由8253的方波信号产生（频率1Hz）,所以8259的IR2端与8253的OUT1端相连。

此芯片实验2,3要用到

（3）8253的片选地址为A000H，所以CS接Y2。

CLK0接入1MHz信号，GATE0与GATE1同时接电源，OUT0输入到CLK1,OUT1经分频输出1HZ的信号。

此芯片实验2,3要用到

（4）对于实验1,8255的连接如下图

8255的片选地址为8000H,所以CS接Y0。

端口A中的两个端口接开关，作为输入（开关弹起时为高电平，按下时为低电平）；端口B中的全部端口作为输出显示。

（5）对于实验2，8255的连接如下图

端口A中的全部端口作为输出显示

（6）对于实验3，8255的连接如下图

端口A中的全部端口连接一二级数码管，且该数码管为共阴极连接。

二．软件

(1)实验1流程图

（2）实验2,3流程图

（3）8255初始化（实验1.2.3都需要）

方式选择控制字，对于实验1：A口方式0输入,B口方式0输出（10010000B=90H）

对于实验2,3：A口方式0输出 (10000000B=80H) 程序如下：

实验1：

MOV DX,8006H ;控制器地址为8006H

MOV AL,90H ;A口方式0输入,B口方式0输出

（10010000B=90H）

OUT DX,AL

实验2,3：

MOV DX,8006H ;控制器地址为8006H

MOV AL,80H ;A口方式0输出(10000000B=80H)

OUT DX,AL

（4）8253及8259的初始化（实验2,3需要）

关于8253及8259的设置与前面的两个实验均相同，其程序如下所示;8253A初始化

MOV DX,0A006H ;控制器地址0A006H

MOV AL,37H ;计数器0选用工作方式3(方波) OUT DX,AL ;写计数器0工作方式

MOV DX,0A000H ;计数器0地址0A000H

MOV AL,00H ;计数器0预置值为1000H

OUT DX,AL

MOV AL,10H

OUT DX,AL

MOV DX,0A006H ;控制器地址0A006H

MOV AL,75H ;计数器1选用工作方式2(分频) OUT DX,AL ;写计数器1工作方式

MOV DX,0A002H ;计数器1地址0A002H

MOV AL,00H ;计数器1预置值1000H

OUT DX,AL

MOV AL,10H

OUT DX,AL

;8259A初始化

MOV DX,9000H ;ICW1=13H(0001 0011B) MOV AL,13H

OUT DX,AL

MOV DX,9002H ;ICW2=80H(1000 0000B) MOV AL,80H

OUT DX,AL

MOV DX,9002H ;ICW3=01H(0000 0001B) MOV AL,01H

OUT DX,AL

（5）关于实验1的设置

读A口信息，判断K0K1

若K1=1,则灯全灭

若K0K1=00,则灯由低到高循环亮灭

若K0K1=10,则灯由高到低循环亮灭

其中还需要两个延时子程序

程序如下

MOV BL,01H

IN AL,DX ;读A口中的数据

AND AL,03H ;AL=AL&00000011B

CMP AL,02H

JAE MIE ;K1=1转到MIE

CMP AL,01H

JE NEXT1 ;K1K0=01B转HL,K1K0=00B做LH

NEXT1:MOV DX,8002H ;B口地址为8002H

MOV AL,BL

OUT DX,AL

CALL DELAY1

ROL BL,1

JMP NEXT1

NEXT2:MOV DX,8002H ;B口地址为8002H

MOV AL,BL

OUT DX,AL

CALL DELAY2

ROR BL,1

JMP NEXT2

MOV AL,00H ;LED全灭

OUT DX,AL

MOV DX,8000H ;A口地址为8000H IN AL,DX ;读A口中的数据

AND AL,03H ;AL=AL&00000011B

CMP AL,01H

JE NEXT2 ;K1K0=01B转到NEXT2 CMP AL,00H

JE NEXT1 ;K1K0=00B转到NEXT1 JMP MIE

延时子程序如下：

DELAY1 PROC NEAR

MOV CX,0A000H

L2: MOV DX,8000H ;A口地址为8000H IN AL,DX ;读A口中的数据

AND AL,03H ;AL=AL&00000011B

CMP AL,02H

JAE MIE ;K1=1转到MIE

CMP AL,01H

JE NEXT2 ;K1K0=01B转到NEXT2

LOOP L2

RET

DELAY1 ENDP

DELAY2 PROC NEAR

MOV CX,0A000H

L4: MOV DX,8000H ;A口地址为8000H

IN AL,DX ;读A口中的数据

AND AL,03H ;AL=AL&00000011B

CMP AL,02H

JAE MIE ;K1=1转到MIE

CMP AL,00H

JE NEXT1 ;K1K0=00B转到NEXT1

LOOP L4

RET

DELAY2 ENDP

（6）关于中断的设置（实验2,3需要）(其中红色字体为实验2和实验3不同的部份)

【需注意，对于实验3要编写数据段，为了能直接调用共阴极数码管段码表】程序如下：

DATA SEGMENT ;"DATA"逻辑段开始

TABLE DB 3FH,06H,5BH,4FH,66H,6DH,7DH,07H,7FH,6FH ;共阴极数码管段码表[0-9]

;非规范装入中断向量

MOV AX,0

MOV DS,AX

MOV BX,82H*4

MOV AX,OFFSET INTP

MOV [BX],AX

INC BX

INC BX

MOV AX,SEG INTP

MOV [BX],AX

MOV DX,9002H ;允许IRQ2中断

IN AL,DX

AND AL,0FBH

OUT DX,AL

;主程序(高电平点灯)

MOV DX,8000H

MOV AL,00H ;LED灯初始化全灭

OUT DX,AL

MOV BL,1 ;第一次中断后LED0亮

【MOV AX,DATA

MOV DS,AX

LEA SI,TABLE ;第一次中断的状态:显示'0'】

MOV CX,8 ;设置中断次数为8

【实验3设置的中断次数为10】

STI

L1:MOV AL,82H ; Proteus中8086模型有问题，它取得的中断号是最后发到总线上的数据，并不是由8259发出的中断号

MOV DX,0B000H

OUT DX,AL ;所以造成了要在这里执行EOI 的假相，地址使用没有用过的，如b000h，8259的地址是9000h CMP CX,0

JZ QUIT

JMP L1 ;等待中断

QUIT:CLI ;关中断

MOV DX,9002H ;禁止IRQ2中断

IN AL,DX

OR AL,4

OUT DX,AL

STI ;开中断

MOV CX,0FFFFH

L: LOOP L ;为看到第8次中断结束后LED灯亮,需用延时

MOV DX,8000H

MOV AL,00H ;8次中断结束，LED灯全灭

OUT DX,AL

JMP $

INTP PROC ;中断服务子程序

MOV AL,BL [ MOV AL,[SI]]

MOV DX,8000H

OUT DX,AL

ROL BL,1 [INC SI]

DEC CX

MOV DX,9000H

MOV AL,20H ;发中断结束命令EOI

OUT DX,AL

IRET ;中断返回

三．运行结果

实验1：当K1高电平时，LED灯都不亮。

当K1低电平，K0高电平时，LED灯按如下次序亮

当K1为高电平，K0为低电平时，LED灯按如下顺序亮

实验2. 第1次中断，LED0亮。以后每来一个上升沿，LED灯依次亮灭，直到计数到8为止。

最后LED灯全灭。

实验3：第一次中断数码管上显示数字“0”，以后每次中断，数码管上依次显示数字“1”—“9”。

数字信号处理实验报告

一、实验名称：基本信号的产生 二、实验目的：I 利用MATLAB 产生连续信号并作图 II 利用MATLAB 产生离散序列并作图 III 利用MATLAB 进行噪声处理 三、 实验内容： I 利用MATLAB 产生下列连续信号并作图 ①X(t)=-2u(t-1),-1=0); plot(t,x); 图形如右： ② X(t)=-(e^-0.1t)*sin(2/3*t),0

-1.5-1 -0.5 0.5 1 1.5 2 II 利用MATLAB 产生下列离散序列并作图 ① X(t)=1,-5<=t<=5 else 0,-15<=t<=15 MATLAB 程序如下: k= -15: 15; x=[zeros(1,10),ones(1,11),zeros(1,10)]; stem(k,x) 图形如下： ② X(t)=0.9^k*(cos(0.25*pi*k)+sin(0.25*pi*p),-20

《数字信号处理》实验报告

数字信号处理》 实验报告 年级：2011 级班级：信通 4 班姓名：朱明贵学号： 111100443 老师：李娟 福州大学 2013 年11 月

实验一快速傅里叶变换(FFT)及其应用 一、实验目的 1. 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB^的有关函数。 2. 熟悉应用FFT对典型信号进行频谱分析的方法。 3. 了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。 4. 熟悉应用FFT实现两个序列的线性卷积和相关的方法。 二、实验类型 演示型 三、实验仪器 装有MATLA爵言的计算机 四、实验原理 在各种信号序列中，有限长序列信号处理占有很重要地位，对有限长序列，我们可以 使用离散Fouier变换(DFT)。这一变换不但可以很好的反映序列的频谱特性，而且易于用快速算法在计算机上实现，当序列x(n)的长度为N时，它的DFT定义为： JV-1 $生 反变换为： 如-器冃吋 科— 有限长序列的DFT是其Z变换在单位圆上的等距采样，或者说是序列Fourier变换的等 距采样，因此可以用于序列的谱分析。 FFT并不是与DFT不同的另一种变换，而是为了减少DFT运算次数的一种快速算法。它 是对变换式进行一次次分解，使其成为若干小点数的组合，从而减少运算量。常用的FFT 是以2为基数的，其长度A - o它的效率高，程序简单，使用非常方便，当要变换的 序列长度不等于2的整数次方时，为了使用以2为基数的FFT，可以用末位补零的方法，使其长度延长至2的整数次方。 (一)在运用DFT进行频谱分析的过程中可能的产生三种误差 1 .混叠 序列的频谱是被采样信号频谱的周期延拓，当采样速率不满足Nyquist定理时，就会 发生频谱混叠，使得采样后的信号序列频谱不能真实的反映原信号的频谱。避免混叠现象的 唯一方法是保证采样速率足够高，使频谱混叠现象不致出现，即在确定采样频率之前，必须

数字信号处理-实验报告

学生实验报告 （理工类） 课程名称：数字信号处理专业班级：通信（4）班学生学号：学生姓名： 所属院部：网络与通信工程学院指导教师： 20 16 ——20 17 学年第一学期 金陵科技学院教务处制

实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求 实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。

实验项目名称：MATLAB语言工作环境和基本操作实验学时： 同组学生姓名：实验地点：工科楼A205 实验日期：实验成绩： 批改教师：批改时间： 一、实验目的和要求 目的： 1．初步了解MATLAB开发环境和常用菜单的使用方法； 2．熟悉MATLAB常用窗口，包括命令窗口、历史窗口、当前工作窗口、工作空间浏览器窗口、数组编辑器窗口和M文件编辑/调试窗口等； 3．了解MATLAB的命令格式； 4．熟悉MATLAB的帮助系统。 要求： 1. 简述实验原理及目的。 2. 记录调试运行情况及所遇问题的解决方法。 3. 简要回答思考题。 二、实验仪器和设备 微型计算机、Matlab6.5以上版本的编程环境。 三、实验过程 命令窗口(Command Window)： (1) 用于执行MATLAB命令，正常情况下提示符为“>>”，表示MATLAB进入工作状态。 (2) 在提示符后输入运算指令和函数调用等命令（不带“；”），MATLAB将迅速显示出结果并 再次进入准备工作状态。 (3) 若命令后带有“；”，MATLAB执行命令后不显示结果。 (4) 在准备工作状态下，如果按上下键，MATLAB会按顺序依次显示以前输入的命令，若要执 行它，则直接回车即可。 工作空间(Workspace)： (1) 显示计算机内存中现有变量的名称、类型、结构及其占用子节数等。 (2) 如果直接双击某变量，则弹出Array Editor窗口供用户查看及修改变量内容。 (3) 该窗口上有工具条支持用户将某变量存储到文件中或者从文件中载入某变量。 命令历史记录(Command History)： (1) 保存并显示用户在命令窗口中输入过的命令，以及每次启动MATLAB的时间等信息 (2) 若双击某条命令记录，则MATLAB会再次执行该命令。 当前路径窗口(Current Directory)：

信号处理实验报告、

第一题 如何用计算机模拟一个随机事件，并估计随机事件发生的概率以计算圆周率π。 解： （一）蒙特卡洛方法可用于近似计算圆周率：让计算机每次随机生成两个0到1之间的数，看以这两个实数为横纵坐标的点是否在单位圆内。生成一系列随机点，统计单位圆内的点数与总点数，（圆面积和外切正方形面积之比为π：4），当随机点取得越多时，其结果越接近于圆周率。 代码： N=100000000; x=rand(N,1); y=rand(N,1); count=0; for i=1:N if (x(i)^2+y(i)^2<=1) count=count+1; end end PI=vpa(4*count/N,10) PI = 3.1420384

蒙特卡洛法实验结果与试验次数相关，试验次数增加，结果更接近理论值 （二）18世纪，法国数学家布丰和勒可莱尔提出的“投针问题”，记载于布丰1777年出版的著作中：“在平面上画有一组间距为d的平行线，将一根长度为l （l

数字信号处理期末综合实验报告

数字信号处理综合实验报告 实验题目：基于Matlab的语音信号去噪及仿真 专业名称： 学号： 姓名： 日期： 报告内容： 一、实验原理 1、去噪的原理 1.1 采样定理 在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中，最高频率fmax的2倍时，即：fs.max>=2fmax,则采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。1924年奈奎斯特(Nyquist)就推导出在理想低通信道的最高大码元传输速率的公式: 理想低通信道的最高大码元传输速率=2W*log2 N (其中W是理想低通信道的带宽,N是电平强度)为什么把采样频率设为8kHz?在数字通信中，根据采样定理, 最小采样频率为语音信号最高频率的

2倍 频带为F的连续信号f(t)可用一系列离散的采样值f(t1),f(t1±Δt)，f(t1±2Δt)，...来表示,只要这些采样点的时间间隔Δt≤1/2F，便可根据各采样值完全恢复原来的信号f(t)。这是时域采样定理的一种表述方式。 时域采样定理的另一种表述方式是：当时间信号函数f(t)的最高频率分量为fM时,f(t)的值可由一系列采样间隔小于或等于1/2fM的采样值来确定,即采样点的重复频率f≥2fM。图为模拟信号和采样样本的示意图。 时域采样定理是采样误差理论、随机变量采样理论和多变量采样理论的基础。对于时间上受限制的连续信号f(t)（即当│t│>T时,f(t)=0,这里T=T2-T1是信号的持续时间），若其频谱为F（ω）,则可在频域上用一系列离散的采样值 （1-1） 采样值来表示,只要这些采样点的频率间隔 （1-2） 。 1.2 采样频率 采样频率，也称为采样速度或者采样率，定义了每秒从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。采样频率的倒数是采样周期或者叫作采样时间，它是采样之间的时间间隔。通俗的讲采样频率是指计算

视频信号处理实验报告

中南大学 实验报告（实验一） 实验名称 JM代码编译与编解码参数配置 课程名称视频信号处理 姓名：杨慧成绩:__________________ 班级：电子信息工程1301班学号： 0903130117 日期： 2016.6.10 地点：综合实验楼 备注：

1.实验目的 1）掌握常用的编解码器参数及其用法，实现测试序列的编解码 2）初步了解H.264视频编解码的基本原理、熟开发工具的使用 3）学会使用相关的开发工具修改、调试参考软件,掌握使用相应软件实现视频编解码的经验与技巧，锻炼提高分析问题和解决问题的能力 4）调试、编译好相应的实验程序，正确配置测试参数，能预计可能出现的结果2.实验环境（软件、硬件及条件） Windows 7 3.实验方法 1）JM工作目录与文件设置 ①下载并解压JM源代码。 ②在源代码根目录下的bin文件夹中新建backup文件夹，将bin文件夹中所有文件移入该文件夹做备份。 ③在源代码根目录下新建encodtest文件夹，作为编码使用。将编码过程所需要的文件，例如：编码配置文件（encoder_baseline.cfg）、待编码视频序列文件（foreman_part_qcif.yuv,对应为编码配置文件中InputFile参数的值）复制到该文件夹中。 ④在源代码根目录下新建decodtest文件夹，作为解码使用。将解码过程所需要的文件，例如：解码配置文件（decoder.cfg）复制到该文件夹中。 ⑤检查实验用机安装的MS Visual C++版本，根据表3，本实验打开jm_vc10.sln 解决方案。

2）配置、编译、测试编码项目——lencod ①选中lencod项目，打开主菜单“项目——属性”，将所有配置（Debug、Release）和所有平台（Win32、x64）“常规”选项中的“输出目录”设置为 “.\bin\$(Configuration)_$(Platform)\”;将“调试”选项中“工作目录”设置为“.\encodtest”,在“命令参数”中设置要使用的解码配置文件，例如：“-d encoder_baseline.cfg”，然后确定修改。 ②选中lencod工程，选择鼠标右键菜单“设为启动项目”。 ③打开主菜单“生成--批生成”，勾选所有的lencod项目，点击生成后，将会在主目录bin文件夹的Debug_Win32/x64文件夹及Release_Win32/x64文件夹下生成Win32/x64平台的调试版（运行速度慢）和发行版（运行速度快）编码器程序lencod.exe。打开主菜单“生成--配置管理器”，将活动解决方案配置和平台分别设置为Release何Win32，执行调试完成编码。此时会在源代码根目录下的encodtest文件夹中生成几个新文件，其中test.264（对应编码配置文件中OutputFile参数的值）即为压缩码流文件。 3）配置、编译、测试解码项目--ldecod ①选中ldecod项目，打开主菜单“项目——属性”，将所有配置（Debug、Release）和所有平台（Win32、x64）“常规”选项中的“输出目录”设置为 “.\bin\$(Configuration)_$(Platform)\”;将“调试”选项中“工作目录”设置为“.\decodtest”,在“命令参数”中设置要使用的解码配置文件，例如：“ decoder.cfg”，然后确定修改。 ②将编码生成的压缩码流文件test.24复制到decodtest文件夹中。 ③选中lencod工程，选择鼠标右键菜单“设为启动项目”。 ④打开主菜单“生成--批生成”，勾选所有的ldecod项目，点击生成后，将会在主目录bin文件夹的Debug_Win32/x64文件夹及Release_Win32/x64文件夹下生成Win32/x64平台的调试版（运行速度慢）和发行版（运行速度快）编码器程序ldecod.exe。打开主菜单“生成--配置管理器”，将活动解决方案配置和平台分别设置为Release何Win32，执行调试完成编码。此时会在源代码根目录下的decodtest文件夹中生成几个新文件，其中test_dec.yuv（对应解码配置文

哈尔滨工程大学 语音信号处理实验报告

实 验 报 告 实验课程名称： 语音信号处理实验 姓名： 班级： 20120811 学号： 指导教师 张磊 实验教室 21B#293 实验时间 2015年4月12日 实验成绩 实验序号 实验名称 实验过程 实验结果 实验成绩 实验一 语音信号的端点检测 实验二 语音信号的特征提取 实验三 语音信号的基频提取

实验一 语音信号的端点检测 一、实验目的 1、掌握短时能量的求解方法 2、掌握短时平均过零率的求解方法 3、掌握利用短时平均过零率和短时能量等特征，对输入的语音信号进行端点检测。 二、实验设备 HP 计算机、Matlab 软件 三、实验原理 1、短时能量 语音信号的短时能量分析给出了反应这些幅度变化的一个合适的描述方法。对于信号)}({n x ，短时能量的定义如下： ∑ ∑∞ -∞ =∞ -∞ =*=-= -= m m n n h n x m n h m x m n w m x E )()()()()]()([222 2、短时平均过零率 短时平均过零率是指每帧内信号通过零值的次数。对于连续语音信号，可以 考察其时域波形通过时间轴的情况。对于离散信号，实质上就是信号采样点符号变化的次数。过零率在一定程度上可以反映出频率的信息。短时平均过零率的公式为： ∑∑-+=∞ -∞=--= ---=1)] 1(sgn[)](sgn[2 1 ) ()]1(sgn[)](sgn[21N n n m w w m n m x m x m n w m x m x Z 其中，sgn[.]是符号函数，即 ? ? ?<-≥=0)(10)(1 )](sgn[n x n x n x

数字信号处理实验报告

语音信号的数字滤波 一、实验目的： 1、掌握使用FFT进行信号谱分析的方法 2、设计数字滤波器对指定的语音信号进行滤波处理 二、实验内容 设计数字滤波器滤除语音信号中的干扰（4 学时） 1、使用Matlab的fft函数对语音信号进行频谱分析，找出干扰信号的频谱； 2、设计数字滤波器滤除语音信号中的干扰分量，并进行播放对比。 三、实验原理 通过观察原语音信号的频谱，幅值特别大的地方即为噪声频谱分量，根据对称性，发现有四个频率的正弦波干扰，将它们分别滤掉即可。采用梳状滤波器，经过计算可知，梳状滤波器h[n]={1，A，1}的频响|H(w)|=|A+2cos(w)|，由需要滤掉的频率分量的频响w,即可得到A，进而得到滤波器的系统函数h[n]。而由于是在离散频域内进行滤波，所以令w=(2k*pi/N)即可。 对原信号和四次滤波后的信号分别进行FFT变换，可以得到它们的幅度相应。最后，将四次滤波后的声音信号输出。 四、matlab代码 clc;clear;close all; [audio_data,fs]=wavread('SunshineSquare.wav'); %读取未处理声音 sound(audio_data,fs); N = length(audio_data); K = 0:2/N:2*(N-1)/N; %K为频率采样点

%sound(audio_data,fs); %进行一次FFT变换 FFT_audio_data=fft(audio_data); mag_FFT_audio_data = abs(FFT_audio_data); %画图 figure(1) %原信号时域 subplot(2,1,1);plot(audio_data);grid; title('未滤波时原信号时域');xlabel('以1/fs为单位的时间');ylabel('采样值'); %FFT幅度相位 subplot(2,1,2);plot(K,mag_FFT_audio_data);grid; title('原信号幅度');xlabel('以pi为单位的频率');ylabel('幅度'); %构造h[n]={1，A，1}的梳状滤波器，计算A=2cosW，妻子W为要滤掉的频率%由原信号频谱可知要分四次滤波，滤掉频响中幅度大的频率分量 %第一次滤波 a = [1,0,0,0];%y[n]的系数 [temp,k]=max(FFT_audio_data); A1=-2*cos(2*pi*k/N); h1=[1,A1,1]; audio_data_h1 = filter(h1,a,audio_data); FFT_audio_data_h1=fft(audio_data_h1);

信号检测论有无法实验报告剖析

------------------------------------------------------------------------------- 实验报告信息栏 系别心理系年级 13级2班姓名魏晓芹同组成员杨思琪、张彤、韩永超 实验日期 2016年4月学号 120105510215 教师评定 ------------------------------------------------------------------------------- 信号检测论有无法实验报告 摘要本次实验采用信号检测论中的有无法，测定被试在不同先定概率下对呈现信号和刺激的击中率与虚报率，计算其辨别力d′和判定标准β，并绘制出ROC 曲线；检验信号呈现的先定概率发生变化时，被试的击中率、虚报率、辨别力d′和判定标准β是否会受到影响。结果显示：（1）被试在先定概率为0.2、0.5、0.8的条件下，击中率分别为0.8、0.92、0.8625，虚报率分别为0.5125、0.56、0.75，辨别力d′分别为0.592、1.254、0.406，判定标准β分别为0.70、0.38、0.71。 关键词信号检测论；有无法；先定概率；辨别力d′；判定标准β 1引言 传统心理物理学对阈限的理解是有限的，不能将个体客观的感受性和主观的动机、反应偏好等加以区分，从而使研究者渐渐陷入到了由阈限概念本身所引发的僵局之中。而在1954年，坦纳和斯韦茨等人首次应用的信号检测论，正好解决了这个问题。 信号检测论的研究对象是信息传播系统中信号的接收问题。在心理学中，它是借助于数学的形式描述“接收者”在某一观察时间内将掺有噪音的信号从噪音中辨别出来。 信号检测论应用于心理学中的基本原理是：将人的感官、中枢分析综合过程看作是一个信息处理系统，应用信号检测论中的一些概念、原理对它进行分析。信号检测论在心理学中具体应用时，常把刺激变量当作信号，把对刺激变量起干扰作用的因素当作噪音，这样就可以把人接收外界刺激时的分辨问题等效于一个在噪音中检测信号的问题，从而便可以应用信号检测论来处理心理学中的实验结果。 信号检测论的理论基础是统计决策。信号检测论本身就是一个以统计判定为根据的理论。它的基本原理是：根据某一观察到的事件，从两个可选择的方面选

MATLAB实验报告(1-4)

信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn 为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。 举例：计算一个函数并绘制出在对应区间上对应的值。

2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表 示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1： 举例2：

工程信号处理实验报告

( 2011-2012 学年 第二学期) 重庆理工大学研究生课程论文 课程论文题目： 《工程信号处理实验报告》 课程名称 工程信号处理实验 课程类别 □学位课 非学位课 任课教师 谢明 所在学院 汽车学院 学科专业 机械设计及理念 姓名 李文中 学 号 50110802313 提交日期 2012年4月12日

工程信号处理实验报告 姓名：李文中学号：50110802313 实验报告一 实验名称：数据信号采集及采样参数选定 1实验目的 1.1了解信号采集系统的组成，初步掌握信号采集系统的使用。 1.2加深对采样定理的理解，掌握采样参数的选择方法 1.3了解信号采集在工程信号处理中的实际应用，及注意事项。 2 实验原理 2.1 模数转换及其控制 对模拟信号进行采集，就是将模拟信号转换为数字信号，即模/数（A/D）转换，然后送入计算机或专用设备进行处理。模数转换包括三个步骤：（1）采样，（2）量化，（3）编码。采样，是对已知的模拟信号按一定的间隔抽出一个样本数据。若间隔为一定时间 T，则称这种采样为等时间间隔采样。除特别注明外，一般都采用等时间间隔采样；量化，是一种用有限字长的数字量逼近模拟量的过程。编码，是将已经量化的数字量变为二进制数码，因为数字处理器只能接受有限长的二进制数。模拟信号经过这三步转换后，变成了时间上离散、幅值上量化的数字信号。A/D转换器是完成这三个步骤的主要器件。 在信号采集系统中，A/D 转换器与计算机联合使用完成模数转换。用计算机的时钟或用软件产生等间隔采样脉冲控制 A/D 转换器采样。A/D 转换器通过内部电路进行量化与编码，输出有限长的二进制代码。信号采集系统中，通常由以 A/D转换器为核心的接口电路及控制软件，进行信号采集控制。 *注这部分是由本实验所用的信号采集器自动完成的，以上也是实验器材-信号采集器的部分工作原理。以后实验中就不再赘述。 2.2 信号采集的参数选择

信号处理实验报告

数字信号处理 第四次实验报告 一、 实验目的 1.了解离散系统的零极点与系统因果性能和稳定性的关系 2.观察离散系统零极点对系统冲激响应的影响 3.熟悉MATLAB 中进行离散系统零极点分析的常用子函数 4.加深对离散系统的频率响应特性基本概念的理解 5.了解离散系统的零极点与频响特性之间的关系 6.熟悉MATLAB 中进行离散系统分析频响特性的常用子函数，掌握离散系统幅频响应和相频响应的求解方法。 二、实验过程 9.2已知离散时间系统函数分别为 ) 7.05.0)(7.05.0(3 .0)(1j z j z z z H ++-+-= )1)(1(3 .0)() 8.06.0)(8.06.0(3 .0)(32j z j z z z H j z j z z z H ++-+-= ++-+-= 求这些系统的零极点分布图以及系统的冲击响应，并判断系统因果稳定性。 %---------第一式-----------------------------------------------------------------------------% z1=[0.3,0]';p1=[-0.5+0.7j,-0.5-0.7j]';k=1; %z1零点向量矩阵，p1极点向量矩阵,k 系统增益系数---------------------------% [bl,al]=zp2tf(z1,p1,k); %将零极点增益函数转换为系统传递函数 subplot(3,2,1),zplane(bl,al); %zplane 显示离散系统的零极点分布图 ylabel('极点在单位圆内'); subplot(3,2,2),impz(bl,al,20); %impz 绘制系统的冲激响应图 %---------第二式-----------------------------------------------------------------------------% z2=[0,3,0]';p2=[-0.6+0.8j,-0.6-0.8j]'; %z2零点向量矩阵，p2极点向量矩阵---------------------------------------------------% [b2,a2]=zp2tf(z2,p2,k); %将零极点增益函数转换为系统传递函数 subplot(3,2,3),zplane(b2,a2); %zplane 显示离散系统的零极点分布图 ylabel('极点在单位圆上'); subplot(3,2,4),impz(b2,a2,20); %impz 绘制系统的冲激响应图 %---------第三式-----------------------------------------------------------------------------%

数字信号处理综合设计实验报告

数字信号处理实验八 调制解调系统的实现 一、实验目的: （1）深刻理解滤波器的设计指标及根据指标进行数字滤波器设计的过程（2）了解滤波器在通信系统中的应用 二、实验步骤: 1.通过SYSTEMVIEW软件设计与仿真工具，设计一个FIR数字带通滤波器，预先给定截止频率和在截止频率上的幅度值，通过软件设计完后，确认滤波器的阶数和系统函数，画出该滤波器的频率响应曲线，进行技术指标的验证。 建立一个两载波幅度调制与解调的通信系统，将该滤波器作为两个载波分别解调的关键部件，验证其带通的频率特性的有效性。系统框图如下： 规划整个系统，确定系统的采样频率、观测时间、细化并设计整个系统，仿真调整并不断改进达到正确调制、正确滤波、正确解调的目的。（参考文件

zhan3.svu） （1）检查滤波器的波特图，看是否达到预定要求； （2）检查幅度调制的波形以及相加后的信号的波形与频谱是否正常； （3）检查解调后的的基带信号是否正常，分析波形变形的原因和解决措施；（4）实验中必须体现带通滤波器的物理意义和在实际中的应用价值。 2.熟悉matlab中的仿真系统； 3.将1.中设计的SYSTEMVIEW（如zhan3.svu）系统移植到matlab中的仿真环境中，使其达到相同的效果； 4.或者不用仿真环境，编写程序实现该系统，并验证调制解调前后的信号是否一致。 实验总共提供三个单元的时间（6节课）给学生，由学生自行学习和自行设计与移植 三、系统设计 本系统是基于matlab的simulink仿真软件设计的基带信号调制与解调的系统，利用matlab自带的数字信号仿真模块构成其原理框图并通过设置载波、带通滤波器以及低通滤波器等把基带信号经过载波调制后再经乘法器、带通滤波器和低通滤波器等电路系统能解调出基带信号。 1、实验原理框图

数字信号处理实验报告

实验一：信号的表示 1.实现单位采样序列、单位阶跃序列、矩形序列程序及绘图1.1代码部分 subplot(3,1,1); n1=-5:10; y1=[zeros(1,5),1,zeros(1,10)]; stem(n1,y1) axis([-5,10,0,2]); title(' 单位采样序列 ') subplot(3,1,2); n2=-5:10; y2=[zeros(1,5),ones(1,5),zeros(1 ,6)]; stem(n2,y2) axis([-5,10,0,2]) title(' 矩形序列 ') subplot(3,1,3); n3=-5:10; y3=[zeros(1,5),ones(1,11)]; stem(n3,y3,'r') axis([-5,10,0,2]) title(' 单位阶跃序列 ') 1.2仿真结果 2.实现三角波、方波、锯齿波、sinc函数及绘图2.1代码部分 %三角波 subplot(4,1,1); x=0:0.001:0.05; y1=sawtooth(2*pi*50*x,0.5);

plot(x,y1) %锯齿波 subplot(4,1,2); x=0:0.001:0.05; y2=sawtooth(2*pi*50*x); plot(x,y2) %方波 subplot(4,1,3); x=0:0.001:0.05; y3=square(2*pi*50*x,50); plot(x,y3) %sinc函数 subplot(4,1,4); t=-5:0.1:5; y=sinc(t); plot(t,y); xlabel('时间t');ylabel('幅值A'); title('Sa函数') 2.2仿真结果 实验二：FFT频谱分析及应用 1.用FFT函数分析某信号的频率成分和功率谱密度并绘图1.1代码部分 t=0:0.001:0.8; x=sin(2*pi*50*t)+cos(2*pi*120*t) ; y=x+1.5*randn(1,length(t)); subplot(3,1,1); plot(t,x); subplot(3,1,2); plot(t,y); Y=fft(y,512); P=Y.*conj(Y)/512;

数字信号处理实验报告(同名22433)

《数字信号处理》 实验报告 课程名称：《数字信号处理》 学院：信息科学与工程学院 专业班级：通信1502班 学生姓名：侯子强 学号：0905140322 指导教师：李宏 2017年5月28日

实验一 离散时间信号和系统响应 一. 实验目的 1. 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解 2. 掌握时域离散系统的时域特性 3. 利用卷积方法观察分析系统的时域特性 4. 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号及系统响应进行频域分析 二、实验原理 1. 采样是连续信号数字化处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性的变化以及信号信息不丢失的条件，而且可以加深对离散傅里叶变换、Z 变换和序列傅里叶变换之间关系式的理解。 对连续信号()a x t 以T 为采样间隔进行时域等间隔理想采样，形成采样信号： ?()()()a a x t x t p t = 式中()p t 为周期冲激脉冲，$()a x t 为()a x t 的理想采样。 ()a x t 的傅里叶变换为μ ()a X j Ω： 上式表明将连续信号()a x t 采样后其频谱将变为周期的，周期为Ωs=2π/T 。也即采样信 号的频谱μ()a X j Ω是原连续信号xa(t)的频谱Xa(jΩ)在频率轴上以Ωs 为周期，周期延拓而成 的。因此，若对连续信号()a x t 进行采样，要保证采样频率fs ≥2fm ，fm 为信号的最高频率，才可能由采样信号无失真地恢复出原模拟信号 计算机实现时，利用计算机计算上式并不方便，因此我们利用采样序列的傅里叶变换来实现，即 ()() n P t t nT δ∞ =-∞ = -∑μ1()()*() 21 ()n a a a s X j X j P j X j jn T π∞ =-∞ Ω=ΩΩ= Ω-Ω∑μ()()|j a T X j X e ωω=ΩΩ=

数字信号处理实验报告 (实验四)

实验四 离散时间信号的DTFT 一、实验目的 １. 运用MA TLAB 计算离散时间系统的频率响应。 ２. 运用MA TLAB 验证离散时间傅立叶变换的性质。 二、实验原理 （一）、计算离散时间系统的DTFT 已知一个离散时间系统∑∑==-= -N k k N k k k n x b k n y a 00)()(，可以用MA TLAB 函数frequz 非常方便地在给定的L 个离散频率点l ωω=处进行计算。由于)(ωj e H 是ω的连续函数，需要 尽可能大地选取L 的值（因为严格说，在MA TLAB 中不使用symbolic 工具箱是不能分析模拟信号的，但是当采样时间间隔充分小的时候，可产生平滑的图形），以使得命令plot 产生的图形和真实离散时间傅立叶变换的图形尽可能一致。在MA TLAB 中，freqz 计算出序列{M b b b ,,,10 }和{N a a a ,,,10 }的L 点离散傅立叶变换，然后对其离散傅立叶变换值相除 得到L l e H l j ,,2,1),( =ω。为了更加方便快速地运算，应将L 的值选为2的幂，如256或 者512。 例3.1 运用MA TLAB 画出以下系统的频率响应。 y(n)-0.6y(n-1)=2x(n)+x(n-1) 程序： clf; w=-4*pi:8*pi/511:4*pi; num=[2 1];den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1) plot(w/pi,real(h));grid title(‘H(e^{j\omega}的实部’)) xlabel(‘\omega/ \pi ’); ylabel(‘振幅’)； subplot(2,1,1) plot(w/pi,imag(h));grid title(‘H(e^{j\omega}的虚部’)) xlabel(‘\omega/ \pi ’); ylabel(‘振幅’)； （二）、离散时间傅立叶变换DTFT 的性质。 1．时移与频移 设 )]([)(n x FT e X j =ω, 那么

数字信号处理实验报告

南京信息工程大学数字信号处理实验报告 学院：电子与信息工程学院 班级：11通信1班 学号：20111334020 姓名：杨丹 指导教师：乔杰 2013/12/6

目录 实验一Matlab基本知识和信号处理工具箱 (3) 实验二离散信号的产生及分析 (5) 实验三离散傅立叶变换及分析 (8) 实验四IIR数字滤波器的分析与设计 (12) 实验五FIR数字滤波器的分析与设计 (14)

实验一 Matlab 基本知识和信号处理工具箱 一、实验目的 1、 了解Matlab 的基本操作 2、 了解Matlab 工具箱的函数 3、 正确使用Matlab 进行试验仿真 二、实验内容 （一）画出振荡曲线 和它的包络 程序如下： t=0:pi/20:4*pi; yy=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r'); hold on plot(t,yy,':b',t,-yy,':b'); grid on shg 实验结果： （二）sawtooth 产生锯齿波或三角波 在Matlab 命令窗口中输入以下命令，也可以用脚本文件来实现。 程序如下： t=-10:0.1:10; x=sawtooth(t); )3sin(3/t e y t -=3/t e yy -=

plot(t,x) grid on Shg 实验结果： 三、实验小结 在本次实验中，了解了Matlab的基本操作，基本了解了信号处理工具箱的使用方法，能够进行简单的波形绘制。对波形代码参量有了基本的识别能力。

实验二 离散信号的产生及分析 一、实验目的 利用函数、脚本文件或在主命令窗口中编写离散信号命令 二、实验内容 （一）生成离散信号 并画出复序列x(n)的实部，虚部，幅值和相位图。 程序如下： figure(1);clf n=-15:24;alpha=-0.3+j*0.5; x=exp(alpha.*n); subplot(221);stem(n,real(x));title('real part');xlabel('n'); subplot(222);stem(n,imag(x));title('imaginary part');xlabel('n'); subplot(223);stem(n,abs(x));title('magnitude part');xlabel('n'); subplot(224);stem(n,angle(x)*180/pi);title('phase part');xlabel('n'); 实验结果： （二）请给出离散系统 的单位抽样响应和单位阶跃响应 2415,)()5.03.0(<≤-=+-n e n x n j )2(2)()2(7.0)1(3.0)(-+=-+--n x n x n y n y n y

matlab实验报告

实验报告 2. The Branching statements 一、实验目的： 1.To grasp the use of the branching statements; 2.To grasp the top-down program design technique. 二、实验内容及要求： 1．实验内容： 1).编写 MATLAB 语句计算 y(t)的值 (Write the MATLAB program required to calculate y(t) from the equation) ???<+≥+-=0 530 53)(2 2t t t t t y 已知 t 从-5到 5 每隔0.5取一次值。运用循环和选择语句进行计算。 (for values of t between -5 and 5 in steps of 0.5. Use loops and branches to perform this calculation.) 2).用向量算法解决练习 1, 比较这两个方案的耗时。 （tic ，toc 的命令可以帮助你完成的时间计算，请使用'help'函数）。 Rewrite the program 1 using vectorization and compare the consuming time of these two programs. (tic, toc commands can help you to finish the time calculation, please use the …help ? function). 2.实验要求： 在报告中要体现top-down design technique, 对于 3 要写出完整的设计过程。 三、设计思路： 1．用循环和选择语句进行计算： 1).定义自变量t ：t=-5:0.5:5; 2).用循环语句实现对自变量的遍历。 3).用选择语句实现对自变量的判断，选择。 4).将选择语句置入循环语句中，则实现在遍历中对数据的选择，从而实现程序的功能。 2. 用向量法实现： 1).定义自变量t ：t=-5:0.5:5; 2).用 b=t>=0 语句，将t>=0得数据选择出，再通过向量运算y(b)=-3*t(b).^2 + 5; 得出结果。 3).用取反运算，选择出剩下的数据，在进行向量运算，得出结果。 四、实验程序和结果 1.实验程序 实验程序：创建m 文件：y_t.m

《语音信号处理》实验报告

盛年不重来，一日难再晨。及时宜自勉，岁月不待人。 中南大学 信息科学与工程学院 语音信号处理 实验报告 指导老师：覃爱娜 学生班级：信息0704 学生名称：阮光武 学生学好：0903070430 提交日期：2010年6月18日

实验一 语音波形文件的分析和读取 一、实验的任务、性质与目的 本实验是选修《语音信号处理》课的电子信息类专业学生的基础实验。通过实验： （1）掌握语音信号的基本特性理论：随机性，时变特性，短时平稳性，相关性等； （2）掌握语音信号的录入方式和*.WAV音波文件的存储结构； （3）使学生初步掌握语音信号处理的一般实验方法。 二、实验原理和步骤： WAV文件格式简介 WAV文件是多媒体中使用了声波文件的格式之一，它是以RIFF格式为标准。每个WAV文件的头四个字节就是“RIFF”。WAV文件由文件头和数据体两大部分组成，其中文件头又分为RIFF/WAV文件标识段和声音数据格式说明段两部分。常见的WAV声音文件有两种，分别对应于单声道（11.025KHz采样率、8Bit的采样值）和双声道（44.1KHz采样率、16Bit的采样值）。采样率是指声音信号在“模拟→数字”转换过程中，单位时间内采样的次数；采样值是指每一次采样周期内声音模拟信号的积分值。对于单声道声音文件，采样数据为8位的短整数（short int 00H-FFH）；而对于双声道立体声声音文件，每次采样数据为一个16位的整数（int），高八位和低八位分别代表左右两个声道。WAV文件数据块包含以脉冲编码调制（PCM）格式表示的样本。在单声道WAV文件中，道0代表左声道，声道1代表右声道；在多声道WAV文件中，样本是交替出现的。WAV文件的格式见表1。