文档视界 最新最全的文档下载
当前位置:文档视界 › 数字信号处理课程的实验教学设计

数字信号处理课程的实验教学设计

数字信号处理课程的实验教学设计
数字信号处理课程的实验教学设计

数字信号处理期末实验-语音信号分析与处理

语音信号分析与处理 摘要 用MATLAB对语音信号进行分析与处理,采集语音信号后,在MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声,对加入噪声的信号进行频谱分析,设计合适的滤波器滤除噪声,恢复原信号。 数字滤波器是数字信号处理的基础,用来对信号进行过滤、检测和参数估计等处理。IIR数字滤波器最大的优点是给定一组指标时,它的阶数要比相同组的FIR滤波器的低的多。信号处理中和频谱分析最为密切的理论基础是傅立叶变换(FT)。离散傅立叶变换(DFT)和数字滤波是数字信号处理的最基本内容。 关键词:MATLAB;语音信号;加入噪声;滤波器;滤波 1. 设计目的与要求 (1)待处理的语音信号是一个在20Hz~20kHz频段的低频信号。 (2)要求MATLAB对语音信号进行分析和处理,采集语音信号后,在MATLAB平台进行频谱分析;并对所采集的语音信号加入干扰噪声,对加入噪声的信号进行频谱分析,设计合适的滤波器进行滤除噪声,恢复原信号。

2. 设计步骤 (1)选择一个语音信号或者自己录制一段语音文件作为分析对象; (2)对语音信号进行采样,并对语音信号进行FFT频谱分析,画出信号的时域波形图和频谱图; (3)利用MATLAB自带的随机函数产生噪声加入到语音信号中,对语音信号进行回放,对其进行FFT频谱分析; (4)设计合适滤波器,对带有噪声的语音信号进行滤波,画出滤波前后的时域波形图和频谱图,比较加噪前后的语音信号,分析发生的变化; (5)对语音信号进行回放,感觉声音变化。 3. 设计原理及内容 3.1 理论依据 (1)采样频率:采样频率(也称采样速度或者采样率)定义了每秒从连续信号中提取并组成离散信号的采样个数,它用赫兹(Hz)来表示。采样频率只能用于周期性采样的采样器,对于非周期采样的采样器没有规则限制。通俗的讲,采样频率是指计算机每秒钟采集多少个声音样本,是描述声音文件的音质、音调,衡量声卡、声音文件的质量标准。采样频率越高,即采样的间隔时间越短,则在单位之间内计算机得到的声音样本数据就越多,对声音波形的表示也越精确。 (2)采样位数:即采样值或取样值,用来衡量声音波动变化的参数。 (3)采样定理:在进行模拟/数字信号的的转换过程中,当采样频率f s.max大于信号中,最高频率f max的2倍时,即:f s.max>=2f max,则采样之后的数字信号完整的保留了原始信号中的信息,一般实际应用中保证采样频率为信号最高频率的5~10倍;采样频率又称乃奎斯特定理。 (4)时域信号的FFT分析:信号的频谱分析就是计算信号的傅立叶变换。连续信号与系统的傅立叶分析显然不便于直接用计算机进行计算,使其应用受到限制。而FFT是一种时域和频域均离散化的变换,适合数值计算,成为用计算机分析离

课题实验课教学设计新部编版与实施过程的实验报告

教师学科教案[ 20 – 20 学年度第__学期] 任教学科:_____________ 任教年级:_____________ 任教老师:_____________ xx市实验学校

课题实验课教学设计与实施过程的实验报告 ——《和奶奶在一起的时候》一课的教学设计与实施 哈尔滨市道外区新一小学赵长清 一、课题自然情况: 1、课题名称及类别: 准备以《如何解决小学中高年级学生作文难问题的试验与对策的研究》为研究课题开展实验与研究工作。(道外区教育科学规划课题)2、课题简介: 中国是一个具有五千年悠久历史和灿烂文化,古老而文明的国度。在中国古代曾涌现出多少文人墨客,留下多少文化财富。为中国文化的发展,乃至世界文化的传承,奠定了坚实的基础。中国古典四大名著与世界十大名著相提并论,受到世人瞩目。然而随着时间的推移,社会发展到今天,中国的文学创作的在世界文学史上的地位受到撼动。翻开世界近代史,获得诺贝尔文学奖的作品中,没有一部是中国本土作家的大作。这可以说是中国文学界的悲哀。像鲁迅、老舍等一代文学大师香消玉损后,中国文学创作便走向了低谷,创作队伍呈现出了青黄不接的现象。这不能不说明中国的教育,特别是语文教学的悲哀。纵观中国近代的教育,长期以来,一直以应试教育为主,培养出来的学生都是应试教育的牺牲品。作品毫无创新意义,失去了创作本身的灵性,甚至有些作品脱离社会实际,无病呻吟,失去了创作

