用matlab实现数字图像的傅里叶变换_如何用一句诗,解释傅里叶(Fourier)变换?

96fb52b600a92f00569d4b523e155609.png

1768年,著名数学家傅里叶出生于法国欧塞尔,30年后提出傅里叶级数。

a3114183bffb4682b82ef0f060a5a591.png

1965年, FFT(快速傅里叶变换)问世,被评为二十世纪最伟大的几个算法之一(接近170年,终于可用了)。

1084年,苏东坡在江西九江写了一首诗---《题西林璧》。

他们之间究竟有什么联系呢?

......

先放一放,穿越回来,听我发发牢骚:


时至今日,在人工智能着火之前,有多少人仅凭一把傅里叶刀,就横行于各大高科技街头巷尾?

通信、图像处理、语音合成、音乐制作、CT扫描、核磁共振......但凡看起来比较复杂,需要去燥、压缩、提纯的地方,都是他们挥刀横舞的江湖,他们拿着高薪,说着“暗语”,嘘:

- 时频分析
- 时域到频域的变换
- 正交、余弦、卷积、频谱、相位、欧拉表达
- 天王盖地虎,宝塔镇河妖.....

才听得懂?这些该死的家伙!

所以,他们拿着高薪,不似刀客像剑(贱)客!

后来,留学的人多了,能看懂洋文的人也越来越多了。

韩昊同学因对教材恨之入骨,出去搞懂之后回来一篇傅里叶分析之掐死教程不知灌顶了多少迷途少年。

当年,我也始终没能理解时频变换这“黑话”,一心只盯着傅里叶变换的绝对可积条件瞎琢磨,梦想着有一天,凭此一步登天,青史名留 (据说这玩儿至今未解)。

后来在GE,有同事想要从数学角度理解傅里叶变换(其实他已会应用而我不会), 为此我专门做了个步步推导的详细说明Fourier Notes分享给大家,获得过些好评,似乎让人窥探到了图像处理的“核心机密”,可算本人在GE所做的少有几件有意思事情之一。

数学公式总是让大多数人望而生畏,但真的很想把这种美传递给更多的人,所以对这些“黑话”,要不厌其烦的重复解释,专业读者觉得啰嗦的可跳过(见谅)。

时频分析: 指时间频率两个维度的分析,

  • 时间很好理解,符合人的习惯思维,一段语音,第一秒、第二秒、第三秒......分别是什么?可映射成数值,画成函数图像,是一段非常复杂的曲线,这种分析叫按时间序列分析。
  • 单说频率,也好理解,一分钟跳绳您能跳几下,或者跳一次要用多长时间?国际单位是赫兹(HZ),指事件一秒内重复发生的次数。Do、Re、Mi的频率分别是261.6HZ、293.6HZ、329.6HZ 。

频率时间又有什么联系呢?好像互为倒数,没错!一般这样理解都默认单位次单位时,但是对同一事件而言。

傅里叶变换呢?它是指“等价”事件的两个不同观察视角互换,不是同一事件的简单互倒。

举个栗子:

  • 对于一首歌,如果从时间上采样得到了非常复杂的图像,如下:

932b8022c536a14f3a56144f04ece3bb.png

音乐家可能大概知道,它什么时间是什么调(Do、Re、Mi),但我一学数学的五音不全人士(KTV里从来就是陪笑的主),怎么能知道它是什么调呢?来,傅里叶变换大刀一挥,掐指一算,我可能比最顶级的歌唱家还更精准的知道它的调调!如下图:

763afb5e10c91aa4c2272d46042f1e18.png

这图什么意思?整首歌最多的部分在50HZ左右,其它各个频率的占有量一目了然,频率在261.6HZ的Do发音,没多少! 而知道了这个又有什么用呢?假如您唱歌不太上调或录音时有些杂音,通过傅里叶变换一看,锁定500HZ是真凶,260HZ欠缺,直接砍掉、添加还原,再听,杂音消失,歌王再现!放大一想,美颜也是同样的道理(自己骗自已),再深想,您的世界观就要发生变化了.......比如

  • 任何材料,只要您能让它发生261.6HZ的振荡,它必发出一个Do的音,就像唱歌时喉咙发声一样,至于如何做到精准频率控制,那是搞物理和材料的事......一旦找准,人还能唱过机器吗?焦虑,捶足顿胸! 原来数学音乐竟有如此曼妙关系,悔不当初没练练嗓子!

57249fa27ed445eb50573d5ba1dc3601.png

时域到频域的变换:有人喜欢这样说(装),所谓,可能来自数学的实数域、有理数域、群、环、域叫法,其实就是指一集合,让人费解,和时频分析是一个意思!

正交: 是相互垂直的高维扩展。

正(余)弦:是初中学的

,数值上它可表示为直角三形
直角边斜边的比,而
则可表示两种不同
频率的波,这是理解 傅里叶变换意义最为关键的地方。

卷积:把它理解为一种操作(物理算子),看数学公式别怕:

它只是一个用积分形式定义的函数,基于两个其它函数。

频谱:就是以频率为横轴画出来的图像,还有人还喜欢叫连续谱(装),就是说频率是连续的,频率还有无理数的吗?一想就懵!别想多了~。

相位:指正弦波的位置,相对零点的平移。

欧拉表达式:用公式

替换,是让表达更简捷的一种表示,这需要点耐心去研究,而看懂之后您
一定会跪倒在欧拉肖像面前(我保证),他告诉人们,这个世界都是各种不同半径的圆圈圈在以不同的频率不停的旋转。

好,有了这些“黑话”的概念之后,您再去看数学公式,网上搜或直接看我的Fourier Notes,应该就不会太难了,推导一遍,看看会不会有世界观的刷新。


若有豁然开朗,喜悦难于言表之情,那就太好了,先恭喜您。

好,平静之后,再来看看900多年前,苏轼在我的老家庐山留下的诗:

横看成岭侧成峰,远近高低各不同。

不识庐山真面目,只缘身在此山中。

光看第一句,横看成岭侧成峰

就感觉我前面说的那一堆

全是废话

废话......!

有没有这种感觉?

苏轼怕不是傅里叶的灵魂先知吧?

严重怀疑是傅里叶到过庐山,剽窃了老祖宗的学术成果才提出的傅里叶变换!

(别说全诗了,只第一句就是傅里叶变换的全部, 全部!)

数学家的那一堆零乱的公式,

诗人寥寥几字就说清楚了!

真是天外有天,

神外有神!

......


数学之水,更多干货,敬请点击详查~

0b5f2fe319c6f644552171df34f105fd.png
数学之水

文章同步微信公众号:数学之水,欢迎关注指导。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值