【傅里叶变换】

1、为什么要进行傅里叶变换,其物理意义是什么?
傅立叶变换是数字信号处理领域一种很重要的算法。要知道傅立叶变换算法的意义,首先要了解傅立叶原理的意义。傅立叶原理表明:任何连续测量的时序或信号,都可以表示为不同频率的正弦波信号的无限叠加。而根据该原理创立的傅立叶变换算法利用直接测量到的原始信号,以累加方式来计算该信号中不同正弦波信号的频率、振幅和相位。

和傅立叶变换算法对应的是反傅立叶变换算法。该反变换从本质上说也是一种累加处理,这样就可以将单独改变的正弦波信号转换成一个信号。

因此,可以说,傅立叶变换将原来难以处理的时域信号转换成了易于分析的频域信号(信号的频谱),可以利用一些工具对这些频域信号进行处理、加工。最后还可以利用傅立叶反变换将这些频域信号转换成时域信号。

从现代数学的眼光来看,傅里叶变换是一种特殊的积分变换。它能将满足一定条件的某个函数表示成正弦基函数的线性组合或者积分。在不同的研究领域,傅里叶变换具有多种不同的变体形式,如连续傅里叶变换和离散傅里叶变换。

在数学领域,尽管最初傅立叶分析是作为热过程的解析分析的工具,但是其思想方法仍然具有典型的还原论和分析主义的特征。"任意"的函数通过一定的分解,都能够表示为正弦函数的线性组合的形式,而正弦函数在物理上是被充分研究而相对简单的函数类:1. 傅立叶变换是线性算子,若赋予适当的范数,它还是酉算子;2. 傅立叶变换的逆变换容易求出,而且形式与正变换非常类似;3. 正弦基函数是微分运算的本征函数,从而使得线性微分方程的求解可以转化为常系数的代数方程的求解.在线性时不变杂的卷积运算为简单的乘积运算,从而提供了计算卷积的一种简单手段;5.
离散形式的傅立叶的物理系统内,频率是个不变的性质,从而系统对于复杂激励的响应可以通过组合其对不同频率正弦信号的响应来获取;4. 著名的卷积定理指出:傅立叶变换可以化复变换可以利用数字计算机快速的算出(其算法称为快速傅立叶变换算法(FFT))。

正是由于上述的良好性质,傅里叶变换在物理学、数论、组合数学、信号处理、概率、统计、密码学、声学、光学等领域都有着广泛的应用。

2、图像傅立叶变换的物理意义

图像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度。如:大面积的沙漠在图像中是一片灰度变化缓慢的区域,对应的频率值很低;而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域,对应的频率值较高。傅立叶变换在实际中有非常明显的物理意义,设f是一个能量有限的模拟信号,则其傅立叶变换就表示f的谱。从纯粹的数学意义上看,傅立叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看,傅立叶变换是将图像从空间域转换到频率域,其逆变换是将图像从频率域转换到空间域。换句话说,傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数,傅立叶逆变换是将图像的频率分布函数变换为灰度分布函数。

傅立叶变换以前,图像(未压缩的位图)是由对在连续空间(现实空间)上的采样得到一系列点的集合,我们习惯用一个二维矩阵表示空间上各点,则图像可由z=f(x,y)来表示。由于空间是三维的,图像是二维的,因此空间中物体在另一个维度上的关系就由梯度来表示,这样我们可以通过观察图像得知物体在三维空间中的对应关系。为什么要提梯度?因为实际上对图像进行二维傅立叶变换得到频谱图,就是图像梯度的分布图,当然频谱图上的各点与图像上各点并不存在一一对应的关系,即使在不移频的情况下也是没有。傅立叶频谱图上我们看到的明暗不一的亮点,实际上图像上某一点与邻域点差异的强弱,即梯度的大小,也即该点的频率的大小(可以这么理解,图像中的低频部分指低梯度的点,高频部分相反)。一般来讲,梯度大则该点的亮度强,否则该点亮度弱。这样通过观察傅立叶变换后的频谱图,也叫功率图,我们首先就可以看出,图像的能量分布,如果频谱图中暗的点数更多,那么实际图像是比较柔和的(因为各点与邻域差异都不大,梯度相对较小),反之,如果频谱图中亮的点数多,那么实际图像一定是尖锐的,边界分明且边界两边像素差异较大的。对频谱移频到原点以后,可以看出图像的频率分布是以原点为圆心,对称分布的。将频谱移频到圆心除了可以清晰地看出图像频率分布以外,还有一个好处,它可以分离出有周期性规律的干扰信号,比如正弦干扰,一副带有正弦干扰,移频到原点的频谱图上可以看出除了中心以外还存在以某一点为中心,对称分布的亮点集合,这个集合就是干扰噪音产生的,这时可以很直观的通过在该位置放置带阻滤波器消除干扰。