的本身的意义。评价也以分数高低论英雄,这种制度对于学生身心发展和教育教学质量的提高,并没有取得预期的理想的效果。有时甚至扼杀了学生创作的灵感。致使学生失去了创作的欲望,对作文厌恶。创作过程中应付了事,甚至抄袭。现在学生作文难已经成为了中国学生学习的杀手,严重制约了教育教学质量的提升,阻碍了我国的文学创作史的发展。那么如何解决当代小学阶段中高年级学生作文难的问题便显得尤为重要,针对这一现状。我们经过了大量的查阅国内外有关资料,及对我校三至五学年部分学生进行的随机抽取的调查问卷的分析,我们认为学生作文难主要的原因有: ⑴、学生虽然生活条件优越,见闻丰富广博,可学生没有养成良好的观察习惯,经历的事都是过眼云烟,毫无印象,没有积累意识,缺乏教师的有效引导,没有形成观察积累的习惯。 ⑵、部分学生虽然也有观察积累的意识,虽然经历了艰苦的创作过程,却得不到家长、教师的赏识。久而久之就对作文失去了信心,失去了创作的激情。 ⑶、还有一少部分学生虽有着极高的创作热情,易于动笔,但由于他们受到生活阅历,写作经验,创作技巧的制约,习作水平自然可想而知了。 根据掌握的有关数据,我们拟定以《如何解决小学中高年级学生作文难问题的试验与对策的研究》为研究课题开展实验与研究工作,通过行之有效的教学策略,解决当代小学中高年级学生作文难的问题。通过今天的实验课的调研分析,听取有关专家与同仁们的意见与

数字信号处理实验一

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

数字信号处理实验教案剖析

数字信号处理实验教案 信息工程学院-通信工程教研室

数字信号处理是一门理论和实际密切结合的课程,为深入掌握课程内容,最好在学习理论的同时,做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论,而且能锻炼同学们的独立解决问题的能力。本讲义在第三版的基础上编写了五个实验,前2个实验属基础性的验证性实验,第3、4、5个实验属基本应用综合性实验。 实验一离散时间信号的MA TLAB实现 实验二线性卷积与循环卷积的原理及应用 实验三频率采样定理 实验四离散系统的因果性和稳定性及频率响应特性 实验五基于MATLAB的快速傅里叶变换 根据教学进度,理论课结束后进行相关实验。

