族谱网 头条 人物百科

离散傅里叶变换矩阵

2020-10-16
出处:族谱网
作者:阿族小谱
浏览:1199
转发:1
评论:0
定义N点的离散傅里叶变换可以用一个n××-->m{displaystylentimesm}的矩阵乘法来表示,即X=Wx{displaystyleX=Wx},其中x{

定义

N点的离散傅里叶变换可以用一个n× × -->m{\displaystyle n\times m}的矩阵乘法来表示,即X=Wx{\displaystyle X=Wx},其中x{\displaystyle x}是原始的输入信号,X{\displaystyle X}是经过离散傅里叶变换得到的输出信号。 一个n× × -->n{\displaystyle n\times n}的变换矩阵W{\displaystyle W}可以定义成W=(ω ω -->ij)i,j=0,… … -->,N− − -->1/N{\displaystyle W=(\omega ^{ij})_{i,j=0,\ldots ,N-1}/{\sqrt {N}}},或等效如下:

其中ω ω -->{\displaystyle \omega }是1{\displaystyle 1}的n{\displaystyle n}次方根的主值(primitive nth root of unity),大小为e− − -->2π π -->iN{\displaystyle e^{\frac {-2\pi i}{N}}}。需要注意的是在总和前面的正规化因数1N{\displaystyle {\frac {1}{\sqrt {N}}}},还有ω ω -->{\displaystyle \omega }中指数的正负号是依据惯例,并且会因为处理的方法有所不同。以下所有的讨论考虑到大多数的细节变动且不论是否为一般惯例均适用之。唯一重要的是,正变换和逆变换有相反的指数正负号标志,而其正规化因数乘积为1N{\displaystyle {\frac {1}{N}}}。然而,这里为了使得最后的离散傅里叶变换矩阵结果正规化所选择的因数 ,在许多情况下都是通用的。

快速傅里叶变换算法利用矩阵的对称性与W的周期性,以减少乘法所需要的时间(把计算复杂度从O(N2){\displaystyle O(N^{2})}降为O(Nlog⁡ ⁡ -->N){\displaystyle O(N\log N)})。类似的方法也可适用于其他矩阵乘法如阿达马矩阵和Walsh matrix(英语:Walsh matrix)。

特殊情况

3点的离散傅里叶变换具有特殊的意义。例如:Charles Legeyt Fortescue于1918 所发表的对称分量变换(Symmetrical Components Transform, SCT),它定义了三相平衡(three phase balance),即3点离散傅里叶变换可分解成一个直流成分,以及两个交流成分(一个是顺时针相位,另一个为逆时针相位)。

例子

两点离散傅里叶变换矩阵

两点的离散傅里叶变换是一个很简单的例子,其第一列代表是直流成分(总和)和第二列是交流成分(差异)。

第一列处理总和的部分,第二列处理相差的部分。 因数1/2{\displaystyle 1/{\sqrt {2}}}致使整个矩阵规一化(见下文)。

四点离散傅里叶变换矩阵

四点的离散傅里叶变换矩阵如下:

八点离散傅里叶变换矩阵

八点的离散傅里叶变换矩阵如下:

其中

以下用图片来解说离散傅里叶变换的矩阵乘法概念:

图中实部(余弦波)是由实线代表,虚部(正弦波)由虚线代表。 最上面一行全为1,(透过乘上1/8{\displaystyle 1/{\sqrt {8}}}来规一化),因此这个部分代表输入信号的直流分量。下一行是8个负一次循环的复指数取样(samples of negative one cycle of complex exponential),即分频(fractional frequency)为−1/8倍频率的信号。因此,这一行代表在分频+1/8的信号强度。再下一行是8个负二次循环的复指数取样,所以它代表-1/4倍的分频。因此,这一行代表在分频+1/4的信号强度。 以下总结了八点离散傅里叶变换代表的意义,依行排序,以分频表示:

0代表直流信号成分

-1/8代表分频为+1/8 的信号强度

-1/4代表分频为+1/4 的信号强度

-3/8代表分频为+3/8 的信号强度

-1/2代表分频为+1/2 的信号强度

-5/8代表分频为+5/8 的信号强度

-3/4代表分频为+3/4 的信号强度

-7/8代表分频为+7/8 的信号强度

等效上最后一行,可以当作是分频为+1/8即代表分频-1/8的信号强度。如此一来,则可以说这个矩阵的上面列是信号的正频率部分的强度而下面列是信号负频率部分的强度。

规一化变换(unitary transform)

离散傅里叶变换(或可能是透过适当的尺度选择)是一个规一化的变换,即符合能量保留(preserves energy)。可以达到规一化的合适尺度是1/N{\displaystyle 1/{\sqrt {N}}},这使得能量物理意义上跟在傅里叶定义上是一样的,即满足Parseval定理(Parseval"s theorem)。(其他未规一化的尺度,也普遍被使用以方便计算;例如,折积定理(convolution theorem)需较简单的形式与尺度选择,详述于离散傅里叶变换条目中) 。

其他性质

其他离散傅里叶变换矩阵的性质,包括其特征值(特征向量),与折积的关系,应用等,请参见离散傅里叶变换条目。

限制:傅里叶运算(Fourier operator)

如果我们作出一个非常大的矩阵,其中列元素为复指数(即,余弦实部和正弦虚部),并增加分辨率而不考虑边界,我们可近似第二型Fredholm积分方程(由傅里叶运算定义连续傅里叶变换)的”核”(kernal)。此连续傅里叶变换的一部分类似于离散傅里叶变换矩阵,如图所示。其中灰阶像素值的数值是指数量。

参考

The Transform and Data Compression Handbook by P. C. Yip, K. Ramamohan Rao- See chapter 2 for a treatment of the DFT based largely on the DFT matrix


免责声明:以上内容版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。感谢每一位辛勤著写的作者,感谢每一位的分享。

——— 没有了 ———
编辑:阿族小谱
发表评论
写好了,提交
{{item.label}}
{{commentTotal}}条评论
{{item.userName}}
发布时间:{{item.time}}
{{item.content}}
回复
举报
点击加载更多
打赏作者
“感谢您的打赏,我会更努力的创作”
— 请选择您要打赏的金额 —
{{item.label}}
{{item.label}}
打赏成功!
“感谢您的打赏,我会更努力的创作”
返回

更多文章

更多精彩文章
打赏
私信

推荐阅读

· 离散傅里叶变换
定义对于N点序列{x[n]}0≤≤-->n<N{\displaystyle\left\{x[n]\right\}_{0\leqn傅里叶的离散傅里叶变换(DFT)为其中e{\displaystylee}是自然对数的底数,i{\displaystylei}是虚数单位。通常以符号F{\displaystyle{\mathcal{F}}}表示这一变换,即离散傅里叶变换的逆变换(IDFT)为:可以记为:实际上,DFT和IDFT变换式中和式前面的归一化系数并不重要。在上面的定义中,DFT和IDFT前的系数分别为1和1/N。有时会将这两个系数都改成1/N{\displaystyle1/{\sqrt{N}}}。从连续到离散连续时间信号x(t)以及对应的连续傅里叶变换x^^-->(ωω-->){\displaystyle{\hat{x}}(\omega连续函数都是连续函数。由于数字系...
· 离散时间傅里叶变换
定义一组离散的实数或复数:x[n](n为所有整数)的离散时间傅里叶变换是产生以频率为变量的周期函数的一个傅里叶级数。当频率变量ω的单位是归一化的弧度/样本时,周期为2π,而傅里叶级数为:此频率域函数的性质源于泊松求和公式(英语:Poissonsummationformula)。令X(f)为任意函数x(t)的傅里叶变换,采样间隔为T(秒),等价于序列x[n](或与之成正比),即T⋅⋅-->x(nT)=x[n]{\displaystyleT\cdotx(nT)=x[n]}。则以傅里叶级数表示的周期函数是X(f)的周期求和。赫兹以赫兹(周期/秒)为单位的频率f{\displaystyle\textstylef}的话就会是:图一.傅立叶变换(左上)和左下的其周期求和(DTFT)的图示。右下角显示了用离散傅里叶变换(DFT)计算DTFT的采样。整数k的单位为转/样本,采样频率是1/T,fs(样本/秒...
· 离散盘
形成阋神星,已知最大的黄道离散天体,和它的卫星阋卫一(迪丝诺美亚,中央偏左的小光点)。我们对离散盘的所知非常有限,虽然天文学的主流观点认为它是太阳系形成的早期过程中,因为海王星向外迁徙造成的引力扰动才被从柯伊伯带散射入高倾斜和高离心率的轨道内。相比之下,柯伊伯带像是一个相对“圆”和“平坦”的甜甜圈,以平和的圆轨道(QB1天体)和略为椭圆的共振轨道(类冥天体),将天体约束在30至44天文单位的圆环内;离散盘内的黄道离散天体轨道环境就比较怪异了。黄道离散天体,就以矮行星阋神星为例,在垂直黄道方向上的距离几乎和平行方向上与太阳的距离一样远;轨道模拟也显示黄道离散天体的轨道是怪异且不稳定的,并且最终会从太阳的核心区域抛至奥尔特云甚至更遥远的地方。有些迹象显示半人马群只是单纯的从柯伊伯带被向内抛射,而不是向外抛射的天体,可以称为“内海王星天体”(cis-Neptunianobject)。事实上,有些...
· 离散信号
参见采样采样定理信号连续信号离散时间信号
· 离散群
性质因为拓扑群是齐次的,你只需要查看一个单一的点就能确定这个群是否为离散的。特别是,拓扑群是离散的,当且仅当包含单位元的单元素集合是开集。离散群是和零维李群同样的东西(不可数离散群不是第二可数的,所以要求李群满足这个公理的作者不把这些群认做李群)。离散群的单位元单元就是平凡子群而单元的群同构于这个群自身。因为只有在有限集合上的豪斯多夫拓扑是离散拓扑,有限豪斯多夫拓扑群必然是离散群。可得出所有的豪斯多夫群的有限子群是离散群。G的离散子群H是馀紧致(cocompact)的,如果有G的紧子集K使得HK=G。离散正规子群在覆盖群和局部同构群的理论中扮演重要角色。连通群G的离散正规子群必然位于G的中心并因此是阿贝尔群。其他性质:所有离散群的子群都是离散群。所有离散群的商群都是离散群。有限个离散群的乘积是离散群。离散群是紧群当且仅当它是有限的。所有离散群都是局部紧群。所有豪斯多夫群的离散子群都是闭合的...

关于我们

关注族谱网 微信公众号,每日及时查看相关推荐,订阅互动等。

APP下载

下载族谱APP 微信公众号,每日及时查看
扫一扫添加客服微信