实验二快速傅里叶变换算法实验

快速傅里叶变换算法实验

院系：电子与信息工程学院

班级：电信10-2

姓名：张玥

学号：10071201026

指导老师：郜参观

快速傅里叶变换算法实验

一、实验目的

1、加深对DFT算法原理和基本性质的理解；

2、熟悉FFT算法原理和FFT子程序的应用；

3、学习用FFT对连续信号和时域信号进行谱分析的方法，了解可能

出现的分析误差及其原因，以便在实际中正确应用FFT。

二、实验设备

计算机 CCS2.0版软件实验箱 DSP仿真器

三、基本原理

1、离散傅里叶变换DFT的定义：将时域的采样变换成频域的周期性离散函数，频域的采样也可以变换成时域的周期性离散函数，这样的变换称为离散傅里叶变换，简称DFT。

2、FFT是DFT的一种快速算法，将DFT的步运算减少为

步，极大地提高了运算速度。

3、旋转因子的变化规律。

4、蝶形运算规律。

5、基2FFT算法。

四、实验步骤

复习DFT的定义、性质和用DFT作谱分析的有关内容；

1、复习FFT算法原理与编程思想，并对照DIT-FFT运算流程图

和程序框图，了解本试验提供的FFT子程序；

2、阅读本实验所提供的样例子程序；

3、运行CCS软件，对样例程序进行跟踪，分析结果；记录

必要的参数。

4、填写实验报告。

5、提供样例程序实验操作说明。

A、实验前准备：

⑴、开关K9拨到右边，即仿真器选择连接右边的CPU：CPU2；

⑵、“A/D转换单元”的拨码开关设置；

⑶、拨码开关设置：

JP3拨码开关：

SW2拨码开关：

S23拨码开关：

⑷、检查：计算机、DSP仿真器、实验箱是否正确连接。

正确完成计算机、DSP仿真器和实验箱的连接后，系统上电。

B、实验

⑴、启动CCS2.0，project/open打开“algorithm”目录“exp01_cpu2”中子目录下“exp01.pjt”工程文件；双击“exp01.pjt”及“source”可查看各源程序；加载“exp01.out”；在中断子程序中，flag=0处设置断点；单击“run”运行程序，程序将运行至断点处停止，如图（一）所示；

图（一）

⑵、用view/graph/time/frequency打开一个图形观察窗口；设置该观察图形窗口变量及参数；采用双踪观察在起始地址分别为x和mo，长度为128的单元中数值的变化，数值类型为32位浮点数，这两段存储单元中分别存放的是经A/D转换后的输入信号和对该信号进行FFT变换的结果如图（二）所示；

图（二）

单击“animate”运行程序，或按F10运行；调整模拟信号源的频率、幅值、种类观察窗口并观察输入信号波形及其FFT变换结果，如图（三）、图（四）所示；单击“Halt”暂停程序运行，关闭窗口，本实验结束。

图（三）

图（四）五、子程序流程图：

常用函数傅里叶变换

信号与系统的基本思想：把复杂的信号用简单的信号表示，再进行研究。怎么样来分解信号？任何信号可以用Delta 函数的移位加权和表示。只有系统是线性时不变系统，才可以用单位冲激函数处理，主要讨论各个单位冲激函数移位加权的响应的叠加能得到总的响应。线性系统（齐次性，叠加定理）时不变系统对一个系统输入单位冲激函数，得到的响应为h(t).表征线性时不变系统的非常重要的东西，只要知道了系统对单位冲击函数的响应，就知道了它对任何信号的响应，因为任何信号都可以表示为单位冲激函数的移位加权和。例如：d(t)__h(t) 那么a*d(t-t0)__a*h(t-t0) -()= ()(t-)d f t f τδττ∝∝? 的响应为-y()=()(-)t f h t d τττ∝ ∝ ? 记为y(t)=f(t)*h(t),称为f(t)和h(t)的卷积总结为两点：对于现行时不变系统，任何信号可以用单位冲激信号的移位加权和表示，任何信号的响应可以用输入函数和单位冲激函数响应的卷积来表示连续时间信号和系统的频域分析时域分析的重点是把信号分解为单位冲激函数的移位加权和，只讨论系统对单位冲激函数的响应。而频域的分析是把信号分解为各种不同频率的正弦函数的加权和，只讨论系统对sinwt 的响应。都是把信号分解为大量单一信号的组合。

