CN116227324A - A Fractional Order Memristive Neural Network Estimation Method Under Variance Constraint - Google Patents

A Fractional Order Memristive Neural Network Estimation Method Under Variance Constraint Download PDF

Info

Publication number
CN116227324A
CN116227324A CN202211559637.5A CN202211559637A CN116227324A CN 116227324 A CN116227324 A CN 116227324A CN 202211559637 A CN202211559637 A CN 202211559637A CN 116227324 A CN116227324 A CN 116227324A
Authority
CN
China
Prior art keywords
matrix
kth moment
row
kth
neural network
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
CN202211559637.5A
Other languages
Chinese (zh)
Other versions
CN116227324B (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.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and 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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN202211559637.5A priority Critical patent/CN116227324B/en
Publication of CN116227324A publication Critical patent/CN116227324A/en
Application granted granted Critical
Publication of CN116227324B publication Critical patent/CN116227324B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/06Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
    • G06N3/063Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using electronic means
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Biophysics (AREA)
  • Computational Mathematics (AREA)
  • Computing Systems (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Biomedical Technology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Algebra (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Neurology (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Complex Calculations (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明公开了一种方差受限下的分数阶忆阻神经网络估计方法,所述包括如下步骤:步骤一、建立分数阶忆阻神经网络动态模型;步骤二、放大转发协议下对分数阶忆阻神经网络动态模型进行状态估计;步骤三、计算分数阶忆阻神经网络的误差协方差矩阵的上界及H性能约束条件;步骤四、利用随机分析方法,并通过解线性矩阵不等式求解出估计器增益矩阵Kk的解,实现对放大转发协议下分数阶忆阻神经网络动态模型的状态估计,判断k+1是否达到总时长N,若k+1<N,则执行步骤二,反之结束。本发明解决了现有状态估计方法不能同时处理放大转发协议下具有H性能约束及方差受限分数阶忆阻神经网络的状态估计导致的估计性能准确率低的问题,从而提高了估计性能的准确率。

Figure 202211559637

The invention discloses a method for estimating a fractional-order memristive neural network under variance limitation, which comprises the following steps: Step 1, establishing a dynamic model of a fractional-order memristive neural network; The dynamic model of the memristive neural network is used to estimate the state; Step 3, calculate the upper bound of the error covariance matrix of the fractional order memristive neural network and the H performance constraints; Step 4, use the stochastic analysis method and solve the linear matrix inequality to obtain The solution of the estimator gain matrix K k realizes the state estimation of the dynamic model of the fractional-order memristive neural network under the amplification and forwarding protocol, and judges whether k+1 reaches the total duration N. If k+1<N, execute step 2, otherwise Finish. The invention solves the problem that the existing state estimation method cannot simultaneously deal with the low accuracy of the estimation performance caused by the state estimation of the fractional-order memristive neural network with H performance constraints and variance constraints under the amplification and forwarding protocol, thereby improving the estimation performance. Accuracy.

Figure 202211559637

Description

一种方差受限下的分数阶忆阻神经网络估计方法A variance-constrained fractional-order memristor neural network estimation method

技术领域Technical Field

本发明涉及一种神经网络的状态估计方法,具体涉及一种放大转发协议下具有H性能约束及方差受限的分数阶忆阻神经网络的状态估计方法。The invention relates to a state estimation method of a neural network, and in particular to a state estimation method of a fractional-order memristor neural network with H∞ performance constraints and limited variance under an amplify-and-forward protocol.

背景技术Background Art

神经网络是根据人脑中的神经细胞结构和功能模拟出来的信息处理系统,具有较强的联想能力、自适应性和容错能力等优势。在现实的许多网络中,这类网络能高效地解决模式识别、信号处理和图像识别等实际系统建模和分析方面。Neural networks are information processing systems that simulate the structure and function of nerve cells in the human brain. They have the advantages of strong associative ability, adaptability and fault tolerance. In many real networks, this type of network can efficiently solve practical system modeling and analysis aspects such as pattern recognition, signal processing and image recognition.

在过去的几十年里,递归神经网络的状态估计问题已成为一个引人关注的课题,它已成功地应用于联想记忆、模式识别和组合优化等广泛领域。然而,在实际应用过程中,神经元的信息往往是不完全可测的,因此需要使用有效的估计方法来估计它们。到目前为止,已经研究了许多不同类型的神经网络状态估计问题。但值得注意的是,目前的结果仅适用于定常的情况下,这可能会导致应用存在的局限性。In the past few decades, the state estimation problem of recurrent neural networks has become an interesting topic, which has been successfully applied to a wide range of fields such as associative memory, pattern recognition and combinatorial optimization. However, in practical applications, the information of neurons is often not completely measurable, so effective estimation methods are needed to estimate them. So far, many different types of neural network state estimation problems have been studied. But it is worth noting that the current results are only applicable to steady-state cases, which may lead to limitations in application.

目前已有的状态估计方法不能同时处理在方差受限下具有H性能约束及放大转发协议的分数阶忆阻神经网络的状态估计问题,导致估计性能准确率低。The existing state estimation methods cannot simultaneously handle the state estimation problem of fractional-order memristor neural networks with H∞ performance constraints and amplify-and-forward protocols under variance constraints, resulting in low estimation performance accuracy.

发明内容Summary of the invention

本发明针对时变系统进行研究,提供了一种方差受限下的分数阶忆阻神经网络估计方法。该方法解决了现有状态估计方法不能同时处理放大转发协议下具有H性能约束的分数阶忆阻神经网络的状态估计问题,从而导致估计精度准确率低,以及在放大转发协议下存在信息无法接收到其他时刻信息的情况下,导致估计性能准确率低的问题,可用于忆阻神经网络状态估计领域。The present invention studies time-varying systems and provides a fractional-order memristor neural network estimation method under variance constraints. The method solves the problem that the existing state estimation method cannot simultaneously process the state estimation problem of the fractional-order memristor neural network with H performance constraints under the amplification and forwarding protocol, resulting in low estimation accuracy, and the problem that when there is information under the amplification and forwarding protocol that cannot receive information at other times, the estimation performance accuracy is low, and can be used in the field of memristor neural network state estimation.

本发明的目的是通过以下技术方案实现的:The objective of the present invention is achieved through the following technical solutions:

一种方差受限下的分数阶忆阻神经网络估计方法,包括如下步骤:A variance-constrained fractional-order memristor neural network estimation method comprises the following steps:

步骤一、建立放大转发协议下的分数阶忆阻神经网络动态模型;Step 1: Establish a dynamic model of fractional-order memristor neural network under amplification and forwarding protocol;

步骤二、放大转发协议下对步骤一建立的分数阶忆阻神经网络动态模型进行状态估计;Step 2: performing state estimation on the dynamic model of the fractional-order memristor neural network established in step 1 under the amplification and forwarding protocol;

步骤三、给定H性能指标γ、半正定矩阵一号

Figure BDA0003984088630000021
半正定矩阵二号
Figure BDA0003984088630000022
及初始条件
Figure BDA0003984088630000023
计算分数阶忆阻神经网络的误差协方差矩阵的上界及H性能约束条件;Step 3: Given the H∞ performance index γ and the semi-positive definite matrix No.
Figure BDA0003984088630000021
Semi-positive definite matrix II
Figure BDA0003984088630000022
and initial conditions
Figure BDA0003984088630000023
Calculate the upper bound of the error covariance matrix and H∞ performance constraints of fractional-order memristor neural networks;

步骤四、利用随机分析方法,并通过解线性矩阵不等式求解出估计器增益矩阵Kk的解,实现对放大转发协议下的分数阶忆阻神经网络动态模型的状态估计,判断k+1是否达到总时长N,若k+1<N,则执行步骤二,反之结束。Step 4: Use the random analysis method and solve the linear matrix inequality to solve the solution of the estimator gain matrix Kk to realize the state estimation of the dynamic model of the fractional-order memristor neural network under the amplification and forwarding protocol, and judge whether k+1 reaches the total duration N. If k+1<N, execute step 2, otherwise end.

本发明中,所述神经网络可以为车辆悬挂构成的网络、质点弹簧构成的网络、航天器构成的网络或雷达构成的网络,在生物学、数学、计算机以及联想记忆、模式识别、组合优化和图像处理等多学科领域具有重要应用。In the present invention, the neural network can be a network composed of vehicle suspension, a network composed of mass springs, a network composed of spacecraft or a network composed of radar, and has important applications in multiple disciplines such as biology, mathematics, computers, associative memory, pattern recognition, combinatorial optimization and image processing.

相比于现有技术,本发明具有如下优点:Compared with the prior art, the present invention has the following advantages:

1、本发明同时考虑了放大转发协议下具有H性能约束和方差受限对状态估计性能的影响,利用不等式处理技术以及随机分析方法,全面考虑了估计误差协方差矩阵的有效信息,与现有的神经网络状态估计方法相比,本发明的分数阶忆阻神经网络状态估计方法同时考虑在放大转发协议下具有H性能约束及方差受限的分数阶忆阻神经网络的状态估计问题,得到了误差系统同时满足估计误差协方差有上界和给定的H性能要求的分数阶忆阻神经网络状态估计方法,同时达到了抑制扰动,并且提高了估计精度的目的,而且目前的结果仅适用于定常的情况下,这可能会导致应用存在的局限性,本发明考虑了时变神经网络,更能接近实际。1. The present invention simultaneously considers the influence of H∞ performance constraint and variance constraint on state estimation performance under the amplification and forwarding protocol, and comprehensively considers the effective information of the estimation error covariance matrix by using inequality processing technology and random analysis method. Compared with the existing neural network state estimation method, the fractional-order memristor neural network state estimation method of the present invention simultaneously considers the state estimation problem of the fractional-order memristor neural network with H∞ performance constraint and variance constraint under the amplification and forwarding protocol, and obtains the fractional-order memristor neural network state estimation method in which the error system simultaneously satisfies the upper bound of the estimation error covariance and the given H∞ performance requirements, and achieves the purpose of suppressing disturbances and improving the estimation accuracy. Moreover, the current results are only applicable to steady-state situations, which may lead to limitations in application. The present invention considers time-varying neural networks and is closer to reality.

2、本发明解决了现有状态估计方法不能同时处理放大转发协议下具有H性能约束及方差受限分数阶忆阻神经网络的状态估计导致的估计性能准确率低的问题,从而提高了估计性能的准确率。从仿真图可以看出,功率越小,分数阶忆阻神经网络的状态估计性能逐渐降低,估计误差相对较大。此外,验证了本发明所提出的状态估计方法的可行性和有效性。2. The present invention solves the problem that the existing state estimation method cannot simultaneously handle the state estimation of the fractional-order memristor neural network with H∞ performance constraints and variance constrained under the amplification and forwarding protocol, resulting in low estimation performance accuracy, thereby improving the estimation performance accuracy. It can be seen from the simulation diagram that the smaller the power, the state estimation performance of the fractional-order memristor neural network gradually decreases, and the estimation error is relatively large. In addition, the feasibility and effectiveness of the state estimation method proposed in the present invention are verified.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1是本发明放大转发协议下的分数阶忆阻神经网络状态估计方法的流程图;FIG1 is a flow chart of a method for estimating a state of a fractional-order memristor neural network under an amplify-and-forward protocol of the present invention;

图2是分数阶忆阻神经网络实际状态轨迹zk在两个不同情形下状态估计轨迹

Figure BDA0003984088630000031
的对比图,zk为神经网络在第k时刻的状态变量;其中
Figure BDA0003984088630000032
是系统状态轨迹,
Figure BDA0003984088630000033
是情形一下的状态估计轨迹,
Figure BDA0003984088630000034
是情形二下的状态估计轨迹;Figure 2 is the actual state trajectory z k of the fractional-order memristor neural network under two different situations.
Figure BDA0003984088630000031
Comparison chart, z k is the state variable of the neural network at the kth moment;
Figure BDA0003984088630000032
is the system state trajectory,
Figure BDA0003984088630000033
is the state estimation trajectory for the case,
Figure BDA0003984088630000034
is the state estimation trajectory under case 2;

图3是神经网络控制输出估计误差轨迹图在两个不同情形下的误差对比图;其中

Figure BDA0003984088630000035
是情形一下的控制输出估计误差轨迹,
Figure BDA0003984088630000036
是情形二下的控制输出估计误差轨迹;FIG3 is a comparison diagram of the error trajectory of the neural network control output estimation error in two different situations;
Figure BDA0003984088630000035
is the control output estimation error trajectory for the case 1,
Figure BDA0003984088630000036
is the control output estimation error trajectory under case 2;

图4是神经网络实际状态误差协方差和误差协方差上界第一个分量的轨迹图;其中

Figure BDA0003984088630000037
是方差约束的轨迹,
Figure BDA0003984088630000038
是实际误差协方差的轨迹;FIG4 is a trajectory diagram of the actual state error covariance of the neural network and the first component of the upper bound of the error covariance;
Figure BDA0003984088630000037
is the variance-constrained trajectory,
Figure BDA0003984088630000038
is the locus of the actual error covariance;

图5是神经网络实际状态误差协方差和误差协方差上界第二个分量的轨迹图;其中

Figure BDA0003984088630000039
是方差约束的轨迹,
Figure BDA00039840886300000310
是实际误差协方差的轨迹。FIG5 is a trajectory diagram of the actual state error covariance of the neural network and the second component of the upper bound of the error covariance;
Figure BDA0003984088630000039
is the variance-constrained trajectory,
Figure BDA00039840886300000310
is the locus of the actual error covariance.

具体实施方式DETAILED DESCRIPTION

下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。The technical solution of the present invention is further described below in conjunction with the accompanying drawings, but is not limited thereto. Any modification or equivalent replacement of the technical solution of the present invention without departing from the spirit and scope of the technical solution of the present invention should be included in the protection scope of the present invention.

本发明提供了一种方差受限下的分数阶忆阻神经网络估计方法,所述方法利用随机分析方法和不等式处理技术,首先,分别考虑估计误差系统满足H性能约束条件及误差协方差有上界的充分条件;然后,再同时得到估计误差系统满足H性能约束条件及误差协方差有上界的判别条件;最后,通过求解一系列线性矩阵不等式得到估计器增益矩阵的值,实现了在放大转发协议下具有H性能约束以及方差受限同时发生的情况下性能估计不受影响,从而提高了估计准确率。如图1所示,具体包括如下步骤:The present invention provides a fractional-order memristor neural network estimation method under variance constraint. The method utilizes a random analysis method and an inequality processing technique. First, sufficient conditions for the estimated error system to satisfy the H∞ performance constraint and the error covariance to have an upper bound are considered respectively; then, the judgment conditions for the estimated error system to satisfy the H∞ performance constraint and the error covariance to have an upper bound are obtained simultaneously; finally, the value of the estimator gain matrix is obtained by solving a series of linear matrix inequalities, so that the performance estimation is not affected when the H∞ performance constraint and variance constraint occur simultaneously under the amplification and forwarding protocol, thereby improving the estimation accuracy. As shown in FIG1 , the method specifically comprises the following steps:

步骤一、建立放大转发协议下分数阶忆阻神经网络动态模型。具体步骤如下:Step 1: Establish a dynamic model of fractional-order memristor neural network under the amplify-and-forward protocol. The specific steps are as follows:

首先,介绍Grunwald-Letnikov分数阶导数定义,这是一种适合于数值实现和应用的形式。该定义的离散形式表示为:First, we introduce the Grunwald-Letnikov fractional derivative definition, which is a form suitable for numerical implementation and application. The discrete form of this definition is expressed as:

Figure BDA0003984088630000041
Figure BDA0003984088630000041

Figure BDA0003984088630000042
Figure BDA0003984088630000042

式中,Δα表示α阶的Grunwald-Letnikov分数阶导数定义,h为相应的采样间隔,假设采样间隔为1,k为采样时刻,

Figure BDA0003984088630000045
表示h→0的所有极限值,i!表示的是i的所有阶层,
Figure BDA0003984088630000043
表示i=0到k的所有求和值。In the formula, Δ α represents the definition of the Grunwald-Letnikov fractional derivative of order α, h is the corresponding sampling interval, assuming that the sampling interval is 1, and k is the sampling time.
Figure BDA0003984088630000045
represents all the extreme values of h→0, i! represents all the classes of i,
Figure BDA0003984088630000043
represents all the summed values from i=0 to k.

根据Grunwald-Letnikov分数阶导数定义,分数阶忆阻神经网络动态模型的状态空间形式为:According to the definition of Grunwald-Letnikov fractional-order derivative, the state space form of the dynamic model of the fractional-order memristor neural network is:

Figure BDA0003984088630000044
Figure BDA0003984088630000044

式中:Where:

Figure BDA0003984088630000051
Figure BDA0003984088630000051

Figure BDA0003984088630000052
Figure BDA0003984088630000052

Figure BDA0003984088630000053
Figure BDA0003984088630000053

这里,

Figure BDA0003984088630000054
表示微分算子,
Figure BDA0003984088630000055
为分数阶(j=1,2,…,n),n为维数,
Figure BDA0003984088630000056
是在第k时刻的分数阶忆阻神经网络的状态向量,
Figure BDA0003984088630000057
是在第k-ι+1时刻的分数阶忆阻神经网络的状态向量,
Figure BDA0003984088630000058
是在第k-d时刻的分数阶忆阻神经网络的状态向量,
Figure BDA0003984088630000059
是在第k+1时刻的分数阶忆阻神经网络的状态向量,
Figure BDA00039840886300000510
为神经网络动态模型状态的实数域且其维数为n;
Figure BDA00039840886300000511
为在第k时刻的被控测量输出,
Figure BDA00039840886300000512
为神经网络动态模型被控输出状态的实数域且其维数为r;
Figure BDA00039840886300000513
是给定的初始序列,d为离散固定的网络时滞;A(xk)=diagn{ai(xi,k)}为在第k时刻的神经网络自反馈对角矩阵,n为维数,diag{·}表示的是对角矩阵,ai(xi,k)为A(xk)的第i个分量,n为维数;Ad(xk)={aij,d(xi,k)}n*n为在第k时刻的已知维数且与时滞相关的系统矩阵,aij,d(xi,k)为在第k时刻Ad(xk)的第i个分量形式;B(xk)={bij(xi,k)}n*n为在第k时刻的已知的连接激励函数的权重矩阵,bij(xi,k)为在第k时刻B(xk)的第i个分量形式;f(xk)为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;v1k为在第k时刻均值为零并且协方差为V1>0的高斯白噪声序列,v2k为在第k时刻均值为零并且协方差为V2>0的高斯白噪声序列,
Figure BDA00039840886300000514
表示的是ι=1到k+1求和的值。here,
Figure BDA0003984088630000054
represents the differential operator,
Figure BDA0003984088630000055
is a fractional order (j=1,2,…,n), n is the dimension,
Figure BDA0003984088630000056
is the state vector of the fractional-order memristor neural network at the kth moment,
Figure BDA0003984088630000057
is the state vector of the fractional-order memristor neural network at the k-ι+1th moment,
Figure BDA0003984088630000058
is the state vector of the fractional-order memristor neural network at the kdth time,
Figure BDA0003984088630000059
is the state vector of the fractional-order memristor neural network at the k+1th time,
Figure BDA00039840886300000510
is the real number field of the states of the neural network dynamic model and its dimension is n;
Figure BDA00039840886300000511
is the controlled measured output at the kth moment,
Figure BDA00039840886300000512
is the real number domain of the controlled output state of the neural network dynamic model and its dimension is r;
Figure BDA00039840886300000513
is a given initial sequence, d is a discrete fixed network time delay; A(x k ) = diag n {a i (xi ,k )} is the neural network self-feedback diagonal matrix at the kth moment, n is the dimension, diag{·} represents a diagonal matrix, a i (xi ,k ) is the i-th component of A(x k ), and n is the dimension; A d (x k ) = {a ij,d (xi ,k )} n*n is the system matrix of known dimension and time delay at the kth moment, a ij,d (xi ,k ) is the i-th component form of A d (x k ) at the kth moment; B(x k ) = {b ij (xi ,k )} n*n is the weight matrix of the known connection activation function at the kth moment, b ij ( xi,k ) is the i-th component form of B(x k ) at the kth moment; f(x k ) is the nonlinear activation function at the kth moment; C C 1k is the noise distribution matrix of the first component known system at the kth moment, C 2k is the noise distribution matrix of the second component known system at the kth moment, H k is the adjustment matrix of the known measurement at the kth moment; D k is the measurement matrix of the known measurement at the kth moment; v 1k is a Gaussian white noise sequence with zero mean and covariance V 1 >0 at the kth moment, v 2k is a Gaussian white noise sequence with zero mean and covariance V 2 >0 at the kth moment,
Figure BDA00039840886300000514
It represents the sum of the values from ι=1 to k+1.

状态依赖矩阵参数ai(xi,k)、aij,d(xi,k)和bij(xi,k)满足:The state dependence matrix parameters ai (xi ,k ), aij ,d (xi ,k ) and bij (xi ,k ) satisfy:

Figure BDA00039840886300000515
Figure BDA00039840886300000515

式中,ai(xi,k)、aij,d(xi,k)和bij(xi,k)分别为A(xk),Ad(xk)和B(xk)的第i个分量,Ωi>0为已知的切换阈值,

Figure BDA0003984088630000061
为第i个已知的上存储变量矩阵,
Figure BDA0003984088630000062
为第i个已知的下存储变量矩阵,
Figure BDA0003984088630000063
为第ij,d个已知的左存储变量矩阵,
Figure BDA0003984088630000064
为第ij,d个已知的右存储变量矩阵,
Figure BDA0003984088630000065
为第ij个已知的内存储变量矩阵,
Figure BDA0003984088630000066
为第ij个已知的外存储变量矩阵。Where ai (xi ,k ), aij ,d (xi ,k ) and bij (xi ,k ) are the i-th components of A( xk ), Ad ( xk ) and B( xk ), respectively, Ωi >0 is the known switching threshold,
Figure BDA0003984088630000061
is the i-th known upper storage variable matrix,
Figure BDA0003984088630000062
is the i-th known storage variable matrix,
Figure BDA0003984088630000063
is the ij,dth known left storage variable matrix,
Figure BDA0003984088630000064
is the ij,dth known right storage variable matrix,
Figure BDA0003984088630000065
is the ijth known internal storage variable matrix,
Figure BDA0003984088630000066
is the ijth known external storage variable matrix.

定义:definition:

Figure BDA0003984088630000067
Figure BDA0003984088630000067

Figure BDA0003984088630000068
Figure BDA0003984088630000068

Figure BDA0003984088630000069
Figure BDA0003984088630000069

式中,

Figure BDA00039840886300000610
为第i个最小存储的第一号度量矩阵,
Figure BDA00039840886300000611
为第i个已知的上存储区间变量矩阵,
Figure BDA00039840886300000612
为第i个已知的下存储区间变量矩阵,min{·}表示两个存储矩阵中取最小值,max{·}表示两个存储矩阵中取最大值,
Figure BDA00039840886300000613
为第i个最大存储的第一号度量矩阵,
Figure BDA00039840886300000614
为第ij,d个最小存储的第二号度量矩阵,
Figure BDA00039840886300000615
为第i个最大存储的第二号度量矩阵,
Figure BDA00039840886300000616
为第ij,d个已知的左存储变量矩阵,
Figure BDA00039840886300000617
为第ij,d个已知的右存储变量矩阵,
Figure BDA00039840886300000618
为第ij个最小存储的第三号度量矩阵,
Figure BDA00039840886300000619
为第ij个最大存储的第三号度量矩阵,
Figure BDA00039840886300000620
为第ij个已知的内存储变量矩阵,
Figure BDA00039840886300000621
为第ij个已知的外存储变量矩阵,diag{·}为对角矩阵,A-为定义的第一号对角矩阵,A+为定义的第二号对角矩阵,
Figure BDA00039840886300000622
为定义的第三号对角矩阵,
Figure BDA00039840886300000623
为定义的第四号对角矩阵,B-为定义的第五号对角矩阵,B+为定义的第六号对角矩阵,n为维数。In the formula,
Figure BDA00039840886300000610
is the first metric matrix with the smallest storage of i,
Figure BDA00039840886300000611
is the i-th known upper storage interval variable matrix,
Figure BDA00039840886300000612
is the i-th known lower storage interval variable matrix, min{·} means taking the minimum value between the two storage matrices, and max{·} means taking the maximum value between the two storage matrices.
Figure BDA00039840886300000613
is the first metric matrix with the largest storage in the i-th place,
Figure BDA00039840886300000614
is the second metric matrix with the ij,dth smallest storage,
Figure BDA00039840886300000615
is the second largest metric matrix stored in the i-th place,
Figure BDA00039840886300000616
is the ij,dth known left storage variable matrix,
Figure BDA00039840886300000617
is the ij,dth known right storage variable matrix,
Figure BDA00039840886300000618
is the third smallest metric matrix stored for the ijth time.
Figure BDA00039840886300000619
is the third largest metric matrix stored in the ijth order,
Figure BDA00039840886300000620
is the ijth known internal storage variable matrix,
Figure BDA00039840886300000621
is the ijth known external storage variable matrix, diag{·} is a diagonal matrix, A - is the first diagonal matrix defined, A + is the second diagonal matrix defined,
Figure BDA00039840886300000622
is the third diagonal matrix defined,
Figure BDA00039840886300000623
is the fourth diagonal matrix defined, B - is the fifth diagonal matrix defined, B + is the sixth diagonal matrix defined, and n is the dimension.

容易得出A(xk)∈[A-,A+]、

Figure BDA00039840886300000624
和B(xk)∈[B-,B+]。令
Figure BDA00039840886300000625
Figure BDA00039840886300000626
Figure BDA00039840886300000627
则有:It is easy to derive A(x k )∈[A - ,A + ],
Figure BDA00039840886300000624
and B(x k )∈[B ,B + ]. Let
Figure BDA00039840886300000625
Figure BDA00039840886300000626
and
Figure BDA00039840886300000627
Then we have:

Figure BDA00039840886300000628
Figure BDA00039840886300000628

式中,

Figure BDA00039840886300000629
为定义的左右区间的第一号矩阵,
Figure BDA00039840886300000630
为定义的左右区间的第二号矩阵,
Figure BDA0003984088630000071
为定义的左右区间的第三号矩阵,
Figure BDA0003984088630000072
Figure BDA0003984088630000073
满足范数有界不确定性:In the formula,
Figure BDA00039840886300000629
is the first matrix of the defined left and right intervals,
Figure BDA00039840886300000630
is the second matrix of the defined left and right intervals,
Figure BDA0003984088630000071
is the third matrix of the defined left and right intervals,
Figure BDA0003984088630000072
and
Figure BDA0003984088630000073
Satisfies norm-bounded uncertainty:

Figure BDA0003984088630000074
Figure BDA0003984088630000074

式中,ΔAk为满足范数有界不确定性的第一号矩阵,ΔAdk为满足范数有界不确定性的第二号矩阵,ΔBk为满足范数有界不确定性的第三号矩阵,

Figure BDA0003984088630000075
Figure BDA0003984088630000076
均为已知的实值权重矩阵,
Figure BDA0003984088630000077
是未知矩阵且满足
Figure BDA0003984088630000078
Where ΔA k is the first matrix satisfying the norm bounded uncertainty, ΔA dk is the second matrix satisfying the norm bounded uncertainty, and ΔB k is the third matrix satisfying the norm bounded uncertainty.
Figure BDA0003984088630000075
and
Figure BDA0003984088630000076
are all known real-valued weight matrices,
Figure BDA0003984088630000077
is an unknown matrix and satisfies
Figure BDA0003984088630000078

步骤二、放大转发协议下对步骤一建立的分数阶忆阻神经网络动态模型进行状态估计。具体步骤如下:Step 2: Perform state estimation on the dynamic model of the fractional-order memristor neural network established in step 1 under the amplification and forwarding protocol. The specific steps are as follows:

步骤二一、为了顺利完成远程数据传输的任务,在无线网络信道中安置一种放大-转发中继器以此来补充数据传输所消耗的能量。令ps,k和ns,k分别表示传感器和放大-转发中继器具有的随机能量,放大-转发中继器的输出信号由

Figure BDA0003984088630000079
表示,其满足如下方程:Step 21: In order to successfully complete the task of long-distance data transmission, an amplifier-forward repeater is placed in the wireless network channel to supplement the energy consumed by data transmission. Let ps,k and ns ,k represent the random energy of the sensor and the amplifier-forward repeater respectively. The output signal of the amplifier-forward repeater is given by
Figure BDA0003984088630000079
It means that it satisfies the following equation:

Figure BDA00039840886300000710
Figure BDA00039840886300000710

式中,

Figure BDA00039840886300000711
表示在第k时刻已知信道衰减矩阵,
Figure BDA00039840886300000712
为已知信道的衰减矩阵的m个信道的分量形式,diag{·}表示的是对角矩阵,yk是第k时刻的理想测量输出,
Figure BDA00039840886300000713
是第k时刻的实际测量输出,
Figure BDA00039840886300000722
是第k时刻在传感器-中继器信道的白噪声序列且满足
Figure BDA00039840886300000714
Figure BDA00039840886300000715
表示的是数学期望,
Figure BDA00039840886300000721
是在第k时刻
Figure BDA00039840886300000720
的转置,ps,k表示在第k时刻传感器具备的随机能量,满足如下统计特性:In the formula,
Figure BDA00039840886300000711
represents the known channel attenuation matrix at the kth time,
Figure BDA00039840886300000712
is the component form of the attenuation matrix of the known channel for m channels, diag{·} represents a diagonal matrix, y k is the ideal measurement output at the kth moment,
Figure BDA00039840886300000713
is the actual measured output at the kth moment,
Figure BDA00039840886300000722
is a white noise sequence in the sensor-repeater channel at the kth moment and satisfies
Figure BDA00039840886300000714
Figure BDA00039840886300000715
represents the mathematical expectation,
Figure BDA00039840886300000721
At the kth moment
Figure BDA00039840886300000720
The transpose of , ps,k represents the random energy of the sensor at the kth moment, satisfying the following statistical characteristics:

Figure BDA00039840886300000716
Figure BDA00039840886300000716

式中,Pr{·}表示的是数学概率,

Figure BDA00039840886300000717
表示所有概率的求和值为1,并且概率满足区间
Figure BDA00039840886300000718
Figure BDA00039840886300000719
为第k时刻的传感器具备的随机能量的期望值,φ表示的是所有信道的数量。In the formula, Pr{·} represents the mathematical probability,
Figure BDA00039840886300000717
Indicates that the sum of all probabilities is 1, and the probability satisfies the interval
Figure BDA00039840886300000718
Figure BDA00039840886300000719
is the expected value of the random energy of the sensor at the kth moment, and φ represents the number of all channels.

放大-转发中继器的输出值可表示为:The output value of the amplify-and-forward repeater can be expressed as:

Figure BDA0003984088630000081
Figure BDA0003984088630000081

式中,χk>0表示在第k时刻的放大系数,

Figure BDA0003984088630000082
是第k时刻已知信道的衰减矩阵,
Figure BDA0003984088630000083
为已知信道的衰减矩阵的m个信道的分量形式,m表示的是第m个信道,ns,k为第k时刻的传输随机能量的变量,
Figure BDA0003984088630000084
是第k时刻的实际测量输出,
Figure BDA00039840886300000819
是第k时刻在中继器-估计器信道的白噪声信号且满足
Figure BDA0003984088630000085
Figure BDA0003984088630000086
表示的是数学期望,
Figure BDA00039840886300000820
是在第k时刻
Figure BDA00039840886300000821
的转置。同样地,随机能量ns,k具有如下统计特性:In the formula, χ k >0 represents the amplification factor at the kth moment,
Figure BDA0003984088630000082
is the attenuation matrix of the known channel at the kth moment,
Figure BDA0003984088630000083
is the component form of the attenuation matrix of the known channel, m represents the mth channel, ns,k is the variable of the transmission random energy at the kth moment,
Figure BDA0003984088630000084
is the actual measured output at the kth moment,
Figure BDA00039840886300000819
is the white noise signal in the repeater-estimator channel at the kth moment and satisfies
Figure BDA0003984088630000085
Figure BDA0003984088630000086
represents the mathematical expectation,
Figure BDA00039840886300000820
At the kth moment
Figure BDA00039840886300000821
Similarly, the random energy n s,k has the following statistical properties:

Figure BDA0003984088630000087
Figure BDA0003984088630000087

式中,

Figure BDA0003984088630000088
表示所有概率的求和值为1,并且概率满足区间
Figure BDA0003984088630000089
Figure BDA00039840886300000810
为第k时刻的传输随机能量的期望值,ψ表示的是所有的信道的数量。In the formula,
Figure BDA0003984088630000088
Indicates that the sum of all probabilities is 1, and the probability satisfies the interval
Figure BDA0003984088630000089
Figure BDA00039840886300000810
is the expected value of the random energy transmitted at the kth moment, and ψ represents the number of all channels.

非线性函数f(s)满足如下扇形有界条件:The nonlinear function f(s) satisfies the following fan-shaped bounded conditions:

Figure BDA00039840886300000811
Figure BDA00039840886300000811

式中,

Figure BDA00039840886300000812
是第1个分量在k时刻的已知适当维数的第一号实矩阵,
Figure BDA00039840886300000813
是第2个分量的在k时刻的已知适当维数的第二号实矩阵。In the formula,
Figure BDA00039840886300000812
is the first real matrix of known appropriate dimension of the first component at time k,
Figure BDA00039840886300000813
is the second real matrix of known appropriate dimension of the second component at time k.

步骤二二、基于可获得的测量信息,构造如下的时变状态估计器:Step 22: Based on the available measurement information, construct the following time-varying state estimator:

Figure BDA00039840886300000814
Figure BDA00039840886300000814

式中,

Figure BDA00039840886300000815
是神经网络在第k时刻的状态估计,
Figure BDA00039840886300000816
是神经网络在第k时刻的状态估计,
Figure BDA00039840886300000817
是神经网络在第k-d时刻的状态估计,
Figure BDA00039840886300000818
为神经网络动态模型状态的实数域且其维数为n;χk表示在第k时刻的放大系数,d为一个固定的网络时滞,
Figure BDA0003984088630000091
为在第k时刻的被控输出的状态估计,
Figure BDA0003984088630000092
为神经网络动态模型被控输出状态的实数域且其维数为r,
Figure BDA0003984088630000093
为定义的左右区间的第一号矩阵,
Figure BDA0003984088630000094
为定义的左右区间的第二号矩阵,
Figure BDA0003984088630000095
为定义的左右区间的第三号矩阵,
Figure BDA0003984088630000096
为在第k时刻的非线性激励函数,Hk为在第k时刻的已知测量的调节矩阵,Dk是在第k时刻的已知测量的量度矩阵,
Figure BDA0003984088630000097
是第k时刻解码器的测量输出,Kk是在第k时刻估计器增益矩阵,
Figure BDA0003984088630000098
为传感器具备的随机能量期望的求和,
Figure BDA0003984088630000099
表示在第k时刻传感器具备的随机能量的期望,
Figure BDA00039840886300000910
为传感器具备的随机能量期望的求和,
Figure BDA00039840886300000911
表示在第k时刻传感器具备的随机能量的期望,
Figure BDA00039840886300000912
为所有的二项式组成的对角矩阵,
Figure BDA00039840886300000913
为分数阶(j=1,2,…,n),n为维数,diag{·}表示的是对角矩阵,χk表示在第k时刻的放大系数。In the formula,
Figure BDA00039840886300000815
is the state estimate of the neural network at the kth moment,
Figure BDA00039840886300000816
is the state estimate of the neural network at the kth moment,
Figure BDA00039840886300000817
is the state estimate of the neural network at the kdth moment,
Figure BDA00039840886300000818
is the real number domain of the dynamic model state of the neural network and its dimension is n; χ k represents the amplification factor at the kth moment, d is a fixed network delay,
Figure BDA0003984088630000091
is the state estimate of the controlled output at the kth moment,
Figure BDA0003984088630000092
is the real number domain of the controlled output state of the neural network dynamic model and its dimension is r,
Figure BDA0003984088630000093
is the first matrix of the defined left and right intervals,
Figure BDA0003984088630000094
is the second matrix of the defined left and right intervals,
Figure BDA0003984088630000095
is the third matrix of the defined left and right intervals,
Figure BDA0003984088630000096
is the nonlinear excitation function at the kth time, Hk is the adjustment matrix of the known measurement at the kth time, Dk is the measurement matrix of the known measurement at the kth time,
Figure BDA0003984088630000097
is the measured output of the decoder at time k, Kk is the estimator gain matrix at time k,
Figure BDA0003984088630000098
is the expected sum of random energies possessed by the sensor,
Figure BDA0003984088630000099
represents the expected random energy of the sensor at the kth moment,
Figure BDA00039840886300000910
is the expected sum of random energies possessed by the sensor,
Figure BDA00039840886300000911
represents the expected random energy of the sensor at the kth moment,
Figure BDA00039840886300000912
is the diagonal matrix of all binomials,
Figure BDA00039840886300000913
is a fractional order (j=1,2,…,n), n is the dimension, diag{·} represents a diagonal matrix, and χ k represents the amplification factor at the kth moment.

步骤二三、定义估计误差

Figure BDA00039840886300000914
和控制输出估计误差
Figure BDA00039840886300000915
进一步,可以得到估计误差系统:Step 2: Define the estimated error
Figure BDA00039840886300000914
and control output estimation error
Figure BDA00039840886300000915
Furthermore, we can get the estimated error system:

Figure BDA00039840886300000916
Figure BDA00039840886300000916

Figure BDA00039840886300000917
Figure BDA00039840886300000917

式中,

Figure BDA00039840886300000918
为第k时刻的激励函数,
Figure BDA00039840886300000919
为在第k时刻的非线性激励函数,
Figure BDA00039840886300000920
是神经网络在第k时刻的状态估计,
Figure BDA00039840886300000921
是神经网络在第k-d时刻的状态估计,
Figure BDA00039840886300000922
是神经网络在第k-ι+1时刻的状态估计,
Figure BDA00039840886300000923
为神经网络动态模型状态的实数域,n为维数,χk表示在第k时刻的放大系数,
Figure BDA00039840886300000924
表示的是开根号的值,Kk是在第k时刻估计器增益矩阵,
Figure BDA00039840886300000925
表示在第k时刻传感器具备的随机能量的期望,
Figure BDA0003984088630000101
表示在第k时刻传感器具备的随机能量的期望,ΔAk为满足范数有界不确定性的第一号矩阵,ΔAdk为满足范数有界不确定性的第二号矩阵,ΔBk为满足范数有界不确定性的第三号矩阵,
Figure BDA0003984088630000102
为定义的左右区间的第一号矩阵,
Figure BDA0003984088630000103
为定义的左右区间的第二号矩阵,
Figure BDA0003984088630000104
为定义的左右区间的第三号矩阵,ek是在第k时刻的估计误差,ek+1是在第k+1时刻的估计误差,ek-d是在第k-d时刻的估计误差,
Figure BDA0003984088630000105
是在第k时刻的被控输出估计误差,A(xk)=diagn{ai(xik)}为在第k时刻的神经网络自反馈对角矩阵,diag{·}表示的是对角矩阵,ai(xik)为A(xk)的第i个分量,n为维数;Ad(xk)为在第k时刻的已知维数且与时滞相关的系统矩阵,B(xk)为在第k时刻已知的连接激励函数的权重矩阵;f(xk)为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;v1k为在第k时刻均值为零并且协方差为V1>0的高斯白噪声序列,v2k为在第k时刻均值为零并且协方差为V2>0的高斯白噪声序列,
Figure BDA0003984088630000106
表示的是ι=1到k+1求和的值,
Figure BDA00039840886300001020
是第k时刻在传感器-中继器信道的白噪声序列且满足
Figure BDA0003984088630000107
Figure BDA0003984088630000108
表示的是数学期望,
Figure BDA00039840886300001021
是在第k时刻
Figure BDA00039840886300001016
的转置,
Figure BDA00039840886300001017
是第k时刻在中继器-估计器信道的白噪声信号且满足
Figure BDA0003984088630000109
Figure BDA00039840886300001010
表示的是数学期望,
Figure BDA00039840886300001018
是在第k时刻
Figure BDA00039840886300001019
的转置,
Figure BDA00039840886300001011
表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,
Figure BDA00039840886300001012
是第k时刻已知信道衰减矩阵,m表示的是第m个信道。In the formula,
Figure BDA00039840886300000918
is the activation function at the kth moment,
Figure BDA00039840886300000919
is the nonlinear activation function at the kth moment,
Figure BDA00039840886300000920
is the state estimate of the neural network at the kth moment,
Figure BDA00039840886300000921
is the state estimate of the neural network at the kdth moment,
Figure BDA00039840886300000922
is the state estimate of the neural network at the k-ι+1th moment,
Figure BDA00039840886300000923
is the real number domain of the dynamic model state of the neural network, n is the dimension, χ k represents the amplification factor at the kth moment,
Figure BDA00039840886300000924
represents the value of the square root, K k is the estimator gain matrix at the kth moment,
Figure BDA00039840886300000925
represents the expected random energy of the sensor at the kth moment,
Figure BDA0003984088630000101
represents the expectation of the random energy possessed by the sensor at the kth moment, ΔA k is the first matrix satisfying the norm bounded uncertainty, ΔA dk is the second matrix satisfying the norm bounded uncertainty, ΔB k is the third matrix satisfying the norm bounded uncertainty,
Figure BDA0003984088630000102
is the first matrix of the defined left and right intervals,
Figure BDA0003984088630000103
is the second matrix of the defined left and right intervals,
Figure BDA0003984088630000104
is the third matrix of the defined left and right intervals, e k is the estimated error at the kth moment, e k+1 is the estimated error at the k+1th moment, e kd is the estimated error at the kdth moment,
Figure BDA0003984088630000105
is the estimated error of the controlled output at the kth moment, A( xk ) = diagn {a i ( xik )} is the diagonal matrix of the neural network self-feedback at the kth moment, diag{·} represents a diagonal matrix, a i ( xik ) is the i-th component of A( xk ), and n is the dimension; Ad ( xk ) is the system matrix of known dimension and time delay related at the kth moment, B( xk ) is the weight matrix of the connection activation function known at the kth moment; f( xk ) is the nonlinear activation function at the kth moment; C1k is the noise distribution matrix of the system with the first component known at the kth moment, C2k is the noise distribution matrix of the system with the second component known at the kth moment, Hk is the adjustment matrix of the known measurement at the kth moment; Dk is the measurement matrix of the known measurement at the kth moment; v1k is a Gaussian white noise sequence with zero mean and covariance V1 >0 at the kth moment, and v2k is a Gaussian white noise sequence with zero mean and covariance V2 at the kth moment. >0 Gaussian white noise sequence,
Figure BDA0003984088630000106
It represents the sum of ι=1 to k+1.
Figure BDA00039840886300001020
is a white noise sequence in the sensor-repeater channel at the kth moment and satisfies
Figure BDA0003984088630000107
Figure BDA0003984088630000108
represents the mathematical expectation,
Figure BDA00039840886300001021
At the kth moment
Figure BDA00039840886300001016
The transpose of
Figure BDA00039840886300001017
is the white noise signal in the repeater-estimator channel at the kth moment and satisfies
Figure BDA0003984088630000109
Figure BDA00039840886300001010
represents the mathematical expectation,
Figure BDA00039840886300001018
At the kth moment
Figure BDA00039840886300001019
The transpose of
Figure BDA00039840886300001011
represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix,
Figure BDA00039840886300001012
is the known channel attenuation matrix at the kth moment, and m represents the mth channel.

本步骤的主要目的是设计一个基于放大转发协议下的时变状态估计器(2),使得估计误差系统同时满足以下两个性能约束要求:The main purpose of this step is to design a time-varying state estimator (2) based on the amplify-and-forward protocol so that the estimation error system satisfies the following two performance constraints at the same time:

(1)令扰动衰减水平γ>0,半正定矩阵一号和半正定矩阵二号分别为

Figure BDA00039840886300001013
Figure BDA00039840886300001014
对于初始状态e0,控制输出估计误差
Figure BDA00039840886300001015
满足如下的H性能约束条件:(1) Let the disturbance attenuation level γ>0, and the semi-positive definite matrix No. 1 and the semi-positive definite matrix No. 2 are respectively
Figure BDA00039840886300001013
and
Figure BDA00039840886300001014
For the initial state e 0 , the control output estimation error
Figure BDA00039840886300001015
The following H∞ performance constraints are met:

Figure BDA0003984088630000111
Figure BDA0003984088630000111

式中,N为有限的节点个数,

Figure BDA0003984088630000112
表示的是数学期望,
Figure BDA0003984088630000113
是第一号权重矩阵,
Figure BDA0003984088630000114
是第一号权重矩阵,e0是在第0时刻的估计误差,γ>0是给定的扰动衰减水平,
Figure BDA0003984088630000115
是噪声v1k和v2k增广的向量,
Figure BDA0003984088630000116
是在第k时刻ek的转置,·表示的是范数形式,·2表示的是范数平方的形式。Where N is the finite number of nodes.
Figure BDA0003984088630000112
represents the mathematical expectation,
Figure BDA0003984088630000113
is the first weight matrix,
Figure BDA0003984088630000114
is the first weight matrix, e 0 is the estimation error at time 0, γ>0 is the given disturbance attenuation level,
Figure BDA0003984088630000115
is the vector augmented by the noise v 1k and v 2k ,
Figure BDA0003984088630000116
is the transpose of e k at the kth moment, · represents the norm form, and · 2 represents the squared norm form.

(2)估计误差协方差满足如下的上界约束条件:(2) The estimated error covariance satisfies the following upper bound constraints:

Figure BDA0003984088630000117
Figure BDA0003984088630000117

式中,

Figure BDA0003984088630000118
是在第k时刻ek的转置,Πk(0≤k<N)是在第k时刻的一系列预先给定的可接受的估计精度矩阵。In the formula,
Figure BDA0003984088630000118
is the transpose of e k at the kth moment, and Π k (0≤k<N) is a series of pre-given acceptable estimation accuracy matrices at the kth moment.

步骤三、给定H性能指标γ、半正定矩阵一号

Figure BDA0003984088630000119
半正定矩阵二号
Figure BDA00039840886300001110
及初始条件
Figure BDA00039840886300001111
计算分数阶忆阻神经网络的误差协方差矩阵的上界及H性能约束条件。具体步骤如下:Step 3: Given the H∞ performance index γ and the semi-positive definite matrix No.
Figure BDA0003984088630000119
Semi-positive definite matrix II
Figure BDA00039840886300001110
and initial conditions
Figure BDA00039840886300001111
Calculate the upper bound of the error covariance matrix and H∞ performance constraints of the fractional-order memristor neural network. The specific steps are as follows:

步骤三一、按照下式证明出H性能分析问题并给出相应的易于求解的判别准则:Step 3: 1. Prove the H∞ performance analysis problem according to the following formula and give the corresponding easy-to-solve judgment criteria:

Figure BDA00039840886300001112
Figure BDA00039840886300001112

式中:Where:

Figure BDA0003984088630000121
Figure BDA0003984088630000121

Figure BDA0003984088630000122
Figure BDA0003984088630000122

Figure BDA0003984088630000123
Figure BDA0003984088630000123

Figure BDA0003984088630000124
Figure BDA0003984088630000124

Figure BDA0003984088630000125
Figure BDA0003984088630000125

Figure BDA0003984088630000126
Figure BDA0003984088630000126

Figure BDA0003984088630000127
Figure BDA0003984088630000127

Figure BDA0003984088630000128
Figure BDA0003984088630000128

Figure BDA0003984088630000129
Figure BDA0003984088630000129

Figure BDA00039840886300001210
Figure BDA00039840886300001210

Figure BDA00039840886300001211
Figure BDA00039840886300001211

Figure BDA00039840886300001212
Figure BDA00039840886300001212

Figure BDA00039840886300001213
Figure BDA00039840886300001213

Figure BDA00039840886300001214
Figure BDA00039840886300001214

式中,γ为给定的正标量;

Figure BDA00039840886300001215
为半正定矩阵一号,
Figure BDA00039840886300001216
Figure BDA00039840886300001217
分别为
Figure BDA00039840886300001218
Dk、Kk、Et,k、Ct,k、ΔAk、Hk
Figure BDA00039840886300001219
ΔBk
Figure BDA00039840886300001220
ΔAk
Figure BDA00039840886300001221
Ek、Kk、Ck、R3k的转置;
Figure BDA00039840886300001222
为半正定矩阵;Y11是Y的第1行第1列分块矩阵,Y12是Y的第1行第2列分块矩阵,Y22是Y的第2行第2列分块矩阵,Y33是Y的第3行第3列分块矩阵,Y44是Y的第4行第4列分块矩阵,Y55是Y的第5行第5列分块矩阵,Y66是Y的第6行第6列分块矩阵,Y77是Y的第7行第7列分块矩阵,Y88是Y的第8行第8列分块矩阵,Y99是Y的第9行第9列分块矩阵,
Figure BDA00039840886300001223
表示在第k时刻传感器具备的随机能量的期望,
Figure BDA00039840886300001224
表示在第k时刻传感器具备的随机能量的期望,ΔAk为满足范数有界不确定性的第一号矩阵,ΔAdk为满足范数有界不确定性的第二号矩阵,ΔBk为满足范数有界不确定性的第三号矩阵,
Figure BDA00039840886300001225
为定义的左右区间的第一号矩阵,
Figure BDA00039840886300001226
为定义的左右区间的第二号矩阵,
Figure BDA00039840886300001227
为定义的左右区间的第三号矩阵,
Figure BDA0003984088630000131
为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;
Figure BDA0003984088630000132
表示的是ι=1到k+1求和的值,
Figure BDA0003984088630000133
表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,
Figure BDA0003984088630000134
是第k时刻已知信道衰减矩阵,m表示的是第m个信道,
Figure BDA0003984088630000135
Figure BDA0003984088630000136
分别是第一个、第二个、第三个、第四个和第五个相关比例系数,0表示的是矩阵块的元素均为0。Where γ is a given positive scalar;
Figure BDA00039840886300001215
is a semi-positive definite matrix number one,
Figure BDA00039840886300001216
Figure BDA00039840886300001217
They are
Figure BDA00039840886300001218
D k , K k , E t,k , C t,k , ΔA k , H k ,
Figure BDA00039840886300001219
ΔB k ,
Figure BDA00039840886300001220
ΔA k ,
Figure BDA00039840886300001221
Transpose of E k , K k , C k , R 3k ;
Figure BDA00039840886300001222
is a semi-positive definite matrix; Y 11 is the 1st row and 1st column block matrix of Y, Y 12 is the 1st row and 2nd column block matrix of Y, Y 22 is the 2nd row and 2nd column block matrix of Y, Y 33 is the 3rd row and 3rd column block matrix of Y, Y 44 is the 4th row and 4th column block matrix of Y, Y 55 is the 5th row and 5th column block matrix of Y, Y 66 is the 6th row and 6th column block matrix of Y, Y 77 is the 7th row and 7th column block matrix of Y, Y 88 is the 8th row and 8th column block matrix of Y, Y 99 is the 9th row and 9th column block matrix of Y,
Figure BDA00039840886300001223
represents the expectation of the random energy possessed by the sensor at the kth moment,
Figure BDA00039840886300001224
represents the expectation of the random energy possessed by the sensor at the kth moment, ΔA k is the first matrix satisfying the norm bounded uncertainty, ΔA dk is the second matrix satisfying the norm bounded uncertainty, ΔB k is the third matrix satisfying the norm bounded uncertainty,
Figure BDA00039840886300001225
is the first matrix of the defined left and right intervals,
Figure BDA00039840886300001226
is the second matrix of the defined left and right intervals,
Figure BDA00039840886300001227
is the third matrix of the defined left and right intervals,
Figure BDA0003984088630000131
is the nonlinear excitation function at the kth moment; C 1k is the noise distribution matrix of the first component known system at the kth moment, C 2k is the noise distribution matrix of the second component known system at the kth moment, H k is the adjustment matrix of the known measurement at the kth moment; D k is the measurement matrix of the known measurement at the kth moment;
Figure BDA0003984088630000132
It represents the sum of ι=1 to k+1.
Figure BDA0003984088630000133
represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix,
Figure BDA0003984088630000134
is the known channel attenuation matrix at the kth moment, m represents the mth channel,
Figure BDA0003984088630000135
and
Figure BDA0003984088630000136
They are the first, second, third, fourth and fifth relevant proportional coefficients respectively, and 0 means that all elements of the matrix block are 0.

步骤三二、探讨协方差矩阵χk的上界约束问题,并给出如下充分条件:Step 3.2: Discuss the upper bound constraint of the covariance matrix χk and give the following sufficient conditions:

Sk+1≥Ω(Sk), (4)S k+1 ≥Ω(S k ), (4)

式中,In the formula,

Figure BDA0003984088630000137
Figure BDA0003984088630000137

Figure BDA0003984088630000138
Figure BDA0003984088630000138

式中,ek为在第k时刻的误差矩阵;

Figure BDA0003984088630000139
为在第k时刻的状态估计,ρ∈(0,1)为已知的调节正常数;Sk为在第k时刻的误差协方差矩阵的上界;
Figure BDA00039840886300001310
Θ1k T
Figure BDA00039840886300001311
Figure BDA00039840886300001312
分别为
Figure BDA00039840886300001313
Θ1k
Figure BDA00039840886300001314
C1k、Φι、Ct,k、Et,k的转置;ζ为调节系数,Ω(Sk)为在第k时刻求解出的上界矩阵;Sk-d为在第k-d时刻的误差协方差矩阵的上界矩阵;tr(Sk)为在第k时刻的误差协方差矩阵上界的迹;tr()为矩阵的迹,χk=ekek T为在第k时刻的误差上界,ek为在第k时刻的误差矩阵,I为单位矩阵,
Figure BDA0003984088630000141
是第1个分量在k时刻的已知适当维数的第一号实矩阵,
Figure BDA0003984088630000142
是第2个分量在k时刻的已知适当维数的第二号实矩阵。Where, e k is the error matrix at the kth moment;
Figure BDA0003984088630000139
is the state estimate at the kth moment, ρ∈(0,1) is a known normal constant; S k is the upper bound of the error covariance matrix at the kth moment;
Figure BDA00039840886300001310
Θ 1k T ,
Figure BDA00039840886300001311
Figure BDA00039840886300001312
They are
Figure BDA00039840886300001313
Θ 1k
Figure BDA00039840886300001314
is the transpose of C 1k , Φ ι , C t,k , and E t,k ; ζ is the adjustment coefficient, Ω(S k ) is the upper bound matrix solved at the kth moment; S kd is the upper bound matrix of the error covariance matrix at the kdth moment; tr(S k ) is the trace of the upper bound of the error covariance matrix at the kth moment; tr() is the trace of the matrix, χ k = ek e k T is the upper bound of the error at the kth moment, e k is the error matrix at the kth moment, I is the identity matrix,
Figure BDA0003984088630000141
is the first real matrix of known appropriate dimension of the first component at time k,
Figure BDA0003984088630000142
is the second real matrix of known appropriate dimension of the second component at time k.

通过对上述两个结果的分析,得到了保证估计误差系统满足给定的H性能要求和误差协方差有界性的充分条件。By analyzing the above two results, sufficient conditions are obtained to ensure that the estimation error system meets the given H performance requirements and the error covariance is bounded.

步骤四、利用随机分析方法,并通过解一系列线性矩阵不等式求解出估计器增益矩阵Kk的解,实现对放大转发协议下的分数阶忆阻神经网络动态模型进行状态估计;判断k+1是否达到总时长N,若k+1<N,则执行步骤二,反之结束。Step 4: Use the random analysis method and solve a series of linear matrix inequalities to solve the solution of the estimator gain matrix Kk to realize the state estimation of the dynamic model of the fractional-order memristor neural network under the amplification and forwarding protocol; determine whether k+1 reaches the total duration N. If k+1<N, execute step 2, otherwise end.

本步骤中,通过求解(5)~(7)一系列递推线性矩阵不等式,给出估计误差系统同时满足H性能要求和误差协方差有上界的充分条件,即可计算出估计器增益矩阵的值:In this step, by solving a series of recursive linear matrix inequalities (5) to (7), sufficient conditions are given for the estimation error system to simultaneously meet the H∞ performance requirements and the upper bound of the error covariance, and the value of the estimator gain matrix can be calculated:

Figure BDA0003984088630000143
Figure BDA0003984088630000143

Figure BDA0003984088630000144
Figure BDA0003984088630000144

Sk+1k+1≤0 (7)S k+1k+1 ≤0 (7)

更新矩阵为:The updated matrix is:

Figure BDA0003984088630000145
Figure BDA0003984088630000145

式中:Where:

Ω22=diag{-ε1,kI,-ε2,kI,-ε2,kI,-ε3,kI,-ε3,kI},Ω 22 =diag{-ε 1,k I,-ε 2,k I,-ε 2,k I,-ε 3,k I,-ε 3,k I},

Ω33=diag{-ε4,kI,-ε4,kI,-ε5,kI,-ε5,kI},Ω 33 =diag{-ε 4,k I,-ε 4,k I,-ε 5,k I,-ε 5,k I},

Figure BDA0003984088630000146
Figure BDA0003984088630000146

Figure BDA0003984088630000151
Figure BDA0003984088630000151

Figure BDA0003984088630000152
Figure BDA0003984088630000152

Figure BDA0003984088630000153
Figure BDA0003984088630000153

Figure BDA0003984088630000154
Figure BDA0003984088630000154

Figure BDA0003984088630000155
Figure BDA0003984088630000155

Figure BDA0003984088630000156
Figure BDA0003984088630000156

Figure BDA0003984088630000157
Figure BDA0003984088630000157

Figure BDA0003984088630000158
Figure BDA0003984088630000158

Figure BDA0003984088630000159
Figure BDA0003984088630000159

Figure BDA00039840886300001510
Figure BDA00039840886300001510

Figure BDA00039840886300001511
Figure BDA00039840886300001511

Figure BDA00039840886300001512
Figure BDA00039840886300001512

Figure BDA00039840886300001513
Figure BDA00039840886300001513

Figure BDA00039840886300001514
Figure BDA00039840886300001514

Figure BDA00039840886300001515
Figure BDA00039840886300001515

Figure BDA0003984088630000161
Figure BDA0003984088630000161

Figure BDA0003984088630000162
Figure BDA0003984088630000162

Figure BDA0003984088630000163
Figure BDA0003984088630000163

Figure BDA0003984088630000164
Figure BDA0003984088630000164

Figure BDA0003984088630000165
Figure BDA0003984088630000165

Figure BDA0003984088630000166
Figure BDA0003984088630000166

Figure BDA0003984088630000167
Figure BDA0003984088630000167

Figure BDA0003984088630000168
Figure BDA0003984088630000168

式中,Ω11是第1行第1列分块矩阵,Ω12是第1行第2列分块矩阵,Ω13是第1行第3列分块矩阵,Ω22是第2行第2列分块矩阵,Ω33是第3行第3列分块矩阵,

Figure BDA0003984088630000169
是第1行第1列分块矩阵,
Figure BDA00039840886300001610
是第1行第2列分块矩阵,
Figure BDA00039840886300001611
是第1行第3列分块矩阵,
Figure BDA00039840886300001612
是第1行第4列分块矩阵,L15是第1行第5列分块矩阵,L16是第1行第6列分块矩阵,L22是第2行第2列分块矩阵,L33是第3行第3列分块矩阵,L44是第4行第4列分块矩阵,L55是第5行第5列分块矩阵,L66是第6行第6列分块矩阵,
Figure BDA00039840886300001613
是第1行第1列分块矩阵,G12是第1行第2列分块矩阵,G14是第1行第4列分块矩阵,
Figure BDA00039840886300001614
是第1行第5列分块矩阵,G22是第2行第2列分块矩阵,G24是第2行第4列分块矩阵,
Figure BDA00039840886300001615
是第2行第6列分块矩阵,
Figure BDA00039840886300001616
是第2行第7列分块矩阵,
Figure BDA00039840886300001617
是第2行第8列分块矩阵,G33是第3行第3列分块矩阵,G39是第3行第9列分块矩阵,
Figure BDA00039840886300001618
是第4行第10列分块矩阵,
Figure BDA00039840886300001619
是第4行第4列分块矩阵,
Figure BDA00039840886300001620
是第5行第5列分块矩阵,
Figure BDA00039840886300001621
是第6行第6列分块矩阵,
Figure BDA00039840886300001622
是第7行第7列分块矩阵,
Figure BDA00039840886300001623
是第8行第8列分块矩阵,
Figure BDA0003984088630000171
是第9行第9列分块矩阵,
Figure BDA0003984088630000172
是第10行第10列分块矩阵,
Figure BDA0003984088630000173
Wherein, Ω 11 is the 1st row and 1st column block matrix, Ω 12 is the 1st row and 2nd column block matrix, Ω 13 is the 1st row and 3rd column block matrix, Ω 22 is the 2nd row and 2nd column block matrix, Ω 33 is the 3rd row and 3rd column block matrix,
Figure BDA0003984088630000169
is the 1st row and 1st column block matrix,
Figure BDA00039840886300001610
is a block matrix with row 1 and column 2,
Figure BDA00039840886300001611
is a block matrix with row 1 and column 3,
Figure BDA00039840886300001612
is the 1st row and 4th column block matrix, L 15 is the 1st row and 5th column block matrix, L 16 is the 1st row and 6th column block matrix, L 22 is the 2nd row and 2nd column block matrix, L 33 is the 3rd row and 3rd column block matrix, L 44 is the 4th row and 4th column block matrix, L 55 is the 5th row and 5th column block matrix, L 66 is the 6th row and 6th column block matrix,
Figure BDA00039840886300001613
is the 1st row and 1st column block matrix, G12 is the 1st row and 2nd column block matrix, G14 is the 1st row and 4th column block matrix,
Figure BDA00039840886300001614
is the 1st row and 5th column block matrix, G 22 is the 2nd row and 2nd column block matrix, G 24 is the 2nd row and 4th column block matrix,
Figure BDA00039840886300001615
is the 2nd row and 6th column block matrix,
Figure BDA00039840886300001616
is the 2nd row and 7th column block matrix,
Figure BDA00039840886300001617
is the 2nd row and 8th column block matrix, G 33 is the 3rd row and 3rd column block matrix, G 39 is the 3rd row and 9th column block matrix,
Figure BDA00039840886300001618
is the 4th row and 10th column block matrix,
Figure BDA00039840886300001619
is the 4th row and 4th column block matrix,
Figure BDA00039840886300001620
is the 5th row and 5th column block matrix,
Figure BDA00039840886300001621
is the 6th row and 6th column block matrix,
Figure BDA00039840886300001622
is the 7th row and 7th column block matrix,
Figure BDA00039840886300001623
is the 8th row and 8th column block matrix,
Figure BDA0003984088630000171
is the 9th row and 9th column block matrix,
Figure BDA0003984088630000172
is the 10th row and 10th column block matrix,
Figure BDA0003984088630000173

分别为

Figure BDA0003984088630000174
Dk,Kk,Et,k,Ct,k,ΔAk,Hk
Figure BDA0003984088630000175
ΔBk
Figure BDA0003984088630000176
ΔAk
Figure BDA0003984088630000177
Ek,Kk,Ck,R3k的转置,
Figure BDA0003984088630000178
为定义的左右区间的第一号矩阵,
Figure BDA0003984088630000179
为定义的左右区间的第二号矩阵,
Figure BDA00039840886300001710
为定义的左右区间的第三号矩阵,
Figure BDA00039840886300001711
为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;
Figure BDA00039840886300001712
表示的是ι=1到k+1求和的值,
Figure BDA00039840886300001713
表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,
Figure BDA00039840886300001714
是第k时刻已知信道衰减矩阵,m表示的是第m个信道,ρ∈(0,1)为已知的调节正常数;Sk为在第k时刻的误差协方差矩阵的上界;
Figure BDA00039840886300001715
Θ1k T
Figure BDA00039840886300001716
Figure BDA00039840886300001717
分别为
Figure BDA00039840886300001718
Θ1k
Figure BDA00039840886300001719
C1k,Φι,Ct,k,Et,k的转置;ζ为调节系数,Ω(Sk)为在第k时刻求解出的上界矩阵;Sk-d为在第k-d时刻的误差协方差矩阵的上界矩阵;tr(Sk)为在第k时刻的误差协方差矩阵上界的迹;tr()为矩阵的迹,I为单位矩阵;
Figure BDA00039840886300001720
为在第k时刻的第一号权重矩阵;
Figure BDA00039840886300001721
为在第k时刻的第二号权重矩阵;
Figure BDA00039840886300001722
为在第k时刻的第三号权重矩阵;
Figure BDA00039840886300001723
是在第k时刻R3k的转置;
Figure BDA00039840886300001724
是第1个分量在k时刻的已知适当维数的第一号实矩阵,
Figure BDA00039840886300001725
是第2个分量在k时刻的已知适当维数的第二号实矩阵;
Figure BDA00039840886300001726
为在第k时刻的非线性激励函数的状态估计;
Figure BDA00039840886300001727
是第1个分量在k时刻的已知适当维数的第一号量度矩阵;
Figure BDA00039840886300001728
是第2个分量在k时刻的已知适当维数的第二号量度矩阵;
Figure BDA00039840886300001729
是第3个分量在k时刻的已知适当维数的第三号量度矩阵;
Figure BDA00039840886300001730
是第4个分量在k时刻的已知适当维数的第三号量度矩阵;N5是第5个分量在k时刻的已知适当维数的第三号量度矩阵;M1,M2,M3,M4和M5分别是第一号,第二号,第三号,第四号和第五号的量度矩阵,
Figure BDA0003984088630000181
为在第k时刻的神经元状态估计,
Figure BDA0003984088630000182
为第k时刻的半正定矩阵;
Figure BDA0003984088630000183
为第k时刻的半正定矩阵;
Figure BDA0003984088630000184
为第k-d时刻的半正定矩阵;
Figure BDA0003984088630000185
为在第k+1时刻的第一更新矩阵,Sk为估计误差的上界矩阵,tr(Sk)为在第k时刻估计误差上界矩阵Sk的迹;Sk-d为在k-d时刻的上界矩阵,κ为调节的权重系数,
Figure BDA0003984088630000186
Figure BDA0003984088630000187
均为已知的实值权重矩阵,
Figure BDA0003984088630000188
是未知矩阵且满足
Figure BDA0003984088630000189
Figure BDA00039840886300001810
Figure BDA00039840886300001811
的转置,γ为给定的正标量;
Figure BDA00039840886300001812
为给定的半正定矩阵一号;
Figure BDA00039840886300001813
分别是Ω12,Ω13
Figure BDA00039840886300001814
的转置;
Figure BDA00039840886300001815
Figure BDA00039840886300001816
分别是,G12,G14
Figure BDA00039840886300001817
G24,G410
Figure BDA00039840886300001818
的转置;
Figure BDA00039840886300001819
分别是M1,M2,M3,M4,M5的转置;N1,N2,N3,N4,N5分别是
Figure BDA00039840886300001820
的转置;
Figure BDA00039840886300001821
Figure BDA00039840886300001822
分别是第一个、第二个、第三个、第四个和第五个相关比例系数,0表示的是矩阵块的元素均为0。They are
Figure BDA0003984088630000174
D k , K k , E t,k , C t,k , ΔA k , H k ,
Figure BDA0003984088630000175
ΔB k
Figure BDA0003984088630000176
ΔA k
Figure BDA0003984088630000177
The transpose of E k , K k , C k , R 3k ,
Figure BDA0003984088630000178
is the first matrix of the defined left and right intervals,
Figure BDA0003984088630000179
is the second matrix of the defined left and right intervals,
Figure BDA00039840886300001710
is the third matrix of the defined left and right intervals,
Figure BDA00039840886300001711
is the nonlinear excitation function at the kth moment; C 1k is the noise distribution matrix of the first component known system at the kth moment, C 2k is the noise distribution matrix of the second component known system at the kth moment, H k is the adjustment matrix of the known measurement at the kth moment; D k is the measurement matrix of the known measurement at the kth moment;
Figure BDA00039840886300001712
It represents the sum of ι=1 to k+1.
Figure BDA00039840886300001713
represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix,
Figure BDA00039840886300001714
is the known channel attenuation matrix at the kth moment, m represents the mth channel, ρ∈(0,1) is the known normal constant; S k is the upper bound of the error covariance matrix at the kth moment;
Figure BDA00039840886300001715
Θ 1k T ,
Figure BDA00039840886300001716
Figure BDA00039840886300001717
They are
Figure BDA00039840886300001718
Θ 1k ,
Figure BDA00039840886300001719
C 1k , Φ ι , C t,k , the transpose of E t,k ; ζ is the adjustment coefficient, Ω(S k ) is the upper bound matrix solved at the kth time; S kd is the upper bound matrix of the error covariance matrix at the kdth time; tr(S k ) is the trace of the upper bound of the error covariance matrix at the kth time; tr() is the trace of the matrix, and I is the identity matrix;
Figure BDA00039840886300001720
is the first weight matrix at the kth moment;
Figure BDA00039840886300001721
is the second weight matrix at the kth moment;
Figure BDA00039840886300001722
is the third weight matrix at the kth moment;
Figure BDA00039840886300001723
is the transpose of R 3k at the kth moment;
Figure BDA00039840886300001724
is the first real matrix of known appropriate dimension of the first component at time k,
Figure BDA00039840886300001725
is the second real matrix of known appropriate dimension of the second component at time k;
Figure BDA00039840886300001726
is the state estimate of the nonlinear activation function at the kth moment;
Figure BDA00039840886300001727
is the first metric matrix of known appropriate dimension of the first component at time k;
Figure BDA00039840886300001728
is the second metric matrix of known appropriate dimension of the second component at time k;
Figure BDA00039840886300001729
is the third metric matrix of known appropriate dimension of the third component at time k;
Figure BDA00039840886300001730
is the third metric matrix of known appropriate dimension of the 4th component at time k; N5 is the third metric matrix of known appropriate dimension of the 5th component at time k; M1 , M2 , M3 , M4 and M5 are the first, second, third, fourth and fifth metric matrices respectively.
Figure BDA0003984088630000181
is the estimated neuron state at the kth moment,
Figure BDA0003984088630000182
is the semi-positive definite matrix at the kth moment;
Figure BDA0003984088630000183
is the semi-positive definite matrix at the kth moment;
Figure BDA0003984088630000184
is the semi-positive definite matrix at the kdth moment;
Figure BDA0003984088630000185
is the first update matrix at the k+1th time, Sk is the upper bound matrix of the estimation error, tr( Sk ) is the trace of the upper bound matrix Sk of the estimation error at the kth time; Skd is the upper bound matrix at the kdth time, κ is the adjusted weight coefficient,
Figure BDA0003984088630000186
and
Figure BDA0003984088630000187
are all known real-valued weight matrices,
Figure BDA0003984088630000188
is an unknown matrix and satisfies
Figure BDA0003984088630000189
Figure BDA00039840886300001810
yes
Figure BDA00039840886300001811
The transpose of , γ is a given positive scalar;
Figure BDA00039840886300001812
For the given semi-positive definite matrix one;
Figure BDA00039840886300001813
They are Ω 12 , Ω 13 ,
Figure BDA00039840886300001814
The transpose of
Figure BDA00039840886300001815
Figure BDA00039840886300001816
They are G12 , G14 ,
Figure BDA00039840886300001817
G24 , G410 ,
Figure BDA00039840886300001818
The transpose of
Figure BDA00039840886300001819
are the transposes of M 1 , M 2 , M 3 , M 4 , and M 5 respectively; N 1 , N 2 , N 3 , N 4 , and N 5 are
Figure BDA00039840886300001820
The transpose of
Figure BDA00039840886300001821
and
Figure BDA00039840886300001822
They are the first, second, third, fourth and fifth relevant proportional coefficients respectively, and 0 means that all elements of the matrix block are 0.

本发明中,步骤三与步骤四中所述理论为:In the present invention, the theory described in step 3 and step 4 is:

首先,证明出H性能分析问题并给出相应的易于求解的判别准则;其次,探讨协方差矩阵Xk的上界约束问题,并给出如下充分条件;通过对上述两个结果的分析,得到了保证估计误差系统满足给定的H性能要求和误差协方差有界性的充分条件,通过解一系列线性矩阵不等式求解出估计器增益矩阵的解,并计算出估计器增益矩阵Kk的解。Firstly, the H∞ performance analysis problem is proved and the corresponding easy-to-solve judgment criterion is given; secondly, the upper bound constraint problem of the covariance matrix Xk is discussed, and the following sufficient conditions are given; through the analysis of the above two results, sufficient conditions are obtained to ensure that the estimation error system meets the given H∞ performance requirements and the boundedness of the error covariance. The solution of the estimator gain matrix is solved by solving a series of linear matrix inequalities, and the solution of the estimator gain matrix Kk is calculated.

实施例:Example:

本实施例以具有H性能约束及方差约束的分数阶忆阻神经网络为例,此外还可以应用在联想记忆、模式识别和组合优化中,采用本发明所述方法针对语音识别案例进行仿真:This embodiment takes a fractional-order memristor neural network with H∞ performance constraints and variance constraints as an example. In addition, it can also be applied to associative memory, pattern recognition and combinatorial optimization. The method of the present invention is used to simulate a speech recognition case:

放大转发协议下具有H性能约束及方差约束的分数阶忆阻神经网络状态模型、测量输出模型及被控输出模型的相关系统参数选取如下:The relevant system parameters of the fractional-order memristor neural network state model, measurement output model and controlled output model under the amplification and forwarding protocol with H∞ performance constraints and variance constraints are selected as follows:

根据人的声音的状态给定相应的调节矩阵为:According to the state of human voice, the corresponding adjustment matrix is given as:

Figure BDA0003984088630000191
Figure BDA0003984088630000191

Figure BDA0003984088630000192
Figure BDA0003984088630000192

Figure BDA0003984088630000193
Figure BDA0003984088630000193

Figure BDA0003984088630000194
Figure BDA0003984088630000194

Figure BDA0003984088630000195
Figure BDA0003984088630000195

Figure BDA0003984088630000196
Figure BDA0003984088630000196

C1k=[-1.2-0.35sin(2k)]T,C 1k =[-1.2-0.35sin(2k)] T ,

测量调节矩阵为:The measurement adjustment matrix is:

C2k=[-0.2-0.1sin(3k)]T,

Figure BDA0003984088630000197
C 2k = [-0.2-0.1sin(3k)] T ,
Figure BDA0003984088630000197

被控输出调节矩阵为:The controlled output adjustment matrix is:

Hk=[-0.01-0.01sin(2k)]H k =[-0.01-0.01sin(2k)]

状态权重矩阵为:The state weight matrix is:

Figure BDA0003984088630000198
Figure BDA0003984088630000198

Figure BDA0003984088630000199
Figure BDA0003984088630000199

非线性函数的权重矩阵和调节参数为:The weight matrix and adjustment parameters of the nonlinear function are:

Figure BDA00039840886300001910
Figure BDA00039840886300001910

Case I:给出以下传输功率的概率分布:Case I: Given the following probability distribution of transmission power:

Prob{pt,k=1}=0.1,Prob{pt,k=1.5}=0.3,Prob{p t, k = 1} = 0.1, Prob {p t, k = 1.5} = 0.3,

Prob{pt,k=2}=0.6,Prob{nt,k=1}=0.2,Prob{p t, k = 2} = 0.6, Prob {n t, k = 1} = 0.2,

Prob{nt,k=1.5}=0.4,Prob{nt,k=2}=0.4,Prob{n t,k =1.5}=0.4,Prob{n t,k =2}=0.4,

Case II:给出以下传输功率的概率分布:Case II: Given the following probability distribution of transmission power:

Prob{pt,k=1}=0.6,Prob{pt,k=1.5}=0.3,Prob{p t, k = 1} = 0.6, Prob {p t, k = 1.5} = 0.3,

Prob{pt,k=2}=0.1,Prob{nt,k=1}=0.4,Prob{p t, k = 2} = 0.1, Prob {n t, k = 1} = 0.4,

Prob{nt,k=1.5}=0.4,Prob{nt,k=2}=0.2.Prob{n t,k =1.5}=0.4, Prob{n t,k =2}=0.2.

激励函数取为:The activation function is taken as:

Figure BDA0003984088630000201
Figure BDA0003984088630000201

式中,xk=[x1,kx2,k]T是神经元的状态向量,放大系数为χs=1,x1,k为在第k时刻xk的第一个分量比重矩阵,x2,k为在第k时刻xk的第二个分量比重矩阵。Where xk = [ x1, k x2 , k ] T is the state vector of the neuron, the amplification factor is χs = 1, x1, k is the first component weight matrix of xk at the kth moment, and x2 , k is the second component weight matrix of xk at the kth moment.

其它仿真初始值选取如下:Other simulation initial values are selected as follows:

扰动衰减水平γ=0.7,半正定矩阵一号

Figure BDA0003984088630000202
上界矩阵{Ωk}1≤k≤N=diag{0.2,0.2}和协方差
Figure BDA0003984088630000203
初始状态
Figure BDA0003984088630000204
通道参数Ct,s=0.38和Et,s=0.12,分别取传感器到中继信道和中继到估计器信道的噪声协方差
Figure BDA0003984088630000205
Figure BDA0003984088630000206
Perturbation attenuation level γ = 0.7, semi-positive definite matrix No.
Figure BDA0003984088630000202
The upper bound matrix {Ω k } 1≤k≤N = diag{0.2,0.2} and the covariance
Figure BDA0003984088630000203
Initial state
Figure BDA0003984088630000204
The channel parameters C t,s = 0.38 and E t,s = 0.12, which are the noise covariances of the sensor-to-relay channel and the relay-to-estimator channel, respectively.
Figure BDA0003984088630000205
and
Figure BDA0003984088630000206

利用递推线性矩阵不等式,求解线性矩阵不等式(5)~(7),部分数值如下:情形一(Case I):Using the recursive linear matrix inequality, we solve linear matrix inequalities (5) to (7). Some numerical values are as follows: Case I:

Figure BDA0003984088630000207
Figure BDA0003984088630000207

Figure BDA0003984088630000208
Figure BDA0003984088630000208

Figure BDA0003984088630000209
Figure BDA0003984088630000209

情形二(Case II):Case II:

Figure BDA00039840886300002010
Figure BDA00039840886300002010

Figure BDA0003984088630000211
Figure BDA0003984088630000211

Figure BDA0003984088630000212
Figure BDA0003984088630000212

状态估计器效果:State Estimator Effects:

由图2可见,针对放大转发协议下具有H性能约束及方差约束的分数阶忆阻神经网络,所发明的状态估计器设计方法可有效地估计出目标状态。As can be seen from FIG2 , for the fractional-order memristor neural network with H performance constraint and variance constraint under the amplify-and-forward protocol, the invented state estimator design method can effectively estimate the target state.

由图3、图4、图5可以看出,针对每个时刻,随着功率的减少,估计误差效果变差。It can be seen from FIG3 , FIG4 , and FIG5 that, at each moment, as the power decreases, the estimation error effect becomes worse.

Claims (9)

1.一种方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述方法包括如下步骤:1. A variance-constrained fractional-order memristor neural network estimation method, characterized in that the method comprises the following steps: 步骤一、建立放大转发协议下的分数阶忆阻神经网络动态模型;Step 1: Establish a dynamic model of fractional-order memristor neural network under amplification and forwarding protocol; 步骤二、放大转发协议下对步骤一建立的分数阶忆阻神经网络动态模型进行状态估计;Step 2: performing state estimation on the dynamic model of the fractional-order memristor neural network established in step 1 under the amplification and forwarding protocol; 步骤三、给定H性能指标γ、半正定矩阵一号
Figure FDA0003984088620000011
半正定矩阵二号
Figure FDA0003984088620000012
及初始条件
Figure FDA0003984088620000013
计算分数阶忆阻神经网络的误差协方差矩阵的上界及H性能约束条件;
Step 3: Given the H∞ performance index γ and the semi-positive definite matrix No.
Figure FDA0003984088620000011
Semi-positive definite matrix II
Figure FDA0003984088620000012
and initial conditions
Figure FDA0003984088620000013
Calculate the upper bound of the error covariance matrix and H∞ performance constraints of fractional-order memristor neural networks;
步骤四、利用随机分析方法,并通过解线性矩阵不等式求解出估计器增益矩阵Kk的解,实现对放大转发协议下的分数阶忆阻神经网络动态模型的状态估计,判断k+1是否达到总时长N,若k+1<N,则执行步骤二,反之结束。Step 4: Use the random analysis method and solve the linear matrix inequality to solve the solution of the estimator gain matrix Kk to realize the state estimation of the dynamic model of the fractional-order memristor neural network under the amplification and forwarding protocol, and judge whether k+1 reaches the total duration N. If k+1<N, execute step 2, otherwise end.
2.根据权利要求1所述的方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述步骤一中,根据Grunwald-Letnikov分数阶导数定义,分数阶忆阻神经网络动态模型的状态空间形式为:2. The variance-constrained fractional-order memristor neural network estimation method according to claim 1, characterized in that in the step 1, according to the Grunwald-Letnikov fractional-order derivative definition, the state space form of the fractional-order memristor neural network dynamic model is:
Figure FDA0003984088620000014
Figure FDA0003984088620000014
式中:Where:
Figure FDA0003984088620000015
Figure FDA0003984088620000015
Figure FDA0003984088620000016
Figure FDA0003984088620000016
Figure FDA0003984088620000017
Figure FDA0003984088620000017
这里,
Figure FDA0003984088620000021
表示微分算子,
Figure FDA0003984088620000022
为分数阶(j=1,2,…,n),n为维数,
Figure FDA0003984088620000023
是在第k时刻的分数阶忆阻神经网络的状态向量,
Figure FDA0003984088620000024
是在第k-ι+1时刻的分数阶忆阻神经网络的状态向量,
Figure FDA0003984088620000025
是在第k-d时刻的分数阶忆阻神经网络的状态向量,
Figure FDA0003984088620000026
是在第k+1时刻的分数阶忆阻神经网络的状态向量,
Figure FDA0003984088620000027
为神经网络动态模型状态的实数域且其维数为n;
Figure FDA0003984088620000028
为在第k时刻的被控测量输出,
Figure FDA0003984088620000029
为神经网络动态模型被控输出状态的实数域且其维数为r;
Figure FDA00039840886200000210
是给定的初始序列,d为离散固定的网络时滞;A(xk)=diagn{ai(xi,k)}为在第k时刻的神经网络自反馈对角矩阵,n为维数,diag{·}表示的是对角矩阵,ai(xi,k)为A(xk)的第i个分量,n为维数;Ad(xk)={aij,d(xi,k)}n*n为在第k时刻的已知维数且与时滞相关的系统矩阵,aij,d(xi,k)为在第k时刻Ad(xk)的第i个分量形式;B(xk)={bij(xi,k)}n*n为在第k时刻的已知的连接激励函数的权重矩阵,bij(xi,k)为在第k时刻B(xk)的第i个分量形式;f(xk)为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;v1k为在第k时刻均值为零并且协方差为V1>0的高斯白噪声序列,v2k为在第k时刻均值为零并且协方差为V2>0的高斯白噪声序列,
Figure FDA00039840886200000211
表示的是ι=1到k+1求和的值。
here,
Figure FDA0003984088620000021
represents the differential operator,
Figure FDA0003984088620000022
is a fractional order (j=1,2,…,n), n is the dimension,
Figure FDA0003984088620000023
is the state vector of the fractional-order memristor neural network at the kth moment,
Figure FDA0003984088620000024
is the state vector of the fractional-order memristor neural network at the k-ι+1th moment,
Figure FDA0003984088620000025
is the state vector of the fractional-order memristor neural network at the kdth moment,
Figure FDA0003984088620000026
is the state vector of the fractional-order memristor neural network at the k+1th time,
Figure FDA0003984088620000027
is the real number field of the states of the neural network dynamic model and its dimension is n;
Figure FDA0003984088620000028
is the controlled measured output at the kth moment,
Figure FDA0003984088620000029
is the real number domain of the controlled output state of the neural network dynamic model and its dimension is r;
Figure FDA00039840886200000210
is a given initial sequence, d is a discrete fixed network time delay; A(x k ) = diag n {a i (xi ,k )} is the neural network self-feedback diagonal matrix at the kth moment, n is the dimension, diag{·} represents a diagonal matrix, a i (xi ,k ) is the i-th component of A(x k ), and n is the dimension; A d (x k ) = {a ij,d (xi ,k )} n*n is the system matrix of known dimension and time delay at the kth moment, a ij,d (xi ,k ) is the i-th component form of A d (x k ) at the kth moment; B(x k ) = {b ij (xi ,k )} n*n is the weight matrix of the known connection activation function at the kth moment, b ij ( xi,k ) is the i-th component form of B(x k ) at the kth moment; f(x k ) is the nonlinear activation function at the kth moment; C C 1k is the noise distribution matrix of the first component known system at the kth moment, C 2k is the noise distribution matrix of the second component known system at the kth moment, H k is the adjustment matrix of the known measurement at the kth moment; D k is the measurement matrix of the known measurement at the kth moment; v 1k is a Gaussian white noise sequence with zero mean and covariance V 1 >0 at the kth moment, v 2k is a Gaussian white noise sequence with zero mean and covariance V 2 >0 at the kth moment,
Figure FDA00039840886200000211
It represents the sum of the values from ι=1 to k+1.
3.根据权利要求2所述的方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述ai(xi,k)、aij,d(xi,k)和bij(xi,k)满足:3. The variance-constrained fractional-order memristor neural network estimation method according to claim 2, wherein the a i (xi ,k ), a ij,d (xi ,k ) and b ij (xi ,k ) satisfy:
Figure FDA00039840886200000212
Figure FDA00039840886200000212
式中,ai(xi,k)、aij,d(xi,k)和bij(xi,k)分别为A(xk),Ad(xk)和B(xk)的第i个分量,Ωi>0为已知的切换阈值,
Figure FDA00039840886200000213
为第i个已知的上存储变量矩阵,
Figure FDA00039840886200000214
为第i个已知的下存储变量矩阵,
Figure FDA00039840886200000215
为第ij,d个已知的左存储变量矩阵,
Figure FDA00039840886200000216
为第ij,d个已知的右存储变量矩阵,
Figure FDA0003984088620000031
为第ij个已知的内存储变量矩阵,
Figure FDA0003984088620000032
为第ij个已知的外存储变量矩阵。
Where ai (xi ,k ), aij ,d (xi ,k ) and bij (xi ,k ) are the i-th components of A( xk ), Ad ( xk ) and B( xk ), respectively, Ωi >0 is the known switching threshold,
Figure FDA00039840886200000213
is the i-th known upper storage variable matrix,
Figure FDA00039840886200000214
is the i-th known storage variable matrix,
Figure FDA00039840886200000215
is the ij,dth known left storage variable matrix,
Figure FDA00039840886200000216
is the ij,dth known right storage variable matrix,
Figure FDA0003984088620000031
is the ijth known internal storage variable matrix,
Figure FDA0003984088620000032
is the ijth known external storage variable matrix.
4.根据权利要求1所述的方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述步骤二的具体步骤如下:4. The variance-constrained fractional-order memristor neural network estimation method according to claim 1, wherein the specific steps of step 2 are as follows: 步骤二一、令ps,k和ns,k分别表示传感器和放大-转发中继器具有的随机能量,放大-转发中继器的输出信号由
Figure FDA0003984088620000033
表示,其满足如下方程:
Step 21: Let ps,k and ns,k represent the random energy of the sensor and the amplifier-forward repeater respectively. The output signal of the amplifier-forward repeater is given by
Figure FDA0003984088620000033
It means that it satisfies the following equation:
Figure FDA0003984088620000034
Figure FDA0003984088620000034
式中,
Figure FDA0003984088620000035
表示在第k时刻已知信道衰减矩阵,
Figure FDA0003984088620000036
为已知信道的衰减矩阵的m个信道的分量形式,diag{·}表示的是对角矩阵,yk是第k时刻的理想测量输出,
Figure FDA0003984088620000037
是第k时刻的实际测量输出,θs1,k是第k时刻在传感器-中继器信道的白噪声序列且满足
Figure FDA0003984088620000038
Figure FDA0003984088620000039
表示的是数学期望,(θs1,k)T是在第k时刻θs1,k的转置,ps,k表示在第k时刻传感器具备的随机能量;
In the formula,
Figure FDA0003984088620000035
represents the known channel attenuation matrix at the kth time,
Figure FDA0003984088620000036
is the component form of the attenuation matrix of the known channel for m channels, diag{·} represents a diagonal matrix, y k is the ideal measurement output at the kth moment,
Figure FDA0003984088620000037
is the actual measured output at the kth moment, θ s1,k is the white noise sequence in the sensor-repeater channel at the kth moment and satisfies
Figure FDA0003984088620000038
Figure FDA0003984088620000039
represents the mathematical expectation, (θ s1,k ) T is the transpose of θ s1,k at the kth moment, and p s,k represents the random energy possessed by the sensor at the kth moment;
放大-转发中继器的输出值表示为:The output value of the amplify-and-forward repeater is expressed as:
Figure FDA00039840886200000310
Figure FDA00039840886200000310
式中,χk>0表示在第k时刻的放大系数,
Figure FDA00039840886200000311
是第k时刻已知信道的衰减矩阵,
Figure FDA00039840886200000312
为已知信道的衰减矩阵的m个信道的分量形式,m表示的是第m个信道,ns,k为第k时刻的传输随机能量的变量,
Figure FDA00039840886200000313
是第k时刻的实际测量输出,θs2,k是第k时刻在中继器-估计器信道的白噪声信号且满足
Figure FDA00039840886200000314
Figure FDA00039840886200000315
表示的是数学期望,(θs2,k)T是在第k时刻θs2,k的转置;
In the formula, χ k >0 represents the amplification factor at the kth moment,
Figure FDA00039840886200000311
is the attenuation matrix of the known channel at the kth moment,
Figure FDA00039840886200000312
is the component form of the attenuation matrix of the known channel, m represents the mth channel, ns,k is the variable of the transmission random energy at the kth moment,
Figure FDA00039840886200000313
is the actual measured output at the kth moment, θ s2,k is the white noise signal in the repeater-estimator channel at the kth moment and satisfies
Figure FDA00039840886200000314
Figure FDA00039840886200000315
represents the mathematical expectation, (θ s2,k ) T is the transpose of θ s2,k at the kth moment;
步骤二二、基于可获得的测量信息,构造如下的时变状态估计器:Step 22: Based on the available measurement information, construct the following time-varying state estimator:
Figure FDA00039840886200000316
Figure FDA00039840886200000316
式中,
Figure FDA00039840886200000317
是神经网络在第k时刻的状态估计,
Figure FDA00039840886200000318
是神经网络在第k时刻的状态估计,
Figure FDA00039840886200000319
是神经网络在第k-d时刻的状态估计,
Figure FDA00039840886200000320
为神经网络动态模型状态的实数域且其维数为n;χk表示在第k时刻的放大系数,d为一个固定的网络时滞,
Figure FDA0003984088620000041
为在第k时刻的被控输出的状态估计,
Figure FDA0003984088620000042
为神经网络动态模型被控输出状态的实数域且其维数为r,
Figure FDA0003984088620000043
为定义的左右区间的第一号矩阵,
Figure FDA0003984088620000044
为定义的左右区间的第二号矩阵,
Figure FDA0003984088620000045
为定义的左右区间的第三号矩阵,
Figure FDA0003984088620000046
为在第k时刻的非线性激励函数,Hk为在第k时刻的已知测量的调节矩阵,Dk是在第k时刻的已知测量的量度矩阵,
Figure FDA0003984088620000047
是第k时刻解码器的测量输出,Kk是在第k时刻估计器增益矩阵,
Figure FDA0003984088620000048
为传感器具备的随机能量期望的求和,
Figure FDA0003984088620000049
表示在第k时刻传感器具备的随机能量的期望,
Figure FDA00039840886200000410
为传感器具备的随机能量期望的求和,
Figure FDA00039840886200000411
表示在第k时刻传感器具备的随机能量的期望,
Figure FDA00039840886200000412
为所有的二项式组成的对角矩阵,
Figure FDA00039840886200000413
为分数阶(j=1,2,…,n),n为维数,diag{·}表示的是对角矩阵,χk表示在第k时刻的放大系数;
In the formula,
Figure FDA00039840886200000317
is the state estimate of the neural network at the kth moment,
Figure FDA00039840886200000318
is the state estimate of the neural network at the kth moment,
Figure FDA00039840886200000319
is the state estimate of the neural network at the kdth moment,
Figure FDA00039840886200000320
is the real number domain of the dynamic model state of the neural network and its dimension is n; χ k represents the amplification factor at the kth moment, d is a fixed network delay,
Figure FDA0003984088620000041
is the state estimate of the controlled output at the kth moment,
Figure FDA0003984088620000042
is the real number domain of the controlled output state of the neural network dynamic model and its dimension is r,
Figure FDA0003984088620000043
is the first matrix of the defined left and right intervals,
Figure FDA0003984088620000044
is the second matrix of the defined left and right intervals,
Figure FDA0003984088620000045
is the third matrix of the defined left and right intervals,
Figure FDA0003984088620000046
is the nonlinear excitation function at the kth time, Hk is the adjustment matrix of the known measurement at the kth time, Dk is the measurement matrix of the known measurement at the kth time,
Figure FDA0003984088620000047
is the measured output of the decoder at time k, Kk is the estimator gain matrix at time k,
Figure FDA0003984088620000048
is the expected sum of random energies possessed by the sensor,
Figure FDA0003984088620000049
represents the expected random energy of the sensor at the kth moment,
Figure FDA00039840886200000410
is the expected sum of random energies possessed by the sensor,
Figure FDA00039840886200000411
represents the expected random energy of the sensor at the kth moment,
Figure FDA00039840886200000412
is the diagonal matrix of all binomials,
Figure FDA00039840886200000413
is a fractional order (j=1,2,…,n), n is the dimension, diag{·} represents a diagonal matrix, and χ k represents the magnification factor at the kth moment;
步骤二三、定义估计误差
Figure FDA00039840886200000414
和控制输出估计误差
Figure FDA00039840886200000415
得到估计误差系统:
Step 2: Define the estimated error
Figure FDA00039840886200000414
and control output estimation error
Figure FDA00039840886200000415
The estimated error system is obtained:
Figure FDA00039840886200000416
Figure FDA00039840886200000416
Figure FDA00039840886200000417
Figure FDA00039840886200000417
式中,
Figure FDA00039840886200000418
为第k时刻的激励函数,
Figure FDA00039840886200000419
为在第k时刻的非线性激励函数,
Figure FDA00039840886200000420
是神经网络在第k时刻的状态估计,
Figure FDA00039840886200000421
是神经网络在第k-d时刻的状态估计,
Figure FDA00039840886200000422
是神经网络在第k-ι+1时刻的状态估计,
Figure FDA00039840886200000423
为神经网络动态模型状态的实数域,n为维数,χk表示在第k时刻的放大系数,
Figure FDA00039840886200000424
表示的是开根号的值,Kk是在第k时刻估计器增益矩阵,
Figure FDA0003984088620000051
表示在第k时刻传感器具备的随机能量的期望,
Figure FDA0003984088620000052
表示在第k时刻传感器具备的随机能量的期望,ΔAk为满足范数有界不确定性的第一号矩阵,ΔAdk为满足范数有界不确定性的第二号矩阵,ΔBk为满足范数有界不确定性的第三号矩阵,
Figure FDA0003984088620000053
为定义的左右区间的第一号矩阵,
Figure FDA0003984088620000054
为定义的左右区间的第二号矩阵,
Figure FDA0003984088620000055
为定义的左右区间的第三号矩阵,ek是在第k时刻的估计误差,ek+1是在第k+1时刻的估计误差,ek-d是在第k-d时刻的估计误差,
Figure FDA0003984088620000056
是在第k时刻的被控输出估计误差,A(xk)=diagn{ai(xik)}为在第k时刻的神经网络自反馈对角矩阵,diag{·}表示的是对角矩阵,ai(xik)为A(xk)的第i个分量,n为维数;Ad(xk)为在第k时刻的已知维数且与时滞相关的系统矩阵,B(xk)为在第k时刻已知的连接激励函数的权重矩阵;f(xk)为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;v1k为在第k时刻均值为零并且协方差为V1>0的高斯白噪声序列,v2k为在第k时刻均值为零并且协方差为V2>0的高斯白噪声序列,
Figure FDA0003984088620000057
表示的是ι=1到k+1求和的值,θs1,k是第k时刻在传感器-中继器信道的白噪声序列且满足
Figure FDA0003984088620000058
Figure FDA0003984088620000059
表示的是数学期望,(θs1,k)T是在第k时刻θs1,k的转置,θs2,k是第k时刻在中继器-估计器信道的白噪声信号且满足
Figure FDA00039840886200000510
Figure FDA00039840886200000511
表示的是数学期望,(θs2,k)T是在第k时刻θs2,k的转置,
Figure FDA00039840886200000512
表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,
Figure FDA00039840886200000513
是第k时刻已知信道衰减矩阵,m表示的是第m个信道。
In the formula,
Figure FDA00039840886200000418
is the activation function at the kth moment,
Figure FDA00039840886200000419
is the nonlinear activation function at the kth moment,
Figure FDA00039840886200000420
is the state estimate of the neural network at the kth moment,
Figure FDA00039840886200000421
is the state estimate of the neural network at the kdth moment,
Figure FDA00039840886200000422
is the state estimate of the neural network at the k-ι+1th moment,
Figure FDA00039840886200000423
is the real number domain of the dynamic model state of the neural network, n is the dimension, χ k represents the amplification factor at the kth moment,
Figure FDA00039840886200000424
represents the value of the square root, K k is the estimator gain matrix at the kth moment,
Figure FDA0003984088620000051
represents the expected random energy of the sensor at the kth moment,
Figure FDA0003984088620000052
represents the expectation of the random energy of the sensor at the kth moment, ΔA k is the first matrix satisfying the norm bounded uncertainty, ΔA dk is the second matrix satisfying the norm bounded uncertainty, ΔB k is the third matrix satisfying the norm bounded uncertainty,
Figure FDA0003984088620000053
is the first matrix of the defined left and right intervals,
Figure FDA0003984088620000054
is the second matrix of the defined left and right intervals,
Figure FDA0003984088620000055
is the third matrix of the defined left and right intervals, e k is the estimated error at the kth moment, e k+1 is the estimated error at the k+1th moment, e kd is the estimated error at the kdth moment,
Figure FDA0003984088620000056
is the estimated error of the controlled output at the kth moment, A( xk ) = diagn {a i ( xik )} is the diagonal matrix of the neural network self-feedback at the kth moment, diag{·} represents a diagonal matrix, a i ( xik ) is the i-th component of A( xk ), and n is the dimension; Ad ( xk ) is the system matrix of known dimension and time delay related at the kth moment, B( xk ) is the weight matrix of the connection activation function known at the kth moment; f( xk ) is the nonlinear activation function at the kth moment; C1k is the noise distribution matrix of the system with the first component known at the kth moment, C2k is the noise distribution matrix of the system with the second component known at the kth moment, Hk is the adjustment matrix of the known measurement at the kth moment; Dk is the measurement matrix of the known measurement at the kth moment; v1k is a Gaussian white noise sequence with zero mean and covariance V1 >0 at the kth moment, and v2k is a Gaussian white noise sequence with zero mean and covariance V2 at the kth moment. >0 Gaussian white noise sequence,
Figure FDA0003984088620000057
represents the sum of ι=1 to k+1, θ s1,k is the white noise sequence in the sensor-repeater channel at the kth moment and satisfies
Figure FDA0003984088620000058
Figure FDA0003984088620000059
represents the mathematical expectation, (θ s1,k ) T is the transpose of θ s1,k at the kth time, θ s2,k is the white noise signal in the repeater-estimator channel at the kth time and satisfies
Figure FDA00039840886200000510
Figure FDA00039840886200000511
represents the mathematical expectation, (θ s2,k ) T is the transpose of θ s2,k at the kth moment,
Figure FDA00039840886200000512
represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix,
Figure FDA00039840886200000513
is the known channel attenuation matrix at the kth moment, and m represents the mth channel.
5.根据权利要求4所述的方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述ps,k满足如下统计特性:5. The variance-constrained fractional-order memristor neural network estimation method according to claim 4, wherein p s,k satisfies the following statistical characteristics:
Figure FDA00039840886200000514
Figure FDA00039840886200000514
式中,Pr{·}表示的是数学概率,
Figure FDA0003984088620000061
表示所有概率的求和值为1,并且概率满足区间
Figure FDA0003984088620000062
Figure FDA0003984088620000063
为第k时刻的传感器具备的随机能量的期望值,φ表示的是所有信道的数量。
In the formula, Pr{·} represents the mathematical probability,
Figure FDA0003984088620000061
Indicates that the sum of all probabilities is 1, and the probability satisfies the interval
Figure FDA0003984088620000062
Figure FDA0003984088620000063
is the expected value of the random energy of the sensor at the kth moment, and φ represents the number of all channels.
6.根据权利要求4所述的方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述随机能量ns,k具有如下统计特性:6. The variance-constrained fractional-order memristor neural network estimation method according to claim 4, wherein the random energy ns,k has the following statistical characteristics:
Figure FDA0003984088620000064
Figure FDA0003984088620000064
式中,
Figure FDA0003984088620000065
表示所有概率的求和值为1,并且概率满足区间
Figure FDA0003984088620000066
Figure FDA0003984088620000067
为第k时刻的传输随机能量的期望值,ψ表示的是所有的信道的数量。
In the formula,
Figure FDA0003984088620000065
Indicates that the sum of all probabilities is 1, and the probability satisfies the interval
Figure FDA0003984088620000066
Figure FDA0003984088620000067
is the expected value of the random energy transmitted at the kth moment, and ψ represents the number of all channels.
7.根据权利要求4所述的方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述估计误差系统同时满足以下两个性能约束要求:7. The variance-constrained fractional-order memristor neural network estimation method according to claim 4, characterized in that the estimation error system satisfies the following two performance constraints at the same time: (1)令扰动衰减水平γ>0,半正定矩阵一号和半正定矩阵二号分别为
Figure FDA0003984088620000068
Figure FDA0003984088620000069
对于初始状态e0,控制输出估计误差
Figure FDA00039840886200000610
满足如下的H性能约束条件:
(1) Let the disturbance attenuation level γ>0, and the semi-positive definite matrix No. 1 and the semi-positive definite matrix No. 2 are respectively
Figure FDA0003984088620000068
and
Figure FDA0003984088620000069
For the initial state e 0 , the control output estimation error
Figure FDA00039840886200000610
The following H∞ performance constraints are met:
Figure FDA00039840886200000611
Figure FDA00039840886200000611
式中,N为有限的节点个数,
Figure FDA00039840886200000612
表示的是数学期望,
Figure FDA00039840886200000613
是第一号权重矩阵,
Figure FDA00039840886200000614
是第一号权重矩阵,e0是在第0时刻的估计误差,γ>0是给定的扰动衰减水平,
Figure FDA00039840886200000615
是噪声v1k和v2k增广的向量,
Figure FDA00039840886200000616
是在第k时刻ek的转置,||·||表示的是范数形式,||·||2表示的是范数平方的形式;
Where N is the finite number of nodes.
Figure FDA00039840886200000612
represents the mathematical expectation,
Figure FDA00039840886200000613
is the first weight matrix,
Figure FDA00039840886200000614
is the first weight matrix, e 0 is the estimation error at time 0, γ>0 is the given disturbance attenuation level,
Figure FDA00039840886200000615
is the vector augmented by the noise v 1k and v 2k ,
Figure FDA00039840886200000616
is the transpose of e k at the kth moment, ||·|| represents the norm form, and ||·|| 2 represents the squared norm form;
(2)估计误差协方差满足如下的上界约束条件:(2) The estimated error covariance satisfies the following upper bound constraints:
Figure FDA00039840886200000617
Figure FDA00039840886200000617
式中,
Figure FDA00039840886200000618
是在第k时刻ek的转置,Πk(0≤k<N)是在第k时刻的一系列预先给定的可接受的估计精度矩阵。
In the formula,
Figure FDA00039840886200000618
is the transpose of e k at the kth moment, and Π k (0≤k<N) is a series of pre-given acceptable estimation accuracy matrices at the kth moment.
8.根据权利要求1所述的方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述步骤三的具体步骤如下:8. The variance-constrained fractional-order memristor neural network estimation method according to claim 1, wherein the specific steps of step three are as follows: 步骤三一、按照下式证明出H性能分析问题并给出相应的易于求解的判别准则:Step 3: 1. Prove the H∞ performance analysis problem according to the following formula and give the corresponding easy-to-solve judgment criteria:
Figure FDA0003984088620000071
Figure FDA0003984088620000071
式中:Where:
Figure FDA0003984088620000072
Figure FDA0003984088620000072
Figure FDA0003984088620000073
Figure FDA0003984088620000073
Figure FDA0003984088620000074
Figure FDA0003984088620000074
Figure FDA0003984088620000075
Figure FDA0003984088620000075
Figure FDA0003984088620000076
Figure FDA0003984088620000076
Figure FDA0003984088620000077
Figure FDA0003984088620000077
Figure FDA0003984088620000078
Figure FDA0003984088620000078
Figure FDA0003984088620000079
Figure FDA0003984088620000079
Figure FDA00039840886200000710
Figure FDA00039840886200000710
Figure FDA00039840886200000711
Figure FDA00039840886200000711
Figure FDA00039840886200000712
Figure FDA00039840886200000712
Figure FDA00039840886200000713
Figure FDA00039840886200000713
Figure FDA00039840886200000714
Figure FDA00039840886200000714
Figure FDA00039840886200000715
Figure FDA00039840886200000715
式中,γ为给定的正标量;
Figure FDA00039840886200000716
为半正定矩阵一号,
Figure FDA00039840886200000717
Figure FDA0003984088620000081
分别为
Figure FDA0003984088620000082
Dk、Kk、Et,k、Ct,k、ΔAk、Hk
Figure FDA0003984088620000083
ΔBk
Figure FDA0003984088620000084
ΔAk
Figure FDA0003984088620000085
Ek、Kk、Ck、R3k的转置;
Figure FDA0003984088620000086
为半正定矩阵;Y11是Y的第1行第1列分块矩阵,Y12是Y的第1行第2列分块矩阵,Y22是Y的第2行第2列分块矩阵,Y33是Y的第3行第3列分块矩阵,Y44是Y的第4行第4列分块矩阵,Y55是Y的第5行第5列分块矩阵,Y66是Y的第6行第6列分块矩阵,Y77是Y的第7行第7列分块矩阵,Y88是Y的第8行第8列分块矩阵,Y99是Y的第9行第9列分块矩阵,
Figure FDA0003984088620000087
表示在第k时刻传感器具备的随机能量的期望,
Figure FDA0003984088620000088
表示在第k时刻传感器具备的随机能量的期望,ΔAk为满足范数有界不确定性的第一号矩阵,ΔAdk为满足范数有界不确定性的第二号矩阵,ΔBk为满足范数有界不确定性的第三号矩阵,
Figure FDA0003984088620000089
为定义的左右区间的第一号矩阵,
Figure FDA00039840886200000810
为定义的左右区间的第二号矩阵,
Figure FDA00039840886200000811
为定义的左右区间的第三号矩阵,
Figure FDA00039840886200000812
为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;
Figure FDA00039840886200000813
表示的是ι=1到k+1求和的值,
Figure FDA00039840886200000814
表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,
Figure FDA00039840886200000815
是第k时刻已知信道衰减矩阵,m表示的是第m个信道,
Figure FDA00039840886200000816
Figure FDA00039840886200000817
分别是第一个、第二个、第三个、第四个和第五个相关比例系数,0表示的是矩阵块的元素均为0;
In the formula, γ is a given positive scalar;
Figure FDA00039840886200000716
is a semi-positive definite matrix number one,
Figure FDA00039840886200000717
Figure FDA0003984088620000081
They are
Figure FDA0003984088620000082
D k , K k , E t,k , C t,k , ΔA k , H k ,
Figure FDA0003984088620000083
ΔB k ,
Figure FDA0003984088620000084
ΔA k ,
Figure FDA0003984088620000085
Transpose of E k , K k , C k , R 3k ;
Figure FDA0003984088620000086
is a semi-positive definite matrix; Y 11 is the 1st row and 1st column block matrix of Y, Y 12 is the 1st row and 2nd column block matrix of Y, Y 22 is the 2nd row and 2nd column block matrix of Y, Y 33 is the 3rd row and 3rd column block matrix of Y, Y 44 is the 4th row and 4th column block matrix of Y, Y 55 is the 5th row and 5th column block matrix of Y, Y 66 is the 6th row and 6th column block matrix of Y, Y 77 is the 7th row and 7th column block matrix of Y, Y 88 is the 8th row and 8th column block matrix of Y, Y 99 is the 9th row and 9th column block matrix of Y,
Figure FDA0003984088620000087
represents the expected random energy of the sensor at the kth moment,
Figure FDA0003984088620000088
represents the expectation of the random energy of the sensor at the kth moment, ΔA k is the first matrix satisfying the norm bounded uncertainty, ΔA dk is the second matrix satisfying the norm bounded uncertainty, ΔB k is the third matrix satisfying the norm bounded uncertainty,
Figure FDA0003984088620000089
is the first matrix of the defined left and right intervals,
Figure FDA00039840886200000810
is the second matrix of the defined left and right intervals,
Figure FDA00039840886200000811
is the third matrix of the defined left and right intervals,
Figure FDA00039840886200000812
is the nonlinear excitation function at the kth moment; C 1k is the noise distribution matrix of the first component known system at the kth moment, C 2k is the noise distribution matrix of the second component known system at the kth moment, H k is the adjustment matrix of the known measurement at the kth moment; D k is the measurement matrix of the known measurement at the kth moment;
Figure FDA00039840886200000813
It represents the sum of ι=1 to k+1.
Figure FDA00039840886200000814
represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix,
Figure FDA00039840886200000815
is the known channel attenuation matrix at the kth moment, m represents the mth channel,
Figure FDA00039840886200000816
and
Figure FDA00039840886200000817
They are the first, second, third, fourth and fifth related proportional coefficients, and 0 means that all elements of the matrix block are 0;
步骤三二、探讨协方差矩阵
Figure FDA00039840886200000818
的上界约束问题,并给出如下充分条件:
Step 3.2: Explore the covariance matrix
Figure FDA00039840886200000818
The upper bound constraint problem is given, and the following sufficient conditions are given:
Sk+1≥Ω(Sk), (4)S k+1 ≥Ω(S k ), (4) 式中,In the formula,
Figure FDA0003984088620000091
Figure FDA0003984088620000091
Figure FDA0003984088620000092
Figure FDA0003984088620000092
式中,ek为在第k时刻的误差矩阵;
Figure FDA0003984088620000093
为在第k时刻的状态估计,ρ∈(0,1)为已知的调节正常数;Sk为在第k时刻的误差协方差矩阵的上界;
Figure FDA0003984088620000094
Θ1k T
Figure FDA0003984088620000095
Figure FDA0003984088620000096
分别为
Figure FDA0003984088620000097
Θ1k
Figure FDA0003984088620000098
C1k、Φι、Ct,k、Et,k的转置;ζ为调节系数,Ω(Sk)为在第k时刻求解出的上界矩阵;Sk-d为在第k-d时刻的误差协方差矩阵的上界矩阵;tr(Sk)为在第k时刻的误差协方差矩阵上界的迹;tr()为矩阵的迹,
Figure FDA0003984088620000099
为在第k时刻的误差上界,ek为在第k时刻的误差矩阵,I为单位矩阵,
Figure FDA00039840886200000910
是第1个分量在k时刻的已知适当维数的第一号实矩阵,
Figure FDA00039840886200000911
是第2个分量在k时刻的已知适当维数的第二号实矩阵。
Where, e k is the error matrix at the kth moment;
Figure FDA0003984088620000093
is the state estimate at the kth moment, ρ∈(0,1) is a known normal constant; S k is the upper bound of the error covariance matrix at the kth moment;
Figure FDA0003984088620000094
Θ 1k T ,
Figure FDA0003984088620000095
Figure FDA0003984088620000096
They are
Figure FDA0003984088620000097
Θ 1k
Figure FDA0003984088620000098
C 1k , Φ ι , C t,k , E t,k are transposed; ζ is the adjustment coefficient, Ω(S k ) is the upper bound matrix solved at the kth time; S kd is the upper bound matrix of the error covariance matrix at the kdth time; tr(S k ) is the trace of the upper bound of the error covariance matrix at the kth time; tr() is the trace of the matrix,
Figure FDA0003984088620000099
is the upper bound of the error at the kth moment, e k is the error matrix at the kth moment, I is the unit matrix,
Figure FDA00039840886200000910
is the first real matrix of known appropriate dimension of the first component at time k,
Figure FDA00039840886200000911
is the second real matrix of known appropriate dimension of the second component at time k.
9.根据权利要求1所述的方差受限下的分数阶忆阻神经网络估计方法,其特征在于所述步骤四中,通过求解(5)~(7)一系列递推线性矩阵不等式,给出估计误差系统同时满足H性能要求和误差协方差有上界的充分条件,即可计算出估计器增益矩阵的值:9. The variance-constrained fractional-order memristor neural network estimation method according to claim 1, characterized in that in the step 4, by solving a series of recursive linear matrix inequalities (5) to (7), sufficient conditions are given for the estimation error system to simultaneously meet the H∞ performance requirements and the upper bound of the error covariance, and the value of the estimator gain matrix can be calculated:
Figure FDA00039840886200000912
Figure FDA00039840886200000912
Figure FDA0003984088620000101
Figure FDA0003984088620000101
Sk+1k+1≤0 (7)S k+1k+1 ≤0 (7) 更新矩阵为:The updated matrix is:
Figure FDA0003984088620000102
Figure FDA0003984088620000102
式中:Where:
Figure FDA0003984088620000103
Figure FDA0003984088620000103
Ω22=diag{-ε1,kI,-ε2,kI,-ε2,kI,-ε3,kI,-ε3,kI},Ω 22 =diag{-ε 1,k I,-ε 2,k I,-ε 2,k I,-ε 3,k I,-ε 3,k I}, Ω33=diag{-ε4,kI,-ε4,kI,-ε5,kI,-ε5,kI},Ω 33 =diag{-ε 4,k I,-ε 4,k I,-ε 5,k I,-ε 5,k I},
Figure FDA0003984088620000104
Figure FDA0003984088620000104
Figure FDA0003984088620000105
Figure FDA0003984088620000105
Figure FDA0003984088620000106
Figure FDA0003984088620000106
Figure FDA0003984088620000107
Figure FDA0003984088620000107
Figure FDA0003984088620000108
Figure FDA0003984088620000108
Figure FDA0003984088620000109
Figure FDA0003984088620000109
Figure FDA00039840886200001010
Figure FDA00039840886200001010
Figure FDA0003984088620000111
Figure FDA0003984088620000111
Figure FDA0003984088620000112
Figure FDA0003984088620000112
Figure FDA0003984088620000113
Figure FDA0003984088620000113
Figure FDA0003984088620000114
Figure FDA0003984088620000114
Figure FDA0003984088620000115
Figure FDA0003984088620000115
Figure FDA0003984088620000116
Figure FDA0003984088620000116
Figure FDA0003984088620000117
Figure FDA0003984088620000117
Figure FDA0003984088620000118
Figure FDA0003984088620000118
Figure FDA0003984088620000119
Figure FDA0003984088620000119
Figure FDA00039840886200001110
Figure FDA00039840886200001110
Figure FDA00039840886200001111
Figure FDA00039840886200001111
Figure FDA00039840886200001112
Figure FDA00039840886200001112
Figure FDA00039840886200001113
Figure FDA00039840886200001113
Figure FDA00039840886200001114
Figure FDA00039840886200001114
Figure FDA00039840886200001115
Figure FDA00039840886200001115
Figure FDA00039840886200001116
Figure FDA00039840886200001116
式中,Ω11是第1行第1列分块矩阵,Ω12是第1行第2列分块矩阵,Ω13是第1行第3列分块矩阵,Ω22是第2行第2列分块矩阵,Ω33是第3行第3列分块矩阵,
Figure FDA00039840886200001117
Wherein, Ω 11 is the 1st row and 1st column block matrix, Ω 12 is the 1st row and 2nd column block matrix, Ω 13 is the 1st row and 3rd column block matrix, Ω 22 is the 2nd row and 2nd column block matrix, Ω 33 is the 3rd row and 3rd column block matrix,
Figure FDA00039840886200001117
是第1行第1列分块矩阵,
Figure FDA0003984088620000121
是第1行第2列分块矩阵,
Figure FDA0003984088620000122
是第1行第3列分块矩阵,
Figure FDA0003984088620000123
是第1行第4列分块矩阵,L15是第1行第5列分块矩阵,L16是第1行第6列分块矩阵,L22是第2行第2列分块矩阵,L33是第3行第3列分块矩阵,L44是第4行第4列分块矩阵,L55是第5行第5列分块矩阵,L66是第6行第6列分块矩阵,
Figure FDA0003984088620000124
是第1行第1列分块矩阵,G12是第1行第2列分块矩阵,G14是第1行第4列分块矩阵,
Figure FDA0003984088620000125
是第1行第5列分块矩阵,G22是第2行第2列分块矩阵,G24是第2行第4列分块矩阵,
Figure FDA0003984088620000126
是第2行第6列分块矩阵,
Figure FDA0003984088620000127
是第2行第7列分块矩阵,
Figure FDA0003984088620000128
是第2行第8列分块矩阵,G33是第3行第3列分块矩阵,G39是第3行第9列分块矩阵,
Figure FDA0003984088620000129
是第4行第10列分块矩阵,
Figure FDA00039840886200001210
是第4行第4列分块矩阵,
Figure FDA00039840886200001211
是第5行第5列分块矩阵,
Figure FDA00039840886200001212
是第6行第6列分块矩阵,
Figure FDA00039840886200001213
是第7行第7列分块矩阵,
Figure FDA00039840886200001214
是第8行第8列分块矩阵,
Figure FDA00039840886200001215
是第9行第9列分块矩阵,
Figure FDA00039840886200001216
是第10行第10列分块矩阵,
Figure FDA00039840886200001217
is a 1-row, 1-column block matrix,
Figure FDA0003984088620000121
is a block matrix with row 1 and column 2,
Figure FDA0003984088620000122
is a block matrix with row 1 and column 3,
Figure FDA0003984088620000123
is the 1st row and 4th column block matrix, L 15 is the 1st row and 5th column block matrix, L 16 is the 1st row and 6th column block matrix, L 22 is the 2nd row and 2nd column block matrix, L 33 is the 3rd row and 3rd column block matrix, L 44 is the 4th row and 4th column block matrix, L 55 is the 5th row and 5th column block matrix, L 66 is the 6th row and 6th column block matrix,
Figure FDA0003984088620000124
is the 1st row and 1st column block matrix, G12 is the 1st row and 2nd column block matrix, G14 is the 1st row and 4th column block matrix,
Figure FDA0003984088620000125
is the 1st row and 5th column block matrix, G 22 is the 2nd row and 2nd column block matrix, G 24 is the 2nd row and 4th column block matrix,
Figure FDA0003984088620000126
is the 2nd row and 6th column block matrix,
Figure FDA0003984088620000127
is the 2nd row and 7th column block matrix,
Figure FDA0003984088620000128
is the 2nd row and 8th column block matrix, G 33 is the 3rd row and 3rd column block matrix, G 39 is the 3rd row and 9th column block matrix,
Figure FDA0003984088620000129
is the 4th row and 10th column block matrix,
Figure FDA00039840886200001210
is the 4th row and 4th column block matrix,
Figure FDA00039840886200001211
is the 5th row and 5th column block matrix,
Figure FDA00039840886200001212
is the 6th row and 6th column block matrix,
Figure FDA00039840886200001213
is the 7th row and 7th column block matrix,
Figure FDA00039840886200001214
is the 8th row and 8th column block matrix,
Figure FDA00039840886200001215
is the 9th row and 9th column block matrix,
Figure FDA00039840886200001216
is the 10th row and 10th column block matrix,
Figure FDA00039840886200001217
分别为
Figure FDA00039840886200001218
Dk,Kk,Et,k,Ct,k,ΔAk,Hk
Figure FDA00039840886200001219
ΔBk
Figure FDA00039840886200001220
ΔAk
Figure FDA00039840886200001221
Ek,Kk,Ck,R3k的转置,
Figure FDA00039840886200001222
为定义的左右区间的第一号矩阵,
Figure FDA00039840886200001223
为定义的左右区间的第二号矩阵,
Figure FDA00039840886200001224
为定义的左右区间的第三号矩阵,
Figure FDA00039840886200001225
为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;
Figure FDA00039840886200001226
表示的是ι=1到k+1求和的值,
Figure FDA00039840886200001227
表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,
Figure FDA00039840886200001228
是第k时刻已知信道衰减矩阵,m表示的是第m个信道,ρ∈(0,1)为已知的调节正常数;Sk为在第k时刻的误差协方差矩阵的上界;
Figure FDA00039840886200001229
Θ1k T
Figure FDA00039840886200001230
Figure FDA00039840886200001231
分别为
Figure FDA00039840886200001232
Θ1k
Figure FDA00039840886200001233
C1k,Φι,Ct,k,Et,k的转置;ζ为调节系数,Ω(Sk)为在第k时刻求解出的上界矩阵;Sk-d为在第k-d时刻的误差协方差矩阵的上界矩阵;tr(Sk)为在第k时刻的误差协方差矩阵上界的迹;tr()为矩阵的迹,I为单位矩阵;
Figure FDA0003984088620000131
为在第k时刻的第一号权重矩阵;
Figure FDA0003984088620000132
为在第k时刻的第二号权重矩阵;
Figure FDA0003984088620000133
为在第k时刻的第三号权重矩阵;
Figure FDA0003984088620000134
是在第k时刻R3k的转置;
Figure FDA0003984088620000135
是第1个分量在k时刻的已知适当维数的第一号实矩阵,
Figure FDA0003984088620000136
是第2个分量在k时刻的已知适当维数的第二号实矩阵;
Figure FDA0003984088620000137
为在第k时刻的非线性激励函数的状态估计;
Figure FDA0003984088620000138
是第1个分量在k时刻的已知适当维数的第一号量度矩阵;
Figure FDA0003984088620000139
是第2个分量在k时刻的已知适当维数的第二号量度矩阵;
Figure FDA00039840886200001310
是第3个分量在k时刻的已知适当维数的第三号量度矩阵;
Figure FDA00039840886200001311
是第4个分量在k时刻的已知适当维数的第三号量度矩阵;
Figure FDA00039840886200001312
是第5个分量在k时刻的已知适当维数的第三号量度矩阵;M1,M2,M3,M4和M5分别是第一号,第二号,第三号,第四号和第五号的量度矩阵,
Figure FDA00039840886200001313
为在第k时刻的神经元状态估计,
Figure FDA00039840886200001314
为第k时刻的半正定矩阵;
Figure FDA00039840886200001315
为第k时刻的半正定矩阵;
Figure FDA00039840886200001316
为第k-d时刻的半正定矩阵;
Figure FDA00039840886200001317
为在第k+1时刻的第一更新矩阵,Sk为估计误差的上界矩阵,tr(Sk)为在第k时刻估计误差上界矩阵Sk的迹;Sk-d为在k-d时刻的上界矩阵,κ为调节的权重系数,
Figure FDA00039840886200001318
Figure FDA00039840886200001319
均为已知的实值权重矩阵,
Figure FDA00039840886200001320
是未知矩阵且满足
Figure FDA00039840886200001321
Figure FDA00039840886200001322
Figure FDA00039840886200001323
的转置,γ为给定的正标量;
Figure FDA00039840886200001324
为给定的半正定矩阵一号;
Figure FDA00039840886200001325
分别是Ω12,Ω13
Figure FDA00039840886200001326
的转置;
Figure FDA00039840886200001327
Figure FDA00039840886200001328
分别是,G12,G14
Figure FDA00039840886200001329
G24,G410
Figure FDA00039840886200001330
的转置;
Figure FDA00039840886200001331
分别是M1,M2,M3,M4,M5的转置;N1,N2,N3,N4,N5分别是
Figure FDA00039840886200001332
的转置;
Figure FDA00039840886200001333
Figure FDA00039840886200001334
分别是第一个、第二个、第三个、第四个和第五个相关比例系数,0表示的是矩阵块的元素均为0。
They are
Figure FDA00039840886200001218
D k , K k , E t,k , C t,k , ΔA k , H k ,
Figure FDA00039840886200001219
ΔB k
Figure FDA00039840886200001220
ΔA k
Figure FDA00039840886200001221
The transpose of E k , K k , C k , R 3k ,
Figure FDA00039840886200001222
is the first matrix of the defined left and right intervals,
Figure FDA00039840886200001223
is the second matrix of the defined left and right intervals,
Figure FDA00039840886200001224
is the third matrix of the defined left and right intervals,
Figure FDA00039840886200001225
is the nonlinear excitation function at the kth moment; C 1k is the noise distribution matrix of the first component known system at the kth moment, C 2k is the noise distribution matrix of the second component known system at the kth moment, H k is the adjustment matrix of the known measurement at the kth moment; D k is the measurement matrix of the known measurement at the kth moment;
Figure FDA00039840886200001226
It represents the sum of ι=1 to k+1.
Figure FDA00039840886200001227
represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix,
Figure FDA00039840886200001228
is the known channel attenuation matrix at the kth moment, m represents the mth channel, ρ∈(0,1) is the known normal constant; S k is the upper bound of the error covariance matrix at the kth moment;
Figure FDA00039840886200001229
Θ 1k T ,
Figure FDA00039840886200001230
Figure FDA00039840886200001231
They are
Figure FDA00039840886200001232
Θ 1k ,
Figure FDA00039840886200001233
C 1k , Φ ι , C t,k , the transpose of E t,k ; ζ is the adjustment coefficient, Ω(S k ) is the upper bound matrix solved at the kth time; S kd is the upper bound matrix of the error covariance matrix at the kdth time; tr(S k ) is the trace of the upper bound of the error covariance matrix at the kth time; tr() is the trace of the matrix, and I is the identity matrix;
Figure FDA0003984088620000131
is the first weight matrix at the kth moment;
Figure FDA0003984088620000132
is the second weight matrix at the kth moment;
Figure FDA0003984088620000133
is the third weight matrix at the kth moment;
Figure FDA0003984088620000134
is the transpose of R 3k at the kth moment;
Figure FDA0003984088620000135
is the first real matrix of known appropriate dimension of the first component at time k,
Figure FDA0003984088620000136
is the second real matrix of known appropriate dimension of the second component at time k;
Figure FDA0003984088620000137
is the state estimate of the nonlinear activation function at the kth moment;
Figure FDA0003984088620000138
is the first metric matrix of known appropriate dimension of the first component at time k;
Figure FDA0003984088620000139
is the second metric matrix of known appropriate dimension of the second component at time k;
Figure FDA00039840886200001310
is the third metric matrix of known appropriate dimension of the third component at time k;
Figure FDA00039840886200001311
is the third metric matrix of known appropriate dimension of the fourth component at time k;
Figure FDA00039840886200001312
is the third metric matrix of known appropriate dimension of the fifth component at time k; M 1 , M 2 , M 3 , M 4 and M 5 are the first, second, third, fourth and fifth metric matrices respectively,
Figure FDA00039840886200001313
is the estimated neuron state at the kth moment,
Figure FDA00039840886200001314
is the semi-positive definite matrix at the kth moment;
Figure FDA00039840886200001315
is the semi-positive definite matrix at the kth moment;
Figure FDA00039840886200001316
is the semi-positive definite matrix at the kdth moment;
Figure FDA00039840886200001317
is the first update matrix at the k+1th time, Sk is the upper bound matrix of the estimation error, tr( Sk ) is the trace of the upper bound matrix Sk of the estimation error at the kth time; Skd is the upper bound matrix at the kdth time, κ is the adjusted weight coefficient,
Figure FDA00039840886200001318
and
Figure FDA00039840886200001319
are all known real-valued weight matrices,
Figure FDA00039840886200001320
is an unknown matrix and satisfies
Figure FDA00039840886200001321
Figure FDA00039840886200001322
yes
Figure FDA00039840886200001323
The transpose of , γ is a given positive scalar;
Figure FDA00039840886200001324
For the given semi-positive definite matrix one;
Figure FDA00039840886200001325
They are Ω 12 , Ω 13 ,
Figure FDA00039840886200001326
The transpose of
Figure FDA00039840886200001327
Figure FDA00039840886200001328
They are G12 , G14 ,
Figure FDA00039840886200001329
G24 , G410 ,
Figure FDA00039840886200001330
The transpose of
Figure FDA00039840886200001331
are the transposes of M 1 , M 2 , M 3 , M 4 , and M 5 respectively; N 1 , N 2 , N 3 , N 4 , and N 5 are
Figure FDA00039840886200001332
The transpose of
Figure FDA00039840886200001333
and
Figure FDA00039840886200001334
They are the first, second, third, fourth and fifth relevant proportional coefficients respectively, and 0 means that all elements of the matrix block are 0.
CN202211559637.5A 2022-12-06 2022-12-06 Fractional order memristor neural network estimation method under variance limitation Active CN116227324B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211559637.5A CN116227324B (en) 2022-12-06 2022-12-06 Fractional order memristor neural network estimation method under variance limitation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211559637.5A CN116227324B (en) 2022-12-06 2022-12-06 Fractional order memristor neural network estimation method under variance limitation

Publications (2)

Publication Number Publication Date
CN116227324A true CN116227324A (en) 2023-06-06
CN116227324B CN116227324B (en) 2023-09-19

Family

ID=86584853

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211559637.5A Active CN116227324B (en) 2022-12-06 2022-12-06 Fractional order memristor neural network estimation method under variance limitation

Country Status (1)

Country Link
CN (1) CN116227324B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117077748A (en) * 2023-06-15 2023-11-17 盐城工学院 Coupling synchronous control method and system for discrete memristor neural network
CN117949897A (en) * 2024-01-09 2024-04-30 哈尔滨理工大学 Multifunctional radar working mode identification method based on time sequence segmentation and clustering

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107436411A (en) * 2017-07-28 2017-12-05 南京航空航天大学 Battery SOH On-line Estimation methods based on fractional order neural network and dual-volume storage Kalman
CN109088749A (en) * 2018-07-23 2018-12-25 哈尔滨理工大学 The method for estimating state of complex network under a kind of random communication agreement
CN111025914A (en) * 2019-12-26 2020-04-17 东北石油大学 Method and device for remote state estimation of neural network system based on limited communication
US11449754B1 (en) * 2021-09-12 2022-09-20 Zhejiang University Neural network training method for memristor memory for memristor errors

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107436411A (en) * 2017-07-28 2017-12-05 南京航空航天大学 Battery SOH On-line Estimation methods based on fractional order neural network and dual-volume storage Kalman
CN109088749A (en) * 2018-07-23 2018-12-25 哈尔滨理工大学 The method for estimating state of complex network under a kind of random communication agreement
CN111025914A (en) * 2019-12-26 2020-04-17 东北石油大学 Method and device for remote state estimation of neural network system based on limited communication
US11449754B1 (en) * 2021-09-12 2022-09-20 Zhejiang University Neural network training method for memristor memory for memristor errors

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117077748A (en) * 2023-06-15 2023-11-17 盐城工学院 Coupling synchronous control method and system for discrete memristor neural network
CN117077748B (en) * 2023-06-15 2024-03-22 盐城工学院 Coupling synchronous control method and system for discrete memristor neural network
CN117949897A (en) * 2024-01-09 2024-04-30 哈尔滨理工大学 Multifunctional radar working mode identification method based on time sequence segmentation and clustering

Also Published As

Publication number Publication date
CN116227324B (en) 2023-09-19

Similar Documents

Publication Publication Date Title
CN116227324B (en) Fractional order memristor neural network estimation method under variance limitation
Bai Common breaks in means and variances for panel data
Bolzern et al. On almost sure stability of continuous-time Markov jump linear systems
Du et al. A fast multi-output RBF neural network construction method
Mead et al. A Newton root-finding algorithm for estimating the regularization parameter for solving ill-conditioned least squares problems
Antunes et al. Stochastic networked control systems with dynamic protocols
CN112115419A (en) System state estimation method and system state estimation device
Hammouri et al. Constant gain observer for continuous-discrete time uniformly observable systems
CN106972949A (en) A kind of fractional order network system situation method of estimation based on adaptive equalization technology
Guo et al. Identification for Wiener‐Hammerstein systems under quantized inputs and quantized output observations
CN109217844B (en) Hyperparameter optimization method based on pre-trained stochastic Fourier feature kernel LMS
CN105808962A (en) Assessment method considering voltage probabilities of multiple electric power systems with wind power output randomness
Horváth et al. Sample autocovariances of long-memory time series
CN115935787B (en) A Memristive Neural Network State Estimation Method Based on Encoding and Decoding Mechanism
CN109474258B (en) Nuclear parameter optimization method of random Fourier feature kernel LMS (least mean square) based on nuclear polarization strategy
CN115879634A (en) Weather environment factor and agricultural drought cause and effect inference method
Sverchkov A new approach to estimation of response probabilities when missing data are not missing at random
Xu et al. Differential game with mixed strategies: A weak approximation approach
Kawada et al. An application of a generalized inverse regression estimator to Taguchi's T-Method
Allende et al. Robust estimator for the learning process in neural networks applied in time series
Zhu et al. State estimation for a class of discrete-time BAM neural networks with multiple time-varying delays
CN117875032B (en) Optimal filtering method for state-saturated mass-spring system under dynamic-static hybrid protocol
Víšek Instrumental weighted variables
Liu et al. A unified recursive identification algorithm with quantized observations based on weighted least-squares type criteria
Xu et al. Quantised innovation Kalman filter: performance analysis and design of quantised level

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
GR01 Patent grant
GR01 Patent grant