实验一时域离散信号的产生 一实验目的 (1)了解常用的时域离散信号及其特点。 (2)掌握MATLAB产生常用时域离散信号的方法。 二实验内容 (1) 编写程序,产生下列离散序列: A.f(n)=δ(n) (-3> n1=-3;n2=4;n0=0; %在起点n1、终点n2的范围内,于n0处产生冲激 >> n=n1:n2; %生成离散信号的时间序列 >> x=[n==n0]; %生成离散信号x(n) >> stem(n,x,'filled'); %绘制杆状图,且圆心处用实心圆表示 >> title('单位脉冲序列'); >> xlabel('时间(n)');ylabel('幅度x(n)'); 在上述语句输入完成之后,敲击回车键,弹出图形窗口,显示出如下图形,即已经满足题干所述条件,产生了f(n)=δ(n),(-3> n1=16;a=0.1;w=1.6*pi; >> n=0:n1; >> x=exp((a+j*w)*n); >>subplot(2,1,1),stem(n,real(x)); %在指定位置描绘图像 >> title('复指数序列的实部'); >> subplot(2,1,2),stem(n,imag(x)); >> title('复指数序列的虚部'); 在上述语句输入完成之后,敲击回车键,弹出图形窗口,显示出如下图形,即已经满足题干所述条件,产生了f(n)=e(0.1+j1.6π)n,(0> f=50;Um=1;nt=2; %输入信号频率、振幅、显示周期 >> N=16;T=1/f; %N为信号一个采样周期的采样点数,T为信号周期 >> dt=T/N; %采样时间间隔 >> n=0:nt*N-1; %建立离散时间的时间序列 >> tn=n*dt; %确定时间序列样点在时间轴上的位置 >> f=Um*sawtooth(2*f*pi*tn)+1;

数字信号处理实验作业

实验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)

数字信号处理实验

实验一 自适应滤波器 一、实验目的 1、掌握功率谱估计方法 2、会用matlab 对功率谱进行仿真 二、实验原理 功率谱估计方法有很多种,一般分成两大类,一类是经典谱估计;另一类是现代谱估计。经典谱估计可以分成两种,一种是BT 法,另一种是周期法;BT 法是先估计自相关函数,然后将相关函数进行傅里叶变换得到功率谱函数。相应公式如下所示: ||1 *0 1 ?()()()(11) ??()(12) N m xx n jwn BT xx m r m x n x n m N P r m e --=∞ -=-∞ =+-=-∑ ∑ 周期图法是采用功率谱的另一种定义,但与BT 法是等价的,相应的功率谱估计如下所示: 2 1 1? ()()01 (13)N jw jwn xx n P e x n e n N N --==≤≤--∑ 其计算框图如下所示: ) (jw xx e ∧ 图1.1周期图法计算用功率谱框图

由于观测数据有限,所以周期图法估计分辨率低,估计误差大。针对经典谱估计的缺点,一般有三种改进方法:平均周期图法、窗函数法和修正的周期图平均法。 三、实验要求 信号是正弦波加正态零均值白噪声,信噪比为10dB,信号频率为2kHZ,取样频率为100kHZ。 四、实验程序与实验结果 (1)用周期图法进行谱估计 A、实验程序: %用周期法进行谱估计 clear all; N1=128;%数据长度 N2=256; N3=512; N4=1024; f=2;%正弦波频率,单位为kHZ fs=100;%抽样频率,单位为kHZ n1=0:N1-1; n2=0:N2-1; n3=0:N3-1; n4=0:N4-1; a=sqrt(20);%由信噪比为10dB计算正弦信号的幅度

数字信号处理实验作业

实验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的围,频谱出现了镜像对称的部分。

开题报告.有效教学各学科教学设计研究》课题

一、研究背景 促进教师专业化发展是教师教育改革的重大主题之一,也是教师教育研究的核心课题。2006年12月全国教育科学规划领导小组下达了由中央教科所承担的全国教育科学“十一五”规划教育部重点课题“有效推进区域教师专业化发展”,之后,中央教科所在全国遴选课题实验区。2007年4月,合肥市教育局组成课题组积极向中央教科所申报“有效推进合肥地区教师专业化发展”的研究课题,2007年4月13日,该课题经中央教科所“有效推进区域教师专业化发展”总课题组批准立项。为了全面提升我校教师素质,促进教师专业化发展,加快教师教育一体化进程,结合我校教师教育实际,2007年11月,我校整合研究力量,组成有研究能力的课题组,积极向合肥市教育局、瑶海区教育局申报子课题《有效教学各学科教学设计研究》,很快得到了合肥市教育局、瑶海区教育局课题组批准立项,这是上级领导对我校的信任,同时也体现了市、区课题组对我校教师的关心和爱护。 二、研究目的和意义 伴随着我国基础教育新课程的全面实施,基础教育的改革与发展跨入了一个崭新的时代。同时也使我们教师面临着一场严峻的考验。能否在思想上、行动上尽快适应新课程的要求,是这次基础教育课程改革成败的关键。因此,切实搞好教师培训工作,促进教师专业化发展,努力建设一支高素质的教师队伍,是完成本次课程改革任务、全面推进素质教育的保证。本课题的研究,就是面对基础教育发展的新形势,深入研究、探索并最终构建适合我校实际、有效促进教师专业

化发展的新思路、新方式,提高整个教师队伍的专业素质,增强广大教师的创新精神和实践能力。 三、研究目标 (一)全面提高我校教师专业素质。通过课题研究,在我校催生一批具有高专业素质的教育教学能手;骨干教师、学科带头人、特级教师等拔尖人才队伍不断加强。 (二)提升学校办学理念。学校的首要目标是促进全体学生的全面发展,形成“教书育人、管理育人、服务育人、环境育人”的氛围;树立科学的发展观、以人为本的教育观、全面发展的质量观;明确学校的科学发展方向和目标定位,要关注每一个学生的发展,关注每一位教师的成长,研究促进教师专业发展的有效途径,努力丰富校园文化,促进学校办学理念、教育教学管理和教育发展水平的再提升。通过课题研究,努力成为学习型学校。 四、研究内容 什么样的教学才是有效教学?有效教学的指标是什么?有效教 学如何设计与实施?有效教学如何评价?有效教学的反思如何进 行?如何重建有效教学的管理制度等等,这些被大家所关注的问题,是每一位教师专业成长的核心问题,在课题研究中要深入研讨。 五、研究方法 (1)文献借鉴法。多角度、多渠道开展对情报资料的比较研究,借鉴已有研究成果和经验教训,找到新的生长点,防止重复研究,避免和少走弯路。

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

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

数字信号处理上机实验答案(全)1

第十章 上机实验 数字信号处理是一门理论和实际密切结合的课程,为深入掌握课程内容,最好在学习理论的同时,做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论,而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验,前五个实验属基础理论实验,第六个属应用综合实验。 实验一 系统响应及系统稳定性。 实验二 时域采样与频域采样。 实验三 用FFT 对信号作频谱分析。 实验四 IIR 数字滤波器设计及软件实现。 实验五 FIR 数字滤波器设计与软件实现 实验六 应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度,安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一;在学习完第三、四章后作实验二和实验三;实验四IIR 数字滤波器设计及软件实现在。学习完第六章进行;实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行;实验六为综合实验,在学习完本课程后再进行。 10.1 实验一: 系统响应及系统稳定性 1.实验目的 (1)掌握 求系统响应的方法。 (2)掌握时域离散系统的时域特性。 (3)分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中,描写系统特性的方法是差分方程和单位脉冲响应,在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应,本实验仅在时域求解。在计算机上适合用递推法求差分方程的解,最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积,求出系统的响应。 系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性,包括观察系统的暂态响应和稳定响应。 系统的稳定性是指对任意有界的输入信号,系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。 实际中检查系统是否稳定,不可能检查系统对所有有界的输入信号,输出是否都是有界输出,或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列,如果系统的输出趋近一个常数(包括零),就可以断定系统是稳定的[19]。系统的稳态输出是指当∞→n 时,系统的输出。如果系统稳定,信号加入系统后,系统输出的开始一段称为暂态效应,随n 的加大,幅度趋于稳定,达到稳态输出。 注意在以下实验中均假设系统的初始状态为零。 3.实验内容及步骤 (1)编制程序,包括产生输入信号、单位脉冲响应序列的子程序,用filter 函数或conv 函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。 (2)给定一个低通滤波器的差分方程为

数字信号处理实验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)’);

