Movatterモバイル変換


[0]ホーム

URL:


CN114741657A - Fixed time algebraic parameter estimation algorithm of Chua circuit system - Google Patents

Fixed time algebraic parameter estimation algorithm of Chua circuit system
Download PDF

Info

Publication number
CN114741657A
CN114741657ACN202210435371.7ACN202210435371ACN114741657ACN 114741657 ACN114741657 ACN 114741657ACN 202210435371 ACN202210435371 ACN 202210435371ACN 114741657 ACN114741657 ACN 114741657A
Authority
CN
China
Prior art keywords
formula
fixed
circuit system
parameter estimation
estimation algorithm
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.)
Granted
Application number
CN202210435371.7A
Other languages
Chinese (zh)
Other versions
CN114741657B (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHUfiledCriticalHohai University HHU
Priority to CN202210435371.7ApriorityCriticalpatent/CN114741657B/en
Publication of CN114741657ApublicationCriticalpatent/CN114741657A/en
Application grantedgrantedCritical
Publication of CN114741657BpublicationCriticalpatent/CN114741657B/en
Activelegal-statusCriticalCurrent
Anticipated expirationlegal-statusCritical

Links

Images

Classifications

Landscapes

Abstract

The invention discloses a fixed time algebraic parameter estimation algorithm of a Chua circuit system, which is used for predicting and outputting a fixed time algebraic parameter based on a constructed fixed time algebraic parameter estimation algorithm; constructing a fixed time algebraic parameter estimation algorithm: establishing a nonlinear mathematical model according to the Chua circuit system; determining a kernel function by using a Volterra integral operator; determining the stability of the fixed time; and constructing and obtaining a fixed time algebraic parameter estimation algorithm according to the kernel function and the nonlinear mathematical model. The invention can realize the accurate estimation of the parameters in the fixed time which does not depend on the initial value of the system under the condition that only the output signal is known; by skillfully selecting the kernel function in the Volterra integral operator, the influence of the initial value of the system can be effectively eliminated, and the calculation of the output derivative of the Chua circuit system is avoided; different from the traditional adaptive parameter estimation algorithm, the algebraic parameter estimation algorithm provided by the invention can realize the rapid and accurate estimation of parameters.

Description

Translated fromChinese
Chua电路系统的固定时间代数参数估计算法A Fixed-Time Algebraic Parameter Estimation Algorithm for Chua Circuit Systems

技术领域technical field

本发明涉及Chua电路系统的固定时间代数参数估计算法,属于Chua电路系统参数估计技术领域。The invention relates to a fixed-time algebraic parameter estimation algorithm of a Chua circuit system, and belongs to the technical field of parameter estimation of a Chua circuit system.

背景技术Background technique

混沌现象作为微分动力系统中十分重要的非线性动力学现象,因其在安全通信、激光系统、电子化学和神经生理学等领域的应用,受到学者们的广泛关注。自混沌学的发展逐渐步入正轨以来,人们一直尝试用肉眼可见的方法去观察混沌这一令人捉摸不透的动力学现象。由此,电学中非线性电路的研究逐渐进入人们的眼球,同时成为揭开分支与混沌机理过程中的重要研究课题。关于非线性电路中混沌现象的研究有三十年以上的历史,在研究者们的不懈努力下,许多以混沌机理研究为目的的电路被构造出来。同时,研究者们还对一部分实用电路中产生的混沌现象进行了密切研究。Chaos, as a very important nonlinear dynamical phenomenon in differential dynamical systems, has received extensive attention from scholars due to its applications in secure communications, laser systems, electronic chemistry and neurophysiology. Since the development of chaos theory has been on the right track, people have been trying to observe the elusive dynamic phenomenon of chaos with the naked eye. As a result, the study of nonlinear circuits in electricity has gradually entered people's attention, and has become an important research topic in the process of uncovering the mechanism of branching and chaos. The research on the chaotic phenomenon in nonlinear circuits has a history of more than 30 years. With the unremitting efforts of researchers, many circuits have been constructed for the purpose of studying the chaotic mechanism. At the same time, researchers have also closely studied the chaotic phenomenon generated in some practical circuits.

1983年,加利福尼亚大学的Leon Chua教授设计出了第一个可以通过实验模拟混沌现象发生的Chua电路。这是一个具有里程碑意义的设计,它开辟了混沌行为应用于现实世界的新道路。Chua电路是一个可以产生双涡卷混沌的三阶自治电路,通过赋予一个分段线性负电阻不同形式的参数组合,从而产生尤为丰富的分支与混沌行为。Chua电路是能引起混沌现象出现的自治电路中构造最为简单的一个,任何在三阶自治系统中产生的混沌现象,都可以通过Chua电路系统模拟展现出来。Shilnikov定理中对Chua系统产生的混沌形态做出了严格证明,基于此,人们得以通过电路这一形式来对分支和混沌的各种机理进行观察和研究。因此对于Chua电路的研究引起了很多研究者的兴趣,Chua系统也逐渐成为了混沌现象研究中的标准模型。In 1983, Professor Leon Chua of the University of California designed the first Chua circuit that can simulate the occurrence of chaotic phenomena experimentally. This is a landmark design that opens up new avenues for chaotic behavior to be applied to the real world. Chua circuit is a third-order autonomous circuit that can generate double scroll chaos. By giving a piecewise linear negative resistance different forms of parameter combinations, it can produce particularly rich branching and chaotic behaviors. The Chua circuit is the simplest one of the autonomous circuits that can cause chaotic phenomena. Any chaotic phenomenon generated in the third-order autonomous system can be simulated by the Chua circuit system. Shilnikov's theorem has made a strict proof of the chaotic form generated by Chua system. Based on this, people can observe and study various mechanisms of branching and chaos in the form of circuits. Therefore, the study of Chua circuit has aroused the interest of many researchers, and the Chua system has gradually become the standard model in the study of chaotic phenomena.

最初的关于Chua电路系统的研究需要假设系统的参数是精确已知的。考虑到在实际的电路设计中,由于受到元器件老化、不确定干扰以及测量成本等各种内部外部因素的影响,电路中的某些参数往往很难精确或直接获取。因此,如何解决具有未知参数的Chua电路系统的参数估计问题就变得尤为重要,当系统仅有部分状态变量可用时,问题将变得更加困难。为了解决这一问题,学者们提出了各种各样不同的方法,如基于延迟嵌入的方法和基于控制理论的方法等。其中基于控制理论的方法可以通过观测器设计和系统辨识等方法从隐藏变量中估计未知参数,从而得到越来越广泛的应用。Initial studies of Chua's circuit system assumed that the parameters of the system were precisely known. Considering that in the actual circuit design, due to the influence of various internal and external factors such as component aging, uncertain interference and measurement cost, some parameters in the circuit are often difficult to obtain accurately or directly. Therefore, how to solve the parameter estimation problem of the Chua circuit system with unknown parameters becomes particularly important, and the problem becomes more difficult when only part of the state variables of the system are available. In order to solve this problem, scholars have proposed a variety of different methods, such as methods based on delayed embedding and methods based on control theory. Among them, the methods based on control theory can estimate unknown parameters from hidden variables through methods such as observer design and system identification, which are more and more widely used.

估计参数的收敛速度是在参数估计的过程中非常关键的一个性能指标,前文中所涉及的文献大多仅能实现参数的渐近估计。渐近估计仅能保证参数的估计误差在时间趋向于无穷时趋向于零,因此无法得到更快的估计速度。从时间优化的角度来看,使估计误差有限时间收敛的估计方法才是时间最优的估计算法。此外,有限时间估计算法通常带有分数幂项,使得有限时间参数估计和传统的渐近估计相比,往往具有更好的鲁棒性和抗扰动性。正是由于有限时间参数估计的诸多优点,关于混沌系统的有限时间参数估计研究受到越来越多研究者的关注,并取得了许多积极的成果。The convergence speed of the estimated parameters is a very critical performance index in the process of parameter estimation, and most of the literatures involved in the above can only realize the asymptotic estimation of the parameters. Asymptotic estimation can only guarantee that the estimation error of the parameter tends to zero when the time tends to infinity, so it cannot obtain a faster estimation speed. From the point of view of time optimization, the estimation method that makes the estimation error converge in limited time is the time-optimized estimation algorithm. In addition, the finite-time estimation algorithm usually has fractional power terms, which makes the finite-time parameter estimation have better robustness and anti-disturbance than the traditional asymptotic estimation. Because of the many advantages of finite-time parameter estimation, the research on finite-time parameter estimation of chaotic systems has been paid more and more attention by researchers, and many positive results have been achieved.

不管是传统的渐近参数估计算法还是近年来被广泛研究的有限时间估计算法,其收敛时间都是依赖于参数的初始估计误差的,并随初始估计误差的增大而增大。为了克服这一缺陷,有关文献提出了固定时间稳定性的概念。固定时间稳定性在保留了有限时间稳定收敛速度快、收敛精度高等诸多优点的同时,还具有收敛时间不依赖初始值的优良特性,因此成为国内外的研究热点。近年来,混沌系统的固定时间同步和固定时间控制等问题被大量研究,但关于混沌系统的固定时间参数估计问题的研究却鲜有报道。Whether it is the traditional asymptotic parameter estimation algorithm or the finite-time estimation algorithm that has been widely studied in recent years, the convergence time depends on the initial estimation error of the parameters, and increases with the increase of the initial estimation error. In order to overcome this defect, the related literature proposes the concept of fixed time stability. Fixed-time stability retains the advantages of fast convergence speed and high convergence accuracy of finite-time stability, and also has the excellent characteristic that the convergence time does not depend on the initial value, so it has become a research hotspot at home and abroad. In recent years, fixed-time synchronization and fixed-time control of chaotic systems have been studied a lot, but the research on fixed-time parameter estimation of chaotic systems is rarely reported.

发明内容SUMMARY OF THE INVENTION

本发明所要解决的技术问题是克服现有技术的缺陷,提供基于Volterra积分算子的Chua电路系统的固定时间代数参数估计算法,该算法能够做到在仅有输出信号已知的情况下,在不依赖于系统初始值的固定时间内实现对参数的精确估计;通过在Volterra积分算子中巧妙的选取核函数,本发明能够有效消除系统初始值的影响,同时避免了对Chua电路系统输出导数的计算;此外,不同于传统的自适应参数估计算法,本发明提出的代数参数估计算法能够实现参数的快速、准确估计。The technical problem to be solved by the present invention is to overcome the defects of the prior art, and to provide a fixed-time algebraic parameter estimation algorithm of the Chua circuit system based on the Volterra integral operator, which can achieve a Accurate estimation of parameters is realized in a fixed time independent of the initial value of the system; by cleverly selecting the kernel function in the Volterra integral operator, the present invention can effectively eliminate the influence of the initial value of the system, and avoid the output derivative of the Chua circuit system. In addition, different from the traditional self-adaptive parameter estimation algorithm, the algebraic parameter estimation algorithm proposed by the present invention can realize fast and accurate estimation of parameters.

为达到上述目的,本发明提供Chua电路系统的固定时间代数参数估计算法,包括:In order to achieve the above object, the present invention provides a fixed-time algebraic parameter estimation algorithm of the Chua circuit system, including:

基于构建的固定时间代数参数估计算法,预测输出固定时间代数参数;Based on the constructed fixed-time algebraic parameter estimation algorithm, predict and output fixed-time algebraic parameters;

构建固定时间代数参数估计算法,通过以下步骤实现:Construct a fixed-time algebraic parameter estimation algorithm, implemented through the following steps:

根据Chua电路系统,建立非线性数学模型;According to the Chua circuit system, establish a nonlinear mathematical model;

利用Volterra积分算子,确定核函数K(t,τ);Use the Volterra integral operator to determine the kernel function K(t,τ);

确定固定时间稳定性;determine fixed time stability;

根据核函数K(t,τ)和非线性数学模型,构建获得固定时间代数参数估计算法。According to the kernel function K(t,τ) and the nonlinear mathematical model, a fixed-time algebraic parameter estimation algorithm is constructed.

优先地,根据Chua电路系统,构建简化非线性模型,并建立无量纲非线性数学模型,包括:Chua电路系统的简化非线性模型的表达式为:Preferentially, according to the Chua circuit system, a simplified nonlinear model is constructed, and a dimensionless nonlinear mathematical model is established, including: the expression of the simplified nonlinear model of the Chua circuit system is:

Figure BDA0003612689220000031
Figure BDA0003612689220000031

式中,R为线性电阻,

Figure BDA0003612689220000032
是电容器C1的电压,
Figure BDA0003612689220000033
是电容器C2的电压,il是通过电感L的电流;
Figure BDA0003612689220000034
是关于电容器电压的函数,
Figure BDA0003612689220000035
表达式为:where R is the linear resistance,
Figure BDA0003612689220000032
is the voltage across capacitorC1 ,
Figure BDA0003612689220000033
is the voltage across capacitor C2 andi lis the current through inductor L;
Figure BDA0003612689220000034
is a function of capacitor voltage,
Figure BDA0003612689220000035
The expression is:

Figure BDA0003612689220000036
Figure BDA0003612689220000036

公式(2)中m0,m1和Bp为二极管的三个固定常数;In formula (2), m0 , m1 and Bp are three fixed constants of the diode;

利用无量纲形式,构建非线性数学模型;Use dimensionless form to construct nonlinear mathematical model;

非线性数学模型的表达式为:The expression of the nonlinear mathematical model is:

Figure BDA0003612689220000037
Figure BDA0003612689220000037

式中,f(x1)=ax1+b(|x1+1|-|x1-1|) (4),In the formula, f(x1 )=ax1 +b(|x1 +1|-|x1 -1|) (4),

Figure BDA0003612689220000038
Figure BDA0003612689220000038

式中,y1和y2表示公式(3)的输出;m0、m1和Bp为Chua二极管的三个固定常数,

Figure BDA0003612689220000039
是电容C1的电压,
Figure BDA00036126892200000310
是电容C2的电压,il是通过电感L的电流。In the formula, y1 and y2 represent the output of formula (3); m0 , m1 and Bp are the three fixed constants of the Chua diode,
Figure BDA0003612689220000039
is the voltage across capacitorC1 ,
Figure BDA00036126892200000310
is the voltage across capacitor C2 andi lis the current through inductor L.

优先地,利用Volterra积分算子,确定核函数K(t,τ),包括:Preferentially, using the Volterra integral operator, determine the kernel function K(t,τ), including:

定义

Figure BDA00036126892200000311
为R≥0到R的希尔伯特空间内的局部平方可积函数,R表示实数,R≥0表示非负实数,局部平方可积函数f(t)经过Volterra积分算子进行映射定义:definition
Figure BDA00036126892200000311
is a local square integrable function in the Hilbert space from R≥ 0 to R, R represents a real number, R≥ 0 represents a non-negative real number, and the local square integrable function f(t) is defined by the Volterra integral operator mapping:

Figure BDA00036126892200000312
Figure BDA00036126892200000312

式中,K(·,·):R×R→R希尔伯特-施密特核函数;In the formula, K(·,·): R×R→R Hilbert-Schmidt kernel function;

基于莱布尼兹微分规则,[VKf](t)由如下微分方程产生:Based on Leibniz's differential rule, [VK f](t) is generated by the following differential equation:

Figure BDA0003612689220000041
Figure BDA0003612689220000041

式中,K(t,t)为核函数,ξ0为预设的初始值常数;In the formula, K(t, t) is the kernel function, and ξ0 is the preset initial value constant;

定义f(i)(t),i∈Z+为f(t)的i阶导数,Z+为正整数,则[VKf(i)](t)展开为:Define f(i) (t), i∈Z+ is the i-order derivative of f(t), and Z+ is a positive integer, then [VK f(i) ](t) is expanded as:

Figure BDA0003612689220000042
Figure BDA0003612689220000042

式中,K(i)(·,·)表示K(·,·)关于第二个变量的i阶导数;In the formula, K(i) (·,·) represents the i-order derivative of K(·,·) with respect to the second variable;

核函数K(t,τ)的选取决定了Volterra积分算子的映射[VKf](t)效果,选取的核函数为:

Figure BDA0003612689220000043
The selection of the kernel function K(t,τ) determines the mapping [VK f](t) effect of the Volterra integral operator, and the selected kernel function is:
Figure BDA0003612689220000043

公式(9)中,ωh,ω>0均为预设的参数;In formula (9), ωh , ω>0 are all preset parameters;

公式(9)展开得到:Equation (9) is expanded to get:

Figure BDA0003612689220000044
Figure BDA0003612689220000044

其中fp(τ),p=0,1,2满足:Where fp (τ), p = 0, 1, 2 satisfy:

Figure BDA0003612689220000045
Figure BDA0003612689220000045

式中,q={0,1,2};In the formula, q={0,1,2};

基于公式(11),计算获得

Figure BDA0003612689220000046
Based on formula (11), the calculation is obtained
Figure BDA0003612689220000046

Figure BDA0003612689220000051
Figure BDA0003612689220000051

基于公式(10)和公式(12),得到核函数关于的导数:Based on formula (10) and formula (12), the derivative of the kernel function with respect to:

Figure BDA0003612689220000052
Figure BDA0003612689220000052

得到对于

Figure BDA0003612689220000053
有:get for
Figure BDA0003612689220000053
Have:

K(i)(t,t)=0,K(i)(t,0)=0(14)。K(i) (t,t)=0, K(i) (t,0)=0(14).

优先地,确定固定时间稳定性,包括:Preferentially, determine fixed-time stability, including:

Chua电路系统为:The Chua circuit system is:

Figure BDA0003612689220000056
Figure BDA0003612689220000056

式中,f(0)=0,x(t)∈Rn,f(x(t)):Rn→Rn为定义域n维空间Rn到值域n维空间Rn中的一个连续函数,x0=x(t0)表示Chua电路系统的初始值,t0表示初始时刻;In the formula, f(0)=0, x(t)∈Rn , f(x(t)): Rn →Rn is one of the definition domain n-dimensional space Rn to the value domain n-dimensional space Rn Continuous function, x0 =x(t0 ) represents the initial value of the Chua circuit system, and t0 represents the initial moment;

将公式(15)的解记为x(t,x0);Denote the solution of formula (15) as x(t,x0 );

若存在一个不依赖于Chua电路系统初始值x0的固定时间T,使得对任意初始状态x0,有

Figure BDA0003612689220000054
成立,则判定Chua电路系统的平衡点x=0是全局固定时间稳定的。If there is a fixed time T that does not depend on the initial value x0 of the Chua circuit system, such that for any initial state x0 , we have
Figure BDA0003612689220000054
If it is established, it is determined that the balance point x=0 of the Chua circuit system is stable in the global fixed time.

优先地,在构建固定时间代数参数估计算法之前,对Chua电路系统做如下处理:Preferentially, before constructing a fixed-time algebraic parameter estimation algorithm, the Chua circuit system is processed as follows:

将公式(3)变换为关于Chua电路系统输出及其导数形式:Transform equation (3) into the form of the output of the Chua circuit system and its derivative:

Figure BDA0003612689220000055
Figure BDA0003612689220000055

其中,y3=|y1+1|-|y1-1|,where y3 =|y1 +1|-|y1 -1|,

参数θ1,θ2和θ3满足:The parameters θ1 , θ2 and θ3 satisfy:

θ1=-β(1+a),θ2=β,θ3=-bβ (17);θ1 =-β(1+a), θ2 =β, θ3 =-bβ (17);

将对β,a,b,γ的估计转变为对θ123,γ的估计;Convert the estimates of β, a, b, γ into estimates of θ1 , θ2 , θ3 , γ;

对公式(16)两边同时进行积分算子VK运算:Perform the integral operator VK operation on both sides of formula (16) at the same time:

Figure BDA0003612689220000061
Figure BDA0003612689220000061

式中,y1(t),y2(t)为公式(3)的输出;In the formula, y1 (t), y2 (t) is the output of formula (3);

将公式(14)带入公式(18),得到:Substituting equation (14) into equation (18), we get:

Figure BDA0003612689220000062
Figure BDA0003612689220000062

将公式(19)带入公式(18),得到:Substituting equation (19) into equation (18), we get:

Figure BDA0003612689220000063
Figure BDA0003612689220000063

基于公式(20),设计参数估计器;Based on formula (20), design a parameter estimator;

Figure BDA0003612689220000064
构造微分方程:right
Figure BDA0003612689220000064
Construct the differential equation:

Figure BDA0003612689220000065
Figure BDA0003612689220000065

其中,

Figure BDA0003612689220000066
由公式(13)获得;in,
Figure BDA0003612689220000066
Obtained by formula (13);

若核函数K(t,τ)满足公式(9),则

Figure BDA0003612689220000067
Figure BDA0003612689220000068
成立;If the kernel function K(t,τ) satisfies formula (9), then
Figure BDA0003612689220000067
right
Figure BDA0003612689220000068
established;

由微分方程得到

Figure BDA0003612689220000069
因此公式(20)改写成:obtained from the differential equation
Figure BDA0003612689220000069
Therefore formula (20) is rewritten as:

Figure BDA00036126892200000610
Figure BDA00036126892200000610

其中,in,

Figure BDA0003612689220000071
Figure BDA0003612689220000071

其中,v1(t)和v2(t)满足持续可激励假设。Among them, v1 (t) and v2 (t) satisfy the continuous stimuli hypothesis.

优先地,在构建固定时间代数参数估计算法之前,设定如下持续可激励假设:Preferentially, before building a fixed-time algebraic parameter estimation algorithm, set the following persistent motivating assumptions:

假设1:函数v1(t)和v2(t)满足持续可激励条件,存在常数r>0,T0>0使得如下不等式对

Figure BDA0003612689220000072
成立:Assumption 1: The functions v1 (t) and v2 (t) satisfy the continuous excitation condition, and there are constants r > 0 and T0 > 0 such that the following inequalities
Figure BDA0003612689220000072
Established:

Figure BDA0003612689220000073
Figure BDA0003612689220000073

式中,I3×3为3×3的单位矩阵;In the formula, I3×3 is a 3×3 unit matrix;

对公式(23)分别左乘v1T(t)和v2(t),得到:Left-multiplying v1T (t) and v2 (t) to formula (23), respectively, yields:

Figure BDA0003612689220000074
Figure BDA0003612689220000074

其中in

Figure BDA0003612689220000075
Figure BDA0003612689220000075

对于公式(26)两边同时进行积分算子

Figure BDA0003612689220000076
运算,其中Kg=e-g(t-τ),g>0,得到:For both sides of equation (26), the integral operator is performed simultaneously
Figure BDA0003612689220000076
Operation, where Kg =e-g(t-τ) , g>0, we get:

Figure BDA0003612689220000077
Figure BDA0003612689220000077

其中,

Figure BDA0003612689220000078
in,
Figure BDA0003612689220000078

Chua电路系统产生变量s1,f(t)、s2,f(t)、v1,f(t)和v2,f(t):The Chua circuitry produces variables s1,f (t), s2,f (t), v1,f (t) and v2,f (t):

Figure BDA0003612689220000079
Figure BDA0003612689220000079

优先地,根据核函数K(t,τ)和非线性数学模型,构建获得固定时间代数参数估计算法,包括:固定时间代数参数估计算法的表达式为:Preferentially, according to the kernel function K(t,τ) and the nonlinear mathematical model, a fixed-time algebraic parameter estimation algorithm is constructed and obtained, including: the expression of the fixed-time algebraic parameter estimation algorithm is:

Figure BDA0003612689220000081
Figure BDA0003612689220000081

参数ε,g满足:

Figure BDA0003612689220000082
Parameters ε, g satisfy:
Figure BDA0003612689220000082

其中,

Figure BDA0003612689220000083
表示参数
Figure BDA0003612689220000084
的上一个值,
Figure BDA0003612689220000085
表示参数
Figure BDA0003612689220000086
的上一个值,min{eig(v1,f(t))}和min{eig(v2,f(t))}分别表示v1,f(t)和v2,f(t)特征值最小值。in,
Figure BDA0003612689220000083
Indicates parameters
Figure BDA0003612689220000084
the previous value of ,
Figure BDA0003612689220000085
Indicates parameters
Figure BDA0003612689220000086
The previous value of , min{eig(v1,f (t))} and min{eig(v2,f (t))} denote v1,f (t) and v2,f (t) features respectively value minimum.

一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现上述任一项所述方法的步骤。An electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of any one of the methods described above when the processor executes the program.

一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现上述任一项所述方法的步骤。A computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, implements the steps of any one of the above-mentioned methods.

本发明所达到的有益效果:Beneficial effects achieved by the present invention:

1、首先,和已有的基于渐近收敛和有限时间收敛的参数估计策略相比较,本发明所提出的固定时间参数估计算法在保证了收敛速度快和收敛精度高的同时,还具有收敛时间不依赖于初始估计误差的特点;1. First, compared with the existing parameter estimation strategies based on asymptotic convergence and finite-time convergence, the fixed-time parameter estimation algorithm proposed in the present invention not only ensures fast convergence speed and high convergence accuracy, but also has a convergence time. Does not depend on the characteristics of the initial estimation error;

2、此外,本发明所使用的Volterra积分算子,能够有效消除系统初始值的影响,同时避免对系统输出导数的计算,使得所提算法只需用到系统的输出信息,因此能够有效减少传感器的使用,从而有效降低算法的实现成本;2. In addition, the Volterra integral operator used in the present invention can effectively eliminate the influence of the initial value of the system, and at the same time avoid the calculation of the output derivative of the system, so that the proposed algorithm only needs to use the output information of the system, so it can effectively reduce the number of sensors. It can effectively reduce the implementation cost of the algorithm;

3、最后,在现有的文献中,为了得到足够多的方程来消去未知变量的影响,需要选取几个不同的核函数,而本发明所提算法除了能实现固定时间参数估计外,只需要选取一个核函数,从而使得算法复杂度大大降低。3. Finally, in the existing literature, in order to obtain enough equations to eliminate the influence of unknown variables, it is necessary to select several different kernel functions. In addition to realizing fixed-time parameter estimation, the algorithm proposed in the present invention only needs to A kernel function is selected, which greatly reduces the complexity of the algorithm.

附图说明Description of drawings

图1为本发明的Chua电路系统的原理图;1 is a schematic diagram of the Chua circuit system of the present invention;

图2为本发明的Chua电路系统的混沌吸引子三维图;2 is a three-dimensional diagram of a chaotic attractor of the Chua circuit system of the present invention;

图3为本发明的Chua电路系统在x1Ox2平面上投影的混沌吸引子图;Fig. 3 is the chaotic attractor diagram of the Chua circuit system of the present invention projected on the x1 Ox2 plane;

图4为本发明的Chua电路系统在x2Ox3平面上投影的混沌吸引子图;Fig. 4 is the chaotic attractor diagram of the Chua circuit system of the present invention projected on the x2 Ox3 plane;

图5为本发明的Chua电路系统在x1Ox3平面上投影的混沌吸引子图;Fig. 5 is the chaotic attractor diagram of the Chua circuit system of the present invention projected on the x1 Ox3 plane;

图6为本发明在自适应初始值

Figure BDA0003612689220000087
下β真实值和估计值曲线图;Fig. 6 is the adaptive initial value of the present invention
Figure BDA0003612689220000087
Lower β true value and estimated value curve;

图7为本发明在自适应初始值

Figure BDA0003612689220000091
下γ真实值和估计值曲线图;Fig. 7 is the adaptive initial value of the present invention
Figure BDA0003612689220000091
Lower γ true value and estimated value curve;

图8为本发明在自适应初始值

Figure BDA0003612689220000092
下a真实值和估计值曲线图;Fig. 8 is the adaptive initial value of the present invention
Figure BDA0003612689220000092
Below a graph of true and estimated values;

图9为本发明在自适应初始值

Figure BDA0003612689220000093
下b真实值和估计值曲线图;Fig. 9 is the adaptive initial value of the present invention
Figure BDA0003612689220000093
Lower b true value and estimated value curve;

图10为本发明在自适应初始值

Figure BDA0003612689220000094
下β真实值和估计值曲线图;Fig. 10 shows the adaptive initial value of the present invention
Figure BDA0003612689220000094
Lower β true value and estimated value curve;

图11为本发明在自适应初始值

Figure BDA0003612689220000095
下γ真实值和估计值曲线图;Fig. 11 is the adaptive initial value of the present invention
Figure BDA0003612689220000095
Lower γ true value and estimated value curve;

图12为本发明在自适应初始值

Figure BDA0003612689220000096
下a真实值和估计值曲线图;Fig. 12 is the adaptive initial value of the present invention
Figure BDA0003612689220000096
Below a graph of true and estimated values;

图13为本发明在自适应初始值

Figure BDA0003612689220000097
下n真实值和估计值曲线图;Fig. 13 shows the adaptive initial value of the present invention
Figure BDA0003612689220000097
Lower n real value and estimated value curve;

图14为本发明的流程图。Figure 14 is a flow chart of the present invention.

具体实施方式Detailed ways

以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。The following examples are only used to illustrate the technical solutions of the present invention more clearly, and cannot be used to limit the protection scope of the present invention.

实施例一Example 1

Chua电路系统的固定时间代数参数估计算法,包括以下步骤:A fixed-time algebraic parameter estimation algorithm for the Chua circuit system, including the following steps:

步骤1、根据Chua电路系统建立非线性数学模型;Step 1. Establish a nonlinear mathematical model according to the Chua circuit system;

步骤2、利用Volterra积分算子,确定核函数K(t,τ);Step 2. Use the Volterra integral operator to determine the kernel function K(t, τ);

步骤3、根据核函数K(t,τ)和非线性数学模型,确定固定时间稳定性,构建固定时间代数参数估计算法;Step 3, according to the kernel function K(t, τ) and the nonlinear mathematical model, determine the fixed-time stability, and construct a fixed-time algebraic parameter estimation algorithm;

基于构建的固定时间代数参数估计算法,预测输出固定时间代数参数。Based on the constructed fixed-time algebraic parameter estimation algorithm, the fixed-time algebraic parameters are predicted and output.

步骤1包括:Step 1 includes:

Chua电路系统包括三个储能元件(一个电感L和两个电容C1、C2)、一个线性电阻R和一个称为Chua二极管的非线性电阻器,Chua电路系统简化非线性数学模型如下:The Chua circuit system includes three energy storage elements (an inductor L and two capacitors C1 , C2 ), a linear resistor R and a nonlinear resistor called a Chua diode. The simplified nonlinear mathematical model of the Chua circuit system is as follows:

Figure BDA0003612689220000098
Figure BDA0003612689220000098

公式(1)中,R为线性电阻,

Figure BDA0003612689220000099
是C1的电压,
Figure BDA00036126892200000910
是C2的电压,il是通过电感L的电流;
Figure BDA0003612689220000101
是电压
Figure BDA0003612689220000102
的函数,用来表示通过非线性电阻的电流,该非线性函数
Figure BDA0003612689220000103
由如下的奇对称分段线性函数描述:In formula (1), R is the linear resistance,
Figure BDA0003612689220000099
is the voltage ofC1 ,
Figure BDA00036126892200000910
is the voltage of C2 and il is the current through the inductor L;
Figure BDA0003612689220000101
is the voltage
Figure BDA0003612689220000102
A function of , used to represent the current through a nonlinear resistor, the nonlinear function
Figure BDA0003612689220000103
It is described by the following oddly symmetric piecewise linear function:

Figure BDA0003612689220000104
Figure BDA0003612689220000104

公式(2)中m0,m1和Bp为Chua二极管的三个固定常数;In formula (2), m0 , m1 and Bp are three fixed constants of Chua diode;

假设1:在实验过程中,面对变量多或者无法完全复原模型的情况下,可以通过无量纲化的方式,通过一个合适的变量替代,将一个涉及物理量的方程的部分或全部的单位移除,以求简化实验或者计算的目的。Assumption 1: During the experiment, when there are too many variables or the model cannot be completely restored, part or all of the units of an equation involving a physical quantity can be removed by a suitable variable substitution in a dimensionless way. , in order to simplify the experiment or calculation purpose.

公式(1)转化为无量纲形式:Equation (1) is transformed into dimensionless form:

Figure BDA0003612689220000105
Figure BDA0003612689220000105

公式(3)中,y1和y2表示Chua电路系统的输出;f(x1)具体表示为公式(4);β,γ,a,b,x1,x2,x3具体表示公式(5);In formula (3), y1 and y2 represent the output of the Chua circuit system; f(x1 ) is specifically represented by formula (4); β, γ, a, b, x1 , x2 , and x3 specifically represent the formula (5);

f(x1)=ax1+b(|x1+1|-|x1-1|) (4)f(x1 )=ax1 +b(|x1 +1|-|x1 -1|) (4)

Figure BDA0003612689220000106
Figure BDA0003612689220000106

步骤2具体为:Step 2 is specifically:

定义

Figure BDA0003612689220000107
为R≥0到R的希尔伯特空间内的局部平方可积函数,R表示实数,R≥0表示非负实数,局部平方可积函数f(t)经过Volterra积分算子进行映射定义:definition
Figure BDA0003612689220000107
is a local square integrable function in the Hilbert space from R≥ 0 to R, R represents a real number, R≥ 0 represents a non-negative real number, and the local square integrable function f(t) is defined by the Volterra integral operator mapping:

Figure BDA0003612689220000108
Figure BDA0003612689220000108

公式(6)中K(·,·):R×R→R希尔伯特-施密特核函数;K(·,·) in formula (6): R×R→R Hilbert-Schmidt kernel function;

对于公式(6),用微分方程的形式能更好的进行实际计算,基于莱布尼兹微分规则,[VKf](t)由如下微分方程产生:For formula (6), the actual calculation can be better performed in the form of a differential equation. Based on Leibniz’s differential rule, [VK f](t) is generated by the following differential equation:

Figure BDA0003612689220000111
Figure BDA0003612689220000111

式中,K(t,t)为核函数,ξ0为预设的初始值常数;In the formula, K(t, t) is the kernel function, and ξ0 is the preset initial value constant;

定义f(i)(t),i∈Z+为f(t)的i阶导数,Z+为正整数,则[VKf(i)](t)展开为:Define f(i) (t), i∈Z+ is the i-order derivative of f(t), and Z+ is a positive integer, then [VK f(i) ](t) is expanded as:

Figure BDA0003612689220000112
Figure BDA0003612689220000112

其中K(i)(·,·)表示K(·,·)关于第二个变量的i阶导数;where K(i) (·,·) represents the i-order derivative of K(·,·) with respect to the second variable;

核函数K(t,τ)的选取决定了Volterra积分算子的映射[VKf](t)效果,本发明选取核函数为:The selection of the kernel function K(t,τ) determines the mapping [VK f](t) effect of the Volterra integral operator, and the kernel function selected in the present invention is:

Figure BDA0003612689220000113
Figure BDA0003612689220000113

公式(9)中ωh,ω>0均为预设的参数;In formula (9), ωh and ω>0 are all preset parameters;

公式(9)展开得到:Equation (9) is expanded to get:

Figure BDA0003612689220000114
Figure BDA0003612689220000114

其中fp(τ),p=0,1,2满足:Where fp (τ), p = 0, 1, 2 satisfy:

Figure BDA0003612689220000115
Figure BDA0003612689220000115

式中,q={0,1,2};In the formula, q={0,1,2};

基于公式(11),计算获得

Figure BDA0003612689220000116
Based on formula (11), the calculation is obtained
Figure BDA0003612689220000116

Figure BDA0003612689220000121
Figure BDA0003612689220000121

基于公式(10)和公式(12),得到核函数关于的导数:Based on formula (10) and formula (12), the derivative of the kernel function with respect to:

Figure BDA0003612689220000122
Figure BDA0003612689220000122

基于公式(12)与公式(13),得到对于

Figure BDA0003612689220000123
有:Based on formula (12) and formula (13), we get for
Figure BDA0003612689220000123
Have:

K(i)(t,t)=0,K(i)(t,0)=0 (14)。K(i) (t,t)=0, K(i) (t,0)=0 (14).

步骤3中确定固定时间稳定性,具体包括以下步骤:Instep 3, the fixed time stability is determined, which specifically includes the following steps:

性质1:考虑如下Chua电路系统:Property 1: Consider the following Chua circuit system:

Figure BDA0003612689220000124
Figure BDA0003612689220000124

其中f(0)=0,x(t)∈Rn,f(x(t)):Rn→Rn为定义域n维空间Rn到值域n维空间Rn中的一个连续函数,x0=x(t0)表示Chua电路系统的初始值,t0表示初始时刻,为一个定值;where f(0)=0, x(t)∈Rn , f(x(t)): Rn →Rn is a continuous function from the definition domain n-dimensional space Rn to the value domain n-dimensional space Rn , x0 =x(t0 ) represents the initial value of the Chua circuit system, and t0 represents the initial moment, which is a fixed value;

性质2:公式(15)解记为x(t,x0):如果存在一个不依赖于Chua电路系统初始值x0的固定时间T,使得对任意初始状态x0,有

Figure BDA0003612689220000125
成立,那么公式(15)的平衡点x=0是全局固定时间稳定的。Property 2: Formula (15) can be solved as x(t, x0 ): if there is a fixed time T that does not depend on the initial value x0 of the Chua circuit system, such that for any initial state x0 , there is
Figure BDA0003612689220000125
is established, then the equilibrium point x=0 of formula (15) is globally fixed time stable.

在构建固定时间代数参数估计算法之前,需要对Chua电路系统做如下处理:Before constructing a fixed-time algebraic parameter estimation algorithm, the Chua circuit system needs to be processed as follows:

为了方便参数估计器的设计,将公式(3)重新写成关于Chua电路系统输出及其导数形式:In order to facilitate the design of the parameter estimator, formula (3) is rewritten as the output of the Chua circuit system and its derivative form:

Figure BDA0003612689220000126
Figure BDA0003612689220000126

其中y3=|y1+1|-|y1-1|,参数θ1,θ2和θ3满足:where y3 =|y1 +1|-|y1 -1|, the parameters θ1 , θ2 and θ3 satisfy:

θ1=-β(1+a),θ2=β,θ3=-bβ (17)θ1 =-β(1+a), θ2 =β, θ3 =-bβ (17)

对β,a,b,γ的估计转变为对θ123,γ的估计;The estimation of β, a, b, γ is transformed into the estimation of θ1 , θ2 , θ3 , γ;

对公式(16)两边同时进行积分算子VK运算:Perform the integral operator VK operation on both sides of formula (16) at the same time:

Figure BDA0003612689220000131
Figure BDA0003612689220000131

公式(18)中,y1(t),y2(t)是已知的,但其导数

Figure BDA0003612689220000132
是未知的,因此
Figure BDA0003612689220000133
是未知的;将公式(18)展开,并将公式(14)带入公式(18),得到:In formula (18), y1 (t), y2 (t) are known, but their derivatives
Figure BDA0003612689220000132
is unknown, so
Figure BDA0003612689220000133
is unknown; expanding Equation (18) and substituting Equation (14) into Equation (18) yields:

Figure BDA0003612689220000134
Figure BDA0003612689220000134

将公式(19)带入公式(18),得到:Substituting equation (19) into equation (18), we get:

Figure BDA0003612689220000135
Figure BDA0003612689220000135

公式(20)中,除了参数θ123,γ外,其他变量都是已知的,因此可以用来设计参数估计器。然而直接根据公式(6)去计算公式(20)中的Volterra积分算子映射非常繁琐,用微分方程对公式(20)中的所需变量进行计算,更方便实际。In formula (20), except the parameters θ1 , θ2 , θ3 , γ, other variables are known, so they can be used to design the parameter estimator. However, it is very cumbersome to directly calculate the Volterra integral operator mapping in formula (20) according to formula (6), and it is more convenient and practical to use differential equations to calculate the required variables in formula (20).

性质1:对

Figure BDA00036126892200001312
构造微分方程:Property 1: Yes
Figure BDA00036126892200001312
Construct the differential equation:

Figure BDA0003612689220000136
Figure BDA0003612689220000136

其中

Figure BDA0003612689220000137
由公式(13)获得;若核函数K(t,τ)满足公式(9),则
Figure BDA0003612689220000138
Figure BDA0003612689220000139
成立;in
Figure BDA0003612689220000137
Obtained by formula (13); if the kernel function K(t,τ) satisfies formula (9), then
Figure BDA0003612689220000138
right
Figure BDA0003612689220000139
established;

由性质1得到

Figure BDA00036126892200001310
因此公式(20)改写成:obtained fromproperty 1
Figure BDA00036126892200001310
Therefore formula (20) is rewritten as:

Figure BDA00036126892200001311
Figure BDA00036126892200001311

为了方便参数估计器设计,公式(22)重新写为:For the convenience of parameter estimator design, Equation (22) is rewritten as:

Figure BDA0003612689220000141
Figure BDA0003612689220000141

其中in

Figure BDA0003612689220000142
Figure BDA0003612689220000142

其中v1(t),v2(t)满足持续可激励假设。where v1 (t), v2 (t) satisfy the continual motivating hypothesis.

在构建固定时间代数参数估计算法之前,设定如下持续可激励假设:Before building a fixed-time algebraic parameter estimation algorithm, the following persistent motivating assumptions are made:

假设1:函数v1(t)和v2(t)满足持续可激励条件,即存在常数r>0,T0>0使得如下不等式对

Figure BDA0003612689220000143
成立:Assumption 1: The functions v1 (t) and v2 (t) satisfy the continuous excitation condition, that is, there are constants r > 0, T0 > 0 such that the following inequalities
Figure BDA0003612689220000143
Established:

Figure BDA0003612689220000144
Figure BDA0003612689220000144

式中,I3×3为3×3的单位矩阵;In the formula, I3×3 is a 3×3 unit matrix;

根据假设1,公式(23)的两个等式中一个是向量形式,另一个是标量形式,因此对两式分别左乘v1T(t)和v2(t),得到:According toAssumption 1, one of the two equations of formula (23) is in vector form and the other is in scalar form, so the two equations are left-multiplied by v1T (t) and v2 (t), respectively, to get:

Figure BDA0003612689220000145
Figure BDA0003612689220000145

其中in

Figure BDA0003612689220000146
Figure BDA0003612689220000146

对于公式(26)两边同时进行积分算子

Figure BDA0003612689220000147
运算,其中Kg=e-g(t-τ),g>0,得到:For both sides of equation (26), the integral operator is performed simultaneously
Figure BDA0003612689220000147
Operation, where Kg =e-g(t-τ) , g>0, we get:

Figure BDA0003612689220000148
Figure BDA0003612689220000148

其中

Figure BDA0003612689220000149
变量s1,f(t),s2,f(t),v1,f(t)和v2,f(t)由如下Chua电路系统产生:in
Figure BDA0003612689220000149
The variables s1,f (t), s2,f (t), v1,f (t) and v2,f (t) are generated by the following Chua circuit system:

Figure BDA0003612689220000151
Figure BDA0003612689220000151

步骤3中的固定时间代数参数估计算法具体为:The fixed-time algebraic parameter estimation algorithm instep 3 is specifically:

定理1:参数θ,γ的估计值

Figure BDA0003612689220000152
满足:Theorem 1: Estimation of parameters θ, γ
Figure BDA0003612689220000152
Satisfy:

Figure BDA0003612689220000153
Figure BDA0003612689220000153

其中

Figure BDA0003612689220000154
表示参数
Figure BDA0003612689220000155
的上一个值,
Figure BDA0003612689220000156
表示参数
Figure BDA0003612689220000157
的上一个值,其中各变量由公式(21)、(24)、(27)和(29)给出,min{eig(v1,f(t))}和min{eig(v2,f(t))}分别表示v1,f(t)和v2,f(t)特征值最小值。in
Figure BDA0003612689220000154
Indicates parameters
Figure BDA0003612689220000155
the previous value of ,
Figure BDA0003612689220000156
Indicates parameters
Figure BDA0003612689220000157
, where the variables are given by equations (21), (24), (27), and (29), min{eig(v1,f (t))} and min{eig(v2,f (t))} represent the minimum eigenvalues of v1,f (t) and v2,f (t), respectively.

若满足假设1,且参数ε,g满足Ifassumption 1 is satisfied, and the parameters ε, g satisfy

Figure BDA0003612689220000158
Figure BDA0003612689220000158

那么对于任意t≥T0,都满足Then for any t≥T0 , all satisfy

Figure BDA0003612689220000159
Figure BDA0003612689220000159

Figure BDA00036126892200001510
能够在固定时间T0实现对参数θ,γ的估计。which is
Figure BDA00036126892200001510
The estimation of the parameters θ, γ can be achieved at a fixed time T0 .

定理1的证明如下:The proof ofTheorem 1 is as follows:

由v1,f(t)的定义以及假设1可得,当t≥T0From the definition of v1,f (t) andassumption 1, when t ≥ T0 :

Figure BDA00036126892200001511
Figure BDA00036126892200001511

同理可得:The same can be obtained:

Figure BDA00036126892200001512
Figure BDA00036126892200001512

由(33)和(34)可知当t≥T0时v1,f(t)和v2,f(t)可逆,并满足From (33) and (34), it can be known that v1,f (t) and v2,f (t) are invertible when t≥T0 , and satisfy

Figure BDA0003612689220000161
Figure BDA0003612689220000161

因此由公式(28)可知,当t≥T0时,

Figure BDA0003612689220000162
成立,由此定理1得证。Therefore, according to formula (28), when t≥T0 ,
Figure BDA0003612689220000162
is established, andtheorem 1 is proved.

实施例二Embodiment 2

Chua系统电路根据说明书附图1进行建模,参数选择为:The Chua system circuit is modeled according to Figure 1 of the specification, and the parameters are selected as:

x1(0)=-0.9,x2(0)=-0.15,x3(0)=1.47x1 (0)=-0.9, x2 (0)=-0.15, x3 (0)=1.47

γ=27,β=15.6,a=-5/7,b=-3/14γ=27, β=15.6, a=-5/7, b=-3/14

系统的混沌吸引子及其在各坐标平面的投影如说明书附图2所示。固定时间代数参数估计算法按照定理1进行设计,参数选择按照说明书附表1所示。The chaotic attractor of the system and its projection on each coordinate plane are shown in Figure 2 of the description. The fixed-time algebraic parameter estimation algorithm is designed according toTheorem 1, and the parameter selection is shown in Table 1 of the specification.

表1Table 1

Figure BDA0003612689220000163
Figure BDA0003612689220000163

为了验证算法的固定时间收敛性,分别选取两个不同的自适应初始值:In order to verify the fixed-time convergence of the algorithm, two different adaptive initial values are selected respectively:

(C.1):

Figure BDA0003612689220000164
(C.1):
Figure BDA0003612689220000164

(C.2):

Figure BDA0003612689220000165
(C.2):
Figure BDA0003612689220000165

来进行算法的仿真验证。for the simulation verification of the algorithm.

当选取k=10-12,T0=0.6s,可以算出参数估计器的收敛时间Tmax≈1.4s,选取自适应初始值

Figure BDA0003612689220000166
仿真结果如附图6-图9所示。由图6-图9可得,本发明提出算法能够在1.4秒内实现对系统参数的准确估计。当选取一个大的自适应初始值
Figure BDA0003612689220000167
仿真结果如附图10-图13所示,由附图10-图13可得,即使在大的初始状态下,本发明算法依然能够在1.4秒内实现对系统的准确估计,由此验证了算法的固定时间收敛性质。When k=10-12 , T0 =0.6s, the convergence time of the parameter estimator Tmax ≈ 1.4s can be calculated, and the adaptive initial value is selected
Figure BDA0003612689220000166
The simulation results are shown in Figures 6-9. As can be seen from Fig. 6-Fig. 9, the algorithm proposed by the present invention can realize accurate estimation of system parameters within 1.4 seconds. When choosing a large adaptive initial value
Figure BDA0003612689220000167
The simulation results are shown in Figures 10-13, which can be obtained from Figures 10-13. Even in a large initial state, the algorithm of the present invention can still achieve accurate estimation of the system within 1.4 seconds, which verifies The fixed-time convergence properties of the algorithm.

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the present application. It will be understood that each flow and/or block in the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to the processor of a general purpose computer, special purpose computer, embedded processor or other programmable data processing device to produce a machine such that the instructions executed by the processor of the computer or other programmable data processing device produce Means for implementing the functions specified in a flow or flow of a flowchart and/or a block or blocks of a block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory result in an article of manufacture comprising instruction means, the instructions The apparatus implements the functions specified in the flow or flow of the flowcharts and/or the block or blocks of the block diagrams.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded on a computer or other programmable data processing device to cause a series of operational steps to be performed on the computer or other programmable device to produce a computer-implemented process such that The instructions provide steps for implementing the functions specified in the flow or blocks of the flowcharts and/or the block or blocks of the block diagrams.

以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。The above are only the preferred embodiments of the present invention. It should be pointed out that for those skilled in the art, without departing from the technical principle of the present invention, several improvements and modifications can also be made. These improvements and modifications It should also be regarded as the protection scope of the present invention.

Claims (9)

Translated fromChinese
1.Chua电路系统的固定时间代数参数估计算法,其特征在于,包括:1. The fixed time algebraic parameter estimation algorithm of Chua circuit system, is characterized in that, comprises:基于构建的固定时间代数参数估计算法,预测输出固定时间代数参数;Based on the constructed fixed-time algebraic parameter estimation algorithm, predict and output fixed-time algebraic parameters;构建固定时间代数参数估计算法,通过以下步骤实现:Construct a fixed-time algebraic parameter estimation algorithm, implemented through the following steps:根据Chua电路系统,建立非线性数学模型;According to the Chua circuit system, establish a nonlinear mathematical model;利用Volterra积分算子,确定核函数K(t,τ);Use the Volterra integral operator to determine the kernel function K(t,τ);确定固定时间稳定性;determine fixed time stability;根据核函数K(t,τ)和非线性数学模型,构建获得固定时间代数参数估计算法。According to the kernel function K(t,τ) and the nonlinear mathematical model, a fixed-time algebraic parameter estimation algorithm is constructed.2.根据权利要求1所述的Chua电路系统的固定时间代数参数估计算法,其特征在于,2. The fixed-time algebraic parameter estimation algorithm of Chua circuit system according to claim 1, is characterized in that,根据Chua电路系统,构建简化非线性模型,并建立无量纲非线性数学模型,包括:According to the Chua circuit system, a simplified nonlinear model is constructed, and a dimensionless nonlinear mathematical model is established, including:Chua电路系统的简化非线性模型的表达式为:The expression for the simplified nonlinear model of the Chua circuit system is:
Figure FDA0003612689210000011
Figure FDA0003612689210000011
式中,R为线性电阻,
Figure FDA0003612689210000012
是电容器C1的电压,
Figure FDA0003612689210000013
是电容器C2的电压,il是通过电感L的电流;
Figure FDA0003612689210000014
是关于电容器电压的函数,
Figure FDA0003612689210000015
表达式为:
where R is the linear resistance,
Figure FDA0003612689210000012
is the voltage across capacitorC1 ,
Figure FDA0003612689210000013
is the voltage across capacitor C2 andi lis the current through inductor L;
Figure FDA0003612689210000014
is a function of capacitor voltage,
Figure FDA0003612689210000015
The expression is:
Figure FDA0003612689210000016
Figure FDA0003612689210000016
公式(2)中m0,m1和Bp为二极管的三个固定常数;In formula (2), m0 , m1 and Bp are three fixed constants of the diode;利用无量纲形式,构建非线性数学模型;Use dimensionless form to construct nonlinear mathematical model;非线性数学模型的表达式为:The expression of the nonlinear mathematical model is:
Figure FDA0003612689210000017
Figure FDA0003612689210000017
式中,f(x1)=ax1+b(|x1+1|-|x1-1|) (4),In the formula, f(x1 )=ax1 +b(|x1 +1|-|x1 -1|) (4),
Figure FDA0003612689210000021
Figure FDA0003612689210000021
式中,y1和y2表示公式(3)的输出;m0、m1和Bp为Chua二极管的三个固定常数,
Figure FDA0003612689210000022
是电容C1的电压,
Figure FDA0003612689210000023
是电容C2的电压,il是通过电感L的电流。
In the formula, y1 and y2 represent the output of formula (3); m0 , m1 and Bp are the three fixed constants of the Chua diode,
Figure FDA0003612689210000022
is the voltage across capacitorC1 ,
Figure FDA0003612689210000023
is the voltage across capacitor C2 andi lis the current through inductor L.
3.根据权利要求2所述的Chua电路系统的固定时间代数参数估计算法,其特征在于,3. The fixed-time algebraic parameter estimation algorithm of Chua circuit system according to claim 2, is characterized in that,利用Volterra积分算子,确定核函数K(t,τ),包括:Using the Volterra integral operator, determine the kernel function K(t,τ), including:定义
Figure FDA0003612689210000024
为R≥0到R的希尔伯特空间内的局部平方可积函数,R表示实数,R≥0表示非负实数,局部平方可积函数f(t)经过Volterra积分算子进行映射定义:
definition
Figure FDA0003612689210000024
is a local square integrable function in the Hilbert space from R≥ 0 to R, R represents a real number, R≥ 0 represents a non-negative real number, and the local square integrable function f(t) is defined by the Volterra integral operator mapping:
Figure FDA0003612689210000025
Figure FDA0003612689210000025
式中,K(·,·):R×R→R希尔伯特-施密特核函数;In the formula, K(·,·): R×R→R Hilbert-Schmidt kernel function;基于莱布尼兹微分规则,[VKf](t)由如下微分方程产生:Based on Leibniz's differential rule, [VK f](t) is generated by the following differential equation:
Figure FDA0003612689210000026
Figure FDA0003612689210000026
式中,K(t,t)为核函数,ξ0为预设的初始值常数;In the formula, K(t, t) is the kernel function, and ξ0 is the preset initial value constant;定义f(i)(t),i∈Z+为f(t)的i阶导数,Z+为正整数,则[VKf(i)](t)展开为:Define f(i) (t), i∈Z+ is the i-order derivative of f(t), and Z+ is a positive integer, then [VK f(i) ](t) is expanded as:
Figure FDA0003612689210000027
Figure FDA0003612689210000027
式中,K(i)(·,·)表示K(·,·)关于第二个变量的i阶导数;In the formula, K(i) (·,·) represents the i-order derivative of K(·,·) with respect to the second variable;核函数K(t,τ)的选取决定了Volterra积分算子的映射[VKf](t)效果,选取的核函数为:The selection of the kernel function K(t,τ) determines the mapping [VK f](t) effect of the Volterra integral operator, and the selected kernel function is:
Figure FDA0003612689210000028
Figure FDA0003612689210000028
公式(9)中,ωh,ω>0均为预设的参数;In formula (9), ωh , ω>0 are all preset parameters;公式(9)展开得到:Equation (9) is expanded to get:
Figure FDA0003612689210000031
Figure FDA0003612689210000031
其中fp(τ),p=0,1,2满足:Where fp (τ), p = 0, 1, 2 satisfy:
Figure FDA0003612689210000032
Figure FDA0003612689210000032
式中,q={0,1,2};In the formula, q={0,1,2};基于公式(11),计算获得
Figure FDA0003612689210000033
Based on formula (11), the calculation is obtained
Figure FDA0003612689210000033
Figure FDA0003612689210000034
Figure FDA0003612689210000034
基于公式(10)和公式(12),得到核函数关于的导数:Based on formula (10) and formula (12), the derivative of the kernel function with respect to:
Figure FDA0003612689210000035
Figure FDA0003612689210000035
得到对于
Figure FDA0003612689210000036
有:
get for
Figure FDA0003612689210000036
Have:
K(i)(t,t)=0,K(i)(t,0)=0 (14)。K(i) (t,t)=0, K(i) (t,0)=0 (14).
4.根据权利要求3所述的Chua电路系统的固定时间代数参数估计算法,其特征在于,4. The fixed-time algebraic parameter estimation algorithm of Chua circuit system according to claim 3, is characterized in that,确定固定时间稳定性,包括:Determine fixed-time stability, including:Chua电路系统为:The Chua circuit system is:
Figure FDA0003612689210000037
Figure FDA0003612689210000037
式中,f(0)=0,x(t)∈Rn,f(x(t)):Rn→Rn为定义域n维空间Rn到值域n维空间Rn中的一个连续函数,x0=x(t0)表示Chua电路系统的初始值,t0表示初始时刻;In the formula, f(0)=0, x(t)∈Rn , f(x(t)): Rn →Rn is one of the definition domain n-dimensional space Rn to the value domain n-dimensional space Rn Continuous function, x0 =x(t0 ) represents the initial value of the Chua circuit system, and t0 represents the initial moment;将公式(15)的解记为x(t,x0);Denote the solution of formula (15) as x(t,x0 );若存在一个不依赖于Chua电路系统初始值x0的固定时间T,使得对任意初始状态x0,有
Figure FDA0003612689210000041
成立,则判定Chua电路系统的平衡点x=0是全局固定时间稳定的。
If there is a fixed time T that does not depend on the initial value x0 of the Chua circuit system, such that for any initial state x0 , we have
Figure FDA0003612689210000041
If it is established, it is determined that the balance point x=0 of the Chua circuit system is stable in the global fixed time.
5.根据权利要求4所述的Chua电路系统的固定时间代数参数估计算法,其特征在于,5. The fixed-time algebraic parameter estimation algorithm of Chua circuit system according to claim 4, is characterized in that,在构建固定时间代数参数估计算法之前,对Chua电路系统做如下处理:Before constructing a fixed-time algebraic parameter estimation algorithm, the Chua circuit system is processed as follows:将公式(3)变换为关于Chua电路系统输出及其导数形式:Transform equation (3) into the form of the output of the Chua circuit system and its derivative:
Figure FDA0003612689210000042
Figure FDA0003612689210000042
其中,y3=|y1+1|-|y1-1|,where y3 =|y1 +1|-|y1 -1|,参数θ1,θ2和θ3满足:The parameters θ1 , θ2 and θ3 satisfy:θ1=-β(1+a),θ2=β,θ3=-bβ (17);θ1 =-β(1+a), θ2 =β, θ3 =-bβ (17);将对β,a,b,γ的估计转变为对θ123,γ的估计;Convert the estimates of β, a, b, γ into estimates of θ1 , θ2 , θ3 , γ;对公式(16)两边同时进行积分算子VK运算:Perform the integral operator VK operation on both sides of formula (16) at the same time:
Figure FDA0003612689210000043
Figure FDA0003612689210000043
式中,y1(t),y2(t)为公式(3)的输出;In the formula, y1 (t), y2 (t) is the output of formula (3);将公式(14)带入公式(18),得到:Substituting equation (14) into equation (18), we get:
Figure FDA0003612689210000044
Figure FDA0003612689210000044
将公式(19)带入公式(18),得到:Substituting equation (19) into equation (18), we get:
Figure FDA0003612689210000045
Figure FDA0003612689210000045
基于公式(20),设计参数估计器;Based on formula (20), design a parameter estimator;
Figure FDA0003612689210000051
构造微分方程:
right
Figure FDA0003612689210000051
Construct the differential equation:
Figure FDA0003612689210000052
Figure FDA0003612689210000052
其中,
Figure FDA0003612689210000053
由公式(13)获得;
in,
Figure FDA0003612689210000053
Obtained by formula (13);
若核函数K(t,τ)满足公式(9),则
Figure FDA0003612689210000054
Figure FDA0003612689210000055
成立;
If the kernel function K(t,τ) satisfies formula (9), then
Figure FDA0003612689210000054
right
Figure FDA0003612689210000055
established;
由微分方程得到
Figure FDA0003612689210000056
因此公式(20)改写成:
obtained from the differential equation
Figure FDA0003612689210000056
Therefore formula (20) is rewritten as:
Figure FDA0003612689210000057
Figure FDA0003612689210000057
其中,in,
Figure FDA0003612689210000058
Figure FDA0003612689210000058
其中,v1(t)和v2(t)满足持续可激励假设。Among them, v1 (t) and v2 (t) satisfy the continuous stimuli hypothesis.
6.根据权利要求5所述的Chua电路系统的固定时间代数参数估计算法,其特征在于,6. The fixed time algebraic parameter estimation algorithm of Chua circuit system according to claim 5, is characterized in that,在构建固定时间代数参数估计算法之前,设定如下持续可激励假设:Before building a fixed-time algebraic parameter estimation algorithm, the following persistent motivating assumptions are made:假设1:函数v1(t)和v2(t)满足持续可激励条件,存在常数r>0,T0>0使得如下不等式对
Figure FDA0003612689210000059
成立:
Assumption 1: The functions v1 (t) and v2 (t) satisfy the continuous excitation condition, and there are constants r > 0 and T0 > 0 such that the following inequalities
Figure FDA0003612689210000059
Established:
Figure FDA00036126892100000510
Figure FDA00036126892100000510
式中,I3×3为3×3的单位矩阵;In the formula, I3×3 is a 3×3 unit matrix;对公式(23)分别左乘v1T(t)和v2(t),得到:Left-multiplying v1T (t) and v2 (t) to formula (23), respectively, yields:
Figure FDA00036126892100000511
Figure FDA00036126892100000511
其中in
Figure FDA0003612689210000061
Figure FDA0003612689210000061
对于公式(26)两边同时进行积分算子
Figure FDA0003612689210000062
运算,其中Kg=e-g(t-τ),g>0,得到:
For both sides of equation (26), the integral operator is performed simultaneously
Figure FDA0003612689210000062
Operation, where Kg =e-g(t-τ) , g>0, we get:
Figure FDA0003612689210000063
Figure FDA0003612689210000063
其中,
Figure FDA0003612689210000064
in,
Figure FDA0003612689210000064
Chua电路系统产生变量s1,f(t)、s2,f(t)、v1,f(t)和v2,f(t):The Chua circuitry produces variables s1,f (t), s2,f (t), v1,f (t) and v2,f (t):
Figure FDA0003612689210000065
Figure FDA0003612689210000065
7.根据权利要求6所述的Chua电路系统的固定时间代数参数估计算法,其特征在于,7. The fixed-time algebraic parameter estimation algorithm of Chua circuit system according to claim 6, is characterized in that,根据核函数K(t,τ)和非线性数学模型,构建获得固定时间代数参数估计算法,包括:According to the kernel function K(t,τ) and the nonlinear mathematical model, a fixed-time algebraic parameter estimation algorithm is constructed, including:固定时间代数参数估计算法的表达式为:The expression for the fixed-time algebraic parameter estimation algorithm is:
Figure FDA0003612689210000066
Figure FDA0003612689210000066
参数ε,g满足:
Figure FDA0003612689210000067
其中,
Figure FDA0003612689210000068
表示参数
Figure FDA0003612689210000069
的上一个值,
Figure FDA00036126892100000610
表示参数
Figure FDA00036126892100000611
的上一个值,min{eig(v1,f(t))}和min{eig(v2,f(t))}分别表示v1,f(t)和v2,f(t)特征值最小值。
Parameters ε, g satisfy:
Figure FDA0003612689210000067
in,
Figure FDA0003612689210000068
Indicates parameters
Figure FDA0003612689210000069
the previous value of ,
Figure FDA00036126892100000610
Indicates parameters
Figure FDA00036126892100000611
The previous value of , min{eig(v1,f (t))} and min{eig(v2,f (t))} denote v1,f (t) and v2,f (t) features respectively value minimum.
8.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至7中任一项所述方法的步骤。8. An electronic device comprising a memory, a processor and a computer program stored on the memory and running on the processor, wherein the processor implements any one of claims 1 to 7 when the processor executes the program the steps of the method described in item.9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至7中任一项所述方法的步骤。9. A computer-readable storage medium on which a computer program is stored, characterized in that, when the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.
CN202210435371.7A2022-04-242022-04-24 Fixed-time Algebraic Parameter Estimation Algorithm for Chua Circuit SystemsActiveCN114741657B (en)

Priority Applications (1)

Application NumberPriority DateFiling DateTitle
CN202210435371.7ACN114741657B (en)2022-04-242022-04-24 Fixed-time Algebraic Parameter Estimation Algorithm for Chua Circuit Systems

Applications Claiming Priority (1)

Application NumberPriority DateFiling DateTitle
CN202210435371.7ACN114741657B (en)2022-04-242022-04-24 Fixed-time Algebraic Parameter Estimation Algorithm for Chua Circuit Systems

Publications (2)

Publication NumberPublication Date
CN114741657Atrue CN114741657A (en)2022-07-12
CN114741657B CN114741657B (en)2025-06-24

Family

ID=82283655

Family Applications (1)

Application NumberTitlePriority DateFiling Date
CN202210435371.7AActiveCN114741657B (en)2022-04-242022-04-24 Fixed-time Algebraic Parameter Estimation Algorithm for Chua Circuit Systems

Country Status (1)

CountryLink
CN (1)CN114741657B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
KR101219545B1 (en)*2011-09-142013-01-09주식회사 파워이십일Optimized parameter estimation method for power system
WO2017097909A1 (en)*2015-12-112017-06-15Bundesrepublik Deutschland, Vertreten Durch Das Bundesministerium Des Innern, Vertreten Durch Das Bundesamt Für Sicherheit In Der InformationstechnikChaotic circuit having variable dynamic states as secure information memory
CN107092746A (en)*2017-04-192017-08-25江西理工大学A kind of circuit design method of the isomery magnetic control memristor model based on Chua circuits
CN111106926A (en)*2019-12-192020-05-05宁波海特技术转移有限公司Working method for discrete chaotic sequence control continuous chaotic system
CN112506045A (en)*2020-09-282021-03-16中国科学院数学与系统科学研究院Adaptive control method and system of non-standard discrete time nonlinear system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
KR101219545B1 (en)*2011-09-142013-01-09주식회사 파워이십일Optimized parameter estimation method for power system
WO2017097909A1 (en)*2015-12-112017-06-15Bundesrepublik Deutschland, Vertreten Durch Das Bundesministerium Des Innern, Vertreten Durch Das Bundesamt Für Sicherheit In Der InformationstechnikChaotic circuit having variable dynamic states as secure information memory
CN107092746A (en)*2017-04-192017-08-25江西理工大学A kind of circuit design method of the isomery magnetic control memristor model based on Chua circuits
CN111106926A (en)*2019-12-192020-05-05宁波海特技术转移有限公司Working method for discrete chaotic sequence control continuous chaotic system
CN112506045A (en)*2020-09-282021-03-16中国科学院数学与系统科学研究院Adaptive control method and system of non-standard discrete time nonlinear system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李建芬 等: "不确定Chua’s电路的参数辨识与自适应同步", 物理学报, no. 12, 15 December 2008 (2008-12-15)*
石尚 等: "不确定Chua电路系统的固定时间自适应参数估计", 控制理论与应用, vol. 39, no. 8, 31 August 2022 (2022-08-31)*

Also Published As

Publication numberPublication date
CN114741657B (en)2025-06-24

Similar Documents

PublicationPublication DateTitle
CN105871356B (en)Adaptive filter method based on maximum mixing cross-correlation entropy criterion
JP6404909B2 (en) How to calculate the output model of a technical system
CN104009696B (en)Interactive model reference adaptive speed and stator resistance identification method based on sliding-mode control
CN107181474A (en)A kind of kernel adaptive algorithm filter based on functional expansion
CN105425587A (en)Hysteresis nonlinear motor identification and control method
CN111881764B (en)Target detection method and device, electronic equipment and storage medium
JP7090734B2 (en) Control system, control method and storage medium
CN116987878A (en)Anti-interference temperature control method for heat treatment of aerospace complex cast parts
CN117252102A (en)Nonlinear wave propagation simulation method and system based on physical information neural network
KR20210115863A (en)Method and appartus of parallel processing for neural network model
CN114037052A (en)Training method and device for detection model, electronic equipment and storage medium
CN115330579B (en) Construction method, device, equipment and storage medium of model watermark
CN110362881B (en)Microwave power device nonlinear model method based on extreme learning machine
CN114741657A (en)Fixed time algebraic parameter estimation algorithm of Chua circuit system
CN113361719B (en)Incremental learning method and image processing method based on image processing model
CN119046887A (en)Complex electromechanical equipment fault diagnosis method and system based on confidence rule base
CN111159956B (en)Feature-based flow field discontinuity capturing method
Michalowsky et al.Model-based extremum seeking for a class of nonlinear systems
LiPeriodic solutions of non-autonomous cellular neural networks with impulses and delays on time scales
Song et al.Ricci Planner: Zero-Shot Transfer for Goal-Conditioned Reinforcement Learning via Geometric Flow
Chen et al.Modeling of nonlinear dynamical systems based on deterministic learning and structural stability
Bechtold et al.Nonlinear model order reduction in nanoelectronics: combination of POD and TPWL
ZhangNonlinear system identification with output error model through stabilized simulation
ShuEfficient algorithms for solving partial differential equations with discontinuous solutions
US20230142773A1 (en)Method and system for real-time simulations using convergence stopping criterion

Legal Events

DateCodeTitleDescription
PB01Publication
PB01Publication
SE01Entry into force of request for substantive examination
SE01Entry into force of request for substantive examination
GR01Patent grant
GR01Patent grant

[8]ページ先頭

©2009-2025 Movatter.jp