周期函数可以展开为傅里叶级数，将矩形脉冲展开成傅里叶级数，得到傅里叶级数的系数 n A sin F = T x x τ 其中0=2 nw x τ。取样函数sin ()=x S a x 。产生一种震荡，0点的值最大，然后渐渐衰减直至0 第一：对于傅里叶级数的系数，n 是离散的，所以频谱也是离散状的每条谱线都出现在基波频率的整数倍上，其包络是取样函数。第二：谱线的间距是0w .。零点是0=2nw x τ，02w =T π是谱的基波频率。如果τ不变，T 增大，那么0w 减小，当T 非常大的时候，0w 非常小，谱线近似连续，越来越密，幅度越来越小。傅里叶变换：非周期函数正变换：--F jw)= ()iwt f t e dt ∝ ∝?（反变换：-1()=()2jnwt f t F jw e dw π ∝∝ ? 常用函数的傅里叶变换（典型非周期信号的频谱）

阿贝成像原理实验报告

佛山科学技术学院实验报告课程名称近代物理实验实验项目阿贝成像原理和空间滤波专业班级 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];

离散傅里叶变换及其快速算法

第五章离散傅里叶变换及其快速算法 1 离散傅里叶变换(DFT)的推导 (1) 时域抽样：目的：解决信号的离散化问题。效果：连续信号离散化使得信号的频谱被周期延拓。 (2) 时域截断：原因：工程上无法处理时间无限信号。方法：通过窗函数(一般用矩形窗)对信号进行逐段截取。结果：时域乘以矩形脉冲信号，频域相当于和抽样函数卷积。 (3) 时域周期延拓：目的：要使频率离散，就要使时域变成周期信号。 ! 方法：周期延拓中的搬移通过与)(s nT t -δ的卷积来实现。表示：延拓后的波形在数学上可表示为原始波形与冲激串序列的卷积。结果：周期延拓后的周期函数具有离散谱。 (4) 1。图1 DFT 推导过程示意图 (5) 处理后信号的连续时间傅里叶变换：∑∑ ∞ -∞=-=π--δ???? ? ????= k N n N kn j s kf f e nT h f H )()()(~ 010/2

(i) )(~ f H 是离散函数，仅在离散频率点S NT k T k kf f == =00处存在冲激，强度为k a ，其余各点为0。 (ii) )(~ f H 是周期函数，周期为s s T NT N T N Nf 100===，每个周期内有N 个不同的幅值。 (iii) 时域的离散时间间隔(或周期)与频域的周期(或离散间隔)互为倒数。 2 DFT 及IDFT 的定义 (1) ， (2) DFT 定义：设()s nT h 是连续函数)(t h 的N 个抽样值1,,1,0-=N n ，这N 个点的宽度为N 的DFT 为：[])1,...,1,0(,)()(1 /2-=???? ??==? -=π-∑N k NT k H e nT h nT h DFT s N n N nk j s s N (3) IDFT 定义：设??? ? ??s NT k H 是连续频率函数)(f H 的N 个抽样值1,,1,0-=N k ，这N 个点的宽度为N 的IDFT 为： ())1,...,1,0(,11 0/21 -==??? ? ? ?=??? ???? ?? ??? ???-=π--∑ N k nT h e NT k H N NT k H DFT s N k N nk j s s N (4) N nk j e /2π-称为N 点DFT 的变换核函数，N nk j e /2π称为N 点IDFT 的变换核函数。它们互为共轭。 (5) 同样的信号，宽度不同的DFT 会有不同的结果。DFT 正逆变换的对应关系是唯一的，或者说它们是互逆的。 (6) 引入N j N e W /2π-= (i) 用途： (a) 正逆变换的核函数分别可以表示为nk N W 和nk N W -。 (b) 核函数的正交性可以表示为：() )(* 1 0r n N W W kr N N k kn N -δ=∑-= (c) DFT 可以表示为：)1,,1,0(,)(10 -==? ??? ??∑ -=N k W nT h NT k H N n nk N s s (d) IDFT 可以表示为：)1,,1,0(,1 )(10 -=??? ? ??=∑-=-N n W NT k H N nT h N k nk N s s (ii) ） (iii) 性质：周期性和对称性： (a) 12==π-j N N e W (b) 12 /-==π-j N N e W (c) r N r N N N r N N W W W W ==+ (d) r N r N N N r N N W W W W -=-=+2/2/ (e) )(1Z m W m N ∈?= (f) ),(/2/2Z n m W e e W n N N n j mN mn j mn mN ∈?===π-π- 3 离散谱的性质 (1) 离散谱定义：称)(Z k NT k H H S k ∈??? ? ??=? 为离散序列)0)((N n nTs h <≤的DFT 离散谱，简称离散谱。 (2) 性质：

傅里叶光学实验

傅里叶光学的空间频谱与空间滤波实验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．信号、信息与消息的差别？信号：随时间变化的物理量；消息：待传送的一种以收发双方事先约定的方式组成的符号，如语言、文字、图像、数据等信息：所接收到的未知内容的消息，即传输的信号是带有信息的。 2．什么是奇异信号？函数本身有不连续点或其导数或积分有不连续点的这类函数统称为奇异信号或奇异函数。例如：单边指数信号（在t =0点时，不连续），单边正弦信号（在t =0时的一阶导函数不连续）。较为重要的两种奇异信号是单位冲激信号δ(t )和单位阶跃信号u(t )。 3．单位冲激信号的物理意义及其取样性质？冲激信号：它是一种奇异函数，可以由一些常规函数的广义极限而得到。它表达的是一类幅度很强，但作用时间很短的物理现象。其重要特性是筛选性，即： ()()()(0)(0)t x t dt t x dt x δδ∞ ∞ -∞ -∞ ==? ? 4．什么是单位阶跃信号？单位阶跃信号也是一类奇异信号，定义为： 10()00t u t t >?=?

12()()()x t ax t bx t =+，其中a 和b 是任意常数时，输出信号()y t 是1()y t 和2()y t 的线性叠加，即：12()()()y t ay t by t =+；且当输入信号()x t 出现延时，即输入信号是0()x t t -时，输出信号也产生同样的延时，即输出信号是0()y t t -。其中，如果当12()()()x t x t x t =+时，12()()()y t y t y t =+，则称系统具有叠加性；如果当1()()x t ax t =时，1()()y t ay t =则称系统具有均匀性。线性时不变系统是最基本的一类系统，是研究复杂系统，如非线性、时变系统的基础。 6．线性时不变系统的意义与应用？线性时不变系统是我们本课程分析和研究的主要对象，对线性时不变性进行推广，可以得到线性时不变系统具有微分与积分性质，假设系统的输入与输出信号分别为()x t 和()y t ，则当输入信号为 ()dx t dt 时，输出信号则为() dy t dt ；或者当输入信号为()t x d ττ-∞ ?时，输出信号则为()t y d ττ-∞ ?。另外，线性时不变系统对信号的处理作用可以用冲激响应（或单位脉冲响应）、系统函数或频率响应进行描述。而且多个系统可以以不同的方式进行连接，基本的连接方式为：级联和并联。假设两个线性时不变系统的冲激响应分别为：1()h t 和2()h t ，当两个系统级联后，整个系统的冲激响应为：12()()*()h t h t h t =；当两个系统并联后，整个系统的冲激响应为：12()()()h t h t h t =+；当0t <时，若()0h t =，则此系统为因果系统；若|()|h t dt ∞ -∞<∞?，则此系统为稳定系统。第二章连续时间系统的时域分析 1．如何获得系统的数学模型？数学模型是实际系统分析的一种重要手段，广泛应用于各种类型系统的分析和控制之中。不同的系统，其数学模型可能具有不同的形式和特点。对于线性时不变系统，其数学模型

傅里叶光学实验报告

实验原理：（略）实验仪器：光具座、氦氖激光器、白色像屏、作为物的一维、二维光栅、白色像屏、傅立叶透镜、小透镜实验内容与数据分析 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 =

按频率抽取基2-快速傅里叶逆变换算法_MATLAB代码

function x=MyIFFT_FB(y) %MyIFFT_TB:My Inverse Fast Fourier Transform Time Based %按频率抽取基2-傅里叶逆变换算法 %input: % y -- 傅里叶正变换结果,1*N的向量 %output: % x -- 逆变换结果，1*N的向量 %参考文献： % https://www.360docs.net/doc/2c15068411.html,/view/fea1e985b9d528ea81c779ee.html N=length(y); x=conj(y); %求共轭 x=MyFFT_FB(x);%求FFT x=conj(x);%求共轭 x=x./N;%除以N end %% 内嵌函数====================================================== function y=MyFFT_FB(x,n) %MYFFT_TB:My Fast Fourier Transform Frequency Based %按频率抽取基2-fft算法 %input: % x -- 输入的一维样本 % n -- 变换长度，缺省时n=length(x) 当n小于x数据长度时，x数据被截断到第n个数据% 当n大于时，x数据在尾部补0直到x 含n个数据 %output: % y -- 1*n的向量，快速傅里叶变换结果 %variable define: % N -- 一维数据x的长度 % xtem -- 临时储存x数据用 % m,M -- 对N进行分解N=2^m*M,M为不能被2整除的整数 % two_m -- 2^m % adr -- 变址，1*N的向量 % l -- 当前蝶形运算的级数 % W -- 长为N/2的向量，记录W(0,N),W(1,N),...W(N/2-1,N) % d -- 蝶形运算两点间距离 % t -- 第l级蝶形运算含有的奇偶数组的个数 % mul -- 标量，乘数 % ind1,ind2 -- 标量，下标 % tem -- 标量，用于临时储存 %参考文献： % https://www.360docs.net/doc/2c15068411.html,/view/fea1e985b9d528ea81c779ee.html %% 输入参数个数检查

常用傅立叶变换表

时域信号弧频率表示的傅里叶变换注释 1 线性 2 时域平移 3 频域平移, 变换2的频域对应4 如果值较大，则会收缩到原点附近，而会扩散并变得扁平. 当 | a | 趋向无穷时，成为 Delta函数。 5 傅里叶变换的二元性性质。通过交换时域变量和频域变量得到. 6 傅里叶变换的微分性质 7 变换6的频域对应 8 表示和的卷积—这

9 矩形脉冲和归一化的sinc 函数 10 变换10的频域对应。矩形函数是理想的低通滤波器，sinc 函数是这类滤波器对反因果冲击的响应。 11 tri 是三角形函数 12 变换12的频域对应 13 高斯函数 exp( ? αt 2) 的傅里叶变换是他本身. 只有当 Re(α) > 0时，这是可积的。 14 15 16 a>0 17 变换本身就是一个公式

18 δ(ω) 代表狄拉克δ函数分布. 这个变换展示了狄拉克δ函数的重要性：该函数是常函数的傅立叶变换 19 变换23的频域对应 20 由变换3和24得到. 21 由变换1和25得到，应用了欧拉公式: cos(at) = (e iat + e?iat) / 2. 22 由变换1和25得到 23 这里, n是一个自然数. δ(n)(ω) 是狄拉克δ函数分布的n阶微分。这个变换是根据变换7和24得到的。将此变换与1结合使用，我们可以变换所有多项式。 24 此处sgn(ω)为符号函数；注意此变换与变换7和24是一致的. 25 变换29的推广. 26 变换29的频域对应. 27 此处u(t)是单位阶跃函数; 此变换根据变换1和31得到.

光学仪器实验报告

常用光电仪器原理及使用实验报告班级：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 。

第七章傅里叶变换.

第七章傅里叶变换 1．求下列函数的傅氏变换： (1)1,10, ()1, 01,0,; t f t t --<? 解: (1)[()]()j t F f t f t e dt ω+∞--∞ =? 1 101 10 1 1 22sin cos | 2(1cos ).j t j t j t j t e dt e dt e dt e dt j i tdt t j ωωωωωωω ωω -----=-+=-+=-= =- -????? (2) ()()j t F f t e dt ωω+∞--∞ =? 0(1)(1)0 11|.11t j t j t j t e e dt e dt e j j ωωωωω ---∞ -∞ --∞====--?? 6．求下列函数的傅氏变换 (1) 1,0,sgn 1,0;t t t -? (2) ()sin(5).3f t t π =+ 解: (1)已知 1 [()](),[1]2(),F u t F j πδωπδωω = +=由sgn 2()1t u t =-有 12[sgn ]2( ())2().F t j j πδωπδωωω =+-= (2) 由于 1()sin(5)sin 5cos5,322f t t t t π=+=+ 故 [()][(5)(5)](5)(5)].2j F f t πδωδωδωδω= +--++- 7．已知00()[()()]F ωπδωωδωω=++-为函数()f t 的傅氏变换，求().f t

快速傅里叶变换实验报告

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

快速傅里叶变换实验报告机械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('傅立叶逆变换后时域函数'); 运行结果

实验二 快速傅里叶变换算法实验