信号与系统实验报告

1、信号与线性系统实验报告 班 级: 电科122 学 号: 124633224 姓 名: 纳扎尔库尔曼别克 2015年10月计算机与信息工程学院 班级 ： 电科122班 实验名称：实验一 常见信号的MATLAB表示及运算 指导教师： 张婧婧 姓名 ： 纳扎尔库尔曼别克 学号 ： 124633224 一、实验目的1熟悉常见信号的意义、特性及波形2学会使用MATLAB表示信号的方法并绘制信号波形3. 掌握使用MATLAB进行信号基本运算的指令4. 熟悉用MATLAB实现卷积积分的方法二、实验原理信号一般是随时间而变化的某些物理量。按照自变量的取值是否连续，信号分为连续时间信号和离散时间信号，一般用和来

2、表示。若对信号进行时域分析，就需要绘制其波形，如果信号比较复杂，则手工绘制波形就变得很困难，且难以精确。MATLAB强大的图形处理功能及符号运算功能，为实现信号的可视化及其时域分析提供了强有力的工具。根据MATLAB的数值计算功能和符号运算功能，在MATLAB中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MATLAB语句表示出信号后，就可以利用MATLAB中的绘图命令绘制出直观的信号波形了。下面分别介绍连续时间信号和离散时间信号的MATLAB表示及其波形绘制方法。作业题自编程序：1、已知两信号，求卷积积分，并与例题比较。t1=-1:0.01:0;t2=0:

