信号与系统实验DOC

合集下载

信号与系统实验

实验一信号与系统认知一、实验目的1、了解实验室的规章制度、强化安全教育、说明考核方法。

2、学习示波器、实验箱的使用、操作知识；3、学习常用连续周期信号的波形以及常用系统的作用。

二、实验仪器1、信号与系统实验箱（本次实验使用其自带的简易信号源，以及实验箱上的“信号通过系统”部分。

）2、示波器三、实验原理1、滤波器滤波器是一种常用的系统，它的作用为阻止某些频率信号通过，或只允许某些频率的信号通过。

滤波器主要有四种：这是四种滤波器的理想状态，实际上的滤波器只能接近这些效果，因此通常的滤波器有一些常用的参数：如带宽、矩形系数等。

通带范围：与滤波器最低衰减处比，衰减在3dB以下的频率范围。

2、线性系统线性系统是现实中广泛应用的一种系统，线性也是之后课程中默认为系统都具有的一种系统性质。

系统的线性表现在可加性与齐次性上。

齐次性：输入信号增加为原来的a倍时，输出信号也增加到原来的a倍。

四、预习要求1、复习安全操作的知识。

2、学习或复习示波器的使用方法。

3、复习典型周期信号的波形及其性质。

4、复习线性系统、滤波器的性质。

5、撰写预习报告。

五、实验内容及步骤1、讲授实验室的规章制度、强化安全教育、说明考核方法2、通过示波器，读出实验箱自带信号源各种信号的频率范围（1）测试信号源1的各种信号参数，并填入表1-1。

（2）测试信号源2的各种信号参数，并填入表1-2。

3、测量滤波器根据相应测量方法，用双踪示波器测出实验箱自带的滤波器在各频率点的输入输出幅度（先把双踪示波器两个接口都接到所测系统的输入端，调节到都可以读出输入幅度值，并把两侧幅度档位调为一致，记录下这个幅度值；之后，将示波器的一侧改接入所测系统的输出端，再调节用于输入的信号源，将信号频率其调至表1-3中标示的值，并使输入信号幅度保持原幅度值不变。

观察输出波形幅度的变化，并与原来的幅度作比较，记录变化后的幅度值。

），并将相应数据计入表1-3中。

4、测量线性系统（1）齐次性的验证自选一个输入信号，观察输出信号的波形并记录输入输出信号的参数，将输入信号的幅度增强为原信号的一定倍数后，再对输入输出输出参数进行记录，对比变化前后的输出。

《信号与系统》课程实验报告

《信号与系统》课程实验报告《信号与系统》课程实验报告一图1-1 向量表示法仿真图形2．符号运算表示法若一个连续时间信号可用一个符号表达式来表示，则可用ezplot命令来画出该信号的时域波形。

上例可用下面的命令来实现（在命令窗口中输入，每行结束按回车键）。

t=-10:0.5:10;f=sym('sin((pi/4)*t)');ezplot(f,[-16,16]);仿真图形如下：图1-2 符号运算表示法仿真图形三、实验内容利用MATLAB实现信号的时域表示。

三、实验步骤该仿真提供了7种典型连续时间信号。

用鼠标点击图0-3目录界面中的“仿真一”按钮，进入图1-3。

图1-3 “信号的时域表示”仿真界面图1-3所示的是“信号的时域表示”仿真界面。

界面的主体分为两部分：1) 两个轴组成的坐标平面（横轴是时间，纵轴是信号值）；2) 界面右侧的控制框。

控制框里主要有波形选择按钮和“返回目录”按钮，点击各波形选择按钮可选择波形，点击“返回目录”按钮可直接回到目录界面。

图1-4 峰值为8V，频率为0.5Hz，相位为180°的正弦信号图1-4所示的是正弦波的参数设置及显示界面。

在这个界面内提供了三个滑动条，改变滑块的位置，滑块上方实时显示滑块位置代表的数值，对应正弦波的三个参数：幅度、频率、相位；坐标平面内实时地显示随参数变化后的波形。

在七种信号中，除抽样函数信号外，对其它六种波形均提供了参数设置。

矩形波信号、指数函数信号、斜坡信号、阶跃信号、锯齿波信号和抽样函数信号的波形分别如图1-5～图1-10所示。

图1-5 峰值为8V，频率为1Hz，占空比为50%的矩形波信号图1-6 衰减指数为2的指数函数信号图1-7 斜率=1的斜坡信号图1-8 幅度为5V，滞后时间为5秒的阶跃信号图1-9 峰值为8V，频率为0.5Hz的锯齿波信号图1-10 抽样函数信号仿真途中，通过对滑动块的控制修改信号的幅度、频率、相位，观察波形的变化。

《信号与系统》实验三

实验记录及个人小结（包括：实验源程序、注释、结果分析与讨论等）
三：
源程序：
（1）：τ/T=1/4时的周期矩形脉冲的幅度谱和相位谱：
n=-20:20;
F=zeros(size(n));
forii=-20:20
F(ii+21)= sin(ii*pi/4)/(ii*pi+eps);
end
F(21)=1/4;
实验
内容
1.求图1所示周期信号（，）的傅里叶级数，用Matlab做出其前3、9、21、45项谐波的合成波形与原信号作比较，并做出其单边幅度谱和相位谱。
图1 周期为2的三角脉冲信号
2. 求图2所示的单个三角脉冲（）的傅里叶变换，并做出其幅度谱和相位谱。
图2 单个三角脉冲
3. 求不同占空比下周期矩形脉冲的幅度谱和相位谱，例如、。
y=1/4;
forn=1:m
y=y+4/(n*n*pi*pi)*(1-cos(n*pi/2)).*cos(n*pi.*t);
end
源代码：
t=-6:0.01:6;
d=-6:2:6;
fxx=pulstran(t,d,'tripuls');
f1=fourierseries(3,t);
f2=fourierseries(9,t);
n=1:10;
a=zeros(size(n));
fori=1:10
a(i)=angle(4/(i*i*pi*pi)*(1-cos(i*pi/2)))
end
n=0:pi:9*pi
stem(n,a,'fill','linewidth',2);
axis([0,9*pi,-0.2,0.2])

信号与系统实验报告

