从时域和频域来解析傅里叶变换(含代码和性质)

冒泡~十二月啦!艰难的十一月总算是熬过去,十二月希望少熬一点吧(不可能)。
【题外话:傅里叶变换讲道理应该是大一高数就学习,然而当时的老师因为考试不考就放弃了教学,于是乎现在的我学起来真是恶补的痛苦】

傅里叶变换

傅里叶变换实质涉及的是频域函数和时域函数的转换。

概念解释

a.先引入时域和频域这两个概念的解释。

时域

时域是真实世界,是惟一实际存在的域。可以这样理解,从我们出生开始,所接触的这个世界就是随着时间在变化的,是在运动的。

频域

频域它不是真实的,而是一个数学构造。如果说时域是惟一客观存在的域,那么频域是一个遵循特定规则的数学范畴,频域也被一些学者称为上帝视角。结合上面对时域的理解,如果时域是运动永不停止的,那么频域就是静止的。
正弦波是频域中唯一存在的波形,这是频域中最重要的规则,即正弦波是对频域的描述,因为频域中的任何波形都可用正弦波合成
通过图片来直观解释:
在时域里面,一段音乐是什么?是一个随着时间变化的震动(我们可以观察到钢琴的琴弦一会上一会下的摆动)。
相比较,在频域里面,一段音乐又是什么?是一个个音符,是乐谱。音符的个数是有限且固定的,但可以组合出无限多的乐曲。

参考来源:https://zhuanlan.zhihu.com/p/19763358

b.接下来再讲解两个概念 :频谱和相位谱

频谱

在这里先引入一个通俗的解释:


如上图所示:
第一幅图是一个余弦波 cos(x)
第二幅图是 2 个余弦波的叠加 cos (x) +a.cos (3x)
第三幅图是 4 个余弦波的叠加
第四幅图是 10 个余弦波的叠加
那么问题来了:
随着正弦波数量逐渐的增长,他们最终会叠加成一个标准的矩形,那么可以推出什么?
不仅仅是矩形,你能想到的任何波形都是可以如此方法用正弦波叠加起来的。

接下来我们换一个角度来看看:
从侧面看就是矩形波在频域的样子,是不是完全认不出来了?教科书一般就给到这里然后留给了读者无穷的遐想,以及无穷的吐槽,其实教科书只要补一张图就足够了:频域图像,也就是俗称的频谱。

对比展示:

image.png

可以看出,从时域来看,我们会看到一个近似为矩形的波,而我们知道这个矩形的波可以被差分为一些正弦波的叠加。而从频域方向来看,我们就看到了每一个正弦波的幅值,可以发现,在频谱中,偶数项的振幅都是0,也就对应了图中的彩色直线。振幅为 0 的正弦波。
正面--时域图像 侧面--频域图像
动图展示:

相位谱

频谱只代表了一个正弦函数的幅值,而要准确描述一个正弦函数,我们不仅需要幅值,还需要相位,不同相位决定了波的位置,所以对于频域分析,仅仅有频谱(振幅谱)是不够的,我们还需要一个相位谱。
频谱的重点是侧面看,相位谱的重点则是从下面看。


如上图所示:投影点我们用粉色点来表示,红色的点表示离正弦函数频率轴最近的一个峰值,而相位差就是粉色点和红色点水平距离除以周期。将相位差画到一个坐标轴上就形成了相位谱。
【这里需要纠正一个概念:时间差并不是相位差。如果将全部周期看作2Pi或者360度的话,相位差则是时间差在一个周期中所占的比例。我们将时间差除周期再乘2Pi,就得到了相位差。】

傅里叶变换

为什么需要傅里叶变换?
如下图所示:上面我们能看到的仅仅是一个类似正弦波的波形,其幅值在按照一定的规律变化。如何记载这个波形的信息呢?尤其是量化的记载!是很困难的。那么这个时候引入傅里叶变换就可以得到一个频谱(幅值谱),主要包括3、5、7、9次谐波,一目了然!


傅里叶变换是一种信号分析方法,让我们对信号的构成和特点进行深入的、定量的研究。
把信号通过频谱的方式(包括幅值谱、相位谱和功率谱)进行准确的、定量的描述。
这就是傅里叶变换的主要目的。
(具体分析可参考:(https://www.sohu.com/a/224533027_99907714))

公式定义:


可以这样理解,
①为频率域的函数F(ω)表示为时间域的函数f(t)的积分形式
②为时间域的函数f(t)表示为频率域的函数F(ω)的积分

傅里叶变换过程图解:


具体看参考(https://blog.csdn.net/qq_33414271/article/details/79117586

其他性质
参考:(https://www.cnblogs.com/AndyJee/p/5073281.html)

参考:百度百科

Matlab实现(频谱和相位谱)

%% 二维傅立叶变换
%读取图片
I=imread('Naruto.jpg');
img=rgb2gray(I);
%幅值图和相位谱
%进行傅立叶变换
f=fft2(img);        %傅里叶变换
f=fftshift(f);      %使图像对称,中心化
r=real(f);          %图像频域实部
i=imag(f);          %图像频域虚部
margin=log(sqrt(r.^2+i.^2));      %图像幅度谱,加log便于显示
phase=real(angle(f)*180/pi);     %图像相位谱
figure
subplot(1,3,1),imshow(img),title('源图像');
subplot(1,3,2),imshow(margin,[]),title('图像幅度谱');
subplot(1,3,3),imshow(phase,[]),title('图像相位谱');

参考(https://blog.csdn.net/qq_33414271/article/details/79117586
https://ww2.mathworks.cn/help/matlab/math/fouriertransforms.html#responsive_offcanvas

Ending~在自我理解的基础上完成了资料的搬运是很不容容易了!
十二月加油鸭!

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 215,723评论 6 498
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,003评论 3 391
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 161,512评论 0 351
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,825评论 1 290
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,874评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,841评论 1 295
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,812评论 3 416
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,582评论 0 271
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,033评论 1 308
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,309评论 2 331
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,450评论 1 345
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,158评论 5 341
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,789评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,409评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,609评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,440评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,357评论 2 352

推荐阅读更多精彩内容

  • 定义 一、傅里叶级数 法国数学家傅里叶发现,任何周期函数都可以用正弦函数和余弦函数构成的无穷级数来表示(选择正弦函...
    MiracleJQ阅读 6,074评论 1 7
  • 作 者:韩 昊知 乎:Heinrich 微 博:@花生油工人 知乎专栏:与时间无关的故事 谨以此文献给大连海事大学...
    泉泉酱阅读 663评论 0 1
  • 最近圈里朋友们兴起一股写作热,而好友每日一篇坚持近一年累计10万以上的字数更是鞭策了我!休养4年的脑子该动起...
    Tracy_H阅读 215评论 3 1
  • “立春一日,百草回芽”,尽管严寒的日子并未离去,但温暖不会太远。一候东风解冻、二候蛰虫始振、三候鱼陟负冰,大自然万...
    追梦人船長阅读 269评论 0 0
  • iOS切圆角的方式有三种 1.通过设置layer的属性 最简单的一种,但是很影响性能,一般在正常的开发中使用很少....
    weicyNO_1阅读 2,170评论 0 2