3、0.01:1; t3=-1:0.01:1; f1=ones(size(t1); f2=ones(size(t2); g=conv(f1,f2); subplot(3,1,1),plot(t1,f1); subplot(3,1,2),plot(t2,f2); subplot(3,1,3),plot(t3,g);2. 已知，求两序列的卷积和 a=1,1,1,2; b=1,2,3,4,5; g=conv(a,b); h=0:1:7; stem(g,h,filled) stem(h,g,filled)班级 ： 电科122班 实验名称：实验二 LTI系统的响应 指导教师： 张婧婧 姓名 ： 纳扎尔库尔曼

4、别克 学号 ： 124633224 一、实验目的1.熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法2.熟悉连续（离散）时间系统在任意信号激励下响应的求解方法3.熟悉应用MATLAB实现求解系统响应的方法二、实验原理1.连续时间系统对于连续的LTI系统，当系统输入为f(t)，输出为y(t)，则输入与输出之间满足如下的线性常系数微分方程：，当系统输入为单位冲激信号(t)时产生的零状态响应称为系统的单位冲激响应，用h(t)表示。若输入为单位阶跃信号(t)时，系统产生的零状态响应则称为系统的单位阶跃响应，记为g(t)，如下图所示。系统的单位冲激响应h(t)包含了系统的固有特性，它是由系统本身

5、的结构及参数所决定的，与系统的输入无关。我们只要知道了系统的冲激响应，即可求得系统在不同激励下产生的响应。因此，求解系统的冲激响应h(t)对我们进行连续系统的分析具有非常重要的意义。在MATLAB中有专门用于求解连续系统冲激响应和阶跃响应， 并绘制其时域波形的函数impulse( ) 和step( )。如果系统输入为f(t)，冲激响应为h(t)，系统的零状态响应为y(t)，则有：。若已知系统的输入信号及初始状态，我们便可以用微分方程的经典时域求解方法，求出系统的响应。但是对于高阶系统，手工计算这一问题的过程非常困难和繁琐。在MATLAB中，应用lsim( )函数很容易就能对上述微分方程所描述的

6、系统的响应进行仿真，求出系统在任意激励信号作用下的响应。lsim( )函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解，而且还能同时绘制出系统响应的时域波形图。作业题自编程序：1. 已知描述系统的微分方程和激励信号e(t) 分别如下，试用解析方法求系统的单位冲激响应h(t)和零状态响应r(t)，并用MATLAB绘出系统单位冲激响应和系统零状态响应的波形，验证结果是否相同。； a=1 4 4;b=1 3; subplot(2,1,1), impulse(b,a,4) %冲激响应 subplot(2,1,2), step(b,a,4) %阶跃响应 a=1 4 4;b=1 3; p1=

7、0.01; t1=0:p1:5; x1=exp(-t1); lsim(b,a,x1,t1), hold on; p2=0.5; t2=0:p2:5; x2=exp(-t2); lsim(b,a,x2,t2), hold off2. 请用MATLAB分别求出下列差分方程所描述的离散系统，在020时间范围内的单位函数响应、阶跃响应和系统零状态响应的数值解，并绘出其波形。另外，请将理论值与MATLAB仿真结果在对应点上的值作比较。；a=1,2,1; b=1; impz(b,a,-3:10), title(单位响应) a=1 2 1;b=1;subplot(2,1,1), dstep(b,a,20)

8、%阶跃响应a=1,2,1; b=1,0,0;N=20;n=0:N-1;hn=impz(b,a,n);gn=dstep(b,a,n);subplot(2,1,1),stem(n,hn);title(单位冲击响应);xlabel(n);ylabel(h(n);subplot(2,1,2),stem(n,gn);title(单位阶跃响应);xlabel(n);ylabel(g(n);a=1,2,1;b=1k=0:20; %定义输入序列取值范围x=(1/4).k; %定义输入序列表达式y=filter(b,a,x) %求解零状态响应样值subplot(2,1,1),stem(k,x) %绘制输入序列的

9、波形 title(输入序列)subplot(2,1,2),stem(k,y) %绘制零状态响应的波形title(输出序列)a=1 2 1;b=1;p1=0.01; %定义取样时间间隔为0.01t1=0:p1:5; %定义时间范围x1=exp(-2*t1); %定义输入信号lsim(b,a,x1,t1), %对取样间隔为0.01时系统响应进行仿真hold on; %保持图形窗口以便能在同一窗口中绘制多条曲线p2=0.5; %定义取样间隔为0.5t2=0:p2:5; %定义时间范围x2=exp(-2*t2); %定义输入信号lsim(b,a,x2,t2), hold off %对取样间隔为0.5时

10、系统响应进行仿真并解除保持班级 ： 电科122班 实验名称：实验三 连续时间信号的频域分析 指导教师： 张婧婧 姓名 ： 纳扎尔库尔曼别克 学号 ： 124633224 一、实验目的1熟悉傅里叶变换的性质2熟悉常见信号的傅里叶变换3了解傅里叶变换的MATLAB实现方法二、实验原理傅里叶变换是信号分析 的最重要的内容之一。从已知信号求出相应的频谱函数的数学表示为：的傅里叶变换存在的充分条件是在无限区间内绝对可积，即满足下式：但上式并非傅里叶变换存在的必要条件。在引入广义函数概念之后，使一些不满足绝对可积条件的函数也能进行傅里叶变换。 傅里叶反变换的定义为：。在这一部分的学习中，大家都体会到了这种

11、数学运算的麻烦。在MATLAB语言中有专门对信号进行正反傅里叶变换的语句，使得傅里叶变换很容易在MATLAB中实现。在MATLAB中实现傅里叶变换的方法有两种，一种是利用MATLAB中的Symbolic Math Toolbox提供的专用函数直接求解函数的傅里叶变换和傅里叶反变换，另一种是傅里叶变换的数值计算实现法。作业题自编程序：1.编程实现求下列信号的幅度频谱 求出的频谱函数F1(j)，请将它与上面门宽为2的门函数的频谱进行比较，观察两者的特点，说明两者的关系。syms t w %定义两个符号变量t,wGt=sym(Heaviside(2*t+1)-Heaviside(2*t-1); %产

12、生门宽为2的门函数Fw=fourier(Gt,t,w); %对门函数作傅氏变换求F(jw)FFw=maple(convert,Fw,piecewise); %数据类型转换，转为分段函数，此处可以去掉FFP=abs(FFw); %求振幅频谱| F(jw)|ezplot(FFP,-10*pi 10*pi);grid; %绘制函数图形，并加网格axis(-10*pi 10*pi 0 2.2) FFw= -i*exp(i*w)/w+i*exp(-i*w)/w % FFw为复数FFP= abs(-i*exp(i*w)/w+i*exp(-i*w)/w)2.利用ifourier( ) 函数求下列频谱函数的傅

13、氏反变换 syms t w %定义两个符号变量t,wFw=sym(-j*2*w/(16+w2); %定义频谱函数F(jw)ft=ifourier(Fw,w,t); %对频谱函数F(jw)进行傅氏反变换班级 ： 电科122班 实验名称：实验四 系统的零极点及频率响应特性 指导教师： 张婧婧 姓名 ： 纳扎尔库尔曼别克 学号 ： 124633224 一、实验目的1掌握系统函数零极点的定义2熟悉零极点与频率响应的关系3掌握极点与系统稳定性的关系4状态方程与系统函数的关系5在MATLAB中实现系统函数与状态方程间的转换二、实验原理描述连续系统的系统函数H(s)的一般表示形式为：其对应的零极点形式的系统

14、函数为： 共有n个极点：p1，p2，pn和m个零点：z1，z2，zm。把零极点画在S平面中得到的图称为零极点图，人们可以通过零极点分布判断系统的特性。当系统的极点处在S的左半平面时系统稳定；处在虚轴上的单阶极点系统稳定；处在S的右半平面的极点及处在虚轴上的高阶极点，系统是不稳定的。描述系统除了可以用系统函数和零极图以外，还可以用状态方程。对应上述用系统函数H (s) 描述的系统，其状态方程可用相变量状态方程和对角线变量状态方程描述，作业题自编程序：1已知下列系统函数H (s)或状态方程，求其零极点，并画出零极点图。 num = 1 0 1; %分子系数，按降幂顺序排列。den = 1 2 5;

15、 %分母系数，按降幂顺序排列。z,p = tf2zp(num,den); %求零点z和极点pzplane (z,p) %作出零极点图 2已知下列系统函数H (s)，求其频率特性。 num = 0 2 0; %分子系数，按降幂顺序排列。den = 1 sqrt(2) 1; %分母系数，按降幂顺序排列。 w =logspace (-1,1); %频率范围freqs(num,den,w) %画出频率响应曲线 3. 已知系统函数H (s)，求其频率特性和零极点图。num = 1 35 291 1093 1700; %分子系数，按降幂顺序排列。den = 1 9 66 294 1029 2541 468

16、4 5856 4629 1700; %分母系数，按降幂顺序排列。z,p = tf2zp(num,den); %求零点z和极点pzplane (z,p) %作出零极点图 w =logspace (-1,1); %频率范围freqs(num,den,w) %画出频率响应曲线班级 ： 电科122班 实验名称：实验五 连续信号与系统的S域分析 指导教师： 张婧婧 姓名 ： 纳扎尔库尔曼别克 学号 ： 124633224 一、实验目的1. 熟悉拉普拉斯变换的原理及性质2. 熟悉常见信号的拉氏变换3. 了解正/反拉氏变换的MATLAB实现方法和利用MATLAB绘制三维曲面图的方法4. 了解信号的零极点分布

17、对信号拉氏变换曲面图的影响及续信号的拉氏变换与傅氏变换的关系二、实验原理拉普拉斯变换是分析连续时间信号的重要手段。对于当t时信号的幅值不衰减的时间信号，即在f(t)不满足绝对可积的条件时，其傅里叶变换可能不存在，但此时可以用拉氏变换法来分析它们。连续时间信号f(t)的单边拉普拉斯变换F(s)的定义为：拉氏反变换的定义为： 显然，上式中F(s)是复变量s的复变函数，为了便于理解和分析F(s)随s的变化规律，我们将F(s)写成模及相位的形式：。其中，|F(s)|为复信号F(s)的模，而为F(s)的相位。由于复变量s=+j，如果以为横坐标（实轴），j为纵坐标（虚轴），这样，复变量s就成为一个复平面，

18、我们称之为s平面。从三维几何空间的角度来看，和分别对应着复平面上的两个曲面，如果绘出它们的三维曲面图，就可以直观地分析连续信号的拉氏变换F(s)随复变量s的变化情况，在MATLAB语言中有专门对信号进行正反拉氏变换的函数，并且利用 MATLAB的三维绘图功能很容易画出漂亮的三维曲面图。作业题自编程序：1. 求出下列函数的拉氏变换式，并用MATLAB绘制拉氏变换在s平面的三维曲面图 Syms t s %定义符号变量ft=sym(2*(e(-t)*Heaviside(t)+5*(e(-3*t)*Heaviside(t); %定义时间函数f(t)的表达式Fs=laplace(ft) %求f(t)的拉

19、氏变换式F(s)syms x y se=2.718281829s=x+i*y; %产生复变量sFFs=2/(s+log(e)+5/(s+3*log(e); %将F(s)表示成复变函数形式FFss=abs(FFs); %求出F(s)的模ezmesh(FFss); %画出拉氏变换的网格曲面图ezsurf(FFss); %画出带阴影效果的三维曲面图colormap(hsv); 2. 已知信号的拉氏变换如下，请用MATLAB画出其三维曲面图，观察其图形特点，说出函数零极点位置与其对应曲面图的关系，并且求出它们所对应的原时间函数f (t) syms t s %定义符号变量Fs =sym(2*(s-3)*

20、(s+3)/(s-5)*(s2+16); %定义F(s)的表达式ft=ilaplace(Fs) %求F(s)的拉氏反变换式f(t)ft=400/41*cos(t)4-400/41*cos(t)2+500/41*sin(t)*cos(t)3-250/41*sin(t)*cos(t)+2*exp(5/2*t)*(cosh(5/2*t)-9/41*sinh(5/2*t)syms t s %定义符号变量ft=sym(400/41*cos(t)4-400/41*cos(t)2+500/41*sin(t)*cos(t)3-250/41*sin(t)*cos(t)+2*exp(5/2*t)*(cosh(5/

21、2*t)-9/41*sinh(5/2*t)*Heaviside(t); %定义时间函数f(t)的表达式Fs=laplace(ft) Fs=9600/41/s/(s2+4)/(s2+16)*(1+1/2*s2+1/24*s2*(s2+4)-800/41/s/(s2+4)*(1+1/2*s2)+125/328/(1/16*s2+1)+125/82/(1/4*s2+1)-250/41/(s2+4)+8/25*(s-5/2)/(4/25*(s-5/2)2-1)-36/205/(4/25*(s-5/2)2-1)syms x y ss=x+i*y; %产生复变量sFFs=9600/41/s/(s2+4)/

22、(s2+16)*(1+1/2*s2+1/24*s2*(s2+4)-800/41/s/(s2+4)*(1+1/2*s2)+125/328/(1/16*s2+1)+125/82/(1/4*s2+1)-250/41/(s2+4)+8/25*(s-5/2)/(4/25*(s-5/2)2-1)-36/205/(4/25*(s-5/2)2-1); %将F(s)表示成复变函数形式FFss=abs(FFs); %求出F(s)的模ezmesh(FFss); %画出拉氏变换的网格曲面图ezsurf(FFss); %画出带阴影效果的三维曲面图colormap(hsv); %设置图形中多条曲线的颜色顺序 班级 ： 电

23、科122班 实验名称：实验六 离散信号与系统的Z变换分析 指导教师： 张婧婧 姓名 ： 纳扎尔库尔曼别克 学号 ： 124633224 一、实验目的1.熟悉离散信号Z变换的原理及性质2.熟悉常见信号的Z变换3.了解正/反Z变换的MATLAB实现方法4.了解离散信号的Z变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系5.了解利用MATLAB实现离散系统的频率特性分析的方法二、实验原理正/反Z变换Z变换分析法是分析离散时间信号与系统的重要手段。如果以时间间隔对连续时间信号f(t)进行理想抽样，那么，所得的理想抽样信号为：理想抽样信号的双边拉普拉斯变换Fd (s)为：若令 ，那么的双边拉普拉斯变换F (s)为：则离散信号f(k)的Z变换定义为： 从上面关于Z变换的推导过程中可知，离散信号f(k)的Z变换F(z)与其对应的理想抽样信号的拉氏变换Fd (s)之间存在以下关系： 同理，可以