电气学科大类2012 级《信号与控制综合实验》课程实验报告(基本实验一:信号与系统基本实验)姓名丁玮学号U201216149 专业班号水电1204 同组者1 余冬晴学号U201216150 专业班号水电1204 同组者2 学号专业班号指导教师日期实验成绩评阅人实验评分表基本实验实验编号名称/内容实验分值评分实验一常用信号的观察实验二零输入响应、零状态相应及完全响应实验五无源滤波器与有源滤波器实验六LPF、HPF、BPF、BEF间的变换实验七信号的采样与恢复实验八调制与解调设计性实验实验名称/内容实验分值评分创新性实验实验名称/内容实验分值评分教师评价意见总分目录1.实验一常用信号的观察 (1)2.实验二零输入响应、零状态响应及完全响应 (4)3.实验五无源滤波器与有源滤波器 (7)4.实验六 LPF、HPF、BPF、BEF间的转换 (14)5.实验七信号的采样与恢复 (19)6.实验八调制与解调 (29)7.实验心得与自我评价 (33)8.参考文献 (34)实验一常用信号的观察一．任务与目标1.了解常见信号的波形和特点；2.了解常见信号有关参数的测量，学会观察常见信号组合函数的波形；3.学会使用函数发生器和示波器，了解所用仪器原理与所观察信号的关系；4.掌握基本的误差观察与分析方法。

二．总体方案设计1.实验原理描述信号的方法有许多种，可以用数学表达式（时间的函数），也可以使用函数图形（信号的波形）。

信号可以分为周期信号和非周期信号两种。

普通示波器可以观察周期信号，具有暂态拍摄功能的示波器可以观察到非周期信号的波形。

目前，常用的数字示波器可以方便地观察周期信号及非周期信号的波形。

2.总体设计⑴观察常用的正弦波、方波、三角波、锯齿波等信号及一些组合函数的波形，如y=sin(nx)+cos(mx)。

⑵用示波器测量信号，读取信号的幅值与频率。

三．方案实现与具体设计1.用函数发生器产生正弦波，并且设定波形的峰值及频率，用示波器观察并记录波形，测量和读取信号的幅值与频率；2.用函数发生器产生方波，并且设定波形的峰值及频率，用示波器观察并记录波形，测量和读取信号的幅值与频率；3.用函数发生器产生三角波，并且设定波形的峰值及频率，用示波器观察并记录波形，测量和读取信号的幅值与频率；4.用函数发生器产生锯齿波，并且设定波形的峰值及频率，用示波器观察并记录波形，测量和读取信号的幅值与频率；5.用函数发生器产生两个不同频率的正弦波，分别设定波形的峰值及频率，用示波器叠加波形，并观察组合函数的波形。

信号与系统实验报告

信号与系统实验报告一、实验目的(1) 理解周期信号的傅里叶分解，掌握傅里叶系数的计算方法；(2)深刻理解和掌握非周期信号的傅里叶变换及其计算方法；(3) 熟悉傅里叶变换的性质，并能应用其性质实现信号的幅度调制；(4) 理解连续时间系统的频域分析原理和方法，掌握连续系统的频率响应求解方法，并画出相应的幅频、相频响应曲线。

二、实验原理、原理图及电路图(1) 周期信号的傅里叶分解设有连续时间周期信号()f t ，它的周期为T ，角频率22fT，且满足狄里赫利条件，则该周期信号可以展开成傅里叶级数，即可表示为一系列不同频率的正弦或复指数信号之和。

傅里叶级数有三角形式和指数形式两种。

1）三角形式的傅里叶级数：01212011()cos()cos(2)sin()sin(2)2cos()sin()2n n n n a f t a t a t b t b t a a n t b n t 式中系数n a ，n b 称为傅里叶系数，可由下式求得：222222()cos(),()sin()T T T T nna f t n t dtb f t n t dtTT2）指数形式的傅里叶级数：()jn tn nf t F e式中系数n F 称为傅里叶复系数，可由下式求得：221()T jn tT nF f t edtT周期信号的傅里叶分解用Matlab进行计算时，本质上是对信号进行数值积分运算。

Matlab中进行数值积分运算的函数有quad函数和int函数。

其中int函数主要用于符号运算，而quad函数（包括quad8，quadl）可以直接对信号进行积分运算。

因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算，也可以采用符号运算方法。

quadl函数(quad系)的调用形式为：y＝quadl(‘func’,a,b)或y＝quadl(@myfun,a,b)。

其中func是一个字符串，表示被积函数的.m文件名（函数名）；a、b分别表示定积分的下限和上限。

《信号与系统实验》信号的采样与恢复(抽样定理)实验

《信号与系统实验》信号的采样与恢复（抽样定理）实验一、实验目的1、了解电信号的采样方法与过程以及信号恢复的方法。

2、验证抽样定理。

二、实验设备1、信号与系统实验箱2、双踪示波器三、原理说明1、离散时间信号可以从离散信号源获得，也可以从连续时间信号抽样而得。

抽样信号f s(t)可以看成连续f(t)和一组开关函数s (t)的乘积。

s (t)是一组周期性窄脉冲，见实验图5-1，T s(t)称为抽样周期，其倒数f s(t)= 1/T s称为抽样频率。

图5-1 矩形抽样脉冲对抽样信号进行傅立叶分析可知，抽样信号的频率包括了原连续信号以及无限个经过平移的信号频率。

平移的频率等于抽样频率f s(t)及其谐波频率2f s、3f s》》》》》》。

当抽样信号是周期性窄脉冲时，平移后的频率幅度(sinx)/x规律衰减。

抽样信号的频谱是原信号频谱周期的延拓，它占有的频带要比原信号频谱宽得多。

2、正如测得了足够的实验数据以后，我们可以在坐标纸上把一系列数据点连起来，得到一条光滑的曲线一样，抽样信号在一定条件下也可以恢复到原信号。

只要用一截止频率等于原信号频谱中最高频率f n的低通滤波器，滤除高频分量，经滤波后得到的信号包含了原信号频谱的全部内容，故在低通滤波器输出可以得到恢复后的原信号。

3、但原信号得以恢复的条件是f s 2，其中f s为抽样频率，为原信号占有的频带宽度。

而f min=2 为最低抽样频率又称“柰奎斯特抽样率”。

当f s<2 时，抽样信号的频谱会发生混迭，从发生混迭后的频谱中我们无法用低通滤波器获得原信号频谱的全部内容。

在实际使用中，仅包含有限频率的信号是及少的，因此即使f s=2 ，恢复后的信号失真还是难免的。

图5-2画出了当抽样频率f s>2 （不混叠时）f s<2 （混叠时）两种情况下冲激抽样信号的频谱。

t f(t)0F()t 0m ωm ω-(a)连续信号的频谱Ts t 0f s (t)F()t0m ωm ω-s ω-s ω（）(b)高抽样频率时的抽样信号及频谱不混叠图5-2 冲激抽样信号的频谱实验中f s >2 、f s =2 、f s <2 三种抽样频率对连续信号进行抽样，以验证抽样定理——要使信号采样后能不失真地还原，抽样频率f s 必须大于信号频率中最高频率的两倍。

