信号分析与处理综合作业

作业一：信号分析与处理在某一方面的应用综述

要求：

1、查阅学校图书馆的纸质和电子图书与期刊，撰写一篇信号分析与处理在某一方面的应用的小论文，要求至少查阅5篇以上的论文；

2、论文的格式要按照期刊论文发表的格式，由以下几个部分组成；

题目，中文摘要，英文摘要，中文关键词，英文关键词，前言，正文，结论，参考文献；

3、论文中引用的部分要注明；

4、论文字数不少于1500字，要求手写，不能打印。

作业二：声音信号的分析与处理

1、利用windows里的录音机等工具，记录声音“专业班级姓名”，如“勘查083班张三”，

最好存储为wav的单声道格式；

2、在公用邮箱中下载地震数据文件和测井数据文件；

3、对上述信号的频谱进行分析，写出其频谱特征；

4、对上述信号进行时域简单处理，如反褶，尺度变换，试听处理后的结果；；

5、分别设计低通，高通和带通滤波器对此信号进行处理，试听处理后的结果；

6、（选做）对此信号进行纯相位滤波或其他的处理；

7、写出相应的分析与处理的报告，报告中要插入有关的图件（要有分析）和分析与处理的

程序；

8、将报告、程序和处理前后的声音信号压缩成一个文件，发送到邮箱upcsignal@https://www.360docs.net/doc/4718386945.html,，

邮件的标题格式为“专业班级学号姓名”。

五款信号完整性仿真工具介绍

现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 (一)Ansoft公司的仿真工具现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 Ansoft的信号完整性工具采用一个仿真可解决全部设计问题： SIwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何数量的过孔和信号引线条构成。仿真结果采用先进的3D图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿真器。 (二)SPECCTRAQuest Cadence的工具采用Sun的电源层分析模块： Cadence Design Systems的SpecctraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI。该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具：（1）SigXplorer可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在PCB详细设计前使用此工具，对互连线的不同情况进行仿真，把仿真结果存为拓扑结构模板，在后期详细设计中应用这些模板进行设计。（2）DF/Signoise工具是信号仿真分析工具，可提供复杂的信号延时和信号畸变分析、IBIS 模型库的设置开发功能。SigNoise是SPECCTRAQUEST SI Expert和SQ Signal Explorer Expert进行分析仿真的仿真引擎，利用SigNoise可以进行反射、串扰、SSN、EMI、源同步及系统级的仿真。（3）DF/EMC工具——EMC分析控制工具。（4）DF/Thermax——热分析控制工具。 SPECCTRAQuest中的理想高速PCB设计流程：由上所示，通过模型的验证、预布局布线的space分析、通过floorplan制定拓朴规则、由规

信号分析与处理答案第二版完整版

信号分析与处理答案第二版 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】

第二章习题参考解答求下列系统的阶跃响应和冲激响应。 (1) 解当激励为时，响应为，即：由于方程简单，可利用迭代法求解：，， …，由此可归纳出的表达式：利用阶跃响应和冲激响应的关系，可以求得阶跃响应： (2) 解 (a)求冲激响应，当时，。特征方程，解得特征根为。所以： …(2.1.2.1) 通过原方程迭代知，，，代入式(2.1.2.1)中得：解得，代入式(2.1.2.1)： …(2.1.2.2) 可验证满足式(2.1.2.2)，所以： (b)求阶跃响应通解为特解形式为，，代入原方程有，即完全解为通过原方程迭代之，，由此可得解得，。所以阶跃响应为： (3)

解 (4) 解当t>0时，原方程变为：。 …(2.1.3.1) …(2.1.3.2) 将(2.1.3.1)、式代入原方程，比较两边的系数得：阶跃响应：求下列离散序列的卷积和。 (1) 解用表格法求解 (2) 解用表格法求解 (3) 和如题图2.2.3所示解用表格法求解

(4) 解 (5) 解 (6) 解参见右图。当时：当时：当时：当时：当时： (7) ，解参见右图：当时：当时：当时：当时：当时： (8) ，解参见右图

当时：当时：当时：当时： (9) ，解 (10) ，解或写作：

求下列连续信号的卷积。 (1) ，解参见右图：当时：当时：当时：当时：当时：当时： (2) 和如图2.3.2所示解当时：当时：当时：当时：当时： (3) ，解 (4) ，解 (5) ，解参见右图。当时：当时：当时：

《测试技术与信号处理》习题答案-华科版

