一种新的定点DSP Q格式数乘法运算方法

中南大学数字信号处理实验三

实验报告 实验名称用双线性变换法设计IIR数字滤波器课程名称数字信号处理 姓名成绩 班级学号 日期 2014年5月24号地点综合实验楼机房备注：

1.实验目的 （1）熟悉用双线性变换法设计IIR 数字滤波器的原理与方法； （2）掌握数字滤波器的计算机仿真方法； （3）通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。 2.实验环境 应用MATLAB 6.5软件 操作系统：windows XP 3.实验内容及原理 （1）用双线性变换法设计一个巴特沃斯低通IIR 数字滤波器。设计指标参数为：在通带内截止频率低于0.2π时，最大衰减小于1dB ；在阻带内[0.3π,π]频率区间上，最小衰减大于15dB 。 （2）以0.02π为采样间隔，打印出数字滤波器在频率区间[0，π/2]上的幅频响应特性曲线。 （3）用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，并分别打印出滤波前后的心电图信号波形图，观察总结滤波作用与效果。 教材例中已求出满足本实验要求的数字滤波系统函数： ∏==3 1)()(k k z H z H , 3,2,1,1)21()(2121=--++=----k z C z B z z A z H k k k 式中 A=0.09036, 2155 .0,9044.03583 .0,0106.17051 .0,2686.1332211-==-==-==C B C B C B 4.实验结果 心电图信号采样序列 0510152025 303540455055-100-50 50 n x (n ) 心电图信号采样序列x(n)

一级滤波后的心电图信号： 01020 30405060-100-80 -60 -40 -20 20 40 n y 1(n ) 一级滤波后的心电图信号 二级滤波后的心电图信号： 01020 30405060-100-80 -60 -40 -20 20 40 n y 2(n ) 二级滤波后的心电图信号 三级滤波后的心电图信号： 0102030 405060-80-60 -40 -20 20 40 n y 3(n )三级滤波后的心电图信号

数字信号处理实验一

实验一 离散时间信号分析 班级 信息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）完成实验报告。 五、实验结果 六、实验总结

中南大学RFID课程设计报告

CENTRAL SOUTH UNIVERSITY 课程设计报告 课程： RFID课程设计 班级：物联网工程1201班 学号： 0909120316 姓名：王兆岳 指导教师：李刚 日期： 2015年4月25日

第一节课程设计选题 (1) 1.1选题背景 (1) 1.2课程设计目标 (1) 1.3课程设计使用的相关语言及数据库 (2) 1.4测试环境 (2) 第二节总体设计 (2) 2.1处理流程概要 (2) 2.2总体架构设计 (3) 2.3总体处理流程 (4) 第三节 PC端具体设计 (4) 3.1PC端模块划分 (4) 3.2出入库控制模块 (5) 3.3信息查询模块 (6) 3.4账号注册模块 (8) 3.5充值缴费模块 (8) 3.6硬件通讯中间件 (10) 第四节移动端具体设计 (11) 4.1剩余车位展示 (11) 4.2停车场线路导航 (12) 4.3个人记录、余额查询 (13) 第五节主要算法 (13) 6.1避免刷卡同时激活入库和出库 (13) 6.2多张卡同时在区域内时的屏蔽 (14) 6.3屏蔽偶发错误 (15) 第六节实验总结 (15)

第一节课程设计选题 1.1选题背景 近几年随着我国高速发展，我国的机动车保有量也在不断攀升，因此楼宇、社区和商业区构建停车场及管理系统就显得十分迫切，构建一套包含车辆进出、停车泊位、缴费结算、资料查询、信息提示等功能的相对完善的管理系统，已成为停车场管理部门的共同愿望，同时由于传统停车场并没有与互联网实现对接，经常造成停车位的浪费或是由于驾驶员不能及时获知停车位已满的消息而导致能源的极大浪费、加剧交通拥堵的状况，基于此我选择停车场管理系统作为本次RFID课程设计的题目。 1.2课程设计目标 在本方案中，效率、正确率、信息的整合、以及便捷性是重点追求的目标。 效率读取后数据应及时进行处理，并写入数据库备查 正确率保证每次读取信息的准确性，避免“漏读”或“重读” 信息的整合不同功能模块要实现良好的整合 便捷性尽可能减少人员手动操作，尽量实现自动化

数字信号处理实验作业

实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。

数字信号处理实验

实验一 离散傅里叶变换（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）

(完整版)《中南大学数字信号处理》2014试卷及答案

中南大学考试试卷 2013-- 2014学年 下 学期期末考试试题 时间100分钟 数字信号处理 课程 48 学时 3 学分 考试形式： 闭 卷 专业年级： 电子信息、通信2012级 总分100分，占总评成绩 70 % 注：此页不作答题纸，请将答案写在答题纸上 一、填空题（本题20分，每空2分） 1. 系统稳定的充要条件是系统的单位脉冲响应满足: ∞<∑+∞ -∞=|)(|n n h 。 （p17） 2.若()a x t 是频带宽度有限的，要想抽样后()()a x n x nT =能够不失真地还原出原始信号()a x t ，则抽样频率必须 大于或等于 两倍信号谱的最高频率，这就是奈奎斯特抽样定理。P24 3. 如果序列)(n x 是一长度为64点的有限长序列)630(≤≤n ，序列)(n h 是一长度为128点的有限长序列)1270(≤≤n ，记)()()(n h n x n y *=（线性卷积），则)(n y 为 64+128-1＝191 点的序列，如果采用基FFT 2算法以快速卷积的方式实现线性卷积，则FFT 的点数至少为 256 点。P12、p111 4. 设序列()x n 傅立叶变换为()jw X e ，则0()x n n -（0n 为任意实整数）的傅立叶变换是 0)(jwn jw e e X -? 。P35 5. 序列()(3)x n n δ=-的傅里叶变换是 3jw e - 。P35 6.某DFT 的表达式是1 0()()N kn N n X k x n W -==∑，则变换后数字频域上相邻两个频率样点之间的 间隔是 2/N π 。 p76 7.用DFT 对模拟信号进行谱分析，会有 频谱混叠、截断效应、栅栏效应 三种误差来源。 P103 二、单项选择题(10分，每题2分) 1. 序列()(1)n x n a u n =---，则()X z 的收敛域为（ A ）。P48列 2.5.4 A. ||||z a < B. ||||z a ≤ C. ||||z a > D. ||||z a ≥ 2．下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？( D )p11 A.5()()y n x n = B.()()(2)y n x n x n =+ C.()()2y n x n =+ D. 2 ()()y n x n = 3. 直接计算N 点DFT 所需的复数乘法次数与( B )成正比。P110 A.N B.N 2 C.N 3 D.Nlog 2N 4.ZT[2()]n u n --=__B____。P46，例2.5.1

数字信号处理实验作业

实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示：

原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。

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

《测试技术与信号处理》习题答案 第二章 信号分析基础 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

数字信号处理实验答案完整版

数字信号处理实验答案 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),绘出四个周期。