信号与系统试验----信号卷积

一、实验目的1. 理解卷积的概念及物理意义；2. 通过实验的方法加深对卷积运算的图解方法及结果的理解。

二、实验设备1.信号与系统实验箱 1台2.双踪示波器1台三、实验原理卷积积分的物理意义是将信号分解为冲激信号之和，借助系统的冲激响应，求解系统对任意激励信号的零状态响应。

设系统的激励信号为)t (x ，冲激响应为)t (h ，则系统的零状态响应为)(*)()(t h t x t y =⎰∞∞--=ττd t h t x )()(。

对于任意两个信号)t (f 1和)t (f 2，两者做卷积运算定义为：⎰∞∞--=ττd t f t f t f )(2)(1)(=)t (f 1*)t (f 2=)t (f 2*)t (f 1。

1. 两个矩形脉冲信号的卷积过程两信号)t (x 与)t (h 都为矩形脉冲信号，如图9-1所示。

下面由图解的方法（图9-1）给出两个信号的卷积过程和结果，以便与实验结果进行比较。

0≤<∞-t210≤≤t 1≤≤t 41≤≤t ∞<≤t 2124τ（b)(a)(c)(d)(e)(f)(g)(h)(i)2卷积结果2. 矩形脉冲信号与锯齿波信号的卷积信号)t (f 1为矩形脉冲信号，)t (f 2为锯齿波信号，如图9-2所示。

根据卷积积分的运算方法得到)t (f 1和)t (f 2的卷积积分结果)t (f ，如图9-2(c)所示。

图9-2 矩形脉冲信号与锯齿脉冲信号的卷积积分的结果3. 本实验进行的卷积运算的实现方法在本实验装置中采用了DSP 数字信号处理芯片，因此在处理模拟信号的卷积积分运算时，是先通过A/D 转换器把模拟信号转换为数字信号，利用所编写的相应程序控制DSP 芯片实现数字信号的卷积运算，再把运算结果通过D/A 转换为模拟信号输出。

结果与模拟信号的直接运算结果是一致的。

数字信号处理系统逐步和完全取代模拟信号处理系统是科学技术发展的必然趋势。

图9-3为信号卷积的流程图。

信号与系统实验

实验一抽样定理与信号恢复一、实验目的1. 观察离散信号频谱，了解其频谱特点；2. 验证抽样定理并恢复原信号。

二、实验原理1. 离散信号不仅可从离散信号源获得，而且也可从连续信号抽样获得。

抽样信号 Fs （t ）=F （t ）·S （t ）。

其中F （t ）为连续信号（例如三角波），S （t ）是周期为Ts 的矩形窄脉冲。

Ts 又称抽样间隔，Fs=1Ts 称抽样频率，Fs （t ）为抽样信号波形。

F （t ）、S （t ）、Fs （t ）波形如图1-1。

t-4T S -T S 0T S 4T S8T S 12T S tt02/1τ1τ2/31τ2/1τ1τ2/31τ2/1τ-(a)(b)(c)图1-1 连续信号抽样过程将连续信号用周期性矩形脉冲抽样而得到抽样信号，可通过抽样器来实现，实验原理电路如图1-2所示。

2. 连续周期信号经周期矩形脉冲抽样后，抽样信号的频谱()∑∞∞--∙=m s s m m SaTsA j )(22s F ωωπδτωτω 它包含了原信号频谱以及重复周期为fs （f s =πω2s 、幅度按ST A τSa （2τωs m ）规律变化的原信号频谱，即抽样信号的频谱是原信号频谱的周期性延拓。

因此，抽样信号占有的频带比原信号频带宽得多。

以三角波被矩形脉冲抽样为例。

三角波的频谱 F （j ω）=∑∞-∞=-K k k sa E )2()2(12τπωδππ抽样信号的频谱Fs （j ω）=式中取三角波的有效带宽为31ω18f f s =作图，其抽样信号频谱如图1-3所示。

图1-2 信号抽样实验原理图)(2(212s m k s m k k Sa m Sa TS EA ωωωδπτωτπ--∙∙∑∞-∞=-∞=111112ττπω==f 或(b) 抽样信号频谙图1-3 抽样信号频谱图如果离散信号是由周期连续信号抽样而得，则其频谱的测量与周期连续信号方法相同，但应注意频谱的周期性延拓。

信号与系统实验报告

信号与系统实验报告目录1. 内容概要 (2)1.1 研究背景 (3)1.2 研究目的 (4)1.3 研究意义 (4)2. 实验原理 (5)2.1 信号与系统基本概念 (7)2.2 信号的分类与表示 (8)2.3 系统的分类与表示 (9)2.4 信号与系统的运算法则 (11)3. 实验内容及步骤 (12)3.1 实验一 (13)3.1.1 实验目的 (14)3.1.2 实验仪器和设备 (15)3.1.4 实验数据记录与分析 (16)3.2 实验二 (16)3.2.1 实验目的 (17)3.2.2 实验仪器和设备 (18)3.2.3 实验步骤 (19)3.2.4 实验数据记录与分析 (19)3.3 实验三 (20)3.3.1 实验目的 (21)3.3.2 实验仪器和设备 (22)3.3.3 实验步骤 (23)3.3.4 实验数据记录与分析 (24)3.4 实验四 (26)3.4.1 实验目的 (27)3.4.2 实验仪器和设备 (27)3.4.4 实验数据记录与分析 (29)4. 结果与讨论 (29)4.1 实验结果汇总 (31)4.2 结果分析与讨论 (32)4.3 结果与理论知识的对比与验证 (33)1. 内容概要本实验报告旨在总结和回顾在信号与系统课程中所进行的实验内容，通过实践操作加深对理论知识的理解和应用能力。

实验涵盖了信号分析、信号处理方法以及系统响应等多个方面。

实验一：信号的基本特性与运算。

学生掌握了信号的表示方法，包括连续时间信号和离散时间信号，以及信号的基本运算规则，如加法、减法、乘法和除法。

实验二：信号的时间域分析。

在本实验中，学生学习了信号的波形变换、信号的卷积以及信号的频谱分析等基本概念和方法，利用MATLAB工具进行了实际的信号处理。

实验三：系统的时域分析。

学生了解了线性时不变系统的动态响应特性，包括零状态响应、阶跃响应以及脉冲响应，并学会了利用MATLAB进行系统响应的计算和分析。

信号与系统实验报告

信号与系统实验报告一、信号的时域基本运算1．连续时间信号的时域基本运算两实验之一实验分析：输出信号值就等于两输入信号相加（乘）。

由于b=2，故平移量为2时，实际是右移1，符合平移性质。

两实验之二心得体会：时域中的基本运算具有连续性，当输入信号为连续时，输出信号也为连续。

平移，伸缩变化都会导致输出结果相对应的平移伸缩。

2.离散时间信号的时域基本运算两实验之一实验分析：输出信号的值是对应输入信号在每个n值所对应的运算值，当进行拉伸变化后，n值数量不会变，但范围会拉伸所输入的拉伸系数。

两实验之二心得体会：离散时间信号可以看做对连续时间信号的采样，而得到的输出信号值，也可以看成是连续信号所得之后的采样值。

二、连续信号卷积与系统的时域分析1.连续信号卷积积分两实验之一实验分析：当两相互卷积函数为冲激函数时，所卷积得到的也是一个冲激函数，且该函数的冲激t值为函数x，函数y冲激t值之和。

两实验之二心得体会：连续卷积函数每个t值所对应的卷积和可以看成其中一个在k值取得的函数与另外一个函数相乘得到的一个分量函数，并一直移动k值直至最后，最后累和出来的最终函数便是所得到的卷积函数。

3.RC电路时域积分两实验之一实验分析：全响应结果正好等于零状态响应与零输入响应之和。

两实验之二心得体会：具体学习了零状态，零输入，全响应过程的状态及变化，与之前所学的电路知识联系在一起了。

三、离散信号卷积与系统的时域分析1.离散信号卷积求和两实验之一实验分析：输出结果的n值是输入结果的k号与另一个n-k的累和两实验之二心得体会：直观地观察到卷积和的产生，可以看成连续卷积的采样形式，从这个方面去想，更能深入地理解卷积以及采样的知识。

2.离散差分方程求解两实验之一实验分析：其零状态响应序列为0 0 4 5 7.5，零输入响应序列为2 4 5 5.5 5.75，全状态响应序列为2 4 9 10.5 13.25，即全状态=零输入+零状态。

两实验之二心得体会：求差分方程时，可以根据全状态响应是由零输入输入以及零状态相加所得，分开来求，同时也加深了自己对差分方程的求解问题的理解。

信号与系统课程实验报告

合肥工业大学宣城校区《信号与系统》课程实验报告专业班级学生姓名《信号与系统》课程实验报告一实验名称一阶系统的阶跃响应姓名系院专业班级学号实验日期指导教师成绩一、实验目的1．熟悉一阶系统的无源和有源电路；2．研究一阶系统时间常数T的变化对系统性能的影响；3．研究一阶系统的零点对系统响应的影响。

二、实验原理1．无零点的一阶系统无零点一阶系统的有源和无源电路图如图2-1的(a)和(b)所示。

它们的传递函数均为：10.2s1G(s)＝+(a) 有源(b) 无源图2-1 无零点一阶系统有源、无源电路图2．有零点的一阶系统(|Z|<|P|)图2-2的(a)和(b)分别为有零点一阶系统的有源和无源电路图，它们的传递函数为：10.2s1)0.2(sG(s)++=，⎪⎪⎪⎪⎭⎫⎝⎛++=S611S161G(s）(a) 有源(b) 无源图2-2 有零点(|Z|<|P|)一阶系统有源、无源电路图3．有零点的一阶系统(|Z|>|P|)图2-3的(a)和(b)分别为有零点一阶系统的有源和无源电路图，它们的传递函数为：1s10.1sG(s)＝++(a) 有源(b) 无源图2-3 有零点(|Z|>|P|)一阶系统有源、无源电路图三、实验步骤1．打开THKSS-A/B/C/D/E型信号与系统实验箱，将实验模块SS02插入实验箱的固定孔中，利用该模块上的单元组成图2-1(a)(或(b))所示的一阶系统模拟电路。

2．实验线路检查无误后，打开实验箱右侧总电源开关。

3．将“阶跃信号发生器”的输出拨到“正输出”，按下“阶跃按键”按钮，调节电位器RP1，使之输出电压幅值为1V，并将“阶跃信号发生器”的“输出”端与电路的输入端“Ui”相连，电路的输出端“Uo”接到双踪示波器的输入端，然后用示波器观测系统的阶跃响应，并由曲线实测一阶系统的时间常数T。

4．再依次利用实验模块上相关的单元分别组成图2-2(a)(或(b))、2-3(a)(或(b))所示的一阶系统模拟电路，重复实验步骤3，观察并记录实验曲线。

信号与系统实验报告一-连续时间信号

实验一连续时间信号§1.2 连续时间复指数信号基本题1．对下面信号创建符号表达式()()t t t x ππ2c o s2sin )(= 这两个信号应分别创建，然后用symmul 组合起来。

对于T=4，8和16，利用ezplot 画出320≤≤t 内的信号。

什么是)(t x 的基波周期？x(t) =cos((pi*t)/2)*sin((pi*t)/2)=1/2sin(pi*t) (T=4)若令f1=1 /T1=1/2,很容易得到其基波分量：1/2sin(pi*t)同理可得：x(t)=cos((pi*t)/4)*sin((pi*t)/4)=1/2sin((pi*t)/2) (T=8)其基波分量为1/2sin((pi*t)/2),基频为f1=1/T1=1/4x(t)= cos((pi*t)/8)*sin((pi*t)/8)=1/2sin((pi*t)/4) (T=16)其基波分量为1/2sin((pi*t)/4)，基频为f1=1/T1=1/8 中等题2．对下面信号创建一个符号表达式()t e t x at π2cos )(-=对于81,41,21=a ，利用ezplot 确定d t ，d t 为)(t x 最后跨过0.1的时间，将d t 定义为该信号的消失的时间。

利用ezplot 对每一个a 值确定在该信号消失之前，有多少个完整的余弦周期出现，周期数目是否正比于品质因素a T Q 2)2(π=？1）当a=1/2时： x(t)= cos(2*pi*t)/exp(t/2)利用Tool菜单中的data cursor项目可大致确定d t=4.548在该信号消失之前，有个约4(4.5)完整的余弦周期出现，对应的品质因数为6.28。

2）当a=1/4时： x(t)= cos(2*pi*t)/exp(t/4)利用Tool菜单中的data cursor项目可大致确定d t=9.053在该信号消失之前，有个约9完整的余弦周期出现，对应的品质因数为12.57。

信号与系统实验五（docX页）

实验五低通滤波系统的频率特性分析实验报告一、实验名称低通滤波系统的频率特性分析二、实验目的（1）观察理想低通滤波器的单位冲击响应与频谱图；（2）观察RC 低通网络的单位冲击响应与频谱图。

三、实验原理RC 低通滤波电路如图其系统函数为()()()()12211tan 11j RC RC H RC H RC RC ωωωωωω-==∠-∠∂++式中()()2211RC H RC ωω=+称为幅频特性； ()()1tan RC H ωωω-∂=- 称为相频特性。

当0ω=，()()1,;H ωω=∂当11RC ωτ==时，()12H ω=，()45ω︒∂=-；当ω→+∞时，()0H ω→，()90ω︒∂→-。

电路的幅频特性表明，对于同样大小的输入信号，频率越高，输出信号衰减越大；频率越低，输出信号衰减越小或者可以认为无衰减。

也就是说，对该电路而言，低频信号比较容易通过，而高频信号则不容易通过，因此这个电路称为低通滤波器。

（1）理想低通的单位冲击响应为()0Sa t t - 函数，幅频特性在通带内为常数，阻带内为零。

在截止频率点存在阶越性跳变。

相频特性为通过原点斜率为 0t ω- 的直线。

（2）实际物理可实现的RC 低通网络通带阻带存在过渡时间，与RC 时间常数有关，通带阻带也不在完全是常数。

相频特性为通过原点的曲线（在原点附近近似为直线）。

四、实验步骤（1）打开MA TLAB 软件，建立一个M 文件。

（2）MA TLAB 所在目录的\work 子目录下建立一个名为heaviside 的M 文件，创建子程序函数。

（3）建立一个新的M 文件，编写主程序并保存。

（4）运行主程序，观察理想低通滤波器及实际RC 低通滤波电路的单位冲击响应与频谱图并记录试验结果。

五、实验结果（1）实验程序1.子程序（定义阶越函数）function f=heaviside(t)f=(t>0);2.主程序%理想低通滤波器的单位冲击响应、幅频特性、相频特性syms t f w;figure(1)f=sin(t-1)/(t-1); Fw=fourier(f); %傅里叶变换x=[-20:0.05:20]; fx=subs(f,t,x);subplot(2,1,1);plot(x,fx); %波形图grid;W=[-4:0.01:4];FW=subs(Fw,w,W);subplot(2,2,3);plot(W,abs(FW)); %幅频特性grid;xlabel('频率');ylabel('幅值');subplot(2,2,4);plot(W,angle(FW)); %相频特性grid;xlabel('频率');ylabel('相位');%RC低通网络的单位冲击响应、幅频特性、相频特性figure(2)f=exp(-2*t)*sym('heaviside(t)');Fw=fourier(f); %傅里叶变换x=[-4:0.02:4]; fx=subs(f,t,x);subplot(2,1,1);plot(x,fx); %波形图grid;W=[-4:0.02:4];FW=subs(Fw,w,W);subplot(2,2,3);plot(W,abs(FW)); %幅频特性grid;xlabel('频率');ylabel('幅值');subplot(2,2,4);plot(W,angle(FW)); %相频特性grid;xlabel('频率');ylabel('相位');（2）运行结果理想低通滤波器的单位冲击响应及频率特性RC低通滤波电路的单位冲击响应及频率特性六、思考题（1）理想低通滤波电路的幅频曲线和相频曲线有什么特点？（2）实际RC低通与理想低通滤波器的频谱有何不同？为什么？（3）在实验中的低通网络RC时间常数是多少？对低通滤波器有何影响？。

信号与系统实验报告

实验三常见信号的MATLAB 表示及运算一、实验目的1．熟悉常见信号的意义、特性及波形2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法二、实验原理根据MATLAB 的数值计算功能和符号运算功能,在MATLAB 中,信号有两种表示方法,一种是用向量来表示,另一种则是用符号运算的方法;在采用适当的MATLAB 语句表示出信号后,就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了;1.连续时间信号从严格意义上讲,MATLAB 并不能处理连续信号;在MATLAB 中,是用连续信号在等时间间隔点上的样值来近似表示的,当取样时间间隔足够小时,这些离散的样值就能较好地近似出连续信号;在MATLAB 中连续信号可用向量或符号运算功能来表示; ⑴ 向量表示法对于连续时间信号()f t ,可以用两个行向量f 和t 来表示,其中向量t 是用形如12::t t p t =的命令定义的时间范围向量,其中,1t 为信号起始时间,2t 为终止时间,p 为时间间隔;向量f 为连续信号()f t 在向量t 所定义的时间点上的样值; ⑵ 符号运算表示法如果一个信号或函数可以用符号表达式来表示,那么我们就可以用前面介绍的符号函数专用绘图命令ezplot 等函数来绘出信号的波形; ⑶ 常见信号的MATLAB 表示单位阶跃信号单位阶跃信号的定义为：10()0t u t t >⎧=⎨<⎩方法一：调用Heavisidet 函数首先定义函数Heavisidet 的m 函数文件,该文件名应与函数名同名即;%定义函数文件,函数名为Heaviside,输入变量为x,输出变量为y function y= Heavisidety=t>0; %定义函数体,即函数所执行指令%此处定义t>0时y=1,t<=0时y=0,注意与实际的阶跃信号定义的区别;方法二：数值计算法在MATLAB 中,有一个专门用于表示单位阶跃信号的函数,即stepfun 函数,它是用数值计算法表示的单位阶跃函数()u t ;其调用格式为：stepfunt,t0其中,t 是以向量形式表示的变量,t0表示信号发生突变的时刻,在t0以前,函数值小于零,t0以后函数值大于零;有趣的是它同时还可以表示单位阶跃序列()u k ,这只要将自变量以及取样间隔设定为整数即可; 符号函数符号函数的定义为：10sgn()1t t t >⎧=⎨-<⎩在MATLAB 中有专门用于表示符号函数的函数sign ,由于单位阶跃信号 t 和符号函数两者之间存在以下关系：1122()sgn()t t ε=+,因此,利用这个函数就可以很容易地生成单位阶跃信号;2.离散时间信号离散时间信号又叫离散时间序列,一般用()f k 表示,其中变量k 为整数,代表离散的采样时间点采样次数;在MATLAB 中,离散信号的表示方法与连续信号不同,它无法用符号运算法来表示,而只能采用数值计算法表示,由于MATLAB 中元素的个数是有限的,因此,MATLAB 无法表示无限序列；另外,在绘制离散信号时必须使用专门绘制离散数据的命令,即stem 函数,而不能用plot 函数; 单位序列()k δ单位序列()k δ的定义为10()0k k k δ=⎧=⎨≠⎩单位阶跃序列()u k单位阶跃序列()u k 的定义为10()0k u k k ≥⎧=⎨<⎩3.卷积积分两个信号的卷积定义为：MATLAB 中是利用conv 函数来实现卷积的;功能：实现两个函数1()f t 和2()f t 的卷积;格式：g=convf1,f2说明：f1=f 1t,f2=f 2t 表示两个函数,g=gt 表示两个函数的卷积结果;三、实验内容1.分别用MATLAB 的向量表示法和符号运算功能,表示并绘出下列连续时间信号的波形： ⑴ 2()(2)()tf t e u t -=- ⑵[]()cos()()(4)2tf t u t u t π=--1 t=-1::10;t1=-1::; t2=0::10;f1=zeros1,lengtht1,ones1,lengtht2;f=2-exp-2t.f1; plott,faxis-1,10,0, syms t;f=sym'2-exp-2theavisidet'; ezplotf,-1,10;2t=-2::8;f=0.t<0+cospit/2.t>0&t<4+0.t>4; plott,f syms t;f=sym'cospit/2heavisidet-heavisidet-4 '; ezplotf,-2,8;2.分别用MATLAB 表示并绘出下列离散时间信号的波形： ⑵ []()()(8)f t k u k u k =-- ⑶()sin()()4k f k u k π= 2 t=0:8; t1=-10:15;f=zeros1,10,t,zeros1,7; stemt1,faxis-10,15,0,10; 3 t=0:50; t1=-10:50;f=zeros1,10,sintpi/4; stemt1,faxis-10,50,-2,23.已知两信号1()(1)()f t u t u t =+-,2()()(1)f t u t u t =--,求卷积积分12()()()g t f t f t =*,并与例题比较;t1=-1::0; t2=0::1; t3=-1::1;f1=onessizet1; f2=onessizet2; g=convf1,f2;subplot3,1,1,plott1,f1; subplot3,1,2,plott2,f2; subplot3,1,3,plott3,g;与例题相比较,gt 的定义域不同,最大值对应的横坐标也不同;4.已知{}{}12()1,1,1,2,()1,2,3,4,5f k f k ==,求两序列的卷积和 ;N=4; M=5; L=N+M-1; f1=1,1,1,2;f2=1,2,3,4,5; g=convf1,f2; kf1=0:N-1; kf2=0:M-1; kg=0:L-1;subplot1,3,1,stemkf1,f1,'k';xlabel'k'; ylabel'f1k';grid onsubplot1,3,2,stemkf2,f2,'k';xlabel'k'; ylabel'f2k';grid onsubplot1,3,3;stemkg,g,'k';xlabel'k'; ylabel'gk';grid on 实验心得：第一次接触Mutlab 这个绘图软件,觉得挺新奇的,同时 ,由于之前不太学信号与系统遇到一些不懂的问题,结合这些图对信号与系统有更好的了解;实验四连续时间信号的频域分析一、实验目的1．熟悉傅里叶变换的性质 2．熟悉常见信号的傅里叶变换3．了解傅里叶变换的MATLAB 实现方法二、实验原理从已知信号()f t 求出相应的频谱函数()F j ω的数学表示为：()F j ω()j t f t e dt ω∞--∞=⎰傅里叶反变换的定义为：1()()2j t f t F j e d ωωωπ∞-∞=⎰在MATLAB 中实现傅里叶变换的方法有两种,一种是利用MATLAB 中的Symbolic Math Toolbox 提供的专用函数直接求解函数的傅里叶变换和傅里叶反变换,另一种是傅里叶变换的数值计算实现法;1.直接调用专用函数法①在MATLAB 中实现傅里叶变换的函数为：F=fourier f 对ft 进行傅里叶变换,其结果为Fw F =fourierf,v 对ft 进行傅里叶变换,其结果为Fv F=fourier f,u,v 对fu 进行傅里叶变换,其结果为Fv ②傅里叶反变换f=ifourier F 对Fw 进行傅里叶反变换,其结果为fx f=ifourierF,U 对Fw 进行傅里叶反变换,其结果为fu f=ifourier F,v,u 对Fv 进行傅里叶反变换,其结果为fu 注意：1在调用函数fourier 及ifourier 之前,要用syms 命令对所有需要用到的变量如t,u,v,w 等进行说明,即要将这些变量说明成符号变量;对fourier 中的f 及ifourier 中的F 也要用符号定义符sym 将其说明为符号表达式;2采用fourier 及fourier 得到的返回函数,仍然为符号表达式;在对其作图时要用ezplot 函数,而不能用plot 函数;3fourier 及fourier 函数的应用有很多局限性,如果在返回函数中含有δω等函数,则ezplot 函数也无法作出图来;另外,在用fourier 函数对某些信号进行变换时,其返回函数如果包含一些不能直接表达的式子,则此时当然也就无法作图了;这是fourier 函数的一个局限;另一个局限是在很多场合,尽管原时间信号ft 是连续的,但却不能表示成符号表达式,此时只能应用下面介绍的数值计算法来进行傅氏变换了,当然,大多数情况下,用数值计算法所求的频谱函数只是一种近似值;2、傅里叶变换的数值计算实现法严格说来,如果不使用symbolic 工具箱,是不能分析连续时间信号的;采用数值计算方法实现连续时间信号的傅里叶变换,实质上只是借助于MATLAB 的强大数值计算功能,特别是其强大的矩阵运算能力而进行的一种近似计算;傅里叶变换的数值计算实现法的原理如下：对于连续时间信号ft,其傅里叶变换为：其中τ为取样间隔,如果ft 是时限信号,或者当|t|大于某个给定值时,ft 的值已经衰减得很厉害,可以近似地看成是时限信号,则上式中的n 取值就是有限的,假定为N,有：若对频率变量ω进行取样,得：通常取：02k k k MM ωπωτ==,其中0ω是要取的频率范围,或信号的频带宽度;采用MATLAB 实现上式时,其要点是要生成ft 的N 个样本值()f n τ的向量,以及向量k j n eωτ-,两向量的内积即两矩阵的乘积,结果即完成上式的傅里叶变换的数值计算;注意：时间取样间隔τ的确定,其依据是τ必须小于奈奎斯特Nyquist 取样间隔;如果ft 不是严格的带限信号,则可以根据实际计算的精度要求来确定一个适当的频率0ω为信号的带宽;三、实验内容1.编程实现求下列信号的幅度频谱1 求出1()(21)(21)f t u t u t =+--的频谱函数F 1jω,请将它与上面门宽为2的门函数()(1)(1)f t u t u t =+--的频谱进行比较,观察两者的特点,说明两者的关系;2 三角脉冲21||||1()0||1t t f t t -≤⎧=⎨>⎩3 单边指数信号3()()tf t e t ε-=4 高斯信号23()t f t e -=1 syms t w Gt=sym'Heaviside2t+1-Heaviside2t-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi 0与()(1)(1)f t u t u t =+--的频谱比较,1()(21)(21)f t u t u t =+--的频谱函数F 1jω最大值是其的1/2; 2syms t w;Gt=sym'1+tHeavisidet+1-Heavisidet+1-tHeavisidet-Heavisidet-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi 0 3syms t w Gt=sym'exp-tHeavisidet';Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi -1 2 4syms t w Gt=sym'exp-t^2';Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; ezplotFFw,-30 30;grid; axis-30 30 -1 22.利用ifourier 函数求下列频谱函数的傅氏反变换122()16F j j ωωω=-+ 222()58()()65j j F j j j ωωωωω+-=++1syms t w Fw=sym'-i2w/16+w^2'; ft=ifourierFw,w,t; ft运行结果： ft =-exp4theaviside-t+exp-4theavisidet 2syms t wFw=sym'iw^2+5iw-8/iw^2+6iw+5'; ft=ifourierFw,w,t; ft运行结果： ft =diract+-3exp-t+2exp-5theavisidet实验心得matlab 不但具有数值计算能力,还能建模仿真,能帮助我们理解不同时间信号的频域分析;实验五连续时间系统的频域分析一、实验目的1. 学习由系统函数确定系统频率特性的方法;2. 学习和掌握连续时间系统的频率特性及其幅度特性、相位特性的物理意义;3.通过本实验了解低通、高通、带通、全通滤波器的性能及特点;二、实验原理及方法频域分析法与时域分析法的不同之处主要在于信号分解的单元函数不同;在频域分析法中,信号分解成一系列不同幅度、不同频率的等幅正弦函数,通过求取对每一单元激励产生的响应,并将响应叠加,再转换到时域以得到系统的总响应;所以说,频域分析法是一种变域分析法;它把时域中求解响应的问题通过 Fourier 级数或 Fourier 变换转换成频域中的问题；在频域中求解后再转换回时域从而得到最终结果;在实际应用中,多使用另一种变域分析法：复频域分析法,即 Laplace 变换分析法;所谓频率特性,也称频率响应特性,是指系统在正弦信号激励下稳态响应随频率变化的情况,包括幅度随频率的响应和相位随频率的响应两个方面;利用系统函数也可以确定系统频率特性,公式如下：幅度响应用()ωj H 表示,相位响应用)(ωϕH 表示;本实验所研究的系统函数Hs 是有理函数形式,也就是说,分子、分母分别是m 、n 阶多项式; 要计算频率特性,可以写出为了计算出()ωj H 、)(ωϕH 的值,可以利用复数三角形式的一个重要特性：而⎥⎦⎤⎢⎣⎡+=2sin 2cosππωωj j ,则()⎥⎦⎤⎢⎣⎡+=2sin 2cos ππωωn j n j n n利用这些公式可以化简高次幂,因此分子和分母的复数多项式就可以转化为分别对实部与虚部的实数运算,算出分子、分母的实部、虚部值后,最后就可以计算出幅度()ωj H 、相位)(ωϕH 的值了;三、实验内容a)sm m ms H )(1)(2-+=,m 取值区间 0,1,绘制一组曲线 m=,,,,； b) 绘制下列系统的幅频响应对数曲线和相频响应曲线,分析其频率特性; a %figurealpha=,,,,;colorn='r' 'g' 'b' 'y' 'k'; % r g b y m c k 红,绿,蓝,黄,品红,青,黑 for n=1:5b=0 alphan; % 分子系数向量a=alphan-alphan^2 1; % 分母系数向量 printsysb,a,'s' Hz,w=freqsb,a; w=w./pi; magh=absHz;zerosIndx=findmagh==0; maghzerosIndx=1; magh=20log10magh; maghzerosIndx=-inf; angh=angleHz;angh=unwrapangh180/pi; subplot1,2,1plotw,magh,colornn;hold onsubplot1,2,2plotw,angh,colornn;hold onendsubplot1,2,1hold offxlabel'特征角频率\times\pi rad/sample' title'幅频特性曲线 |Hw| dB';subplot1,2,2hold offxlabel'特征角频率 \times\pi rad/sample' title'相频特性曲线 \thetaw degrees';b1 %b=1,0; % 分子系数向量a=1,1; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';2 %b=0,1,0; % 分子系数向量a=1,3,2; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';3 %b=1,-1; % 分子系数向量a=1,1; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';实验心得：虽然之前用公式转换到频域上分析,但是有时会觉得挺抽象的,不太好理解;根据这些图像结合起来更进一步对信号的了解;同时,这个在编程序时,虽然遇到一些问题,但是总算解决了;实验六离散时间系统的Z域分析一、实验目的1. 学习和掌握离散系统的频率特性及其幅度特性、相位特性的物理意义;2. 深入理解离散系统频率特性和对称性和周期性;3. 认识离散系统频率特性与系统参数之间的系统4.通过阅读、修改并调试本实验所给源程序,加强计算机编程能力; 二、实验原理及方法对于离散时间系统,系统单位冲激响应序列)(n h 的 Fourier 变换)(ωj e H 完全反映了系统自身的频率特性,称)(ωj eH 为离散系统的频率特性,可由系统函数)(z H 求出,关系式如下：ωωj j e z z H e H ==)()( 6 – 1由于ωj e是频率的周期函数,所以系统的频率特性也是频率的周期函数,且周期为π2,因此研究系统频率特性只要在πωπ≤≤-范围内就可以了;∑∑∑∞-∞=∞-∞=∞-∞=--==n n n j j n n h j n n h en h e H )sin()()cos()()()(ωωωω6 – 2容易证明,其实部是ω的偶函数,虚部是ω的奇函数,其模ωj e H (的ω的偶函数,相位[])(arg ωj e H 是ω的奇函数;因此研究系统幅度特性)(ωj e H 、相位特性[])(arg ωj e H ,只要在πω≤≤0范围内讨论即可;综上所述,系统频率特性)(ωj eH 具有周期性和对称性,深入理解这一点是十分重要的;当离散系统的系统结构一定,它的频率特性)(ωj e H 将随参数选择的不同而不同,这表明了系统结构、参数、特性三者之间的关系,即同一结构,参数不同其特性也不同; 例如,下图所示离散系统,其数学模型由线性常系数差分方程描述：)()1()(n x n ay n y +-=系统函数：a z az z z H >-=,)(系统函数频率特性：ωωωωωsin )cos 1(1)(ja a a e e e H j j j +-=-=幅频特性：ωωcos 211)(2a a eH j -+=相频特性：[]ωωωcos 1sin arctan)(arg a a eH j --= 容易分析出,当10<<a 时系统呈低通特性,当01<<-a 时系统呈高通特性；当0=a 时系统呈全通特性;同时说明,在系统结构如图所示一定时,其频率特性随参数a 的变化而变化;三、实验内容a 2281.011)(----=z z z H ;b 1.04.06.01.03.03.01.0)(2323+++-+-=z z z z z z z Hc 2181.011)(--+-=zz z H a %b=1,0,-1; % 分子系数向量a=1,0,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';带通b %b=,,,; % 分子系数向量a=1,,,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';高通c %b=1,-1,0; % 分子系数向量a=1,0,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';带通实验心得：本来理论知识不是很强的,虽然已经编出程序得到相关图形,但是不会辨别相关通带,这让我深刻地反省;。

《信号与系统》实验报告

信号与系统实验报告班级:姓名:信息与通信工程学院实验一系统的卷积响应实验性质：提高性实验级别：必做开课单位：信息与通信工程学院学时：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 )(0t h t h t P t =→=∆→∆∆→∆δ若输入为f(t):∆∆-∆=≈∑∞-∞=∆∆)()()()(k t P k f t f t f k得输出：∆∆-∆=∑∞-∞=∆∆)()()(k t hk 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 hk f t y t y k )()()()(lim)(lim )(0所以：∆∆-∆=-==∑⎰→∆)()(lim)()()(*)()(212121k t f k fd t f f t f t f t s τττ如果只求离散点上的f 值)(n f ∆])[()()()()(2121∑∑∞-∞=∞-∞=∆-∆∆=∆∆-∆∆=∆k k k n f k f k n f k fn f所以，可以用离散卷积和CONV （）求连续卷积，只需∆足够小以及在卷积和的基础上乘以∆。

信号与系统实验报告

《信号与系统》实验报告湖南工业大学电气与信息工程学院实验一用同时分析法观测50Hz非正弦周期信号的分解与合成一、实验目的1、用同时分析法观测50Hz非正弦周期信号的频谱，并与傅立叶级数各项的频率与系数作比较。

2、观测基波和其谐波的合成。

二、实验设备1、信号与系统实验箱：TKSS －Ａ型或TKSS －B 型TKSS －C 型；2、双踪示波器三、实验原理1、一个非正弦周期函数可以用一系列频率成整数倍的正弦函数来表示，其中与非正弦具有相同频率的成分称为基波或一次谐波，其他成分则根据其频率为基波频率的2、3、4、…、n 等倍数分别称为二次、三次、四次、…、n 次谐波，其幅度将随着谐波次数的增加而减小，直至无穷小。

2、不同频率的谐波可以合成一个非正弦周期波，反过来，一个非正弦周期波也可以分解为无限个不同频率的谐波成分，3、一个非正弦周期函数可以用傅立叶级数来表示，级数各项系数之间的关系可用一个频谱来表示，不同的非正弦周期函数具有不同的频谱图，各种不同波形及其傅氏级数表达式见表2-1，方波频谱图如图2-1表示Um1351/91/51/71/3790ωωωωωω图1-1 方波频谱图表2-1 各种不同波形的傅立叶级数表达式UmtTU 2τ方波UmTU 2τ正弦整流全波UmTU 2τ三角波Um0T2τ正弦整流半波t tUm0tT U 2τ矩形波U1、方波 ())7sin 715sin 513sin 31(sin 4 ++++=t t t t u t u m ωωωωπ 2、三角波())5sin 2513sin 91(sin 82++-=t t t u t u mωωωπ3、半波())4cos 1512cos 31sin 421(2 +--+=t t t u t u m ωωωππ 4、全波 ())6cos 3514cos 1512cos 3121(4 +---=t t t u t u m ωωωπ5、矩形波())3cos 3sin 312cos 2sin 21cos (sin 2 ++++=t T t T t T U T U t u m m ωτπωτπωτππτ实验装置的结构如图1-2所示DC20f f f f f f 3456图1-2信号分解于合成实验装置结构框图图中LPF 为低通滤波器，可分解出非正弦周期函数的直流分量。

信号与系统实验报告-1(常用信号的分类与观察)

信号与系统实验报告-1(常用信号的分类与观察)实验一：信号的时域分析一、实验目的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)为单位阶跃函数。

其信号如下图所示：f(t)…………0 t图1－6 脉冲信号7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示U(t)…………0 t图1－7 方波信号四、实验内容及主要步骤下列实验中信号产生器的工作模式为111、指数信号观察通过信号选择键1，设置A组输出为指数信号（此时信号输出指示灯为000000）。

信号与系统实验-周期信号的合成与分解

武汉大学教学实验报告电子信息学院专业年月日实验名称周期信号的合成与分解指导教师姓名年级学号成绩一、预习部分1.实验目的2.实验基本原理3.主要仪器设备（含必要的元器件、工具）一、实验目的1．在理论学习的基础上，通过实验深刻领会周期信号傅里叶级数分解的物理意义。

2．理解实际应用中通常采用有限项级数来逼近无限项级数，此时方均误差随项数的增加而减小。

3．观察并初步了解 Gibbs 现象。

4．深入理解周期信号的频谱特点，比较不同周期信号频谱的差异。

二、实验原理满足 Dirichlet 条件的周期信号f(t)可以分解成三角函数形式的傅里叶级数，表达式为：式中n为正整数；角频率ω1由周期T1决定：112=Tπω。

该式表明：任何四、实验内容1．周期对称方波信号的合成图示方波既是一个奇对称信号，又是一个奇谐信号。

根据函数的对称性与傅里叶系数的关系可知，它可以用无穷个奇次谐波分量的傅里叶级数来表示:选取奇对称周期方波的周期T=0.02s，幅度E =6，请采用有限项级数替代无限项级数来逼近该函数。

分别取前 1、2、5 和 100 项有限级数来近似，编写程序并把结果显示在一幅图中，观察它们逼近方波的过程。

MATLAB 程序如下：%奇对称方波合成f(t)=2*E/pi*( sum(sin(2*pi*f*t*(2*i+1)/(2*i+1));i=0,1,2¡-¡-)t=0:0.001:0.1;E=6;A=2*E/pi;T=0.02;f=1/T;y=A*sin(2*pi*f*t);subplot(221)plot(t,y);显示结果如图 4-2 所示。

图 4-2 奇对称方波信号的合成2．观察 Gibbs 现象分别取前 5,6,7和8项有限级数来逼近奇对称方波，观察 Gibbs 现象。

MATLAB 程序如下：%观察 Gibbs 现象%奇对称方波合成，观察 Gibbs 现象t=0:0.001:0.04;sishu=12/pi;y=0;for i=1:5y=y+sishu*(sin((2*i-1)*100*pi*t)/(2*i-1));endsubplot(221)plot(t,y);axis([0,0.04,-4,4]);xlabel('time');ylabel('前 5 项有限级数);y=0;图 4-3 Gibbs 现象3．周期对称三角信号的合成设计采用有限项级数逼近偶对称周期三角信号的实验，编制程序并显示结果。