《测试技术与信号处理》习题答案第二章信号分析基础 1、请判断下列信号是功率信号还是能量信号：（1）)()(10cos 2 ∞<<-∞=t e t x t π （2）)()(||10∞<<-∞=-t e t x t 【解】（1）该信号为周期信号，其能量无穷大，但一个周期内的平均功率有限，属功率信号。（2）信号能量：? ∞ ∞ -= =10 1 )(2dt t x E ，属于能量信号。 2、请判断下列序列是否具有周期性，若是周期性的，请求其周期。)8 ()(π-=n j e n x 【解】设周期为N ，则有：8 )8 8()()(N j N n j e n x e N n x ?==+-+π 若满足)()(n x N n x =+，则有1)8/sin()8/cos(8/=-=-N j N e jN 即：k N π28/=，k N π16=，k = 0，1，2，3，… N 不是有理数，故序列不是周期性的。 3、已知矩形单脉冲信号x 0(t)的频谱为X 0(ω)=A τsinc(ωτ/2) ，试求图示三脉冲信号的频谱。【解】三脉冲信号的时域表达式为：)()()()(000T t x t x T t x t x -+++= 根据Fourier 变换的时移特性和叠加特性，可得其频谱： )]cos(21)[2 ( sin )()()()(000T c A e X X e X X T j T j ωωτ τωωωωωω+=++=- 4、请求周期性三角波（周期为T ，幅值为0—A ）的概率分布函数F(x)与概率密度函数p(x) 。【解】在一个周期T 内，变量x (t )小于某一特定值x 的时间间隔平均值为：T A x t i = ? 取n 个周期计算平均值，当∞→n 时，可有概率分布函数：A x nT t n x F i n =?=∞→lim )( 概率密度函数：A dx x dF x p 1 )()(== t -τ/2 0 τ/2 -T T

现代信号处理大作业

现代信号处理大作业姓名：潘晓丹学号：0140349045 班级：A1403492

作业1 LD 算法实现AR 过程估计 1.1 AR 模型 p 阶AR 模型的差分方程为： )()()(1 n w i n x a n x p i i =-+ ∑=，其中)(n w 是均值为0的白噪声。 AR 过程的线性预测方法为：先求得观测数据的自相关函数，然后利用Yule-Walker 方程递推求得模型参数，再根据公式求得功率谱的估计。 Yule-Walker 方程可写成矩阵形式： ??????? ? ????????= ??????? ? ?? ???? ????????????? ??? ??--+-+--000)()2()1(1) 0() 2()1()()2()0()1()2()1()1()0() 1()() 2()1()0(2 σp a a a r p r p r p r p r r r r p r r r r p r r r r p p p xx xx xx xx xx xx xx xx xx xx xx xx xx xx xx xx 1.2 LD 算法介绍 Levinson-Durbin 算法可求解上述问题，其一般步骤为： 1) 计算观测值各自相关系数p j j r xx ,,1, 0),( =；)0(0xx r =ρ；i=1； 2) 利用以下递推公式运算： ) 1(1,...,2,1),()()()() ()()(2 1111 1 1 1 i i i i i i i i i i i j xx i xx i k i j j i a k j a j a k i a j i r j a i r k -=-=--==-?+ -=-----=-∑ρρρ 3) i=i+1，若i>p ，则算法结束；否则，返回(2)。 1.3 matlab 编程实现以AR 模型：xn=12xn-1-12xn-2+w(n)为例，Matlab 程序代码如下： clear; clc;

Altium Designer中进行信号完整性分析

在高速数字系统中，由于脉冲上升/下降时间通常在10到几百p秒，当受到诸如内连、传输时延和电源噪声等因素的影响，从而造成脉冲信号失真的现象；在自然界中，存在着各种各样频率的微波和电磁干扰源，可能由于很小的差异导致高速系统设计的失败；在电子产品向高密和高速电路设计方向发展的今天，解决一系列信号完整性的问题，成为当前每一个电子设计者所必须面对的问题。业界通常会采用在PCB制板前期，通过信号完整性分析工具尽可能将设计风险降到最低，从而也大大促进了EDA设计工具的发展…… 信号完整性（Signal Integrity，简称SI）问题是指高速数字电路中，脉冲形状畸变而引发的信号失真问题，通常由传输线不阻抗匹配产生的问题。而影响阻抗匹配的因素包括信号源的架构、输出阻抗(output impedance)、走线的特性阻抗、负载端的特性、走线的拓朴(topology)架构等。解决的方式可以采用端接(termination)与调整走线拓朴的策略。信号完整性问题通常不是由某个单一因素导致的，而是板级设计中多种因素共同作用的结果。信号完整性问题主要表现形式包括信号反射、信号振铃、地弹、串扰等； 1，Altium Designer信号完整性分析（机理、模型、功能）在Altium Designer设计环境下，您既可以在原理图又可以在PCB编辑器内实现信号完整性分析，并且能以波形的方式在图形界面下给出反射和串扰的分析结果。 Altium Designer的信号完整性分析采用IC器件的IBIS模型，通过对版图内信号线路的阻抗计算，得到信号响应和失真等仿真数据来检查设计信号的可靠性。Altium Designer的信号完整性分析工具可以支持包括差分对信号在内的高速电路信号完整性分析功能。 Altium Designer仿真参数通过一个简单直观的对话框进行配置，通过使用集成的波形观察仪，实现图形显示仿真结果，而且波形观察仪可以同时显示多个仿真数据图像。并且可以直接在标绘的波形上进行测量，输出结果数据还可供进一步分析之用。 Altium Designer提供的集成器件库包含了大量的的器件IBIS模型，用户可以对器件添加器件的IBIS模型，也可以从外部导入与器件相关联的IBIS模型，选择从器件厂商那里得到的IBIS 模型。 Altium Designer的SI功能包含了布线前（即原理图设计阶段）及布线后（PCB版图设计阶段）两部分SI分析功能；采用成熟的传输线计算方法，以及I/O缓冲宏模型进行仿真。基于快速反射和串扰模型，信号完整性分析器使用完全可靠的算法，从而能够产生出准确的仿真结果。布线前的阻抗特征计算和信号反射的信号完整性分析，用户可以在原理图环境下运行SI仿真功能，对电路潜在的信号完整性问题进行分析，如阻抗不匹配等因素。更全面的信号完整性分析是在布线后PCB版图上完成的，它不仅能对传输线阻抗、信号反射和信号间串扰等多种设计中存在的信号完整性问题以图形的方式进行分析，而且还能利用规则检查发现信号完整性问题，同时，Altium Designer还提供一些有效的终端选项，来帮助您选择最好的解决方案。 2，分析设置需求在PCB编辑环境下进行信号完整性分析。为了得到精确的结果，在运行信号完整性分析之前需要完成以下步骤：

