信号处理实验 基本信号
实验结果如下所示：
1.]5[9.0][1-=n n x δ 1≤n ≤20
L=20;
nn=1:20;
imp=zeros(L,1);
imp(5)=0.9;
stem(nn,imp)
][8.0][2n n x δ= -15
≤n ≤15
L=31;
nn=-15:15;
imp=zeros(L,1);
imp(16)=0.8;
stem(nn,imp)
]333[5.1][3
-=n n x δ 300≤n ≤350 L=51;
nn=300:350;
imp=zeros(L,1);
imp(34)=1.5;
stem(nn,imp)
]7[45][4+=n n x δ -10
≤n ≤0
L=11;
nn=-10:0;
imp=zeros(L,1);
imp(4)=4.5;
stem(nn,imp)
2.正弦信号
)17/sin(][1n n x π= 0≤n ≤25
L=26;
nn=0:25;
sinus=sin(pi*nn/17);
stem(nn,sinus) 0510152025
-1-0.8
-0.6
-0.4
-0.2
0.2
0.4
0.6
0.8
1
)17/sin(][2n n x π= -15≤n ≤25
L=41;
nn=-15:25;
sinus=sin(pi*nn/17);
stem(nn,sinus)
-15-10-50510152025
-1-0.8
-0.6
-0.4
-0.2
0.2
0.4
0.6
0.8
1
)2/3sin(][3ππ+=n n x -10≤n≤10
L=21;
nn=-10:10;
sinus=sin(3*pi*nn+pi/2)
stem(nn,sinus)
sinus =
Columns 1 through 12
1 -1 1 -1 1 -1 1 -1 1 -1 1
-1
Columns 13 through 21
1 -1 1 -1 1 -1 1 -1 1
-10-8-6-4-20246810
-1-0.8
-0.6
-0.4
-0.2
0.2
0.4
0.6
0.8
1
)23/cos(][4n n x π= 0≤n ≤50
L=51;
nn=0:50;
cosus=cos(pi/sqrt(23)*nn);
stem(nn,cosus) 05101520253035404550
-1-0.8
-0.6
-0.4
-0.2
0.2
0.4
0.6
0.8
1
3.指数信号
(a).x[n]=(0.9)n
function y=genexp(b,n0,L)
if(L<=0)
error('GENEXP:length not positive') end
nn=n0+[1:L]'-1;
y=b.^nn;
end
genexp(0.9,0,21)
ans =
1.0000
0.9000
0.8100
0.7290
0.6561
0.5905
0.5314
0.4783
0.4305
0.3874
0.3487
0.3138
0.2824
0.2542
0.2288
0.2059
0.1853
0.1668
0.1501
0.1351
0.1216
nn=0:20;
stem(nn,ans)
02
46810121416182000.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
(b). x[n]=(0.9)
n genexp(0.9,0,20)
genexp(0.9,0,21)
ans =
1.0000
0.9000
0.8100
0.7290
0.6561
0.5905
0.5314
0.4783
0.4305
0.3874
0.3487
0.3138
0.2824
0.2542
0.2288
0.2059
0.1853
0.1668
0.1501
0.1351
0.1216
sum(ans)
ans =
8.9058
A=(1-0.9^21)/(1-0.9)
A =
8.9058
( c ).
L=20;
N =2:L;
x1=0.9.^n;
stem(n,x1);
246810121416182000.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
>> L=20;
>> n=1:L-1;
>> x2=0.9*0.9.^n;
>> stem(n,x2)
02468101214161820
00.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
( d ). >> L=21;
>> nn=0:L-1;
>> imp=zeros(L,1);
>> imp(1)=1;
>> yout=filter(1,[1,-0.9],imp);
>> stem(nn,yout);
0246810121416182000.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
数字信号处理基础实验指导书
《数字信号处理》实验指导书 光电工程学院二○○九年十月
实验一离散时间信号分析 一、实验目的 1．掌握各种常用的序列，理解其数学表达式和波形表示。 2．掌握在计算机中生成及绘制数字信号波形的方法。 3．掌握序列的相加、相乘、移位、反转等基本运算及计算机实现与作用。 4．掌握线性卷积软件实现的方法。 5．掌握计算机的使用方法和常用系统软件及应用软件的使用。 6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列来表示，其中代表序列的第n个数字，n代表时间的序列，n的取值范围为的整数，n取其它值没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号进行等间隔采样，采样间隔为T，得到一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反转、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将和的变量换成，变成和，再将以纵轴为对称轴反褶成。 （2）移位：将移位，得。当为正数时，右移位；当为负数时，左
移位。 （3）相乘：将和的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得。 三、主要实验仪器及材料 微型计算机、Matlab软件6.5或更高版本。 四、实验内容 1．知识准备 认真复习以上基础理论，理解本实验所用到的实验原理。 2．离散时间信号（序列）的产生 利用MATLAB或C语言编程产生和绘制下列有限长序列： （1）单位脉冲序列 （2）单位阶跃序列 （3）矩形序列 （4）正弦型序列 （5）任意序列 3．序列的运算 利用MATLAB编程完成上述两序列的移位、反转、加法、乘法等运算，并绘制运算后序列的波形。 4．卷积运算 利用MATLAB编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制卷积后序列的波形。 5．上机调试并打印或记录实验结果。 6．完成实验报告。 五、实验报告要求 1. 简述实验原理及目的。 2. 给出上述序列的实验结果。 3. 列出计算卷积的公式，画出程序框图，并列出实验程序清单 （可略）（包括必要的程序说明）。 4. 记录调试运行情况及所遇问题的解决方法。 5. 给出实验结果，并对结果做出分析。 6. 简要回答思考题。 1 如何产生方波信号序列和锯齿波信号序列? 2 实验中所产生的正弦序列的频率是多少?是否是周期序列?
数字信号处理实验(吴镇扬)答案-2
（1） 观察高斯序列的时域和幅频特性，固定信号)(n x a 中参数p=8,改变q 的 值，使q 分别等于2、4、8，观察他们的时域和幅频特性，了解当q 取不同值时，对信号序列的时域和幅频特性的影响；固定q=8，改变p,使p 分别等于8、13、14，观察参数p 变化对信号序列的时域和幅频特性的影响，注意p 等于多少时会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。 ()() ?????≤≤=-其他0150,2n e n x q p n a 解：程序见附录程序一： P=8,q 变化时： t/T x a (n ) k X a (k ) t/T x a (n ) p=8 q=4 k X a (k ) p=8 q=4 t/T x a (n ) p=8 q=8 k X a (k ) p=8 q=8 幅频特性 时域特性
t/T x a (n ) p=8 q=8 k X a (k ) p=8 q=8 t/T x a (n ) 5 10 15 k X a (k ) p=13 q=8 t/T x a (n ) p=14 q=8 5 10 15 k X a (k ) p=14 q=8 时域特性幅频特性 分析： 由高斯序列表达式知n=p 为期对称轴； 当p 取固定值时，时域图都关于n=8对称截取长度为周期的整数倍，没有发生明显的泄漏现象；但存在混叠，当q 由2增加至8过程中，时域图形变化越来越平缓，中间包络越来越大，可能函数周期开始增加，频率降低，渐渐小于fs/2,混叠减弱； 当q 值固定不变，p 变化时，时域对称中轴右移，截取的时域长度渐渐地不再是周期的整数倍，开始无法代表一个周期，泄漏现象也来越明显，因而图形越来越偏离真实值， p=14时的泄漏现象最为明显，混叠可能也随之出现；
数字信号处理实验一
实验一 离散时间信号分析 班级 信息131班 学号 201312030103 姓名 陈娇 日期 一、实验目的 掌握两个序列的相加、相乘、移位、反褶、卷积等基本运算。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为+∞<<∞-n 的整数，n 取其它值)(n x 没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t x a 进行等间隔采样，采样间隔为T ，得到)}({nT x a 一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）） （n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 ∑∞ -∞==-= m n h n x m n h m x n y )(*)()()()( 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。
（2）移位：将)(m h -移位n ，得)(m n h -。当n 为正数时，右移n 位；当n 为负数时，左移n 位。 （3）相乘：将)(m n h -和)(m x 的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得)(n y 。 三、主要实验仪器及材料 微型计算机、Matlab6.5 教学版、TC 编程环境。 四、实验内容 （1）用Matlab 或C 语言编制两个序列的相加、相乘、移位、反褶、卷积等的程序； （2）画出两个序列运算以后的图形； （3）对结果进行分析； （4）完成实验报告。 五、实验结果 六、实验总结
数字信号处理-实验报告
学生实验报告 （理工类） 课程名称：数字信号处理专业班级：通信（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)：
数字信号处理实验报告(实验1_4)
实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B
（2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；
（5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；
（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π]
（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear;
数字信号处理实验
实验一 离散傅里叶变换（DFT ）对确定信号进行谱分析 一．实验目的 1．加深对DFT 算法原理和基本性质的理解。 2．熟悉DFT 算法和原理的编程方法。 3．学习用DFT 对信号进行谱分析的方法，了解可能出现的误差及其原因，以便在实际中正确利用。 二．实验原理 一个连续信号)(t x a 的频谱可以用其傅里叶变换表示，即 dt e t x j X t j a a Ω-∞ ∞ -? = Ω)()( 若对)(t x a 进行理想采样可得采样序列 )(|)()(nT x t x n x a nT t a === 对)(n x 进行DTFT ，可得其频谱为： ∑∞ -∞ =-= n n j j e n x e X ωω )()( 其中数字频率ω与模拟频率Ω的关系为： s f T Ω = Ω=ω )(n x 的DFT 为∑∞ -∞ =-= n nk N j e n x k X π 2)()( 若)(t x a 是限带信号，且在满足采样定理的条件下，)(ω j e X 是)(Ωj X a 的周期延拓， )(k X 是)(ωj e X 在单位圆上的等间隔采样值，即k N j e X k X πωω2| )()(= =。 为在计算机上分析计算方便，常用)(k X 来近似)(ω j e X ，这样对于长度为N 的有限 长序列（无限长序列也可用有限长序列来逼近），便可通过DFT 求其离散频谱。 三．实验内容 1．用DFT 对下列序列进行谱分析。 （1）)()04.0sin(3)(100n R n n x π=
1 （2）]0,0,0,0,0,0,0,0,1,1,1,1[)(=n x 2．为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 )52.0cos()48.0cos()(n n n x ππ+= （1）当0≤n ≤10时，确定并画出x(n)的离散傅里叶变换。 （2）当0≤n ≤100时，确定并画出x(n)的离散傅里叶变换。 四．实验结果 1. （1） （2）
视频信号处理实验报告
中南大学 实验报告（实验一） 实验名称 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（对应解码配置文
数字信号处理实验答案完整版
数字信号处理实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】
实验一熟悉Matlab环境 一、实验目的 1.熟悉MATLAB的主要操作命令。 2.学会简单的矩阵输入和数据读写。 3.掌握简单的绘图命令。 4.用MATLAB编程并学会创建函数。 5.观察离散系统的频率响应。 二、实验内容 认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉了MATLAB基本命令的基础上，完成以下实验。 上机实验内容： （1）数组的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。 clear all; a=[1 2 3 4]; b=[3 4 5 6]; c=a+b; d=a-b; e=a.*b; f=a./b; g=a.^b; n=1:4; subplot(4,2,1);stem(n,a); xlabel('n');xlim([0 5]);ylabel('A'); subplot(4,2,2);stem(n,b); xlabel('n');xlim([0 5]);ylabel('B'); subplot(4,2,3);stem(n,c); xlabel('n');xlim([0 5]);ylabel('C'); subplot(4,2,4);stem(n,d); xlabel('n');xlim([0 5]);ylabel('D'); subplot(4,2,5);stem(n,e); xlabel('n');xlim([0 5]);ylabel('E'); subplot(4,2,6);stem(n,f); xlabel('n');xlim([0 5]);ylabel('F'); subplot(4,2,7);stem(n,g); xlabel('n');xlim([0 5]);ylabel('G'); （2）用MATLAB实现下列序列： a) x(n)= 0≤n≤15 b) x(n)=e+3j)n 0≤n≤15 c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15 d) 将c)中的x(n)扩展为以16为周期的函数x(n)=x(n+16),绘出四个周期。
信号检测论有无法实验报告剖析
------------------------------------------------------------------------------- 实验报告信息栏 系别心理系年级 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年，坦纳和斯韦茨等人首次应用的信号检测论，正好解决了这个问题。 信号检测论的研究对象是信息传播系统中信号的接收问题。在心理学中，它是借助于数学的形式描述“接收者”在某一观察时间内将掺有噪音的信号从噪音中辨别出来。 信号检测论应用于心理学中的基本原理是：将人的感官、中枢分析综合过程看作是一个信息处理系统，应用信号检测论中的一些概念、原理对它进行分析。信号检测论在心理学中具体应用时，常把刺激变量当作信号，把对刺激变量起干扰作用的因素当作噪音，这样就可以把人接收外界刺激时的分辨问题等效于一个在噪音中检测信号的问题，从而便可以应用信号检测论来处理心理学中的实验结果。 信号检测论的理论基础是统计决策。信号检测论本身就是一个以统计判定为根据的理论。它的基本原理是：根据某一观察到的事件，从两个可选择的方面选
数字信号处理实验(吴镇扬)答案-4
实验四 有限长单位脉冲响应滤波器设计 朱方方 0806020433 通信四班 (1) 设计一个线性相位FIR 高通滤波器，通带边界频率为0.6π，阻带边界频率为0.4π，阻 带衰减不小于40dB 。要求给出h(n)的解析式，并用MATLAB 绘出时域波形和幅频特性。 解： (1) 求数字边界频率: 0.6 , .c r ωπωπ== (2) 求理想滤波器的边界频率: 0.5n ωπ= (3) 求理想单位脉冲响应: []d s i n ()s i n [()] () ()1n n n n n n h n n παωαα παωα π?-- -≠??-=? ? -=?? (4) 选择窗函数。阻带最小衰减为-40dB ，因此选择海明窗(其阻带最小衰减为-44dB)；滤 波器的过渡带宽为0.6π-0.4π=0.2π，因此 6.21 0.231 , 152 N N N ππα-=?=== (5) 求FIR 滤波器的单位脉冲响应h(n): []31d sin (15)sin[0.5(15)] 1cos ()15()()()15(15)1 15 n n n R n n h n w n h n n n ππππ?---????-? ?≠? ???==-???? ? ?=? 程序： clear; N=31; n=0:N-1; hd=(sin(pi*(n-15))-sin(0.5*pi*(n-15)))./(pi *(n-15)); hd(16)=0.5; win=hanning(N); h=win'.*hd; figure; stem(n,h); xlabel('n'); ylabel('h(n)'); grid; title('FIR 高通滤波单位脉冲响应h(n)'); [H,w]=freqz(h,1); H=20*log10(abs(H)); figure;3 plot(w/pi,H); axis([0 1 -100 10]); xlabel('\omega/\pi'); ylabel('幅度/dB'); grid; title('FIR 高通滤波器，hanning 窗，N=31');
数字信号处理实验三
实验三：离散LSI 系统的频域分析 一、实验内容 2、求以下各序列的z 变换： 12030() ()sin() ()sin()n an x n na x n n x n e n ωω-=== 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z 变换 0 312342 1 1() () () ()() 1j z z z z X z X z X z X z z a z a z e z ω---= = = = ---- 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a 课程名称 数字信号 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期
x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） (0.3)()(1)(1) z z H z z j z j -= +-++ z1=[0,0.3]';p1=[-1+j,-1-j]';k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title('极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 -1 -0.5 00.51 -2 -1.5-1-0.500.511.5 2Real Part I m a g i n a r y P a r t 极点在单位圆外 n (samples) A m p l i t u d e Impulse Response
信号处理实验报告、
第一题 如何用计算机模拟一个随机事件，并估计随机事件发生的概率以计算圆周率π。 解： （一）蒙特卡洛方法可用于近似计算圆周率：让计算机每次随机生成两个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 的数字信号处理基本操作 一、 实验目的：学会运用MA TLAB 表示的常用离散时间信号；学会运用MA TLAB 实现离 散时间信号的基本运算。 二、 实验仪器：电脑一台，MATLAB6.5或更高级版本软件一套。 三、 实验内容： （一） 离散时间信号在MATLAB 中的表示 离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用)(n x 来表示，自变量必须是整数。 离散时间信号的波形绘制在MATLAB 中一般用stem 函数。stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill ”、“filled ”，或者参数“.”。由于MATLAB 中矩阵元素的个数有限，所以MA TLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。 1. 单位取样序列 单位取样序列)(n δ，也称为单位冲激序列，定义为 ) 0() 0(0 1)(≠=?? ?=n n n δ 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n =0处是取确定的值1。在MATLAB 中，冲激序列可以通过编写以下的impDT .m 文件来实现，即 function y=impDT(n) y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n 必须为整数或整数向量。 【实例1-1】 利用MATLAB 的impDT 函数绘出单位冲激序列的波形图。 解：MATLAB 源程序为 >>n=-3:3; >>x=impDT(n); >>stem(n,x,'fill'),xlabel('n'),grid on >>title('单位冲激序列') >>axis([-3 3 -0.1 1.1]) 程序运行结果如图1-1所示。 图1-1 单位冲激序列 数字信处理上机实验答 案全 Document number【SA80SAB-SAA9SYT-SAATC-SA6UT-SA18】 第十章 上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一系统响应及系统稳定性。 实验二时域采样与频域采样。 实验三用FFT对信号作频谱分析。 实验四 IIR数字滤波器设计及软件实现。 实验五 FIR数字滤波器设计与软件实现 实验六应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。也可以用MATLAB语言的工具箱函数conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。 系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。 系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。 实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。系统的稳态输出是指当∞ n时，系统的输出。如果系统稳定，信号加入 → 系统后，系统输出的开始一段称为暂态效应，随n的加大，幅度趋于稳定，达到稳态输出。 注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤 数字信号处理实验四 第一题结果： （1）没有增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 %H(3,13) = 0.75;H(5,11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 （2）增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 H(3) = 0.75;H(13) = 0.75;H(5) = 0.25;H(11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果： 实验一：信号的表示 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; 本科生实验报告 实验课程数字信号处理基础 学院名称地球物理学院 专业名称地球物理学 学生姓名 学生学号 指导教师王山山 实验地点5417 实验成绩 二〇一四年十一月二〇一四年十二月 填写说明 1、适用于本科生所有的实验报告（印制实验报告册除外）； 2、专业填写为专业全称，有专业方向的用小括号标明； 3、格式要求： ①用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 ②打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm， 左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用小五号字底端居中。 ③具体要求： 题目（二号黑体居中）； 摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）； 关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)； 正文部分采用三级标题； 第1章××(小二号黑体居中，段前0.5行) 1.1 ×××××小三号黑体×××××（段前、段后0.5行） 1.1.1小四号黑体（段前、段后0.5行） 参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照《参考文献著录规则（GB/T 7714－2005）》。 实验一生成离散信号并计算其振幅谱 并将信号进行奇偶分解 一、实验原理 单位脉冲响应h(t)=exp(-a*t*t)*sin(2*3.14*f*t)进行离散抽样，分别得到t=0.002s，0.009s，0.011s采样的结果。用Excel软件绘图显示计算结果。并将信号进行奇偶分解，分别得到奇对称信号h(n)-h(-n)与偶对称信号h(n)+h(-n)。用Excel 软件绘图显示计算结果。 二、实验程序代码 （1）离散抽样 double a,t; a=2*f*f*log(m); int i; for(i=0;i 南京信息工程大学数字信号处理实验报告 学院：电子与信息工程学院 班级：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的数字信号处理基本
数字信处理上机实验答案全
数字信号处理实验4
数字信号处理实验报告
数字信号处理基础实验报告_
数字信号处理实验报告