短时傅里叶变换(STFT)及其逆变换实验报告

短时傅里叶变换（STFT）及其逆变换实验报告

1.实验目的

1、掌握STFT及其逆变换的matlab程序。

2、能够在STFT的基础上编写matlab函数实现低通滤波。

3、分析滤波前后的效果变化和使用不同的语音质量评估法对恢复的语音信

号的效果。

2.实验概述

基于STFT的程序，编写matlab的函数，实现低通滤波。测试2个48k的music和2个8k的speech的音频通过不同低通截止频率，分析低通滤波前后的音频频谱图的变化，分别使用主观语音质量评估法和客观语音质量评估法对恢复的语音测试得分，并分析得到的结果。

3.1实验环境和配置

本次实验在matlab2012a下运行，在matlab工具箱中有STFT变换的函数，STFT的逆变换调用的外部函数，另外需要给出主观和客观语音质量评估法的程序。

3.2实验步骤

1、在matlab工具箱中找到STFT变换的spectrogram.m函数,熟悉程序，尤其分析输入输出的变量，然后了解STFT逆变换的overlapadd.m函数。

2、通过调用STFT和逆变换的函数编写一个对语音信号进行STFT及其逆变换的matlab程序，并能够显示原始语音和恢复出的语音信号的波形。

3、在对语音信号进行STFT和逆变换的matlab程序的基础上编写低通滤波的函数，并且能够显示信号在低通滤波前后的频谱图。

4、在matlab程序中采用fwSNRseg对所有的语音进行客观语音质量评估，对48k的music采用PqevalAudio，8k的speech采用PESQ进行主观质量评估。

5、对2个48k的music和2个8k的speech的音频在1k、2k、4k、8k、12k 的低通截止频率下进行测试，观察滤波前后频谱图的变化，记录在不同截止频率

下音频的主客观语音质量评估的得分，并分析得到的结果。

3.3实验结果

分别对2个48k的music和2个8k的speech语音，在1k、2k、4k、8k、12k 的低通截至频率的条件下恢复语音信号，对于语音检测的客观语音质量评估法均采用fwSNRseg，主观语音质量评估法48k的music采用PqevalAudio，8k的speech 采用PESQ。实验表格如下：

3.4实验总结

总结结论

1、对STFT及其逆变换有了比较深入的了解，尤其是对STFT的过程有了更加深入的了解。

2、在STFT的基础上编写了简单的滤波器的matlab函数，并能有很好的滤波效果。

3、fwSNRseg得分最高为35db,并且越高越好，PESQ的得分在1.0和4.5之间，同样越高越好，PqevalAudio的得分小于0，并且越接近0越好。

注：提交该报告的同时，提交测试所用的用例，必要时可给结果文件，以便可以重现实验。

阿贝成像原理实验报告

佛山科学技术学院 实验报告 课程名称近代物理实验实验项目阿贝成像原理和空间滤波 专业班级 10物师姓名邓新炬学号 02 仪器组号 指导教师朱星成绩日期 2013年月日

2、关于阿贝成像原理 成像的这两个步骤本质上就是两次傅里叶变换。第一步把物面光场的空间分布()y x g ,变为频谱面上空间频率分布() y x f f G ,，第二步则是再作一次变换，又将() y x f f G ,还原到空间分布()y x g ,。 3、空间滤波 空间函数变为频谱函数，再变回到空间函数（忽略放大率）。显然如果我们在频谱面（即透镜的后焦面）上放一些不同结构的光阑，以提取（或摒弃）某些频段的物信息，则必然使像面上的图像发生相应的变化，这样的图像处理称为空间滤波，频谱面上这种光阑称为滤波器。滤波器使频谱面上一个或一部分频率分量通过，而挡住其它频率分量，从而改变了像面上图像的频率成分。例如光轴上的圆孔光栏可以作为一个低通滤波器，而圆屏就可以用作为高通滤波器。 四 实验步骤 1、实验光路调节 在光具座上将小圆孔光阑靠近激光管的输出端，上下左右调节激光管，使激光束能穿过小孔；然后移远小孔，如光束偏离光阑，调节激光管的仰俯，再使激光能穿过小孔，重新将光阑移近，反复调节，直至小孔光阑在光具座上平移时，激光束能通过小孔光阑。 2、阿贝成像原理实验 如实验光路图在物平面上放上一维光栅，用激光器发出的细锐光束垂直照到光栅上，用一短焦距薄透镜（6~10cm ）组装一个放大的成像系统，调节透镜位置，使光栅狭缝清晰地成像在像平面屏上，那么在频谱面上的衍射点如图所示。在频谱面上放上可调狭缝或滤波模板，使通过的衍射点如下图所示：（a ）全部；（b ）零级；（c ）零和±1级；分别记录图片信息。 3、阿贝一波特实验（方向滤波） （1）光路不变，将一维光栅的物换成二维正交光栅，在频谱面上可以观察到二维分立的光点阵（频谱），像面上可以看到放大了的正交光栅像，测出像面上的网格间距。 （2）在频谱面放上可旋转狭缝光阑（方向滤波器），在下述情况：（a ）只让光轴上水平的一行频谱分量通过；（b ）只让光轴上垂直的一行频谱分量通过；（c ）只让光轴上45°的一行频谱分量通过。记录像面上的图像变化、像面上条纹间距，并做出适当的解释。 五 实验数据和数据处理 1. 1解释阿贝成像实验

MAtlab傅里叶变换实验报告