信号分析与处理习题

2.1 有一个理想采样系统，其采样角频率Ωs =6π，采样后经理想低通滤波器H a (j Ω)还原，其中 ?? ???≥Ω<Ω=Ωππ 3032 1 )(，，j H a 现有两个输入，x 1(t )=cos2πt ，x 2(t )=cos5πt 。试问输出信号y 1(t )，y 2(t )有无失真？为什么？分析：要想时域采样后能不失真地还原出原信号，则采样角频率Ωs 必须大于等于信号谱最高角频率Ωh 的2倍，即满足Ωs ≥2Ωh 。解：已知采样角频率Ωs =6π，则由香农采样定理，可得因为x 1(t )=cos2πt ，而频谱中最高角频率ππ π32621=< =Ωh ，所以y 1(t )无失真；因为x 2(t )=cos5πt ，而频谱中最高角频率ππ π32 652=>=Ωh ，所以y 2(t )失真。 3.2 设x (n )的傅里叶变换为X (e j ω)，试利用X (e j ω )表示下列序列的傅里叶变换：（1） )1()1()(1n x n x n x --+-= （2） )]()([2 1 )(2n x n x n x -+= * 分析：利用序列翻褶后的时移性质和线性性质来求解，即 )()(ωj e X n x ?，)()(ωj e X n x -?- )()(ωωj m j e X e n m x --?- 解：（1）由于)()]([ω j e X n x DTFT =，)()]([ωj e X n x DTFT -=-，则 )()]1([ωωj j e X e n x DTFT --=- )()]1([ωωj j e X e n x DTFT -=-- 故ωωωωω cos )(2])[()]([1j j j j e X e e e X n x DTFT ---=+= （2）由于)()]([ω j e X n x DTFT * * =- 故)](Re[2 ) ()()]([2ωωωj j j e X e X e X n x DTFT =+= * 3.7 试求下列有限长序列的N 点离散傅里叶变换（闭合形式表达式）：

华中科技大学工程测试与信息处理作业及答案

1 请给出3种家用电器中的传感器及其功能。洗衣机：水位传感器冰箱：温度传感器彩电：亮度传感器热水器：温度传感器空调：温度传感器 2 请给出智能手机中用到的测试传感器。重力传感器、三维陀螺仪、GPS 、温度传感器、亮度传感器、摄像头等。 3 系统地提出（或介绍）你了解的（或设想的）与工程测试相关的某一（小）问题。考虑通过本门课程的学习，你如何来解决这一问题。（注意：本题的给出的答案将于课程最后的综合作业相关联，即通过本课程的学习，给出详细具体可行的解决方案）第二章普通作业1 请写出信号的类型 1) 简单周期信号 2) 复杂周期信号 3) 瞬态信号准周期信号 4) 平稳噪声信号 5) 非平稳噪声信号第二章信号分析基础测试题 1. 设时域信号x(t)的频谱为X(f)，则时域信号（C ）的频谱为X(f ＋fo ）。 A . B. C. D. 2. 周期信号截断后的频谱必是（A ）的。 A. 连续 B. 离散 C. 连续非周期 D. 离散周期 3. 不能用确定的数学公式表达的信号是 (D) 信号。 A 复杂周期 B 非周期 C 瞬态 D 随机 4. 信号的时域描述与频域描述通过 (C) 来建立关联。 A 拉氏变换 B 卷积 C 傅立叶变换 D 相乘 5. 以下 (B) 的频谱为连续频谱。 A 周期矩形脉冲 B 矩形窗函数 C 正弦函数 D 周期方波 6. 单位脉冲函数的采样特性表达式为（A ）。 A )(d )()(00t x t t t t x =-?∞ ∞-δ B )()(*)(00t t x t t t x -=-δ C )()(*)(t x t t x =δ D 1)(?t δ 思考题 1) 从下面的信号波形图中读出其主要参数。

