用虚数i与欧拉公式来解释分数阶微积分

article/2025/10/22 3:51:01

ps:如果研究又遇到分数阶相关课题的童鞋可以与博主联系获取分数阶相关的源代码,不图什么,就图个点赞

一、问题简介

在最近的课题中接触到了分数阶导数相关的东西,比如函数 f ( t ) f(t) f(t)对t的0.5阶导等等,网上有一些对于分数阶导数的处理办法是将其放到频域中进行处理,根据大学《自动控制原理》等课程我们可以知道,通过拉氏变换我们可以将微积分运算变为代数运算,对一个输入信号进行一阶微分就是乘以s,一阶积分就是乘以 1 / s 1/s 1/s,因此对于0.5阶的微分运算可以看作 s 0.5 s^{0.5} s0.5,因此分数阶微积分运算可以看作是一个传递函数,通过整数阶的以 s s s为自变量的方程可以逼近s的分数阶方程,现在通过这种逼近模型得到了 s 0.5 s^{0.5} s0.5 s 0.3 s^{0.3} s0.3两种运算的传递函数,利用Matlab观察到两种传递函数的bode图如下:
在这里插入图片描述
在这里插入图片描述
上述两幅图中的上半部分是传递函数的幅频特性曲线,下半部分是传递函数的相频特性曲线。在课程《自动控制原理》频率响应一章中我们可以知道将拉氏变换中的 s s s替换为 j w jw jw j j j是虚数)之后得到一个新的与 w w w传递函数F(W),对于一个角频率为 w w w的输入正弦信号,系统会对这个正弦信号产生幅值和相位的变化,这就对应着F(W)的模值与相位,同样也可以得到幅频曲线和相频曲线。在分数阶导数中,以0.5阶导 s 0.5 s^{0.5} s0.5为例,将 s s s换成 j w jw jw得到 ( j w ) 0.5 = j 0.5 ∗ w 0.5 (jw)^{0.5}=j^{0.5}*w^{0.5} (jw)0.5=j0.5w0.5,此时出现了 j 0.5 j^{0.5} j0.5,这样分数阶导数的问题就变成了复数的分数阶次方的问题,但是出现虚数的0.5次方时应该如何求解其幅值和相角来验证由整数阶逼近模型得到的bode图是否正确呢?通过一系列的资料查阅后得到了两种解释虚数的分数阶次方 j a j^a ja(a为实数)的方法。

二、从虚数的定义说起

1、虚数的起源与定义
虚数i的定义其实十分简单,在我们的中学校园里就有学过,虚数i的平方等于-1,虚数就是-1的平方根。但是我们在学习过程中一直都不知道虚数究竟代表着什么,在虚数出现前的数学中有4中基本运算:+-×÷,假如对于一个实轴上的点,加减运算意味着点的向左向右的平移,而乘除可以看做是对轴线的拉伸与压缩(也可以看做是一种累积的连续作用),但是这些运算都没有解决旋转的表示,在虚数发明之前人们对于旋转的处理都显得十分麻烦,而虚数与复数的出现就使得旋转的处理大大简化。

前面说到 i 2 = − 1 i^2=-1 i21,既可以写作 1 × i × i = − 1 1×i×i=-1 1×i×i1,乘法可以看做是一种连续的变化,那么可以想象,对于实轴上的点1经过了两次乘以i的变化变为了-1,相信大家已经知道了i代表着什么,没错,答案就是旋转,逆时针的旋转!如下图所示:
在这里插入图片描述
一个i就代表着逆时针旋转90°,而1经过两次乘以i的变换既逆时针旋转了180°,那么1就变成了-1。那么总结一下就是i代表了旋转,i是一个旋转量,i并不是一个数。
2、复数的引入
此时已经知道了i是一个旋转量,此时将实数作为横轴,虚数作为纵轴,并用加号将两种数连接起来就形成了复数 a + b i a+bi abi(也可以表示为 A ∠ a A∠a Aa),复数拥有着和虚数很相似的性质。比如两个复数的相乘就等于其模值相乘,相角求和,其实就是模值相乘同时在第一个复数的相角基础上逆时针继续旋转相应的角度。
3.虚数的分数阶解释
在理解了虚数和复数的定义及意义后我们再来看 s 0 . 5 s^0.5 s0.5次方,将 s s s替换为 j w jw jw后得到 ( j w ) 0.5 = j 0.5 ∗ w 0.5 (jw)^{0.5}=j^{0.5} *w^{0.5} (jw)0.5j0.5w0.5,由于j是一个虚数,是一个旋转量,因此 j 0.5 j^{0.5} j0.5的模值为1,因此 j w jw jw…的模值就等于 w 0.5 w^{0.5} w0.5。又由于乘以j表示旋转了90°,又 j = j 0.5 × j 0.5 j=j^{0.5}×j^{0.5} jj0.5×j0.5,既经过连续的两次 j 0.5 j^{0.5} j0.5的变换后会逆时针旋转90°,既 j 0.5 j^{0.5} j0.5代表了旋转45°。因此我们可以得出结果, s 0.5 s^{0.5} s0.5的相角为恒定的45°,与上图中的44.8°基本一致,其幅值增益为 w 0.5 w^{0.5} w0.5,将增益写作对数形式则表示为 10 l o g w 10logw 10logw,因此如图一,当 w = 100 w=100 w100是增益为20,与计算结果一致,最终判断逼近模型是正确的。

