從時(shí)域和頻域來(lái)解析傅里葉變換(含代碼和性質(zhì))

冒泡~十二月啦!艱難的十一月總算是熬過(guò)去,十二月希望少熬一點(diǎn)吧(不可能)。
【題外話:傅里葉變換講道理應(yīng)該是大一高數(shù)就學(xué)習(xí),然而當(dāng)時(shí)的老師因?yàn)榭荚嚥豢季头艞壛私虒W(xué),于是乎現(xiàn)在的我學(xué)起來(lái)真是惡補(bǔ)的痛苦】

傅里葉變換

傅里葉變換實(shí)質(zhì)涉及的是頻域函數(shù)和時(shí)域函數(shù)的轉(zhuǎn)換。

概念解釋

a.先引入時(shí)域和頻域這兩個(gè)概念的解釋。

時(shí)域

時(shí)域是真實(shí)世界,是惟一實(shí)際存在的域。可以這樣理解,從我們出生開(kāi)始,所接觸的這個(gè)世界就是隨著時(shí)間在變化的,是在運(yùn)動(dòng)的。

頻域

頻域它不是真實(shí)的,而是一個(gè)數(shù)學(xué)構(gòu)造。如果說(shuō)時(shí)域是惟一客觀存在的域,那么頻域是一個(gè)遵循特定規(guī)則的數(shù)學(xué)范疇,頻域也被一些學(xué)者稱為上帝視角。結(jié)合上面對(duì)時(shí)域的理解,如果時(shí)域是運(yùn)動(dòng)永不停止的,那么頻域就是靜止的。
正弦波是頻域中唯一存在的波形,這是頻域中最重要的規(guī)則,即正弦波是對(duì)頻域的描述,因?yàn)轭l域中的任何波形都可用正弦波合成。
通過(guò)圖片來(lái)直觀解釋:
在時(shí)域里面,一段音樂(lè)是什么?是一個(gè)隨著時(shí)間變化的震動(dòng)(我們可以觀察到鋼琴的琴弦一會(huì)上一會(huì)下的擺動(dòng))。
相比較,在頻域里面,一段音樂(lè)又是什么?是一個(gè)個(gè)音符,是樂(lè)譜。音符的個(gè)數(shù)是有限且固定的,但可以組合出無(wú)限多的樂(lè)曲。

參考來(lái)源:https://zhuanlan.zhihu.com/p/19763358

b.接下來(lái)再講解兩個(gè)概念 :頻譜和相位譜

頻譜

在這里先引入一個(gè)通俗的解釋:


如上圖所示:
第一幅圖是一個(gè)余弦波 cos(x)
第二幅圖是 2 個(gè)余弦波的疊加 cos (x) +a.cos (3x)
第三幅圖是 4 個(gè)余弦波的疊加
第四幅圖是 10 個(gè)余弦波的疊加
那么問(wèn)題來(lái)了:
隨著正弦波數(shù)量逐漸的增長(zhǎng),他們最終會(huì)疊加成一個(gè)標(biāo)準(zhǔn)的矩形,那么可以推出什么?
不僅僅是矩形,你能想到的任何波形都是可以如此方法用正弦波疊加起來(lái)的。

接下來(lái)我們換一個(gè)角度來(lái)看看:
從側(cè)面看就是矩形波在頻域的樣子,是不是完全認(rèn)不出來(lái)了?教科書(shū)一般就給到這里然后留給了讀者無(wú)窮的遐想,以及無(wú)窮的吐槽,其實(shí)教科書(shū)只要補(bǔ)一張圖就足夠了:頻域圖像,也就是俗稱的頻譜。

對(duì)比展示:

image.png

可以看出,從時(shí)域來(lái)看,我們會(huì)看到一個(gè)近似為矩形的波,而我們知道這個(gè)矩形的波可以被差分為一些正弦波的疊加。而從頻域方向來(lái)看,我們就看到了每一個(gè)正弦波的幅值,可以發(fā)現(xiàn),在頻譜中,偶數(shù)項(xiàng)的振幅都是0,也就對(duì)應(yīng)了圖中的彩色直線。振幅為 0 的正弦波。
正面--時(shí)域圖像 側(cè)面--頻域圖像
動(dòng)圖展示:

相位譜

頻譜只代表了一個(gè)正弦函數(shù)的幅值,而要準(zhǔn)確描述一個(gè)正弦函數(shù),我們不僅需要幅值,還需要相位,不同相位決定了波的位置,所以對(duì)于頻域分析,僅僅有頻譜(振幅譜)是不夠的,我們還需要一個(gè)相位譜。
頻譜的重點(diǎn)是側(cè)面看,相位譜的重點(diǎn)則是從下面看。