信号与系统大作业

中北大学信号与系统综合性报告学院：仪器与电子学院专业：电子科学与技术学号姓名：王鹏学号姓名：张艺超学号姓名：郭靖锋学号姓名：蔡宪庆学号姓名：指导教师: 张晓明 2019年5 月13 日

1 设计题目时频域语音信号的分析与处理 2 设计目标对语音信号进行时频域分析和处理的基本方法 3 设计要求 1）分别录制一段男生和女生语音文件及相应有明显高频或低频干扰的语音文件*.wav，并将文件导入Matlab中； 2）分别分析各段语音的频谱，绘制其频谱图，分析语音信号和干扰信号的频段； 3）设计相应的滤波器，剔除含干扰的语音段的干扰信号，并分析滤波信号的频谱； 4）生成滤波后的语音文件，分析听觉效果。 4 理论分析声音作为一种波，频率在20 Hz~20 kHz之间的声音是可以被人耳识别的通过查阅资料显示，实际人声频率范围男：低音82～392Hz，基准音区64～523Hz 男中音123～493Hz，男高音164～698Hz 女：低音82～392Hz，基准音区160～1200Hz 女低音123～493Hz，女高音220～1.1KHz 声音作为波的一种，频率和振幅就成了描述波的重要属性，频率的大小与我们通常所说的音高对应，而振幅影响声音的大小。声音可以被分解为不同频率不同强度正弦波的叠加。这种变换（或分解）的过程，称为傅立叶变换(Fourier Transform)。傅里叶变换之后可以得到男女声的频谱，从而分析男女声的特点，观察男女声频率集中的区域，在声音中加入高频噪声，分析高频噪声频率的分布，从而设计巴特沃斯滤波器进行滤波。 5 实验内容及步骤 5.1 获取音频文件 5.1.1 通过手机录音可直接获取wav音频文件，对于噪声的添加，我们选择单独录制高频件，读取音频数据，在时域领域上相加，便获取到含有高频噪声的音频 5.2 音频的时域处理 5.2.1 wav属于无损音乐格式的一种，其文件包含采样频率，左右声道数据，在处理时，由于我们使用的是matlab2012a，且录制时只有一个声道，可使用函数wavread()读取到一个一维数组，使用plot函数即可获取其音频时域图像 5.3 音频的频域处理 5.3.1 对于音频数组，我们采用fft函数进行傅里叶变换，获取到的是对称的复数数组，数组的前一半即为其频域，同样使用plot将其画出。 5.3.2 观察频域图，分析男女声特点。 5.4 噪声的去除 5.4.1 分析高频噪声频谱，找到合适的截止频率，设计巴特沃斯滤波器对高频噪声进行过滤。 5.4．2 将去除噪声的数组转换成音频文件

于博士信号完整性分析入门(修改)

于博士信号完整性分析入门于争博士 https://www.360docs.net/doc/4718386945.html, for more information,please refer to https://www.360docs.net/doc/4718386945.html, 电设计网欢迎您

什么是信号完整性? 如果你发现，以前低速时代积累的设计经验现在似乎都不灵了，同样的设计，以前没问题，可是现在却无法工作，那么恭喜你，你碰到了硬件设计中最核心的问题：信号完整性。早一天遇到，对你来说是好事。在过去的低速时代，电平跳变时信号上升时间较长，通常几个ns。器件间的互连线不至于影响电路的功能，没必要关心信号完整性问题。但在今天的高速时代，随着IC输出开关速度的提高，很多都在皮秒级，不管信号周期如何，几乎所有设计都遇到了信号完整性问题。另外，对低功耗追求使得内核电压越来越低，1.2v内核电压已经很常见了。因此系统能容忍的噪声余量越来越小，这也使得信号完整性问题更加突出。广义上讲，信号完整性是指在电路设计中互连线引起的所有问题，它主要研究互连线的电气特性参数与数字信号的电压电流波形相互作用后，如何影响到产品性能的问题。主要表现在对时序的影响、信号振铃、信号反射、近端串扰、远端串扰、开关噪声、非单调性、地弹、电源反弹、衰减、容性负载、电磁辐射、电磁干扰等。信号完整性问题的根源在于信号上升时间的减小。即使布线拓扑结构没有变化，如果采用了信号上升时间很小的IC芯片，现有设计也将处于临界状态或者停止工作。下面谈谈几种常见的信号完整性问题。反射：图1显示了信号反射引起的波形畸变。看起来就像振铃，拿出你制作的电路板，测一测各种信号，比如时钟输出或是高速数据线输出，看看是不是存在这种波形。如果有，那么你该对信号完整性问题有个感性的认识了，对，这就是一种信号完整性问题。很多硬件工程师都会在时钟输出信号上串接一个小电阻，至于为什么，他们中很多人都说不清楚，他们会说，很多成熟设计上都有，照着做的。或许你知道，可是确实很多人说不清这个小小电阻的作用，包括很多有了三四年经验的硬件工程师，很惊讶么?可这确实是事实，我碰到过很多。其实这个小电阻的作用就是为了解决信号反射问题。而且随着电阻的加大，振铃会消失，但你会发现信号上升沿不再那么陡峭了。这个解决方法叫阻抗匹配，奥，对了，一定要注意阻抗匹配，阻抗在信号完整性问题中占据着极其重要的

