最近在知乎上发现了一篇超厉害的关于傅里叶变换、拉普拉斯变换和Z变换的 帖子,赶紧整理一下学习笔记~作者是DBinary
傅里叶牛逼!
- 他说
- “任何”周期信号都可以用一系列成谐波关系的正弦曲线表示
- 等会,啥是谐波呀?
- 首先介绍一下基波:
- 基波是指在复杂的周期性振荡中与该振荡最长周期相等的正弦波分量
- 大概的意思就是在傅里叶变换后周期最长的那个分量波
- 百百说的是
- 谐波是指电流中所含有的频率为基波的整数倍的电量
- 一般是指对周期性的非正弦电量进行傅立叶级数分解,除了基波频率的电量,其余大于基波频率的电流产生的电量,称为谐波。
- 首先介绍一下基波:
- 当然,这需要满足一定的条件——狄里赫利条件
- 函数在任意有限区间内连续,或只有有限个第一类间断点(当t从左或右趋于这个间断点时,函数有有限的左极限和右极限)
- 在一个周期内,函数有有限个极大值或极小值
- x(t)在单个周期内绝对可积,即 ∫ 0 T ∣ x ( t ) ∣ d t < ∞ \int_0^T|x(t)|dt<\infty ∫0T∣x(t)∣dt<∞
这个结论超牛的!这说明下面一些奇形怪状的周期信号最后都能用一个式子表示
-
作者还用了一个gif图表示:
-
-
意义:如果可以将信号分解为正弦函数的累加和,就等于知道了这个信号是由哪些频率的正弦波构成,那么,也能知道对应频率的波在信号中的能量和相位信息了
-
所以,如何将信号 f ( t ) f(t) f(t)分解成的形式呢?
先介绍一下三角函数的正交性:
三角函数的正交性
- 对正弦波从 + ∞ 到 − ∞ +\infty到-\infty +∞到−∞积分=0,这里的 ∞ \infty ∞默认为 lim n → ∞ n T \lim\limits_{n\to \infty}nT n→∞limnT
- 不同频率的正弦余项波相乘的无穷积分也=0
- 这又叫做波的相干性,在信号里有重要应用
- 如果你想检验信号里有没有频率为100HZ的正弦波信号,只需将这个信号与100HZ的正弦波相乘,然后在周期内积分,若 ≠ 0 \ne 0 =0则信号中有含100HZ的信号
以上方法可以确定频率f了,那么如何求公式中的 c n 和 φ c_n和\varphi cn和φ呢?
- 这里才发现所有的 φ \varphi φ都是一个值哎!
求 c n 和 φ c_n和\varphi cn和φ
-
首先利用三角函数的变换公式,上式可变成
-
令 a n = c n cos φ , b n = − c n sin φ a_n=c_n\cos\varphi,b_n=-c_n\sin\varphi an=cncosφ,bn=−cnsinφ
-
上式变成
-
现在确定 a n , b n a_n,b_n an,bn,先求 b n b_n bn好了:只需让 f ( t ) ∗ sin ( n w t ) f(t)*\sin(nwt) f(t)∗sin(nwt)再积分, 即
-
同样的办法可得: a n = 2 T ∫ 0 T f ( t ) cos ( n w t ) d t a_n=\frac2T\int_0^Tf(t)\cos(nwt)dt an=T2∫0Tf(t)cos(nwt)dt
-
所以
-
令人疑惑的是,这与傅里叶的标准公式好像有点不大一样
-
这是欧拉公式呀,小笨猪
-
将欧拉代入,得
-
多出来的复数就是为了在计算里将sin与cos区分开,没有别的用处
-
是不是感觉好像还是不大对,那是因为之前那个式子是求和,这里只有一个,如果想变成求和的形式,要好多个 F ( f ( t ) ) F(f(t)) F(f(t))加起来,乘的 e − i w t e^{-iwt} e−iwt就变成对应的 e − i k w t e^{-ikwt} e−ikwt
拉普拉斯变换
- 公式:
- 发现,这就是在原傅里叶级数的基础上×了个 e − σ e^{-\sigma} e−σ呀!那为啥要这样啊
- 像下面这种增长太快准备上天的函数这怎么用上天的正弦函数去拟合呀,高人说只要“掰弯”它即可哈哈哈哈哈哈,×个 e − σ e^{-\sigma} e−σ之后就成这熊样了
- 不过为啥呀,为啥×了之后就掉下来了?还是说本来就是要掉下来的,只不过×上 e − σ e^{-\sigma} e−σ之后加快了它掉下来的速度?难道这,就是 e − σ e^{-\sigma} e−σ又称为衰减因子的原因?
离散型傅里叶变换
- 以上都是说的连续型傅里叶变换,当你需要在计算机中处理信号时就会发愁了,计算机只能将连续型采点成离散型,那么离散型傅里叶变化是啥样的呢?
- 再回忆一下连续型
- 要变成离散的,首先, f ( t ) f(t) f(t)就要变成 x [ n ] x[n] x[n](对应的第n个采样点)
- 积分要变成求和
- e − i w t e^{-iwt} e−iwt中的t也要变成n
- 还有 w w w这个动态的角速度值换成具体的角度 ϕ \phi ϕ?这里有点不太懂
- 反正最后离散型傅里叶公式就变成这样了
- 再令 Z = A e j ϕ Z=Ae^{j\phi } Z=Aejϕ,上式可写成 ∑ − ∞ ∞ x [ n ] Z − n \sum\limits_{-\infty}^{\infty}x[n]Z^{-n} −∞∑∞x[n]Z−n
- 妈耶,原来这玩意就是Z变换,我当时看着就头大所以拖到现在都没看的Z变换呀!