另外我还想说明以下几点:
1、图像经过二维傅立叶变换后,其变换系数矩阵表明:
若变换矩阵Fn原点设在中心,其频谱能量集中分布在变换系数短阵的中心附近(图中阴影区)。若所用的二维傅立叶变换矩阵Fn的原点设在左上角,那么图像信号能量将集中在系数矩阵的四个角上。这是由二维傅立叶变换本身性质决定的。同时也表明一股图像能量集中低频区域。
2 、变换之后的图像在原点平移之前四角是低频,最亮,平移之后中间部分是低频,最亮,亮度大说明低频的能量大(幅角比较大)

另一个人的说法
傅里叶变换是一种解决问题的方法,一种工具,一种看待问题的角度。
  理解的关键是:一个连续的信号可以看作是一个个小信号的叠加,从时域叠加与从频域叠加都可以组成原来的信号,将信号这么分解后有助于处理。
  我们原来对一个信号其实是从时间的角度去理解的,不知不觉中,其实是按照时间把信号进行分割,每一部分只是一个时间点对应一个信号值,一个信号是一组这样的分量的叠加。傅里叶变换后,其实还是个叠加问题,只不过是从频率的角度去叠加,只不过每个小信号是一个时间域上覆盖整个区间的信号,但他确有固定的周期,或者说,给了一个周期,我们就能画出一个整个区间上的分信号,那么给定一组周期值(或频率值),我们就可以画出其对应的曲线,就像给出时域上每一点的信号值一样,不过如果信号是周期的话 ,频域的更简单,只需要几个甚至一个就可以了,时域则需要整个时间轴上每一点都映射出一个函数值。
  
  傅里叶变换就是将一个信号的时域表示形式映射到一个频域表示形式;逆傅里叶变换恰好相反。这都是一个信号的不同表示形式。它的公式会用就可以,当然把证明看懂了更好。
  
  对一个信号做傅立叶变换,可以得到其频域特性,包括幅度和相位两个方面。幅度是表示这个频率分量的大小,那么相位呢,它有什么物理意义?频域的相位与时域的相位有关系吗?信号前一段的相位(频域)与后一段的相位的变化是否与信号的频率成正比关系?
  傅立叶变换就是把一个信号,分解成无数的正弦波(或者余弦波)信号。也就是说,用无数的正弦波,可以合成任何你所需要的信号。
  想一想这个问题,给你很多正弦信号,你怎样才能合成你需要的信号呢。
  答案是要两个条件,一个是每个正弦波的幅度,另一个就是每个正弦波之间的相位差。
  所以现在应该明白了吧,频域上的相位,就是每个正弦波之间的相位。
  
  傅立叶变换用于信号的频率域分析,一般我们把电信号描述成时间域的数学模型,而数字信号处理对信号的频率特性更感兴趣,而通过傅立叶变换很容易得到信号的频率域特性。
  
  傅里叶变换简单通俗理解就是把看似杂乱无章的信号考虑成由一定振幅、相位、频率的基本正弦(余弦)信号组合而成,傅里叶变换的目的就是找出这些基本正弦(余弦)信号中振幅较大(能量较高)信号对应的频率,从而找出杂乱无章的信号中的主要振动频率特点。
  如减速机故障时,通过傅里叶变换做频谱分析,根据各级齿轮转速、齿数与杂音频谱中振幅大的对比,可以快速判断哪级齿轮损伤。

时间: 2024-08-17 17:05:13

【傅里叶变换】的相关文章

二维傅里叶变换转一维傅里叶变换

问题描述 二维傅里叶变换转一维傅里叶变换 看了一个论文说二维傅里叶变换图谱经圆积分计算后可以得到一维傅里叶变换,求大神推导公式或者过程 解决方案 // ==============================================================================// 快速离散傅里叶变换和功率谱// 一维快速傅里叶变换FFT1和二维快速傅里叶变换FFT2// 测试环境 C++ builder 2010// MinGW (GCC) 4.5// wuxup

傅立叶变换-是否有支持定点化,支持NEON, 支持任意长度的的傅里叶变换的库

问题描述 是否有支持定点化,支持NEON, 支持任意长度的的傅里叶变换的库 现在正在做实时信号处理,在ARM上对性能要求比较高,所以需要比较快的FFT库 需要能同时支持定点化,NEON指令优化,并且支持非2的N次方的傅里叶变换库.希望有谁有用过的,了解过的,介绍一下给点指导. PS:现在找到大部分库都都不能满足上面三个特性,例如FFTWKISSFFTNET10 等都是只能满足两个特性.

