《信号与系统及MATLAB实现》实验指导书
前言 长期以来，《信号与系统》课程一直采用单一理论教学方式，同学们依靠做习题来巩固和理解教学内容，虽然手工演算训练了计算能力和思维方法，但是由于本课程数学公式推导较多，概念抽象，常需画各种波形，作题时难免花费很多时间，现在，我们给同学们介绍一种国际上公认的优秀科技应用软件MA TLAB，借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。 MA TLAB的功能非常强大，我们此处仅用到它的一部分，在后续课程中我们还会用到它，在未来地科学研究和工程设计中有可能继续用它，所以有兴趣的同学，可以对MA TLAB 再多了解一些。 MA TLAB究竟有那些特点呢？ 1．高效的数值计算和符号计算功能，使我们从繁杂的数学运算分析中解脱出来； 2．完备的图形处理功能，实现计算结果和编程的可视化； 3．友好的用户界面及接近数学表达式的自然化语言，易于学习和掌握； 4．功能丰富的应用工具箱，为我们提供了大量方便实用的处理工具； MA TLAB的这些特点，深受大家欢迎，由于个人电脑地普及，目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景，我们编写了这本《信号与系统及MA TLAB实现》指导书，内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习，同学们在学习《信号与系统》的同时，掌握MA TLAB的基本应用，学会应用MA TLAB的数值计算和符号计算功能，摆脱烦琐的数学运算，从而更注重于信号与系统的基本分析方法和应用的理解与思考，将课程的重点、难点及部分习题用MA TLAB进行形象、直观的可视化计算机模拟与仿真实现，加深对信号与系统的基本原理、方法及应用的理解，为学习后续课程打好基础。另外同学们在进行实验时，最好事先预习一些MA TLAB的有关知识，以便更好地完成实验，同时实验中也可利用MA TLAB的help命令了解具体语句以及指令的使用方法。
学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟
20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求
实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。
实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1;
数字信号处理作业提交日期：2016年7月15日
实验一 维纳滤波器的设计 第一部分 设计一维纳滤波器。 (1)产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪(信噪比分别为20,10,6dB dB dB )，得到观测数据123(),(),()x n x n x n 。 (2)估计()i x n ，1,2,3i =的AR 模型参数。假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。 1 实验原理 滤波技术是信号分析、处理技术的重要分支，无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传递是至关重要的。信号分析检测与处理的一个十分重要的内容就是从噪声中提取信号，实现这种功能的有效手段之一是设计一种具有最佳线性过滤特性的滤波器，当伴有噪声的信号通过这种滤波器的时候，它可以将信号尽可能精确地重现或对信号做出尽可能精确的估计，而对所伴随噪声进行最大限度地抑制。维纳滤波器就是这种滤波器的典型代表之一。 维纳（Wiener ）是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。 设一线性系统的单位样本响应为()h n ，当输入以随机信号()x n ，且 ()() () x n s n v n =+，其中()s n 表示原始信号，即期望信号。()v n 表示噪声，则输出()y n 为()=()()m y n h m x n m -∑，我们希望信号()x n 经过线性系统()h n 后得到的()y n 尽可能接近 于()s n ，因此称()y n 为估计值，用?()s n 表示。 则维纳滤波器的输入-输出关系可用下面表示。 设误差信号为()e n ,则?()()()e n s n s n =-,显然)(n e 可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即 它的平方的统计期望最小：222?[|()|][|()()|][|()()|]E e n E s n s n E s n y n =-=-=min 。而要使均方误差最小，则需要满足2[|()|]j E e n h ?=0. 进一步导出维纳-霍夫方程为：()()()()*(),0,1,2...xs xx xx i R m h i R m i R m h m m =-==∑ 写成矩阵形式为：xs xx R R h =，可知：1xs xx h R R -=。表明已知期望信号与观测数据的互相关函数以及观测信号的自相关函数时，可以通过矩阵求逆运算，得到维纳滤波器的
实验报告 实验二连续时间系统的时域分析 一、实验目的： 1、掌握用Matlab进行卷积运算的数值方法和解析方法,加深对卷积积分的理解。 2、学习利用Matlab实现LTI系统的冲激响应、阶跃响应和零状态响应。 二、实验内容及步骤 实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。
1、 编写程序Q2_1，完成)(1t f 与)(2t f 两函数的卷积运算。 2、 编写程序Q2_2，完成)(1t f 与)(2t f 两函数的卷积运算。 3、编写程序Q2_3。利用程序Q2_1，验证卷积的相关性质。 (a) 验证性质：)()(*)(t x t t x =δ (b) 验证性质： )()(*)(00t t x t t t x -=-δ 4、编写程序Q2_4。某线性时不变系统的方程为 )(8)(2)(6)(5)(t f t f t y t y t y +'=+'+''， (a)系统的冲激响应和阶跃响应。 (b)输入()()t f t e u t -=，求系统的零状态响应)(t y zs 。 三. 实验结果 一： dt=0.01 t1=0:dt:2 f1=0.5*t1 t2=0:dt:2 f2=0.5*t2 f=dt*conv(f1,f2) t=0:0.01:4 plot(t,f);axis([-1 5 0 0.8])
二： dt=0.01 t=-3:dt:3 t1=-6:dt:6 ft1=2*rectpuls(t,2) ft2=rectpuls(t,4) y=dt*conv(ft1,ft2) plot(t1,y) axis([-4 4 0 5]) 以上两题出现错误点：（1）最开始模仿例1的写法用function [f,k]=sconv，总提示出现 错误 （2）t0+t2 ≤ t ≤ t1+t3 不大能理解的运用个特点，在编写的时候总是被忽略。导致t和t1设置的长度总出错。 三： （a） dt=0.01 t=0:dt:2 t0=0 t1=0:dt:2t2=0:dt: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所示）转换为级联型和并联型。
信号与系统实验报告 班级: 姓名: 信息与通信工程学院
实验一 系统的卷积响应 实验性质：提高性 实验级别：必做 开课单位：信息与通信工程学院 学 时：2 一、实验目的：深刻理解卷积运算，利用离散卷积实现连续卷积运算；深刻理解信号与系统的关系，学习MATLAB 语言实现信号通过系统的仿真方法。 二、实验设备： 计算机，MATLAB 软件 三、实验原理： 1、 离散卷积和： 调用函数：conv （） ∑∞ -∞ =-= =i i k f i f f f conv S )()(1)2,1(为离散卷积和， 其中，f1(k), f2 (k) 为离散序列，K=…-2, -1, 0 , 1, 2, …。但是，conv 函数只给出纵轴的序列值的大小，而不能给出卷积的X 轴序号。为得到该值，进行以下分析： 对任意输入：设)(1k f 非零区间n1~n2，长度L1=n2-n1+1；)(2k f 非零区间m1~m2，长度L2=m2-m1+1。则：)(*)()(21k f k f k s =非零区间从n1+m1开始，长度为L=L1+L2-1，所以S （K ）的非零区间为：n1+m1~ n1+m1+L-1。 2、 连续卷积和离散卷积的关系： 计算机本身不能直接处理连续信号，只能由离散信号进行近似： 设一系统（LTI ）输入为)(t P ?，输出为)(t h ?，如图所示。 )t )()(t h t P ??→
)()(lim )(lim )(0 t h t h t P t =→=?→??→?δ 若输入为f(t): ??-?= ≈∑∞ -∞ =? ?)()()()(k t P k f t f t f k 得输出： ??-?= ∑∞ -∞ =? ?)()()(k t h k f t y k 当0→?时：?∑∞ ∞-∞ -∞ =? →??→?-=??-?==ττδτd t f k t P k f t f t f k )()()()(lim )(lim )(0 ?∑∞ ∞ -∞ -∞ =? →??→?-= ??-?==τττd t h f k t h k f t y t y k )()()()(lim )(lim )(0 所以： ? ?-?=-==∑?→?)()(lim )()()(*)()(21 2121k t f k f d t f f t f t f t s τ ττ 如果只求离散点上的f 值)(n f ? ] )[()()()()(2121 ∑∑∞ -∞ =∞ -∞=?-??=? ?-??= ?k k k n f k f k n f k f n f 所以，可以用离散卷积和CONV （）求连续卷积，只需?足够小以及在卷积和的基础上乘以?。 3、 连续卷积坐标的确定： 设)(1t f 非零值坐标范围：t1~t2，间隔P )(2t f 非零值坐标范围：tt1~tt2，间隔P )(*)()(21t f t f t s =非零值坐标：t1+tt1~t2+tt2+1 根据给定的两个连续时间信号x(t) = t[u(t)-u(t-1)]和h(t) = u(t)-u(t-1)，编写程序，完成这
实验一MATLAB语言的基本使用方法 实验类别：基础性实验 实验目的： （1）了解MATLAB程序设计语言的基本方法，熟悉MATLAB软件运行环境。 （2）掌握创建、保存、打开m文件的方法，掌握设置文件路径的方法。 （3）掌握变量、函数等有关概念，具备初步的将一般数学问题转化为对应计算机模型并进行处理的能力。 （4）掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理。 实验内容和步骤： 1、打开MATLAB，熟悉MATLAB环境。 2、在命令窗口中分别产生3*3全零矩阵，单位矩阵，全1矩阵。 3、学习m文件的建立、保存、打开、运行方法。 4、设有一模拟信号f(t)=1.5sin60πt,取?t=0.001,n=0,1,2,…,N-1进行抽样，得到 序列f(n),编写一个m文件sy1_1.m,分别用stem,plot,subplot等命令绘制32 点序列f(n)(N=32)的图形，给图形加入标注，图注，图例。 5、学习如何利用MATLAB帮助信息。 实验结果及分析： 1）全零矩阵 >> A=zeros(3,3) A = 0 0 0 0 0 0 0 0 0 2）单位矩阵 >> B=eye(3) B = 1 0 0 0 1 0 0 0 1 3）全1矩阵 >> C=ones(3) C = 1 1 1 1 1 1 1 1 1 4)sy1_1.m N=32; n=0:N-1; dt=0.001; t=n*dt; y=1.5*sin(60*pi*t); subplot(2,1,1), plot(t,y); xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('正弦函数'); title('二维图形'); subplot(2,1,2), stem(t,y) xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('序列函数'); title('条状图形'); 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 二维图形 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 条状图形
信号与系统实验实验四：周期信号的傅里叶级数 小组成员： 黄涛13084220 胡焰焰13084219 洪燕东13084217
一、实验目的 1、分析典型的矩形脉冲信号，了解矩形脉冲信号谐波分量的构成。 2、观察矩形脉冲信号通过多个数字滤波器后，分解出各谐波分量的情况。 3、掌握用傅里叶级数进行谐波分析的方法。 4、观察矩形脉冲信号分解出的各谐波分量可以通过叠加合成出原矩形脉冲信号。 二、预习内容 1、周期信号的傅里叶级数分解及其物理意义。 2、典型信号傅里叶级数计算方法。 三、实验原理 1. 信号的时间特性与频率特性 信号可以表示为随时间变化的物理量，比如电压)(t u 和电流)(t i 等，其特性主要表现为随时间的变化，波形幅值的大小、持续时间的长短、变化速率的快慢、波动的速度及重复周期的大小等变化，信号的这些特性称为时间特性。 信号还可以分解为一个直流分量和许多不同频率的正弦分量之和。主要表现在各频率正弦分量所占比重的大小不同；主要频率分量所占的频率范围也不同，信号的这些特性称为信号的频率特性。 无论是信号的时间特性还是频率特性都包含了信号的全部信息量。 2. 信号的频谱 信号的时间特性和频率特性是对信号的两种不同的描述方式。根据傅里叶级数原理，任意一个时域的周期信号)t (f ，只要满足狄利克莱(Dirichlet)条件，就可以将其展开成三角形式或指数形式的傅里叶级数。例如，对于一个周期为T 的时域周期信号)t (f ，可以用三角形式的傅里叶级数求出它的各次分量，在区间),(11T t t +内表示为 ()∑∞ =Ω+Ω+=10sin cos )(n n n t n b t n a a t f 即将信号分解成直流分量及许多余弦分量和正弦分量，研究其频谱分布情况。 A 0t A n 0A 0t （a）（b） Ω（c）ωΩ 5Ω3Ω Ω3Ω5 3. 信号的时间特性与频率特性关系 信号的时域特性与频域特性之间有着密切的内在联系，这种联系可以用图4-1来形象地表示。其中图4-1(a)是信号在幅度--时间--频率三维坐标系统中的图形；图4-1(b)是信号在幅度--时间坐标系统中的图形即波形图；把周期信号分解得到的各次谐波分量按频率的高低排列，就可以得到频谱图。反映各频率分量幅度的频谱称为振幅频谱。图4-1(c)是信号在幅度--频率坐标系统中的图形即振幅频谱图。反映各分量相位的频谱称为相位频谱。 4. 信号频谱的测量 在本实验中只研究信号振幅频谱。周期信号的振幅频谱有三个性质：离散性、谐波性、收敛
实验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的围，频谱出现了镜像对称的部分。
信号与系统实验（六） 班级11083415 章仕波（11081522） 刘贺洋（11081515） 实验内容 1离散时间傅里叶变换 (1)下面参考程序是如下序列在范围44πωπ-≤≤的离散时间傅里叶变换 ()210.6j j j e F e e ω ω ω --+=- %计算离散时间傅里叶变换的频率样本 clear all; w=-4*pi;8*pi/511;4*pi; num=[2 1]; den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1) plot(w/pi,real(h)); grid; title(‘实部’) xlabel(‘omega/\pi ’); ylabel(‘振幅’); subplot(2,1,2) plot(w/pi, imag(h)); grid; title(‘虚部’) xlabel(‘omega/\pi ’); ylabel(‘振幅’); figure; subplot(2,1,1) plot(w/pi, abs(h)); grid; title(‘幅度谱’) xlabel(‘omega/\pi ’); ylabel(‘振幅’); subplot(2,1,2) plot(w/pi, angle (h)); grid; title(‘相位谱’) xlabel(‘omega/\pi ’); ylabel(‘以弧度为单位的相位’); 修改程序，在范围0ωπ≤≤内计算如下有限长序列的离散时间傅里叶变换 h[n]=[1 2 3 4 5 6 7 8 9] (2)利用(1)的程序，通过比较结果的幅度谱和相位谱，验证离散时间傅里叶变换的时移
武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02
一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分
实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列');
实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列');
正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列');
随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列');
实验一：信号的时域分析 一、实验目的 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)
实验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）’）；
实验二 常用信号分类与观察 一、实验目的 1、观察常用信号的波形特点及产生方法。 2、学会使用示波器对常用波形参数的测量。 二、实验内容 1、信号的种类相当的多，这里列出了几种典型的信号，便于观察。 2、这些信号可以应用到后面的“基本运算单元”和“无失真传输系统分析”中。 三、实验仪器 1、信号与系统实验箱一台（主板）。 2、20MHz 双踪示波器一台。 四、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、抽样信号、钟形信号、脉冲信号等。 1、正弦信号：其表达式为)sin()(θω+=t K t f ，其信号的参数：振幅K 、角频率ω、与初始相位θ。其波形如下图所示： 图 1-5-1 正弦信号 2、指数信号：指数信号可表示为at Ke t f =)(。对于不同的a 取值，其波形表现为不同的形式，如下图所示：
图 1-5-2 指数信号 3、指数衰减正弦信号：其表达式为 ?? ???><=-)0()sin()0(0)(t t Ke t t f at ω 其波形如下图： 图 1-5-3 指数衰减正弦信号 4、抽样信号：其表达式为： sin ()t Sa t t = 。)(t Sa 是一个偶函数，t = ±π,±2π,…,±n π时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：
实验一 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;
实验1 信号的时域描述与运算 一、实验目的 1. 掌握信号的MATLAB表示及其可视化方法。 2. 掌握信号基本时域运算的MA TLAB实现方法。 3. 利用MA TLAB分析常用信号，加深对信号时域特性的理解。 二、实验原理与方法 1. 连续时间信号的MATLAB表示 连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。 从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： >> t=0:0.01:10; >> x=sin(t); 利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。 如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： >> x=sin(t); >> ezplot(X); 利用ezplot(x)命令可以绘制上述信号的时域波形 Time(seconds) 图1 利用向量表示连续时间信号
t 图 2 利用符号对象表示连续时间信号 sin(t) 2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。 1）相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。 2）微分和积分 对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量[N t t t ,,,21?]和采样值向量[N x x x ,,,21?]表示的连续时间信号，其微分可以通过下式求得 1,,2,1,|)('1-?=?-≈ +=N k t x x t x k k t t k 其中t ?表示采样间隔。MA TLAB 中用diff 函数来计算差分 k k x x -+1。 连续时间信号的定积分可以由MATLAB 的qud 函数实现，调用格式为 quad ('function_name',a,b) 其中，function_name 为被积函数名，a 、b 为积分区间。
一.实验目的 1.复习采样定理 2.掌握应用matlab 函数设计模拟滤波器的方法 3.掌握系统性能分析的方法 4.结合实际综合应用信号与系统的基础理论 二.实验原理 在数字语音系统中，需首先对语音信号（模拟信号）采样，语音信号频率范围[-fh ，fh]，信号中一般含有干扰噪声，其频带宽度远大于fh 。本次实验以电话系统中的语音信号采样系统为对象，设计语音信号采样前滤波器。数字电话系统结构框图如图8.1，电话系统中一般要保证4kHz 的音频带宽，即取fh ＝4kHz ，但送话器发出的信号的带宽比fh 大很多。因此在A/D 转换之前需对其进行模拟预滤波，以防止采样后发生频谱混叠失真。为使信号采集数量尽量少，设模数转换器的采样频率为8kHz 。 图8.1 数字电话系统结构框图 滤波器的定义 在信号处理时，通常都会遇到有用信号中混入(叠加)噪声的问题，消除或减弱噪声对信号的干扰，是信号处理中的一种最基本且重要的技术。根据有用信号与噪声不同的特性，抑制不需要的噪声或干扰， 提取出有用信号的过程称为滤波，实现滤波功能的装置称为滤波器。 在A/D 变换前，常常需要设置一个模拟滤波器进行预滤波以限 制信号带宽，去掉高于1/2抽样频率以上的高频分量，防止频谱 混叠现象的发生，称为抗混叠滤波器或预抽样滤波器 模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。 典型的模拟滤波器
巴特沃斯 Butterworth 滤波器 幅频特性单调下降 切比雪夫 Chebyshev 滤波器 幅频特性在通带或者在阻带有波动 贝塞尔 Bessel 滤波器 通带内有较好的线性相位持性 椭圆 Ellipse 滤波器 以这些数学函数命名的滤波器是低通滤波器的原型 模拟滤波器按幅度特性可分成低通、高通、带通和带阻滤波器，它们的理想幅度特性如图所示。 模拟低通滤波器的设计指标有αp, Ωp,αs 和Ωs 。 Ωp ；通带截止频率 Ωs ：阻带截止频率 αp ：通带中最大衰减系数 αs ；阻带最小衰减系数 αp 和αs 一般用dB 数表示。对于单调下降的幅度特性，可表示成： 222 2 (0) (0) 10lg 10lg () () a a p s a s a p H j H j H j H j αα==ΩΩ 三.实验内容
信息科学与技术学院本科三年级 数字信号处理实验报告 2011 年12 月21日
实验一 序列的傅立叶变换 实验目的 进一步加深理解DFS,DFT 算法的原理;研究补零问题;快速傅立叶变换 （FFT ）的应用。 实验步骤 1. 复习DFS 和DFT 的定义，性质和应用； 2. 熟悉MATLAB 语言的命令窗口、编程窗口和图形窗口的使用；利用提供的 程序例子编写实验用程序；按实验内容上机实验，并进行实验结果分析；写出完整的实验报告，并将程序附在后面。 实验内容 1. 周期方波序列的频谱试画出下面四种情况下的的幅度频谱,并分析补零后，对信号频谱的影响。 实验结果： 60 ,7)4(;60,5)3(; 40,5)2(;20,5)1()] (~[)(~,2,1,01 )1(,01,1)(~=========±±=???-+≤≤+-+≤≤=N L N L N L N L n x DFS k X m N m n L m N L m N n m N n x ) 52.0cos()48.0cos()(n n n x ππ+=
2. 有限长序列x(n)的DFT （1） 取x(n)(n=0:10)时，画出x(n)的频谱X(k) 的幅度； （2） 将（1）中的x(n)以补零的方式，使x(n)加长到(n:0~100)时，画出 x(n)的频谱X(k) 的幅度； （3） 取x(n)(n:0~100)时，画出x(n)的频谱X(k) 的幅度。利用FFT 进行谱分析 已知：模拟信号 以t=0.01n(n=0:N-1)进行采样，求N 点DFT 的幅值谱。 请分别画出N=45; N=50;N=55;N=60时的幅值曲线。 实验结果： ) 8cos(5)4sin(2)(t t t x ππ+=
信号与系统实验 实验1 阶跃响应与冲激响应 一、实验目的 1、观察和测量RLC串联电路的阶跃响应与冲激响应的波形和有关参数，并 研究其电路元件参数变化对响应状态的影响； 2、掌握有关信号时域的测量方法。 二、实验原理说明 实验如图1－1所示RLC串联电路的阶跃响应与冲激响应的电路连接图，图1
用周期方波通过微分电路后得到的尖顶脉冲代替冲激信号。 三、实验内容 1、阶跃响应波形观察与参数测量 设激励信号为方波，其幅度为1.5V 峰峰值，频率为500Hz 。 实验电路连接图如图1-1（a ）所示。 ① 连接如图1-1所示 ② 调整激励源信号为方波，调节频率旋钮，使f=500Hz ，调节幅度旋钮， 使信号幅度为1.5V 。（注意：实验中，在调整信号源的输出信号的参数时，需连接上负载后调节） ③ 示波器CH1接于TP909，调节滑动变阻器，使电路分别工作于欠阻尼、 临界和过阻尼三种状态，并将实验数据填入表格1-1中。 ④ TP908为输入信号波形的测量点，可把示波器的CH ·接于TP908上，便 于波形比较。 表1-1 注:描绘波形要使三状态的X 轴坐标(扫描时间)一致。 2、冲激响应的波形观察 冲激信号是由阶跃信号经过微分电路而得到。 实验电路如图1—1（b ）所示。 参数测量 波形观察 欠阻尼状态 临界状态 过阻尼状态 状态 参数测量 R< Tr= Ts= δ= R= Tr= R>
①将信号输入接于P905。（频率与幅度不变）； ②将示波器的CH1接于TP906，观察经微分后响应波形（等效为冲激激 励信号）； ③连接如图1-1（b）所示 ④将示波器的CH2接于TP909，调整滑动变阻器，使电路分别工作于欠 阻尼、临界和过阻尼三种状态 ④观察TP909端三种状态波形，并填于表1-2中。 表1-2 表中的激励波形为在测量点TP906观察到的波形（冲激激励信号）。 四、实验报告要求 1、描绘同样时间轴阶跃响应与冲激响应的输入、输出电压波形时， 要标明信号幅度A、周期T、方波脉宽T1以及微分电路的τ值。 2、分析实验结果，说明电路参数变化对状态的影响。 五、实验设备 双踪示波器 1 台 信号系统实验箱 1台 上升时间t r :y(t)从0.1到第一次达到0.9所需时间。 峰值时间t p :y(t)从0上升y max 所需的时间。 调节时间t s :y(t)的振荡包络线进入到稳态值的% 5 误差范围所需的时间。 激励波形 响应波形 欠阻尼状态临界状态过阻尼状态