一种变阶分数阶微积分时域分析方法与流程

本发明涉及电子信息
技术领域：
，尤其涉及一种变阶分数阶微积分时域分析方法。
背景技术：
：分数阶微积分学的发展起源于19世纪60年代，但其理论在近年来才被应用到工程技术中。目前主要应用领域包括了：涵盖粘弹性分析、反馈放大器、电容理论、控制领域、生物系统电导率分析、分数阶神经系统建模等领域。尤其是在一些自动化生产线中，会采用分数阶微积分优化自动控制流程，从而提高生产效率以及控制精度。在现阶段的方案中，通过irid_fod(impulseresponseinvariantdiscretization_fractionalorderderivative)方法得到分数阶微积分离散传递函数，其方法是利用拉普拉斯逆变换得到分数阶微积分算子的单位脉冲响应函数，并取响应函数的n个离散点，用prony方法来拟合n个离散点，获得z离散传递函数，但是在这种计算传递函数的方式，实时性较低，在对于实时性要求较高的应用领域中(比如自动控制领域)，限制了该方案的使用范围。技术实现要素：本发明的实施例提供一种变阶分数阶微积分时域分析方法，能够提高传递函数的计算效率。为达到上述目的，本发明的实施例采用如下技术方案：所述方法用于一种分析系统，所述分析系统用于为外部系统提供分数阶微积分计算，所述方法用于一种分析系统，所述分析系统用于为外部系统提供分数阶微积分计算，所述分析系统包括：计算机设备、利用pci接口连接所述计算机设备的quanser板卡(qpide)、所述quanser板卡与模拟传感器适配器电源和信号路由器组成的通信模块连接二级倒立摆；提供一种基于两点差值的变阶分数阶微积分时域分析方法包括：步骤(1)：对于分数阶积分算子所述计算机设备令阶次λ在(01)范围内取m个均匀点，并记为λi，其中i＝1,2,…,m；并且，对于分数阶微分算子所述计算机设备令阶次λ在(-10)范围内取m个均匀点，记为λi，其中i＝1,2,…,m，分数阶积分算子和分数阶微分算子的符号都为其中积分的运用范围是(0,1)，微分的运用范围是(-1,0)，s表示拉普拉斯变换后的算子；步骤(2)：离散传递函数将分数阶积分作为滤波器的时域离散传递函数,并求得单位脉冲输入δ(t)的情况下，所述滤波器输出的时域步骤(3)：根据所述滤波器的输入δ(t)和输出利用prony算法求出所述滤波器的时域离散传递函数；步骤(4)：根据所述时域离散传递函数，构造分母系数矩阵a；步骤(5)：根据所述时域离散传递函数，构造分子系数矩阵b；步骤(6)：以离散点阶次x＝[λ0,λ1,…λm]为输入，以分母系数矩阵a的列向量转置为输出，采用两点差值法进行拟合计算，得到n+1个分母的系数，其中，n表示大于等于0的整数，a[n,1],a[n,2],…a[n,m]代表矩阵a第一列到最后一列的数值；步骤(7)：以离散点阶次x＝[λ0,λ1,…λm]为输入，以分子系数矩阵b的列向量转置为输出，采用两点差值法进行拟合计算，得到n+1个分子的系数；步骤(8)：由所述分母系数的n+1多项式和所述分子系数的n+1多项式，得到分数阶微积分的时域整数阶传递函数。可选的，还提供一种基于多项式拟合的变阶分数阶微积分时域分析方法包括：步骤(1)：对于分数阶积分算子所述计算机设备令阶次λ在(01)范围内取m个均匀点，并记为λi，其中i＝1,2,…,m；并且，对于分数阶微分算子所述计算机设备令阶次λ在(-10)范围内取m个均匀点，记为λi，其中i＝1,2,…,m，分数阶积分算子和分数阶微分算子的符号都为其中积分的运用范围是(0,1)，微分的运用范围是(-1,0)，s表示拉普拉斯变换后的算子；步骤(2)：离散传递函数将分数阶积分作为滤波器的时域离散传递函数,并求得单位脉冲输入δ(t)的情况下，所述滤波器输出的时域步骤(3)：根据所述滤波器的输入δ(t)和输出利用prony算法求出所述滤波器的时域离散传递函数；步骤(4)：根据所述时域离散传递函数，构造分母系数矩阵a；步骤(5)：根据所述时域离散传递函数，构造分子系数矩阵b；步骤(6)：以离散点阶次x＝[λ0,λ1,…λm]为输入，以分母系数矩阵a的列向量转置为输出，采用最小二乘方法进行拟合计算，得到所述分母系数的n+1多项式，其中，n表示大于等于0的整数，a[n,1],a[n,2],…a[n,m]代表矩阵a第一列到最后一列的数值；步骤(7)：以离散点阶次x＝[λ0,λ1,…λm]为输入，以分子系数矩阵b的列向量转置为输出，采用最小二乘方法进行拟合计算，得到所述分子系数的n+1多项式；步骤(8)：由n+1个分母系数和n+1个分子系数，得到分数阶微积分的时域整数阶传递函数。本发明实施例提供的变阶分数阶微积分时域分析方法，采用脉冲响应不变原理，通过时域拟合获得多个等间隔阶次时分数阶微积分算子的时域离散传递函数；然后构造离散传递函数同位置系数矩阵，并以系数矩阵的列向量进行多项式拟合，获得阶次为输入、系数为输出的拟合函数组；最后对于任意阶次的分数阶微积分算子，可快速解算并构建出其离散传递函数，从而实现变阶次分数阶微积分的数值计算。仿真结果表明该方法具有较高的精度和快速性。因此，本发明提供的方案实现了快速求解分数阶传递函数，提高了传递函数的计算效率。也间接提高了在工业控制时的精度以及效率。相对于现有的irid_fod(impulseresponseinvariantdiscretization_fractionalorderderivative)irid_fod方法，本实施例提出的方法能够更快地运算出其传递函数。从而可以克服irid_fod方法中无法实时获得传递函数的缺点，提高了传递函数的计算效率。附图说明为了更清楚地说明本发明实施例中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。图1为本发明实施例提供的系统架构示意图；图2为本发明实施例提供的逻辑原理示意图；图3至图14为本发明实施例提供的y相对于x的散点图，其中离散的阶次x＝[α1,α2,…,αm]为横坐标,以为纵坐标；图15为本发明实施例提供的具体实例中给定的正弦波输入；图16通过irid_fod方法得到的波形图；图17为本发明实施例提供的两点法得到的波形图。具体实施方式为使本领域技术人员更好地理解本发明的技术方案，下面结合附图和具体实施方式对本发明作进一步详细描述。下文中将详细描述本发明的实施方式，所述实施方式的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。本
技术领域：
技术人员可以理解，除非特意声明，这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是，本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件，但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解，当我们称元件被“连接”或“耦接”到另一元件时，它可以直接连接或耦接到其他元件，或者也可以存在中间元件。此外，这里使用的“连接”或“耦接”可以包括无线连接或耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的任一单元和全部组合。本
技术领域：
技术人员可以理解，除非另外定义，这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样定义，不会用理想化或过于正式的含义来解释。本发明实施例提供一种变阶分数阶微积分时域分析方法，所述方法用于一种如图1所示的分析系统所述分析系统用于为外部系统提供分数阶微积分计算，所述分析系统包括：计算机设备、利用pci接口连接所述计算机设备的quanser板卡(qpide，一种多功能数据采集卡)、所述quanser板卡与模拟传感器适配器电源和信号路由器组成的通信模块连接二级倒立摆；其中辅助控制二级倒立摆的硬件设备还有多功能数据采集卡、模拟传感器适配器电源和信号路由器。所述计算机设备上运行的分析系统是依靠quarc通用快速控制实时仿真系统，其通过matalb中simulink设计的控制机，直接生成实时代码，并将工业级的实时应用程序下载到电脑操作系统中。其中，具体可以在计算机设备利用matlab仿真建模，将外部系统全部仿真到计算机设备中，建立控制系统的仿真模型。在matlab平台上编写步骤1到步骤7的程序(一种用于外部系统的控制算法)，在给定一个输入变量后，通过quanser板卡控制二级倒立摆。quanser板卡具体可以采用quanser公司的qpide卡。在本实施例中，二级倒立摆具体可以采用目前的倒立摆控制系统。倒立摆控制系统是一个一个复杂的、不稳定的、非线性系统的实验平台，主要用于控制理论教学及开展各种控制实验。对倒立摆系统的研究能有效的反映控制中的许多典型问题:如非线性问题、鲁棒性问题、镇定问题、随动问题以及跟踪问题等。通过对倒立摆的控制，用来检验新的控制方法是否有较强的处理非线性和不稳定性问题的能力。同时，其控制方法在军工、航天、机器人和一般工业过程领域中都有着广泛的用途，如机器人行走过程中的平衡控制、火箭发射中的垂直度控制和卫星飞行中的姿态控制等。所述方法的执行流程包括：步骤(1)：对于分数阶积分算子所述计算机设备令阶次λ在(01)范围内取m个均匀点，并记为λi，其中i＝1,2,…,m；并且，对于分数阶微分算子所述计算机设备令阶次λ在(-10)范围内取m个均匀点，记为λi，其中i＝1,2,…,m，分数阶积分算子和分数阶微分算子的符号都为其中积分的运用范围是(0,1)，微分的运用范围是(-1,0)。具体的，s表示拉普拉斯变换后的算子，例如设一个系统的输入函数为x(t)，输出函数为y(t)，则y(t)的拉氏变换y(s)与x(t)的拉氏变换x(s)的商:w(s)＝y(s)/x(s)称为这个系统的传递函数。，后面的斜体和正体是一样的。；例如：通过quanser板卡的多功能采集器，从外部系统采集输入变量给到了外部二级倒立摆。步骤(2)：离散传递函数将分数阶积分作为滤波器的时域离散传递函数,并求得单位脉冲输入δ(t)的情况下，所述滤波器输出的时域步骤(3)：根据所述滤波器的输入δ(t)和输出利用prony算法求出所述滤波器的时域离散传递函数；其中，γ()表示伽玛函数(gamma函数)，也叫欧拉第二积分，是阶乘函数在实数与复数上扩展的一类函数，伽玛函数作为阶乘的延拓，是定义在复数范围内的亚纯函数；步骤(4)：根据所述时域离散传递函数，构造分母系数矩阵a；步骤(5)：根据所述时域离散传递函数，构造分子系数矩阵b；步骤(6)：以离散点阶次x＝[λ0,λ1,…λm]为输入，以分母系数矩阵a的列向量转置为输出，采用两点差值法进行拟合计算，得到n+1个分母的系数，其中，n表示大于等于0的整数，a[n,1],a[n,2],…a[n,m]代表矩阵a第一列到最后一列的数值；步骤(7)：以离散点阶次x＝[λ0,λ1,…λm]为输入，以分子系数矩阵b的列向量转置为输出，采用两点差值法进行拟合计算，得到n+1个分子的系数；步骤(8)：由n+1个分母系数和n+1个分子系数，得到分数阶微积分的时域整数阶传递函数。所得到的分数阶微积分的时域整数阶传递函数表示为：由计算机设备中的程序，计算得到传递函数的完整表达式后，可以输入quanser板卡，计算机设备的显示模块(比如屏幕)显示quanser板卡在输入传递函数的运行状态，其中，步骤(1)到步骤(7)可以采用matlab平台编写，步骤(8)可以采用目前主流的程序语言编写。在本实施例中，所述步骤(3)，具体包括：根据所述滤波器的输入δ(t)和输出，利用prony算法求出所述滤波器的时域离散传递函数其中，z代表复数域，在所述时域离散传递函数中作为自变量。在本实施例中，所述步骤(4)，具体包括：包括：根据m个离散的时域离散传递函数，获取分母系数矩阵a：其中an,…,a0为n+1个m×1维列向量。在本实施例中，具体包括：所述步骤(5)，包括：根据m个离散的时域离散传递函数，获取分子系数矩阵b：其中bn,…,b0为n+1个m×1维列向量。在本实施例中，所述步骤(6)中，具体包括：得到n+1个分母的系数：a0＝(λ*a[0,p]-λ*a[0,q]+λp*a[0,q]-λq*a[0,p])/(λp-λq)a1＝(λ*a[1,p]-λ*a[1,q]+λp*a[1,q]-λq*a[1,p])/(λp-λq)...an＝(λ*a[n,p]-λ*a[n,q]+λp*a[n,q]-λq*a[n,p])/(λp-λq)先求出两点坐标(λp，a[n,p])与(λq，a[n,q])之间的直线方程，再代入所求点λ，得到分母系数，其中λ为所求点，其介于λp与λq之间，a[n,p]是λp所对应的输出，a[n,q]是λq所对应的输出；p，q两点是在x＝[λ0,λ1,…λm]中任意取值的x横坐标，以便得到相对应的y纵坐标，从而满足差值方程的条件。所述步骤(7)中，包括：得到n+1个分子的系数：b0＝(λ*b[0,p]-λ*b[0,q]+λp*b[0,q]-λq*b[0,p])/(λp-λq)b1＝(λ*b[1,p]-λ*b[1,q]+λp*b[1,q]-λq*b[1,p])/(λp-λq)...bn＝(λ*b[n,p]-λ*b[n,q]+λp*b[n,q]-λq*b[n,p])/(λp-λq)先求出两点坐标(λp，b[n,p])与(λq，b[n,q])之间的直线方程，再代入所求点λ，得到分子系数，其中λ为所求点，其介于λp与λq之间，b[n,p]是λp所对应的输出，b[n,q]是λq所对应的输出。所述步骤(8)中，所得到的分数阶微积分的时域整数阶传递函数表示为：在现有方案中，通常通过irid_fod(impulseresponseinvariantdiscretization_fractionalorderderivative)方法得到分数阶微积分离散传递函数，其方法是利用拉普拉斯逆变换得到分数阶微积分算子的单位脉冲响应函数，并取响应函数的n个离散点，用prony方法来拟合n个离散点，获得z离散传递函数如式：其中，α为分数阶微积分的阶次，n表示αi传递函数分子分母的阶次，a0…an为分母系数,bn…b0分子系数。可以看出现有技术所设计的积分仅仅能在知道α和n的情况下，得出传递函数，如取α＝-0.5，n＝5得传递函数如式：此方法无法做到实时获得传递函数。具体的计算逻辑如图2所示，包括：本发明公开了一种快速得到任意阶次分数阶微积分的时域传递函数实现方法。该方法首先在(01)或(-10)区间内取m个均匀阶次点，然后运用拉普拉斯逆变换得到m个离散阶次分数阶微积分的时域传递函数组；然后利用上述传递函数组的分子系数、分母系数利用最小二乘方法构造拟合多项式；最后对于任意可变阶次，可迅速通过多项式函数获得该阶次的对应分数阶微积分的时域传递函数。本发明所要解决的技术问题是提供一种快速得到任意阶次的离散s传递函数的方法。基于irid_fod分数阶数值方法的定阶次分析，基于此提出并设计了一种变阶次微积分α(t)的数值实现方法。并依据此对变阶次分数阶微积分数值计算进行时域二次拟合，方便了变阶次分数阶的数值方法实现。在本实施例中，提出一种变阶次分数阶微积分数值计算的时域二次拟合法。该方法通过时域拟合获得多个等间隔阶次时分数阶微积分算子的时域连续传递函数；然后构造时域传递函数同位置系数矩阵，并以系数矩阵的列向量进行多项式拟合，获得阶次为输入、系数为输出的拟合函数组；最后对于任意阶次的分数阶微积分算子，可快速解算并构建出其连续传递函数，从而实现变阶次分数阶微积分的时域数值计算。通过多项式计算的方法快速获得变阶分数阶微积分的时域传递函数,且具有很高的精度。在所述步骤(8)中，所得到的分数阶微积分的时域整数阶传递函数表示为：参照图3至图14所示的y相对于x的散点图。其中，离散的阶次x＝[α1,α2,…,αm]为横坐标,以为纵坐标。从图3-图14中可以看出，阶次在(-1,0)上时得到的z传递函数系数可以拟合成连续的曲线。通过图16和图17对比可以发现，本实施例所述的两点差值法得到的波形与irid_fod方法得到的波形误差不大，充分证明了两点差值法的可行性。相对于现有方案中，通常通过irid_fod方法得到分数阶微积分离散传递函数，本实施例提出的方法能够更快地运算出其传递函数，从表1-3中任意取的三个阶次的运算时间如表4(计算环境：cpu主频2ghz)：表1α＝-0.23时系数分析表表2α＝-0.57时系数分析表分母、分子系数iridfod方法本文方法绝对误差b50.007220.00742.49％b4-0.01808-0.01852.32％b30.015210.01562.56％b2-0.00448-0.00462.68％b10.000055010.00007260531.99％b00.00008030.0000807410.55％a51.00001.00000.00％a4-3.1947-3.19470.00％a33.80113.80110.00％a2-2.0387-2.03880.00％a10.46230.46240.02％a0-0.02993-0.0299-0.10％表3α＝-0.72时系数分析表相对于irid_fod方法，本文提出的方法能够更快地运算出其传递函数，上文中任意取的三个阶次的运算时间如下表(计算环境：cpu主频2ghz)：表4运算时间分析表由此可见，相对于irid_fod方法，本实施例提出的方法能够更快地运算出其传递函数；该方法采用脉冲响应不变原理，通过时域拟合获得多个等间隔阶次时分数阶微积分算子的时域离散传递函数；然后构造离散传递函数同位置系数矩阵，并以系数矩阵的列向量进行多项式拟合，获得阶次为输入、系数为输出的拟合函数组；最后对于任意阶次的分数阶微积分算子，可快速解算并构建出其离散传递函数，从而实现变阶次分数阶微积分的数值计算。仿真结果表明该方法具有较高的精度和快速性irid_fod。因此，本发明提供的方案实现了快速求解分数阶传递函数，提高了传递函数的计算效率。也间接提高了在工业控制时的精度以及效率。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于设备实施例而言，由于其基本相似于方法实施例，所以描述得比较简单，相关之处参见方法实施例的部分说明即可。以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本
技术领域：
的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。当前第1页12