我对信息与信号处理的理解

我对信息与信号处理的理解 --电气1031班肖斯诺 第一次认识到信号这个概念是在小学自然课上的一个小实验：用两个杯 子和一根很长的线远距离聊天。我现在都还记得当时我们几个最先做完的小伙 伴得意的表情，像是吃到了最甜的糖果。后来我慢慢知道，其实信号充斥着我 周围的每一个角落，电视，空调，微波炉等等……信号几乎无处不在。而第一 次深入了解和学习信号是在大一的这门信息与信号课程上，接下来说说我对这 门课程的理解吧。 先说说何为信号，信号是运载消息的工具，是消息的载体。从广义上讲， 它包含光信号、声信号和电信号等。例如，古代人利用点燃烽火台而产生的滚 滚狼烟，向远方军队传递敌人入侵的消息，这属于光信号；当我们说话时，声 波传递到他人的耳朵，使他人了解我们的意图，这属于声信号；遨游太空的各 种无线电波、四通八达的电话网中的电流等，都可以用来向远方表达各种消息，这属电信号。人们通过对光、声、电信号进行接收，才知道对方要表达的消息。 总的来说，信息的具体表现形式是信号，信息是信号包含的内容。没有信息，信号将毫无意义，这是两个分不开，却又完全不同的概念。 为了充分地获取信息和有效利用信息，必须对信号进行分析和处理。其中 包括两个方面，即信号分析和信息处理。而信息处理则指按某种需要或目的， 对信号进行特定的加工，操作或修改。信号处理涉及的领域非常广泛，包括信 号滤波，信号中的干扰/噪声抑制或滤除、信号平滑、信号锐减、信号增强、信 号的数字化、信号的恢复和重建、信号的编译和译码、信号的调制和调解、信 号加密和解密、信号均衡或校正、信号的特征提取、信号的辨识或目标识别、

