北京理工大学信号与系统实验报告2 LTI系统的时域分析
北京理工大学信号和系统实验报告材料
本科实验报告实验名称：信号与系统实验实验一信号的时域描述与运算一、实验目的①掌握信号的MATLAB表示及其可视化方法。
②掌握信号基本时域运算的MATLAB实现方法。
③利用MATLAB分析常用信号，加深对信号时域特性的理解。
二、实验原理与方法1. 连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。
在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。
从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。
表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。
例如一个正弦信号可以表示如下：>> t=0:0.01:10;>> x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。
如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。
例如对于上述正弦信号，可以用符号对象表示如下：>> x=sin(t);>> ezplot(X);利用ezplot(x)命令可以绘制上述信号的时域波形常用的信号产生函数2.连续时间信号的时域运算-1-0.8-0.6-0.4-0.200.20.40.60.81Time(seconds)图1 利用向量表示连续时间信号-1-0.50.51t图 2 利用符号对象表示连续时间信号sin(t)对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。
1）相加和相乘信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。
实验2离散时间LTI系统的时域分析
实验二 离散时间LTI 系统的时域分析一 实验目的（1） 学会运用MATLAB 求解离散时间系统的零状态响应；（2） 学会运用MATLAB 求解离散时间系统的单位取样响应；（3） 学会运用MATLAB 求解离散时间系统的卷积和。
二 实验原理及实例分析1、离散时间系统的响应离散时间LTI 系统可用线性常系数差分方程来描述，即∑∑==-=-Mj jN i i j n x b i n y a 00)()( （1） 其中，i a （0=i ，1，…，N ）和j b （0=j ，1，…，M ）为实常数。
MATLAB 中函数filter 可对式（1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter 的语句格式为y = filter (b , a , x )其中，x 为输入的离散序列；y 为输出的离散序列；y 的长度与x 的长度一样；b 与a 分别为差分方程右端与左端的系数向量。
【实例1】 已知某LTI 系统的差分方程为)1(2)()2(2)1(4)(3-+=-+--n x n x n y n y n y试用MATLAB 命令绘出当激励信号为)()2/1()(n u n x n=时，该系统的零状态响应。
解：MATLAB 源程序为>>a=[3 -4 0 2];>>b=[1 2]; >>n=0:30;>>x=(1/2).^n;>>y=filter(b,a,x);>>stem(n,y,'fill'),grid on>>xlabel('n'),title('系统响应y(n)')程序运行结果如图1所示。
2、离散时间系统的单位取样响应系统的单位取样响应定义为系统在)(n δ激励下系统的零状态响应，用)(n h 表示。
MATLAB 求解单位取样响应可利用函数filter ，并将激励设为前面所定义的impDT 函数。
实验二 离散时间LTI系统的时域分析
实验二离散时间LTI系统的时域分析1．实验目的通过本实验，要求学生学会运用MATLAB求解离散时间系统的零状态响应；学会运用MATLAB求解离散时间系统的单位冲激响应；学会运用MATLAB求解离散时间系统的卷积和。
2．实验原理MATLAB中函数filter可对上式的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter的语句格式为 y=filter(b,a,x) 其中，x为输入的离散序列；y 为输出的离散序列；y的长度与x的长度一样；b与a分别为差分方程右端和左端的系数向量。
系统的单位冲激响应定义为系统在δ(n)激励下系统的零状态响应，用h(n）表示。
MATLAB求解单位冲激响应可利用函数filter，并设激励为δ(n)函数。
系统的单位阶跃响应定义为系统在u (n)激励下系统的零状态响应，用g(n）表示。
MATLAB 求解单位阶跃响应可利用函数filter，并设激励为u(n)函数。
系统的零状态响应是激励与系统的单位冲激响应的卷积。
离散时间信号的卷积运算是求和运算，因而常称为卷积和。
MATLAB求离散时间信号卷集和的命令为conv，其语句格式为 y=conv(x,h)。
例如，利用MATLAB的conv命令求两个长为4的矩形序列的卷积和，即g(n)=[u(n)- u(n-4)]*[ u(n)- u(n-4)]，其结果应是长为7(4+4-1=7）的三角序列。
运行结果如图10-7所示。
图10-7 矩形序列3．实验内容(1) 已知x(n)=u(n-“学号后两位”)-u(n-(“学号后两位”+4))、h(n)=R5(n)。
(2) 编制程序求解下列两个系统的单位冲激响应和单位阶跃响应，并绘出其图形。
(3) 已知某系统的单位冲激响应为h(n)=0.5n［u(n)-u(n-4)］，试分别利用MATLAB卷积和两种方法求当激励信号为x(n)=u(n)-u(n-4)时，系统的零状态响应。
4．实验设备计算机，MATLAB软件。
LTI系统的时域分析
实验二：LTI系统的时域分析一实验目的1．掌握利用计算机进行连续时间信号卷积运算和离散时间信号卷积运算的原理和方法；2．掌握用MATLAB求解LTI系统的时域响应。
二、涉及的MATLAB函数1、conv（）函数2、filter（）函数3、impulse（）函数4、dsolve（）函数三、实验内容与方法1、信号的卷积运算p=0.001; % 取样时间间隔nf=0:p:2; % f(t)对应的时间向量f=2*((nf>=0)-(nf>=2)); % 序列f(n)的值nh=0:p:5; % h(t)对应的时间向量h=exp(-nh); % 序列h(n)的值y=conv(f,h); % 计算序列f(n)与h(n)的卷积和y(n)y=y*p; % y(n)变成y(t)left=nf(1)+nh(1) % 计算序列y(n)非零样值的起点位置right=length(nf)+length(nh)-2 % 计算序列y(n)非零样值的终点位置k=p*(left:right);subplot(3,1,1),stairs(nf,f); % 绘制f(t)的波形title('f(t)');axis([0 3 0 2.1]);subplot(3,1,2),stairs(nh,h); % 绘制h(t)的波形title('h(t)');axis([0 3 0 2.1]);subplot(3,1,3),plot(k,y); % 绘制y(t)=f(t)*h(t)的波形title('y(t)=f(t)*h(t)');axis([0 3 0 3]);2、系统的时域分析f(k)=(k),b=[1 1];a=[1 -0.25 0.5]; % 差分方程的系数n=0:15; % 序列的个数fn=0.5.^n; % 输入序列y1=filter(b,a,fn); % 零状态响应y2=impz(b,a,16); % 单位响应subplot(1,2,1),stem(n,y1,'filled');title('零状态响应');grid onsubplot(1,2,2),stem(n,y2,'filled');title('单位响应');grid on五、实验报告要求1、实验目的概述；2、实验原理分析；3、将源代码及其运行结果附在报告中并加以分析说明；4、实验收获及体会；。
北京理工大学信号与系统实验实验报告材料
实验1 信号的时域描述与运算一、实验目的1. 掌握信号的MATLAB 表示及其可视化方法。
2. 掌握信号基本时域运算的MATLAB 实现方法。
3. 利用MATLAB 分析常用信号，加深对信号时域特性的理解。
二、实验原理与方法1. 连续时间信号的MATLAB 表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。
在MATLAB 中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。
从严格意义上来说，MATLAB 并不能处理连续时间信号，在MATLAB 中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。
表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。
例如一个正弦信号可以表示如下：>> t=0:0.01:10; >> x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。
如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。
例如对于上述正弦信号，可以用符号对象表示如下：>> x=sin(t); >> ezplot(X);利用ezplot(x)命令可以绘制上述信号的时域波形012345678910-1-0.8-0.6-0.4-0.200.20.40.60.81Time(seconds)图1 利用向量表示连续时间信号-6-4-20246-1-0.50.51t图 2 利用符号对象表示连续时间信号sin(t)常用的信号产生函数 函数名 功能 函数名 功能 heaviside 单位阶跃函数 rectpuls 门函数 sin 正弦函数 tripuls 三角脉冲函数 cos 余弦函数 square 周期方波sinc sinc 函数 sawtooth 周期锯齿波或三角波 exp指数函数2.连续时间信号的时域运算对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。
信号与系统实验二 连续LTI系统的时域研究分析
信号与系统实验二连续LTI 系统的时域分析————————————————————————————————作者：————————————————————————————————日期：2实验二 连续LTI 系统的时域分析一. 实验目的1. 加深对线性时不变系统中零状态响应概念的理解，掌握其求解方法；2. 掌握求解给定连续系统的冲激响应和阶跃响应的方法。
二. 实验原理1.连续系统零状态响应的数值解线性时不变 (LTI) 连续时间系统用常系数线性微分方程进行描述，系统的零状态响应就是在系统初始状态为零条件下微分方程的解。
MATLAB 控制系统工具箱提供了一个lism 函数来求解连续时间系统的零状态响应，其调用格式为y = lism(sys,f,t)其中t 表示计算系统响应的时间抽样点向量，f 是系统输入信号向量，sys 是LTI 系统模型，用来表示微分方程、状态方程。
在求解微分方程时，微分方程的LTI 系统模型sys 要借助MATLAB 中的tf 函数来获得，其调用格式为sys = tf(b,a)其中a 、b 分别为微分方程左端和右端各项的系统向量。
例如系统方程 (3)(2)(1)(2)(1)2210210()()()()()()()a y t a y t a y t a y t b f t b f t b f t +++=++该方程左边、右边的系数向量分别为3210[,,,]a a a a a =，210[,,]b b b b =。
例1：描述某线性时不变系统的方程为"()4'()4()'()3()y t y t y t f t f t ++=+试求：当()()tf t e t ε-=时，系统的零状态响应()zs y t 。
解：实现所要求运算的m 文件如下，a = [1 4 4]; %将y （t ）各阶导数的系数放在向量a 中b = [1 3]; %将f （t ）各阶导数的系数放在向量b 中sys = tf(b, a); %求系统模型systd = 0.01; %定义时间间隔t = 0 : td : 10; %定义时间向量f = exp(-t); %将f （t ）表示出来y = lsim(sys, f, t); %求系统的零状态响应plot(t, y); %绘出零状态响应的波形xlabel('t(sec)'); %给出x 坐标的标签ylabel('y(t)'); %给出y 坐标的标签grid on %在图上显示方格程序运行结果见图1。
北京理工大学信号与系统实验报告
实验1 信号的时域描述与运算一、实验目的1、掌握信号的MATLAB表示及其可视化方法。
2、掌握信号基本时域运算的MATLAB实现方法。
3、利用MATLAB分析常用信号，加深对信号时域的理解。
二、实验原理1、连续时间的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除若干个不连续点外，在任何信号都有意义。
在MATLAB中，连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。
向量表示法：严格意义上来说，MATLAB并不能处理连续时间信号，都必须是用信号等时间间隔采样后的采样值来近似表示的，采样时间间隔足够小的时候，这些采样值就可以近似地表示出连续时间信号。
例如：>>t=0:0.01:10;>>x=sin(t);此时利用plot(t,x)命令即可绘制上述信号的时域波形。
符号对象表示法：连续时间信号先用表达式表示出来，然后采用符号表达式来表示信号。
例如：>>sym t;>>x=xin(t);此时利用ezplot(x)命令即可绘制上述信号的时域波形。
常用的信号产生函数：2、连续时间信号的时域运算对连续时间信号的运算包括量信号想家、相乘、微分、积分以及位移反转、尺度变换（尺度伸缩）等1）相加和相乘信号的相加和相乘指两个信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“•”来计算，此时要求表示两信号的向量时间范围和采样间隔相同，采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。
2）微分和积分对于向量表示发表示的连续时间信号，可以用过数值计算的方法计算信号的微分和积分。
这里由时间向量[t1,t2,…,t N]和采样值向量[x1,x2,…,x N]表示的连续信号的微分是利用差分来近似求取的。
MATLAB里用diff来计算差分x(k+1)-x(k)。
连续信号的定积分可以由MATLAB的quad函数实现，调用格式为quad（‘functions_name’,a,b）其中，functions_name为被积函数名，a、b为积分区间。
《信号与系统》第二章 北京理工大学
2.1 引言
时域分析方法：
输入信号（激励）
时域数学模型 输出信号（响应） （输入输出法）
LTI系统零输入和零状态的求取： ➢解微分方程法 ➢卷积积分法
2.2 LTI系统的微分方程表示及其响应
2 微分方程的求解 ▪微分方程的完全解：由齐次解 yh (t)和特解 yp (t)两部分 组成。 ▪齐次解的求法 表2－1 ▪特解的求法 表2－2 ▪完全解:其中齐次解中的待定系数应在完全解中由给 定的附加初始条件确定。
将x(t ) e t 代入方程的右端，得 e t ,其中a 1与1＝－1相重，
故特解 y p (t ) A1te t
代入微分方程，得 A1＝1，因此 y p (t ) te t
（3）完全解
y(t ) yh (t ) y p (t ) c1e t c2e 2t te t ,由初始条件
P47，例2－8
2.4 卷积积分
1 用冲激函数表示任意信号
u(t) 0 (t )d
2 卷积积分
x(t) x( ) (t )d
y(t) x( )h(t )d y(t) x(t) h(t)
卷积积分公式
LTI系统对于任意信号 x(t) 的零状态响应，可以由 该信号与系统的单位冲激响应的卷积积分得到。
y(0) 1, y(0) 3 得 y(0) c1 c2 0
解得c1 2,
y(0) c1 2c2 1 3
c2 2
y(t ) 2e t 2e 2t te t , t 0
2.3零输入响应和零状态响应
➢零输入响应（ y0(t)）：由微分方程的齐次解得到，齐次解中 的系数由给定的初始条件确定； ➢零状态响应 ( yx (t) ) ：由方程的全解得到，其中齐次解的系 数应在全解中由初始条件确定（y(0)＝y(1)(0)=…= y(n-1)(0)= 0。 ？完全解和完全响应 完全解＝齐次解＋特解 完全响应＝零输入响应＋零状态响应
实验2-LTI系统地时域分析报告
一，实验目的作为根底性实验局部，实验2使我们了解和掌握利用MATLAB对系统进展时域分析的方法，掌握了连续时间系统和离散时间系统下对零状态响应、单位抽样响应的方法，以与求卷积积分和卷积和的方法。
二，实验原理〔1〕连续时间系统时域分析的MATLAB实现。
①连续时间系统的MATLAB表示。
用系统微分方程描述LTI连续系统，然后在matlab中建立模型：b=[b1,b2,……]a=[a1,a2,……]sys=tf(b,a)②连续时间系统的零状态响应。
调用函数lsim(sys,x,t)绘出信号与响应的波形。
③连续时间系统的冲击响应与阶跃响应。
描述系统的单位冲击响应调用impulse函数：impulse(sys)在默认时间X围内绘出系统冲激响应的时域波形。
impulse(sys，T)绘出系统在0~TX围内冲激响应的时域波形。
impulse(sys，ts:tp:te)绘出系统在ts~teX围内，以tp为时间间隔取样的冲击响应波形。
描述系统的单位阶跃响应调用step函数：impulse(sys)impulse(sys，T)impulse(sys，ts:tp:te)〔2〕离散时间系统时域分析的MATLAB实现。
①离散时间系统的MATLAB表示。
用向量b=[b1,b2,……]，a=[a1,a2,……]可以表示系统。
②离散时间系统对任意输入的响应。
可以调用函数filter〔b,a,x〕③离散时间系统的单位抽样响应。
可以调用函数impz:impz〔b,a〕在默认时间X围内绘出系统单位抽样响应的时域波形。
impz (b,a,N绘出系统在0~NX围内单位抽样响应的时域波形。
impz(b,a,ns:ne)绘出系统在ns~neX围内的单位抽样响应波形。
〔3〕卷积与卷积积分①离散时间序列的卷积和可以调用函数conv求得两个离散序列的卷积和。
②连续时间信号的卷积积分在取样间隔足够小的情况下，由卷积和近似求得卷积积分。
三，实验内容〔1〕描述模拟低通、高通、带通和带阻滤波器的微分方程如下，试采用MATLAB绘出各系统的单位冲激响应和单位阶跃响应波形。
北京理工大学信号与系统实验实验报告
北京理工大学信号与系统实验实验报告信号与系统实验报告姓名:肖枫学号:1120111431班号:05611102专业:信息对抗技术学院:信息与电子学院12实验1 信号的时域描述与运算一、实验目的1. 掌握信号的MATLAB表示及其可视化方法。
2. 掌握信号基本时域运算的MATLAB实现方法。
3. 利用MATLAB分析常用信号，加深对信号时域特性的理解。
二、实验原理与方法1. 连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。
在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。
从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。
表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。
例如一个正弦信号可以表示如下:>> t=0:0.01:10;>> x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。
如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。
例如对于上述正弦信号，可以用符号对象表示如下:>> x=sin(t);>> ezplot(X);利用ezplot(x)命令可以绘制上述信号的时域波形10.80.60.40.2-0.2-0.4-0.6-0.8-1012345678910Time(seconds)图1 利用向量表示连续时间信号3sin(t)10.5-0.5-1-6-4-20246t图 2 利用符号对象表示连续时间信号常用的信号产生函数函数名功能函数名功能 heaviside 单位阶跃函数 rectpuls 门函数 sin 正弦函数 tripuls 三角脉冲函数 cos 余弦函数 square 周期方波 sinc sinc函数 sawtooth 周期锯齿波或三角波 exp 指数函数2.连续时间信号的时域运算对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换(尺度伸缩)等。
实验2-LTI系统的时域分析报告
一，实验目的作为基础性实验部分，实验2使我们了解和掌握利用MATLAB对系统进行时域分析的方法，掌握了连续时间系统和离散时间系统下对零状态响应、单位抽样响应的方法，以及求卷积积分和卷积和的方法。
二，实验原理（1）连续时间系统时域分析的MATLAB实现。
①连续时间系统的MATLAB表示。
用系统微分方程描述LTI连续系统，然后在matlab中建立模型：b=[b1,b2,……]a=[a1,a2,……]sys=tf(b,a)②连续时间系统的零状态响应。
调用函数lsim(sys,x,t)绘出信号及响应的波形。
③连续时间系统的冲击响应与阶跃响应。
描述系统的单位冲击响应调用impulse函数：impulse(sys)在默认时间围绘出系统冲激响应的时域波形。
impulse(sys，T)绘出系统在0~T围冲激响应的时域波形。
impulse(sys，ts:tp:te)绘出系统在ts~te围，以tp为时间间隔取样的冲击响应波形。
描述系统的单位阶跃响应调用step函数：impulse(sys)impulse(sys，T)impulse(sys，ts:tp:te)（2）离散时间系统时域分析的MATLAB实现。
①离散时间系统的MATLAB表示。
用向量b=[b1,b2,……]，a=[a1,a2,……]可以表示系统。
②离散时间系统对任意输入的响应。
可以调用函数filter（b,a,x）③离散时间系统的单位抽样响应。
可以调用函数impz:impz（b,a）在默认时间围绘出系统单位抽样响应的时域波形。
impz(b,a,N绘出系统在0~N围单位抽样响应的时域波形。
impz(b,a,ns:ne)绘出系统在ns~ne围的单位抽样响应波形。
（3）卷积与卷积积分①离散时间序列的卷积和可以调用函数conv求得两个离散序列的卷积和。
②连续时间信号的卷积积分在取样间隔足够小的情况下，由卷积和近似求得卷积积分。
三，实验容（1）已知描述模拟低通、高通、带通和带阻滤波器的微分方程如下，试采用MATLAB绘出各系统的单位冲激响应和单位阶跃响应波形。
北京理工大学信号与系统实验报告2 LTI系统的时域分析
实验2 LTI 系统的时域分析（基础型实验）一．实验目的1. 掌握利用MATLAB 对系统进行时域分析的方法。
2. 掌握连续时间系统零状态响应、冲击响应和阶跃响应的求解方法。
3. 掌握求解离散时间系统响应、单位抽样响应的方法。
4.加深对卷积积分和卷积和的理解。
掌握利用计算机进行卷积积分和卷积和计算的方法。
二． 实验原理与方法1. 连续时间系统时域分析的MATLAB 实现1） 连续时间系统的MA TLAB 表示LTI 连续系统通常可以由系统微分方程描述，设描述系统的微分方程为：(N)(N 1)(M)(M 1)1010(t)(t)...(t)b (t)b (t)...b (t)N N M M a y a y a y x x x ----++=++则在MATLAB 中可以建立系统模型如下：1010[b ,b ,...,b ];a [a ,a ,...,a ];sys tf(b,a);M M N N b --=== 其中，tf 是用于创建系统模型的函数，向量a 和b 的元素是以微分方程求导的降幂次序来排列的，如果有缺项，应用0补齐，例如由微分方程2''(t)y'(t)3y(t)x(t)y ++=描述的系统可以表示为： >> b=[1]; >> a=[2 1 3]; >> sys=tf(b,a); 而微分方程由''(t)y'(t)y(t)x''(t)x(t)y ++=-描述的系统则要表示成 >> b=[1 0 -1]; >> a=[1 1 1]; >> sys=tf(b,a);2） 连续时间系统的零状态响应零状态响应指系统的初始状态为零，仅由初始信号所引起的响应。
MATLAB 提供了一个用于求解零状态响应的函数lism ，其调用格式如下：lism （sys,x,t ）绘出输入信号及响应的波形，x 和t 表示输入信号数值向量及其时间向量。
实验二-LTI系统时域分析(guo)
实验二——LTI 系统的时域分析实验性质：提高性 实验级别：必做开课单位：机械电子工程学院 学 时：2一、实验目的1、深刻理解卷积运算，利用离散卷积实现连续卷积运算；2、深刻理解信号与系统的关系，学习MA TLAB 语言实现信号通过系统的仿真方法。
二、实验设备计算机，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 )(00t 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 )(00⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==τττd t h f k t h k f t y t y k )()()()(lim )(lim )(00所以：∆∆-∆=-==∑⎰→∆)()(lim )()()(*)()(2102121k 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 （）求连续卷积，只需∆足够小以及在卷积和的基础上乘以∆。
北京理工大学信号与系统实验报告
实验1 信号的时域描述与运算一、实验目的1、掌握信号的MATLAB表示及其可视化方法。
2、掌握信号基本时域运算的MATLAB实现方法。
3、利用MATLAB分析常用信号，加深对信号时域的理解。
二、实验原理1、连续时间的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除若干个不连续点外，在任何信号都有意义。
在MATLAB中，连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。
向量表示法：严格意义上来说，MATLAB并不能处理连续时间信号，都必须是用信号等时间间隔采样后的采样值来近似表示的，采样时间间隔足够小的时候，这些采样值就可以近似地表示出连续时间信号。
例如：>>t=0:0.01:10;>>x=sin(t);此时利用plot(t,x)命令即可绘制上述信号的时域波形。
符号对象表示法：连续时间信号先用表达式表示出来，然后采用符号表达式来表示信号。
例如：>>sym t;>>x=xin(t);此时利用ezplot(x)命令即可绘制上述信号的时域波形。
常用的信号产生函数：2、连续时间信号的时域运算对连续时间信号的运算包括量信号想家、相乘、微分、积分以及位移反转、尺度变换（尺度伸缩）等1）相加和相乘信号的相加和相乘指两个信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“•”来计算，此时要求表示两信号的向量时间范围和采样间隔相同，采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。
2）微分和积分对于向量表示发表示的连续时间信号，可以用过数值计算的方法计算信号的微分和积分。
这里由时间向量[t1,t2,…,t N]和采样值向量[x1,x2,…,x N]表示的连续信号的微分是利用差分来近似求取的。
MATLAB里用diff来计算差分x(k+1)-x(k)。
连续信号的定积分可以由MATLAB的quad函数实现，调用格式为quad（‘functions_name’,a,b）其中，functions_name为被积函数名，a、b为积分区间。
北京理工大学信号与系统实验报告
实验一信号的时域描述与运算一、实验目的①掌握信号的MATLAB表示及其可视化方法。
②掌握信号基本时域运算的MATLAB实现方法。
③利用MATLAB分析常用信号，加深对信号时域特性的理解。
二、实验原理与方法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 利用向量表示连续时间信号常用的信号产生函数sin(t)t图 2 利用符号对象表示连续时间信号2.连续时间信号的时域运算对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。
1）相加和相乘信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。
采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。
信号与系统第2章 信号通过LTI系统的时域分析
因此，f(t)的第n个分段可近似表示为
f n (t ) f (tn )[ (t tn ) (t tn )]（2-3）
图2-1
使用矩形脉冲逼近f(t)
而f(t)就可近似表示为这个分段之和 ，即
f (t )
fn (t ) f (tn )[ (t tn ) (t tn )] n 0 n 0 （2-4） N 1 (t tn ) (t tn ) f (tn )
对式（2-8）中的积分变量作变量置换， d dt1 ，得到 令 t t1 ，因此 t t1 ，
y(t )
∞ ∞
x(t t1 )h(t1 )dt1
∞ ∞
x(t )h( )d h(t ) x(t ) （2-9）
比较式（2-8）、式（2-9）可知，卷 积服从交换律。 这个分解表达式及其物理意义
首先考察下面的数学表达式
∞ ∞
f ( )δ(t )d f (t )
（2-1）
表达式（2-1）在前面1.3.2小节介绍 (t)性质时已经指出，这个表达式的物理 意义是指任何一个连续时间信号可以分 解为单位冲激信号的线性组合。 下面对此进行展开说明。
n 0
N 1
N 1
2.3 信号通过LTI系统的时域分析与卷积积分
2.3.1 分析
如图2-2所示，假设LTI系统处于初始 松弛状态，输入信号为x(t)，则利用LTI系 统的线性和时不变性，输出信号为
y (t ) T [ x(t )] T[
∞ ∞ ∞ ∞ ∞ ∞ ∞ ∞
初始松弛时，LTI系统输出y(t)是输 入x(t)与系统单位冲激响应h(t)的卷积这 一结果表明，对于LTI系统，h(t)已经给 出了系统的全部信息，也即表征了系统 的全部性质。 因此，LTI系统现已可用图2-4所示 的框图来表示。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验2 LTI 系统的时域分析（基础型实验）一．实验目的1. 掌握利用MATLAB 对系统进行时域分析的方法。
2. 掌握连续时间系统零状态响应、冲击响应和阶跃响应的求解方法。
3. 掌握求解离散时间系统响应、单位抽样响应的方法。
4.加深对卷积积分和卷积和的理解。
掌握利用计算机进行卷积积分和卷积和计算的方法。
二． 实验原理与方法1. 连续时间系统时域分析的MATLAB 实现1） 连续时间系统的MA TLAB 表示LTI 连续系统通常可以由系统微分方程描述，设描述系统的微分方程为：(N)(N 1)(M)(M 1)1010(t)(t)...(t)b (t)b (t)...b (t)N N M M a y a y a y x x x ----++=++则在MATLAB 中可以建立系统模型如下：1010[b ,b ,...,b ];a [a ,a ,...,a ];sys tf(b,a);M M N N b --=== 其中，tf 是用于创建系统模型的函数，向量a 和b 的元素是以微分方程求导的降幂次序来排列的，如果有缺项，应用0补齐，例如由微分方程2''(t)y'(t)3y(t)x(t)y ++=描述的系统可以表示为： >> b=[1]; >> a=[2 1 3]; >> sys=tf(b,a); 而微分方程由''(t)y'(t)y(t)x''(t)x(t)y ++=-描述的系统则要表示成 >> b=[1 0 -1]; >> a=[1 1 1]; >> sys=tf(b,a);2） 连续时间系统的零状态响应零状态响应指系统的初始状态为零，仅由初始信号所引起的响应。
MATLAB 提供了一个用于求解零状态响应的函数lism ，其调用格式如下：lism （sys,x,t ）绘出输入信号及响应的波形，x 和t 表示输入信号数值向量及其时间向量。
y= lism （sys,x,t ）这种调用格式不绘出波形，而是返回响应的数值向量。
3） 连续时间系统的冲激响应与阶跃响应MATLAB 提供了函数impulse 来求指定时间范围内，由模型sys 描述的连续时间系统的单位冲激响应impulse 函数基本调用格式如下： impulse(sys)在默认时间范围内绘出系统冲激响应impulse(sys,T)绘出系统在0-T 范围内冲激响应的时域波形。
Impulse(sys,ts:tp:te)绘出系统在ts-te 范围内以tp 为时间间隔取样的冲激响应波形。
[y,t]=impulse(...)该调用格式不绘出冲激响应波形，而是返回冲激响应的数值向量及其对应的时间向量。
求解单位阶跃响应用到函数step ，几种调用格式如下： step(sys) step(sys,T)step(sys,ts:tp:te) [y,t]=step(...)2. 离散时间系统时域分析的MATLAB 实现1） 离散时间系统的MA TLAB 表示LTI 离散系统通常可以由系统差分方程描述，设描述系统的差分方程为：0101(n)(n 1)...(n M)b (n)b (n 1)...b (n N)N M a y a y a y x x x +-++-=+-++-则在MATLAB 里，我们可以用如下两个向量来表示这个系统：0101[,,...,];[,,...,];M N b b b b a a a a ==2） 离散时间系统对任意输入的响应MATLAB 提供了求LTI 离散系统响应的专用函数filter ，该函数用于求取由差分方程所描述的离散时间系统在指定时间范围内对输入序列所产生的响应，该函数的基本调用格式为 y=filter(b,a,x)其中x 为输入序列，y 为输出序列，输出序列y 对应的时间区间与x 对应的时间区间相同。
3） 离散时间系统的单位抽样响应MATLAB 提供了函数impz 来求指定时间范围内，由向量b 和a 描述的离散时间系统的单位抽样响应，具体调用格式如下：impz(b,a)在默认的时间范围内绘出系统抽样响应的时域波形。
impz(b,a,N)绘出系统在0~N 范围内单位抽样响应的时域波形。
impz(b,a,ns:ne)绘出系统在ns~ne 范围内的单位抽样响应波形。
[y,t]=impz(...)该调用格式不绘出单位抽样响应波形，而是返回单位抽样响应的数值向量及其对应的时间向量。
3. 卷积和与卷积积分1） 离散时间序列的卷积和卷积和是离散系统时域分析的基本方法之一，离散时间序列1(n)x 和2(n)x 的卷积和(n)x 定义如下:1212(n)(n)*(n)(k)(n k)k x x x x x +∞=-∞==-∑对已离散LTI 系统，设其输入信号为(n)x ，单位抽样响应为h(n)，则其零状态响应(n)y 为(n)(n)*h(n)y x =即离散LTI 系统的零状态响应可以表示成输入信号x(n)和单位抽样响应h(n)的卷积。
MATLAB 中conv 函数可以用来求两个离散序列的卷积和，调用格式为conv(x1,x2)。
例如： >> x1=ones(1,3); >> x2=[1 2 3 4]; >> x=conv(x1,x2) x =1 3 6 9 7 4此例中x1、x2、x 中都没有时间信息，实际上要考察的信号中还需要知道各时刻对应的时间序列，所以还需要根据序列x1、x2对应的时间序列确定卷积结果x 对应的时间序列。
设x1、x2为两个在有限时间区间内非零的离散时间序列，即序列x1在区间n1~n2内非零，序列x2在区间m1~m2内非零，则序列x1的时域宽度为L1=n2-n1+1，序列x2的时域宽度为L2=m2-m1+1。
由卷积和的定义可知，卷积和序列序列x 的时域宽度为L=L1+L2-1，且只在区间(n1+m1)~(n1+m1)+(L1+L2-2)非零。
2） 连续时间信号的卷积积分对连续LTI 系统，输入x(t)，单位冲击响应为h(t)，零状态响应y(t)，则有(t)(t)*h(t)y x =利用MATLAB 可以采用数值计算的方法近似计算卷积积分，卷积积分可以用求和运算来实现121212(t)(t)*(t)()(t )d lim (k )(t k )k x x x x x x x τττ+∞+∞∆→=-∞-∞==-=∆-∆⋅∆∑⎰现在考虑只求t n =∆时(t)x 的值(n )x ∆，则由上式可得12(n )(k )((n k))k x x x +∞=-∞∆=∆⋅∆-∆∑当∆足够小，(n )x ∆就是()x t 的数值近似。
可以利用计算离散序列卷积和的conv 来计算卷积积分，具体步骤如下：1） 将连续时间信号1()x t 和2()x t 以时间间隔∆进行取样，得到离散序列1(n )x ∆和2(n )x ∆；2） 构造离散序列1(t)x 和2(t)x 对应的时间向量1t 和2t ；3） 调用该函数conv 计算卷积积分在t n =∆时的近似采样值(n )x ∆； 4） 构造离散序列(n )x ∆对应的时间向量n.。
三．实验内容(1)已知描述模拟低通、高通、带通和带阻滤波器的微分方程如下，试采用MA TLAB 绘制出各系统的单位冲激响应和单位阶跃响应的波形。
1) ''(t)'(t)y(t)x(t)y += >> b=[0 0 1]; >> a=[1 sqrt(2) 1] >> sys=tf(b,a); >> subplot(121)>> impulse(sys) >> subplot(122) >> step(sys)Time (seconds)A m p l i t u deTime (seconds)A m p l i t u d e2). ''(t)'(t)y(t)x''(t)y += >> b=[2 0 0];>> a=[1 sqrt(2) 1]; >> sys=tf(b,a); >> subplot(121)>> impulse(sys) >> subplot(122)>>step(sys)Time (seconds)A m p l i t u d eTime (seconds)A m p l i t u d e3) ''(t)'(t)y(t)x'(t)y y ++= >> b=[0 1 0]; >> a=[1 1 1]; >> sys=tf(b,a); >> subplot(121)>> impulse(sys) >> subplot(122) >> step(sys)Time (seconds)A m p l i t u d eTime (seconds)A m p l i t u d e4) ''(t)'(t)y(t)x''(t)x(t)y y ++=+ >> b=[1 1 0]; >> a=[1 1 1]; >> sys=tf(b,a); >> subplot(121)>> impulse(sys) >> subplot(122) >> step(sys)Time (seconds)A m p l i t u deTime (seconds)A m p l i t u d e(2)已知某系统可以由如下微分方程描述''(t)'(t)6y(t)x(t)y y ++= 1)利用MATLAB 绘制出该系统冲激响应和阶跃响应的时域波形。
>> b=[0 0 1]; >> a=[1 1 6]; >> sys=tf(b,a); >> subplot(121)>> impulse(sys) >> subplot(122) >> step(sys)Impulse ResponseTime (seconds)A m p l i t u d eStep ResponseTime (seconds)A m p l i t u d e2)根据冲激响应的时域波分析系统的稳定性。