大学计算机实验报告信息的表示与转换,大连理工大学城市学院信号与系统上机实验报告.doc...

大连理工大学城市学院信号与系统上机实验报告.doc

大连理工大学城市学院 信号与系统实验报告 学 院电子与自动化学院专 业 电子信息工程 学 生 秦俊 学 号 201115128 班 级 电子1105班 19目 录实验一 信号的Matlab表示31.1 实验原理31.2 实验内容3实验二 连续系统的时域分析82.1 实验原理82.2 实验内容8实验三 连续系统的频域分析123.1 实验原理123.2 实验内容12实验四 连续系统的S域分析174.1 实验原理174.2 实验内容18实验一 信号的Matlab表示1.1 实验原理信号是随时间变换的物理量。时域信号是指将信号表示成时间的函数，信号的时间特性是指信号的波形出现的先后、持续时间的长短、随时间变化的快慢等。离散信号是只在某些不连续的时间上有信号值，而在其他时间点上信号没有定义的一类信号。离散信号一般可以利用模数转换由连续信号来得到，计算机所能处理的只是离散信号。常用的连续信号有直流信号、正弦信号、单位阶跃信号、单位冲激信号、抽样信号等。常用的离散信号有正弦信号序列、单位阶跃序列、单位冲激序列等。信号的运算包括信号的基本运算，信号的时域变换等。1.2 实验内容1、 单位阶跃信号第一种实现方式的程序及图程序如下 t-20.016;ut0;plott,yaxis-2,6,0,1.2;结果如下图第二种实现方式的程序及图程序如下clcclearclose allt-20.016;ustepfunt,0;plott,uaxis-2,6,0.1,2;结果如下第三种实现方式的程序及图程序如下clcclearclose allt-20.016; u1.*heavisidet;plott,uaxis-2,6,0,1.2;结果如下三种方式得到的结果都一致的2、 u k的程序和图如下3、 正弦序列的程序和图形如下程序clcclearclose allk -2020;fsink*pi/6;stemk,f,filled;图形 4、 ftt*ut-ut-1 f-2t1clcclearclose allt-1.50.0011.5;fbht;f1bh-2*t1;subplot1,2,1plott,fsubplot1,2,2plotf,f1function fbht;ft.*t0-t1;图形 实验二 连续系统的时域分析一、实验原理系统的响应一般包括两个部分，即由当前输入所产生的响应(零状态响应)和由历史输入(初始状态)所产生的响应(零输入响应)。对于低阶系统，一般可以通过解析的方法得到响应。但对于高阶系统，手工计算就比较困难，这时MATLAB强大的计算功能就能比较容易地确定系统的各种响应，如冲激响应、阶跃响应、零输入响应、零状态响应、全响应等。二 连续系统时域分析实验内容1、 完成卷积的计算。要求(1)应用MATLAB符号计算方法 (2)编程实现，并同时显示xt、ft及yt的图形 (3)对于结论要有验证程序如下clcclearclose allsyms t tao;xexp-t;fexp-2*t;ysubsx,t,tao*subsf,t,t-tao;y1inty,tao,0,t;subplot2,2,1;ezplotx,-8,8;grid onsubplot2,2,2;ezplotf,-8,8;grid on;subplot2,2,3;ezploty1,-8,8;grid on;图形如下2、求系统的冲激响应和阶跃响应。冲激响应程序如下clcclearclose allb3,9;a1,6,8;systfb,a;t00.110;yimpulsesys,t;plott,y结果图 阶跃响应程序如下clcclearclose allb3,9;a1,6,8;systfb,a;t00.110;ystepsys,t;plott,y结果图如下 3、 求系统的零状态响应。clcclear 结果图close all b1;a1,0,1;systfb,a;t00.110;xcost;ylsimsys,x,t;plott,y实验三 连续系统的频域分析一实验原理根据傅里叶级数的原理，任何周期信号都可以分解为三角级数的组合，称为的傅里叶级数。例如一个方波信号可以分解为傅里叶变换傅里叶逆变换求解傅里叶变换，可以调用fourier函数，调用格式为Ffourierf,u,v，是关于u的函数f的傅里叶变换，返回函数F是关于v的函数。求解傅里叶逆变换，可以调用ifourier函数，调用格式为fifourierF,u,v，是关于v的函数F的傅里叶变换，返回函数f是关于u的函数。二实验内容1 信号的分解用正弦信号的叠加近似合成一个频率为50Hz，幅度为3的方波clcclearclose allfs10000; t01/fs0.1;f050;sum0;subplot2,1,1for n129plott,4/pi*1/n*sin2*pi*n*f0*t,k;hold on;endtitle信号叠加前;subplot2,1,2for n129sumsum4/pi*1/n*sin2*pi*n*f0*t; 三个波叠加后的波形endplott,sum,k;title信号叠加后;图形如下2、傅里叶变换和傅里叶逆变换例2 程序如下clcclear close allsyms t; Ffourierexp-2*abst;ezplotF 结果如图 例3 已知连续信号，通过程序求其傅里叶逆变换。程序如下clcclear close allsyms t w;ifourier1/1w2,t;结果图如下3 理解例2及例3中对于傅里叶变换和逆变换的命令应用，结合教材，完成对傅里叶变换对称性的验证。程序如下clcclearclose allsyms w1 w2;ftsymft;Fw1fourierft,w1Fw2fourierFw1,w1,w2图形如下4、已知系统微分方程。编程绘制系统的幅频响应，相频响应，频率响应的实部和虚部。b1; 分母的系数a1,3,2; 分子的系数H,wfreqsb,a; 求解微分方程HmabsH; 相频响应phaiangleH; 幅频响应HrrealH; 频率响应的实部HiimagH; 频率响应的虚部subplot2,2,1plotw,Hm 输出相频响应的图像title相频响应grid onsubplot2,2,2plotw,Hrtitle(相频响应的实部)grid onsubplot2,2,3plotw,phaititle(幅频响应)grid onsubplot2,2,4plotw,Hititle(频率响应的虚部)grid on实验结果 实验四 连续系统的S域分析一实验原理连续时间信号的拉普拉斯变换定义为 拉普拉斯逆变换定义为 。考虑到实际问题，人们用物理手段和实验方法所能记录和处理的一切信号都是有起始时刻的，对于这类单边信号或因果信号，我们引入单边拉普拉斯变换，定义为。如果连续时间信号可用符号表达式表达，则可利用MATLAB的符号数学工具箱中的laplace函数来实现其单边拉普拉斯变换，其语句格式为Llaplacef式中L返回的是默认符号为自变量s的符号表达式；f则为时域符号表达式，可通过sym函数来定义。如果连续时间信号可用符号表达式表达，则可利用MATLAB的符号数学工具箱中的ilaplace函数来实现其单边拉普拉斯变换，其语句格式为filaplaceL式中f返回的是默认符号为自变量t的符号表达式；L则为时域符号表达式，可通过sym函数来定义。用MATLAB函数residue可得到复杂有理分式Fs的部分分式展开式，其语句格式为 r,p,kresidueB,A其中，B、A分别表示Fs的分子和分母多项式的系数向量；r为部分分式的系数；p为极点；k为Fs中整式部分的系数。若Fs为有理真分式，则k为0。二实验内容【 1 】试用MATLAB的laplace函数求的拉普拉斯变换。 MATLAB程序如下 clcclearclose allfsymexp-t*sina*t;Llaplacef图形 【 2 】试用MATLAB的ilaplace函数求的拉普拉斯变换。 MATLAB程序如下 Fsyms2/s21;ftilaplaceF【 3 】利用MATLAB部分分式展开法求的反变换。 【4】分别用手算和MATLAB部分分式展开法求的反变换。【 5 】已知系统函数为，试用MATLAB命令画出其零极点分布图，并判断该系统是否稳定。【 6 】已知系统函数为，试用MATLAB命令求出零极点，并画出其零极点分布图，同时判断该系统是否稳定。b1,4,3; a1,1,7,2;systfb,a;pzmapsysaxis-3,3,-2,2