分数阶傅立叶变换 matlab,【综述】分数阶傅里叶变换(FRFT)
原标题:【综述】分数阶傅里叶变换(FRFT)
作者:WTT整理
傅里叶级数(傅里叶变换)几乎在所有科学和工程领域发挥着重要作用。黎曼积分和勒贝格积分均起源于对傅里叶级数的研究,傅里叶级数(傅里叶变换)在齿轮箱故障诊断分析上也有广泛应用,是目前齿轮箱技术研究的热点之一,特别是分数阶傅里叶变换(FRFT,fractional Fourier transform),在齿轮箱早期故障的诊断上,具有良好的发展潜力。
1. 历史回顾
三角函数展开起源于18世纪。1808年,傅里叶在他著名的热力学论文集“热的分析理论”中详细的研究了三角级数,并用三角级数成功的解决了许多热传导问题的偏微分方程。由于当时傅里叶的结论并无确凿证据及数学家对傅里叶的观点还很陌生,所以傅里叶的研究当时饱受争议,并且人们还认为他夸大了结论的适用范围。事实上,目前数学家已经证明了傅里叶级数足以表示绝大多数曲线。
将傅里叶级数的定义域拓展到整个实轴上,就得到了傅里叶变换。傅里叶变换是分析和处理平稳信号的一种标准和有力的工具,但随着研究范围的扩大逐渐也显示出了它的局限性:主要是傅里叶变换是全局变换,得到的是信号的整体频谱,因而无法表述时频的局部特征,因此人们进一步提出了分数阶傅里叶变换,短时傅里叶变换。
分数阶傅里叶变换最早是纳米亚(Namias)于1980年提出,其目的是用于求解量子力学中出现的线性时变偏微分方程。随后麦克布莱德(McBride)等对FRFT理论从数学上以积分形式给出了严格定义。1993年洛曼(Lohmann)阐述了FRFT的物理意义,即可理解为时频平面的旋转,洛曼(Lohmann)开创性的工作使得FRFT首先在光学中得到应用。
FRFT实际上是一种统一的时频变换,同时反映了信号在时域和频域上的信息,与传统的傅里叶变换相比,它适用于处理非平稳信号,尤其是线性调频信号,这使得FRFT在齿轮箱故障诊断上也得到越来越多的应用。
以下我们通过定义来了解它们之间的关系。
2. 傅里叶级数定义:
傅里叶级数就是将某一函数改写成正弦(sin)函数和余弦(cos)函数的和的形式。
2.1 傅里叶级数在[-π,л]区间上的形式为:
2.2傅里叶级数拓展到任意周期区间[-a,a]上的形式为:
2.3 将以上形式改写成复数形式为
如果我们将[-a.a]区间拓展到整个实轴上,就引出以下傅里叶变换定义。
3. 傅里叶变换定义:
我们用F表示傅里叶算子,如果对信号f(t)连续进行傅里叶变换,则有:
FF[f(t)]=F2[f(t)]=f(-t);
F3[f(t)]=f(-λ);
F4[f(t)]=f(t);
如果我们将t轴(时间轴)与λ(频率轴)轴构成一个直角坐标系,将每次的傅里叶变换看作坐标轴的90度旋转,在每次旋转的同时改变f(t)的表达形式,则这种坐标系定义的t-λ平面(可理解为时间-频率平面)满足以上傅里叶变换的性质。
图1傅里叶变换t-λ平面(时-频平面)表示
显然,时频平面的90度旋转对应一次傅里叶变换,那么时频平面的任意角度α旋转的数学定义和物理意义是什么呢?这就引出了以下的分数阶傅里叶变换的定义。
4. 分数阶傅里叶变换定义:
j为复数符号,p为傅里叶阶数。当p=1时,f1(u)就是f(t)的普通的傅里叶变换;当p=-1时,f-1(u)就是f(t)的普通的傅里叶逆变换;同时,FRFT满足阶数可加性(指数可加性,或旋转可加性),即FP1FP2=FP1+P2。所以可以认为分数阶傅里叶变换是一种广义的傅里叶变换。
5. 矩形波的FRFT展开:
通过下面矩形波不同阶次的展开图,能形象的理解信号从时域逐渐过渡到频域的过程。
f(t)=1,-4≤t<4;f(t)=0,t是其它值时。
以下图2到图6中,实线表示FRFT的实部,虚线表示FRFT的虚部。当p很小时,FRFT的实部非常接近时域的矩形方波(见图2);当p值接近于1时,FRFT的结果就接近矩形波的傅里叶变换,即sinc函数(见图6)。
图2 p=0.005时的FRFT
图3 p=0.1时的FRFT
图4 p=0.4时的FRFT
图5 p=0.6时的FRFT
图6 p=0.995时的FRFT
参考文献:
梅检民,肖云魁,齿轮箱早期故障精细诊断技术—分数阶傅立叶变换原理和 应用[M],高等教育出版社,2016
申永军,杨绍谱,齿轮系统的非线性动力学与故障诊断[M],科学出版社,2014
Albert Boggess,Francis J.Narcowich,小波与傅里叶分析基础(第二板)[M],芮国胜,等译,电子工业出版社,2017
E-mail:
Suzhou Walnut Transmission Technology Co.,ltd
Perfect Design Comes From Profession And Experience返回搜狐,查看更多
责任编辑:
分数阶傅立叶变换 matlab,【综述】分数阶傅里叶变换(FRFT)相关推荐
- 信号与系统仿真实验——实验二 傅立叶变换MATLAB的实现及傅里叶变换性质的分析
[ 实验目的] 1.利用MATLAB分析非周期信号的频谱 2.观察信号频谱变化验证傅里叶变换性质 [ 实验内容] [ 实验报告要求] (1)记录实验一和实验三中的波形: (2)总结实验二中频谱特性曲线 ...
- 分数阶傅立叶变换程序汇总(转载)
转载链接http://forum.vibunion.com/forum.php?mod=viewthread&action=printable&tid=41561 声振论坛 标题: 分 ...
- 分数S(Stockwell)变换Matlab代码
分数S变换(FrST) 离散的FrST Matlab代码 实验结果 离散的FrST Matlab代码 function FRST= FrST(t,sig,freqlow,freqhigh,alpha, ...
- 分数阶傅立叶变换中午matlab,怎么做短时分数阶傅里叶变换
我的分数阶傅里叶变换如下,怎么将它和短时傅里叶变换结合起来,变成短时分数阶傅里叶变换,有没有大佬指点一下下!!! frft: function Faf = frft(f, a) % The fast ...
- 高阶系统降阶matlab程序,高阶系统的降阶控制研究
引言近年来,主动控制技术在很多方面都得到了很大的发展,尤其在飞机设计中,这门新技术越来越为人们所重视.但是,对飞机实施主动控制很关键的一个问题就是其实时性,离开实时性,主动控制就没韦现实性可言.计算机 ...
- 数字信号处理FFT快速傅立叶变换MATLAB实现——实例
今天做作业的时候发现要对一个信号进行FFT变换,在网上找了半天也没找到个能看懂的(因为我太菜了),后来自己研究了一下,感觉一知半解的 起因是这道作业题 例题-满足奈奎斯特 我画了两个图,一个是原信号经 ...
- 分数傅立叶变换的性质
主要内容: 回顾--分数傅立叶变换基于积分核形式的推导 当P=1的时候,可以把分数傅里叶变换的积分核退化为傅里叶变换的积分核形式. 在坐标系当中,当旋转角度到90°的时候就退化为傅里叶变换,当旋转角度 ...
- 线性代数28——复矩阵和快速傅立叶变换
原文 | https://mp.weixin.qq.com/s/YzPoPnRb-gEm_EiV9et0TA 实矩阵也可能碰到复特征值,因此无可避免地在矩阵运算中碰到复数. 矩阵当然也有可能包含复数, ...
- 为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换?
https://www.douban.com/note/164400821/ 写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,绝大部分内容非我所原创.在此向多位原创作者致敬!!! ...
最新文章
- squid与varnish
- php 支付宝支付 简书,支付宝支付接口
- 交换机出现电源故障怎么解决?
- 程序员外包避坑指南?
- 线程间通讯《代码》pthread_cond_t
- Intel异常的分类:错误,陷阱,终止
- c++ 随机字符串_第3章 别碰白块(《C和C++游戏趣味编程》配套教学视频)
- Linux基础-5进程管理工具
- CCF201609-1 最大波动(100分)
- Kotlin教程(一):走进Kotlin的世界
- WLAN/WIFI信道列表 2020-11-23
- 2021年山东省安全员C证考试内容及山东省安全员C证模拟试题
- CMMB 全国各地市频率规划表
- 云存储可靠性的计算方法
- File Juicer for Mac(文件内容提取工具)
- 商品抢购时间倒计时结束后加入购物清单demo
- C#实现的系统内存清理
- 信号分析与处理——拉普拉斯变化
- SCOM2019导入管理包
- 泰尔解说Google云计算业务:关键在于垄断
热门文章
- python基于opencv的手势识别_求python opencv 手势识别源代码
- CCF的基于UEBA的用户上网异常行为分析baseline(线上0.9263)
- java-net-php-python-jsp学生社团信息演示录像2019计算机毕业设计程序
- U盘出现拒绝访问怎么解决
- 计算机二级讲座宣传语,计算机二级讲座策划书.doc
- ThinkPad E460如何进入bios
- 中秋节快到了,一起用MATLAB绘制一款2.5D月饼叭
- 推荐使用Windows10企业版LTSC的理由
- (三)JVM成神路之全面详解执行引擎子系统与JIT即时编译原理
- layui追加或动态修改表单元素不生效