班级信工142 学号 22 姓名何岩实验组别实验日期室温报告日期成绩报告内容：(目的和要求，原理，步骤，数据，计算，小结等) 1.求信号的离散时间傅立叶变换并分析其周期性和对称性； 给定正弦信号x(t)=2*cos(2*pi*10*t),fs=100HZ,求其DTFT。 (a)代码： f=10;T=1/f;w=-10:0.2:10; t1=0:0.0001:1;t2=0:0.01:1; n1=-2;n2=8;n0=0;n=n1:0.01:n2; x5=[n>=0.01]; x1=2*cos(2*f*pi*t1); x2=2*cos(2*f*pi*t2); x3=(exp(-j).^(t2'*w)); x4=x2*x3; subplot(2,2,1);plot(t1,x1); axis([0 1 1.1*min(x2) 1.1*max(x2)]); xlabel('x(n)');ylabel('x(n)'); title('原信号x1'); xlabel('t');ylabel('x1'); subplot(2,2,3);stem(t2,x2); axis([0 1 1.1*min(x2) 1.1*max(x2)]); title('原信号采样结果x2'); xlabel('t');ylabel('x2'); subplot(2,2,2);stem(n,x5); axis([0 1 1.1*min(x5) 1.1*max(x5)]); xlabel('n');ylabel('x2'); title('采样函数x2'); subplot(2,2,4);stem(t2,x4); axis([0 1 -0.2+1.1*min(x4) 1.1*max(x4)]); xlabel('t');ylabel('x4'); title('DTFT结果x4'); （b）结果： 2.用以下两个有限长序列来验证DTFT的线性、卷积和共轭特性； (n) x1(n)=[1 2 3 4 5 6 7 8 9 10 11 12];x2(n)=R 10 (1)线性:（a）代码： w=linspace(-8,8,10000); nx1=[0:11]; nx2=[0:9]; x1=[1 2 3 4 5 6 7 8 9 10 11 12];

傅里叶光学实验

傅里叶光学的空间频谱与空间滤波实验11系09级姓名张世杰日期2011年3月30日学号PB09210044 实验目的： 1.了解傅里叶光学中基本概念，如空间频率，空间频谱，空间滤波和卷积 2.理解透镜成像的物理过程 3.通过阿贝尔成像原理，了解透镜孔径对分辨率的影响 实验原理： 一、基本概念 频谱面：透镜的后焦面 空间函数：实质即光波照明图形时从图形反射或透射出来的光波可用空间两维复变函数 空间频谱：一个复变函数f(x,y)的傅立叶变换为 ??+ ) exp[ , F)] ( ( (π , u ) { , ( )} v =dxdy vy ? = f ux - y x 2i f x y F(u,v)叫作f(x,y)的变换函数或频谱函数 空间滤波：在频谱面上放一些光栅以提取某些频段的物信息的过程 滤波器：频谱面上的光阑 二、阿贝尔成像原理 本质就是经过两次傅里叶变换，先是使单色平行光照在光栅上，经衍射分解成不同方向的很多束平行光，经过透镜分别在后焦面上形成点阵，然后代表不同空间频率的光束又在向面上复合而成像。 需要提及的是，由于透镜的大小有限，总有一部分衍射角度大的高频成分不 能进入到透镜而被丢弃了，因此像平面上总是可能会丢失一些高频的信息，即在 透镜的后焦平面上得到的不是物函数的严格的傅立叶变换（频谱），不过只有一 个位相因子的差别，对于一般情况的滤波处理可以不考虑。这个光路的优点是光 路简单，而且可以得到很大的像以便于观察。

三、空间滤波器 在频谱面上放置特殊的光阑，以滤去特定的光信号(1)单透镜系统 (2)双透镜系统 (3)三透镜系统

四、空间滤波器的种类 a ．低通滤波：在频谱面上放如图2.4-3(1）所示的光阑，只允许位于频谱面中心及附近的低频分量通过，可以滤掉高频噪音。 b ．高通滤波：在频谱面上放如图2.4-3(2)所示的光阑，它阻挡低频分量而让高频分量通过，可以实现图像的衬度反转或边缘增强。 c ． 带通滤波：在频谱面上放如图2.4-3(3)所示的光阑，它只允许特定区域的频谱通过，可以去除随机噪音。 d ．方向滤波：在频谱面上放如图2.4-3(4)或（5）所示的光阑，它阻挡或允许特定方向上的频谱分量通过，可以突出图像的方向特征。 以上滤波光阑因透光部分是完全透光，不透光部分是将光全 部挡掉，所以称作“二元振幅滤波器”。还有各种其它形式的滤波器，如：“振幅 滤波器”、“相位滤波器”和“复数滤波器”等。 e ．相幅滤波器：是将位相转变为振幅的滤波器,它的重要应用就是把”位相物体”显现出来，所谓位相物体是指那些只有空间的位相结构而透明度却一样的透明物体。如生物切片、油膜、热塑等，它们只改变入射光的位相而不影响其振幅。所以人眼不能直接看到透明体中的位相分布也就是它们的形状和结构，利用相幅转换技术就能使人眼看到透明体的形状和结构，从而扩展了人眼的视觉功能。 图 3 图2.4-3 各种形式的空间滤波器

MAtlab 傅里叶变换 实验报告

陕西科技大学实验报告 班级信工142 学号22 姓名何岩实验组别实验日期__________ 室温_____________ 报告日期________________ 成绩报告内容：(目的和要求，原理，步骤，数据，计算，小结等) 1.求信号的离散时间傅立叶变换并分析其周期性和对称性； 给定正弦信号x(t)=2*cos(2*pi*10*t),fs=100HZ, 求其DTFT (a)代码： f=10；T=1/f；w=-10:0.2:10; t1=0:0.0001:1;t2=0:0.01:1; n1=-2; n2=8; n0=0; n=n 1:0.01: n2; x5=[ n>=0.01]; x1=2*cos(2*f*pi*t1); x2=2*cos(2*f*pi*t2); x3=(exp(-j)4(t2'*w)); x4=x2*x3; subplot(2,2,1);plot(t1,x1); axis([0 1 1.1*mi n(x2) 1.1*max(x2)]); xlabel('x( n)');ylabel('x( n)'); title('原信号x1'); xlabel('t');ylabel('x1'); subplot(2,2,3);stem(t2,x2); axis([0 1 1.1*mi n(x2) 1.1*max(x2)]); title(' 原信号采样结果x2'); xlabel('t');ylabel('x2'); subplot(2,2,2);stem( n, x5); axis([0 1 1.1*mi n(x5) 1.1*max(x5)]); xlabel(' n');ylabel('x2'); title(' 采样函数x2'); subplot(2,2,4);stem(t2,x4); axis([0 1 -0.2+1.1*mi n(x4) 1.1*max(x4)]); xlabel('t');ylabel('x4'); title('DTFT 结果x4'); (b)结果：

傅里叶光学实验报告

实验原理：（略） 实验仪器： 光具座、氦氖激光器、白色像屏、作为物的一维、二维光栅、白色像屏、傅立叶透镜、小透镜 实验内容与数据分析 1．测小透镜的焦距f 1 （付里叶透镜f 2=45.0CM ） 光路：激光器→望远镜（倒置）（出射应是平行光）→小透镜→屏 操作及测量方法：打开氦氖激光器，在光具座上依次放上扩束镜，小透镜和光屏，调节各光学元件的相对位置是激光沿其主轴方向射入，将小透镜固定，调节光屏的前后位置，观察光斑的会聚情况，当屏上亮斑达到最小时，即屏处于小透镜的焦点位置，测量出此时屏与小透镜的距离，即为小透镜的焦距。 112.1913.2011.67 12.3533 f cm ++= = 0.7780cm σ= = 1.320.5929 p A p t t cm μ=== 0.68P = 0.0210.00673 B p B p t k cm C μ?==?= 0.68P = 0.59cm μ== 0.68P = 1(12.350.59)f cm =± 0.68P =

2．利用弗朗和费衍射测光栅的的光栅常数 光路：激光器→光栅→屏（此光路满足远场近似） 在屏上会观察到间距相等的k 级衍射图样，用锥子扎孔或用笔描点，测出衍射图样的间距，再根据sin d k θλ=测出光栅常数d （1）利用夫琅和费衍射测一维光栅常数； 衍射图样见原始数据； 数据列表： sin || i k Lk d x λλ θ= ≈ 取第一组数据进行分析： 2105 13 43.0910******* 4.00106.810d m ----????==?? 210 523 43.0910******* 3.871014.110d m ----????==?? 2105 33 43.0910******* 3.95106.910d m ----????==?? 210 543 43.0910******* 4.191013.010 d m ----????==?? 554.00 3.87 3.95 4.19 10 4.0025104 d m m --+++= ?=? 61.3610d m σ-=? 忽略b 类不确定度：

快速傅里叶变换实验报告..

快速傅里叶变换实验报告 班级： 姓名： 学号:

快速傅里叶变换 一．实验目的 1.在理论学习的基础上，通过本实验加深对快速傅立叶变换的理解； 2.熟悉并掌握按时间抽取FFT 算法的程序； 3.了解应用FFT 进行信号频谱分析过程中可能出现的问题，例如混淆、泄漏、栅栏效应等，以便在实际中正确应用FFT 。 二．实验内容 1.仔细分析教材第六章‘时间抽取法FFT ’的算法结构，编制出相应的用FFT 进行信号分析的C 语言（或MATLAB 语言）程序； 2.用FFT 程序分析正弦信号 ()sin(2)[()(*)],(0)1y t f t u t u t N T t u π=---∞<<+∞=设 分别在以下情况进行分析并讨论所得的结果： a ） 信号频率f ＝50Hz ，采样点数N=32,采样间隔T=0.000625s b ） 信号频率f ＝50Hz ，采样点数N=32,采样间隔T=0.005s c ） 信号频率f ＝50Hz ，采样点数N=32,采样间隔T=0.0046875s d ） 信号频率f ＝50Hz ，采样点数N=32,采样间隔T=0.004s e ） 信号频率 f ＝50Hz ，采样点数N=64,采样间隔T=0.000625s f ） 信号频率f ＝250Hz ，采样点数N=32,采样间隔T=0.005s g ） 将c ） 信号后补32个0，做64点FFT 三．实验要求 1.记录下实验内容中各种情况下的X (k)值，做出频谱图并深入讨论结果，说明参数的变化对信号频谱产生哪些影响。频谱只做模特性，模的最大值＝1，全部归一化；

2.打印出用C 语言（或MATLAB 语言）编写的FFT 源程序，并且在每一小段处加上详细的注释说明； 3.用C 语言（或MATLAB 语言）编写FFT 程序时，要求采用人机界面形式： N , T , f 变量均由键盘输入，补零或不补零要求设置一开关来选择。 四．实验分析 对于本实验进行快速傅里叶变换，依次需要对信号进行采样，补零（要求补零时），码位倒置，蝶形运算，归一化处理并作图。 此外，本实验要求采用人机界面形式，N,T,F 变量由键盘输入，补零或不补零设置一开关来选择。 1.采样 本实验进行FFT 运算，给出的是正弦信号，需要先对信号进行采样，得到有限 长序列()n x ， N n ...... 2,1,0= Matlab 实现： t=0:T:T*(N-1); x=sin(2*pi*f*t); 2.补零 根据实验要求确定补零与否，可以用if 语句做判断，若为1，再输入补零个数， 并将补的零放到采样得到的序列的后面组成新的序列，此时新的序列的元素个数等于原采样点个数加上补零个数，并将新的序列个数赋值给N 。 Matlab 实现： a=input('是否增加零点? 是请输入1 否请输入0\n'); if (a) ZeroNum=input('请输入增加零点的个数:\n'); else ZeroNum=0; end if (a) x=[x zeros(1, ZeroNum)];%%指令zeros(a,b)生成a 行b 列全0矩阵，在单行矩阵x 后补充0 end N=N+ZeroNum; 3.码位倒置 本实验做FFT 变换的级数为M ，N M 2log =

光学仪器实验报告

常用光电仪器原理及使用 实验报告 班级：11级光信息1班 姓名：姜萌萌 学号：110104060016 指导老师：李炳新

数字存储示波器 一、实验目的 1、熟悉数字存储示波器的使用方法； 2、测量数字存储示波器产生方波的上升时间； 二、实验仪器 数字存储示波器 三、实验步骤 1、产生方波波形 ⑴、打开示波器电源阅读探头警告，然后按下OK。按下“DEFAULT SETUP”按钮，默认的电压探头衰减选项是10X。 ⑵、在P2200探头上将开关设定到10X并将探头连接到示波器的通道1上，然后向右转动将探头锁定到位，将探头端部和基线导线连接到“PROBE COMP”终端上。 ⑶、按下“AUTOSET”按钮，在数秒钟内，看到频率为1KHz 电压为5V峰峰值得方波。按两次CH1BNC按钮删除通道1，

按下CH2BNC按钮显示通道2，重复第二步和第三步。 2、自动测量 ⑴、按下“MUASURE”按钮，查看测量菜单。 ⑵、按下顶部的选项按钮，显示“测量1菜单”。 ⑶、按下“类型”“频率”“值”读书将显示测量结果级更新信息。 ⑷、按下“后退”选项按钮。 ⑸、按下顶部第二个选项按钮；显示“测量2菜单”。 ⑹、按下“类型”“周期”“值”读数将显示测量结果与更新信息。 ⑺、按下“后退”选项按钮。 ⑻、按下中间选项按钮；显示“测量3菜单”。 ⑼、按下“类型”“峰-峰值”“值”读数将显示测量结果与更新信息。 ⑽、按下“后退”选项按钮。 ⑾、按下底部倒数第二个按钮；显示“测量4菜单”。⑿、按下“类型”“上升时间”“值”读数将显示测量结果与更新信息。

LCR测试仪 一、实验目的 1、熟悉LCR测试仪的使用方法； 2、了解LCR测试仪的工作原理； 3、精确测量一些电阻，电感，电容的值； 二、实验仪器 LCR测试仪，电阻，电容，电感等元件 三、LCR测试原理 根据待测元器件实际使用的条件和组合上的差别，LCR 测量仪有两种检测模式，串联模式和并联模式。串联模式以检测元器件Z为基础，并联模式以检测元器件的导纳Y为基础，当用户将测出流过待测元件的电流I，数字电压表将测出待测元件两端的电压V，数字鉴相器将测出电压V和电流I 之间的相位角 。检测结果被储存在仪器内部微型计算机的

信号与系统实验报告3实验3 傅里叶变换及其性质

信息工程学院实验报告 课程名称： 实验项目名称：实验3 傅里叶变换及其性质 实验时间：2015/11/17 班级：通信141 姓名： 学号： 一、实 验 目 的: 学会运用MATLAB 求连续时间信号的傅里叶（Fourier ）变换；学会运用MATLAB 求连续时间信号的频谱图；学会运用MATLAB 分析连续时间信号的傅里叶变换的性质。 二、实 验 设 备 与 器 件 软件：Matlab 2008 三、实 验 原 理 3.1傅里叶变换的实现 信号()f t 的傅里叶变换定义为： ()[()]()j t F F f t f t e dt ωω∞ --∞ ==? ， 傅里叶反变换定义为：1 1()[()]()2j t f t F F f e d ωωωωπ ∞ --∞ == ? 。 信号的傅里叶变换主要包括MATLAB 符号运算和MATLAB 数值分析两种方法，下面分别加以探讨。同时，学习连续时间信号的频谱图。 3.1.1 MATLAB 符号运算求解法 MATLAB 符号数学工具箱提供了直接求解傅里叶变换与傅里叶反变换的函数fourier( )和ifourier( )。Fourier 变换的语句格式分为三种。 （1）F=fourier(f)：它是符号函数f 的Fourier 变换，默认返回是关于ω的函数。 （2）F=fourier(f,v)：它返回函数F 是关于符号对象v 的函数，而不是默认的 ω，即 ()()j v t F v f t e d t ∞ --∞ =? 。 （3）F=fourier(f,u,v)：是对关于u 的函数f 进行变换，返回函数F 是关于v 的函数，即 ()()jvu F v f t e du ∞ --∞ =?。 傅里叶反变换的语句格式也分为三种。 （1）f=ifourier(F)：它是符号函数F 的Fourier 反变换，独立变量默认为ω，默认返回是关于x 的函数。 （2）f=ifourier(F,u)：它返回函数f 是u 的函数，而不是默认的x 。 （3）f=ifourier(F,u,v)：是对关于v 的函数F 进行反变换，返回关于u 的函数f 。

实验三傅里叶变换及其性质

信息工程学院实验报告 课程名称：信号与系统 实验项目名称：实验3 傅里叶变换及其性质实验时间：2013-11-29 班级： 姓名： 学号： 一、实验目的: 1、学会运用MATLAB 求连续时间信号的傅里叶（Fourier ）变换； 2、学会运用MATLAB 求连续时间信号的频谱图； 3、学会运用MATLAB 分析连续时间信号的傅里叶变换的性质。 二、实验环境: 1、硬件：在windows 7 操作环境下； 2、软件：Matlab 版本7.1 三、实验原理： 3.1傅里叶变换的实现 信号()f t 的傅里叶变换定义为： ()[()]()j t F F f t f t e dt ωω∞ --∞ == ? ， 傅里叶反变换定义为：1 1 ()[()]()2j t f t F F f e d ωωωωπ ∞ --∞ == ? 。 信号的傅里叶变换主要包括MATLAB 符号运算和MATLAB 数值分析两种方法，下面分别加以探讨。同时，学习连续时间信号的频谱图。 3.1.1 MATLAB 符号运算求解法 MATLAB 符号数学工具箱提供了直接求解傅里叶变换与傅里叶反变换的函数fourier( )和ifourier( )。Fourier 变换的语句格式分为三种。 （1）F=fourier(f)：它是符号函数f 的Fourier 变换，默认返回是关于ω的函数。 （2）F=fourier(f,v)：它返回函数F 是关于符号对象v 的函数，而不是默认的ω，即()()jvt F v f t e dt ∞ --∞ = ? 。 （3）F=fourier(f,u,v)：是对关于u 的函数f 进行变换，返回函数F 是关于v 的函数，即 ()()jvu F v f t e du ∞ --∞ =? 。 傅里叶反变换的语句格式也分为三种。 （1）f=ifourier(F)：它是符号函数F 的Fourier 反变换，独立变量默认为ω，默认返回是关于x 的函数。 （2）f=ifourier(F,u)：它返回函数f 是u 的函数，而不是默认的x 。 （3）f=ifourier(F,u,v)：是对关于v 的函数 F 进行反变换，返回关于u 的函数f 。

快速傅里叶变换实验报告

快速傅里叶变换实验报告

————————————————————————————————作者：————————————————————————————————日期: ?

快速傅里叶变换实验报告 机械34班 刘攀 2０1３0105５8 一、 基本信号（函数)的FF Ｔ变换 1. 000()sin()sin 2cos36x t t t t π ωωω=+++ 1) 采样频率08s f f =,截断长度N ＝１６； 取02ωπ=ｒad/ｓ，则0f =1Hz ，s f =8Ｈz ，频率分辨率 f ?=s f f N ?= =0.5Hz 。 最高频率c f ＝30f =３Hz ，s f ＞2c f ,故满足采样定理,不会发生混叠现象。 截断长度02T T =,整周期截取,不会发生栅栏效应。理论上有一定的泄漏,但在整周期 截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下:

幅值误差0A ?=,相位误差0??=。 2) 采样频率08s f f =，截断长度N=3２; 取02ωπ=rad/s ，则0f =１Hz,s f ＝８Hz ，频率分辨率f ?=s f f N ?==０.2５Hz 。 最高频率c f =30f =3H ｚ,s f >２c f ,故满足采样定理,不会发生混叠现象。 截断长度04T T =,整周期截取,不会发生栅栏效应。理论上有一定的泄漏，但在整周期 截取的情况下,旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下：

幅值误差0A ?=，相位误差0??=。 2. 00()sin()sin116x t t t π ωω=++ 1) 采样频率08s f f =,截断长度N=16； 取02ωπ=ra ｄ/s,则0f =１Hz ，s f =８Hz,频率分辨率f ?=s f f N ?==0．5H ｚ。 最高频率c f =110f =11H ｚ,s f <２c f ，故不满足采样定理，会发生混叠现象。 截断长度02T T =,整周期截取，不会发生栅栏效应。理论上有一定的泄漏,但在整周期 截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图：

大学物理实验傅里叶分析实验报告

脉搏、语音及图像信号的傅里叶分析 一、实验简介 任何波形的周期信号均可用傅里叶级数来表示。傅里叶级数的各项代表了不同频率的正弦或余弦信号，即任何波形的周期信号都可以看作是这些信号（谐波）的叠加。利用不同的方法，可以从周期信号中分解出它的各次谐波的幅值和相位。也可依据信号的傅里叶级数表达式，将各次谐波按表达式的要求叠加得到所期望的信号。 二、实验目的 1、了解常用周期信号的傅里叶级数表示。 2、了解周期脉搏信号、语音信号及图像信号的傅里叶分析过程 3、理解体会傅里叶分析的理论及现实意义 三、实验仪器 脉搏语音实验仪器，数字信号发生器，示波器 四、实验原理 1、周期信号傅里叶分析的数学基础 任意一个周期为T 的函数f(t)都可以表示为傅里叶级数： 00010000 1 ()(cos sin ) 21()() 1 ()cos()()1 ()sin()()n n n n n f t a a n t b n t a f t d t a f t n t d t b f t n t d t π π π ππ πωωωωπ ωωωπ ωωωπ ∞ =-- - =++=== ∑??? 其中0ω为角频率，称为基频，0a 为常数，n a 和n b 称为第n 次谐波的幅值。任何

周期性非简谐交变信号均可用上述傅里叶级数进行展开，即分解为一系列不同次谐波的叠加。 对于如图1所示的方波，一个周期内的函数表达式为: (0t<)2() (-t 0) 2 h f t h ππ? ≤??=? ?-≤

傅里叶变换实验报告

南昌大学实验报告 学生姓名：学号：6100209228 班级：电子093班 实验类型：□验证□综合■设计□创新实验日期：2011-04-8 实验成绩： 傅里叶变换 （一）实验目的 1、掌握对不同的函数进行傅里叶变换的程序编写； 2、熟悉生成联系周期信号的方法； 3、练习matlab编程。 (二) 实验内容 1.请编写函数F=fsana(t,f,,N),计算周期信号f的前N个指数形式的傅立叶级数系数，t表示f对应的抽样时间（均为一个周期）；再编写函数f=fssyn（F，t），由傅立叶级数系数F合成抽样时间t对应的函数。设计信号验证这两个是否正确。 定义F=fsana（t，f，N）。 function F=fsana(t,f,N) omg1=2*pi/(max(t)-min(t)); k=[0:N]'; F=1/length(t)*exp(-j*kron(k*omg1,t.'))*f 定义f=fssyn（F,t） function f=fssyn(F,t) omg1=2*pi/(max(t)-min(t)); N=floor(length(F)/2); k=[0:N]; f=exp(j*kron(t,k*omg1))*F; 运行所定义的函数 T1=2*pi; %一个周期时域范围 N1=300; %时域抽样点数