如上圖所示:投影點(diǎn)我們用粉色點(diǎn)來(lái)表示,紅色的點(diǎn)表示離正弦函數(shù)頻率軸最近的一個(gè)峰值,而相位差就是粉色點(diǎn)和紅色點(diǎn)水平距離除以周期。將相位差畫(huà)到一個(gè)坐標(biāo)軸上就形成了相位譜。
【這里需要糾正一個(gè)概念:時(shí)間差并不是相位差。如果將全部周期看作2Pi或者360度的話,相位差則是時(shí)間差在一個(gè)周期中所占的比例。我們將時(shí)間差除周期再乘2Pi,就得到了相位差。】

傅里葉變換

為什么需要傅里葉變換?
如下圖所示:上面我們能看到的僅僅是一個(gè)類似正弦波的波形,其幅值在按照一定的規(guī)律變化。如何記載這個(gè)波形的信息呢?尤其是量化的記載!是很困難的。那么這個(gè)時(shí)候引入傅里葉變換就可以得到一個(gè)頻譜(幅值譜),主要包括3、5、7、9次諧波,一目了然!


傅里葉變換是一種信號(hào)分析方法,讓我們對(duì)信號(hào)的構(gòu)成和特點(diǎn)進(jìn)行深入的、定量的研究。
把信號(hào)通過(guò)頻譜的方式(包括幅值譜、相位譜和功率譜)進(jìn)行準(zhǔn)確的、定量的描述。
這就是傅里葉變換的主要目的。
(具體分析可參考:(https://www.sohu.com/a/224533027_99907714))

公式定義:


可以這樣理解,
①為頻率域的函數(shù)F(ω)表示為時(shí)間域的函數(shù)f(t)的積分形式
②為時(shí)間域的函數(shù)f(t)表示為頻率域的函數(shù)F(ω)的積分

傅里葉變換過(guò)程圖解:


具體看參考(https://blog.csdn.net/qq_33414271/article/details/79117586

其他性質(zhì)
參考:(https://www.cnblogs.com/AndyJee/p/5073281.html)

參考:百度百科

Matlab實(shí)現(xiàn)(頻譜和相位譜)

%% 二維傅立葉變換
%讀取圖片
I=imread('Naruto.jpg');
img=rgb2gray(I);
%幅值圖和相位譜
%進(jìn)行傅立葉變換
f=fft2(img);        %傅里葉變換
f=fftshift(f);      %使圖像對(duì)稱,中心化
r=real(f);          %圖像頻域?qū)嵅?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~在自我理解的基礎(chǔ)上完成了資料的搬運(yùn)是很不容容易了!
十二月加油鴨!

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書(shū)系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

  • 定義 一、傅里葉級(jí)數(shù) 法國(guó)數(shù)學(xué)家傅里葉發(fā)現(xiàn),任何周期函數(shù)都可以用正弦函數(shù)和余弦函數(shù)構(gòu)成的無(wú)窮級(jí)數(shù)來(lái)表示(選擇正弦函...
    MiracleJQ閱讀 6,386評(píng)論 1 7
  • 作 者:韓 昊知 乎:Heinrich 微 博:@花生油工人 知乎專欄:與時(shí)間無(wú)關(guān)的故事 謹(jǐn)以此文獻(xiàn)給大連海事大學(xué)...
    泉泉醬閱讀 735評(píng)論 0 1
  • 最近圈里朋友們興起一股寫(xiě)作熱,而好友每日一篇堅(jiān)持近一年累計(jì)10萬(wàn)以上的字?jǐn)?shù)更是鞭策了我!休養(yǎng)4年的腦子該動(dòng)起...
    Tracy_H閱讀 274評(píng)論 3 1
  • “立春一日,百草回芽”,盡管嚴(yán)寒的日子并未離去,但溫暖不會(huì)太遠(yuǎn)。一候東風(fēng)解凍、二候蟄蟲(chóng)始振、三候魚(yú)陟負(fù)冰,大自然萬(wàn)...
    銘印船長(zhǎng)閱讀 346評(píng)論 0 0
  • iOS切圓角的方式有三種 1.通過(guò)設(shè)置layer的屬性 最簡(jiǎn)單的一種,但是很影響性能,一般在正常的開(kāi)發(fā)中使用很少....
    weicyNO_1閱讀 2,240評(píng)論 0 2

友情鏈接更多精彩內(nèi)容