抽取之后，原来的频谱被展宽了D 倍，为了避免混叠，需要满足f
f
m
s
D ≤
。
()
3、运行结果运行的结果
（1）D=1时，是原来的音乐，总共是20秒，运行结果 d=1_fs.wav.
（2）D=2时，以fs保存，运行结果d=2_fs,wav，以fs/2保存，运行结果
对时域截短，使频谱变宽拖尾，称为泄露程序分析
clc;close all;
n1=0:130;n2=0:220;
x1=sin(pi*n1/32);
x2=sin(pi*n2/32);
X1=fft(x1);
X2=fft(x2);
窗函数设计法：
窗函数法的设计步骤：
2.MATLAB函数工具
打开新的工程，输入“fdatool”。之后运行，打开滤波器设计应用程序。
通过改变滤波器类型（lowpass、highpass）选择FIR window，并设置Fpass Fstoop来得到一个低通滤波器。
设置完毕之后打开generate matlab rade来完成程序编写
电子工程系 信号与系统课程实验报告 2011-----2012学年第一学期 专业: 电子信息工程技术班级: 学号 : 姓名: 指导教师: 实常用连续时间信号的实现
一、实验目的 （1）了解连续时间信号的特点； （2）掌握连续时间信号表示的向量法和符号法； （3）熟悉MATLAB Plot函数等的应用。 二、实验原理 1、信号的定义 信号是随时间变化的物理量。信号的本质是时间的函数。 2、信号的描述 1）时域法 时域法是将信号表示成时间的函数f(t)来对信号进行描述的方法。信号的时间特性指的是信号的波形出现的先后，持续时间的长短，随时间变化的快慢和大小，周期的长短等。 2）频域（变换域）法 频域法是通过正交变换，将信号表示成其他变量的函数来对信号进行描述的方法。一般常用的是傅立叶变换。信号的频域特性包括频带的宽窄、频谱的分布等。 信号的频域特性与时域特性之间有着密切的关系。 3、信号的分类 按照特性的不同，信号有着不同的分类方法。 （1）确定性信号：可以用一个确定的时间函数来表示的信号。 随机信号：不可以用一个确定的时间函数来表示，只能用统计特性加以描述的信号。 （2）连续信号：除若干不连续的时间点外，每个时间点在t上都有对应的数值信号。离散信号：只在某些不连续的点上有数值，其他时间点上信号没有定义的信号。 （3）周期信号：存在T，使得等式f(t+T)=f(t)对于任意时间t都成立的信号。非周期信号：不存在使得等式f(t+T)=f(t)对于任意时间t都成立的信号。 绝对的周期信号是不存在的，一般只要在很长时间内慢走周期性就可以了。 （4）能量信号：总能量有限的信号。 功率信号:平均功率有限切非零的信号。 （5）奇信号：满足等式f(t)=--f(--t)的信号。偶信号：满足等式f(t)=f(--t)的信号。 三、涉及的MATLAB函数 1、plot函数 功能：在X轴和Y轴方向都按线性比例绘制二维图形。 调用格式： Plot（x,y）:绘出相x对y的函数线性图。 Plot（x1,y1,x2,y2,…..）:会出多组x对y的线性曲线图。 2、ezplot函数 功能：绘制符号函数在一定范围内的二维图形。简易绘制函数曲线。 调用格式： Ezplot (fun):在[-2π，2π]区间内绘制函数。 Ezplot (fun，［min，max］)：在［min,max］区间内绘函数。 Ezplot (funx,funy):定义同一曲面的函数，默认的区间是[0, 2π]。】 3、sym函数 功能：定义信号为符号的变量。 调用格式：sym(fun):fun为所要定义的表达式。 4、subplot函数
一、实验名称：基本信号的产生 二、实验目的：I 利用MATLAB 产生连续信号并作图 II 利用MATLAB 产生离散序列并作图 III 利用MATLAB 进行噪声处理 三、 实验内容： I 利用MATLAB 产生下列连续信号并作图 ①X(t)=-2u(t-1),-1 -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 实验一：信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示： 图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示： 图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图： 图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示： 图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示： 图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t) 学生实验报告 （理工类） 课程名称：数字信号处理专业班级：通信（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 电 子 科 技 大 学 信号与系统实验报告 学生姓名： 杜杰 学 号： 2014030103007 指导教师：张鹰 一、实验室名称：信号与系统实验室 二、实验项目名称：连续系统的幅频特性测量 三、实验原理： 正弦波信号)cos()(0t A t x ω=输入连续LTI 系统，输出)(t y 仍为正弦波信号。 图信号输入连续LTI 系统 图中， )(cos()()(000ωωωj H t j H A t y ∠+=) 通过测量输入)(t x 、输出)(t y 的正弦波信号幅度，计算输入、输出的正弦波信号幅度比值，可以得到系统的幅频特性在0ω处的测量值)(0ωj H 。改变0ω可以测出不同频率处的系统幅频特性。 四、实验目的： 使学生对系统的频率特性有深入了解。 五、实验内容： 实验内容（一）、低通滤波器的幅频特性测量 实验内容（二）、带通滤波器的幅频特性测量 六、实验器材（设备、元器件）： 数字信号处理实验箱、信号与系统实验板的低通滤波器模块U11、高通滤波器模块U21、PC 机端信号与系统实验软件、＋5V 电源 七、实验步骤：打开PC 机端软件SSP.EXE ，在下拉菜单“实验选择”中选择 “实验三”；使用串口电缆连接计算机串口和实验箱串口，打开实验箱电源。 实验内容（一）、低通滤波器的幅频特性测量 (x ) (t y 实验步骤： 1、信号选择：按实验箱键盘“3”选择“正弦波”，再按“＋”或“－”依次 选择一个频率。 2、连接接口区的“输入信号1”和“输出信号”，如下图所示。点击SSP软 件界面上的按钮，观察输入正弦波。将正弦波频率值和幅度值(Vpp/2, Vpp为峰-峰值)记录于表1。 观察输入正弦波的连线示意图 3、按下图的模块连线示意图连接各模块。 实验三实验内容（一）模块连线示意图 4、点击SSP软件界面上的按钮，观察输入正弦波通过连续系统的 响应波形；适当调整X、Y轴的分辨率可得到如下图所示的实验结果。将输出正弦波的幅度值(Vpp/2, Vpp为峰-峰值)记录于表1中。 输入正弦波和响应波形 中南大学 实验报告（实验一） 实验名称 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年，坦纳和斯韦茨等人首次应用的信号检测论，正好解决了这个问题。 信号检测论的研究对象是信息传播系统中信号的接收问题。在心理学中，它是借助于数学的形式描述“接收者”在某一观察时间内将掺有噪音的信号从噪音中辨别出来。 信号检测论应用于心理学中的基本原理是：将人的感官、中枢分析综合过程看作是一个信息处理系统，应用信号检测论中的一些概念、原理对它进行分析。信号检测论在心理学中具体应用时，常把刺激变量当作信号，把对刺激变量起干扰作用的因素当作噪音，这样就可以把人接收外界刺激时的分辨问题等效于一个在噪音中检测信号的问题，从而便可以应用信号检测论来处理心理学中的实验结果。 信号检测论的理论基础是统计决策。信号检测论本身就是一个以统计判定为根据的理论。它的基本原理是：根据某一观察到的事件，从两个可选择的方面选 ( 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．实验目的 ● 学会使用MATLAB 产生各种常见的连续时间信号与离散时间信号； ● 通过MATLAB 中的绘图工具对产生的信号进行观察，加深对常用信号的 理解； ● 熟悉MATLAB 的基本操作，以及一些基本函数的使用，为以后的实验奠 定基础。 2．实验内容 ⑴ 运行以上九个例子程序，掌握一些常用基本信号的特点及其MATLAB 实现方法；改变有关参数，进一步观察信号波形的变化。 ⑵ 在 k [10:10]=- 范围内产生并画出以下信号： a) 1f [k][k]δ=； b) 2f [k][k+2]δ=； c) 3f [k][k-4]δ=； d) 4f [k]2[k+2][k-4]δδ=-。 源程序： k=-10:10; f1k=[zeros(1,10),1,zeros(1,10)]; subplot(2,2,1) stem(k,f1k) title('f1[k]') f2k=[zeros(1,8),1,zeros(1,12)]; subplot(2,2,2) stem(k,f2k) title('f2[k]') f3k=[zeros(1,14),1,zeros(1,6)]; subplot(2,2,3) stem(k,f3k) title('f3[k]') f4k=2*f2k-f3k; subplot(2,2,4) stem(k,f4k) title('f4[k]') ⑶ 在 k [0:31]=范围内产生并画出以下信号： a) ()()k k 144f [k]sin cos π π=； b) ()2k 24f [k]cos π =； c) ()()k k 348f [k]sin cos π π=。 请问这三个信号的基波周期分别是多少？ 源程序： k=0:31; f1k=sin(pi/4*k).*cos(pi/4*k); subplot(3,1,1) stem(k,f1k) title('f1[k]') f2k=(cos(pi/4*k)).^2; subplot(3,1,2) stem(k,f2k) title('f2[k]') f3k=sin(pi/4*k).*cos(pi/8*k); subplot(3,1,3) stem(k,f3k) title('f3[k]') 其中f1[k]的基波周期是4, f2[k]的基波周期是4, f3[k]的基波周期是16。 数字信号处理 第四次实验报告 一、 实验目的 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.实现单位采样序列、单位阶跃序列、矩形序列程序及绘图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; 实验三 常见信号的MATLAB 表示及运算 一、实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二、实验原理 根据MATLAB 的数值计算功能和符号运算功能，在MA TLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MA TLAB 语句表示出信号后，就可以利用MA TLAB 中的绘图命令绘制出直观的信号波形了。 1.连续时间信号 从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t =的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 ⑵ 符号运算表示法 如果一个信号或函数可以用符号表达式来表示，那么我们就可以用前面介绍的符号函数专用绘图命令ezplot()等函数来绘出信号的波形。 ⑶ 常见信号的MATLAB 表示 单位阶跃信号 单位阶跃信号的定义为：10 ()0 t u t t >?=? 0); %定义函数体，即函数所执行指令 《数字信号处理》 实验报告 课程名称：《数字信号处理》 学院：信息科学与工程学院 专业班级：通信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 ωω=ΩΩ= Harbin Institute of Technology 匹配滤波器实验报告 课程名称：信号检测理论 院系：电子与信息工程学院 姓名：高亚豪 学号：14SD05003 授课教师：郑薇 哈尔滨工业大学 1. 实验目的 通过Matlab 编程实现对白噪声条件下的匹配滤波器的仿真，从而加深对匹配滤波器及其实现过程的理解。通过观察输入输出信号波形及频谱图，对匹配处理有一个更加直观的理解，同时验证匹配滤波器具有时间上的适应性。 2. 实验原理 对于一个观测信号()r t ，已知它或是干扰与噪声之和，或是单纯的干扰， 即 0()()()()a u t n t r t n t +?=?? 这里()r t ，()u t ，()n t 都是复包络，其中0a 是信号的复幅度，()u t 是确知的归一化信号的复包络，它们满足如下条件。 2|()|d 1u t t +∞ -∞=? 201||2 a E = 其中E 为信号的能量。()n t 是干扰的均值为0，方差为0N 的白噪声干扰。 使该信号通过一个线性滤波系统，有效地滤除干扰，使输出信号的信噪比在某一时刻0t 达到最大，以便判断信号的有无。该线性系统即为匹配滤波器。 以()h t 代表系统的脉冲响应，则在信号存在的条件下，滤波器的输出为 0000()()()d ()()d ()()d y t r t h a u t h n t h τττττττττ+∞+∞+∞ =-=-+-??? 右边的第一项和第二项分别为滤波器输出的信号成分和噪声成分，即 00()()()d x t a u t h τττ+∞ =-? 0 ()()()d t n t h ?τττ+∞ =-? 则输出噪声成分的平均功率（统计平均）为 2 20E[|()|]=E[|()()d |]t n t h ?τττ+∞ -? **00*000200 =E[()(')]()(')d d '=2()(')(')d d ' 2|()|d n t n t h h N h h N h ττττττδττττττττ+∞+∞+∞+∞+∞ ---=?? ?? ? 而信号成分在0t 时刻的峰值功率为 22 20000|()||||()()d |x t a u t h τττ+∞ =-? 输出信号在0t 时刻的总功率为 22000E[|()|]E[|()()|]y t x t t ?=+ 22**0000002200E[|()||()|()()()()] |()|E[|()|] x t t x t t t x t x t t ????=+++=+ 上式中输出噪声成分的期望值为0，即0E[()]0t ?=，因此输出信号的功率 成分中只包含信号功率和噪声功率。 则该滤波器的输出信噪比为 222000022000|||()()d ||()|E[|()|]2|()|d a u t h x t t N h τττρ?ττ+∞ +∞-==?? 根据Schwartz 不等式有 南京信息工程大学数字信号处理实验报告 学院：电子与信息工程学院 班级：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信号与系统实验报告_1(常用信号的分类与观察)
数字信号处理-实验报告
信号处理实验报告、
信号实验报告
视频信号处理实验报告
哈尔滨工程大学 语音信号处理实验报告
数字信号处理实验报告
信号检测论有无法实验报告剖析
工程信号处理实验报告
信号与系统实验报告
信号处理实验报告
数字信号处理实验报告
信号与系统实验报告汇总
数字信号处理实验报告(同名22433)
信号检测实验报告
数字信号处理实验报告
相关主题
文本预览