《信号分析与处理》复习总结

信号是带有信息（如语音、音乐、图象、数据等）的随时间（和空间）变化的物理或物理现象，其图象称为信号的波形。信号是消息的表现形式，消息则是信号的具体内容。分类：根据不同分类原则，信号可分为：连续时间信号与离散时间信号；确定信号与随机信号；周期信号和非周期信号；功率信号与能量信号等等反因果信号：若当t ≥0时，f (t )=0;当t <0时，f (t )≠0. 系统：由若干相互作用和相互依赖的事物组合而成的具有特定功能的整体。 ???????=???≠=∞=?∞ ∞ -1)()0( 0)0( )(dt t t t t δδ()()t t δδ-= ()t δ为偶对称函数 1()d 2j t t e ωδωπ ∞-∞= ?——()t δ的逆傅立叶变换 ()()d ()() t x t t t t x t t t δε-∞ -=-?) ()()()(000t t t x t t t x -=-δδ)(| |1 )(t a at δδ= )(t δ'是奇对称函数 ) ()(, 0)(t d d t δττδττδ='='? ? ∞ -∞ ∞ -离散时间单位: 0()(), ()()(1) m n n m n n n εδδεε+∞ ==-=--∑稳定性 ∑? +∞-∞ =∞ +∞ -∞ <∞-=-? -z z z z n Z ε

测试信号分析与处理作业实验一二

王锋实验一：利用FFT 作快速相关估计一、实验目的 a.掌握信号处理的一般方法，了解相关估计在信号分析与处理中的作用。 b.熟悉FFT算法程序；熟练掌握用FFT作快速相关估计的算法。 c.了解快速相关估计的谱分布的情况。二、实验内容 a.读入实验数据[1]。 b.编写一利用FFT作相关估计的程序[2]。 c.将计算结果表示成图形的形式，给出相关谱的分布情况图。注[1]：实验数据文件名为“Qjt.dat”。实验数据来源：三峡前期工程 “覃家沱大桥” 实测桥梁振动数据。实验数据采样频率：50Hz。可从数据文件中任意截取几段数据进行分析，数据长度N 自定。注[2]：采用Matlab 编程。三、算法讨论及分析算法为有偏估计，利用FFT计算相关函数 Step 1: 对原序列补N个零，得新序列x2N(n) Step2: 作FFT[x2N(n)]得到X2N(k) Step 3: 取X2N(k)的共轭，得 Step 4: 作 Step 5: 调整与的错位。四、实验结果分析 1. 该信号可以近似为平稳信号么？可以近似为平稳信号，随机过程的统计特性不随样本的采样时刻而发生变化。取N=8192，分别取间隔m=500,m=700,m=1000,所得到的均值均为0.5366，方差为47369，与时间无关。

图1-1 自相关函数图 (上图表示的R0，下图为调整后的R0) 2. 该信号是否具有周期性，信噪比如何？ >> load Qjt.dat; %加载数据 N=32768; %数据长度 i=1:1:N; %提取数据 plot(i,Qjt(i)); 抛去几个极值点，从图1-2可以看出，数据具有一定的周期性，杂音比较少，说明信噪比较高。图1-2 数据图

现代信号处理大作业题目+答案

研究生“现代信号处理”课程大型作业（以下四个题目任选三题做） 1. 请用多层感知器（MLP ）神经网络误差反向传播（BP ）算法实现异或问题（输入为[00;01;10;11]X T =，要求可以判别输出为0或1），并画出学习曲线。其中，非线性函数采用S 型Logistic 函数。 2. 试用奇阶互补法设计两带滤波器组(高、低通互补)，进而实现四带滤波器组；并画出其频响。滤波器设计参数为：F p ＝1.7KHz , F r ＝2.3KHz , F s ＝8KHz , A rmin ≥70dB 。 3. 根据《现代数字信号处理》（姚天任等，华中理工大学出版社，2001）第四章附录提供的数据(pp.352-353)，试用如下方法估计其功率谱，并画出不同参数情况下的功率谱曲线： 1） Levinson 算法 2） Burg 算法 3） ARMA 模型法 4） MUSIC 算法 4. 图1为均衡带限信号所引起失真的横向或格型自适应均衡器(其中横向FIR 系统长M =11), 系统输入是取值为±1的随机序列)(n x ，其均值为零；参考信号)7()(-=n x n d ；信道具有脉冲响应： 1 2(2)[1cos( )]1,2,3()20 n n h n W π-?+=?=???其它式中W 用来控制信道的幅度失真(W = 2～4, 如取W = 2.9,3.1,3.3,3.5等)，且信道受到均值为零、方差001.02 =v σ(相当于信噪比为30dB)的高斯白噪声)(n v 的干扰。试比较基于下列几种算法的自适应均衡器在不同信道失真、不同噪声干扰下的收敛情况(对应于每一种情况,在同一坐标下画出其学习曲线): 1）横向/格-梯型结构LMS 算法 2）横向/格-梯型结构RLS 算法并分析其结果。

西电随机信号分析大作业

随机信号分析大作业学院：电子工程学院班级：021151 学号：02115037 姓名：隋伟哲

第一题：设有随机信号X（t）=5cos(t+a),其中相位a是在区间（0,2π）上均匀分布的随机变量，使用Matlab编程产生其三个样本函数。解：源程序如下： clc;clear; C=2*pi*rand(1,3);%在[0,2π]产生均匀分布的相位角 t=1:.1:80; y1=5*cos(t+C(1)); %将产生的随机相位角逐一代入随机过程中 y2=5*cos(t+C(2)); %将产生的随机相位角逐一代入随机过程中 y3=5*cos(t+C(3)); %将产生的随机相位角逐一代入随机过程中 plot(t,y1,'r-'); hold on; plot(t,y2,'g--'); hold on; plot(t,y3,'k-'); xlabel('t');ylabel('X(t)'); grid on;axis([0 30 -8 8]); title('随机相位的三条样本曲线'); 产生的三条样本曲线：

第二题：利用Matlab程序设计一正弦型信号加高斯白噪声的复合信号。（1）分析复合信号的功率谱密度、幅度分布特性；（2）分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性；（3）分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。解：设定正选信号的频率为10HZ，抽样频率为100HZ x=sin(2*pi*fc*t)

（1）正弦函数加上高斯白噪声： y=awgn(x,10) y 的幅度分布特性可以通过傅里叶变换得到： Y(jw)=fft(y) y 的功率谱密度： G(w)=Y(jw).*conj(Y(jw)/length(Y(jw))) 随机序列自相关函数的无偏估计公式为： 1 01()()()N m xx n R m x n x n m N m --==+-∑ 01m N ≤≤- （2）复合信号 y 通过RC 积分电路后得到信号y2 通过卷积计算可以得到y2 即：y2= conv2(y,b*pi^-b*t) y2的幅度分布特性可以通过傅里叶变换得到： Y2(jw)=fft(y2) y2的功率谱密度： G2(w)=Y2(jw).*conj(Y2(jw)/length(Y2(jw))) （3）复合信号 y 通过理想滤波器电路后得到信号y3 通过卷积计算可以得到y3 即：y3=conv2(y,sin(10*t)/(pi*t)) y3的幅度分布特性可以通过傅里叶变换得到： Y3(jw)=fft(y3) y3的功率谱密度： G3(w)=Y3(jw).*conj(Y3(jw)/length(Y3(jw)))

测试信号分析与处理作业实验五

王锋实验五：多种功率谱估计的比较一、实验目的 a.了解功率谱估计在信号分析中的作用； b.掌握随机信号分析的基础理论，掌握参数模型描述形式下的随机信号的功率谱的计算方法； c.掌握在计算机上产生随机信号的方法； d.了解不同的功率谱估计方法的优缺点。二、实验准备有三个信号源，分别代表三种随机信号（序列）。信号源1： 123()2cos(2)2cos(2)2cos(2)()x n f n f n f n z n πππ=+++ 其中，1230.08,=0.38,0.40f f f == z(n)是一个一阶 AR 过程，满足方程： ()(1)(1)()z n a z n e n =--+ (1)0.823321a =- e(n)是一高斯分布的实白噪声序列，方差20.1σ= 信号源2和信号源3：都是4阶的AR 过程，它们分别是一个宽带和一个窄带过程，满足方程： ()(1)(1)(2)(2)(3)(3)(4)(4)()x n a x n a x n a x n a x n e n =--------+ e(n)是一高斯分布的实白噪声序列，方差2σ，参数如下：三、实验内容 a. 描绘出这三个实验信号的真实功率谱波形。 b. 在计算机上分别产生这个三个信号，令所得到的数据长度 N= 256 。注意：产生信号的时候注意避开起始瞬态点。例如，可以产生长度为512 的信号序列，然后取后面256 个点作为实验数据。 c. 分别用如下的谱估计方法，对三个信号序列进行谱估计。 1、经典谱估计周期图法自相关法平均周期图法（Bartlett 法）

Welch法（可选每段64 点，重叠32 点，用Hamming 窗）2、现代谱估计 Yule - Walker方程（自相关法）最小二乘法注：阶次p可在3－20之间，由自己给定。四、实验结果分析生成的信号源

西南交大现代信号处理作业

现代信号处理作业 1.(5″)证明下面定理：任何一个无偏估计子方差的下界叫作Cramer-Rao 下界定理：令1(,,)N x x x =为一样本向量，(|)f x θ是x 的条件密度，若?θ是θ的一个无偏估计子，且(|)/f x θθ??存在，则 22 1 ??var()()[ln (|)]E E f x θ θθθθ =-≥? ? 式中?ln (|)()()f x K θθθθθ ?=-?。其中()K θ是θ的某个不包含x 的正函数。 2.(10″)Wiener 滤波是信号处理中最常用和基础的波形估计工具之一，对其在自己研究领域的应用情况进行一个简单综述。 3.(5″)二阶滑动平均过程由 2()()1(1)2(2), {()~(0,)}x n w n b w n b w n w n N σ=+-+- 定义，式中2(0,)N σ表示正态分布，其均值为零、方差为2σ。求x(n)的功率谱。 4.(20″)信号的函数表达式为： ()sin(2100)sin(2300)()sin(2200)()()x t t t A t t dn t n t πππ=++++，其中，A(t)为一随时间变化的随机过程，dn(t)为经过390-410Hz 带通滤波器后的高斯白噪声，n(t)为高斯白噪声，采样频率为1kHz ，采样时间为2.048s 。 (1) 利用现代信号处理知识进行信号的谱估计； (2) 利用现代信号处理知识进行信号的频率提取； (3) 分别利用Wiener 滤波和Kalman 滤波进行去噪； (4) 利用Wigner-Ville 分布分析信号的时频特征。 5.(10″)附件中表sheet1 为某地2008年4月28日凌晨12点至2008年5月4日凌晨12点的电力系统负荷数据，采样时间间隔为1小时，利用ARMA 方法预测该地5月5日的电力系统负荷，并给出预测误差（5月5日的实际负荷数据如表sheet2）。

西电数字信号处理大作业

第二章 2.25 已知线性时不变系统的差分方程为若系统的输入序列x(x)={1,2,3,4,2,1}编写利用递推法计算系统零状态响应的MATLAB程序，并计算出结果。代码及运行结果： >> A=[1,-0.5]; >> B=[1,0,2]; >> n=0:5; >> xn=[1,2,3,4,2,1]; >> zx=[0,0,0];zy=0; >> zi=filtic(B,A,zy,zx); >> yn=filter(B,A,xn,zi); >> figure(1) >> stem(n,yn,'.'); >> grid on;

2.28图所示系统是由四个子系统T1、T2、T3和T4组成的，分别用单位脉冲响应或差分方程描述为 T1：其他 T2: 其他 T3： T4: 编写计算整个系统的单位脉冲响应h(n)，0≤n≤99的MATLAB程序，并计算结果。代码及结果如下： >> a=0.25;b=0.5;c=0.25; >> ys=0; >> xn=[1,zeros(1,99)]; >> B=[a,b,c]; >> A=1; >> xi=filtic(B,A,ys); >> yn1=filter(B,A,xn,xi); >> h1=[1,1/2,1/4,1/8,1/16,1/32]; >> h2=[1,1,1,1,1,1]; >> h3=conv(h1,h2); >> h31=[h3,zeros(1,89)]; >> yn2=yn1+h31; >> D=[1,1];C=[1,-0.9,0.81]; >> xi2=filtic(D,C,yn2,xi); >> xi2=filtic(D,C,ys); >> yn=filter(D,C,yn2,xi); >> n=0:99; >> figure(1) >> stem(n,yn,'.'); >> title('单位脉冲响应'); >> xlabel('n');ylabel('yn');

哈工大测试大作业——信号的分析与系统特性——锯齿波

1 题目：写出下列信号中的一种信号的数学表达通式，求取其信号的幅频谱图（单边谱和双边谱）和相频谱图，若将此信号输入给特性为传递函数为)(s H 的系统，试讨论信号参数的取值，使得输出信号的失真小。（选其中一个信号） 000 2=tan ,=45,=1w 2K T s T π ααπ= =假设锯齿波的斜取周期，则圆周率，A=1 2 幅频谱和相频谱 00()(+nT )(

所以0001111 (t)=(sin(w t)+sin(2w t)+sin(3w t)+223 w π-…) 转换为复指数展开式的傅里叶级数： 0000000-20 2 1-0 --1 00-02222 0001= (t)e =e 11 =e e |11 = e (2) T jnw t T n jnw t jnw t jnw t jnw t c w dt T t dt t jnw jnw jnw n w n w w π-??-+? ???+-=? ? 其中当n=0时，01 = =22 A c ，0=0? ； =1,2,3,n ±±±当… 时， 111 222n n c A n π=== ， 1,2,32 =1,2,32 n n n π ?π?=??? ?-=---?? 等等用Matlab 做出其双边频谱图 1锯齿波双边幅频谱 A = 1 T0 = 1

南邮现代信号处理最后大作业4道题目(含答案)

南邮研究生“现代信号处理”期末课程大作业（四个题目任选三题做） 1. 请用多层感知器（MLP ）神经网络误差反向传播（BP ）算法实现异或问题（输入为[00;01;10;11]X T =，要求可以判别输出为0或1），并画出学习曲线。其中，非线性函数采用S 型Logistic 函数。 2. 试用奇阶互补法设计两带滤波器组(高、低通互补)，进而实现四带滤波器组；并画出其频响。滤波器设计参数为：F p ＝1.7KHz , F r ＝2.3KHz , F s ＝8KHz , A rmin ≥70dB 。 3. 根据《现代数字信号处理》（姚天任等，华中理工大学出版社，2001）第四章附录提供的数据(pp.352-353)，试用如下方法估计其功率谱，并画出不同参数情况下的功率谱曲线： 1） Levinson 算法 2） Burg 算法 3） ARMA 模型法 4） MUSIC 算法 4. 图1为均衡带限信号所引起失真的横向或格型自适应均衡器(其中横向FIR 系统长M =11), 系统输入是取值为±1的随机序列)(n x ，其均值为零；参考信号)7()(-=n x n d ；信道具有脉冲响应： 1 2(2)[1cos( )]1,2,3()20 n n h n W π-?+=?=???其它式中W 用来控制信道的幅度失真(W = 2～4, 如取W = 2.9,3.1,3.3,3.5等)，且信道受到均值为零、方差001.02 =v σ(相当于信噪比为30dB)的高斯白噪声)(n v 的干扰。试比较基于下列几种算法的自适应均衡器在不同信道失真、不同噪声干扰下的收敛情况(对应于每一种情况,在同一坐标下画出其学习曲线): 1）横向/格-梯型结构LMS 算法 2）横向/格-梯型结构RLS 算法并分析其结果。

随机信号分析大作业

随机信号分析实验报告信息25班 2120502123 赵梦然

作业题三：利用Matlab 产生一个具有零均值、单位方差的的高斯白噪声随机序列X(n)，并通过一脉冲响应为 (0.8)(0)0 n n h n else =≥??? 的线性滤波器。 (1) 产生一个具有零均值、单位方差的的高斯白噪声随机序列X(n)，检验其一维概率密度函数是否与理论相符。 (2) 绘出输入输出信号的均值、方差、自相关函数及功率谱密度的图形，讨论输出信号服从何种分布。 (3) 试产生在[-1,+1]区间均匀分布的白噪声序列，并将其替换高斯白噪声通过上述系统。画出此时的输出图形，并观察讨论输出信号服从何种分布。作业要求 (1) 用MATLAB 编写程序。最终报告中附代码及实验结果截图。 (2) 实验报告中必须有对实验结果的分析讨论。提示： (1) 可直接使用matlab 中已有函数产生高斯白噪声随机序列。可使用hist 函数画出序列的直方图，并与标准高斯分布的概率密度函数做对比。 (2) 为便于卷积操作，当N 很大时，可近似认为h(N)=0。卷积使用matlab 自带的conv 函数。 (3) 分析均值、方差等时，均可使用matlab 现有函数。功率谱密度和自相关函数可通过傅里叶变换相互获得。傅里叶变换使用matlab 自带的fft 函数。 (4) 作图使用plot 函数。

一、作业分析：本题主要考察的是加性高斯白噪声相关问题，因此构造一个高斯白噪声十分重要，故在本题中使用randn函数随机生成一个个符合高斯分布的数据，并由此构成高斯白噪声；而且由于白噪声是无法完全表示的，故此根据噪声长度远大于信号长度时可视为高斯白噪声，构造了一个长度为2000的高斯白噪声来进行试验。二、作业解答：（1）matlab程序为： x-1000:1:1000; k=1*randn(1,length(x));% 生成零均值单位方差的高斯白噪声。 [f,xi]=ksdensity(x);%利用ksdensity函数估计样本的概率密度。 subplot(1,2,1); plot(x,k); subplot(1,2,2); plot(xi,f); 实验结果为：