实验2 离散傅里叶变换(DFT)
一、实验目的
(1)加深对离散傅里叶变换(DFT)基本概念的理解。
(2)了解有限长序列傅里叶变换(DFT)与周期序列傅里叶级数(DFS)、离散时间傅里叶变换(DTFT)的联系。
(3)掌握用MA TLAB 语言进行离散傅里叶变换和逆变换的方法。
二、实验内容
1.有限长序列的傅里叶变换(DFT)和逆变换(IDFT)
2.有限长序列DFT 与周期序列DFS 的联系
3.有限长序列DFT 与离散时间傅里叶变换DTFT 的联系
三、实验环境
MA TLAB7.0
四、实验原理
1.有限长序列的傅里叶变换(DFT)和逆变换(IDFT)
在实际中常常使用有限长序列。如果有限长序列信号为x(n)，则该序列的离散傅里叶变换对可以表示为
1N ,0,1,k ,
W x (n)DFT[x (n)]X(k)1
N 0n nk
N -===∑-= (2-1)
1
N ,0,1,n ,
W X(k)N 1IDFT[X(k)]x (n)1N 0
k nk N -===∑-=- (2-2)
从离散傅里叶变换定义式可以看出，有限长序列在时域上是离散的，在频域上也是离散
的。式中，N
π2j N e W -=即仅在单位圆上N 个等间距的点上取值，这为使用计算机进行
处理带来了方便。
由有限长序列的傅里叶变换和逆变换定义可知，DFT 和DFS 的公式非常相似，因此在程序编写上也基本一致。
例2-1 已知x(n)＝［0，1，2，3，4，5，6，7］，求x(n)的DFT 和IDFT 。要求： (1)画出序列傅里叶变换对应的|X(k)|和arg ［X(k)］图形。 (2)画出原信号与傅里叶逆变换IDFT ［X(k)］图形进行比较。 解 MA TLAB 程序如下：
>> title('|X(k)|');
运行结果如图2-1所示。
024680
2468x(n)
02468
2468IDFT|X (k)|
2
4
6
8
10
20
30|X (k)|
2
4
6
8
-4-2
24arg|X (k)|
图2-1 例2-1有限长序列的傅里叶变换和逆变换结果
从得到的结果可见，与周期序列不同的是，有限长序列本身是仅有N 点的离散序列，相当于周期序列的主值部分。因此，其频谱也对应序列的主值部分，是含N 点的离散序列。
2.有限长序列DFT 与周期序列DFS 的联系
将周期序列的傅里叶级数变换与有限长序列离散傅里叶变换对(式(2-1)和式(2-2))进行比较，可以看出两者的区别仅仅是将周期序列换成了有限长序列。
例2-2 已知周期序列的主值x(n)＝［0，1，2，3，4，5，6，7］，求x(n)周期重复次数为4次时的DFS 。要求：
(1)画出原主值和信号周期序列信号。
(2)画出序列傅里叶变换对应的(k)X ~
和(k)]X ~
arg[的图形。 解 MA TLAB 程序如下：
运行结果如图2-2所示。
024680
2468原主值信号x(n)
010203040
2468周期序列信号
10
20
30
40
50
100
150|X (k)|
10
20
30
40
-4-2
24arg|X (k)|
图2-2 例2-2周期序列的傅里叶级数(DFS)结果
由这个周期序列的实验我们可以看出，与例2-1相比，有限长序列x(n)可以看成是周期序列(n)x ~的一个周期；反之，周期序列(n)x ~
可以看成是有限长序列x(n)以N 为周期的周期延拓。频域上的情况也是相同的。从这个意义上说，周期序列只有有限个序列值有意义。
3.有限长序列DFT 与离散时间傅里叶变换DTFT 的联系
离散时间傅里叶变换(DTFT)是指信号在时域上为离散的，而在频域上则是连续的。 如果离散时间非周期信号为x(n)，则它的离散傅里叶变换对(DTFT)表示为
∑
∞
-∞
=-=
=n j ωω
j ωe x (n))X(e DTFT[x (n)]d ωe )X(e π
21
x(n))]IDTFT[X(e j ωωπ
π
j ωj ω
?
-=
=
其中X(e jw )称为信号序列的频谱。将频谱表示为)
ω(j ωj ωe )X(e )X(e ?=.|X(e jw )|
称为序列的幅度谱,
[])X(e arg )ω(j ω=?称为序列的相位谱。
从离散时间傅里叶变换的定义可以看出，信号在时域上是离散的、非周期的，而在频域上则是连续的、周期性的。
与有限长序列相比，X(e jw )仅在单位圆上取值，X(k)是在单位圆上N 个等间距的点上取值。因此，连续谱X(e jw )可以由离散谱X(k)经插值后得到。
为了进一步理解有限长序列的傅里叶变换(DFT)与离散时间傅里叶变换(DTFT)的联系，我们举例说明离散时间傅里叶变换的使用方法和结果。
例2-3 求x(n)＝［0，1，2，3，4，5，6，7］，0≤n ≤7的DTFT ，将(－2p ，2p)区间分成500份。要求：
(1)画出原信号。
(2)画出由离散时间傅里叶变换求得的幅度谱X(ejw)和相位谱arg ［X(ejw)］图形。 解 MA TLAB 程序如下：
运行结果如图2-3所示。
01234567
0510
x (n )
-6-4
-2
2
4
6
10
2030
幅度谱
-6
-4
-2
2
4
6
-202
相位谱
图2-3 例2-3离散时间傅里叶变换(DTFT)的结果
由图2-3与DFT 的结果图2-1相比可以看出，两者有一定的差别。主要原因在于，该例进行DTFT 时，X(e jw )在单位圆上取250个点进行分割；而图2-1进行DFT 时，X(k)是在单位圆上N ＝8的等间距点上取值，X(k)的序列长度与X(e jw )相比不够长。
例2-4 仍然用x(n)＝［0，1，2，3，4，5，6，7］，将x(n)的有限长序列后面补足至N ＝100，求其DFT ，并与例2-3进行比较。
解 将例2-1程序的前2行改为 N ＝100；
xn ＝［0，1，2，3，4，5，6，7，zeros (1，N －8)］；
则|X(k)|和arg ［X(k)］的图形接近由离散时间傅里叶变换求得的幅度谱X(e jw )和相位谱arg ［X(ejw)］的图形，如图2-4所示。注意，此图对应［0，2p ］区间。
MATLAB 程序如下：
运
行结果如图2-4所示。
0102030405060708090100
10
20
30|X (k)|
10
20
30
40
50
60
70
80
90
100
-4-2024arg|X (k)|
图2-4 增长有限长序列的长度得到|X(k)|和arg ［X(k)］
五、实验内容
1. 已知有限长序列x(n)＝［7，6，5，4，3，2］，求x(n)的DFT 和IDFT 。要求： ①画出序列傅里叶变换对应的|X(k)|和arg ［X(k)］的图形。 ②画出原信号与傅里叶逆变换IDFT ［X(k)］的图形进行比较。
解：
运行结果：
2. 已知周期序列的主值x(n)＝［7，6，5，4，3，2］，求x(n)周期重复次数为3次时的DFS 和IDFS 。要求：
①画出原信号序列的主值和周期序列的图形。
②画出序列傅里叶变换对应的(k)X ~
和(k)]X ~arg[的图形。 解：
3. 求x(n)＝［7，6，5，4，3，2］，0≤n ≤5的DTFT ，将(－2p ，2p)区间分成500份。要求：
①求出原信号。
②画出由离散时间傅里叶变换求得的幅度谱X(ejw)和相位谱arg ［X(ejw)］的图形。
③求有限长序列x(n)＝［7，6，5，4，3，2］，N ＝100时的DFT ，并与DTFT 的结果进行比较。 解 ：
运行结果：
4. 讨论序列补零及增加数据长度对信号频谱的影响
(1) 求出序列x(n)=cos(0.48 n)+cos(0.52 n)基于有限个样点n=10的频谱 (2) 求n=100时，取x(n)的前10个，后90个设为零，得到x(n)的频谱 (3) 增加x(n)有效的样点数，取100个样点得到x(n)的频谱
六 实验要求
1. 将1-4题中的所有图的title 都以自己名字的全拼开头来命名； 2．根据自己的理解，借鉴例题的主要程序编程完成仿真； 3. 对仿真结果作相应的分析。
快速傅里叶变换实验报告
————————————————————————————————作者：————————————————————————————————日期: ?
快速傅里叶变换实验报告 机械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，泄漏现象没有体现出来。 频谱图：
佛山科学技术学院 实验报告 课程名称近代物理实验实验项目阿贝成像原理和空间滤波 专业班级 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解释阿贝成像实验
班级信工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];
实验八 利用FFT 实现快速卷积 一、 实验目的 （1） 通过这一实验，加深理解FFT 在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT 进行数字信号处理。 （2） 进一步掌握循环卷积和线性卷积两者之间的关系。 二、 实验原理与方法 数字滤波器根据系统的单位脉冲响应h(n)是有限长还是无限长可分为有限长单位脉冲响应（Finite Impulse Response ）系统（简记为FIR 系统）和无限长单位脉冲响应（Infinite Impulse Response ）系统（简记为IIR 系统）。 对于FIR 滤波器来说，除了可以通过数字网络来实现外，也可以通过FFT 的变换来实现。 一个信号序列x(n)通过FIR 滤波器时，其输出应该是x(n)与h(n)的卷积： ∑+∞ -∞ =-= =m m n h m x n h n x n y )()()(*)()( 或 ∑+∞ -∞ =-= =m m n x m h n x n h n y ) ()()(*)()( 当h(n)是一个有限长序列，即h(n)是FIR 滤波器，且10-≤≤N n 时 ∑-=-=1 0) ()()(N m m n x m h n y 在数字网络（见图6.1）类的FIR 滤波器中，普遍使用的横截型结构（见下图6.2 图6.1 滤波器的数字网络实现方法 图6.2 FIR 滤波器横截型结构 y(n) y(n) -1-1-1-1
应用FFT 实现数字滤波器实际上就是用FFT 来快速计算有限长度列间的线性卷积。 粗略地说，这种方法就是先将输入信号x(n)通过FFT 变换为它的频谱采样 值X(k)，然后再和FIR 滤波器的频响采样值H(k)相乘，H(k)可事先存放在存储器中，最后再将乘积H(k)X(k)通过快速傅里叶变换（简称IFFT ）还原为时域序列，即得到输出y(n)如图6.3所示。 图6.3 数字滤波器的快速傅里叶变换实现方法 现以FFT 求有限长序列间的卷积及求有限长度列与较长序列间的卷积为例来讨论FFT 的快速卷积方法。 （1） 序列)(n x 和)(n h 的列长差不多。设)(n x 的列长为1N ，)(n h 的列长为2N ，要求 )()(n x n y =N ∑-=-==1 ) ()()(*)()(N r r n h r x n h n x n h 用FFT 完成这一卷积的具体步骤如下： i. 为使两有限长序列的线性卷积可用其循环卷积代替而不发生混叠，必须选择循环卷积长度121-+≥N N N ，若采用基2-FFT 完成卷积运 算，要求m N 2=（m 为整数）。 ii. 用补零方法使)(n x ，)(n h 变成列长为N 的序列。 ?? ?-≤≤-≤≤=10 10)()(11N n N N n n x n x ?? ?-≤≤-≤≤=10 1 0)()(22N n N N n n h n h iii. 用FFT 计算)(),(n h n x 的N 点离散傅里叶变换 )()(k X n x FFT ??→? )()(k H n h FFT ??→? iv. 做)(k X 和)(k H 乘积，)()()(k H k X k Y ?= v. 用FFT 计算)(k Y 的离散傅里叶反变换得 y(n)
1 / 7 信息工程学院实验报告 课程名称：信号与系统 实验项目名称：实验 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()[()] ()2 j 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 Fv fte 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 。 成 绩： 指导教师（签名）：
傅里叶光学的空间频谱与空间滤波实验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 各种形式的空间滤波器
陕西科技大学实验报告 班级信工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、熟悉并掌握按时间抽取FFT 算法的程序； 3、了解应用FFT 进行信号频谱分析过程中可能出现的问题，例如混淆、泄漏、 栅栏效应等，以便在实际中正确应用FFT 。 二 实验内容 a ） 信号频率F ＝50Hz ，采样点数N=32,采样间隔T= matlab 程序代码为： F=50; T=; N=32; n=0:N-1; t=n*T; A=sin(2*pi*F*t); figure; Y = fft(A,N); h = (abs(Y)); h=h/max(h(1:N)); for n=1:N; string1=strcat('X(',num2str(n-1), ')=',num2str(h(n))); disp(string1); f=(n/T)/N; end stem([0:N-1]/N/T,h); xlabel('?μ?ê/HZ'); ylabel('??·ùX￡¨ejw￡?'); title('·ù?μì?D?'); 上述代码命令中，将FFT 变换后的数字变量K ，在画图时转换成频域中的频率f 。这主 要是根据数字频率与模拟域频率之间的关系： T Ω=ω 其中ω、Ω分别为数字和模拟域中的频率，且N k πω2= f π2=Ω 于是有： NT k f = 运算结果： X(1)=1 X(2)= X(3)= X(4)=
X(5)= X(6)= X(7)= X(8)= X(9)= X(10)= X(11)= X(12)= X(13)= X(14)= X(15)= X(16)= X(17)= X(18)= X(19)= X(20)= X(21)= X(22)= X(23)= X(24)= X(25)= X(26)= X(27)= X(28)= X(29)= X(30)= X(31)=1 b）信号频率F＝50Hz，采样点数N=32,采样间隔T= 同理可将a）中F、N、T，参数改成要求值（以下均是如此），即可得，X(0)= X(1)= X(2)= X(3)= X(4)= X(5)= X(6)= X(7)= X(8)=1 X(9)= X(10)= X(11)= X(12)= X(13)= X(14)= X(15)= X(16)= X(17)= X(18)= X(19)= X(20)= X(21)= X(22)= X(23)= X(24)=1 X(25)= X(26)= X(27)= X(28)= X(29)= X(30)= X(31)=
实验四：离散傅里叶变换 实验原理： DFT的快速算法FFT利用了的三个固有特性：（1）对称性（2）周期性（3）可约性。FFT算法基本上可以分为两大类，即按时间抽选法(DIT,Decimation-In-Time)和按频率抽选法(DIF,Decimation-In-frequency)。 MATLAB中提供了进行快速傅里叶变换的fft函数： X=fft(x)，基2时间抽取FFT算法，x是表示离散信号的向量；X是系数向量； X=fft(x,N),补零或截断的N点DFT，当x得长度小于N时，对补零使其长度为N，当x的长度大于N时，对x截断使其长度为N。 实验内容： =60; n=[0:1:k/2]; xa1=2*sin(10*pi*n/k)+cos(18*pi*n/k); subplot(321) stem(n,xa1) xlabel('N');ylabel('x(n)'); xk1=fft(xa1);xk1=abs(xk1) subplot(322) stem(n,xk1) xlabel('k');ylabel('X(k)'); n=[0:1:k*]; xa1=2*sin(10*pi*n/k)+cos(18*pi*n/k); subplot(323) stem(n,xa1) xlabel('N');ylabel('x(n)'); xk1=fft(xa1);xk1=abs(xk1) subplot(324) stem(n,xk1) xlabel('k');ylabel('X(k)'); n=[0:1:k*2]; xa1=2*sin(10*pi*n/k)+cos(18*pi*n/k); subplot(325) stem(n,xa1) xlabel('N');ylabel('x(n)'); xk1=fft(xa1);xk1=abs(xk1) subplot(326) stem(n,xk1) xlabel('k');ylabel('X(k)');
实验原理：（略） 实验仪器： 光具座、氦氖激光器、白色像屏、作为物的一维、二维光栅、白色像屏、傅立叶透镜、小透镜 实验内容与数据分析 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 =
离散信号的变换（MATLAB 实验） 一、实验目的 掌握用Z 变换判断离散系统的稳定与否的方法，掌握离散傅立叶变换及其基本性质和特点，了解快速傅立叶变换。 二、实验内容 1、已经系统函数为 5147.13418.217.098.2250 5)(2342-++--+=z z z z z z Z H (1) 画出零极点分布图，判断系统是否稳定； (2)检查系统是否稳定； (3) 如果系统稳定，求出系统对于u(n)的稳态输出和稳定时间b=[0,0,1,5,-50];a=[2,-2.98,0.17,2.3418,-1.5147]; subplot(2,1,1);zplane(b,a);title('零极点分布图'); z=roots(a); magz=abs(z) magz = 0.9000 0.9220 0.9220 0.9900 n=[0:1000]; x=stepseq(0,0,1000); s=filter(b,a,x); subplot(2,1,2);stem(n,s);title('稳态输出'); （1）因为极点都在单位园内，所以系统是稳定的。 （2）因为根的幅值（magz ）都小于1，所以这个系统是稳定的。 （3）稳定时间为570。 2、综合运用上述命令，完成下列任务。 (1) 已知)(n x 是一个6点序列： ???≤≤=其它,050,1)(n n x
计算该序列的离散时间傅立叶变换，并绘出它们的幅度和相位。 要求：离散时间傅立叶变换在[-2π,2π]之间的两个周期内取401个等分频率上进行数值求值。 n=0:5;x=ones(1,6); k=-200:200;w=(pi/100)*k; X=x*(exp(-j*pi/100)).^(n'*k); magX=abs(X);angX=angle(X); subplot(2,1,1);plot(w/pi,magX);grid;title('幅度'); subplot(2,1,2);plot(w/pi,angX);grid;title('相位'); (2) 已知下列序列： a. ,1000),52.0cos()48.0cos()(≤≤+=n n n n x ππ； b ．)4sin()(πn n x =是一个N ＝32的有限序列； 试绘制)(n x 及它的离散傅立叶变换 )(k X 的图像。 a ． n=[0:1:100];x=cos(0.48*pi*n)+cos(0.52*pi*n); subplot(2,1,1);plot(n,x);title('x(n)的图像'); X=dft(x,101); magX=abs(X); subplot(2,1,2);plot(n,magX);title('丨X(k)丨的图像');
常用光电仪器原理及使用 实验报告 班级：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 傅里叶变换及其性质 实验时间：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 。
戶幵,戈丿、弟实验报告 课程名称：彳 _____________ 指导老师 _____________ 成绩： ____________________ 实验名称：离散傅里叶变换和快速傅里叶变换 实验类型: _________________ 同组学生姓名： 一、实验目的和要求(必填) 二、实验内容和原理(必填) 三、主要仪器设备(必填) 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析(必填) 七、讨论、心得 一、实验目的和要求 1. 掌握DFT 的原理和实现 2. 掌握FFT 的原理和实现，掌握用 FFT 对连续信号和离散信号进行谱分析的方法。 二、实验内容和原理 2.1 DTFT 和 DFT N 1 如果x(n)为因果有限长序列,n=0,1,...,N-1，则x(n)的DTFT 表示为：X(e j ) x(n)e n 0 序列的N 点DFT 是DTFT 在［0,2 n 上的N 点等间隔采样，采样间隔为2 d N 。通过DFT , 可以完成由一组有限个信号采样值 x(n)直接计算得到一组有限个频谱采样值 X(k)。X(k)的幅 度谱为X(k) v 'x R (k ) X |2(k ) , X R (k)和X i (k)分别为X(k)的实部和虚部。X(k)的相位谱 为(k) 列吩 序列x(n)的离散事件傅里叶变换(DTFT )表示为: X(e j ) x( n)e x(n)的离散傅里叶变换(DFT )表达式为: X(k) x(n)e n 0 j^nk N (k 0,1,…,N 1)
IDFT )定义为 x(n)丄 N \(k)e j_Nnk (n 0,1，…，N 1) N n 0 2.2 FFT 快速傅里叶变换(FFT )是DFT 的快速算法，它减少了 DFT 的运算量，使数字信号的处理 速度大大提高。 三、主要仪器设备 PC 一台，matlab 软件 四、实验内容 4.1第一题 x(n)的离散时间 傅里叶变换(DTFT ) X(e j Q )并绘图。 0 其2他n 2; (2)已知 x(n) 2n 0 n 10。 0其他 4.1.1理论分析 1) 由DTFT 计算式, X (Q)是实数，可以直接作出它的图像。 离散傅里叶反变换 求有限长离散时间信号 (1)已知 x(n) X( ) x(n)e j n e 2j 1 5j e 1 e j e 2? e 2? 0.5j e 0.5 j e sin(2.5 )
信息工程学院实验报告 课程名称：信号与系统 实验项目名称：实验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班刘攀 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，故满足采样定理，不会发生混叠现象。 频谱图： 在忽略旁瓣信号的情况下，可近似认为：
实验四离散傅里叶变换的性质 一、实验目的 1. 熟悉matlab软件中离散傅里叶变换的实现方法及FFT函数的使用方法； 2. 通过软件仿真，加深对离散傅里叶变换性质的理解。 二、实验内容 1. 验证离散傅里叶变换的线性性质； 2. 掌握用matlab实现圆周移位的方法； 3. 验证圆周卷积与线性卷积的关系。 三、实验步骤 1. 验证线性性质 设两个有限长序列分别为xn1=[3,1,-2,2,3,4]，xn2=[1,1,1,1]，做4DFT[xn1]+2DFT[xn2],及DFT[4xn1+2xn2]的运算，比较它们的结果。 代码如下： clear,N=20;n=[0:1:N-1]; xn1=[3,1,-2,2,3,4];n1=0:length(xn1)-1; %定义序列xn1 xn2=[1,1,1,1];n2=0:length(xn2)-1; %定义序列xn2 yn1=4*xn1;yn2=2*xn2;[yn,ny]=seqadd(yn1,n1,yn2,n2); %计算4xn1+2xn2 xk1=fft(xn1,N);xk2=fft(xn2,N); %分别求DFT[xn1] 和DFT[xn2] yk0=4*xk1+2*xk2; %计算4DFT[xn1]+2DFT[xn2] yk=fft(yn,N); %计算DFT[4xn1+2xn2] subplot(2,1,1);stem(n,yk0);title('傅里叶变换之和') %显示4DFT[xn1]+2DFT[xn2] subplot(2,1,2);stem(n,yk);title('序列和之傅里叶变换') %显示DFT[4xn1+2xn2] 运行结果如图1所示，从图中可知，用两种方法计算的DFT完全相等，所以离散傅里叶变换的线性性质得到验证。