【报告】《信号与系统》课程研究性学习报告

1、精品word可编辑资料- - - - - - - - - - - - -信号与系统课程讨论性学习报告2021.11指导老师薛健时间第 10 页，共 7 页- - - - - - - - - -信号与系统的时域分析专题研讨【目的 】(1) 加深对信号与系统时域分析基本原理和方法的懂得；(2) 学会利用 matlab进行信号与系统的分析；(3) 培育同学自主学习才能，以及发觉问题、分析问题和解决问题的才能；【研讨题目 】 题目 1连续信号通过系统的响应一连续 lti 系统满意的微分方程为(1) 已知xte 3t ut ，试求该系统的零状态响应yzs t ；(2) 用 lsim 求出该系统的零状态响

2、应的数值解；利用1所求得的结果，比较不同的抽样间隔对数值解精度的影响；(3) 用命令x,fsam,bits = wavreadyourn；将硬盘上的语音文件 yourn.wav读入运算机；用命令soundx,fsam；播放该语音信号；(4) 用命令load model01将磁盘文件 model01.mat读入运算机后， matlab的 workspace 中将会新增变量 den 和 num；den 表示微分方程左边的系数，变量num表示微分方程右边的系数；写出磁盘文件model01.mat定义的微分方程；(5) 运算 3 中的信号通过 4 中系统的响应，播放系统输出的语音信号；与处理前的信号比

3、较，信号有何不同？能用已学学问说明所得结果吗？【题目目的 】1. 学会用运算机求解信号通过系统响应；2. 熟识用 matlab处理语音信号的基本命令；【仿真结果 】（ 1） 解出 y=-2.5*exp-3*t+3*exp-2*t-0.5*exp-t（ 2）【结果分析 】由图可知，抽样间隔越小，精度越高；【仿真结果 】（ 3）杂 音 音 频0.60.40.20-0.2-0.4-0.6【结果分析 】-0.802468101214161820可以听出与看出高频部分存在杂音；【仿真结果 】（ 4） untitled3 den=den= 1.309536e+04den= 7.076334e+08 den

4、= 6.939120e+12 den= 1.396319e+17 den= 8.396151e+20 den= 5.648432e+24 num= 3.162278e-03 num= 9.235054e-14 num= 1.649476e+07 num= 3.566819e-04 num= 1.646178e+16 num= 1.058969e+05 num= 4.486709e+24【仿真结果 】（ 5）【结果分析 】杂音部分被去掉，音乐复原正常；该系统有去噪的成效相当于一个低通录波器，高频的噪声信号被滤去；【研讨题目 】题目 2连续信号卷积积分的数值近似运算两个连续信号的卷积积分定义为为了

5、能用数值方法进行运算，需对连续信号进行抽样；记xk= xk , hk= hk ,为进行数值运算所选定的抽样间隔，就可证明连续信号卷积积分可近似的表示为yk x k hk 1由式 1可知，可以利用 matlab 供应的 conv 函数近似运算连续信号的卷积积分；设xt= ut ut 1，ht= xt xt，(a) 为了与近似运算的结果作比较，用解析法求出yt= xt ht；(b) 用不同的 运算出卷积积分的数值近似值，并与a 中的理论结果进行比较；(c) 证明 1式成立；(d) 如 xt和 ht不是时限信号，如xt= ut， ht=e tut ,， 就用上面的方法进行近似运算是否会遇到问题？如显

6、现问题请分析显现问题的缘由，并给出一种解决问题的方案；依据提出的方案完成近似运算卷积分的子程序；(e) 卷积函数 conv 新选项的定义与应用讨论在新版 matlab中，卷积函数conv 供应了选项conva, b, valid ，下面将讨论conva,b, valid 的定义及应用；(1) 读 matlab供应的关于 conv 的 help ，给出卷积函数 conva, b, valid 的定义；设计一些简洁的试验，验证你给出的定义；你认为这样定义的卷积有何优缺点？(2) 能否用 conva, b,valid 完成 conva, b ？请给出说明，编程验证你的观点；(3) 探究能否用 con