t=linspace(0,T1-T1/N1,N1)'; %生成抽样时间点 f=cos(t); %生成抽样函数值 subplot(2,2,1) plot(t,f); title ('原函数') N=10; F1=fsana(t,f,N); %调用fsana函数求解前N项傅立叶级数系数 subplot(2,2,2) stem(abs(F1),'s'); %绘制离散的幅度曲线 title('前N项傅立叶级数系数幅度曲线'); f2=fssyn(F1,t); %调用fssyn函数求原时域函数 subplot(2,2,3) plot(t,f2,'k'); title('傅立叶逆变换后时域函数'); 运行结果

傅里叶变换光学

中山大学光信息专业实验报告：傅里叶光学变换系统 一、实验目的和内容 1、了解透镜对入射波前的相位调制原理。 2、加深对透镜复振幅、传递函数、透过率等参量的物理意义的认识。 3、观察透镜的傅氏变换（FT ）图像，观察4f 系统的反傅氏变换（IFT ）图像，并进行比较。 4、在4f 系统的变换平面（T ）插入各种空间滤波器，观察各种试件相应的频谱处理图像。 二、实验原理 1、透镜的FT 性质及常用函数与图形的关学频谱分析 透镜由于本身厚度的不同，使得入射光在通过透镜时， 图1 点的厚度。设原复振幅分布为(,)L U x y 的光通过透镜后，幅分布受到透镜的位相调制，附加了一个位相因子(,)x y ?为(,)L U x y '： 图1 (,)(,)exp[(,)]L L U x y U x y j x y ?'= （1） 若对于任意一点（x ，y ）透镜的厚度为(,)D x y ，透镜的中心厚度为0D 。光线由该点通过透镜时在透镜中的距离为(,)D x y ，空气空的距离为0D －(,)D x y ，透镜折射率为n ，则该点的总的位相差为： 00(,)[(,)](,)(1)(,)x y k D D x y knD x y kD k n D x y ?=-+=+- （2） （2）中的k ＝2π/λ，为入射光波波数。 用位相延迟因子(,)t x y 来表示即为： 0(,)exp()exp[(1)(,)]t x y jkD jk n D x y =- （3） 由此可见只要知道透镜的厚度函数(,)D x y 就可得出其相位调制。在球面镜傍轴区域，用抛物面近似球面，可以得到球面透镜的厚度函数为：

