用虚数i与欧拉公式来解释分数阶微积分
ps:如果研究又遇到分数阶相关课题的童鞋可以与博主联系获取分数阶相关的源代码,不图什么,就图个点赞
一、问题简介
在最近的课题中接触到了分数阶导数相关的东西,比如函数f(t)f(t)f(t)对t的0.5阶导等等,网上有一些对于分数阶导数的处理办法是将其放到频域中进行处理,根据大学《自动控制原理》等课程我们可以知道,通过拉氏变换我们可以将微积分运算变为代数运算,对一个输入信号进行一阶微分就是乘以s,一阶积分就是乘以1/s1/s1/s,因此对于0.5阶的微分运算可以看作s0.5s^{0.5}s0.5,因此分数阶微积分运算可以看作是一个传递函数,通过整数阶的以sss为自变量的方程可以逼近s的分数阶方程,现在通过这种逼近模型得到了s0.5s^{0.5}s0.5和s0.3s^{0.3}s0.3两种运算的传递函数,利用Matlab观察到两种传递函数的bode图如下:
上述两幅图中的上半部分是传递函数的幅频特性曲线,下半部分是传递函数的相频特性曲线。在课程《自动控制原理》频率响应一章中我们可以知道将拉氏变换中的sss替换为jwjwjw(jjj是虚数)之后得到一个新的与www传递函数F(W),对于一个角频率为www的输入正弦信号,系统会对这个正弦信号产生幅值和相位的变化,这就对应着F(W)的模值与相位,同样也可以得到幅频曲线和相频曲线。在分数阶导数中,以0.5阶导s0.5s^{0.5}s0.5为例,将sss换成jwjwjw得到(jw)0.5=j0.5∗w0.5(jw)^{0.5}=j^{0.5}*w^{0.5}(jw)0.5=j0.5∗w0.5,此时出现了j0.5j^{0.5}j0.5,这样分数阶导数的问题就变成了复数的分数阶次方的问题,但是出现虚数的0.5次方时应该如何求解其幅值和相角来验证由整数阶逼近模型得到的bode图是否正确呢?通过一系列的资料查阅后得到了两种解释虚数的分数阶次方jaj^aja(a为实数)的方法。
二、从虚数的定义说起
1、虚数的起源与定义
虚数i的定义其实十分简单,在我们的中学校园里就有学过,虚数i的平方等于-1,虚数就是-1的平方根。但是我们在学习过程中一直都不知道虚数究竟代表着什么,在虚数出现前的数学中有4中基本运算:+-×÷,假如对于一个实轴上的点,加减运算意味着点的向左向右的平移,而乘除可以看做是对轴线的拉伸与压缩(也可以看做是一种累积的连续作用),但是这些运算都没有解决旋转的表示,在虚数发明之前人们对于旋转的处理都显得十分麻烦,而虚数与复数的出现就使得旋转的处理大大简化。
前面说到i2=−1i^2=-1i2=−1,既可以写作1×i×i=−11×i×i=-11×i×i=−1,乘法可以看做是一种连续的变化,那么可以想象,对于实轴上的点1经过了两次乘以i的变化变为了-1,相信大家已经知道了i代表着什么,没错,答案就是旋转,逆时针的旋转!如下图所示:
一个i就代表着逆时针旋转90°,而1经过两次乘以i的变换既逆时针旋转了180°,那么1就变成了-1。那么总结一下就是i代表了旋转,i是一个旋转量,i并不是一个数。
2、复数的引入
此时已经知道了i是一个旋转量,此时将实数作为横轴,虚数作为纵轴,并用加号将两种数连接起来就形成了复数a+bia+bia+bi(也可以表示为A∠aA∠aA∠a),复数拥有着和虚数很相似的性质。比如两个复数的相乘就等于其模值相乘,相角求和,其实就是模值相乘同时在第一个复数的相角基础上逆时针继续旋转相应的角度。
3.虚数的分数阶解释
在理解了虚数和复数的定义及意义后我们再来看s0.5s^0.5s0.5次方,将sss替换为jwjwjw后得到(jw)0.5=j0.5∗w0.5(jw)^{0.5}=j^{0.5} *w^{0.5}(jw)0.5=j0.5∗w0.5,由于j是一个虚数,是一个旋转量,因此j0.5j^{0.5}j0.5的模值为1,因此jwjwjw…的模值就等于w0.5w^{0.5}w0.5。又由于乘以j表示旋转了90°,又j=j0.5×j0.5j=j^{0.5}×j^{0.5}j=j0.5×j0.5,既经过连续的两次j0.5j^{0.5}j0.5的变换后会逆时针旋转90°,既j0.5j^{0.5}j0.5代表了旋转45°。因此我们可以得出结果,s0.5s^{0.5}s0.5的相角为恒定的45°,与上图中的44.8°基本一致,其幅值增益为w0.5w^{0.5}w0.5,将增益写作对数形式则表示为10logw10logw10logw,因此如图一,当w=100w=100w=100是增益为20,与计算结果一致,最终判断逼近模型是正确的。
三、通过欧拉公式的计算来解释
欧拉公式也是我们在大学就学过的一个公式了,但我们对欧拉公式的了解也仅仅局限于公式之上,在下一次的公众号文章中我将会和大家讨论关于欧拉公式的起源与意义。在接下来的推导中我们以sas^asa(a为实数)为例,将s替换为iw后得到(iw)a=ia∗wa(iw)^a=i^a*w^a(iw)a=ia∗wa(这里的j用i来表示),此时通过欧拉公式来导出iai^aia。
首先我们知道欧拉公式的表达式为 eix=cosx+i∗sinxe^{i x}=\cos x+i * \sin xeix=cosx+i∗sinxi表示虚数,此时取x=π/2x=π/2x=π/2,则cosx=0cosx=0cosx=0,sinx=1sinx=1sinx=1,因此得到公式:eπi2=ie^{\frac{\pi i}{2}}=ie2πi=i此时iai^aia可以替换为eπai/2e^{πai/2}eπai/2,又将eπai/2e^{πai/2}eπai/2代入欧拉公式得到:
ia=eπai2=cosπa2+i∗sinπa2=[(cosπa2)2+((sinπa2)2]∠πa2=∠πa2i^{a}=e^{\frac{\pi a i}{2}}=\cos \frac{\pi a}{2}+i * \sin \frac{\pi a}{2}=\left[\left(\cos \frac{\pi a}{2}\right)^{2}+\left(\left(\sin \frac{\pi a}{2}\right)^{2}\right] \angle \frac{\pi a}{2}=\angle \frac{\pi a}{2}\right.ia=e2πai=cos2πa+i∗sin2πa=[(cos2πa)2+((sin2πa)2]∠2πa=∠2πa
当a=0.3时,s0.3=w0.3∗i0.3=w0.3∠27°s^{0.3}=w^{0.3}*i^{0.3}=w^{0.3}∠27°s0.3=w0.3∗i0.3=w0.3∠27°,相角恒定为27°与第二个bode图的相频曲线一致,而对于幅频特性,增益为20log(w0.3)=6log(w)20log(w^{0.3})=6log(w)20log(w0.3)=6log(w)此时取w=100,即增益为6*2=12,与逼近模型得到的幅频曲线一致,最终验证得到分数阶逼近模型正确。
用虚数i与欧拉公式来解释分数阶微积分相关推荐
- 分数阶simulink工具箱_CCDC 2021特别专题:分数阶微积分与分数阶系统
点击蓝字,关注我们 征文通知 第33届中国控制与决策会议(CCDC 2021)将于2021年5月22日~24日在中国昆明举行.会议特别专题分数阶微积分与分数阶系统现向从事相关研究的广大专家.学者征稿. ...
- 基于分数阶的图像边缘细节检测
1 分数阶微分理论 分数阶微分几乎和整数阶微分同时诞生,但由于一直没有常见的物理现象能够解释这一数学表达式的含义,所以其也未被广泛运用.几个世纪以来,虽然分数阶微分理论得到了长足的发展,但至今没有统一 ...
- MATLAB求分数阶微分的数值解,G-L定义,R-L定义,Caputo定义
分数阶微积分学是整数阶微积分学的直接拓展,将一阶导数.二阶导数.一重积分.二重积分等整数阶微积分拓展到0.75阶导数.阶导数等实数甚至是复数阶的导数或积分.这无疑拓展了微积分学的深度. 对于整数阶微积 ...
- matlab 分数阶混沌系统的完全同步控制
1.内容简介 略 625-可以交流.咨询.答疑 2.内容说明 分数阶微积分这一重要的数学分支,其诞生在1695年,几乎和经典微积分同时出现.那一年,德国数学家Leibniz 和法国数学家L'Hopit ...
- caputo分数阶导数程序_基于Caputo导数的分数阶非线性振动系统响应计算
基于 Caputo 导数的分数阶非线性振动系统响应计算 李亚杰 1,2 , 吴志强 *1,2 , 章国齐 1,2 [摘 要] 摘 要 : 研究了含分数阶 Caputo 导数的非线性振动系统响应的数值计 ...
- 图像增强-分数阶微分(vc++)
分数阶微分,Grumwald-Letnikov定义在图像的数值实现中更为准确: Gamma函数: 若s(t)的持续期t [a,t],将函数持续期间[a,t]按单位间隔h=1进行等分,得到: 推到一元函 ...
- 【控制】基于灰狼算法改进分数阶PD滑模控制器附matlab代码
1 内容介绍 分数微积分已经被研究了将近 3 个世纪,并且已 经被科学家广泛应用到科学与控制工程领域中.分 数阶 PID 控制系统是由斯洛伐克学者 Podlubny于 1994 年提出,并应用于分数阶 ...
- 基于Duffing系统的分数阶混沌研究【基于matlab的动力学模型学习笔记_5】
/*仅当作学习笔记,若有纰漏欢迎友好交流指正,此外若能提供一点帮助将会十分荣幸*/ 前面的几篇博文我们提到提到的都是整数阶模型,这里我们将对分数阶模型进行一个简单的研究. 摘要:与整数阶混沌相比,分数 ...
- 基于C语言的分数阶微分方程,基于分数阶C~α空间的积分微分方程的配置解法
摘要: 分数阶积分微分方程是在处理实际问题时应运而生,在工程模型和物理现象中极其常见,比如化学反应扩散,弹性力学,热传导方面,种群生态模型,控制理论和生物化学等.近年来,随着科学技术的迅猛发展,分数阶 ...
- 论文总结2:分数阶傅里叶变换在信号检测与图像处理中的应用研究
分数阶傅里叶变换在信号检测与图像处理中的应用研究(李琼) 发展历史 特性 二维分数阶傅里叶变换 基于分数阶傅里叶变换的图像分析 分数阶变换域中图像的能量分布 分数阶傅里叶变换域中的图像的幅度和相位信息 ...
最新文章
- 一条关于互联网创业公司健康与否的判断法则
- python【蓝桥杯vip练习题库】PREV-52小数第n位(模拟)
- 第十八章 12判断string类型字符串是否为空
- java 蓝桥杯算法训练 整除问题(题解)
- 向前的快捷键_平面设计基础知识:平面设计师应该知道的快捷键。
- uniapp底部导航栏中间凸起_PPT导航栏,真的非常好用!
- xp电脑多少位怎么看_怎么看电脑是32位还是64位
- Liferay 6.2 改造系列之十三:修改用户编辑页面表单内容
- gmm聚类python_聚类算法GMM和KMeans?
- Html中的表格、Table和div,表格页面、员工登记表、Div代码
- 用神经网络实现手写数字识别
- 许愿墙 php源码,许愿网许愿墙 v1 共享版
- UVA 1449 Dominating Patterns
- Python爬虫selenium的使用实例爬取知乎首页数据!
- 象棋大战 v1.075 免费
- Error creating bean with name ‘sqlSessionFactory‘ defined in class path resource [applicationContext
- Linux删除文件,df查看磁盘空间未减少
- 【练习题】第二章--变量,表达式,语句(Think Python)
- 光棍必看的十部电影(ZT)
- excel表格换行快捷键_常州办公excel表格常用技巧, 办公软件常用快捷键【金坛吧】...
热门文章
- java drm_DRM系统工作原理
- 轻松解决U盘中病毒,文件变成.exe执行文件的问题
- 计算机控制课程设计体会,计算机控制技术课程设计心得体会
- android .dex文件,对于Android DEX文件详细说明
- java中 什么意思?比如130
- 读书随笔《36岁, 半熟人生》
- linux分析rna-seq,RNA-seq 分析流程(一)linux部分
- [血泪整理]为何加载EDEM和FLUENT耦合接口显示ERROR126
- 基于react hook的砸金蛋动画
- 10大H5前端ui框架,ui让你开发不愁