数字信号处理实验三

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

数字信号处理实验教学大纲

《数字信号处理》实验教学大纲 一、大纲说明 课程名称(中文):数字信号处理实验 课程名称(英文):Experiment of Digital Signal Processing 适用专业:电子信息工程及电子信息科学与技术 课程属性:专业限选课 课程性质:非独立设课 课程学时:12 先修课程:高等数学、复变函数与积分变换、信号与系统等 二、实验课教学内容和要求 实验一:系统响应及系统稳定性 实验目的和教学要求 1)掌握求系统响应的方法 2)掌握时域离散系统的时域特性 3)分析、观察及检验系统的稳定性 实验二:时域采样与频域采样 实验目的和教学要求 1)掌握模拟信号采样前后频谱的变化 2)掌握如何选择采样频率才能使采样后的信号不丢失信息 3)掌握频率域采样会引起时域周期化的概念 4)掌握频率域采样定理及其对频率采样点数选择的指导作用 实验三:用FFT对信号作频谱分析 实验目的和教学要求 1)了解用FFT对连续信号和时域离散信号进行谱分析的方法 2)了解可能出现的分析误差及其原因,以便正确应用FFT 实验四:IIR数字滤波器设计及软件实现 实验目的和教学要求 1)熟悉用双线性变换法设计IIR数字滤波器的原理与方法 2)熟悉调用MATLAB信号处理工具箱中滤波器设计函数设计各种IIR数字滤波器3)掌握根据滤波需求确定滤波器指标参数 4)掌握IIR数字滤波器的MATLAB实现方法

实验五:FIR数字滤波器设计及软件实现 实验目的和教学要求 1)掌握用窗函数法设计FIR数字滤波器的原理与方法 2)掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法 3)掌握FIR滤波器的快速卷积实现原理 4)熟悉调用MATLAB函数设计与实现FIR滤波器 实验六:数字信号处理在双音多频拨号系统中的应用 实验目的和教学要求 1)了解音频电话中的拨号信号--双音多频信号 2)了解小型信号处理系统,用数字方法产生模拟信号并进行传输 3)了解双音多频信号的产生方法和检测方法 三、主要仪器设备 微型计算机,MATLAB软件等 四、实验方式与要求 1)实验前,教师需向学生讲清实验的性质、任务、要求、实验守则及实验室安全制度等; 2)该课以验证性实验为主,指导书中给出实验内容与方法,也可由学生自主设计实验方法,实验前学生必需进行预习,设计报告经教师批阅后,方可进入实验室进行实验; 3)实验2~3人一组,在规定的时间内,由学生独立完成,出现问题,教师要引导学生独立分析、解决,不得包办替; 4)任课教师要认真上好每一堂课,实验前清点学生人数,实验中按要求做好学生实验情况及结果记录,实验后认真填写实验开出记录。 五、考核办法 实验部分不单独考试,实验成绩占课程总评的10%。 六、学时分配 本课程计划学时数为48学时,其中理论课36学时,实验课12学时。 实验教学学时分配表