三、通过欧拉公式的计算来解释

欧拉公式也是我们在大学就学过的一个公式了,但我们对欧拉公式的了解也仅仅局限于公式之上,在下一次的公众号文章中我将会和大家讨论关于欧拉公式的起源与意义。在接下来的推导中我们以 s a s^a sa(a为实数)为例,将s替换为iw后得到 ( i w ) a = i a ∗ w a (iw)^a=i^a*w^a (iw)aiawa(这里的j用i来表示),此时通过欧拉公式来导出 i a i^a ia
首先我们知道欧拉公式的表达式为 e i x = cos ⁡ x + i ∗ sin ⁡ x e^{i x}=\cos x+i * \sin x eix=cosx+isinxi表示虚数,此时取 x = π / 2 x=π/2 x=π/2,则 c o s x = 0 cosx=0 cosx=0 s i n x = 1 sinx=1 sinx=1,因此得到公式: e π i 2 = i e^{\frac{\pi i}{2}}=i e2πi=i此时 i a i^a ia可以替换为 e π a i / 2 e^{πai/2} eπai/2,又将 e π a i / 2 e^{πai/2} eπai/2代入欧拉公式得到:
i a = e π a i 2 = cos ⁡ π a 2 + i ∗ sin ⁡ π a 2 = [ ( cos ⁡ π a 2 ) 2 + ( ( sin ⁡ π a 2 ) 2 ] ∠ π a 2 = ∠ π a 2 i^{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+isin2πa=[(cos2πa)2+((sin2πa)2]2πa=2πa
当a=0.3时, s 0.3 = w 0.3 ∗ i 0.3 = w 0.3 ∠ 27 ° s^{0.3}=w^{0.3}*i^{0.3}=w^{0.3}∠27° s0.3=w0.3i0.3=w0.327°,相角恒定为27°与第二个bode图的相频曲线一致,而对于幅频特性,增益为 20 l o g ( w 0.3 ) = 6 l o g ( w ) 20log(w^{0.3})=6log(w) 20log(w0.3)=6log(w)此时取w=100,即增益为6*2=12,与逼近模型得到的幅频曲线一致,最终验证得到分数阶逼近模型正确。


http://chatgpt.dhexx.cn/article/i1WJet3E.shtml

相关文章

常用的积分方法讨论(数学表达与代码整理)(龙格-库塔、中值积分、欧拉积分)

积分方法讨论(数学表达与代码整理) 数学原理 1.1 四元数与角速度的关系 在无人机或无人车的导航系统中常常采用四元数代替欧拉角来表示机体的旋转,因为欧拉角在计算过程中容易产生奇异,这与欧拉角的计算需要利用正弦、余弦公式…

常微分方程的解法 (二): 欧拉(Euler)方法

上一节讲了 常微分方程的三种离散化 方法:差商近似导数、数值积分、Taylor 多项式近似。 目录 2 欧拉(Euler)方法 2.1 向前 Euler 公式、向后 Euler 公式 2.2 Euler 方法的误差估计 3 改进的 Euler 方法 3.1 梯形公式 …

微积分 --- 欧拉数e的计算方法(个人学习笔记)

计算方法1: 计算方法2: 对于100来说,分37份的话,其值最接近e,且所有份的乘积最大为9.2944e15。 下面是从微积分的角度去求证,如果要让y为最大值,应该让xc/e,这就是最优份数。 &am…

考研数二第十七讲 反常积分与反常积分之欧拉-泊松(Euler-Poisson)积分

反常积分 反常积分又叫广义积分,是对普通定积分的推广,指含有无穷上限/下限,或者被积函数含有瑕点的积分,前者称为无穷限广义积分,后者称为瑕积分(又称无界函数的反常积分)。 含有无穷上限/下…

#欧拉第二积分(伽马函数)

伽玛函数,也叫欧拉第二积分,是阶乘函数在实数与复数上扩展的一类函数。 伽玛函数作为阶乘函数的延拓,是定义在复数范围内的亚纯函数,通常写成,负整数和0是它的一阶极点。 实数域: 复数域: 在解…

常见的数值积分方法 (欧拉、中值、龙格-库塔,【常用于IMU中】)

1. 积分基本概念 设F(x)为函数f(x)的一个原函数,我们把函数f(x)的所有原函数F(x)C(C为任意常数)叫做函数f(x)的不定积分(indefinite integral)。 非线性微分方程: 在有限的时间间隔Δt积分: 连续时间内积分: 工程上最常见的有三种&#xff…

欧拉积分法

数值积分法是求定积分的近似值的数值方法。即用被积函数的有限个抽样值的离散或加权平均近似值代替定积分的值。 数值积分法也是计算机仿真中常用的一种方法。在已知函数的微分方程时,求解函数下一时刻的值,我们主要有欧拉法、梯形法和龙格库塔法。 欧拉…

欧拉积分

欧拉积分 两个公式: $\Gamma(s)\int_{0}^{\infty}x^{s-1}e^{-x}dx,s>0$ (1) $B(p,q)\int_{0}^{1}x^{p-1}(1-x)^{q-1}dx,p>0,q>0$ (2) 一 、$\Gamma(Gamma)$函数 性质: 1.$\Gamma(s)$在定义域…

Oracle现使用CVSS 3.0对漏洞进行评级

Oracle今年4月关键补丁更新(Critical Patch Update)涉及多款产品中的136个漏洞,其中最大的变化是切换到通用安全漏洞评分系统3.0版本或者说CVSSv3,该版本可更准确反映漏洞带来的影响。 Oracle公司在其补丁公告中指出,这个关键补丁更新中的漏洞…

Cvss v2 complete documentation

通用漏洞计分系统(CVSS)为沟通IT漏洞的特征和影响提供了一个开放的框架。 CVSS由3组组成:基础,时间和环境。 每个组产生的范围从0到10的数字分数,以及Vector,一个反映用于得出分数的值的压缩文本表示。 基础…

东莞dell服务器维修上门服务,CVSS 10分漏洞影响Dell Wyse Thin客户端设备

近日,CyberMDX 研究人员公开了今年6月在Dell Wyse Thin客户端中发现了2个安全漏洞,漏洞CVE编号为CVE-2020-29491 和 CVE-2020-29492,这两个漏洞CVSS 评分都为10分,漏洞影响运行ThinOS v8.6及更低版本的所有设备。攻击者利用这两个…

通用漏洞评估方法CVSS3.0简表

CVSS3.0计算分值共有三种维度: 1. 基础度量。 分为 可利用性 及 影响度 两个子项,是漏洞评估的静态分值。 2. 时间度量。 基础维度之上结合受时间影响的三个动态分值,进而评估该漏洞的动态分值。 3. 环境度量。 根据用户实际环境需求结合时间…

通用漏洞评估方法CVSS3.0详解

CVSS(Common Vulerability Scoring System, 通用漏洞评估方法),是由NIAC 发布、FITST维护的开放式行业标准,CVSS 的发布为信息安全产业从业人员交流网络中所存在的系统漏洞的特点与影响提供了一个开放式的评价方法。 1.度量(Metrics&#xf…

通用漏洞评估方法CVSS3.0介绍

CVSS(Common Vulerability Scoring System, 通用漏洞评估方法),是由NIAC 发布、FITST维护的开放式行业标准,CVSS 的发布为信息安全产业从业人员交流网络中所存在的系统漏洞的特点与影响提供了一个开放式的评价方法。 1.度量(Metrics&#xf…

通用漏洞评估方法CVSS 3.0 计算公式及说明

CVSS 3.0 计算公式及说明 一、基础评价 1. 基础评价公式为: 当 影响度分值 <= 0: 基础分值 = 0 当 0 < 影响度分值 + 可利用度分值 < 10: 作用域 = 固定: 基础分值 = Roundup(影响度分值 + 可利用度分值) 作用域 = 变化: 基础分值 = Roundup[1.08 (影响度分值 + 可…

漏洞评估-CVSS3

详细可以参考ITU-T X.1521 Exploitability Attack Vector&#xff08;AV&#xff09; Attack Complexity(AC) Priviliages required(PR) User Interation&#xff08;UI&#xff09; Scope(S) 范围指的是计算授权主体&#xff08;如应用、操作系统或沙盘环境&#xff09;在授…

基于CVSS3.1的一种评估框架

原文 Vulnerability Modelling for Hybrid Industrial Control System Networks 出版 Springer Nature B.V. 2020June 2020https://doi.org/10.1007/s10723-020-09528-w 文章目录 一、摘要二、相关技术介绍&#xff08;一&#xff09;CVSS3.1计算&#xff08;二&#xff09;…

【CVSS V3.1漏洞评分计算方法】

​ 漏洞计算的官网CVSSV3.1 &#xff1a;Common Vulnerability Scoring System Version 3.1 Calculator 1、使用方法&#xff1a; 用鼠标移到对应的字符即可了解相应的含义。选中后即可得到一个分值。2、严酷度评估&#xff1a; 严重程度分值范围9.0-10致命7.0-8.9严重4.0-6…

【文献翻译】基于CVSS的IT系统网络安全风险定量评估方法-A Quantitative CVSS-Based Cyber Security Risk Assessment Methodology

基于CVSS的IT系统网络安全风险定量评估方法 A Quantitative CVSS-Based Cyber Security Risk Assessment Methodology For IT Systems 摘要 由于我们不断增长的IT系统中网络威胁不断增加&#xff0c;IT系统风险评估是必不可少的。此外&#xff0c;法律法规敦促组织定期进行风…

cvss评分及漏洞矢量

CVSS CVSS全称为Common Vulnerability Scoring System&#xff0c;即“通用漏洞评分系统”&#xff0c;是一个行业公开的标准。其被设计用来评测漏洞的严重程度&#xff0c;并帮助确定所需反应的紧急度和重要度。通过漏洞难易程度以及对机密性、完整性、可用性的影响综合评估后…