信号与系统课程设计应用MATLAB实现连续信号的采样与重构仿真Word格式文档下载.docx

1、plot（x1）;plot（x2,y2）;plot（x3,y3,x4,y4）;其中，x1、x2、y2、x3、y3等符号均代表需要绘制的参数，它门可以是向量、矩阵、复数矩阵等，plot指令将根据不同的参数绘制不同的图形。2.1.3多子图绘制指令-subplotMATLAB为方便用户进行仿真分析，设置了subplot指令，利用它可以在不同的子图下绘制图形，以进行对比分析。subplot的基本指令格式如下：subplot（m,n,k） %作出（mn）幅子图中的第k幅图形subplot（position,left bottom width height） %在人工指定位置作出字图。subplot（m,

2、n,k）指令表示在图形窗口中产生（mn）幅子图，k代表当前绘制子图号。如subplot（2,2,1）就是产生22幅子图，当前在子图1绘制图形。2.1.4sinc序列Sinc函数是Matlab软件中经常使用的函数之一，sinc序列定义为： （2-1）这个信号可以利用Signal Processing Toolbox中的函数sinc来实现。2.2 周期信号的频谱周期信号可以分解成一系列正弦信号和指数信号之和，即f（t）= + cos（n0t+） （2-2）或f（t）= ejn （2-3） 沈 阳 大 学 课程设计说明书 NO.3图1 周期信号的频谱由图1见，周期信号频谱具有以下特点：（1）频谱图由

3、频率离散的谱线组成，每根谱线代表一个谐波分量。即周期信号的频谱是离散谱。（2）频谱图中的谱线只能在基波频率的整数倍频率上出现，即谐波性。（3）频谱中个谱线的高度，随谐波次数的增高而逐渐减小。当谐波次数无限增高时，谐波分量的振幅趋于无穷小，即收敛性。 （2-4）（2-5） 课程设计说明书 NO.4 （2-6） 图5 抽样函数波形 （2-7） 沈 阳 大 学课程设计说明书 NO.5 图2 周期与频谱的关系可见，信号的频带宽度与信号的持续时间成反比，信号持续时间愈长，其频带愈窄；反之，信号脉冲愈窄，其频带愈宽。2.3非周期信号的频谱 门函数可以表示为 （2-8） （2-9）图3 门函数及其频谱课程设

4、计说明书 NO.6（1）非周期矩形脉冲信号的频谱是连续频谱，其形状与周期矩形脉冲的离散频谱的包络线相似，都有Sa（x）的形式。周期信号的离散频谱可以通过对非周期信号的连续频谱等间隔取样求得。（2）信号在时域中的持续时间有限，则在频域其频谱将延续到无限。（3）信号的频谱分量主要集中在零频到第一个过零点之间，工程上往往将此宽度作为有效宽带。即以（0，）频率范围作为门函数的有效宽度。与脉冲宽度成反比，脉宽愈窄，频带愈宽；脉冲愈宽，频带愈窄。3、实验内容运用MATLAB软件，输入周期信号与非周期信号的频谱，并观察输出图形，并做出分析。具体程序以及步骤如下：3.1 周期信号频谱% 周期与频谱的关系实现程

5、序function A_sym,B_sym=CTFSingsym% 采用符号计算求0,T内时间函数的三角级数展开系数% 函数的输入输出都是数值量% Nn 输出数据的准确位数% A_sym 第1元素是直流项，其后元素依次是1,2,3,.次斜波cos项展开系数% B_sym 第2,3,4,.元素依次是1,2,3,.次斜波sin项展开系数% T T=m*tao, 周期信号% Nf 谐波的阶数% m（m=T/tao）周期与脉冲宽度之比，如m=4,8,16,100等% tao 脉宽:tao=T/Msyms t n y课程设计说明书 NO.7if nargin3;Nf=input（pleas Input

6、所需展开的最高谐波次数：Nf=）;endT=input（pleas Input 信号的周期T 5;Nn=32;y=time_fun_s（t）;A0=2*int（y,t,0,T）/T;As=int（2*y*cos（2*pi*n*t/T）/T,t,0,T）;Bs=int（2*y*sin（2*pi*n*t/T）/T,t,0,T）;A_sym（1）=double（vpa（A0,Nn）;for k=1:Nf A_sym（k+1）=double（vpa（subs（As,n,k）,Nn）;B_sym（k+1）=double（vpa（subs（Bs,n,k）,Nn）;if nargout=0 S1=flipl

7、r（A_sym） S1（1,k+1）=A_sym（1） S2=fliplr（1/2*S1） S3=fliplr（1/2*B_sym） S3（1,k+1）=0 S4=fliplr（S3） S5=S2-i*S4; N=Nf*2*pi/T; k2=0:2*pi/T:N; x=time_fun_e（t） subplot 212 stem（k2,abs（S5）; title（连续时间函数周期矩形脉冲的单边幅度谱） axis（0,80,0,0.12） line（0,80,0,0） %-function y=time_fun_s（t）% 该函数是CTFStpshsym.m的子函数。它是由符号变量和表达式写成