小课题研究教学设计

《圆的认识》教学设计 杨英庄小学六年级王志娟 【教学内容】人教课标版六年级上册第55–58页,练习十四第1-4题。 【教学目标】 1、通过学生的实践活动、合作交流,使学生进一步掌握圆的特征。从而培养学生的动手能力、分析能力、抽象能力。 2、让学生自己尝试用圆规画圆,培养学生的独立探索能力及观察能力。 3、使学生认识到圆在日常生活中的存在和作用,体验数学的价值。 【教学重点】使学生充分掌握圆的特征。 【教学难点】理解并掌握同圆中直径与半径的特征及关系。 【教学准备】圆的认识的课件、圆形的物体、圆形纸片、圆规。 【教学过程】 一、创设情境,感受新知。 1、同学们,你们见过平静的的水面吗?如果我们从上面往下面丢进一颗小石子(播放水纹),你们发现了什么? 2、其实这样的现象在大自然中随处可见,让我们一起来看看。(伴随着优美的音乐,阳光下绽放的向日葵、花丛中五颜六色的鲜花、光折射后形成的美妙光环、圆形的瓷盘及剪纸等画面一一展现在学生的眼前),从这些现象中,你们找到圆了吗?

3、在我们的生活周围,你还知道哪些物体的形状是圆形的?谁来说说? 4、有人说,因为有了圆,我们的世界才变得如此美妙而神奇。今天这节课,就让我们一起走进圆的世界,去认识圆!(板书课题:圆的认识) 二、提出学习问题 师:谁愿意和老师合作在黑板上面一个圆?(同学们争先恐后举手) 师:看着黑板上这么大的一个圆,你有什么问题想说吗? 生: ①我想知道什么是圆? ②圆的各部分名称是什么? ③圆的周长怎么求? ④圆的面积怎么求? ⑤圆的性质、特点都是什么? ⑥我真羡慕老师,也想在本上画一个大圆,我想了解它的画法。 师:你们的思维太开阔了,想了解这么多的内容,我们归纳起来,就是什么是圆,各部分名称……那么现在就运用我们的聪明才智一个 一个攻破它。 三、自主探究、汇报交流 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))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果:

_数字信号处理实验一.

姓名:罗格学号:2012302530084 班级:信安3班日期:2014.3.15 实验1基本离散信号的MATLAB产生和图形显示 实验前言: MATLAB 是一套功能强大的工程计算及数据处理软件,在众多领域得到广泛应用。它是一种面向对象的,交互式程序设计语言,其结构完整有优良的可移植性。它在矩阵运算,数字信号处理方面有强大的功能。另外,MATLAB提供了方便的绘图功能,便于用户直观地输出处理结果。 本课程实验要求学生运用MATLAB编程完成一些数字信号处理的基本功能,加深对教学内容的理解。 实验目的: 1.熟悉掌握MATLAB的基本操作 2.通过使用来熟悉常用离散信号; 实验内容: 1.1 G ENERATION OF SEQUENCES Project 1.1 Unit sample and unit step sequences A copy of Program P1_1 is given below. % Program P1_1 % Generation of a Unit Sample Sequence clf; % Generate a vector from -10 to 20 n = -10:20; % Generate the unit sample sequence u = [zeros(1,10) 1 zeros(1,20)]; % Plot the unit sample sequence stem(n,u); xlabel('Time index n');ylabel('Amplitude'); title('Unit Sample Sequence'); axis([-10 20 0 1.2]); < Insert program code here. Copy from m-file(s) and paste. > Answers: Q1.1The unit sample sequence u[n] generated by running Program P1_1 is shown below: < Insert MATLAB figure(s) here. Copy from figure window(s) and paste. > Q1.2The purpose of clf command is–清除当前窗口图形 The purpose of axis command is–设定坐标轴最大值最小值和

相关文档
相关文档 最新文档