图像处理-快速傅里叶变换蝶形公式如何理解

问题描述 快速傅里叶变换蝶形公式如何理解 // 采用蝶形算法进行快速傅里叶变换 for(k = 0; k < r; k++) { for(j = 0; j < 1 << k; j++) { bfsize = 1 << (r-k); for(i = 0; i < bfsize / 2; i++) { p = j * bfsize; X2[i + p] = X1[i + p] + X1[i + p + bfsize / 2]; X2[i + p + bfsize /

如果看了这篇文章你还不懂傅里叶变换,那就过来掐死我吧(下)

上一篇文章发出来之后,为了掐死我,大家真是很下工夫啊,有拿给姐姐看的,有拿给妹妹看的,还有拿给女朋友看的,就是为了听到一句"完全看不懂啊".幸亏我留了个心眼,不然就真的像标题配图那样了.我的文章题目是,如果看了这篇文章你"还"不懂就过来掐死我,潜台词就是在你学了,但是没学明白的情况下看了还是不懂,才过来掐死我. 另外,想跟很多人抱歉,因为评论太多了,时间有限,不能给每个人回复,还望大家谅解.但是很感谢一直在评论区帮忙解答读者问题的各位,就不一一@了. 这里郑重感谢大

跪求-图像增强关于二维离散快速傅里叶变换

问题描述 图像增强关于二维离散快速傅里叶变换 matlab7环境下自己编程实现fft2函数,及ifft2函数,不能调用工具箱函数,求大神们帮忙解决,分享一下 解决方案 http://wenku.baidu.com/link?url=N3wlVVKKz4uoRNzyV8BlnXi4QiV72zCiNE9oMHyis0z4NyZsP4MfC-AXVVCrgnHoWgDNi7Zw8SVqh8wipdmqZNuFUOG65N2K1rck0xXnm3i

快速傅里叶变换FFT的C语言算法彻底研究 看不到插图

问题描述 快速傅里叶变换FFT的C语言算法彻底研究 看不到插图 http://blog.csdn.net/wangyoufeng8889/article/details/8867380 这篇文章里边的插图怎么看不到啊?

图像处理-如何使用opencv的Java接口中的Core.dft()函数来实现图像的傅里叶变换?

问题描述 如何使用opencv的Java接口中的Core.dft()函数来实现图像的傅里叶变换? opencv中的dft()函数网上只找到了c++接口的,部分代码: Mat padded; //expand input image to optimal size int m = getOptimalDFTSize( I.rows ); int n = getOptimalDFTSize( I.cols ); // on the border add zero values copyMakeBor

HDU 1402 快速傅里叶变换FFT

题意:大数乘法. 看了算法导论(第2版)第30章的<多项式与快速傅里叶变换> 多项式有系数表示法和点值表示法. 两个次数界为n的多项式A(x)和B(x)相乘,输入输出均采用系数表示法.(假定n为2的幂) 1)使次数界增加一倍:A(x)和B(x)扩充为次数界为2n的多项式,并构造起系数表示 2)求值:两次应用2n阶FFT,计算出A(x)和B(x)的长度为2n的点值表示 3)点乘:计算多项式C(x)=A(x)*B(x)的点值表示 4)插值:对2n个点值对应用一次FFT计算出其逆DFT,就可以构造出

HDU 4609 快速傅里叶变换

题意:给出10^5个边长,问任取三个组成三角形的概率. 这题没遇到还不知道什么是傅里叶变换FFT,当时怎么想时间复杂度都不能降到n*log(n).实际是这样的,把每个记录长度个数的数组存在num数组中. 直接摘大神的吧 网址http://www.cnblogs.com/kuangbin/archive/2013/07/24/3210565.html 学会了FFT. 这题就很容易了. 其实题目是给了n条线段.问随机取三个,可以组成三角形的概率. 其实就是要求n条线段,选3条组成三角形的选法有多少种

神作:深入浅出傅里叶变换

一.什么是频域 从我们出生,我们看到的世界都以时间贯穿,股票的走势.人的身高.汽车的轨迹都会随着时间发生改变.这种以时间作为参照来观察动态世界的方法我们称其为时域分析.而我们也想当然的认为,世间万物都在随着时间不停的改变,并且永远不会静止下来.但如果我告诉你,用另一种方法来观察世界的话,你会发现世界是永恒不变的,你会不会觉得我疯了?我没有疯,这个静止的世界就叫做频域. 先举一个公式上并非很恰当,但意义上再贴切不过的例子: 在你的理解中,一段音乐是什么呢? 这是我们对音乐最普遍的理解,一个随着时间