几种傅里叶变换之间的关系

几种傅里叶变换之间的关系

首先是傅里叶级数,一般来说,傅里叶级数用于分析周期信号,可以将信号表示成正弦的求和级数:
s N ( x ) = A 0 2 + ∑ n = 1 N A n ⋅ sin ⁡ ( 2 π n x P + ϕ n ) , for integer  N   ≥   1. s_{N}(x)={\frac {A_{0}}{2}}+\sum _{{n=1}}^{N}A_{n}\cdot \sin({\tfrac {2\pi nx}{P}}+\phi _{n}),\quad \scriptstyle {\text{for integer}}\ N\ \geq \ 1. sN(x)=2A0+n=1NAnsin(P2πnx+ϕn),for integer N  1.

也可以使用欧拉公式进而表示成指数函数的求和级数;如果是非周期的,那么就需要在数学意义上进行延拓,并使其满足一定的条件(迪利克雷)才可以计算傅里叶级数,但是似乎信号处理中探讨的不多。但是可能最多用的地方,就是考虑一个周期信号,可以通过傅里叶级数将其分解为有限个指数函数(正余弦函数)乘以不同大小的权重的相加的结果,就是下图这种效果(图来自wikipeida)。
在这里插入图片描述

其次是傅里叶变换,这就又分为离散和连续的两种,之间的换算关系如下图。
在这里插入图片描述

  • 连续时间傅里叶变换(CTFT) X a ( j Ω ) = ∫ x a ( t ) e − j Ω t d t X_a(j\Omega)=\int x_a(t) e^{-j\Omega t}dt Xa(jΩ)=xa(t)ejΩtdt。值得注意的是,这里分析的是连续信号,所以频率采用 Ω \Omega Ω表示,表示角频率,单位 弧度/秒。要想求得信号的CTFT,须满足迪利克雷条件。CTFT其实可以被认为是最接近原本傅里叶为了解决热力学问题而提出的那个运算定义的,其他的各种傅里叶变换都是不同程度的变体。换句话说,CTFT就是老老实实的通过公式,对函数进行积分然后得到计算结果。如果说这玩意它有什么用,只能说它其实没啥用(不是真的没用,只是信号处理里一般不会去计算这个),就只是一种运算的定义,和 卷积 是一个道理。当然也就是这种运算方式,让我们得到了连续信号在频域中的一种描述,我们就可以通过这个运算的结果来对信号进行相关的频域分析。
    同时这也满足一些经验性的结论:
    1. 时域连续非周期,频域连续
    2. 时域连续周期,频域离散
  • 离散时间傅里叶变换(DTFT) X ( e j ω ) = ∑ x [ n ] e − j ω n X(e^{j\omega})=\sum x[n]e^{-j\omega n} X(ejω)=x[n]ejωn。值得注意的是,这里分析的是离散信号 x [ n ] x[n] x[n],所以频率采用小写的数字角频率 ω \omega ω表示。这个数字角频率和CTFT中的模拟角频率之间有一个对应关系: ω = Ω / T s \omega=\Omega/T_s ω=Ω/Ts(Ts是采样周期),可以发现,在这个过程中,数字角频率失去了量纲,意味着这是个没有量纲的物理量,它和真实频率(Hz)之间的关系是通过采样率来传递的。对于给定的数字角频率 ω \omega ω值,采样率越大,其对应到真实世界中的信号频率就越大,这一点要怎么理解?我觉得可以从数字离散信号的时域形式下手,因为其时域形式将 n n n作为自变量,本身就已经模糊了时间的概念,对应到频域上,那么就失去了其本来的频率单位。
    同时还需要注意一点的是DTFT的结果,具有周期性,见下图公式:
    在这里插入图片描述
    从推导中可以了解, X ( e j ω 0 ) X(e^{j\omega_0}) X(ejω0)具有周期性,且周期为 2 π 2\pi 2π。这说明,DTFT之后的结果,由于是周期性的,那么只有一段 ( [ − π , π ] , o r [ 0 , 2 π ] ) ([-\pi,\pi],or [0,2\pi]) ([π,π],or[0,2π])是有用的,其他的只是这一段的周期迭代。对于实信号,其频谱是偶函数对称的,所以甚至这一段中也只有一半 ( [ 0 , π ] ) ([0,\pi]) ([0,π])是有效的。如果再将这一段进行一次归一化为 [ 0 , 1 ] [0,1] [0,1](就是全部除以 π \pi π),那么就得到了Matlab里使用pwelch函数计算功率谱结果的横轴。这也解释了下图的现象:对于同一个信号,使用Fs采样率进行采样时,得到右上的频谱,如果使用2Fs进行采样的话,就得到了右下的频谱,发现频谱被压缩了。可以这么理解:因为 π \pi π对应的真实频率增大了,所以 ω \omega ω轴上的对应真实频率的范围更广了,原本的 [ 0 , π ] [0,\pi] [0,π]对应的是真实频率的 [ 0 , F s ] [0,F_s] [0,Fs], 现在则变成了 [ 0 , π ] [0,\pi] [0,π]对应 [ 0 , 2 F s ] [0,2F_s] [0,2Fs]。所以对于 相同的 真实的 带宽 B B B,展示在新的数字频谱上,就是变窄了。
    在这里插入图片描述

    由于该变换面向的是离散信号,离散信号又意味着这是对原本连续的信号的采样、量化之后得到的结果,而采样的过程,其实就是一个信息压缩的过程。在信息论中学过,熵是描述一个消息信息量的一个指标,对于离散消息,那么其熵一定是个有限值,对于连续消息,其熵是一个无穷的值,所以从连续到离散的过程中,信息量会减少(其实说减少也不对,因为从无穷怎么也减少不到一个有限值,这可能是其中一些地方发生了本质上的变化导致的)。这也可以解释对于一个非周期连续信号,无法用有限的离散点对其采样然后原样恢复(当然对于这个问题最常规的解释还是因为采样频率不够而导致混叠)。而对于计算机来说,永远只能处理那么一段的信号,就算原本是一个周期正弦波,但是计算机处理的只有那一段,所以就相当于是一个非周期信号,在处理的过程中就会损失信息量。这也是为什么用于分析连续周期信号的傅里叶级数没有在信号处理领域得到广泛应用的原因(我乱猜的),因为数学可以很方便通过理论得到无穷的计算,在计算机中就没有无穷的概念。
  • 最后是离散傅里叶变换(DFT):从以上分析也可以发现,尽管数字采样已经对原本信息量进行了压缩,已经进行了妥协,但是仍然有一个问题就是其频域还是连续函数,而连续函数在计算机的处理过程中也是无法实现的,这就导致了DFT的出现。DFT的最简单的理解:它是在DTFT的频谱上以 N N N个等间隔点 ω k = 2 π k / N \omega_k=2\pi k/N ωk=2πk/N进行了采样。在时域采样中,采样点数需要满足一定的要求才不会导致混叠。DFT同理,但是这里需要明确的一个非常重要的点是,对DTFT的频谱进行采样之后得到的DFT序列,将其恢复到时域之后,和进行DTFT之前的时域信号不一样,是 y [ n ] = ∑ x [ n + m N ] , 0 < n < N − 1 y[n]=\sum x[n+mN], 0<n<N-1 y[n]=x[n+mN],0<n<N1这样的形式,说明是时域序列的周期延拓。那么到这就可以很容易的理解为什么N的取值很重要,因为当N较小的时候,小于序列长度时,进行周期延拓就会导致时域序列混叠,但是如果N大于了信号长度,就不会发生混叠。得到了离散点之后计算机就可以快乐的进行工作了。当然,既然又是采样,所以信息量又一次的减少,如果将DFT得到的信号进行IDFT恢复,得到的离散序列一定和原本的不一样,可能进行了周期延拓,也可能因为混叠而完全失真,这都和采样率(其实就是DFT的点数)相关。
  • 其余的类似快速傅里叶变换(FFT),没有特别需要在这里区分的必要,因为这完全是为了解决DFT计算量过大而提出的优化算法,目的、思路和DFT无异,只是实现方式很巧妙,这里不细说。

