matla图像b傅里叶逆变换_【学习/方法】尝试理解傅里叶变换

7ba6bb92b013550697b377e29932dd4e.png

注:本文只是对https://www.bilibili.com/video/av19141078/视频的个人笔记,侵权删

1cd446f472246487bb8278f88e7ce169.png

之前有看过几篇关于傅里叶变换和拉普拉斯变换的科普文。

66048c0cb63859a936b0ff938e1a41a6.png

是,这些文章讲了时域与频域的差别,讲了波叠加后的图像。

但看来看去,总觉得差了点什么,我拿出书本,看着那些公式,依旧不明白其意义,不明白为什么傅里叶变换偏偏就能把一个函数变成无数正弦波的叠加,为什么要有负无穷到正无穷的积分,为什么会有乘以一个e^-jwt?为什么会用冲激函数来表示傅里叶变换的结果?

于是我又打开的B站……

刷完视频,终于能够过自己那关了,虽然理解只是粗浅层面。

下面结合视频说说自己对傅里叶变换的理解:

首先需要知道傅里叶变换的结果,我们都知道傅里叶变换把任一的波形变成了无数正弦波的叠加,但这个说法依旧很抽象,我还是不知道他是怎么滤波的。但你看到下面这个声音波形,大概就能理解了:

53d84b7c8e091904843f9eae065682af.png

在音频处理软件中,如果要进行消除噪声或是实现高通低通的效果,就是利用傅里叶变换将声音进行分析,从上图看到,啸叫声在高频的位置,于是我们就能够直接把这段高频噪声消除,改善音质!

40084c28a1a26120945364dac3eac686.png

知道了傅里叶变换的这个应用,我想我们可以尝试取理解它了。我们现在把上面的波形图缠绕在一个圆上,有什么用等会你就知道了

8c15b3512f70d5d233e441e9211bf6d0.png

缠绕的频率会决定这个缠绕图像的样子

251f5751c3c42e90e8f74a756ebc9c31.png

37107f4ced36e5575371caf897dcf9b5.png

现在我们引入一个质心的概念,只记录x轴,然后就让绕在圆上 的图像随着频率的增加旋转起来吧,同时记录质心的位置

216ae31ef334ba2d0e6502ebf3bd6737.png

73fb285236372fc49eb163cf5b367ae1.png

b9022501baa2d964439b0543efaddb8e.png

(图中的虚线是缠绕周期,比如2cycle/second,你就在原波形图每0.5秒划一条线,这样就代表两条线直接的图像绕这下方圆一圈)

发现没有,只有当缠绕频率等于原本波形的频率的时候,这个质心的横坐标才会明显向右移动!

现在你是不是决定这个特性是不是大加利用,就像用共振来检验物体的固有频率一样,我们也可以通过这个缠绕来求出原波形的频率,可以从sint+sin2t的图像中把sint和sin2t分别抓出来!是不是很神奇,是不是觉得稍微懂了点

55ba9121a222ac8f555f2be397d4d5ab.png

(至于频率的0那点,其实只是代表波形的位置,我们只需要把整个波形向下移动,就可以消除0处那个尖峰)

43b0112c37a9ebc6368e40c9b0ac9a30.png

606805a6451b0f7b21e5a35fcaffed7e.png

好了,有了上面 的缠绕,你现在可以尝试理解真·傅里叶变换了,上面其实是一个近·傅里叶变换,毕竟我们总不可能老是把一个图形当成小仓鼠,在圈里绕来绕去吧,我们需要一个简洁的式子,来表示这个有点绕晕人的图形。

想到旋转,你能想到什么?

……

(给你点时间思考)

你知道吗,复数的本质就是旋转。

如果你想要一个一秒钟转一圈的圆,就可以用e2πit!但要一个任一频率转圈的圆呢?

那我们就可以用e2πfit,或者表示为ejwt。

是不是很眼熟!

好像差一个负号?

好,那就加上吧,e-jwt

(这是因为在傅里叶变换下默认是以顺时针旋转的,但复数的旋转矢量是逆时针。)

2337b5e33a8ec3e69bf0313f1f6d1b1d.png

现在我们就可以把一个函数乘以一个e-jwt来得到一个缠绕图像啦,这和上面我们缠的图本质上一样

77adf428959030cdd5c2e60d40996612.png

既然现在得到图了,我们就该想办法来求出这个图像的质心了,来把波形转换成频谱

质心的公式是什么,取很多个点,相加,然后求其平均

于是我们用积分表示

582236308105e92459e071a1f55b0500.png

49e9c211cb0521e3225de093285c9ecb.png

b68e025032d576428cd7585f893ad92a.png

好了,再看看书,书里的傅里叶正变换公式是什么

是不是只有右边的积分公式,而没有取平均

75f216345948e3f3a343e79ca9c6757c.png

那这样的真·傅里叶变换求出来的是什么,想一下,如果我们取一段3秒的波形,那积分式就是质心的x轴坐标的三倍,取6秒波形,那积分式就是质心的x轴坐标的6倍……

于是我们如果取无穷长的波形呢,那傅里叶变换的模长是不是就变得无穷大了?

5116fea16ce1240e8ae67d435b97fc02.png

话说,为什么要取无穷大呢,可以这样理解:时间越长,这样缠绕的图像就越接近中心,也就是除了等于波形频率外的几点周围其他会接近于0,而等于波形频率那点会是一个明显的峰

5c35d454264741555e1e3b175d19ee6e.png

75ace60609f5149c3e4e5fd273c9a23d.png

现在你是不是已经大概理解了这个公式有多吓人,用复数的指数形式乘以一个函数f(t),意味着画一张缠绕版的图像,至于复函数的积分则是利用质心的思想

现在回想起来是不是觉得这些符号变得异常直观!

217c7f1389a336f845bae44a46a6f76c.png

到最后,尝试解释一下为什么要引入一个单位冲激函数δ(t)

上面的分析我们已经知道了傅里叶变换公式的意义——它把一个无序的波形转化为频谱的,但是傅里叶变换的公式是对复函数进行积分欸,怎么算呀。

于是我们需要的是一个能简单显示频谱的工具,就像用音符来简单表示钢琴曲一样,而不是上来就把钢琴曲的声音波形图丢给你。

上面提到了,傅里叶变换的频谱中,除了等于波形频率外的几点周围其他会接近于0,而等于波形频率那点会是一个明显的峰,趋近于无穷大,现在回想单位冲激函数的定义,是不是就好像懂了?

41769f5e1635289f05b7b562d5b5258e.png

其实我们就是把那一个个波峰用单位冲激函数进行表示

75eea0501ced3d7d63843441aa7713d2.png

利用单位冲激函数的傅里叶变换只是1,得到其逆变换

66de9aa29b9ba6287d965f29a5c188a0.png

于是就可以非常简洁的表示傅里叶变换把这些函数的傅里叶变换转化为单位冲激函数的左右移动和叠加!

c9f3a904d354753e52c678bdb41a59e8.png

好了,我该去刷题准备考试了,不搞这些旁门左道了

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值