信息融合及信号的控制，等等。 这是现代信号处理的过程，而古往至今信息处理是经过了多年的演变才有了今天对信号如此多变的应用。 概括说来，信息与信号处理大致经历了一个这样的发展过程：肢体语言信息语言文字信息远程通讯与信息处理模拟信号与信息处理数字信号与信息处理。 信号其实在人类之前就有了，蜜蜂跳舞就是一种信号，蜜蜂们通过跳舞产生信号，让自己的同伴了解到自己所要表达的信息。而当人类诞生以后，信号的世界才变得丰富多彩。 古时候，大概还在石器时期的时候，类人猿通过吼叫以及各种肢体语言在种群生活中向其他类人猿表达自己的想法，后来，随着人类祖先的不断进化，开始使用各种工具，人类的生活中不只只有寻找食物以谋求生存，于是，语言出现了。这是人类进化史上的一大步。人类文明史上的一个伟大的里程碑就是语言的诞生。语言的诞生让人类之间的信息交流变得更方便，人类文明也因此不断地进化，再之后，人类又发明了沉默的语言——文字，文字的产生让人类的学习能力增加，文字这种信号形式让信息可以长时间的保存，人类的技能和知识通过文字的形式得到保存，后人通过前人保存下来的信息可以直接得到前人总价下来的知识精华，并以此获取更多的知识，这让人类在相对来说短短的几千年来成为世界当之无愧的统治者。 自工业革命以后200来年，人类的文明又得到了一步巨大的跨越，科学知

数字信号处理实验1认识实验

实验1认识实验-MATLAB语言上机操作实践 一、实验目的 ㈠了解MATLAB语言的主要特点、作用。 ㈡学会MATLAB主界面简单的操作使用方法。 ㈢学习简单的数组赋值、运算、绘图、流程控制编程。 二、实验原理 ㈠简单的数组赋值方法 MATLAB中的变量和常量都可以是数组（或矩阵）,且每个元素都可以是复数。 在MATLAB指令窗口输入数组A=[1 2 3；4 5 6；7 8 9]，观察输出结果。然后，键入：A（4，2）= 11 键入：A (5，：) = [-13 -14 -15] 键入：A（4，3）= abs (A(5，1)) 键入：A ([2，5]，：) = [ ] 键入：A/2 键入：A (4，：) = [sqrt(3) (4+5)/6*2 –7] 观察以上各输出结果。将A式中分号改为空格或逗号，情况又如何？请在每式的后面标注其含义。 2．在MATLAB指令窗口输入B=[1+2i，3+4i；5+6i ，7+8i], 观察输出结果。 键入：C=[1，3；5，7]+[2，4；6，8]*i，观察输出结果。 如果C式中i前的*号省略，结果如何？ 键入：D = sqrt (2+3i) 键入：D*D 键入：E = C’, F = conj(C), G = conj(C)’ 观察以上各输出结果, 请在每式的后面标注其含义。 3．在MATLAB指令窗口输入H1=ones(3,2)，H2=zeros(2,3)，H3=eye(4)，观察输出结果。 ㈡、数组的基本运算 1．输入A=[1 3 5]，B= [2 4 6]，求C=A+B，D=A-2，E=B-A 2．求F1=A*3，F2=A.*B，F3=A./B，F4=A.\B, F5=B.\A, F6=B.^A, F7=2./B, F8=B.\2 *3．求B'，Z1=A*B’，Z2=B’*A 观察以上各输出结果,比较各种运算的区别，理解其含义。 ㈢、常用函数及相应的信号波形显示 例1：显示曲线f(t)=2sin(2πt)，（t>0） ⅰ点击空白文档图标（New M-file），打开文本编辑器。 ⅱ键入：t=0:0.01:3; (1) f=2*sin(2*pi*t); (2) plot(t,f); title(‘f（t）-t曲线’)； xlabel（‘t’），ylabel（‘f（t）’）；

2007年华中科技大学计算机专业基础综合(计算机组成原理)真题试卷

