信号与系统实验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转换为模拟信号输出。
结果与模拟信号的直接运算结果是一致的。
信号与系统matlab作业2：对任意两段声音的卷积
题目：录两段不同的声音（语音、音乐、噪声均可），时间自己设定，然后再将这两段声音卷积。
要求：分别播放出每段声音（自己听）、分别显示每段声音的波形，从声音和波形两方面理解卷积的作用。
解：先用MATLAB录制两段不同的声音（即附件里的shengyin1,shengyin2）,其音频时间均为5秒，采样频率为40000Hz。
然后再用MATLAB画出两段声音及卷积后声音的波形图（如图【1】），并保存下两段卷积后的声音（即附件里的juanjishengyin）。
图【1】图【1】的MATALB程序：clc,clear[Y1,fs1]=audioread('shengyin1.wav');%获取音频1文件，并返回采样数据到向量y 中，fs表示采样频率[Y2,fs2]=audioread('shengyin2.wav');%获取音频1文件，并返回采样数据到向量y 中，fs表示采样频率ft1 = Y1(:,1);sigLength1 = length(ft1); %获取声音长度t1=(0:sigLength1-1)/fs1; %求出音频1对应的时间坐标subplot(3,1,1);plot(t1,ft1);%画出音频1的波形图title('音频1的波形图');xlabel('时间(s)');ylabel('振幅');grid;ft2 = Y2(:,1);sigLength2 = length(ft2); %获取声音长度t2=(0:sigLength2-1)/fs2; %求出音频2对应的时间坐标subplot(3,1,2);plot(t2, ft2);%画出音频2的波形图title('音频2的波形图');xlabel('时间(s)');ylabel('振幅');grid;ft=conv(ft1,ft2); %将两段音频进行卷积sigLength = length(ft); %获取声音长度t=(0:sigLength-1)/fs1;subplot(3,1,3);plot(t,ft);%画出卷积后的波形图title('音频1和音频2卷积后的波形图');xlabel('时间(s)');ylabel('振幅');grid;wavwrite(ft,40000,16,'juanjishengyin');%保存卷积后的声音总结：通过这次MATALB作业，观察三个不同的波形图（如图【1】）可知，卷积的作用就是一个信号函数在另一个信号函数上的加权叠加，通俗的说，它就是在输入信号的每个位置，叠加一个单位响应，就得到了输出信号。
2010信号与系统实验2
实验二：连续和离散系统的频域分析一：实验目的1：学习傅里叶正变换和逆变换，理解频谱图形的物理含义2：了解连续和离散时间系统的单位脉冲响应3：掌握连续时间系统的频率特性二：实验原理1. 傅里叶正变换和逆变换公式 正变换：()()j t F f t e dt ωω∞--∞=⎰逆变换：1()()2j t f t F e d ωωωπ∞-∞=⎰2. 频域分析t j tj e d d e t e ωωωπωωωπ⎰⎰∞∞-∞∞-E =E =)(21)(21)(将激励信号分解为无穷多个正弦分量的和。
⎰∞∞-H E =ωωωπωd e t r tj zs )()(21)(，R(ω)为)(t r zs 傅里叶变换;πωωd )(E 各频率分量的复数振幅 激励单位冲激响应时的零状态响应→ )(t δ)(t h单位阶跃响应时的零状态响应激励→)(t u )(t g3 各函数说明：(1)impulse 冲激响应函数：[Y ,X,T]=impulse(num,den);)1()2()1()1()2()1()()()(11++++++++==--n a s a s a m b s b s b s A s B s H n n m m num 分子多项式系数； num=[b(1) b(2) … b(n+1)]; den 分母多项式系数； den=[a(1) a(2) … a(n+1)];Y ,X,T 分别表示输出响应，中间状态变量和时间变量； 如：352)(2+++=s s s s H ，等价于)(2)()(3)(5)(t e t e t r t r t r +=++ 定义den=[1 5 3];num=[1 2]; [Y ,X,T]=impulse(num,den);(2)step 阶跃响应函数：[Y,X,T]=step(num,den);num 分子多项式；den 分母多项式 Y ,X,T 分别表示输出响应，中间状态变量和时间变量；如：352)(2+++=s s s s H ，den=[1 5 3];num=[1 2];[Y ,X,T]= step (num,den);（3）impz 数字滤波器的冲激响应 [h,t] = impz(b,a,n) b 分子多项式系数；a 分母多项式系数；n 采样样本h 离散系统冲激响应；t 冲激时间，其中t=[0:n-1]', n=length(t)时间样本数（4）freqs 频域响应 [h,w] = freqs(b,a,f) b,a 定义同上，f 频率点个数 h 频域响应，w 频域变量)1()2()1()1()2()1()()()(11++++++++==--m a s a s a n b s b s b s A s B s H m m n n三．实验内容1 周期信号傅里叶级数 已知连续时间信号()()2/8cos 3/4coscos )(321ππ++++=t A t A wt A t x ，其中321,,A A A 取值如下：（X 为学号的后两位）]10,1[,5.02321∈⎪⎩⎪⎨⎧===X X A X A X A ]20,11[,55321∈⎪⎩⎪⎨⎧+==-=X X A XA X A ⎪⎩⎪⎨⎧=-=-=X A X A X A 32151020,>X 要求画出信号的时域波形和频域波形（幅度谱和相位谱）。
信号与系统实验2：信号与卷积
Matlab 的函数 conv 函数计算卷积积分 x(t) h(t) ，并画出图形。
(2) 画出函数 f1(t) eatu(t) 和 f2 (t) sin tu(t) 的图形，并利用 Matlab 的函数 conv 函数计算卷积 积分 f1(t) f2 (t) ，并画出图形。
(3) 画出教材 P131 例 3-45 中 x[k] 1, 2,3, 4;k 0,1, 2,3, y[k] 1,1,1,1,1;k 0,1, 2,3, 4 的图形，
4.3 step 函数： 计算并画出系统阶跃响应曲线 调用格式：该函数与函数 impulse() 一样，也有相似的调用格式。
三、实验内容及步骤 对书中的例题进行仿真：
1、连续时间系统的冲激响应、阶跃响应
(1)
利用 impulse 函数画出教材 P127 例 3-42:
LTI
系统
d
2
yzs dt
(t
)
stem(f1,'fill'); title('f1=[1 2 3 4];'); xlabel('t'); ylabel('f1(t)'); subplot(1,3,2); stem(f2,'fill') title('f2=[1 1 1 1]'); xlabel('t'); ylabel('f2(t)'); subplot(1,3,3); stem(f,'fill') title('卷积结果'); xlabel('t(sec)'); ylabel('f1(t)');
调用格式： impulse(b,a) 该调用格式以默认方式绘出向量 a 和 b 定义的连续系统的冲激响应的时域波形。
信号与系统实验实验2常用离散时间信号的实现_OK
信号与系统实验实验2常用离散时间信号的实现_OK实验2:常用离散时间信号的实现实验目的：1.掌握离散时间信号的生成方法；2.学习利用MATLAB进行常用离散时间信号的绘制；3.了解常用离散时间信号的特点。
实验原理：离散时间信号是在离散时间点上取值的信号。
常用离散时间信号包括单位冲激函数(δ[n]), 阶跃函数(u[n]), 单位斜坡函数(ramp[n]), 正弦函数(sin[n]), 余弦函数(cos[n])等。
其中，单位冲激函数是一类特殊的信号，它在n=0时取值为1，其余时刻取值为0。
阶跃函数是另一类重要的信号，它在n>=0时取值为1，n<0时取值为0。
单位斜坡函数是一个连续递增的线性函数。
正弦函数和余弦函数是一类周期性信号，它们的周期为2π。
实验步骤：1.打开MATLAB软件，并新建一个空白的脚本文件；2.定义一个离散时间变量n，并确定它的取值范围；3.根据离散时间变量n生成不同的离散时间信号；4.利用MATLAB的绘图函数绘制离散时间信号的图像；5.对比、分析绘制的不同离散时间信号的特点。
以下是对常用离散时间信号的实现代码和绘制方法的介绍：1.单位冲激函数(δ[n]):单位冲激函数在n=0时取值为1，其余时刻取值为0。
可以使用MATLAB的脉冲函数`dirac(n)`来生成单位冲激函数。
绘制单位冲激函数的图像可以使用MATLAB的`stem(n, dirac(n))`函数。
示例代码：```MATLABn=-10:10;%离散时间变量n的取值范围x = dirac(n); % 生成单位冲激函数stem(n, x); % 绘制单位冲激函数的图像title('单位冲激函数');xlabel('离散时间变量n');ylabel('幅值');```2.阶跃函数(u[n]):阶跃函数在n>=0时取值为1，n<0时取值为0。
可以使用MATLAB的阶跃函数`heaviside(n)`来生成阶跃函数。
信号与系统的实验报告（2）
信号与系统实验报告——连续时间系统的复频域分析班级：05911101学号：**********姓名：***实验五连续时间系统的复频域分析——1120111487 信息工程（实验班）蒋志科一、实验目的①掌握拉普拉斯变换及其反变换的定义，并掌握MA TLAB 实现方法 ②学习和掌握连续时间系统系统函数的定义及其复频域分析方法③掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。
二、实验原理与方法 1、拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为：X s =x (t )e −st dt +∞−∞拉普拉斯反变换为：x t =12πj X (s )e st ds σ+j ∞σ−j ∞在MA TLAB 中可以采用符号数学工具箱中的laplace 函数和ilaplace 函数进行拉氏变换和拉氏反变换。
L=laplace(F)符号表达式F 的拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。
L=laplace(F,t)用t 替换结果中的变量s 。
F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换，返回时间变量t 的结果表达式。
F=ilaplace(L,x)用x 替换结果中的变量t 。
2、连续时间系统的系统函数连续时间系统的系统函数是系统单位冲激响应的拉氏变换H s =ℎ(t )e −st dt +∞−∞此外，连续时间系统的系统函数还可以由系统输入和输出信号的拉氏变换之比得到H s =Y(s)/X(s) 单位冲激响应h(t)反映了系统的固有性质，而H(s)从复频域反映了系统的固有性质。
对于H(s)描述的连续时间系统，其系统函数s 的有理函数H s =b M s M +b M−1s M−1+⋯+b 0a n s n +a n −1s M−1+⋯+a 03、连续时间系统的零极点分析系统的零点指使式H s 的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统函数的值无穷大。
信号与系统实验(MATLAB 西电版)实验2 常用离散时间
实验2 常用离散时间信号的实现 图 2.3 数值法生成的单位阶跃序列
实验2 常用离散时间信号的实现
4) MATLAB
clf; c=-(1/12)+(pi/6)*i; K=2; n=0：40; x=K*exp(c*n); subplot(2，1，1);
实验2 常用离散时间信号的实现
stem(n，real(x)); ylabel(′幅值f(k)′); title(′实部′); subplot(2，1，2); stem(n，imag(x)); xlabel(′时间(k)′); ylabel(′幅值f(k)′); title(′虚部′); 用数值法生成的复指数序列如图2.4
实验2 常用离散时间信号的实现 图 2.6 数值法生成的正弦序列
实验2 常用离散时间信号的实现
7) 单位斜坡序列 MATLAB
clf; k1=-1; k2=20; k0=0; n=［k1：k2］; if k0>=k2 x=zeros(1，length(n)); elseif (k0<k2)&(k0>k1)
实验2 常用离散时间信号的实现 图 2.10 数值法生成的幅值调制序列
实验2 常用离散时间信号的实现
11) MATLAB clf; R=51; d=0.8*(rand(1，R)-0.5); % m=0：R-1; s=2*m.*(0.9.^m); % x=s+d; %
实验2 常用离散时间信号的实现
实验2 常用离散时间信号的实现 图 2.8 数值法生成的随机序列
实验2 常用离散时间信号的实现
9) 扫频正弦序列 MATLAB
n=0：100; a=pi/2/100; b=0; arg=a*n.*n + b*n; x=cos(arg); clf; stem(n，x); axis(［0，100，-1.5，1.5］); grid; axis; title(′扫频正弦序列′) xlabel(′k′); ylabel(′f(k)′); 用数值法生成的扫频正弦序列如图2.9
离散信号与系统的频谱分析实验报告
实验二 离散信号与系统的频谱分析一、实验目的1.掌握离散傅里叶变换（DFT ）及快速傅里叶变换（FFT ）的计算机实现方法。
2.检验序列DFT 的性质。
3.掌握利用DFT （FFT ）计算序列线性卷积的方法。
4.学习用DFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差，以便在实际中正确应用DFT 。
5.了解采样频率对谱分析的影响。
6.了解利用FFT 进行语音信号分析的方法。
二、实验设备1.计算机2.Matlab 软件7.0以上版本。
三、实验内容1.对不同序列进行离散傅里叶变换并进行分析；DFT 共轭对称性质的应用（通过1次N 点FFT 计算2个N 点实序列的DFT ）。
2.线性卷积及循环卷积的关系，以及利用DFT （FFT ）进行线性卷积的方法。
3.比较计算序列的DFT 和FFT 的运算时间。
4.利用FFT 实现带噪信号检测。
5.利用FFT 计算信号频谱及功率谱。
6.扩展部分主要是关于离散系统采样频率、时域持续时间、谱分辨率等参数之间的关系，频谱的内插恢复，对语音信号进行简单分析。
四、实验原理1.序列的离散傅里叶变换及性质离散傅里叶变换的定义：10, )()]([)(102-≤≤==∑-=-N k en x n x DFT k X N n nk Nj π离散傅里叶变换的性质：（1）DFT 的共轭对称性。
若)()()(n x n x n x op ep +=，[])()(n x DFT k X =，则：)()]([k X n x DFT R ep =， )()]([k jX n x DFT I op =。
（2）实序列DFT 的性质。
若)(n x 为实序列，则其离散傅里叶变换)(k X 为共轭对称，即10),()(*-≤≤-=N k k N X k X 。
（3）实偶序列DFT 的性质。
若)(n x 为实偶序列，则其离散傅里叶变换)(k X 为实偶对称，即10),()(-≤≤-=N k k N X k X 。
信号与系统实验一、二
chapter1实验内容：1、画出以下连续时间信号的波形1-0)f(t)=cos(2πt)代码如下：pi=3.14159;t=0:0.01:8;fa=cos(2*pi*t);plot(t,fa);1-1)f (t)=sin(2πt)代码如下：pi=3.14159;t=0:0.01:8;fa=sin(2*t*pi); plot(t,fa);2-0)f (t)=Sa(t/π) 代码如下：pi=3.14159;t=0:0.01:8;fa=sinc(t/pi); plot(t,fa);3-0)f (t)=2[u(t 3)- u(t 5)] 代码如下：t=-1:0.01:10;ft=2*((t>=3)-(t>=5)); plot(t,ft);axis([-1,10,0,3]);4-1)f (t)=e t 代码如下：t=0:0.01:10; ft=exp(t);plot(t,ft);4-2)f (t)=e-t u(t) 代码如下：t=0:0.01:10;f1=(t>=0);f2=exp(-t);plot(t,f1.*f2);5-0）f(t)=2e j(π/4)t，画出实部、虚部、模和相角的波形代码如下：t=0:0.01:10;ft=2*exp(j*(pi/4)*t);h=real(ft); %实部g=imag(ft); %虚部r=abs(ft); %模a=angle(ft); %相角subplot(2,2,1),plot(t,h),title('实部') subplot(2,2,3),plot(t,g),title('虚部') subplot(2,2,2),plot(t,r),title('模')subplot(2,2,4),plot(t,a),title('相角')7）f (t) = u(t)代码如下：t=-1:0.01:5ft=(t>=0);plot(t,ft);axis([-1,5,0,1.5]);8）f (t) =δ(t)代码如下：t=-1:0.01:5;ft=(t>=0)-(t>=0.1); plot(t,ft);axis([-1,1,0,1.1]);9)f9为周期矩形信号，其幅度从-1 到1，占空比为75% 代码如下：pi=3.14159;t=-10:0.01/pi:10;ft=square(t,75);plot(t,ft);2、信号本身运算画出f1(t)为宽度是4，高为1，斜度为0.5 的三角脉冲，然后画出f1(-t)，f1(2t)，f1(2-2t)的波形以及f1(t)的微分和积分波形。
- 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);。