傅里叶变换、拉普拉斯变换、Z 变换的联系是什么?为什么要进行这些变换?
原创 逸珺 嵌入式客栈
[导读] 在知乎上看到一个问题,傅里叶变换、拉普拉斯变换、Z 变换的联系是什么?为什么要进行这些变换?我觉得这是一个非常好的问题,貌似一下子也回答不上来,所以整理学习并分享一下。
什么是数学变换?
要理解这些变换,首先需要理解什么是数学变换!如果不理解什么是数学变换的概念,那么其他的概念我觉得也没有理解。
数学变换是指数学函数从原向量空间在自身函数空间变换,或映射到另一个函数空间,或对于集合X到其自身(比如线性变换)或从X到另一个集合Y的可逆变换函数。比如(图片来源wikipedia):
数学中还有很多其他的数学变换,其本质都可以看成是将函数f(x)利用变换因子进行的一种数学映射,其变换结果是函数的自变量有可能还是原来的几何向量空间,或许会变成其他的几何向量空间,比如傅立叶变换就从时域变换为频域。
而傅立叶变换和拉普拉斯变换的本质都是对连续或有限个第一类间断点函数的一种积分变换,那么什么是积分变换呢?
什么是积分变换?
积分变换通过对原函数对映射函数空间自变量在特定区间进行积分运算,将函数从其原始函数空间映射到另一个函数空间。这样一来,其中原始函数的某些属性在映射函数空间可能比原始函数空间更容易表征或分析。通常可以使用逆变换将变换后的函数映射回到原函数空间,这样的变换称为可逆变换。
假定对于函数为自变量t的函数f(t),通常积分变换都具有如下类似的范式:
(Tf)(u)=∫t1t2f(t)K(t,u)dt(Tf)(u)=\int_{t1}^{t2}{f(t)K(t,u)dt} (Tf)(u)=∫t1t2f(t)K(t,u)dt
函数f(t)是该变换的输入,(Tf)(u)为变换的输出,因此积分变换一般也称为一种特定的数学运算符。而函数K(t,u)称为积分核函数(kernel function)。
这里有一个对称核函数的概念,这是什么意思呢?就是将函数K的两个自变量交换位置仍然相等:
K(t,u)=K(u,t)K(t,u)=K(u,t) K(t,u)=K(u,t)
有的变换可逆,这是什么概念呢?就是变换后通过逆变换,还能还原!
f(t)=∫u1u2(Tf)(u)K−1(u,t)duf(t)=\int_{u_1}^{u_2}{(Tf)(u)K^{-1}(u,t)du} f(t)=∫u1u2(Tf)(u)K−1(u,t)du
观察正变换与逆变换,你会发现:
- 核函数刚好两个自变量交换位置
- 正变换是对原函数f(t)在时间维度上进行积分
- 逆变换是在变换后的函数在u维度上进行积分
什么是傅立叶级数?
在谈傅立叶变换之前,先谈谈傅立叶级数会更容易理解傅立叶变换。在数学中,傅里叶级数(Fourier series)是把类似波的函数表示成简单正弦波的方式。更正式的说法是,它能将任何周期性函数或周期性信号分解成一个(可能由无穷个频率分量组成的)简单振荡函数的集合,即正弦函数和余弦函数(或者,等价地使用复指数),从数学的定义来看:
设f(t)是一周期信号,假定其周期为T。若f(t)在一个周期的能量是有限的,就是:
∫−T2T2∣f(t)∣dt<∞\int_{-\frac{T}{2}}^{\frac{T}{2}}|f(t)|dt\lt\infty ∫−2T2T∣f(t)∣dt<∞
则,可以将f(t)展开为傅立叶级数。怎么展开呢?计算如下:
f(t)=∑k=−∞∞F(kω)ejkωtf(t)=\sum_{k=-\infty}^{\infty}{F(k\omega)e^{jk\omega t}} f(t)=k=−∞∑∞F(kω)ejkωt
而傅立叶级数的系数由下式计算:
F(kω)=1T∫−T2T2f(t)e−jkωtdtF(k\omega)=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-jk\omega t}dt F(kω)=T1∫−2T2Tf(t)e−jkωtdt
对于f(t),利用欧拉公式还可以写成正弦函数与余弦函数的和,这里就不写了。欧拉公式如下:
ejt=cos(t)+jsin(t)e^{jt}=cos(t)+jsin(t) ejt=cos(t)+jsin(t)
公式中的k表示第k次谐波,这是个什么概念呢?不容易理解,看下对于一个方波的前4次谐波合成动图就比较好理解了。这里合成的概念是指时域上的叠加的概念,图片来源wikipedia
从上图可以直观看出,周期性方波,可以看成多次谐波的线性叠加,其幅度谱图,是一根根离散的谱线,且幅度值越来越低,从这个角度可以看出高次谐波的分量,占比越来越小。其谱线的位置为:
- 第一根为:ω1=2πT\omega_1=\frac{2\pi}{T}ω1=T2π
- 第二根为:ω2=2×2πT\omega_2=2 \times \frac{2\pi}{T}ω2=2×T2π
- 第n根为:ωn=n×2πT\omega_n=n \times \frac{2\pi}{T}ωn=n×T2π
其谱线的间隔为:ω=2πT\omega=\frac{2\pi}{T}ω=T2π
应用:这里可以联想到我们的电子系统中的时钟信号,做硬件的朋友或有经验,在做EMC的辐射测试时,发现产品电路板在某些频点超标,有经验的同学会很快定位到辐射源。其实这里大概率就是因为周期性的时钟信号造成的,从频率的角度可以看成是其基频的多次谐波的线性叠加,而某个谐波分量在电路线路尺寸满足辐射条件时,就从电路板上脱逸而出,变为电磁波能量向空间传播。所以反向去查该频率可能对应的周期性时钟信号的基频就能很快定位到辐射源,从而解决问题。
说到傅立叶级数是周期性信号可以用傅立叶级数展开,那么是不是任一周期性信号都可以进行傅立叶级数展开呢?答案是否定的,必须满足著名的狄利克雷(Dirichlet)条件:
- 在一周期内,如果有间断点存在,则间断点的数目需要是有限个数
- 在一周期内,极大值和极小值的数目是有限个数的
- 在一周期内,信号或者函数是绝对可积分的。见前文公式。
什么是傅立叶变换?
前面说了傅立叶级数,接下来再看傅立叶变换。傅立叶变换之所以称为傅立叶变换,是由于1822年,法国数学家傅立叶(J.Fourier) 在研究热传导理论时首次证明了将周期函数展开为傅立叶级数的理论,并进而不断发展成为一个有力的科研分析工具。
假定周期性信号周期T逐渐变大,则谱线间间隔将逐渐变小,如果外推周期T无限放大,变成无穷大,则信号或者函数就变成非周期信号或函数了,此时谱线就变成连续的了,而非一根一根离散的谱线!那么傅立叶变换正是这种一般性的数学定义:
对于连续时间信号f(t),若f(t)在时间维度上可积分,(实际上并不一定是时间t维度,这里可以是任意维度,只需在对应维度空间可积分即可),即:
∫−∞∞∣f(t)∣dt<∞\int_{-\infty}^{\infty}|f(t)|dt\lt\infty ∫−∞∞∣f(t)∣dt<∞
那么,x(t)的傅立叶变换存在,且其计算式为:
F(ω)=F[f(t)]=∫−∞∞f(t)e−jωtdtF(\omega)=\mathcal{F}[f(t)]=\int_{-\infty}^{\infty}f(t)e^{-j\omega t}dt F(ω)=F[f(t)]=∫−∞∞f(t)e−jωtdt
其反变换为:
f(t)=F−1[F(ω)]=12π∫−∞∞F(jω)e−jωtdωf(t)=\mathcal{F}^{-1}[F(\omega)]=\frac{1}{2\pi}\int_{-\infty}^{\infty}F(j\omega)e^{-j\omega t}d\omega f(t)=F−1[F(ω)]=2π1∫−∞∞F(jω)e−jωtdω
前文说傅立叶变换本质上也是一种连续函数的积分变换,那么从上面公式,可以看出傅立叶变换的核函数为:
K(t,ω)=ejωtK(t,\omega)=e^{j\omega t} K(t,ω)=ejωt
其核函数的两个自变量为t,ω\omegaω ,对于ω\omegaω一般称为角速度(可以形象的理解为旋转运动的快慢),是表征频率空间的。
上面这两个公式是啥意思呢?在度量空间可积可以理解成其在度量空间能量有限,也即对其自变量积分(相当于求面积)是一个确定值,那么这样的函数或者信号就可以进行傅立叶变换展开,展开得到的F(jω)F(j\omega)F(jω)就变成是频域的函数了,如果对频率ω\omegaω将函数值绘制出曲线就是我们所说的频谱图,而其逆变换就比较好理解了,如果我们知道一个信号或者函数谱密度函数F(jω)F(j\omega)F(jω),就可以对应还原出其时域的函数,也能绘制出时域的波形图。
傅立叶变换公式,从理解的角度,可以看成无限多无穷小的能量之和,而傅立叶级数也是各谐波分量的加和,所不同的是,前者相对于频率变量是连续的,而后者相对于频率则是离散的!
当然,本文限定讨论时域信号是因为我们电子系统中的应用最为普遍的就是一个时域信号。推而广之,其他的多维度信号也能利用上面定义进行推广,同样在多维空间信号也非常有应用价值,比如2维图像处理、3维图像重建等等。
傅立叶级数与变换的区别?
- 傅立叶级数对应的是周期信号,而傅立叶变换则对应的是一个时间连续可积信号(不一定是周期信号)
- 傅立叶级数要求信号在一个周期内能量有限,而后者则要求在整个区间能量有限
- 傅立叶级数的对应ω\omegaω是离散的,而傅立叶变换则对应ω\omegaω是连续的。
故而,两者的物理含义不同,且其量纲也是不同的,F(jkω)F(jk\omega)F(jkω)代表周期信号的第k次谐波幅度的大小,而F(jω)F(j\omega)F(jω)则是频谱密度的概念。所以答案是这两者从本质上不是一个概念,傅立叶级数是周期信号的另一种时域的表达方式,也就是正交级数,它是不同的频率的波形的时域叠加。而傅立叶变换则是完全的频域分析,傅里叶级数适用于对周期性现象做数学上的分析,傅里叶变换可以看作傅里叶级数的极限形式,也可以看作是对周期现象进行数学上的分析,同时也适用于非周期性现象的分析。
什么是拉普拉斯变换?
1814年法国数学家Pierre-Simon Laplace在研究概率论中给出了拉普拉斯的可靠数学依据,从而发展成拉普拉斯变换理论。对于函数f(t)我们知道其傅立叶变换为:
F(ω)=F[f(t)]=∫−∞∞f(t)e−jωtdtF(\omega)=\mathcal{F}[f(t)]=\int_{-\infty}^{\infty}f(t)e^{-j\omega t}dt F(ω)=F[f(t)]=∫−∞∞f(t)e−jωtdt
那么如果对于函数f(t)e−σf(t)e^{-\sigma}f(t)e−σ其傅立叶变换为:
F(ω)=F[f(t)e−σ]=∫0∞f(t)e−σe−jωtdtF(\omega)=\mathcal{F}[f(t)e^{-\sigma}]=\int_{0}^{\infty}f(t)e^{-\sigma}e^{-j\omega t}dt F(ω)=F[f(t)e−σ]=∫0∞f(t)e−σe−jωtdt
上面的公式整理一下:
F(ω)=F[f(t)e−σ]=∫0∞f(t)e−(σ+jωt)dtF(\omega)=\mathcal{F}[f(t)e^{-\sigma}]=\int_{0}^{\infty}f(t)e^{-(\sigma+j\omega t)}dt F(ω)=F[f(t)e−σ]=∫0∞f(t)e−(σ+jωt)dt
令s=σ+jωs=\sigma+j\omegas=σ+jω,则上面的变换
F(s)=∫0∞f(t)e−stdtF(s)=\int_{0}^{\infty}f(t)e^{-st}dt F(s)=∫0∞f(t)e−stdt
从前文我们知道,拉普拉斯本质上也是一种积分变换,那么上面公式,将e−ste^{-st}e−st看成积分变换的核函数,则其变换核函数为:
K(t,s)=e−stK(t,s)=e^{-st} K(t,s)=e−st
上面引入的因子e−σe^{-\sigma}e−σ,对于函数f(t)e−σf(t)e^{-\sigma}f(t)e−σ函数将变得更容易收敛,傅立叶变换的绝对可积分的限制条件也就更容易满足了。拉普拉斯变换存在的条件为:
limx−>∞f(t)e−σt=0(σ>σ0)\lim_{x->\infty}{f(t)e^{-\sigma t}}=0(\sigma \gt \sigma_0) x−>∞limf(t)e−σt=0(σ>σ0)
傅立叶拉氏变换联系区别
所以傅立叶变换与拉普拉斯变换的联系就比较容易联系了。
- 拉普拉斯变换,将原函数从时间维度(不一定是时间维度,只是方便理解本文以常见的时间维度信号进行描述),映射为复平面s=σ+jωs=\sigma + j\omegas=σ+jω
- 傅立叶变换是拉普拉斯变换的特例,也即变换核函数σ=0\sigma=0σ=0时,拉普拉斯变换就变成傅立叶变换了。相当于只取虚部,实部为0.
- 傅立叶变换是从原维度变换为频率维度,对于信号处理而言相当于将时域信号变换为频域进行分析,为信号处理提供了强大的数学理论基础及工具。
- 拉普拉斯变换,将原维度变换为复频域,在电子电路分析以及控制理论中,为建立系统的数学描述提供了强大的数学理论基础,学过控制理论的一天到晚都与传递函数打交道,其本质就是拉普拉斯变换对系统的一种数学建模描述。为分析系统的稳定性、可控性提供了数学工具。
什么是Z变换?
Z变换本质上是拉普拉斯变换的离散形式。也称为Fisher-Z变换。对于连续信号进行抽样变换就得到了原函数的离散序列:
fs(t)=f(t)⋅δT(t)=∑n=0∞f(nT)δ(t−nT)f_s(t)=f(t)\cdot \delta_T(t)=\sum_{n=0}^{\infty}f(nT)\delta(t-nT) fs(t)=f(t)⋅δT(t)=n=0∑∞f(nT)δ(t−nT)
其中T为采样周期,δT(t)\delta_T(t)δT(t)信号与系统中称为冲激抽样。其实说人话,就是将连续信号,按等间隔理想的转为抽取离散序列样本。看下图就明白了,在电子系统中常用AD转换器进行实现。
对上式进行拉普拉斯变换:
Fs(S)=∫0∞[∑n=0∞f(nT)δ(t−nT)]e−stdtF_s(S)=\int_0^\infty[\sum_{n=0}^{\infty}f(nT)\delta(t-nT)]e^{-st}dt Fs(S)=∫0∞[n=0∑∞f(nT)δ(t−nT)]e−stdt
该公式利用冲激函数的抽样特性,可简化为:
Fs(S)=∑n=0∞f(nT)e−snTF_s(S)=\sum_{n=0}^{\infty}f(nT)e^{-snT} Fs(S)=n=0∑∞f(nT)e−snT
引入z=esTz=e^{sT}z=esT,引入新的自变量Z,则上面的公式就变成这样了:
F(Z)=∑n=0∞f(nT)Z−nF(Z)=\sum_{n=0}^{\infty}f(nT)Z^{-n} F(Z)=n=0∑∞f(nT)Z−n
这就是Z变换了,从上面的过程描述就知道Z变换与拉普拉斯变换的关系了。因此两者的联系也就是Z变换是拉布拉斯变换的离散形式。
那么Z变换的意义在于什么呢?在数字信号处理以及数字控制系统中,Z变换提供了数学基础。利用Z变换很快就能将一个传递函数描述成差分方程形式,这就为编程实现提供了数学依据,比如一个数字滤波器知道其Z变换形式,写代码就是分分钟的事情了,同样知道一个控制算法的Z变换形式,同样编代码也是水到渠成的事情。
这里谈到Z变换的离散形式,那么这里也提一句,傅立叶变换数字落地,也即离散形式是离散傅立叶变换DFT(Discrete Fourier Transform),而大家所熟知的快速傅立叶变换FFT(Fast Fourier Transform)则是DFT的高效率实现。
总结一下
要理解三种变换的联系区别,首先要理解什么是数学变换,什么是积分变换。傅立叶变换以及拉普拉斯变换本质上都是连续或有限个第一类间断点函数的积分变换,而傅立叶变换是拉普拉斯变换的特殊形式,而Z变换是拉普拉斯变换的离散形式。每种变换都有其应用价值,傅立叶变换在信号处理的频域分析中提供了强大的数学工具,而拉普拉斯变换在电子学、控制工程、航空航天等领域提供了建模、分析的数学分析工具;Z变换则将这些变换进而落地为数字实现提供数学理论依据。DFT为FFT的离散化形式,而FFT是DFT的算法优化实现。
傅里叶变换、拉普拉斯变换、Z 变换的联系是什么?为什么要进行这些变换?相关推荐
- 简述计算机三大变换的联系和区别 (傅里叶变换 拉普拉斯变换 z变换)
Q:简述计算机三大变换的联系和区别 (傅里叶变换 拉普拉斯变换 z变换) (1) 傅里叶变换定义: 表示能将满足一定条件的某个函数表示成三角函数(正弦和/或余弦函数)或者它们的积分的线性组合.傅立叶变 ...
- 信号与系统 傅里叶变换 拉普拉斯变换 z变换所有公式和性质 三个变换的联系 整理
这是我考研整理的笔记.基本上涵盖了信号与系统三大变换所有重要的公式. 1.傅里叶变换 2.拉普拉斯变换 3.Z变换 4.三大变换的关系
- 傅里叶变换拉普拉斯变换的物理解释及区别
傅里叶变换在物理学.数论.组合数学.信号处理.概率论.统计学.密码学.声学.光学.海洋学.结构动力学等领域都有着广泛的应用(例如在信号处理中,傅里叶变换的典型用途是将信号分解成幅值分量和频率分量). ...
- laplace变换 z变换本质
laplace变换: 对一个实变量函数作拉普拉斯变换,并在复数域中作各种运算,再将运算结果作拉普拉斯反变换来求得实数域中的相应结果,往往比直接在实数域中求出同样的结果在计算上容易得多. 可把微分方程化 ...
- matlab实现阶跃响应变换,z域的阶跃响应matlab
离散系统的的Z域分析及Matlab实现_计算机软件及应用_IT/计算机_专业资料 暂无评价|0人阅读|0次下载|举报文档 离散系统的的Z域分析及Matlab实现_计算机软件及应用_...... 201 ...
- 计算机图形学三维变换论文,计算机图形学 第5章 三维图形生成和变换技术
计算机图形学 第5章 三维图形生成和变换技术 (63页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 19.90 积分 第五章 目录第五章 三维图形生成和变 ...
- 变换例题_小学语文学习攻略9:句式变换知识点概述+例题讲练
很多家长都会遇到这样的情况:孩子的语文考试卷发下来,作文上老师的评语或者做的记号上总会有几个是句意不通.在平常的交流中也是,支吾了半天也不能完整.通顺的表达自己的意思.这都是孩子在造句方面有障碍!孩子 ...
- c语言二维图形变换程序,C语言图形编程(五、二维图形变换-01)(国外英文资料).doc...
C语言图形编程(五.二维图形变换-01)(国外英文资料) C语言图形编程(五.二维图形变换-01)(国外英文资料) The first user coordinates to the screen c ...
- 多视图几何——变换层次总结(射影变换,仿射变换,相似变换,欧式变换)
书中的2D变换总结: 书中的3D变换总结: 变换层级 欧式->相似->仿射->射影,变换的层级提高,失真越来越严重,不变性质越来越少,变换矩阵的自由度越来越高.后面的变换都兼容前面的 ...
- python 希尔伯特变换_Python在信号与系统中的应用(1)——Hilbert变换,Hilbert在单边带包络检波的应用,FIR_LPF滤波器设计,还有逼格高高的FM(PM)调制...
多谢董老师,董老师是个好老师! 心情久久不能平静,主要是高频这门课的分析方法实在是让我难以理解,公式也背不过,还是放放吧. 最近厌恶了Matlab臃肿的体积和频繁的读写对我的Mac的损害,所以学习了一 ...
最新文章
- java内存分配--引用
- python爬虫,爬取猫眼电影2(xpath和bs4)
- 七、Linux脚本进阶和进程管理
- post json后台处理数据_Python爬虫教程-07-post介绍(百度翻译)(上)
- 山东大学linux实验四CSDN,山东大学操作系统实验报告材料4进程同步实验
- Cisco 4507R+E四引擎VSS故障解决
- web.config中的globalization 标签在将几个不同服务器上的网闸内容以统一面貌集成在自己网站上时的要考虑的问题
- java类型提升_java表达式中类型的自动提升(转)
- oracle 生成 sql语句,Oracle使用SQL语句生成日历的实现方法
- (转)word2vec前世今生
- HTML判断如果并且,将html转换为js,并且其中含有判断语句
- 递归处理汉诺塔问题(c++/python)
- TensorFlow应用实战-18-Policy Gradient算法
- java开发和安卓开发_Android开发和Java开发有什么区别
- 联想Y9000P安装Ubuntu20.04记录
- 唐宇迪学习笔记8:逻辑回归算法
- Visual Studio 2010下载 + 附破解方法
- 联想ghost重装系统_联想笔记本重装系统图文教程
- 完美解决浏览器主页被hao123劫持,打开浏览器时自动进入hao123主页怎么办
- Mac突然中文输入法框不见了??
热门文章
- MATLAB遗传算法求函数最值
- 热烈祝贺小明SEO博客网站正式上线
- IEEE1588v2解析(7)gPTP协议和PTP的关系
- java拆分pdf_itextpdf拆分pdf问题
- python正弦函数幂级数展开_函数展开成正弦级数或余弦级数的奇延拓
- 逻辑回归(LR)算法面试问题汇总
- 三大开源bi工具_不用无限手套,人人都能开发 BI 系统------开源BI软件,可集成的BI工具...
- cad图纸比对lisp_cad图纸怎么找出差异?教你怎么对比CAD图纸版本差异
- 【arp】关于arp和arping命令的使用
- mysql 批量录入 id_Mysql 批量插入数据 提前获得主键id