变速率信号处理的总结

意义

  1. 音频信号采样率和播放速率不同的话,会导致播放出来的音乐和原来相比快了(采样率低于播放采样率),或者慢了(采样率高于播放采样率)。
  2. 边缘陡峭的滤波器设计复杂,处理速度慢,如果抽取之后再滤波就可以加快速度提高效率(相比先处理再抽取)

分类

  1. 抽取:一般按一定整数倍率抽取,如果是3,那么隔2个抽取一个
  2. 内插:一般按一定整数倍率内插,如果是4,那么补三个0

作用效果

  1. 抽取在时域上的表现是:时域被压缩,那么意味着信号频率增大,那么频域就会展宽。所以频域上所有频率全部会乘以抽取倍率,然后进行平移相加,得到新的频谱。换句话说,就是频谱仍然以 2 π 2\pi 2π为周期重复,但是每一个 2 π 2\pi 2π周期内信号占比增大,进一步考虑可知,如果抽取倍率太高,那么频谱就会发生混叠,所以抽取其实就意味着降低采样率
  2. 内插(补零)在时域上的表现是,时域被展宽,那么意味着信号的频率全部降低,所以频谱就会收缩。所以新的频谱是原频谱所有频率除以内插倍率,得到新的频谱。这就导致一个 2 π 2\pi 2π周期内出现了很多个信号的频谱,每一个的宽度都比原来窄了很多
  • 9
    点赞
  • 31
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值