信号与线性系统实验指导书v02.docx

1、信号与线性系统实验指导书v02信号与线性系统实验指导书东华大学信息学院通信与电子信息工程系实验要求及说明一、 实验报告内容实验报告包括原理分析、源程序、执行结果分析及实验总结，其中原理分析和实验总结需要手写,其他可打印。二、 实验成绩实验成绩包括出勤（10）、实验表现（10）、编程（30）和实验报告（50）几部分。三、 其他说明缺席3次及以上取消考试成绩。目 录实验一 连续信号的时域分析 1实验二 连续时间系统的时域分析 3实验三 连续信号的频域分析 9实验四 连续系统的频域分析 12实验五 信号采样与重建 15实验六 离散时间信号和系统分析 17附 录 MATLAB主要命令函数表 20实验一

2、 连续信号的时域分析一、 实验目的1、熟悉MATLAB软件。2、掌握常用连续信号与离散信号的MATLAB表示方法。二、 实验设备安装有matlab6.5以上版本的PC机一台。三、 MATLAB使用说明1、在MATLAB可视化绘图中，对于以t为自变量的连续信号，在绘图时统一用plot函数。例题：绘出t从-10到10的sin(t)的波形。t=-10:0.05:10;f=sin(t);plot(t,f);title(f(t)=sin(t);xlabel(t) ;axis(-10,10,-1,1)grid on可得图1所示图形。图12、此外也可以利用MATLAB的ezplot函数对连续信号画图。例题：

3、绘出t从-10到10的sin(t)的波形。clcclear all close allsyms t f=sin(t)ezplot(f, -10 10);xlabel(t);title ( f(t)=sin(t) ;grid on图2四、 实验内容1、用MATLAB表示连续信号：，。2、用MATLAB表示抽样信号(sinc(t)、矩形脉冲信号(rectpuls(t, width)及三角脉冲信号(tripuls(t, width, skew)。3、编写如图3的函数并用MATLAB绘出满足下面要求的图形。图3实验二 连续时间系统的时域分析一、实验目的1、掌握卷积计算方法。2、掌握函数lsim，imp

4、ulse，step的用法，lsim为求取零状态响应，impulse为求取单位脉冲响应，step为求取单位阶跃响应。3、运用课堂上学到的理论知识，从RC、RL一阶电路的响应中正确区分零输入响应、零状态响应、自由响应与受迫响应。二、实验设备安装有matlab6.5以上版本的PC机一台。三、MATLAB使用说明1、MATLAB求取系统响应的方法MATLAB的内部函数impulse()，step()，initial()，lsim() 可以用来计算并绘制连续时间LTI系统的单位冲激响应，单位阶跃响应，零输入响应和任意信号作用于系统的零状态响应。在调用这些函数时，需要用向量来对连续系统进行表示，设描述连续

5、系统的微分方程为：，则可以用向量a和b来表示该系统，即： ,例如对微分方程，则用MATLAB表示该系统的对应向量应为：a=2 3 6;b=1。这些函数的用法描述如下：1) impulse（）函数函数impulse（）将绘出由向量a和b表示的连续系统在指定时间内的冲激响应h(t)的时域波形图，并能求出指定时间范围内冲激响应的数值解。impulse（）函数有如下几种调用格式：impulse（b,a）该调用格式以默认方式绘出向量a和b定义的连续系统的冲激响应的时域波形。例题，求冲击响应的时域波形运行如下Matalb命令：a=1 3 2;b=1 0 1;impulse(b,a);则绘出系统的冲激响应波

6、形，如图4所示。图4impulse（b,a,t）该调用格式绘出向量a和b定义的连续系统在0t时间范围内冲激响应的时域波形。对上例，运行如下命令：impulse(b,a,10)则绘出系统在010秒范围内冲激响应的时域波形。impulse（b,a,t1:p:t2）该调用格式绘出向量a和b定义的连续系统在t1t2时间范围内,且以时间间隔p均匀取样的冲激响应的时域波形。对上例，运行如下命令：impulse(b,a,1:0.1:2)则绘出系统在12秒范围内, 且以时间间隔0.1秒均匀取样的冲激响应的时域波形。y=impulse（b,a,t1:p:t2）该调用格式并不绘出系统的冲激响应波形，而是求出向量a

7、和b定义的连续系统在t1t2时间范围内,且以时间间隔p均匀取样的冲激响应的数值解。如下命令：y=impulse(b,a,0:0.2:2)则计算出系统在02秒范围内, 且以时间间隔0.2秒均匀取样的冲激响应的数值解。2) step（）函数函数step（）将绘出由向量a和b定义的连续系统的阶跃响应g(t)在指定时间范围内的波形图，并能求出其数值解。和impulse（）函数一样， step（）函数也有四种调用格式：step(b,a)step(b,a,t)step(b,a,t1:p:t2)y=step(b,a,t1:p:t2)上述调用格式的功能与impulse（）函数完全相同。3) Lsim（）函数L

8、sim (b, a, x, t)，用来计算由b和a表示的LTI系统在输入信号x作用下的零状态响应。其中t为指定的时间变化范围，x为输入信号，它们的长度应该是相同的。例题：%系统方程为，的零状态响应。 ts=0;te=5;dt=0.01;t=ts:dt:te;f=10*sin(2*pi*t);a=1 2 100;b=1;y=lsim(b,a,f,t);%求取零状态响应plot(t,y); %画图xlabel(Time(sec) %横坐标定义ylabel(y(t) %纵坐标定义则绘出系统的零状态响应波形，如图5所示。图52、MATLAB中int函数用于卷积计算求积指令int的具体使用格式如下：in

9、f=int(f,v)：给出符号表达式f对指定变量v的(不带积分常数的)不定积分。int f=int(f,v,a,b)：给出符号表达式f对指定变量v的定积分：例题：已知试用int函数计算clearsyms t xe = exp(-3.*(x);h = exp(-2.*(t-x);i = int(e.*h, x, 0, t);ezplot(i, 0 6);xlabel(Time(sec);ylabel(y(t);title(exp(-3t)*exp(-2t) );grid on运行程序，f(t)的结果如图6所示。图6四、实验内容1. 分别用函数lsim和卷积积分两种方法求如图7所示系统的零状态响应

10、。其中L=1，R=2，。图72. 求上述系统的冲激响应与阶跃响应。五、思考题1. 为什么连续时间系统的零状态响应为激励与冲击响应的卷积？2. 利用卷积积分法计算系统响应应从几个方面进行?实验三 连续信号的频域分析一、实验目的1. 掌握周期信号的频谱 Fourier级数的分析方法。2. 深入理解信号频谱的概念，掌握典型信号的频谱以及Fourier变换的主要性质。3. 掌握调制与解调的基本原理及滤波器的使用。二、实验设备安装有matlab6.5以上版本的PC机一台。三、MATLAB使用说明1、MATLAB的符号积分函数int()可以帮助我们求出连续时间周期信号的截断傅立叶级数及傅立叶表示，此外su

11、bs()函数可以为符号变量赋值，例如subs(a+b,a,2)就是把a+b里面的a用2替换，subs(x/y,y,3)就是把x/y里面的y用3替换。2、MATLAB的傅立叶变换求取MATLAB提供了能直接求解傅立叶变换及逆变换的函数fourier()及ifourier()。两者的调用格式如下。(1)Fourier变换F=fourier(f)是符号函数f的Fourier变换，默认返回是关于的函数。即：。F=fourier(f,v)返回函数F是关于符号对象v的函数，而不是默认的，即。F=fouier(f,u,v) 对关于u的函数f进行变换，返回函数F是关于v的函数，即。(2)Fourier逆变换f

12、=ifourier(F)是符号函数F的Fourier变换，默认的独立变量为，默认返回是关于x的函数。f=ifourier(F,u)返回函数f是关于符号对象u的函数，而不是默认的x。f=ifouier(F,v,u) 对关于v的函数F进行变换，返回关于u的函数f 例题：求的傅里叶变换syms t x ;x=fourier(exp(-2*abs(t);ezplot(x);%符号函数作图函数xlabel(Time(sec);ylabel(f(t);grid on运行结果如图8所示。图8四、实验内容1、求如图9所示周期矩形脉冲信号的Fourier级数表达式，画出频谱图，并用前N次谐波合成的信号近似。图9

13、2、试用fourier()函数求下列信号的傅里叶变换,并画出(1) (2) 3、调制信号为一取样信号，利用MATLAB分析幅度调制(AM)产生的信号频谱，比较信号调制前后的频谱并解调已调信号。设载波信号的频率为100Hz。五、思考题1、根据试验1的结果，解释Gibbs现象。2、比较周期信号与非周期信号的频谱。3、调制与解调的基本原理是什么？为什么要进行调制？实验四 连续系统的频域分析一、实验目的1、掌握连续时间系统变换域分析的基本方法。 2、掌握系统无失真传输的基本条件。二、实验设备安装有matlab6.5以上版本的PC机一台。三、MATLAB使用说明MATLAB提供了专门对连续系统频谱响应H

14、(j)进行分析的函数freqs()。该函数可以求出系统频率响应的数值解，并可绘出系统的频率响应的数值解，并可绘出系统的幅频及相频响应曲线。假定系统函数的表达式为 ：Freqs()函数有如下四种调用格式：1、h=freqs(b,a,w)该调用格式中，b为对应于H(j)表达式中的向量b1,b2,b3,bm，a为对应于式中的向量a1,a2,a3,an，w为形如w1:p:w2的冒号运算定义的系统频率响应的频率范围，w1为频率起始值，w2为频率中止值，p为频率取样间隔。向量h则返回在向量w所定义的频率点上，系统频率响应的样值。例如，运行如下命令：a=1 2 1;b=0 1;h=freqs(b,a,0:0

15、.5:2*pi);%计算02频率范围内以间隔0.5取样的系统响应的样值。2、h,w=freqs(b,a)该调用格式将计算默认频率范围内200个频率点的系统频率响应的样值，并赋值给返回变量h,200个频率点记录在w中。3、h,w=freqs(b,a,n)该调用格式将计算默认频率范围内n个频率点上的系统频率响应的样值，并赋值给返回变量h,n个频率点记录在w中。4、freqs(b,a)改革是并不返回系统频率响应的样值，而是以对数坐标的方式绘出系统的幅频响应和相频响应曲线。四、实验内容1、如图10所示系统：图10 (a) 对不同的RC值，用freqs函数画出系统的幅频曲线。(b) 信号包含了一个低频分

16、量和一个高频分量，确定适当的RC值，滤除信号中的高频分量并画出信号和在s范围内的波形。提示：为最大值的处对应的频率为通带截止频率，首先求取并找到和RC关系，然后根据题意选定即可确定RC值。2、信号任选，分析以下几种情况下信号的频谱和波形变化：（1） 系统满足线性不失真条件时；（2） 系统只满足恒定幅值条件时；（3） 系统只满足相位条件时；（4）系统两个条件均不满足时。提示：利用fourier求取信号的傅立叶变换E(j)，然后设计使之满足不同条件，计算R( j)= E(j) H(j)并画频谱图。五、思考题1、连续系统频域与复频域分析的基本方法是什么？2、若信号经过系统不发生失真，则对系统频响有何

17、要求？实验五 信号采样与重建一、实验目的1、 深刻理解采样定理的内容。2、掌握信号采样后的频谱。 二、实验设备安装有matlab6.5以上版本的PC机一台。三、MATLAB使用说明利用fft计算频谱或者根据采样后所得的频谱公式构建采样后频谱。例题：以取样信号为例来分析采样后信号的频谱。syms t w;f=sin(2*t)/(2*t);subplot(3,2,1);ezplot(f,-5,5);title(时域信号);ff=fourier(f);subplot(3,2,2);ezplot(ff,-5,5);title(频谱信号);fs=1;fc=0;for n=-3:3fc=fc+(pi/2)

18、*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs);endsubplot(3,2,3);ezplot(fc,-30,30);title(采样信号频谱，采样频率1);四、实验内容：1、已知，以为采样频率，对进行采样得到，观察随着由小变大，频谱的变化，最后得出与两者频谱之间的关系。2、 由实验1中采样得到的离散信号重建对应的连续时间信号：情况一、；情况二、。提示根据内插公式重建连续时间信号，或者根据采样后频谱通过低通滤波器的形式重建连续时间信号。五、思考题1、随着采样频率从小到大变化，的频谱发生怎样的变化，与频谱之间的关系如何？与理论计算

19、结果之间是否完全一致？如果不一致，请分析可能是什么原因导致的？2、 采样频率分别满足情况一与情况二时，由重建的与原信号是否完全相同？如果不相同，试分析原因。实验六 离散时间信号和系统分析一、 实验目的1. 掌握常用离散信号的MATLAB表示方法。2、掌握用MATLAB计算卷积和及零状态响应的方法。二、 实验设备安装有matlab6.5以上版本的PC机一台。三、 实验说明1、离散信号的表示在MATLAB可视化绘图中， 而对n为自变量的离散序列，在绘图时统一用stem函数。例如：在MATLAB中应该表示为：k=-3,-2, -10,1,2,3 或是k=-3:3; f=1,2,-1,3,2,4,-1

20、; stem(k,f,filled); axis(-4,4, -1.5, 4.5); 2、离散序列卷积的计算MATLAB的conv()函数可以帮助我们快速求出两个离散序列的卷积和。conv函数的调用格式为：f=conv(f1,f2)。其中f1为包含序列f1(k)的非零样值点的行向量，f2为包含序列f2(k)的非零样值点的行向量，向量f则返回序列f(k)=f1(k)*f2(k)的所有非零样值点行向量。3、用离散信号的卷积计算连续时间信号的卷积卷积积分运算实际上可用信号的分段求和来实现，即如果我们只求当t=n（n为整数）时f(t)的值f(n),则由上式可得：当足够小时，就是卷积积分的结果连续时间信

21、号f(t)的较好的数值近似。由于MATLAB运算的特点，两个连续信号的卷积f(t)= f1(t)* f2(t)，用MATLAB实现的过程应为：(1)将连续信号f1(t)、f2(t)以时间间隔进行取样，得到离散序列f1(k )、f2(k )；(2)构造f1(k )、f2(k )与相对应的时间向量k1和k2；(3)调用conv()函数计算离散卷积和f1(k),再计算f(k)= f1 (k)*；(4)构造f(k )对应的时间向量k。四、 实验内容1、用MATLAB表示离散信号：，。2、有两离散序列，用MATLAB绘出它们的波形及f1(k)+f2(k)及f1(k)* f2(k)的波形。3、已知离散序列波形如图11所示，试用MATLAB绘出满足下列要求的序列波形。图11(1)f(k-2)(k) (2)f(-k) (3)f(-k+2) (4)f(k-2)(k-2)4、试用MATLAB的conv()函数计算实验2中第1题的结果。5、假设某系统的单位函数响应，系统激励信号，求系统的零状态响应。