8、syms a a1课程设计说明书 NO.8pleas Input 信号的周期 T=M=input（周期与脉冲宽度之比M=A=1;tao=T/M;a=tao/2;y1=sym（Heaviside（t+a1）*A;y=y1-sym（Heaviside（t-a1）y=subs（y,a1,a）;y=simple（y）;%-function x=time_fun_e（t）% t是时间数组% T是周期duty=dao/Tt=-2*T:0.01:2*T;x=rectpuls（t,tao）;subplot 211plot（t,x）hold onx=rectpuls（t-T,tao）;x=rectpuls（t+

9、T,tao）;title（周期为T,脉宽tao=T/M的矩形脉冲axis（-10-T,10+T,0,1.2）课程设计说明书 NO.9 输出图形如下：图4 T=5,tao=4的周期矩形频谱图5 T=5,tao=8的周期矩形频谱课程设计说明书 NO.10图6 T=5,tao=16的周期矩形频谱课程设计说明书 NO.11图7 T=5,tao=100的周期矩形频谱图8 T=5,tao=4的周期矩形频谱课程设计说明书 NO.123.2非周期信号的频谱 % 门信号的波形及其频谱实现程序R=0.02;t=-2:R:2;f=heaviside（t+0.95）-heaviside（t-0.95）;w1=2*pi

10、*5;N=500;k=0:w=k*w1/N;F=f*exp（-j*t*w）*R;F=real（F）;W=-fliplr（w）,w（2:501）;F=fliplr（F）,F（2:subplot 211;plot（t,f,rxlabel（tylabel（f（t）f（t）=u（t+1）-u（t-1）subplot 212;plot（W,F,bwF（w）f（t）的傅立叶变换F（w）输出图形如下：课程设计说明书 NO.13图9 门函数及其频谱4、运行结果及结果分析由图4到图7可知，信号的频带宽度与信号的持续时间成反比，信号持续时间愈长，其频带愈窄；由图9可知，信号的频谱分量主要集中在零频到第一个过零点之

11、间，工程上往往将此宽度作为有效宽带。5、心得体会通过学习完信号与系统的课程设计，加深了我对该课程的了解，我觉得理论和实际是分不开的，本学期是我们第一次学习信号与系统基础课程，通过该实验巩固了的关于信号的理论，在完整的学习理论课程以后，通过课程设计的事件操作，我们对此门课程有更深的印象。信号抽样的定理是很一个重要的概念，它贯穿整个课程设计的始终。课程设计的草稿完成以后，自以为完成不错，但在老师的指导下发现有很多漏洞与不足，经过再三修改，最终完成令人比较满意的结果由于开始不熟悉MATLAB这个软件，走了很多弯路，后来上图书馆借阅很多相关课程设计说明书 NO.14的材料，经过多日仔细研究，然后在老师

12、的帮助下，最终能熟练运用该软件并顺利运行以上程序，在制作课程设计的过程当中，同时也发现了对此门课程自己还有不尽人意之处，比如对信号线形系统的分析，以及应用MATLAB软件进行仿真的程序，我还不能运用自如。但是，我相信经过日后系统的深入学习再加上我自己的努力，将来我可以将这些知识进行熟练的应用，并在实际生活当中应用这些知识来解决问题。6、参考文献1刘旭.MATLAB的系统分析与设计信号处理M.北京：北京电子科技大学出版社，1998.66-682王涌.MATLAB教程及其应用程序实例M.北京：北京航空航天大学出版社，2001.1-63张昱,周绮敏.信号与系统实验教程M.北京：人民邮电出版社，200

13、5.88-994赵毅，高原等.MATLAB7辅助信号处理技术与应用M.北京：科技工业出版社，2005.1-35吴大正，杨林耀，张永瑞.信号与线性系统分析M（第三版）.北京：高等教育出版社，1998.32-36 吴航行,王模昌. 傅立叶变换光谱仪信号调制度下降的分析J. 红外与毫米波学报, 2004,（05）7 范世福,李艳宁,赵玉春,孙振东. 傅里叶变换光谱经典计算技术J. 分析仪器, 1995,（02）8 杜培胜,郭靖,董前民. 液晶可调谐滤光片及其在光谱仪上的应用J. 红外, 2007,（11）9 相里斌. 傅里叶变换光谱仪中的主要技术环节J. 光子学报, 1997,（06） .10 曾立波,尹邦胜,何冰,龚丹. 傅里叶变换红外光谱仪动镜倾斜误差分析J. 光学精密工程, 2006,（02） .