CN107300862A - 一种变阶分数阶微积分频域分析方法 - Google Patents

一种变阶分数阶微积分频域分析方法 Download PDF

Info

Publication number
CN107300862A
CN107300862A CN201710495644.6A CN201710495644A CN107300862A CN 107300862 A CN107300862 A CN 107300862A CN 201710495644 A CN201710495644 A CN 201710495644A CN 107300862 A CN107300862 A CN 107300862A
Authority
CN
China
Prior art keywords
order
coefficients
fractional
frequency domain
denominator
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710495644.6A
Other languages
English (en)
Inventor
黄家才
舒奇
施昕昕
胡凯
张玎橙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Institute of Technology
Original Assignee
Nanjing Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Institute of Technology filed Critical Nanjing Institute of Technology
Priority to CN201710495644.6A priority Critical patent/CN107300862A/zh
Publication of CN107300862A publication Critical patent/CN107300862A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00—Systems involving the use of models or simulators of said systems
    • G05B17/02—Systems involving the use of models or simulators of said systems electric
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/20—Pc systems
    • G05B2219/23—Pc programming
    • G05B2219/23453—Pc simulates equipment and is connected to sequencer to test program
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/20—Pc systems
    • G05B2219/23—Pc programming
    • G05B2219/23456—Model machine for simulation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Complex Calculations (AREA)

Abstract

本发明实施例公开了一种变阶分数阶微积分频域分析方法，涉及电子信息技术领域，能够提高传递函数的计算效率。本发明包括：通过频域拟合获得多个等间隔阶次时分数阶微积分算子的频域连续传递函数；然后构造频域传递函数同位置系数矩阵，并以系数矩阵的列向量进行多项式拟合，获得阶次为输入、系数为输出的拟合函数组；最后对于任意阶次的分数阶微积分算子，可快速解算并构建出其连续传递函数。本发明适用于工程应用中的传递函数的实时获取。

Description

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

Claims (6)

1.一种变阶分数阶微积分频域分析方法，其特征在于，包括：
所述方法用于一种分析系统，所述分析系统用于为外部系统提供分数阶微积分计算，所述分析系统包括：计算机设备、利用PCI接口连接所述计算机设备的quanser板卡(QPIDe)、所述quanser板卡与模拟传感器适配器电源和信号路由器组成的通信模块连接二级倒立摆；
所述方法包括：
步骤(1)：对于分数阶积分算子所述计算机设备令阶次λ在(0 1)范围内取M个均匀点，并记为λi，其中i＝1,2,…,M；并且，对于分数阶微分算子所述计算机设备令阶次λ在(-1 0)范围内取M个均匀点，记为λi，其中i＝1,2,…,M，分数阶积分算子和分数阶微分算子的符号都为其中积分的运用范围是(0,1)，微分的运用范围是(-1,0)，s表示拉普拉斯变换后的算子；
步骤(2)：所述计算机设备对于各阶次的λi，得到分数阶积分的频域整数阶近似传递函数；
步骤(3)：根据所述频域整数阶近似传递函数，构造分母系数矩阵A；
步骤(4)：根据所述频域整数阶近似传递函数，构造分子系数矩阵B；
步骤(5)：以离散点阶次X＝[λ0,λ1,…λM]为输入，以分母系数矩阵A的列向量转置为输出，由两点(λp,a[n,p])与(λq,a[n,q])构造插值方程，即可求得λ点对应的第n个分母系数，其中，对于时变阶次λ，假设λq≤λ≤λp，其中，p，q两点是在X＝[λ0,λ1,…λM]中任意取值的x横坐标，以便得到相对应的y纵坐标，从而满足差值方程的条件，a[n,p]是λp所对应的第n个分母系数输出，a[n,q]是λq所对应的第n个分母系数输出，λ点对应的第n个分母系数表示为：
步骤(6)：以离散点阶次X＝[λ0,λ1,…λM]为输入，以分子系数矩阵B的列向量转置为输出，由两点(λp,b[n,p])与(λq,b[n,q])构造插值方程，求得λ点对应的第n个分子系数，其中，λ点对应的第n个分子系数表示为：
步骤(7)：由所求点所在间断点区间，得到分数阶微积分的分母系数a0，a1，a2……a2N+1和分子系数b0，b1，b2……b2N+1,并通过所述分子系数和所述分母系数，构建分数阶积分的频域整数阶传递函数。
2.根据权利要求1所述的方法，其特征在于，所述步骤(2)，包括：
所述计算机设备对于任意可变阶次λ，，根据其所在区间可以快速得到分数阶微积分的分母系数a0，a1，a2……a2N+1和分子系数b0，b1，b2……b2N+1,得到分数阶积分的频域整数阶近似传递函数
其中，拉普拉斯变换用f(t)表示实变量t的一个函数，F(s)表示它的拉普拉斯变换，其中s表示是一个复变量，在函数中表示一个自变量。
3.根据权利要求2所述的方法，其特征在于，所述步骤(3)，包括：
根据M个离散的频域整数阶近似传递函数，获取分母系数矩阵A：
其中，A2N+1,…,A0为2N+2个M×1维列向量。
4.根据权利要求2所述的方法，其特征在于，所述步骤(4)，包括：
根据M个离散的频域整数阶近似传递函数，获取分子系数矩阵B：
其中，B2N+1,…,B0为2N+2个M×1维列向量。
5.根据权利要求3或4所述的方法，其特征在于，在所述步骤(5)中，包括：
先求出两点坐标(λp，a[n,p])与(λq，a[n,q])之间的直线方程，再代入所求点λ，得到分母系数，其中，λ为所求点，其介于λp与λq之间，a[n,p]是λp所对应的输出，a[n,q]是λq所对应的输出；
所述步骤(5)中的2N+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)...
a2N+1＝(λ*a[2N+1,p]-λ*a[2N+1,q]+λp*a[2N+1,q]-λq*a[2N+1,p])/(λp-λq)。
6.根据权利要求3或4所述的方法，其特征在于，在所述步骤(6)中，包括：
先求出两点坐标(λp，b[n,p])与(λq，b[n,q])之间的直线方程，再代入所求点λ，得到分子系数，其中，λ为所求点，其介于λp与λq之间，b[n,p]是λp所对应的输出，b[n,q]是λq所对应的输出；
所述步骤(6)中的2N+2个的分子系数表示为：
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)...
b2N+1＝(λ*b[2N+1,p]-λ*b[2N+1,q]+λp*b[2N+1q]-λq*b[2N+1,p])/(λp-λq)。
CN201710495644.6A 2017-06-26 2017-06-26 一种变阶分数阶微积分频域分析方法 Pending CN107300862A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710495644.6A CN107300862A (zh) 2017-06-26 2017-06-26 一种变阶分数阶微积分频域分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710495644.6A CN107300862A (zh) 2017-06-26 2017-06-26 一种变阶分数阶微积分频域分析方法