快速傅里叶变换实验报告

快速傅里叶变换实验报告 快速傅里叶变换实验报告 机械34班刘攀 2019010558 一、基本信号（函数）的FFT变换 1. x(t)=sin(ω0t+)+sin2ω0t+cos3ω0t 6 1) 采样频率fs=8f0，截断长度N=16； 取ω0=2πrad/s，则f0=1Hz，fs=8Hz，频率分辨率?f=?f=fs=0.5Hz。 Nπ最高频率fc=3f0=3Hz，fs>2fc，故满足采样定理，不会发生混叠现象。截断长度T=2T0，整周期截取，不会发生栅栏效应。理论上有一定的泄漏，但在整周期截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下： 幅值误差?A=0，相位误差??=0。 2) 采样频率fs=8f0，截断长度N=32； 取ω0=2πrad/s，则f0=1Hz，fs=8Hz，频率分辨率?f=?f=fs=0.25Hz。 N最高频率fc=3f0=3Hz，fs>2fc，故满足采样定理，不会发生混叠现象。截断长度T=4T0，整周期截取，不会发生栅栏效应。理论上有一定的泄漏，但在整周期截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下： 幅值误差?A=0，相位误差??=0。 2. x(t)=sin(ω0t+π 6)+sin11ω0t 1) 采样频率fs=8f0，截断长度N=16； 取ω0=2πrad/s，则f0=1Hz，fs=8Hz，频率分辨率?f=?f=fs=0.5Hz。 N最高频率 fc=11f0=11Hz，fs 漏，但在整周期截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图：

