离散,离散,不如先说什么是离散时间信号。
就是这样的
离散时间信号是指在离散的时间点上取值的信号。与连续时间信号(如模拟信号)不同,离散时间信号只在特定的、等间隔的时间点上有定义。
-
离散化: 我们可以将连续时间信号想象成一条连续的曲线,而离散时间信号则是这条曲线在某些特定点上的取值。这些点之间的间隔是固定的,称为采样周期。
-
表示: 离散时间信号通常用序列的形式表示,例如
x[n]
,其中n
表示离散时间点,x[n]
表示在时刻n
的信号值。 -
采样: 将连续时间信号在等间隔的时间点上取样,得到离散时间信号。
-
直接产生: 有些信号天生就是离散的,比如数字计数器产生的信号。
-
计算机只能处理离散的数据,因此离散时间信号非常适合计算机处理。
我个人觉得最重要的区别是,这些变换把我们的信号带到了不同的变换域,这就最大的区别。当然我会详细的说明域的不同。
yes
-
时域:直观,容易理解。
-
频域:可以分析信号的频率成分,但不能反映信号的衰减或增长。
-
复频域:可以同时分析信号的频率和衰减特性,更适合系统分析。
-
Z域:适用于离散时间信号,在数字信号处理中应用广泛。
上面是域,下面就是DTFT和DFT的区别了:
-
变换域: Z变换变换到复平面,DTFT变换到连续频域,DFT变换到离散频域。
-
序列长度: Z变换适用于无限长序列,DTFT也适用于无限长序列,而DFT仅适用于有限长序列。
-
频谱: DTFT的频谱是连续的,而DFT的频谱是离散的。
如何区分:
-
看变换域: 如果变换结果是一个关于复变量z的函数,则是Z变换;如果是一个关于连续频率ω的函数,则是DTFT;如果是一个有限长的复数序列,则是DFT。
-
看序列长度: 如果处理的是无限长序列,则只能用Z变换或DTFT;如果处理的是有限长序列,则可以用Z变换、DTFT或DFT。
-
看计算方式: DFT是可以通过快速傅里叶变换(FFT)算法快速计算的。
有一个离散时间信号x[n],我们想分析它的频谱特性。
-
如果想得到关于信号所有频率成分的详细信息,可以对x[n]进行DTFT,得到一个连续的频谱。
-
如果想用计算机进行快速计算,并得到一个离散的频谱,可以对x[n]进行DFT。
-
如果想分析系统的稳定性、因果性等,可以对x[n]进行Z变换。
我随便找了本书看,差不多的
某种程度上说这就是数字处理的全部了
-
Z变换:最通用,变换域为复平面,用于系统分析。
-
DTFT:Z变换的特例,变换域为连续周期函数,用于频谱分析。
-
DFT:DTFT的离散形式,变换域为有限长离散序列,用于计算机实现。
假设有一个离散时间信号x[n],我们可以对其进行Z变换、DTFT和DFT,得到不同的变换结果。
-
Z变换:得到一个关于复变量z的函数X(z),可以分析系统的极点、零点等特性。
-
DTFT:得到一个关于频率ω的连续周期函数X(e^jω),可以观察信号的频谱分布。
-
DFT:得到一个有限长的复数序列X[k],可以利用FFT算法快速计算,并对频谱进行分析。
当然拉普拉斯修改参数值就可以得到傅里叶变换。不知道你有没有发现,其实上面还是变换的是连续的函数,就是说,有衰减因子,离散版本的拉普拉斯转换还没有着落呢!Z变换就是有衰减因子,离散版本的拉普拉斯转换。
骄傲的引用自己
说实话是有一些要区分的点的,有点头大。
在长篇大论之前,wiki的解释很棒!
不可能不放公式的
-
DTFT也是一种将离散时间信号变换到频域的工具。它将时域序列映射到一个周期连续的频谱。
-
DTFT给出了离散时间信号的连续频谱,可以更直观地观察信号的频率成分。
-
DFT是DTFT的离散形式,它将时域序列和频域序列都离散化。
-
DTFT给出了离散时间信号的连续频谱,可以更直观地观察信号的频率成分。
-
DTFT的频谱是连续的,但实际计算时只能得到有限长的频谱。
-
DTFT是Z变换的特例:当Z变换中的z取单位圆上的值时,就得到了DTFT。也就是说,DTFT是Z变换在单位圆上的投影。
-
DFT是DTFT的采样:DFT是对DTFT在频域进行等间隔采样得到的。通过DFT,我们可以得到信号的数字频谱。
DFT(Discrete Fourier Transform)是离散傅里叶变换的缩写,是一种将有限长离散时间信号从时域变换到频域的数学方法。简单来说,就是将一个序列(如声音信号的采样值)分解为不同频率的正弦波的叠加。
-
DTFT(离散时间傅里叶变换):将离散时间信号变换到连续的频域,其频谱是周期性的。
-
DFT:是DTFT的离散形式,将DTFT在频域上进行采样,得到离散的频谱。DFT的频谱也是周期性的。
DFT的计算:
X(k) = Σ(n=0 to N-1) x(n) * exp(-j*2π*k*n/N)
-
X(k) 是频域的第k个样本
-
x(n) 是时域的第n个样本
-
N 是序列的长度
-
j 是虚数单位
吴大正信号与系统-频域分析总结,可以小看一下这个里面的性质。
-
线性性: DFT是线性变换
-
周期性: DFT的频谱是周期性的
-
对称性: 实数序列的DFT具有共轭对称性
-
卷积定理: 时域卷积等于频域乘积
为了提高DFT的计算效率,人们提出了快速傅里叶变换(FFT)算法。FFT算法将DFT的计算复杂度从O(N^2)降低到了O(NlogN),大大提高了计算速度。
FFT,蝶形算法
我觉得就这样应该是解释通了,至于计算上面的细节,那我需要计算的事情。
聊聊复频域吧!
复频域是通过拉普拉斯变换或Z变换将时域信号转换得到的域。在复频域中,信号用一个复变量s(或z)的函数来表示。这个复变量s包含了实部σ和虚部jω两部分:
-
实部σ: 表示信号的衰减或增长。当σ为负时,信号随时间衰减;当σ为正时,信号随时间增长;当σ为零时,信号的幅值不变。
-
虚部jω: 表示信号的频率。ω是角频率,与普通频率f的关系为ω=2πf
连续频域没有人那么多的研究,但是它也是和上面的域有着千丝万缕的联系。
-
连续频域是复频域的一个特例: 当复频域中的s仅取虚轴上的值(即σ=0)时,就退化为连续频域。此时,信号的变换只包含频率信息,不包含衰减或增长信息。
-
复频域包含更丰富的信息: 复频域不仅包含了信号的频率信息,还包含了信号的衰减或增长信息。
假设有一个连续时间信号x(t),它的拉普拉斯变换为X(s)。
-
X(s)表示的是信号x(t)在复频域的表示。
-
当s=jω时,X(jω)就是信号x(t)的傅里叶变换,表示信号在连续频域的频谱。
把复频域想象成一个三维空间,其中横轴表示实部σ,纵轴表示虚部ω。在这个空间中,每个点都对应一个复频率。而连续频域则相当于这个三维空间中σ=0的一个平面。
这里我又要补充一些内容关于离散这个事情:一开始是写的,离散时间的周期信号的傅里叶级数,前面的离散时间没得说,我们的研究对象,周期信号和非周期信号,是要先研究周期信号的,因为把非周期信号可以两边延拓为周期信号,至于傅里叶级数,那只是我们从连续到离散的一个梯子。
离散时间周期信号的傅里叶级数
(连续)离散时间,周期信号的傅里叶级数表示.完全推导版 , 这篇写的很清楚了,但是我还是觉得有必要重新摘取一些有用的。
事实上,最后给出的结果就是这样,最引人注目的地方就是它的周期性。
我们一般是对偶的看这两个
在周期离散信号的基础上我们才开始走向离散时间傅里叶变换(DTFT)
上面是级数,下面是变换:
上面是一个非周期的信号,但是延拓出了周期信号,两边copy
从周期走向非周期,中间的变化是重点
上面也放过,但是这里也要有
OK,我记得推倒的时候较为突兀的定义过一个函数,后来想通了,过了几天又忘了。
在这个公式里面,第一个是周期信号,然后是对一个周期求和。接下来是在这个区间之外都是周期信号等于一个常信号。
第三个成立是因为这样的
文章中说定义这个函数
其实在形式上面看就是第三个公式的样子,不过是,kω变换了,其实就是k没有了,因为上面是系数中表征变几次谐波的参数
三式
下面这个解释就对了
反正就是我们通过变换其实是得到了一个非周期的片段的系数是什么样
然后把这个系数带入:离散时间周期信号的傅里叶级数对中的综合公式
接着看这个,其实还是在对上面的公式做变换,就是系数在非周期的情况下,对角频率公式做变换。带入公式,1/N换成了ω0/2π。
就得到了有些怪异的这个,这个公式是出现了两个e函数
是因为经过变换以后,就是系数上面的变换
就是说这上下两个公式是一样的,这个括号里面的(e)其实是Ak的系数给的。
所以在代入之后,出现了这样的结构
接着处理这个有限项求和
上面公式的可视化就是这个面积
也就是说这个离散非周期时间傅里叶变换,天然的带着周期
一个周期信号 𝑥~[𝑛] 的傅里叶系数 𝑎𝑘 可以用一个有限长序列 𝑥[𝑛] 的傅里叶变换 𝑋(𝑒𝑗𝜔) 的等间隔样本来表示,这个 𝑥[𝑛] 就是 𝑥~[𝑛] 的一个周期,而在其余地方为零。这一点在实际的信号处理和傅里叶分析中极为重要。
离散时间傅里叶变换和连续时间情况相比具有许多类似之处。
两者的主要差别在于离散时间变换 𝑋(𝑒^𝑗𝜔) 的周期性和在综合公式中的有限积分区间。这两者均来自这样一个事实:在频率上相差 2𝜋 的离散时间复指数信号 𝑒^𝑗𝜔𝑛 是完全一样的。
这个是离散时间周期信号
对周期离散时间信号而言,这就意味着傅里叶级数系数也是周期的,并且傅里叶级数表示式是一个有限项的和式。
对非周期信号而言,这就意味着 𝑋(𝑒^𝑗𝜔) 也是周期的(周期为 2𝜋 ),并且综合公式只涉及在一个频率区间内的积分,这个频率区间就是产生不同复指数信号的那个间隔,即任何 2𝜋 长度的间隔。
前面指出离散时间指数信号与正弦信号,𝑒^𝑗𝜔𝑛 作为 𝜔 函数的周期性导致的结果是:𝜔=0 和 𝜔=2𝜋 都得出同一个信号。
因此,位于这些频率值或任何 𝜋 偶数倍的 𝜔 附近都是慢变化的,从而都属于低频率的信号;而靠近 𝜋 的奇数倍的 𝜔 ,在离散时间情况下都属于高的频率。
奇数是高频率
如果单看这个图像,不看轴
这个又是一个连续周期性的图形,似乎应该可以用连续时间傅里叶级数的分析式。和方波差不多。
上面是非周期离散时间的傅里叶级数变换,下面是连续时间傅里叶级数对重写。
可以看到都是中间积分式, 且 𝑥[𝑛] 和 𝑎𝑘 都是离散的序列,中间都是无穷项求和,求和的结果 𝑋(𝑒𝑗𝜔) 和 𝑥(𝑡) 都是连续的周期函数。这真是有趣的现象!
连续时间周期信号的傅里叶级数对,如果换种理解角度,就变成了离散时间傅里叶变换了。
下面是一个例题
a展示出了 𝑎>0 时,𝑋(𝑒𝑗𝜔) 的模和相位;b展示出 𝑎<0 时的模和相位。
应该注意,图中所有这些函数都是周期为 2𝜋 的周期函数。
上述导出离散时间傅里叶变换过程中,将非周期信号 𝑥[𝑛] 看成周期序列 𝑥~[𝑛] 的一个周期,这意味着非周期信号 𝑥[𝑛] 一定是有限长的?
事实上,离散时间傅里叶变换对也适用于无限长的非周期序列 𝑥[𝑛] ,但由于分析式是无限项求和,此时需要考虑收敛性问题。
-
离散时间傅里叶变换导出的基本思路:周期序列傅里叶级数在周期 𝑁→∞ 时的极限情况; -
离散时间傅里叶变换综合式的物理意义:基本序列 𝑒^𝑗𝜔𝑛 的适当线性组合(积分)可以构成非周期序列 𝑥[𝑛] ; -
离散时间傅里叶变换分析式的物理意义:𝑋(𝑒^𝑗𝜔) 提供了非周期序列 𝑥[𝑛] 由基本序列 𝑒^𝑗𝜔𝑛 线性组合(积分)时的“成分表”; -
所有的傅里叶级数/变换式中,如果是无穷项求和,无穷区间积分,则应考虑收敛性问题; -
离散时间傅里叶变换 𝑋(𝑒^𝑗𝜔) 具有 2𝜋 的周期性,但连续时间傅里叶变换 𝑋(𝑗𝜔) 通常是非周期的,其根本原因是基本序列 𝑒^𝑗𝜔𝑛 和基本信号 𝑒^𝑗𝜔𝑡关于 𝜔 的周期性不同; -
离散时间信号的高、低频:𝑋(𝑒^𝑗𝜔) 具有 2𝜋 的周期性,且 𝜔 取 𝜋 的偶数倍附近时,为低频成分, 𝜔 取 𝜋的奇数倍附近时,为高频成分。
https://zhuanlan.zhihu.com/p/640055182
奥本海姆的书确实是好,但是⛓️太长了。
吴大正书籍,其实他的书里面不严格区分周期不周期的事情(区分的)
这章称为序列分析,一开始也是使用周期序列的离散傅里叶级数入手的
给了一道例题
DTFT的推倒很紧凑,值得好好看,但是推倒思想是一样的,符号是不统一的
书里面给出了一个重要的看法,周期和域之间是对偶的
下面是一个DSP的东西,不写了。
https://arm-software.github.io/CMSIS-DSP/latest/
https://arm-software.github.io/CMSIS_6/latest/DSP/index.html#autotoc_md1
https://www.keil.arm.com/packs/?page=2
https://www.arm.com/technologies/cmsis