傅里叶变换和频域滤波基础

最新推荐文章于 2024-07-05 09:36:51 发布

Ian@_@

最新推荐文章于 2024-07-05 09:36:51 发布

阅读量667

点赞数

分类专栏：数学文章标签：数学频域卷积

本文链接：https://blog.csdn.net/weixin_42660942/article/details/115050289

版权

                    
                        
                    
                    傅里叶变换和频域滤波基础 
 
1. 傅里叶变换 
1.1 从傅里叶级数到傅里叶变换的推导 
欧拉公式： e j θ = c o s θ + j s i n θ e^{j\theta}=cos\theta+jsin\theta ejθ=cosθ+jsinθ
从傅里叶级数到傅里叶变换： 
  利用欧拉公式，将傅里叶级数从三角函数形式化为指数形式
利用黎曼和，将其中的累加形式化为积分形式
 
傅里叶级数三角函数形式： 
   f ( t ) f(t) f(t)为具有周期T的周期函数，可以表示为 f ( t ) = a 0 2 + ∑ n = 1 ∞ ( a n c o s ( n ω t ) + b n s i n ( n ω t ) ) f(t)=\frac{a_0}{2}+\sum_{n=1}^{\infin}(a_n cos(n\omega t)+b_n sin(n\omega t)) f(t)=2a0​​+∑n=1∞​(an​cos(nωt)+bn​sin(nωt))
其中， ω = 2 π T , a n = 2 T ∫ − T 2 T 2 f ( t ) c o s ( n ω t )   d t , b n = 2 T ∫ − T 2 T 2 f ( t ) s i n ( n ω t )   d t \omega=\frac{2\pi}{T},a_n=\frac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)cos(n\omega t)\ dt,b_n=\frac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)sin(n\omega t)\ dt ω=T2π​,an​=T2​∫−2T​2T​​f(t)cos(nωt) dt,bn​=T2​∫−2T​2T​​f(t)sin(nωt) dt
 n = 0 n=0 n=0时， a 0 = 2 T ∫ − T 2 T 2 f ( t ) d t a_0=\frac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)dt a0​=T2​∫−2T​2T​​f(t)dt
 
三角函数形式到指数形式的推导： 
  根据欧拉公式 e j θ = c o s θ + j s i n θ e^{j\theta}=cos\theta+jsin\theta ejθ=cosθ+jsinθ有： c o s θ = e j θ + e − j θ 2 , s i n θ = e j θ − e − j θ 2 i = − j ⋅ e j θ − e − j θ 2 cos\theta=\frac{e^{j\theta}+e^{-j\theta}}{2},sin\theta=\frac{e^{j\theta}-e^{-j\theta}}{2i}=-j\cdot\frac{e^{j\theta}-e^{-j\theta}}{2} cosθ=2ejθ+e−jθ​,sinθ=2iejθ−e−jθ​=−j⋅2ejθ−e−jθ​