7、va, b, valid 解决 d 中显现的问题？如行，请给出解决问题的方案及试验结果；如不行，请给出理由；(4) 关于卷积 conva, b, valid 的应用，你仍有什么见解？欢迎发表新想法，胆要大，不要怕犯错；【题目目的 】1. 学会用运算机近似运算连续信号的卷积积分；2. 分析在运算过程中显现的问题并提出解决方案；【题目分析 】（ a） ht xt *xt ut *ut ut1 *ut1 - 2ut *ut1r t r t2 - 2r t1【仿真结果 】绿线是间隔的结果，蓝线是间隔结果，红线是实际结果，可见时间间隔越小，精确度越高；c证明：【仿真结果 】d【结果分析 】卷积结果在 t

8、=10 时刻发生突变， 与实际结果不符， 缘由是 t 只取 0-10 这一范畴，在大于 10 区域 matlab默认信号为 0，造成突变；解决方案，依据实际需要确定t 的范畴，取坐标时只显示有实际含义的范畴内的曲线；【仿真结果 】2【结果分析 】运算卷积只会运算序列完全重合的部分，略去 matlab 对不重合部分默认置0 对结果可能造成的影响；valid 有一个显着缺点， 卷积运算时，只会将b 翻转与 a ，比较， 而不会自动挑选短的序列进行翻转， 当 b 的长度大于 a 时，无法得出卷积结果；【仿真结果 】（ 3）可以；如 b 的长度为 n，就在 a 后面添加 n-1 个 0；（ 3b）在

9、a 前加 n-1 个 0，使用 valid4可以编一个挑选程序， 比较 a， b长度，如 b序列比 a 长就交换 a,b序列，这样就可以算出卷积；附录代码（ 2） t1=0:0.01:5; t2=0:0.1:5; t3=0:0.5:5;sys=tf2 1,1 3 2; x1=exp-3*t1; x2=exp-3*t2; x3=exp-3*t3; y1=lsimsys,x1,t1; y2=lsimsys,x2,t2; y3=lsimsys,x3,t3;y=-2.5*exp-3*t1+3*exp-2*t1-0.5*exp-t1;plott1,y1,b,t2,y2,r,t3,y3,c,t1,y,k;

10、 legenddt=0.01, dt=0.1, dt=0.5;title 零状态响应 ;xlabelt;ylabely;（ 3）x,fsam,bits = wavreadyourn soundx,fsamsl=lengthx; t=0:sl-1/fsam;wavwritex,fsam,yourns.wav; plott,x;title杂音音频 ;4load model01 fprintfden= %sn,den; fprintfnum= %sn,num;5x,fsam,bits = wavreadyourn; load model01;sys = tfnum,den; t=linspace0,

11、20,220500; y=lsimsys,x,t; soundy,fsam;（ b）t=0:0.1:5;t1=0:0.1:5; t2=0:0.01:5 x1=t1=0-t1=1;x2=t2=0-t2=1;ht1=convnx1,x1; y1=convnx1,ht1; ht2=convnx2,x2; y2=convnx2,ht2;y=1/2*t.*t.*t=0-3/2*t-1.*t-1.*t=1+3/2*t-2.*t-2.*t=2-1/ 2*t-3.*t-3.*t=3;n1=lengthy1; n2=lengthy2; plott,y,r;plot0:n2-1*0.01,y2*0.01*0.01

12、,b;plot0:n1-1*0.1,y1*0.1*0.1,g;axis0 5 0 1;d23dt=0.01; t=0:dt:10; x=t=0;y=exp-t.*t=0; y1=convnx,y; l=lengthy1;plot0:l-1*dt,y1*dt*dt;y=1 1 1;x=0 1 1 1 1 0;z=convx,y; z1=convx,y,valid; a=lengthz; a1=lengthz1; t=0:a-1;t1=0:a1-1;subplot211 stemt,z;axis-1 8 -1 6;subplot212 stemt1,z1; axis-1 8 -1 6titleht=y*x;y=1 1 1;x=0 1 1 1 1 0;n=lengthy;x1=zeros1,n-1,x,zeros1,n-1; z=convx1,y,valida=lengthz; t=0:a-1;stemt,z;axis-1 8 -1 6;（ 3 ）y=1 1 1;3bx=0 1 1 1 1 0;n=lengthy;x1=zeros1,n-1,x,zeros1,n-1; z=convx