2007年华中科技大学计算机专业基础综合（计算机组成原理） 真题试卷 (总分：62.00，做题时间：90分钟) 一、填空题(总题数：13，分数：26.00) __________________________________________________________________________________________ 正确答案：(正确答案：1.1101) __________________________________________________________________________________________ 正确答案：(正确答案：111111) 3.磁盘的平均存取时间是指______，该时间由______时间和平均等待时间组成，对于7200转的磁盘而言，其平均等待时间约为______ms(取整数)。 __________________________________________________________________________________________ 正确答案：(正确答案：CPU发出读写命令后，磁头从某一起始位置移动到新的位置，以及从盘片上读写信息所需的平均时间、平均找道、4ms) 4.在微指令格式设计过程中，有8个互斥型的微命令被分配到一组，当该组采用编码方法表示时，微指令格式的相关字段至少需要______位。 __________________________________________________________________________________________ 正确答案：(正确答案：4) 5.CPU中保存当前正在执行指令的寄存器是______。 __________________________________________________________________________________________ 正确答案：(正确答案：指令寄存器LR) 6.CPU的功能包括指令的序列控制_______、_______和_______。 __________________________________________________________________________________________ 正确答案：(正确答案：操作控制、时间控制、数据加工) __________________________________________________________________________________________ 正确答案：(正确答案：1100) 8.设计一个64位全并行的ALU需要______片74181和______片______部件74182。 __________________________________________________________________________________________ 正确答案：(正确答案：16、5、先行进位) 9.响应中断必须具备下列几个条件，它们分别是中断源有中断请求、允许接受中断请求和______。 __________________________________________________________________________________________ 正确答案：(正确答案：现行指令执行结束) __________________________________________________________________________________________ 正确答案：( 11.“数据线双向传输”这句话描述了总线的______特性。 __________________________________________________________________________________________

数字信号处理实验三

实验三：离散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

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

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) 从下面的信号波形图中读出其主要参数。

数字信号处理实验4

数字信号处理实验四 第一题结果： （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))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果：

数字信号处理实验八

实验报告 实验名称：FIR数字滤波器设计及应用 课程名称____数字信号处理________ 院系部:电气与电子工程专业班级：信息1002 学生姓名：王萌学号: 11012000219同组人:实验台号: 指导教师：范杰清成绩： 实验日期: 华北电力大学