由上图可以看出，并未体现出11f0的成分，说明波形出现混叠失真。为了消除混叠 现象，应加大采样频率，使之大于等于 22Hz。 f0处的幅值误差?A=0，11f0处由于出现 了混叠现象，幅值误差没有意义；相位误差??=0。 2) 采样频率fs=32f0，截断长度N=32； 取ω0=2πrad/s，则f0=1Hz，fs=32Hz，频率分辨率?f=?f=fs=1Hz。 N最高频率 fc=11f0=11Hz，fs>2fc，故满足采样定理，不会发生混叠现象。 漏，但在整周期截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图： 该频谱图体现出了f0和11f0的成分，说明未失真，且幅值均为1,。幅值误差?A=0，相位误差??=0。 3. x(t)=0t 1) 采样频率fs=8f0，截断长度N=16； 取ω0=2πrad/s，则f0=1Hz，fs=8Hz，频率分辨率?f=?f=fs=0.5Hz。 N最高频率f cf 0Hz，fs>2fc，故满足采样定理，不会发生混叠现象。 频谱图： 在忽略旁瓣信号的情况下，可近似认为： x(t)≈0.9098cos(3ω0t+56.9520?) 故幅值误差?A=0.9096-1=-0.0904，相位误差??=56.9520?。 2) 采样频率fs=32f0，截断长度N=32； 取ω0=2πrad/s，则f0=1Hz，fs=32Hz，频率分辨率?f=?f=fs=1Hz。N最高频率f cf 0Hz，fs>2fc，故满足采样定理，不会发生混叠现象。 频谱图： 在忽略旁瓣信号的情况下，可近似认为：

