信号与系统实验2
电路与信号系统实验报告2-二阶电路的动态响应
二阶电路的动态响应
一、实验内容
1.Multisim仿真
（1）创建电路：
（2）设置=10mH、=22nF，电容的初始电压为5V，电源电压为10V。
利用Transient Analysis观测电容两端的电压。
（3）用Multisim瞬态分析仿真零输入响应（参数欠阻尼、临界阻尼、过阻尼三
种情况）；在同一张图中画出三条曲线，标出相应阻值。
（4）用Multisim瞬态分析仿真零响应（参数欠阻尼、临界阻尼、过阻尼三种情况）；在同一张图中画出三条曲线，标出相应阻值。
（5）利用Multisim中的函数发生器、示波器和波特仪创建如图所示的电路，观测各种响应。
函数信号发生器设置：方波、频率1kHz、幅度5V、偏置5V。
2.在电路板上按图所示的电路（R1=100Ω、L=10mH、C=47nF）焊接实验电路。
3.调节可变电阻器R2，观察二阶电路的零输入响应和零状态响应由过阻尼过渡到临界阻尼，最后过渡到欠阻尼的变化过渡过程，分别定性地描绘、记录响应的典型变化波形，按表记录所测数据和波形。
欠阻尼：R2=35Ω
临界阻尼：R2=819Ω，衰减时间=100us
过阻尼：R2=1kΩ，衰减时间=160us
4.调节R2使示波器荧光屏上呈现稳定的欠阻尼响应波形，定量测定此时电路的衰减振荡角频率
二、实验结论
电阻越高，响应衰减地越快；电阻越低，响应衰减地越慢。
当电阻为零时，
电路无衰减。
信号与系统实验报告3-(2)
信号与系统实验实验三：信号的卷积（第三次实验）【实验目的】1. 理解卷积的物理意义；2. 掌握运用计算机进行卷积运算的原理和方法；3. 熟悉卷积运算函数conv的应用；【实验内容】给定如下因果线性时不变系统：y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3](1)不用impz函数，使用filter命令，求出以上系统的单位冲激响应h[n]的前20个样本;代码如下：clear all;N=[0:19];num=[0.9 -0.45 0.35 0.002];den=[1 0.71 -0.46 -0.62];h=filter(num,den,N);stem(N,h);xlabel('ʱ¼äÐòºÅ');ylabel('Õñ·ù');title('µ¥Î»³å¼¤ÏìÓ¦');grid;图像如下：(2)得到h[n]后，给定x[n],计算卷积输出y[n]；并用滤波器h[n]对输入x[n]滤波，求得y1[n];代码如下：clear all;N=[0:19];num=[0.9 -0.45 0.35 0.002];den=[1 0.71 -0.46 -0.62];h=filter(num,den,N);x=[1 -2 3 -4 3 2 1];y=conv(h,x);n=0:25;subplot(2,1,1);stem(n,y);xlabel('时间序号n');ylabel('振幅');title('用卷积得到的输出');grid;x1=[x zeros(1,19)];y1=filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel('时间序号n ');ylabel('振幅');title('用滤波得到的输出');grid;图像如下：(3)y[n]和)y1[n]有差别吗？为什么要对x[n]进行补零得到的x1[n]来作为输入来产生y1[n]？(4)思考：设计实验，证明下列结论① 单位冲激信号卷积：)()(*)(t f t f t =δ)()(*)(00t t f t f t t -=-δ代码如下：clc;clear all ;n=[0:20];d=(n==0);f=sin(n);f1=conv(d,f);subplot(3,1,1);f1=f1(1:21);stem(n,f1);title('¦Ä[n]*f[n]');grid;subplot(3,1,2);stem(n,f);title('f[n]');grid;subplot(3,1,3);stem(n,f-f1);title('¦Ä[n]*f[n]-f[n]');grid;图像如下：② 卷积交换律：)(*)()(*)()(1221t f t f t f t f t f ==代码如下：clc;clear all;n=0:30;f1=sin(n);f2=cos(n);y1=conv(f1,f2);y1=y1(1:31);y2=conv(f2,f1);y2=y2(1:31); subplot(3,1,1); stem(n,y1);title('f1*f2'); grid;subplot(3,1,2); stem(n,y2);title('f2*f1'); grid;subplot(3,1,3);y3=(y1-y2)>10^-14; stem(n,y3);grid;图像如下：③卷积分配律：)(*)()(*)()]()([*)(3121321t f t f t f t f t f t f t f +=+代码如下：clc;clear all ;n=1:50;f1=(-1).^n;f2=cos(n);f3=sin(n);y1=conv(f1,(f2+f3));y1=y1(1:50);y2=conv(f1,f2)+conv(f1,f3);y2=y2(1:50);subplot(3,1,1);stem(n,y1);title('f1*[f2+f3]');grid;subplot(3,1,2);stem(n,y2);title('f1*f2+f1*f3');grid;subplot(3,1,3);y3=(y1-y2)>10^-14;stem(n,y3);title('f1*[f2+f3]-f1*f2+f1*f3');grid;图像如下：【实验分析】：1.y[n]和)y1[n]有差别吗？为什么要对x[n]进行补零得到的x1[n]来作为输入来产生y1[n]？答：y[n]和)y1[n]是对同一个系统输入的响应，该系统是因果线性时不变系统，所以y[n]和)y1[n]没有差别；由于y[n]和)y1[n]没有差别，滤波器h[n]对x[n]滤波得到的y1[n]和用卷积计算得到的y[n]是同一个信号；2.卷积分配率程序代码中f1的n时间序号长度n为[1:50]，f2的n时间序号长度为[1:50]，所以输出完整信号的长度为99，而程序中输出长度仅50，说明这只是信号的部分波形。
信号与线性系统实验报告2
实验二连续系统频域分析一、实验目的1．通过观察信号的分解与合成过程，理解利用傅利叶级数进行信号频谱分析的方法。
2．了解波形分解与合成原理。
3．掌握带通滤波器有关特性的设计和测试方法。
4．了解电信号的取样方法与过程以及信号恢复的方法。
5．观察连续时间信号经取样后的波形图，了解其波形特点。
6．验证取样定理并恢复原信号。
二、实验内容1．用示波器观察方波信号的分解，并与方波的傅利叶级数各项的频率与系数作比较。
2．用示波器观察三角波信号的分解，并与三角波的傅利叶级数各项的频率与系数作比较。
3．用示波器观察方波信号基波及各次谐波的合成。
4．用示波器观察三角波信号基波及各次谐波的合成。
5．用示波器观察不同的取样频率抽样得到的抽样信号。
6．用示波器观察各取样信号经低通滤波器恢复后的信号并验证抽样定理。
三、实验仪器1．信号与系统实验箱一台2．信号系统实验平台3．信号的分解与合成模块（DYT3000-69）一块4．信号的取样与恢复模块（DYT3000-68）一块5．同步信号源模块（DYT3000-57）（选用）6．20MHz双踪示波器一台7．连接线若干四、实验原理1、信号的分解与合成任何电信号都是由各种不同频率、幅度和初始相位的正弦波跌加而成的。
对周期信号由它的傅利叶级数展开式可知，各次谐波为基波频率的整数倍。
而非周期信号包含了从零到无穷大的所有频率成份，每一频率成份的幅度均趋向无穷小，但其相对大小是不同的。
通过一个选频网络可以将电信号中所包含的某一频率成份提取出来。
本实验采用性能较好的有源带通滤波器作为选频网络。
对周期信号波形分解的方案框图如图2-1所示。
实验中对周期方波、三角波、锯齿波信号进行信号的分解。
方波信号的傅利叶级数展开式为411()(sin sin 3sin 5)35Af t t t t ωωωπ=+++…；三角波信号的傅利叶级数展开式为2811()(sin sin 3sin 5)925A f t t t t ωωωπ=-+-…；锯齿波信号的傅利叶级数展开式为11()(sin sin 2sin 3)223A A f t t t t ωωωπ=-+++…，其中2T πω=为信号的角频率。
电路、信号与系统(2)实验指导书
描述线性时不变离散系统的差分方程为
编写求解上述方程的通用程序。
[建模]
将方程变形可得（用MATLAB语言表示）
a(1)*y(n)= b(1)*u(n)+…+ b(nb)*u(n-nb+1)- a(2)*y(n-1)-…- a(na)*y(n-na+1)
令us== [u(n),…, u(n-nb+1)]; ys=[y(n-1),…, y(n-na+1)]
x(n)={2,1,-1,3,1,4,3,7}（其中加下划线的元素为第0个采样点）在MATLAB中表示为：
n=[-3,-2,-1,0,1,2,3,4]; x=[2,1,-1,3,1,4,3,7]；
当不需要采样位置信息或这个信息是多余的时候，可以只用x向量来表示。
（一）离散信号的MATLAB表述
[问题]
实验一连续时间信号与系统分析
一、实验目的
1、了解连续时间信号的特点；
2、掌握连续时间信号的MATLAB描述；
3、掌握连续LTI系统单位冲激响应的求解方法；
4、掌握连续LTI系统的零状态响应的求解方法。
二、实验内容
严格说来，只有用符号推理的方法才能分析连续系统，用数值方法是不能表示连续信号的，因为它给出的是各个样点的数据。只有当样本点取得很密时才可看成连续信号。所谓很密，是相对于信号变化的快慢而言的。以下均假定相对于采样点密度而言，信号变化足够慢。
elseif lu<lh nh=0; nu=lh-lu;
else nu=0; nh=0;
end
dt=0.1;
lt=lmax;
u=[zeros(1, lt), uls, zeros(1, nu), zeros(1, lt)];
信号与系统实验实验2常用离散时间信号的实现
信号与系统实验实验2常用离散时间信号的实现信号与系统是电子信息类专业的一门基础课程，是理论与实践相结合的一门课程。
离散时间信号与系统是信号与系统理论的一个重要分支，是实际工程应用中的基础。
本实验主要目的是通过实际操作，实现常用离散时间信号的生成和处理，加深对离散时间信号与系统的理解。
实验一：离散时间单位阶跃信号的生成和显示实验介绍：离散时间单位阶跃信号是离散时间系统的基本信号之一，表示时间从0开始，幅值从0突变到1的信号。
本实验通过编写Matlab程序，实现离散时间单位阶跃信号的生成和显示。
实验步骤：1. 打开Matlab软件，创建一个新的脚本文件。
2.在脚本文件中编写以下程序代码：```matlab%生成离散时间单位阶跃信号n=0:10;%离散时间序列u = ones(1,11); % 生成11个单位阶跃信号的幅值stem(n, u); % 显示离散时间单位阶跃信号title('Unit Step Signal'); % 设置图像标题```3.运行程序，得到离散时间单位阶跃信号的图像及其数值序列。
4.分析实验结果，比较离散时间单位阶跃信号与连续时间单位阶跃信号的区别。
实验二：离散时间指数信号的生成和显示实验介绍：离散时间指数信号是离散时间系统中常见的信号之一，表示时间以指数形式变化的信号。
本实验通过编写Matlab程序，实现离散时间指数信号的生成和显示。
实验步骤：1. 打开Matlab软件，创建一个新的脚本文件。
2.在脚本文件中编写以下程序代码：```matlab%生成离散时间指数信号n=0:10;%离散时间序列a=0.8;%指数信号的衰减系数x=a.^n;%生成离散时间指数信号的幅值stem(n, x); % 显示离散时间指数信号title('Exponential Signal'); % 设置图像标题```3.运行程序，得到离散时间指数信号的图像及其数值序列。
信号与系统-2
Y zi ( z ) =
z2
2z 2z = 1 1 1 1 − z − z − z + 6 6 2 3
2
2
部分分式展开，得
0 .8 z 1 .2 z Y zi ( z ) = + 1 1 z− z+ 3 2
所以，零输入响应为
1 1 y zi ( n ) = 1 . 2 + 0 . 8 − 2 3
s→∞
得到 f (0 + ) = ∞ 的 错误 结论。 由 于 F(s) 的结点全在左半 s 平 面 ， 故可 利用终值定理求终值。
s ( s 3 + s 2 + 2 s + 1) lim f (t ) = lim sF ( s ) = lim =0 t→∞ s→∞ s → ∞ ( s + 1)( s + 2)( s + 3)
t →∞
解：(1)因为 F(s)是真分式，所以
s + 2s + 1 f (0+ ) = lim sF ( s) = lim s =1 s →∞ s →∞ ( s − 1)( s + 2)( s + 3)
2
由于 F(s)在右半 s 平面的极点 s=1,故 f(t)的终值 不存在。
(2) 由 于 F(s) 分子的阶次等于分母的阶次， 故 利用长除法， 使得
−t
已知信号的拉氏变换函数如下， 求原信号的初 值与终值
s + 2s + 1 (1) F ( s ) = ( s − 1)( s + 2)( s + 3)
2
s + s + 2s + 1 (2) F ( s ) = ( s + 1)( s + 2)( s + 3) 1 (3) F ( s ) = s ( s + 1)
信号实验报告 2
信号与系统实验报告实验一信号与系统的时域分析一、实验目的1、熟悉和掌握常用的用于信号与系统时域仿真分析的MA TLAB函数；2、学会用MA TLAB进行信号基本运算的方法；3、掌握连续时间和离散时间信号的MA TLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程。
二、实验内容Q1-1：修改程序Program1_1，将dt改为0.2，再执行该程序，保存图形，看看所得图形的效果如何？dt = 0.01时的程序clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.01; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = sin(2*pi*t); % Generate the signalplot(t,x) % Open a figure window and draw the plot of x(t)title('Sinusoidal signal x(t)')xlabel('Time t (sec)')dt = 0.2时的程序clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.2; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = sin(2*pi*t); % Generate the signalplot(t,x) % Open a figure window and draw the plot of x(t)title('Sinusoidal signal x(t)')xlabel('Time t (sec)')dt = 0.01时的信号波形dt = 0.2时的信号波形这两幅图形有什么区别，哪一幅图形看起来与实际信号波形更像？答：dt = 0.01的图形比dt = 0.2的图形光滑，dt = 0.01看起来与实际信号波形更像。
实验2 信号卷积实验
实验二 信号卷积实验一、实验目的1. 理解卷积的概念及物理意义；2. 通过实验的方法加深对卷积运算的图解方法及结果的理解。
二、实验设备1. 信号与系统实验箱 1台2. 双踪示波器1台3. 铆孔连接线 若干二、实验原理说明卷积积分的物理意义是将信号分解为冲激信号之和，借助系统的冲激响应，求解系统对任意激励信号的零状态响应。
设系统的激励信号为)t (x ，冲激响应为)t (h ，则系统的零状态响应为)(*)()(t h t x t y =()()x t h t d ττ∞-∞=-⎰。
对于任意两个信号)t (f 1和)t (f 2，两者做卷积运算定义为：()()()12f t f t f t d ττ∞-∞=-⎰=)t (f 1*)t (f 2=)t (f 2*)t (f 1。
1. 两个矩形脉冲信号的卷积过程两信号)t (x 与)t (h 都为矩形脉冲信号，如图10-1所示。
下面由图解的方法（图10-1）给出两个信号的卷积过程和结果，以便与实验结果进行比较。
0≤<∞-t210≤≤t 1≤≤t 41≤≤t ∞<≤t 2124τ（b)(a)(c)(d)(e)(f)(g)(h)(i)2卷积结果图10-1 两矩形脉冲的卷积积分的运算过程与结果2. 矩形脉冲信号与锯齿波信号的卷积信号)t(f1为矩形脉冲信号，)t(f2为锯齿波信号，如图10-2所示。
根据卷积积分的运算方法得到)t(f1和)t(f2的卷积积分结果)t(f，如图10-2(c)所示。
(a)(b)(c)图10-2 矩形脉冲信号与锯齿脉冲信号的卷积积分的结果3. 本实验进行的卷积运算的实现方法在本实验装置中采用了DSP数字信号处理芯片，因此在处理模拟信号的卷积积分运算时，是先通过A/D转换器把模拟信号转换为数字信号，利用所编写的相应程序控制DSP芯片实现数字信号的卷积运算，再把运算结果通过D/A转换为模拟信号输出。
结果与模拟信号的直接运算结果是一致的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
信号与系统实验
信息学院通信工程 20101060163 荣华杰
e t-
8.1已知描述连续时间系统的微分方程和激烈信号f（t）为y''(t)+4y'(t)+2y(t)=f''(t)+3f(t),f(t)=)(t u
试用MATLAB的lsim函数求出上述系统在0-10秒时间零状态响应y（t）的值，并绘出系统零状态响应的时域仿真波形。
e t-
y''(t)+4y'(t)+2y(t)=f''(t)+3f(t),f(t)=)(t u
a=[1 4 2];
b=[1 0 3];
sys=tf(b,a);
t=0:1:10;
f=exp(-(t));
y=lsim(sys,f,t)
y = 1.0000
0.1721
0.3629
0.3138
0.2162
0.1356
0.0811
0.0472
0.0271
0.0153
0.0086
波形图：
a=[1 4 2];
b=[1 3];
sys=tf(b,a);
t=0:0.01:10;
f=exp(-2*t);
lsim(sys,f,t)
8.3已知描述系统的微分方程如下，试用MATLAB求系统在0-10秒时间范围内冲激响应和阶跃响应的数值解，并绘出系统冲击响应和阶跃响应的时域波形。
（3）y''(t)+4y'(t)+5y(t)=f'(t)
冲激响应数值解：
a=[1 4 5];
b=[1 0];
y=impulse(b,a,0:1:10)
y =1.0000
-0.1546
-0.0409
-0.0032
0.0003
0.0001
0.0000
-0.0000
-0.0000
-0.0000
0.0000
阶跃响应数值解：
a=[1 4 5];
b=[1 0];
y=step(b,a,0:1:10)
y = 0
0.1139
0.0167
0.0003
-0.0003
-0.0000
-0.0000
0.0000
0.0000
-0.0000
冲激响应、阶跃响应的时域波形：
a=[1 4 5];
b=[1 0];
subplot(1,2,1)
step(b,a,10)
subplot(1,2,2)
impulse(b,a,10)
8.4已知描述离散系统的差分方程和输入系列x（n）如下，试用MATLAB的filter函数求出上述系统在0-20时间采样点范围内零状态响应y（n）的系列样值，并绘出系统零状态响应的时域波形。
（2）y(n)+(1/2)y(n-1)=x(n)+2x(n-1),x(n)=2cos(n*pi/3)*u(n)
零状态响应y（n）的系列样值:
a=[1 1/2];
b=[1 2];
n=0:20;
x=2*cos(n*(pi/3));
y=filter(b,a,x)
y =Columns 1 through 7
2.0000 4.0000 -1.0000 -
3.5000 -3.2500 0.6250 3.6875
Columns 8 through 14
3.1563 -0.5781 -3.7109 -3.1445 0.5723 3.7139 3.1431
Columns 15 through 21
-0.5715 -3.7142 -3.1429 0.5714 3.7143 3.1429 -0.5714
波形图：
a=[1 1/2];
b=[1 2];
n=0:0.1:20;
x=2*cos(n*(pi/3));
y=filter(b,a,x);
stem(n,y,'filled')
响应和阶跃响应的数值解，绘出其序列波形图，并根据单位序列响应的时域波形判断系统的稳定性。
（3）y(n)+y(n-1)+(1/4)y(n-2)=x(n)
单位序列响应和阶跃响应的数值解:
a=[1 1 1/4];
b=[1];
h=impz(b,a,0:20)
y=step(b,a,0:2:20)
h =1.0000
-1.0000
0.7500
-0.5000
0.3125
-0.1875
0.1094
-0.0625
0.0352
-0.0195
0.0107
-0.0059
0.0032
-0.0017
0.0009
-0.0005
0.0003
-0.0001
-0.0000
0.0000
y = 0
1.0570
2.3760
3.2034
3.6337
3.8383
3.9306
3.9708
3.9879
3.9951
3.9980
波形图：
a=[1 1 1/4];
b=[1];
subplot(1,2,1)
impz(b,a,0:20)
subplot(1,2,2)
step(b,a,0:2:20)
8.7已知LTI离散系统的单位序列响应h（n）和激烈x（n）如图8-29（a）所示，试用MATLAB的conv函数求出系统的零状态响应y（n），并绘出其时域波形。
零状态响应y值：
x=[0 1 2 1 0 0];
h=[0 1 1 1 1 0 0];
y=conv(x,h)
y =Columns 1 through 11
0 0 1 3 4 4 3 1 0 0 0
波形图：
n1=-1:5;
n2=-2:3;
n=(n1(1)+n2(1)):(n1(1)+n2(1)+length(n1)+length(n2)-2);
subplot(2,2,1)
stem(n1,x,'filled')
title('x(n)')
xlabel('n')
subplot(2,2,2)
stem(n2,h,'filled')
title('h(n)')
xlabel('n')
subplot(2,2,3)
stem(n,y,'filled')
title('y(n)')
xlabel('n')
绘出卷积
积分信号的时域波形，将其与解析计算结果进行比较。
（2）13f t f t （）*（）
t1=-1:0.01:3;
f1=heaviside(t1)-heaviside(t1-2);
t2=-1:0.01:3;
f2=0.5*t2.*(heaviside(t2)-heaviside(t2-2));
[t,f]=gggfconv(f1,f2,t1,t2);
请输入时间间隔d:5
f = 0 1.2500 3.7500 3.7500 1.2500
（5）34f t f t （）*（）
t3=-1:0.01:3;
f3=t3/2.*(Heaviside(t3)-Heaviside(t3-2));
t4=-3:0.01:3;
f4=(t4/2+1).*(Heaviside(t4+2)-Heaviside(t4))+(1-t4/2).*(Heaviside(t4)-Heaviside(t4-2));
[t,f]=gggfconv(f3,f4,t3,t4);。