信号与系统课程设计报告分析

沈阳理工大学课程设计专用纸1沈阳理工大学成绩评定表学生姓名张丽班级学号1203060101专业通信工程课程设计题目傅里叶变换的对称性和时移特性评语组长签字：成绩日期2014年06月日沈阳理工大学课程设计专用纸2沈阳理工大学课程设计任务书学院信息科学与工程学院专业通信工程专业学生姓名张丽班级学号1203060101课程设计题目傅里叶变换的对称性和时移特性实践教学要求与任务:1、学习Matlab软件及应用；2、学习并研究信号的复频域有关理论；3、利用Matlab编程，完成系统函数的零极点分布决定时域特性课题；4、写出课程设计报告，打印程序，给出运行结果。工作计划与进度安排:第1-2天：1、学习使用Matlab软件、上机练习2、明确课题内容，初步编程第3-5天：1、上机编程、调试2、撰写课程设计报告书3、检查编程、运行结果、答辩4、上交课程设计报告指导教师：2014年06月日专业负责人：2014年06月日学院教学副院长：2014年06月日沈阳理工大学课程设计专用纸3沈阳理工大学摘要本文研究的是傅里叶变换的对称性和时移特性，傅里叶变换的性质有：对称性、线性（叠加性）、奇偶虚实性、尺度变换特性、时移特性、频移特性、微分特性、积分特性、卷积特性（时域和频域）；从信号与系统的角度出发,给出了激励信号的具体模型；应用Matlab软件进行仿真，将研究的信号转化成具体的函数形式，在Matlab得到最终变换结果。使用傅里叶变换的方法、卷积的求解方法以及函数的微分等方法研究题目。关键词:傅里叶变换;对称性;时移特性;Matlab沈阳理工大学课程设计专用纸4沈阳理工大学目录1、Matlab介绍...................................................................................错误!未定义书签。2.利用Matlab实现信号的频域分析—傅里叶变换的对称性与时移特性设计......................................................................................................................52.1.傅里叶变换的定义及其相关性质....................................................................52.2.傅里叶变换的对称性验证编程设计及实现................................................72.3.傅里叶变换的时移特性验证编程设计及实现..........................................113.总结......................................................................................................................................134.参考文献..........................................................................................................................13沈阳理工大学课程设计专用纸5沈阳理工大学1、Matlab介绍MATLAB作为一种功能强大的工程软件，其重要功能包括数值处理、程序设计、可视化显示、图形用户界面和与外部软件的融合应用等方面。MATLAB软件由美国MathWorks公司于1984年推出，经过不断的发展和完善，如今己成为覆盖多个学科的国际公认的最优秀的数值计算仿真软件。MATLAB具备强大的数值计算能力，许多复杂的计算问题只需短短几行代码就可在MATLAB中实现。作为一个跨平台的软件，MATLAB已推出Unix、Windows、Linux和Mac等十多种操作系统下的版本，大大方便了在不同操作系统平台下的研究工作。MATLAB软件具有很强的开放性和适应性。在保持内核不变的情况下，MATLAB可以针对不同的应用学科推出相应的工具箱(toolbox)，目前己经推出了图象处理工具箱、信号处理工具箱、小波工具箱、神经网络工具箱以及通信工具箱等多个学科的专用工具箱，极大地方便了不同学科的研究工作。国内已有越来越多的科研和技术人员认识到MATLAB的强大作用，并在不同的领域内使用MATLAB来快速实现科研构想和提高工作效率。MATLAB提供了20类图像处理函数,涵盖了图像处理的包括近期研究成果在内的几乎所有的技术方法,是学习和研究图像处理的人员难得的宝贵资料和加工工具箱。这些函数按其功能可分为:图像显示;图像文件I/O;图像算术运算;几何变换;图像登记;像素值与统计;图像分析;图像增强;线性滤波;线性二元滤波设计;图像去模糊;图像变换;邻域与块处理;灰度与二值图像的形态学运算;结构元素创建与处理;基于边缘的处理;色彩映射表操作;色彩空间变换;图像类型与类型转换。沈阳理工大学课程设计专用纸6沈阳理工大学2、利用Matlab实现信号的频域分析—傅里叶变换的对称性时移特性的设计2.1傅里叶变换的定义及其相关性质1用周期信号的傅里叶级数通过极限的方法导出的非周期信号频谱称为傅里叶变换。傅里叶正变换为:dttfwFejwt)()(傅里叶逆变换为:dwwFtfejwt)(π21)(式中)(wF是)(tf的频谱函数，它一般是复函数，可以写作φ(w))()(jewFwF。其中)(wF是)(wF的模，它表示信号中各频率分量的相对大小。φ(w)是)(wF的相位函数，它表示信号中各频率分量之间的相位关系。2傅立叶变换的频移性质若)()(wFtf，则)0()(0Ftjetf沈阳理工大学课程设计专用纸7沈阳理工大学结论：将信号)(tf乘以因子tje0，对应于将频谱函数沿轴右移0；将信号)(tf乘以因子tje0，对应于将频谱函数沿轴右移0。3傅立叶变换的尺度变换性质若)()(Ftf，则对于任意实常数a，则有)(1)(aFaatf结论：信号时域波形的压缩，对应其频谱图形的扩展；而时域波形的扩展对应其频谱图形的压缩，且两域内展缩的倍数一致。4傅立叶变换的对称特性若)()(Ftf，则)(2)(ftF5傅立叶变换的时域卷积特性若)(*)()(11tftftf)()(Ftf)()(11Ftf则)()()(11FFF上式表明：如果函数)(tf的频谱为)(F，函数)(1tf的频谱为)(1F，且)(*)()(11tftftf，那么)()()(11FFF沈阳理工大学课程设计专用纸8沈阳理工大学2.2傅里叶变换的对称性验证编程设计及实现傅里叶变换的对称性：若)()(Ftf，则)(2)(ftF上式表明：如果函数)(tf的频谱为)(F，那么时间函数)(tF的频谱函数是)(2f。举例证明：（1）利用matlab画出信号)()(2tgtf及其幅度谱；（2）利用matlab画出信号)()(1tSatf及其幅度谱；并由实验结果验证傅立叶变换的对称特性。分析：)]1(u)1(u[)()(2tttgtf，设)()(Ftf，可知)(2)(SaF；由傅立叶变换的对称特性知:)(2)(2)(2)(2gftSatF，由门函数是偶函数以及傅立叶逆变换的线性性质，得：沈阳理工大学课程设计专用纸9沈阳理工大学)()()(2)()(211gFtSatFtf说明：在matlab中sinc(t)=tttc)sin()(sin，所以)(sinsin)(tctttSa。对称性验证编码如下：N=3001;t=linspace(-15,15,N);f=pi*[heaviside(t+1)-heaviside(t-1)];dt=30/(N-1);M=500;w=linspace(-5*pi,5*pi,M);F=f*exp(-j*t'*w)*dt;subplot(2,2,1),plot(t,f);axis([-2,2,-1,4]);xlabel('t');ylabel('f(t)');subplot(2,2,2),plot(w,real(F));axis([-20,20,-3,7]);xlabel('w');ylabel('F(w)=F[f(t)]');f1=sinc(t/pi);F1=f1*exp(-j*t'*w)*dt;沈阳理工大学课程设计专用纸10沈阳理工大学subplot(2,2,3),plot(t,f1);xlabel('t');ylabel('f1(t)=F(t)/2*pi');subplot(2,2,4),plot(w,real(F1));axis([-2,2,-1,4]);xlabel('w');ylabel('F1(w)=F[f1(t)]=f(w)');程序运行结果如下：图1沈阳理工大学课程设计专用纸11沈阳理工大学2.3傅里叶变换的时移特性验证编程设计及实现傅立叶变换的时移性质：若，则00)(0)()()(tjtjeFeFttf结论:)(tf延时（或超前）0t后，其对应的幅度谱保持不变，但相位谱中一切频率分量的相位均滞后（或超前）0t。举例证明：（1）用matlab画)(u21)(2tetft及频谱（幅度谱及相位谱)程序代码如下：N=256;t=linspace(-2,2,N);f=1/2*exp(-2*t).*heaviside(t);dt=4/(N-1);M=401;w=linspace(-2*pi,2*pi,M);F=f*exp(-j*t'*w)*dt;F1=abs(F);P1=angle(F);subplot(3,1,1);plot(t,f);gridonxlabel('t');ylabel('f(t)');title('f(t)')subplot(3,1,2);plot(w,F1);gridon沈阳理工大学课程设计专用纸12沈阳理工大学xlabel('w');ylabel('abs(F(w))');subplot(3,1,3);plot(w,P1);gridonxlabel('w');ylabel('angle(F(w))')（2）用matlab画)5.0(1tff及频谱（幅度谱及相位谱）。程序代码如下：N=256;t=linspace(-2,2,N);f=1/2*exp(-2*t).*heaviside(t);f1=1/2*exp(-2*(t-0.5)).*heaviside(t-0.5);f(t-0.3)dt=4/(N-1);M=401;w=linspace(-2*pi,2*pi,M);F=f*exp(-j*t'*w)*dt;F1=f1*exp(-j*t'*w)*dt;subplot(3,1,1);plot(t,f,t,f1,'r'),gridonxlabel('t');ylabel('f'),title('f(t),f(t-0.5)')subplot(3,1,2);沈阳理工大学课程设计专用纸13沈阳理工大学plot(w,abs(F),w,abs(F1),'r'),gridonxlabel('w');ylabel('f(t)和f(t-0.5)幅度谱');subplot(3,1,3);plot(w,angle(F),w,angle(F1),'r'),gridonxlabel('w');ylabel('f(t)和f(t-0.5)相位谱')程序运行结果如下所示：)(tf及其频谱图2沈阳理工大学课程设计专用纸14沈阳理工大学图33.总结通过本次综合实践让我们在学习“信号与系统”课程的同时，掌握MATLAB的应用，对MATLA