MAtlab傅里叶变换实验报告

M A t l a b傅里叶变换实 验报告 文档编制序号：[KKIDT-LLE0828-LLETD298-POI08]

班级信工142 学号 22 姓名何岩实验组别 实验日期室温报告日期成绩 报告内容：(目的和要求，原理，步骤，数据，计算，小结等) 1.求信号的离散时间傅立叶变换并分析其周期性和对称性； 给定正弦信号x(t)=2*cos(2*pi*10*t),fs=100HZ,求其DTFT。(a)代码： f=10;T=1/f;w=-10::10; t1=0::1;t2=0::1; n1=-2;n2=8;n0=0;n=n1::n2; x5=[n>=]; x1=2*cos(2*f*pi*t1); x2=2*cos(2*f*pi*t2); x3=(exp(-j).^(t2'*w)); x4=x2*x3; subplot(2,2,1);plot(t1,x1); axis([0 1 *min(x2) *max(x2)]); xlabel('x(n)');ylabel('x(n)'); title('原信号x1'); xlabel('t');ylabel('x1'); subplot(2,2,3);stem(t2,x2); axis([0 1 *min(x2) *max(x2)]); title('原信号采样结果x2'); xlabel('t');ylabel('x2'); subplot(2,2,2);stem(n,x5); axis([0 1 *min(x5) *max(x5)]); xlabel('n');ylabel('x2'); title('采样函数x2'); subplot(2,2,4);stem(t2,x4); axis([0 1 +*min(x4) *max(x4)]); xlabel('t');ylabel('x4'); title('DTFT结果x4'); （b）结果： 2.用以下两个有限长序列来验证DTFT的线性、卷积和共轭特性；

傅立叶光学实验报告

实验报告 陈杨 PB05210097 物理二班 实验题目: 傅里叶光学实验 实验目的: 加深对傅里叶光学中的一些基本概念和理论的理解，验证阿贝成像理论，理解透镜成像过程，掌握光学信息处理的实质，进一步了解透镜孔径对分辨率的影响。 实验原理： 1.傅里叶光学变换 二维傅里叶变换为：??+-=?=dxdy vy ux i y x f v u F )](2exp[),()}y ,x (f {),(π ( 1 ) 1()[(,)]x y g x F a f f -=， ''x y x f f y f f λλ??=????????=???? 复杂的二维傅里叶变换可以用透镜来实现，叫光学傅里叶变换。 2.阿贝成像原理 由于物面与透镜的前焦平面不重合，根据傅立叶光 学的理论可以知换（频谱），不过只有一个位相因 子的差别，对于一般情况的滤波处理可以不考虑。这个光路的优道在透镜的后焦平面上得到的不是物函数的严格的傅立叶变点是光路简单，是显微镜物镜成像的情况—可以得到很大的象以便于观察，这正是阿贝当时要改进显微镜的分辨本领时所用的光路。 3.空间滤波

根据以上讨论：透镜的成像过程可看作是两次傅里叶变换，即从空间函数(,)g x y 变为频谱函数(,)x y a f f ，再变回到空间函数(,)g x y ，如果在频谱面上放一不同结构的光阑，以提取某些频段的信息，则必然使像上发生相应的变化，这样的图像处理称空间滤波。 实验内容: 1.测小透镜的焦距f1 （付里叶透镜f2=）. 光路：直角三棱镜→望远镜（倒置）（出射应是平行光）→小透镜→屏。（思考：如何测焦距） 夫琅和费衍射： 光路：直角三棱镜→光栅→墙上布屏（此光路满足远场近似） （1）利用夫琅和费衍射测一维光栅常数； 光栅方程：dsin θ=k λ 其中,k=0,±1, ±2, ±3,… 请自己选择待测量的量和求光栅常数的方法。（卷尺可向老师索要） 记录一维光栅的衍射图样、可看到哪些级记录 0级、±1级、±2级光斑的位置； （2）记录二维光栅的衍射图样. 3.观察并记录下述傅立叶频谱面上不同滤波条件的图样或特征； 光路：直角三棱镜→光栅→小透镜→滤波模板（位于空间频谱面上）→墙上屏 思考：空间频谱面在距小透镜多远处图样应是何样 （1）一维光栅：（滤波模板自制，一定要注意戴眼镜保护；可用一张纸，一根针扎空来制作，也可用其他方法）. a.滤波模板只让 0级通过;

图像的傅里叶变换实验报告

计算机科学与技术系 实验报告 专业名称计算机科学与技术 课程名称数字图像处理 项目名称Matlab语言、图像的傅里叶变换 班级 14计科2班 学号 1404011023 姓名卢爱胜 同组人员张佳佳、王世兜、张跃文 实验日期 2016.11.30

一、实验目的与要求： （简述本次实验要求达到的目的，涉及到的相关知识点，实验的具体要求。） 实验目的： 1了解图像变换的意义和手段； 2熟悉傅立叶变换的基本性质； 3熟练掌握FFT 变换方法及应用； 4通过实验了解二维频谱的分布特点； 5通过本实验掌握利用MATLAB 编程实现数字图像的傅立叶变换。 6评价人眼对图像幅频特性和相频特性的敏感度。 实验要求： 应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具，它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能，将有助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说，把时间用在学习和掌握博里叶变换上是很有必要的。 二、实验内容 （根据本次实验项目的具体任务和要求，完成相关内容，可包括：实验目的、算法原理、实验仪器、设备选型及连线图、算法描述或流程图、源代码、实验运行步骤、关键技术分析、测试数据与实验结果、其他 ） 1.傅立叶（Fourier ）变换的定义 对于二维信号，二维Fourier 变换定义为： 2()(,)(,)j ux uy F u v f x y e dxdy π∞∞ -+-∞-∞= ?? 逆变换： 2()(,)(,)j ux uy f x y F u v e dudv π∞∞ +-∞-∞= ?? 二维离散傅立叶变换为： 11 2()00 1(,)(,)i k N N j m n N N i k F m n f i k e N π---+===∑∑ 逆变换：

傅里叶光学实验

傅里叶光学的空间频谱与空间滤波实验 11系09级姓名张世杰日期2011年3月30日学号PB0实验目的： 1.了解傅里叶光学中基本概念，如空间频率，空间频谱，空间滤波 和卷积 2.理解透镜成像的物理过程 3.通过阿贝尔成像原理，了解透镜孔径对分辨率的影响 实验原理： 一、基本概念 频谱面：透镜的后焦面 空间函数：实质即光波照明图形时从图形反射或透射出来的光波可用空间两维复变函数 空间频谱：一个复变函数f(x,y)的傅立叶变换为 ??+ ) , u exp[ F)] ( ( )} ) , { ( , (π y =dxdy vy v = - ux ? x f f x y 2i F(u,v)叫作f(x,y)的变换函数或频谱函数 空间滤波：在频谱面上放一些光栅以提取某些频段的物信息的过程 滤波器：频谱面上的光阑 二、阿贝尔成像原理 本质就是经过两次傅里叶变换，先是使单色平行光照在光栅上，经衍射分解成不同方向的很多束平行光，经过透镜分别在后焦面上形成点阵，然后代表不同空间频率的光束又在向面上复合而成像。

需要提及的是，由于透镜的大小有限，总有一部分衍射角度大的高频成分不能进入到透镜而被丢弃了，因此像平面上总是可能会丢失一些高频的信息，即在透镜的后焦平面上得到的不是物函数的严格的傅立叶变换（频谱），不过只有一个位相因子的差别，对于一般情况的滤波处理可以不考虑。这个光路的优点是光路简单，而且可以得到很大的像以便于观察。 物面透镜频谱面 像面 三、空间滤波器 在频谱面上放置特殊的光阑，以滤去特定的光信号 (1)单透镜系统 (2)双透镜系统

(3)三透镜系统 四、空间滤波器的种类 a．低通滤波：在频谱面上放如图2.4-3(1）所示的光阑，只允许位于频谱面中心及附近的低频分量通过，可以滤掉高频噪音。 b．高通滤波：在频谱面上放如图 2.4-3(2)所示的光阑，它阻挡低频分量而让高频分量通过，可以实现图像的衬度反转或边缘增强。 c．带通滤波：在频谱面上放如图2.4-3(3)所示的光阑，它只允许特定区域的频谱通过，可以去除随机噪音。 d．方向滤波：在频谱面上放如图2.4-3(4)或（5）所示的光阑，

大学物理仿真实验傅里叶光学

大学物理仿真实验 ——傅里叶光学实验 实 验 报 告 姓名： 班级： 学号：

实验名称傅里叶光学实验 一、实验目的 1．学会利用光学元件观察傅立叶光学现象。 2．掌握傅立叶光学变换的原理，加深对傅立叶光学中的一些基本概念和基本理论的理解，如空间频率、空间频谱、空间滤波和卷积等。 二、实验所用仪器及使用方法 防震实验台，He-Ne激光器，扩束系统（包括显微物镜，针孔（30μm），水平移动调整器)，全反射镜，透镜及架（f=+150mm，f=+100mm），50线/mm光栅滤波器，白屏 三、实验原理 平面波Ee(x,y)入射到p平面（透过率为）在p平面后Z=0处的光场分布为：E(x,y)= Ee(x,y) 图根据惠更斯原理（Huygens’ Principle）,在p平面后任意一个平面p’处光场的分布可看成p平面上每一个点发出的球面波的组合，也就是基尔霍夫衍射积分（Kirchhoff’s diffraction integral）。 (1) 这里：＝球面波波长； n＝p平面（x,y）的法线矢量；

K＝（波数） 是位相和振幅因子； cos(n,r)是倾斜因子； 在一般的观察成像系统中，cos(n,r)1。 r=Z+，分母项中r z；（1）式可用菲涅尔衍射积分表示：（菲涅尔近似 Fresnel approximation） (2) 当z更大时，即z>>时，公式（2）进一步简化为夫琅和费衍射积分：（Fraunhofer Approximation） 这里： 位相弯曲因子。 如果用空间频率做为新的坐标有： ， 若傅立叶变换为 (4)

(3)式的傅立叶变换表示如下： E(x’,y’,z)＝F[E(x,y)]＝c 图2 空间频率和光线衍射角的关系 tg＝＝，tg＝＝ ＝，＝ 可见空间频率越高对应的衍射角也越大，当z越大时，衍射频谱也展的越宽； 由于感光片和人眼等都只能记录光的强度（也叫做功率谱），所以位相弯曲因子 (5) 理论上可以证明，如果在焦距为f的汇聚透镜的前焦面上放一振幅透过率为g(x,y)的图象作为物，并用波长为的单色平面波垂直照明图象，则在透镜后焦面上的复振幅分布就是g(x,y)的傅立叶变换，其中空间频率，与坐标， 的关系为：，。故面称为频谱面（或傅氏面，由此可见，复杂的二维傅立叶变换可以用一透镜来实现，称为光学傅立叶变换，频谱面上的光强分布，也就是物的夫琅禾费衍射图。 四、实验结果