代入得 f ( t ) = a 0 2 + ∑ n = 1 ∞ ( a n c o s ( n ω t ) + b n s i n ( n ω t ) = a 0 2 + ∑ n = 1 ∞ ( a n e j n ω t + e − j n ω t 2 − j ⋅ b n e j n ω t − e − j n ω t 2 ) = a 0 2 + ∑ n = 1 ∞ ( a n − j b n 2 e j n ω t + a n + j b n 2 e − j n ω t ) f(t)=\frac{a_0}{2}+\sum_{n=1}^{\infin}(a_n cos(n\omega t)+b_n sin(n\omega t)\\=\frac{a_0}{2}+\sum_{n=1}^{\infin}(a_n \frac{e^{jn\omega t}+e^{-jn\omega t}}{2} -j\cdot b_n\frac{e^{jn\omega t}-e^{-jn\omega t}}{2})\\=\frac{a_0}{2}+\sum_{n=1}^{\infin}(\frac{a_n-jb_n}{2}e^{jn\omega t}+\frac{a_n+jb_n}{2}e^{-jn\omega t}) f(t)=2a0​​+∑n=1∞​(an​cos(nωt)+bn​sin(nωt)=2a0​​+∑n=1∞​(an​2ejnωt+e−jnωt​−j⋅bn​2ejnωt−e−jnωt​)=2a0​​+∑n=1∞​(2an​−jbn​​ejnωt+2an​+jbn​​e−jnωt)
令 c 0 = a 0 2 , c n = a n − j b n 2 , d n = a n + j b n 2 , c_0=\frac{a_0}{2},c_n=\frac{a_n-jb_n}{2},d_n=\frac{a_n+jb_n}{2}, c0​=2a0​​,cn​=2an​−jbn​​,dn​=2an​+jbn​​,则 f ( t ) = c 0 + ∑ n = 1 ∞ ( c n e j n ω t + d n e − j n ω t ) f(t)=c_0+\sum_{n=1}^{\infin}(c_ne^{jn\omega t}+d_ne^{-jn\omega t}) f(t)=c0​+∑n=1∞​(cn​ejnωt+dn​e−jnωt)
则 c 0 = 1 T ∫ − T 2 T 2 f ( t ) d t , c n = 1 T ∫ − T 2 T 2 f ( t ) ( c o s ( n ω t ) − j s i n ( n ω t ) ) d t = 1 T ∫ − T 2 T 2 f ( t ) e − j n ω t d t , d n = 1 T ∫ − T 2 T 2 f ( t ) ( c o s ( n ω t ) + j s i n ( n ω t ) ) d t = 1 T ∫ − T 2 T 2 f ( t ) e j n ω t d t c_0=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)dt,\\c_n=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)(cos(n\omega t)-jsin(n\omega t))dt=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-jn\omega t}dt,\\d_n=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)(cos(n\omega t)+jsin(n\omega t))dt=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{jn\omega t}dt c0​=T1​∫−2T​2T​​f(t)dt,cn​=T1​∫−2T​2T​​f(t)(cos(nωt)−jsin(nωt))dt=T1​∫−2T​2T​​f(t)e−jnωtdt,dn​=T1​∫−2T​2T​​f(t)(cos(nωt)+jsin(nωt))dt=T1​∫−2T​2T​​f(t)ejnωtdt
可以看出 d n = c − n , ∑ n = 1 ∞ d n e − j n ω t = ∑ n = 1 ∞ c − n e − j n ω t = ∑ − ∞ − 1 c n e j n ω t d_n=c_{-n},\sum_{n=1}^{\infin}d_ne^{-jn\omega t}=\sum_{n=1}^{\infin}c_{-n}e^{-jn\omega t}=\sum_{-\infin}^{-1}c_{n}e^{jn\omega t} dn​=c−n​,∑n=1∞​dn​e−jnωt=∑n=1∞​c−n​e−jnωt=∑−∞−1​cn​ejnωt
则 f ( t ) = c 0 + ∑ n = 1 ∞ ( c n e j n ω t + d n e − j n ω t ) = c 0 e j 0 ω t + ∑ n = 1 ∞ c n e j n ω t + ∑ n = 1 ∞ d n e − j n ω t = c 0 e j 0 ω t + ∑ n = 1 ∞ c n e j n ω t + ∑ n = − ∞ − 1 c n e j n ω t = ∑ n = − ∞ ∞ c n e j n ω t = 1 T ∑ n = − ∞ ∞ ∫ − T 2 T 2 f ( t ) e − j n ω t d t ⋅ e j n ω t f(t)=c_0+\sum_{n=1}^{\infin}(c_ne^{jn\omega t}+d_ne^{-jn\omega t})\\=c_0e^{j0\omega t}+\sum_{n=1}^{\infin}c_ne^{jn\omega t}+\sum_{n=1}^{\infin}d_ne^{-jn\omega t}\\=c_0e^{j0\omega t}+\sum_{n=1}^{\infin}c_ne^{jn\omega t}+\sum_{n=-\infin}^{-1}c_ne^{jn\omega t}\\=\sum_{n=-\infin}^{\infin}c_n e^{jn\omega t}\\=\frac{1}{T}\sum_{n=-\infin}^{\infin}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-jn\omega t}dt\cdot e^{jn\omega t} f(t)=c0​+∑n=1∞​(cn​ejnωt+dn​e−jnωt)=c0​ej0ωt+∑n=1∞​cn​ejnωt+∑n=1∞​dn​e−jnωt=c0​ej0ωt+∑n=1∞​cn​ejnωt+∑n=−∞−1​cn​ejnωt=∑n=−∞∞​cn​ejnωt=T

                