ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

关于傅里叶变换:时移与频移解读之解读

2022-01-29 19:31:11  阅读:669  来源: 互联网

标签:infty right 频移 t0 解读 信号 时移 omega left


简 介: 本文讨论了傅里叶变换中相互关联的三个性质,时移特性、频移特性和尺度特性。对于信号与系统的理解需要能够在经典理论、物理直观以及数值计算三个方面进行深入挖掘和理解,这样在未来的应用中才能够更好的透过现象看本质。

关键词 FFT时移特性频移特性尺度特性

背 景 目 录
Contents
时移特性 性质证明 物理解释 频移特性 性质证明 物理解释 尺度变换对应
的频率变化
变换总结

 

§00   景

  今天在 西土城山羊卷 的博客看到她在其博文 傅里叶变换:时移与频移性质解读|CSDN创作打卡 中谈到对于傅里叶变换中的时移和频移的理解。其中两个观点还是蛮有新意的。

 

§01 移特性


  里叶变换的时移特性表面了信号 f ( t ) f\left( t \right) f(t) 经过延迟时间 t 0 ,    t 0 ≥ 0 t_0 ,\,\,t_0 \ge 0 t0​,t0​≥0 之后,形成信号 f ( t − t 0 ) f\left( {t - t_0 } \right) f(t−t0​) 。那么它对应的傅里叶变换 F t 0 ( ω ) F_{t0} \left( \omega \right) Ft0​(ω) 与原始信号 f ( t ) f\left( t \right) f(t) 对应的频谱 F ( ω ) F\left( \omega \right) F(ω) 之间的关系为: F t 0 ( ω ) = F ( ω ) ⋅ e − j ω t 0 F_{t0} \left( \omega \right) = F\left( \omega \right) \cdot e^{ - j\omega t_0 } Ft0​(ω)=F(ω)⋅e−jωt0​

1.1 性质证明

  证明这个性质非常容易,只要在傅里叶变换的定义上,通过 变量替换 便可以在三步之内完成证明。

  根据傅里叶变换定义: F ( ω ) = ∫ − ∞ ∞ f ( t ) e − j ω t d t F\left( \omega \right) = \int_{ - \infty }^\infty {f\left( t \right)e^{ - j\omega t} dt} F(ω)=∫−∞∞​f(t)e−jωtdt
  那么对于: F t o ( ω ) = ∫ − ∞ ∞ f ( t − t 0 ) e − j ω t d t F_{to} \left( \omega \right) = \int_{ - \infty }^\infty {f\left( {t - t_0 } \right)e^{ - j\omega t} dt} Fto​(ω)=∫−∞∞​f(t−t0​)e−jωtdt
  利用: t − t 0 → l t - t_0 \to l t−t0​→l ,原来表达式可以写为: F t 0 ( ω ) = ∫ − ∞ + ∞ f ( l ) e − j ω ( l + t 0 ) d l F_{t0} \left( \omega \right) = \int_{ - \infty }^{ + \infty } {f\left( l \right)e^{ - j\omega \left( {l + t_0 } \right)} dl} Ft0​(ω)=∫−∞+∞​f(l)e−jω(l+t0​)dl = [ ∫ − ∞ ∞ f ( l ) e − j ω l d l ] ⋅ e − j ω t 0 = F ( ω ) ⋅ e − j ω t 0 = \left[ {\int_{ - \infty }^\infty {f\left( l \right)e^{ - j\omega l} dl} } \right] \cdot e^{ - j\omega t_0 } = F\left( \omega \right) \cdot e^{ - j\omega t_0 } =[∫−∞∞​f(l)e−jωldl]⋅e−jωt0​=F(ω)⋅e−jωt0​

1.2 物理解释

  从上面的公式可以看出两点:

  • 由于时间延迟,信号中所有频率分量都出现了一个负相位 因子: e − j ω t 0 e^{ - j\omega t_0 } e−jωt0​ ;而信号的幅度谱不变。
  • 每个频率的负相位与该频率成线性关系,比例为 − t 0 - t_0 −t0​ 。

1.2.1 特定信号举例

  下面给出一个由四个频率组成的振荡信号:
f ( t ) = f ( 10 t ) + f ( 12 t ) + f ( 15 t ) + f ( 18 t ) f\left( t \right) = f\left( {10t} \right) + f\left( {12t} \right) + f\left( {15t} \right) + f\left( {18t} \right) f(t)=f(10t)+f(12t)+f(15t)+f(18t)

  它对应的波形为:
▲ 图1.2.1 信号f(t)的四个频率分量

▲ 图1.2.1 信号f(t)的四个频率分量

from headm import *                 # =

fdim = [10,12,15,18]
t = linspace(0,2, 1000)

plt.clf()
plt.figure(figsize=(10,10))

plt.subplot(len(fdim)+1, 1, 1)
sumsin = zeros(len(t))
for f in fdim:
    sumsin = sumsin + sin(f*t)
plt.plot(t, sumsin, label='F:%0.1rad/s')
plt.xlabel("t")
plt.ylabel("sint")
plt.grid(True)
plt.tight_layout()

for id,f in enumerate(fdim):
    plt.subplot(len(fdim)+1, 1, id+2)
    plt.plot(t, sin(f*t), label='F:%0.1rad/s')
    plt.xlabel("t")
    plt.ylabel("sint")
    plt.grid(True)
    plt.tight_layout()

plt.savefig(r"d:\temp\figure1.jpg")
plt.close()
tspshowimage(image=r"d:\temp\figure1.jpg")

  下面假设信号平移了 t 0 = 0.1 t_0 = 0.1 t0​=0.1 ,变成了 t → t − 0.1 t \to t - 0.1 t→t−0.1 。那么对应的波形为:

▲ 图1.2.2 延时信号对应的波形及其各个频率分量

▲ 图1.2.2 延时信号对应的波形及其各个频率分量

  在上面图中,可以明显看到,原来信号延迟的时间 t 0 = 0.1 t_0 = 0.1 t0​=0.1 ,所有频率分量都延迟相同的时间 t 0 = 0.1 t_0 = 0.1 t0​=0.1 。在此基础上, 西土城山羊卷 给出了两点有趣的解释,她认真解释的很到位:

  • 为什么对应的延迟时间 t 0 ≥ 0 t_0 \ge 0 t0​≥0 时,所有的相位都是的?从上图中可以直观看到所有的正弦分量的起始点实际上都基本上对应的是负的相位。
  • 为什么不同频率对应的相位延迟不同?就是对于相同的时间 t 0 t_0 t0​ ,转换到相位时 θ 0 = t 0 ⋅ ω \theta _0 = t_0 \cdot \omega θ0​=t0​⋅ω 与频率成正比。

▲ 图1.2.3  几种不同的典型性好平移后对应的频谱和相位

▲ 图1.2.3 几种不同的典型性好平移后对应的频谱和相位

 

§02 移特性


  际上 西土城山羊卷 博文中有趣的观点在于对于傅里叶变换频移特性的阐述。

2.1 性质证明

  这是在傅里叶变换中与时移特性对偶的一个特性。它讲述了如果把信号 f ( t ) f\left( t \right) f(t) 的频谱 F ( ω ) F\left( \omega \right) F(ω) 在频域进行平移 ω 0 \omega _0 ω0​ ,变成 F ( ω − ω 0 ) F\left( {\omega - \omega _0 } \right) F(ω−ω0​) ,那么对应的时域信号变成: f ω 0 ( t ) = f ( t ) ⋅ e j ω 0 t f_{\omega 0} \left( t \right) = f\left( t \right) \cdot e^{j\omega _0 t} fω0​(t)=f(t)⋅ejω0​t

  对比一下前面的时移特性,你会发现频移特性几乎与时移特性一模一样,只是叙述的时域、频域颠倒了一下,另外在时域信号上增加的相位因子是的。

  证明过程也与前面时移特性一样,只需要从傅里叶反变换出发,利用变量替换也就可以了。

  根据傅里叶反变换: f ( t ) = 1 2 π ∫ − ∞ + ∞ F ( ω ) e j ω t d ω f\left( t \right) = {1 \over {2\pi }}\int_{ - \infty }^{ + \infty } {F\left( \omega \right)e^{j\omega t} d\omega } f(t)=2π1​∫−∞+∞​F(ω)ejωtdω

  那么频率平移之后,对应的信号:
f ω 0 ( t ) = 1 2 π ∫ − ∞ ∞ F ( ω − ω 0 ) e j ω t d ω f_{\omega 0} \left( t \right) = {1 \over {2\pi }}\int_{ - \infty }^\infty {F\left( {\omega - \omega _0 } \right)e^{j\omega t} d\omega } fω0​(t)=2π1​∫−∞∞​F(ω−ω0​)ejωtdω

  利用变量替换: ω − ω 0 → λ \omega - \omega _0 \to \lambda ω−ω0​→λ ,则: f ω 0 ( t ) = ∫ − ∞ + ∞ F ( λ ) e j ( λ + ω 0 ) t d λ f_{\omega 0} \left( t \right) = \int_{ - \infty }^{ + \infty } {F\left( \lambda \right)e^{j\left( {\lambda + \omega _0 } \right)t} d\lambda } fω0​(t)=∫−∞+∞​F(λ)ej(λ+ω0​)tdλ = [ ∫ − ∞ + ∞ F ( λ ) e j λ t d λ ] ⋅ e j ω 0 t = f ( t ) ⋅ e j ω 0 t = \left[ {\int_{ - \infty }^{ + \infty } {F\left( \lambda \right)e^{j\lambda t} d\lambda } } \right] \cdot e^{j\omega _0 t} = f\left( t \right) \cdot e^{j\omega _0 t} =[∫−∞+∞​F(λ)ejλtdλ]⋅ejω0​t=f(t)⋅ejω0​t

2.2 物理解释

  但是在实际中,对于频移特性的物理解释则是从时域的角度来阐明。即在时域中,信号 f ( t ) ⋅ e j ω 0 t f\left( t \right) \cdot e^{j\omega _0 t} f(t)⋅ejω0​t 不是被称为增加了一个相位因子,而是指信号 f ( t ) f\left( t \right) f(t) 对振荡信号 e j ω 0 t e^{j\omega _0 t} ejω0​t 进行 幅度调制 ,调制后的信号 f ( t ) ⋅ e j ω 0 t f\left( t \right) \cdot e^{j\omega _0 t} f(t)⋅ejω0​t 的频谱是原来信号 f ( t ) f\left( t \right) f(t) 的频谱 F ( ω ) F\left( \omega \right) F(ω) 进行平移。通常情况下,是频谱从低频搬移到高频。

  当然,也可以利用这个原理将信号的频谱从高频搬回到低频。从而实现了信号的调制与解调。

  在通常情况下,对于信号调制都是使用正弦波,而不是复振荡信号,因此对应的信号频谱就表现为左右平移的情况。如下图所示。

▲ 图2.2.1  信号频谱的搬移对应的信号调制过程

▲ 图2.2.1 信号频谱的搬移对应的信号调制过程

2.3 尺度变换对应的频率变化

  西土城山羊卷 在解释信号频谱搬移的时候提出了一个有趣的观点:你看看,信号的频谱搬移到高频了,所有的频率都增加了。为什么信号在时域表现不是尺度变化呢?

  也就是: f ( t ) → f ( a ⋅ t ) ,      a > 1 f\left( t \right) \to f\left( {a \cdot t} \right),\,\,\,\,a > 1 f(t)→f(a⋅t),a>1 。比如, f ( t ) = sin ⁡ ( ω t ) f\left( t \right) = \sin \left( {\omega t} \right) f(t)=sin(ωt) ,对应的 f ( 2 t ) = sin ⁡ ( 2 ω t ) f\left( {2t} \right) = \sin \left( {2\omega t} \right) f(2t)=sin(2ωt) ,这个频率也变高了。

  如果这样考虑,这就会引出傅里叶变换中的“尺度变化”性质。

  如果 一个信号在时域发生了尺度变化 f ( t ) → f ( a ⋅ t ) f\left( t \right) \to f\left( {a \cdot t} \right) f(t)→f(a⋅t) ,这个信号有可能被压缩( a > 1 a > 1 a>1 ),也有可能被拉长( a < 1 a < 1 a<1 ),这都会引起其中的频率分量的变化。但这个变化不再是平移,而是呈现相反的尺度变化:


  这个规律与频移是不同的。

  在一般的工程应用中,尺度变化实现起来不太容易,但在“空间换时间,时间换空间”的过程中,可以通过巧妙的设置,完成信号在时域和频域中的尺度变换。

▲ 图2.3.1  傅里叶变幻的尺度性质

▲ 图2.3.1 傅里叶变换的尺度性质

 

换总结 ※


  文讨论了傅里叶变换中相互关联的三个性质,时移特性、频移特性和尺度特性。对于信号与系统的理解需要能够在经典理论、物理直观以及数值计算三个方面进行深入挖掘和理解,这样在未来的应用中才能够更好的透过现象看本质。


■ 相关文献链接:

● 相关图表链接:

标签:infty,right,频移,t0,解读,信号,时移,omega,left
来源: https://blog.csdn.net/zhuoqingjoking97298/article/details/122746177

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有