Publications (1)

Publication Number Publication Date
CN107300862A true CN107300862A (zh) 2017-10-27

Family

ID=60135542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710495644.6A Pending CN107300862A (zh) 2017-06-26 2017-06-26 一种变阶分数阶微积分频域分析方法

Country Status (1)

Country Link
CN (1) CN107300862A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113791353A (zh) * 2021-10-11 2021-12-14 河北工业大学 基于分数阶传递函数的锂电池电压模型构建方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113791353A (zh) * 2021-10-11 2021-12-14 河北工业大学 基于分数阶传递函数的锂电池电压模型构建方法

Similar Documents

Publication Publication Date Title
Chen et al. An enhanced hybrid MobileNet
Shah et al. Review of fractional PID controller
Matouk et al. Dynamical analysis, stabilization and discretization of a chaotic fractional-order GLV model
Zhu et al. Stability analysis for stochastic neural networks of neutral type with both Markovian jump parameters and mixed time delays
Benner et al. Generalised tangential interpolation for model reduction of discrete-time MIMO bilinear systems
Li et al. State-feedback ℋ∞ control for stochastic time-delay nonlinear systems with state and disturbance-dependent noise
CN116151084A (zh) 基于结构网格的模拟方法、装置、终端设备及存储介质
Liu et al. Bifurcation of a predator–prey model with disease in the prey
Ding et al. Slimyolov4: lightweight object detector based on yolov4
CN107220214A (zh) 一种基于多项式拟合的变阶分数阶微积分频域分析方法
CN107300862A (zh) 一种变阶分数阶微积分频域分析方法
Nazemi et al. Solving a class of geometric programming problems by an efficient dynamic model
CN102929842A (zh) 基于fpaa的可重构矢量-矩阵乘法器设计方法
CN101763086A (zh) 一种高效的工业过程动态优化系统及方法
CN111353226B (zh) 多维控制系统模型降阶方法及装置
CN116400598B (zh) 一种惯性cg神经网络的同步控制方法
CN107315347A (zh) 一种变阶分数阶微积分时域分析方法
CN107154802A (zh) 一种校正测量数据的方法及装置
CN117031946A (zh) 一种基于极限学习机的时变系统迭代学习控制方法
Kalangadan et al. PI, PID controller design for interval systems using frequency response model matching technique
Handkiewicz et al. Design automation of a lossless multiport network and its application to image filtering
Thakoor et al. Fast valuation of CEV American options
Angkeaw et al. Design and Implementation of FPAA based LQR Controller for Magnetic Levitation Control System
CN107300854A (zh) 一种基于多项式拟合的变阶分数阶微积分时域分析方法
CN115422860A (zh) 应用于仿真电路的支路电流获取方法、系统、设备和计算机可读存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171027