“快速傅里叶变换的实现”实验报告

“快速傅里叶变换的实现”实验报告
“快速傅里叶变换的实现”实验报告

实验八 利用快速傅里叶变换(FFT)实现快速卷积(精选、)

实验八 利用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)

阿贝成像原理实验报告

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

C语言实现FFT(快速傅里叶变换)

C语言实现FFT(快速傅里叶变换) 函数原型:空快速傅立叶变换(Struct Compx *xin,Intn) 函数函数:对输入复数组执行快速傅立叶变换(FFT)输入参数:*xin复结构组的第一个地址指针。结构输出参数:no * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *结构compx u,w,t。 nv2 =快速傅立叶变换_ N/2;nm1 =快速傅立叶变换_ N-1;(I = 0;i

傅里叶光学实验

傅里叶光学的空间频谱与空间滤波实验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、熟悉并掌握按时间抽取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)=

快速傅里叶变换FFT的FPGA设计与实现--电科1704 郭衡

快速傅里叶变换FFT的FPGA设计与实现 学生姓名郭衡 班级电科1704 学号17419002064 指导教师谭会生 成绩 2020年5 月20 日

快速傅里叶变换FFT 的设计与实现 一、研究项目概述 非周期性连续时间信号x(t)的傅里叶变换可以表示为:= )(?X dt t j e t x ? ∞ ∞ --1 )(?,式中计算出来的是信号x(t)的连续频谱。但是,在实际的控制系统中能够式中计算出来的是信号x(t)的连续频谱。但是,在实际的控制系统中能够算信号x(t)的频谱。 有限长离散信号x(n),n=0,1,…,N-1的DFT 定义为: ∑-=-=-==1 02,1.....10)()(N n N j N kn N e W N k W n x K X π、、。 可以看出,DFT 需要计算大约N2次乘法和N2次加法。当N 较大时,这个计算量是很大的。利用WN 的对称性和周期性,将N 点DFT 分解为两个N /2点的DFT ,这样两个N /2点DFT 总的计算量只是原来的一半,即(N /2)2+(N /2)2=N2/2,这样可以继续分解下去,将N /2再分解为N /4点DFT 等。对于N=2m 点的DFT 都可以分解为2点的DFT ,这样其计算量可以减少为(N /2)log2N 次乘法和Nlog2N 次加法。图1为FFT 与DFT-所需运算量与计算点数的关系曲线。由图可以明显看出FFT 算法的优越性。 图1 FFT 与DFT 所需乘法次数比 较

X[1] 将x(n)分解为偶数与奇数的两个序列之和,即x(n)=x1(n)+x2(n)。 x1(n)和x2(n)的长度都是N /2,x1(n)是偶数序列,x2(n)是奇数序列,则 ∑∑=--=-=+2 )12(120 2)1.....,0()(2)(1)(N n k n N N n km N N k W n x W n x K X 所以)1...,0()(2)(1)(12 22120 -=+=∑∑-=-=N k W n x W W n x K X N n km N k N km N N n 由于km N N j km N j km N W e e W 2/2 /2222===--ππ ,则 )1.....,0)((2)(1)(2)(1)(12 2/120 2/-=+=+=∑∑-=-=N k k X W k X W n x W W n x K X k N N n km N k N N n kn N 其中X1(k)和X2(k)分别为x1(n)和x2(n)的N /2点DFT 。由于X1(k)和X2(k)均以N /2为周期,且WNk+N/2=-WNk ,所以X(k)又可表示为: )12/....,1,0)((2)(1)(-=+=N k k X W k X K X k N )12/....,1,0)((2)(1)2/(-=-=+N k k X W k X N K X k N

傅里叶光学实验报告

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

快速傅里叶变换实验报告

快速傅里叶变换实验报告

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

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

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

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

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

信息工程学院实验报告 课程名称:信号与系统 实验项目名称:实验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,故满足采样定理,不会发生混叠现象。 频谱图: 在忽略旁瓣信号的情况下,可近似认为:

C语言实现FFT(快速傅里叶变换)

#include #include /********************************************************************* 快速福利叶变换C函数 函数简介:此函数是通用的快速傅里叶变换C语言函数,移植性强,以下部分不依赖硬件。此函数采用联合体的形式表示一个复数,输入为自然顺序的复 数(输入实数是可令复数虚部为0),输出为经过FFT变换的自然顺序的 复数 使用说明:使用此函数只需更改宏定义FFT_N的值即可实现点数的改变,FFT_N的应该为2的N次方,不满足此条件时应在后面补0 函数调用:FFT(s); 时间:2010-2-20 版本:Ver1.0 参考文献: **********************************************************************/ #include #define PI 3.1415926535897932384626433832795028841971 //定义圆周率值#define FFT_N 128 //定义福利叶变换的点数 struct compx {float real,imag;}; //定义一个复数结构struct compx s[FFT_N]; //FFT输入和输出:从S[1]开始存放,根据大小自己定义 /******************************************************************* 函数原型:struct compx EE(struct compx b1,struct compx b2) 函数功能:对两个复数进行乘法运算 输入参数:两个以联合体定义的复数a,b 输出参数:a和b的乘积,以联合体的形式输出 *******************************************************************/ struct compx EE(struct compx a,struct compx b) { struct compx c; c.real=a.real*b.real-a.imag*b.imag; c.imag=a.real*b.imag+a.imag*b.real; return(c); } /***************************************************************** 函数原型:void FFT(struct compx *xin,int N)

快速傅里叶变换实验报告

快速傅里叶变换实验报告 机械34班 攀 2013010558 一、 基本信号(函数)的FFT 变换 1. 000()sin()sin 2cos36 x t t t t πωωω=+++ 1) 采样频率08s f f =,截断长度N=16; 取02ωπ=rad/s ,则0f =1Hz ,s f =8Hz ,频率分辨率f ?=s f f N ?==0.5Hz 。 最高频率c f =30f =3Hz ,s f >2c f ,故满足采样定理,不会发生混叠现象。 截断长度02T T =,整周期截取,不会发生栅栏效应。理论上有一定的泄漏,但在整周期 截取的情况下,旁瓣上的采样都约为 0,泄漏现象没有体现出来。 频谱图如下:

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

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

相关文档
最新文档