一、实验目的 加深理解 FIR 数字滤波器的时域特性和频域特性，掌握FIR 数字 滤波器的设计原理与设计方法，以及FIR 数字滤波器的应用。 二、 实验原理 FIR 数字滤波器可以设计成具有线性相位，在数据通信、图像处理、 语音信号处理等实际应用领域得到广泛应用。 M 阶FIR 数字滤波器的系统函数为： FIR 数字滤波器的单位脉冲响应h [k ]是长度为M +1的有限长因果序列。当满足对称条件时，该FIR 数字滤波器具有线性相位。FIR 数字滤波器设计方法主要有窗口法、频率取样法及优化设计法。 MATLAB 中提供的常用FIR 数字滤波器设计函数有： fir1 窗函数法设计FIR 数字滤波器（低通、高通、带通、 带阻、多频带滤波器） fir2 频率取样法设计FIR 数字滤波器：任意频率响应 firls FIR 数字滤波器设计：指定频率响应 firrcos 升余弦型 FIR 数字滤波器设计 intfilt 内插FIR 数字滤波器设计 kaiserord 凯塞(Kaiser)窗函数设计法的阶数估计 firpm Parks-McClellan 算法实现FIR 数字滤波器优化设计 firpmord Parks-McClellan 数字滤波器的阶数选择 cremez 复系数非线性相位FIR 等波纹滤波器设计 1、 窗口法设计FIR 数字滤波器 fir1函数可以很容易地实现FIR 数字滤波器窗口法设计。 可设计低通、高通、带通、带阻滤波器、多频带滤波器。 k M k z k h z H -=∑=][)(0

数字信号处理 实验报告

中南大学信息科学与工程学院 数字信号处理实验报告 学生学院信息科学与工程学院 专业班级通信工程1201 学号 学生姓名____ 指导教师李宏 2014年6月2日

实验一 常见离散时间信号的产生和频谱分析 1）实验指导： 一.实验目的 （1） 熟悉MATLAB 应用环境，常用窗口的功能和使用方法； （2） 加深对常用离散时间信号的理解； （3） 掌握简单的绘图命令； （4） 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离 散信号进行频域分析。 二.实验原理 （1） 常用离散时间信号 a ）单位抽样序列 ???=0 1)(n δ ≠=n n 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： ???=-0 1)(k n δ ≠=n k n b ）单位阶跃序列 ???=0 1 )(n u 00<≥n n c ）矩形序列 ???=01)(n R N 其他1 0-≤≤N n d ）正弦序列 )sin()(?+=wn A n x e ）实指数序列 f ）复指数序列 ()() n x n a u n =

()()jw n x n e σ+= （2）离散傅里叶变换： 设连续正弦信号()x t 为 0()sin()x t A t φ=Ω+ 这一信号的频率为0f ，角频率为002f πΩ=，信号的周期为000 12T f π = = Ω。如果对此连续周期信号()x t 进行抽样，其抽样时间间隔为T ，抽样后信号以()x n 表示，则有0()() sin()t nT x n x t A nT φ===Ω+，如果令w 为数字频率，满足 0000 1 2s s f w T f f π=Ω=Ω=，其中s f 是抽样重复频率，简称抽样频率。 为了在数字计算机上观察分析各种序列的频域特性，通常对)(jw e X 在 []π2,0上进行M 点采样来观察分析。 对长度为N 的有限长序列x(n), 有 ∑-=-=1 )()(N n n jw jw k k e n x e X 其中 1,,1,02-== M k k M w k ，π 通常M 应取得大一些，以便观察谱的细节变化。取模|)(|k jw e X 可绘出幅频特性曲线。 （3）用DFT 进行普分析的三种误差 三种误差：混叠现象、泄露现象、栅栏效应 a) 混叠现象 当采样频率小于两倍信号（这里指是信号）最大频率时，经过采样就会发生频谱混叠，这使得采样后的信号序列频谱不能真实地反映原信号的频谱。所以在利用DFT 分析连续信号的频谱时，必须注意这一问题。避免混叠现象的唯一方法是保证采样速率足够高，使频谱交叠现象不致出现。也就是说，在确定采样频率之前，必须对信号的性质有所了解，一般在采样前，信号通过一个防混叠低通滤波器。

数字信号处理第二章上机题作业

数字信号处理作业实验题报告 第一章16.(1) 实验目的： 求解差分方程所描述的系统的单位脉冲响应和单位阶跃响应。 实验要求： 运用matlab求出y(n)=0.6y(n-1)-0.08y(n-2)+x(n)的单位脉冲响应和单位阶跃响应的示意图。 源程序： B1=1;A1=[1, -0.6, 0.08]; ys=2; %设差分方程 xn=[1, zeros(1, 20)]; %xn=单位脉冲序列，长度N=31 xi=filtic(B1, A1, ys); hn1=filter(B1, A1, xn, xi); %求系统输出信号hn1 n=0:length(hn1)-1; subplot(2, 1, 1);stem(n, hn1, '.') title('单位脉冲响应'); xlabel('n');ylabel('h(n)') xn=ones(1, 20); sn1=filter(B1, A1, xn, xi); %求系统输出信号sn1 n=0:length(sn1)-1; Subplot(2, 1, 2); stem(n, sn1, '.') title('单位阶跃响应'); xlabel('n'); ylabel('s(n)')

运行结果： 实验分析： 单位脉冲响应逐渐趋于0，阶跃响应保持不变，由此可见，是个稳定系统。

第二章31题 实验目的： 用matlab判断系统是否稳定。 实验要求： 用matlab画出系统的极，零点分布图，输入单位阶跃序列u(n)检查系统是否稳定。 源程序： A=[2, -2.98, 0.17, 2.3418, -1.5147]; B=[0, 0, 1, 5, -50]; subplot(2,1,1); zplane(B,A); %求H(z)的极点 p=roots(A); %求H(z)的模 pm=abs(p); if max(pm)<1 disp('系统因果稳定'), else,disp('系统因果不稳定'),end un=ones(1,800); sn=filter(B, A, un); n=0:length(sn)-1; subplot(2, 1, 2);plot(n, sn) xlabel('n');ylabel('s(n)')