信号处理中傅里叶变换简介由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“数字信号处理傅里叶”。
傅里叶变换
一、傅里叶变换的表述
在数学上,对任意函数f(x),可按某一点进行展开,常见的有泰勒展开和傅里叶展开。泰勒展开为各阶次幂函数的线性组合形式,本质上自变量未改变,仍为x,而傅里叶展开则为三角函数的线性组合形式,同时将自变量由x变成ω,且由于三角函数处理比较简单,具有良好的性质,故被广泛地应用在信号分析与处理中,可将时域分析变换到频域进行分析。
信号分析与处理中常见的有CFS(连续时间傅里叶级数)、CFT(连续时间傅里叶变换)、DTFT(离散时间傅里叶变换)、DFS(离散傅里叶级数)、DFT(离散傅里叶变换)。通过对连续非周期信号xc(t)在时域和频域进行各种处理变换,可推导出以上几种变换,同时可得出这些变换之间的关系。以下将对上述变换进行简述,同时分析它们之间的关系。
1、CFS(连续时间傅里叶级数)
在数学中,周期函数f(x)可展开为
由此类比,已知连续周期信号x(t),周期为T0,则其傅里叶级数为
其中,为了简写,有
其中,为了与复数形式联系,先由欧拉公式ejz=cosz+jsinz得
故有
令
则
对于Dn,有
n≤0时同理。故
CFS图示如下:
Figure 1
理论上,CFS对于周期性信号x(t)在任意处展开都可以做到无误差,只要保证n从-∞取到+∞就可以。在实践中,只要n取值范围足够大,就可以保证在某一点附近对x(t)展开都有很高的精度。
2、CFT(连续时间傅里叶变换)
连续非周期信号x(t),可以将其看成一连续周期信号期T0→∞。当然,从时域上将x(t)进行CFS展开,有的周也可以反过来看成x(t)的周期延拓。
若令
则
有
T0→∞使得Ω0→0,则
由此,定义傅里叶变换与其逆变换如下 CFT:
CFT-1:
x(t)是信号的时域表现形式,X(jΩ)是信号的频域表现形式,二者本质上是统一的,相互间可以转换。CFT即将x(t)分解,并按频率顺序展开,使其成为频率的函数。上式中,时域自变量t的单位为秒(s),频域自变量Ω的单位为弧度/秒(rad/s)。
CFS中的Dn与CFT中的X(jΩ)之间有如下关系
即从频域上分析,Dn是对X(jΩ)的采样(可将Figure 1与Figure 2进行对比)。
CFT图示如下:
Figure 23、DTFT(离散时间傅里叶变换)
首先,先从连续信号得到离散信号。用冲激信号序列
对连续非周期信号xc(t)进行采样,采样间隔为Ts,有
此时的xs(t)还不是真正的离散信号,它只是在满足t = nTs的时间点上有值,在其它时间点上值为零。对xs(t)进行进一步处理有
规定
则
其中,x[n]是最终所得的离散信号。xs(t)自变量为t,其单位为秒s,间隔为TS;x[n]自变量为n,其单位为1,间隔为1。
从频域分析上有
其中
。令,定义
以上式为DTFT定义式。DTFT逆变换为
DTFT是在时域上对CFT的采样(图示可见Figure 3与Figure 4),在DTFT中,时域信号x[n]为离散的,而对应的频域表示X(ejω)为连续的,且有周期ωs = 2π。
X(ejω)与Xs(jΩ)之间的关系为
ω = ΩTs
Xs(jΩ)中,自变量Ω单位为弧度/秒(rad/s),周期为Ωs = 2π/Ts;X(ejω)中,自变量ω单位为弧度(rad),周期为ωs = 2π。
CFT时域采样图示如下:
Figure 3
DTFT图示如下:
Figure 44、DFS(离散时间傅里叶级数)
在离散时间信号x[n]基础上,用冲激序列
对DTFT中的X(ejω)进行采样,采样间隔为Δω = 2π/N,则有
而S(ω)的逆DTFT变换为
对Xs(ejω)进行逆DTFT变换,有
xs[n]相当于对x[n]进行了周期延拓,周期为N = 2π/Δω。由上式可得
若延拓周期N大于x[n]的时长,则延拓不会发生混叠,于是
k为任意整数
令周期信号,k为任意整数,则
有
取ω = 2πk/N,令
则有
是以k为自变量的函数,有以下性质
m为任意整数
即的周期为N。为了避免重复计算,我们只考虑一个周期N内的情况,即
同时,为时域表示,为频域表示。故定义DFS为
其逆变换为的自变量n单位为1,周期为N;的自变量k单位为1,周期也为N。DFS应用于离散时间周期性信号中,其相当于在频域中
对DTFT采样,而对应地在时域中相当于对DTFT进行周期延拓(图示见Figure 5与Figure 6)。DFS与DTFT的关系为
DTFT频域采样图示如下:
Figure 5
DFS图示如下:
Figure 65、DFT(离散傅里叶变换)
在DFS基础上,取离散时间周期性信号0,1,2,…N-1这一个周期内的N个点,得
其中,RN[n]表示当n = 0,1,2,…N-1时函数取值为1,当n取其它值时函数取值为0。定义DFT为的基础上,其逆变换为
xd[n]的自变量n单位为1,时长为N;Xd[k]的自变量k单位为1,时长也为N。DFT相当于对DFS的时域及频域都取0,1,2,…N-1这一个周期内的N个点。
6、傅里叶变换之间的关系
傅里叶变换之间的关系主要有两点,一是采样与周期延拓之间的对应关系,二是对自变量的替换关系。(1)采样与周期延拓之间的对应关系
采样与周期延拓之间是一种对应关系,时域中对信号采样相当于在频域中对信号进行周期延拓,同样地,频域中对信号采样相当于在时域中对信号进行周期延拓,二者间是对应与平行的关系,不存在因果关系。
傅里叶变换中的CFS、CFT、DTFT、DFS、DFT可由连续非周期信号xc(t)进行采样及周期延拓处理得到各种变换,它们之间的关系如图Figure 7与Figure 8:
Figure 7
Figure 8
上两图中,蓝色箭头表示在时域或频域中采取的主动措施,白色箭头表示在频域或时域中产生的相应变换。(2)对自变量的替换关系
在对信号进行采样与周期延拓的同时,对自变量进行某种替换,从而完成傅里叶变换类型的转变。
傅里叶变换中对自变量的替换情况如图Figure 9所示。CFS适用于连续周期性信号,其自变量t单位为秒(s),相应的幅频谱|Dn|中,自变量n单位为1。而CFT适用于连续非周期信号xc(t),其自变量t单位为秒(s),对应的频域信号为Xc(jΩ),其自变量Ω单位为弧度/秒(rad/s)。由CFS变成CFT相当于连续周期性信号的周期T0趋于无穷,而在频域中则为自变量的替换,由n变成Ω,替换关系为
DTFT适用于离散时间信号x[n],其自变量n单位为1,对应的频域信号为X(ejω),自变量ω单位为弧度(rad)。由CFT变成DTFT相当于对连续信号xc(t)采样及离散化,自变量由t替换为n,替换关系为t = nTs,而在频域中则为周期延拓及自变量的替换,由Ω替换为ω,替换关系为ω = ΩTs。
DFS适用于离散周期性信号频域信号为,其自变量n单位为1,对应的,自变量k单位为1。由DTFT变成DFS相当于在频
域中对X(ejω)进行采样、离散化与自变量替换,由ω替换为k,替换关系为ω = 2πk/N。
DFT的时域与频域序列长度都为N个点(0,1,2,…N-1),时域自变量n单位为1,频域自变量k单位为1。
由图Figure7、Figure 8和Figure 9可以清楚地研究非相邻变换之间的关系。
Figure 9
二、与相关教材内容的辨析
1、《Signal Proceing and Linear Systems》(B.P.Lathi, Oxford University Pre)
书中首先将高等数学中的向量理论扩展到了信号系统中,引出正交信号空间的定义,指出任意信号x(t)可用正交信号空间的线性组合表示,进而引出三角傅里叶级数,将这种表示用三角函数的线性组合表示。CFS的来源介绍比我对CFS的自述更加详细具体,更有逻辑性,体现了高等数学的延伸,CFS定义部分与我的自述大体相同。
书中由CFS引出CFT,指出连续非周期信号xc(t)相当于将连续
周期性信号的周期T0趋于无穷,然后对xc(t)按照CFS方法展开,中间过程中引出了CFT。这一部分与我的自述大体相同。只是我在对傅里叶变换的总结中将xc(t)进行无混叠的周期性延拓,反向也得出了。这只是对傅里叶变换的又一种理解,但从本源上考虑,还应该是由连续周期性信号
得出连续非周期信号xc(t)。
书中接下来先介绍的是DFS。书中由CFS类比定义了DFS,定义为
其中,这种定义与我对DFS的自述略有差别。书中完全按照CFS的定义模式定义的,书上在此之后也按照CFS的模式给出了Dr的幅频谱与相频谱。而我的自述则采用类似CFT的定义方式,即正变换为从时域变到频域,逆变换为从频域变到时域,其次书中使用的字母表示方式与我的自述略有差异,不过本质上意义是相同的。
紧接着,书中由DFS引出了DFT,指出DFT的时域及频域都为N点有限序列,此处与我对DFT的自述大体相同,但未进行深入说
明。之后,类似于由CFS引出CFT,书中由DFS中的离散时间周期函数引出离散时间非周期函数x[k](令周期N0→∞),然后对x[k]按照DFS的方法展开,在中间推导过程中引出了DTFT。总之,在离散时间信号的傅里叶变换中,书上是类比CFS引出CFT的模式,由DFS引出DTFT,而DFT也由DFS引出,只是未做重点讲解,实质上是从时域角度出发,与连续时间信号进行同等过程的类比。我对离散时间信号傅里叶变换的自述则从频域角度出发,与连续时间信号的时域推导过程进行同等过程的类比。二者分析方向不同,顺序不同,但本质上是相同的。这也从侧面反映出傅里叶变换将单纯的时域分析引向时域与频域的双领域分析,增加了对信号分析与处理的方法与方向,有利于更好地对信号进行理解。
2、《信号与系统》
书中也是首先将高等数学中的向量理论扩展到了信号系统中,引出正交信号空间的定义,指出任意周期为T0的信号x(t)可进行正交分解,而正余弦信号集是比较特殊的正交信号集,并用正余弦信号集表示信号,达到一种分解的目的,从而定义出CFS,并将正余弦信号集进一步扩展为虚指数信号集,从而将指数形式的CFS表示出来。在表示方式上与我的自述基本相同。而书中对三角形式的CFS与指数形式的CFS总结比较清楚,并对各自形式的幅频谱进行了比较,指出指数形式CFS的频谱为双边谱,而三角形式的CFS的频谱为单边谱。而由CFS导出CFT的叙述则基本与我的自述相同,即连续非周期信号xc(t)相当于将连续周期性信号的周期T0趋于无穷,然后对xc(t)按照CFS方法展开,中间过程中引出了CFT。
书中对DFS的描述,类比于对CFS的描述,采用离散形式的虚指数正交信号集对离散时间周期性信号表示,表示方式与上一本书相同。由DFS引出DTFT时类比于由CFS引出CFT的过程,将离散时间周期性信号周期趋于无穷,得出离散时间非周期性信号,按照DFS的方式对信号进行分解表示,在推导过程中引出DTFT的定义,过程与上一本书基本相同。而DTFT也可对离散时间周期性信号进行处理。对DFT并未做重点描述。
总之,两本书对傅里叶变换的描述都是先对连续时间信号进行讨论,然后离散时间信号中的讨论参考连续时间信号中的讨论,层次清晰,可比性强。我的自述主要侧重于对信号的时域或频域进行各种处理,引出傅里叶变换的各种形式,可加深对傅里叶变换各种形式之间关系的理解。
三、傅里叶变换的应用
1、应用
傅里叶变换主要是为了将一般性的信号用较规则的、性质良好的三角函数进行表示,从而可以从频域的角度进行信号分析与处理,扩充了信号分析与处理的分析领域,简化了分析与处理的过程。从理论上,CFS、CFT、DTFT、DFS、DFT在满足相应的条件下都可以使用。而在实际应用中,计算机只能处理离散的、序列长度有限的信号,故实际应用中,DFT具有应用价值,其它形式的傅里叶变换处理的信号
是连续的或无限长的,计算机无法处理,所以只能在理论上进行数学运算。而DFT利用计算机可以快速算出,被称为快速傅立叶变换(FFT)。FFT可以减少计算DFT时乘法的使用次数,简化运算,提高效率。而现代信号分析与处理中必然要对信号进行采样离散化,输入到计算机中进行处理,得到频域形式,所以DFT的实际应用是很广泛的。
2、限制条件及潜在问题
CFS只适用于连续周期性信号,CFT只适用于连续非周期信号,DTFT只适用于离散时间信号,DFS只适用于离散时间周期信号,DFT只适用于有限序列的离散时间信号。CFS、CFT、DTFT、DFS处理的信号具有连续性或无限长特性,适用于在理论上的定性分析,而在实际应用中,我们需要快速高效地处理信号,这必然用到计算机,而计算机只能处理离散的、有限序列长度的信号,故只有DFT有实用意义,CFS、CFT、DTFT、DFS则不行。而DFT计算需要大量的加法与乘法,往往实际应用中不能直接应用,所以实际应用中要根据需要进行优化处理,在提高运算速度与精度之间进行权衡,原始的DFT只是具有实际应用中的象征意义。
(一).信号分析1、编制信号生成程序,产生下述各序列,绘出它们的时域波形1) 单位抽样序列 (n)2)矩形序列 RN(n)3)三角波序列n1,0n3x3(n)8n,4n70,其它4) 反三角波序列4n,0n3x4(n)n3,4n7......
《TMS320C54XDSP语音信号处理》DSP课程设计报告学号:0906030430班级:通信四班姓名:唐溧同组姓名:蒋双福,梁嘉杰指导教师:徐淑芳二〇一二 年四 月十四 日目录《TMS320C54XDSP语音信......
目 录第一章 绪论 .........................................................................................................11.1课程设计的目的及意义 ...................
京信通信系统(广州)有限公司广东分公司直放站及室内分布系统信号覆盖故障现象、产生原因及处理方法目录一、无信号二、覆盖区信号质差三、上行干扰四、掉话五、有信号却不能打......
摘 要语音信号处理是研究数字信号处理技术和语音信号进行处理的一门学科,是一门新型的学科,是在多门学科基础上发展起来的综合性技术,它涉及到数字信号处理、模式识别、语言学......