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 PDFInfo
- 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
Links
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 92
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000011159 matrix material Substances 0.000 claims abstract description 389
- 230000003321 amplification Effects 0.000 claims abstract description 25
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 25
- 238000004458 analytical method Methods 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 claims description 36
- 230000004913 activation Effects 0.000 claims description 15
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000005284 excitation Effects 0.000 claims description 6
- 210000002569 neuron Anatomy 0.000 claims description 5
- 230000003190 augmentative effect Effects 0.000 claims description 2
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 4
- 238000003909 pattern recognition Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 2
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/06—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
- G06N3/063—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using electronic means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic 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∞性能约束及方差受限分数阶忆阻神经网络的状态估计导致的估计性能准确率低的问题,从而提高了估计性能的准确率。
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.
Description
技术领域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∞性能指标γ、半正定矩阵一号半正定矩阵二号及初始条件计算分数阶忆阻神经网络的误差协方差矩阵的上界及H∞性能约束条件;Step 3: Given the H∞ performance index γ and the semi-positive definite matrix No. Semi-positive definite matrix II and initial conditions 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
本发明中,所述神经网络可以为车辆悬挂构成的网络、质点弹簧构成的网络、航天器构成的网络或雷达构成的网络,在生物学、数学、计算机以及联想记忆、模式识别、组合优化和图像处理等多学科领域具有重要应用。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在两个不同情形下状态估计轨迹的对比图,zk为神经网络在第k时刻的状态变量;其中是系统状态轨迹,是情形一下的状态估计轨迹,是情形二下的状态估计轨迹;Figure 2 is the actual state trajectory z k of the fractional-order memristor neural network under two different situations. Comparison chart, z k is the state variable of the neural network at the kth moment; is the system state trajectory, is the state estimation trajectory for the case, is the state estimation trajectory under
图3是神经网络控制输出估计误差轨迹图在两个不同情形下的误差对比图;其中是情形一下的控制输出估计误差轨迹,是情形二下的控制输出估计误差轨迹;FIG3 is a comparison diagram of the error trajectory of the neural network control output estimation error in two different situations; is the control output estimation error trajectory for the case 1, is the control output estimation error trajectory under
图4是神经网络实际状态误差协方差和误差协方差上界第一个分量的轨迹图;其中是方差约束的轨迹,是实际误差协方差的轨迹;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; is the variance-constrained trajectory, is the locus of the actual error covariance;
图5是神经网络实际状态误差协方差和误差协方差上界第二个分量的轨迹图;其中是方差约束的轨迹,是实际误差协方差的轨迹。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; is the variance-constrained trajectory, 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:
式中,Δα表示α阶的Grunwald-Letnikov分数阶导数定义,h为相应的采样间隔,假设采样间隔为1,k为采样时刻,表示h→0的所有极限值,i!表示的是i的所有阶层,表示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. represents all the extreme values of h→0, i! represents all the classes of i, 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:
式中:Where:
这里,表示微分算子,为分数阶(j=1,2,…,n),n为维数,是在第k时刻的分数阶忆阻神经网络的状态向量,是在第k-ι+1时刻的分数阶忆阻神经网络的状态向量,是在第k-d时刻的分数阶忆阻神经网络的状态向量,是在第k+1时刻的分数阶忆阻神经网络的状态向量,为神经网络动态模型状态的实数域且其维数为n;为在第k时刻的被控测量输出,为神经网络动态模型被控输出状态的实数域且其维数为r;是给定的初始序列,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的高斯白噪声序列,表示的是ι=1到k+1求和的值。here, represents the differential operator, is a fractional order (j=1,2,…,n), n is the dimension, is the state vector of the fractional-order memristor neural network at the kth moment, is the state vector of the fractional-order memristor neural network at the k-ι+1th moment, is the state vector of the fractional-order memristor neural network at the kdth time, is the state vector of the fractional-order memristor neural network at the k+1th time, is the real number field of the states of the neural network dynamic model and its dimension is n; is the controlled measured output at the kth moment, is the real number domain of the controlled output state of the neural network dynamic model and its dimension is r; 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, 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:
式中,ai(xi,k)、aij,d(xi,k)和bij(xi,k)分别为A(xk),Ad(xk)和B(xk)的第i个分量,Ωi>0为已知的切换阈值,为第i个已知的上存储变量矩阵,为第i个已知的下存储变量矩阵,为第ij,d个已知的左存储变量矩阵,为第ij,d个已知的右存储变量矩阵,为第ij个已知的内存储变量矩阵,为第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, is the i-th known upper storage variable matrix, is the i-th known storage variable matrix, is the ij,dth known left storage variable matrix, is the ij,dth known right storage variable matrix, is the ijth known internal storage variable matrix, is the ijth known external storage variable matrix.
定义:definition:
式中,为第i个最小存储的第一号度量矩阵,为第i个已知的上存储区间变量矩阵,为第i个已知的下存储区间变量矩阵,min{·}表示两个存储矩阵中取最小值,max{·}表示两个存储矩阵中取最大值,为第i个最大存储的第一号度量矩阵,为第ij,d个最小存储的第二号度量矩阵,为第i个最大存储的第二号度量矩阵,为第ij,d个已知的左存储变量矩阵,为第ij,d个已知的右存储变量矩阵,为第ij个最小存储的第三号度量矩阵,为第ij个最大存储的第三号度量矩阵,为第ij个已知的内存储变量矩阵,为第ij个已知的外存储变量矩阵,diag{·}为对角矩阵,A-为定义的第一号对角矩阵,A+为定义的第二号对角矩阵,为定义的第三号对角矩阵,为定义的第四号对角矩阵,B-为定义的第五号对角矩阵,B+为定义的第六号对角矩阵,n为维数。In the formula, is the first metric matrix with the smallest storage of i, is the i-th known upper storage interval variable matrix, 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. is the first metric matrix with the largest storage in the i-th place, is the second metric matrix with the ij,dth smallest storage, is the second largest metric matrix stored in the i-th place, is the ij,dth known left storage variable matrix, is the ij,dth known right storage variable matrix, is the third smallest metric matrix stored for the ijth time. is the third largest metric matrix stored in the ijth order, is the ijth known internal storage variable matrix, 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, is the third diagonal matrix defined, 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+]、和B(xk)∈[B-,B+]。令 和则有:It is easy to derive A(x k )∈[A - ,A + ], and B(x k )∈[B − ,B + ]. Let and Then we have:
式中,为定义的左右区间的第一号矩阵,为定义的左右区间的第二号矩阵,为定义的左右区间的第三号矩阵,和满足范数有界不确定性:In the formula, is the first matrix of the defined left and right intervals, is the second matrix of the defined left and right intervals, is the third matrix of the defined left and right intervals, and Satisfies norm-bounded uncertainty:
式中,ΔAk为满足范数有界不确定性的第一号矩阵,ΔAdk为满足范数有界不确定性的第二号矩阵,ΔBk为满足范数有界不确定性的第三号矩阵,和均为已知的实值权重矩阵,是未知矩阵且满足 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. and are all known real-valued weight matrices, is an unknown matrix and satisfies
步骤二、放大转发协议下对步骤一建立的分数阶忆阻神经网络动态模型进行状态估计。具体步骤如下: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分别表示传感器和放大-转发中继器具有的随机能量,放大-转发中继器的输出信号由表示,其满足如下方程: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 It means that it satisfies the following equation:
式中,表示在第k时刻已知信道衰减矩阵,为已知信道的衰减矩阵的m个信道的分量形式,diag{·}表示的是对角矩阵,yk是第k时刻的理想测量输出,是第k时刻的实际测量输出,是第k时刻在传感器-中继器信道的白噪声序列且满足 表示的是数学期望,是在第k时刻的转置,ps,k表示在第k时刻传感器具备的随机能量,满足如下统计特性:In the formula, represents the known channel attenuation matrix at the kth time, 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, is the actual measured output at the kth moment, is a white noise sequence in the sensor-repeater channel at the kth moment and satisfies represents the mathematical expectation, At the kth moment The transpose of , ps,k represents the random energy of the sensor at the kth moment, satisfying the following statistical characteristics:
式中,Pr{·}表示的是数学概率,表示所有概率的求和值为1,并且概率满足区间 为第k时刻的传感器具备的随机能量的期望值,φ表示的是所有信道的数量。In the formula, Pr{·} represents the mathematical probability, Indicates that the sum of all probabilities is 1, and the probability satisfies the interval 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:
式中,χk>0表示在第k时刻的放大系数,是第k时刻已知信道的衰减矩阵,为已知信道的衰减矩阵的m个信道的分量形式,m表示的是第m个信道,ns,k为第k时刻的传输随机能量的变量,是第k时刻的实际测量输出,是第k时刻在中继器-估计器信道的白噪声信号且满足 表示的是数学期望,是在第k时刻的转置。同样地,随机能量ns,k具有如下统计特性:In the formula, χ k >0 represents the amplification factor at the kth moment, is the attenuation matrix of the known channel at the kth moment, 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, is the actual measured output at the kth moment, is the white noise signal in the repeater-estimator channel at the kth moment and satisfies represents the mathematical expectation, At the kth moment Similarly, the random energy n s,k has the following statistical properties:
式中,表示所有概率的求和值为1,并且概率满足区间 为第k时刻的传输随机能量的期望值,ψ表示的是所有的信道的数量。In the formula, Indicates that the sum of all probabilities is 1, and the probability satisfies the interval 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:
式中,是第1个分量在k时刻的已知适当维数的第一号实矩阵,是第2个分量的在k时刻的已知适当维数的第二号实矩阵。In the formula, is the first real matrix of known appropriate dimension of the first component at time k, 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:
式中,是神经网络在第k时刻的状态估计,是神经网络在第k时刻的状态估计,是神经网络在第k-d时刻的状态估计,为神经网络动态模型状态的实数域且其维数为n;χk表示在第k时刻的放大系数,d为一个固定的网络时滞,为在第k时刻的被控输出的状态估计,为神经网络动态模型被控输出状态的实数域且其维数为r,为定义的左右区间的第一号矩阵,为定义的左右区间的第二号矩阵,为定义的左右区间的第三号矩阵,为在第k时刻的非线性激励函数,Hk为在第k时刻的已知测量的调节矩阵,Dk是在第k时刻的已知测量的量度矩阵,是第k时刻解码器的测量输出,Kk是在第k时刻估计器增益矩阵,为传感器具备的随机能量期望的求和,表示在第k时刻传感器具备的随机能量的期望,为传感器具备的随机能量期望的求和,表示在第k时刻传感器具备的随机能量的期望,为所有的二项式组成的对角矩阵,为分数阶(j=1,2,…,n),n为维数,diag{·}表示的是对角矩阵,χk表示在第k时刻的放大系数。In the formula, is the state estimate of the neural network at the kth moment, is the state estimate of the neural network at the kth moment, is the state estimate of the neural network at the kdth moment, 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, is the state estimate of the controlled output at the kth moment, is the real number domain of the controlled output state of the neural network dynamic model and its dimension is r, is the first matrix of the defined left and right intervals, is the second matrix of the defined left and right intervals, is the third matrix of the defined left and right intervals, 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, is the measured output of the decoder at time k, Kk is the estimator gain matrix at time k, is the expected sum of random energies possessed by the sensor, represents the expected random energy of the sensor at the kth moment, is the expected sum of random energies possessed by the sensor, represents the expected random energy of the sensor at the kth moment, is the diagonal matrix of all binomials, 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.
步骤二三、定义估计误差和控制输出估计误差进一步,可以得到估计误差系统:Step 2: Define the estimated error and control output estimation error Furthermore, we can get the estimated error system:
式中,为第k时刻的激励函数,为在第k时刻的非线性激励函数,是神经网络在第k时刻的状态估计,是神经网络在第k-d时刻的状态估计,是神经网络在第k-ι+1时刻的状态估计,为神经网络动态模型状态的实数域,n为维数,χk表示在第k时刻的放大系数,表示的是开根号的值,Kk是在第k时刻估计器增益矩阵,表示在第k时刻传感器具备的随机能量的期望,表示在第k时刻传感器具备的随机能量的期望,ΔAk为满足范数有界不确定性的第一号矩阵,ΔAdk为满足范数有界不确定性的第二号矩阵,ΔBk为满足范数有界不确定性的第三号矩阵,为定义的左右区间的第一号矩阵,为定义的左右区间的第二号矩阵,为定义的左右区间的第三号矩阵,ek是在第k时刻的估计误差,ek+1是在第k+1时刻的估计误差,ek-d是在第k-d时刻的估计误差,是在第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的高斯白噪声序列,表示的是ι=1到k+1求和的值,是第k时刻在传感器-中继器信道的白噪声序列且满足 表示的是数学期望,是在第k时刻的转置,是第k时刻在中继器-估计器信道的白噪声信号且满足 表示的是数学期望,是在第k时刻的转置,表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,是第k时刻已知信道衰减矩阵,m表示的是第m个信道。In the formula, is the activation function at the kth moment, is the nonlinear activation function at the kth moment, is the state estimate of the neural network at the kth moment, is the state estimate of the neural network at the kdth moment, is the state estimate of the neural network at the k-ι+1th moment, 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, represents the value of the square root, K k is the estimator gain matrix at the kth moment, represents the expected random energy of the sensor at the kth moment, 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, is the first matrix of the defined left and right intervals, is the second matrix of the defined left and right intervals, 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, 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, It represents the sum of ι=1 to k+1. is a white noise sequence in the sensor-repeater channel at the kth moment and satisfies represents the mathematical expectation, At the kth moment The transpose of is the white noise signal in the repeater-estimator channel at the kth moment and satisfies represents the mathematical expectation, At the kth moment The transpose of represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix, 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,半正定矩阵一号和半正定矩阵二号分别为和对于初始状态e0,控制输出估计误差满足如下的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 and For the initial state e 0 , the control output estimation error The following H∞ performance constraints are met:
式中,N为有限的节点个数,表示的是数学期望,是第一号权重矩阵,是第一号权重矩阵,e0是在第0时刻的估计误差,γ>0是给定的扰动衰减水平,是噪声v1k和v2k增广的向量,是在第k时刻ek的转置,·表示的是范数形式,·2表示的是范数平方的形式。Where N is the finite number of nodes. represents the mathematical expectation, is the first weight matrix, is the first weight matrix, e 0 is the estimation error at
(2)估计误差协方差满足如下的上界约束条件:(2) The estimated error covariance satisfies the following upper bound constraints:
式中,是在第k时刻ek的转置,Πk(0≤k<N)是在第k时刻的一系列预先给定的可接受的估计精度矩阵。In the formula, 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∞性能指标γ、半正定矩阵一号半正定矩阵二号及初始条件计算分数阶忆阻神经网络的误差协方差矩阵的上界及H∞性能约束条件。具体步骤如下:Step 3: Given the H∞ performance index γ and the semi-positive definite matrix No. Semi-positive definite matrix II and initial conditions 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:
式中:Where:
式中,γ为给定的正标量;为半正定矩阵一号, 分别为Dk、Kk、Et,k、Ct,k、ΔAk、Hk、ΔBk、ΔAk、Ek、Kk、Ck、R3k的转置;为半正定矩阵;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列分块矩阵,表示在第k时刻传感器具备的随机能量的期望,表示在第k时刻传感器具备的随机能量的期望,ΔAk为满足范数有界不确定性的第一号矩阵,ΔAdk为满足范数有界不确定性的第二号矩阵,ΔBk为满足范数有界不确定性的第三号矩阵,为定义的左右区间的第一号矩阵,为定义的左右区间的第二号矩阵,为定义的左右区间的第三号矩阵,为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;表示的是ι=1到k+1求和的值,表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,是第k时刻已知信道衰减矩阵,m表示的是第m个信道,和分别是第一个、第二个、第三个、第四个和第五个相关比例系数,0表示的是矩阵块的元素均为0。Where γ is a given positive scalar; is a semi-positive definite matrix number one, They are D k , K k , E t,k , C t,k , ΔA k , H k , ΔB k , ΔA k , Transpose of E k , K k , C k , R 3k ; 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, represents the expectation of the random energy possessed by the sensor at the kth moment, 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, is the first matrix of the defined left and right intervals, is the second matrix of the defined left and right intervals, is the third matrix of the defined left and right intervals, 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; It represents the sum of ι=1 to k+1. represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix, is the known channel attenuation matrix at the kth moment, m represents the mth channel, and 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,
式中,ek为在第k时刻的误差矩阵;为在第k时刻的状态估计,ρ∈(0,1)为已知的调节正常数;Sk为在第k时刻的误差协方差矩阵的上界;Θ1k T、 分别为Θ1k、C1k、Φι、Ct,k、Et,k的转置;ζ为调节系数,Ω(Sk)为在第k时刻求解出的上界矩阵;Sk-d为在第k-d时刻的误差协方差矩阵的上界矩阵;tr(Sk)为在第k时刻的误差协方差矩阵上界的迹;tr()为矩阵的迹,χk=ekek T为在第k时刻的误差上界,ek为在第k时刻的误差矩阵,I为单位矩阵,是第1个分量在k时刻的已知适当维数的第一号实矩阵,是第2个分量在k时刻的已知适当维数的第二号实矩阵。Where, e k is the error matrix at the kth moment; 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; Θ 1k T , They are Θ 1k 、 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, is the first real matrix of known appropriate dimension of the first component at time k, 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
本步骤中,通过求解(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:
Sk+1-Ωk+1≤0 (7)S k+1 -Ω k+1 ≤0 (7)
更新矩阵为:The updated matrix is:
式中: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},
式中,Ω11是第1行第1列分块矩阵,Ω12是第1行第2列分块矩阵,Ω13是第1行第3列分块矩阵,Ω22是第2行第2列分块矩阵,Ω33是第3行第3列分块矩阵,是第1行第1列分块矩阵,是第1行第2列分块矩阵,是第1行第3列分块矩阵,是第1行第4列分块矩阵,L15是第1行第5列分块矩阵,L16是第1行第6列分块矩阵,L22是第2行第2列分块矩阵,L33是第3行第3列分块矩阵,L44是第4行第4列分块矩阵,L55是第5行第5列分块矩阵,L66是第6行第6列分块矩阵,是第1行第1列分块矩阵,G12是第1行第2列分块矩阵,G14是第1行第4列分块矩阵,是第1行第5列分块矩阵,G22是第2行第2列分块矩阵,G24是第2行第4列分块矩阵,是第2行第6列分块矩阵,是第2行第7列分块矩阵,是第2行第8列分块矩阵,G33是第3行第3列分块矩阵,G39是第3行第9列分块矩阵,是第4行第10列分块矩阵,是第4行第4列分块矩阵,是第5行第5列分块矩阵,是第6行第6列分块矩阵,是第7行第7列分块矩阵,是第8行第8列分块矩阵,是第9行第9列分块矩阵,是第10行第10列分块矩阵, 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, is the 1st row and 1st column block matrix, is a block matrix with row 1 and column 2, is a block matrix with row 1 and column 3, 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, 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, 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, is the 2nd row and 6th column block matrix, is the 2nd row and 7th column block matrix, 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, is the 4th row and 10th column block matrix, is the 4th row and 4th column block matrix, is the 5th row and 5th column block matrix, is the 6th row and 6th column block matrix, is the 7th row and 7th column block matrix, is the 8th row and 8th column block matrix, is the 9th row and 9th column block matrix, is the 10th row and 10th column block matrix,
分别为Dk,Kk,Et,k,Ct,k,ΔAk,Hk,ΔBk,ΔAk,Ek,Kk,Ck,R3k的转置,为定义的左右区间的第一号矩阵,为定义的左右区间的第二号矩阵,为定义的左右区间的第三号矩阵,为在第k时刻的非线性激励函数;C1k为在第k时刻第一个分量已知系统的噪声分布矩阵,C2k为在第k时刻第二个分量已知系统的噪声分布矩阵,Hk为在第k时刻的已知测量的调节矩阵;Dk为在第k时刻的已知测量的量度矩阵;表示的是ι=1到k+1求和的值,表示在第k时刻已知信道衰减矩阵,diag{·}表示的是对角矩阵,是第k时刻已知信道衰减矩阵,m表示的是第m个信道,ρ∈(0,1)为已知的调节正常数;Sk为在第k时刻的误差协方差矩阵的上界;Θ1k T, 分别为Θ1k,C1k,Φι,Ct,k,Et,k的转置;ζ为调节系数,Ω(Sk)为在第k时刻求解出的上界矩阵;Sk-d为在第k-d时刻的误差协方差矩阵的上界矩阵;tr(Sk)为在第k时刻的误差协方差矩阵上界的迹;tr()为矩阵的迹,I为单位矩阵;为在第k时刻的第一号权重矩阵;为在第k时刻的第二号权重矩阵;为在第k时刻的第三号权重矩阵;是在第k时刻R3k的转置;是第1个分量在k时刻的已知适当维数的第一号实矩阵,是第2个分量在k时刻的已知适当维数的第二号实矩阵;为在第k时刻的非线性激励函数的状态估计;是第1个分量在k时刻的已知适当维数的第一号量度矩阵;是第2个分量在k时刻的已知适当维数的第二号量度矩阵;是第3个分量在k时刻的已知适当维数的第三号量度矩阵;是第4个分量在k时刻的已知适当维数的第三号量度矩阵;N5是第5个分量在k时刻的已知适当维数的第三号量度矩阵;M1,M2,M3,M4和M5分别是第一号,第二号,第三号,第四号和第五号的量度矩阵,为在第k时刻的神经元状态估计,为第k时刻的半正定矩阵;为第k时刻的半正定矩阵;为第k-d时刻的半正定矩阵;为在第k+1时刻的第一更新矩阵,Sk为估计误差的上界矩阵,tr(Sk)为在第k时刻估计误差上界矩阵Sk的迹;Sk-d为在k-d时刻的上界矩阵,κ为调节的权重系数,和均为已知的实值权重矩阵,是未知矩阵且满足 是的转置,γ为给定的正标量;为给定的半正定矩阵一号;分别是Ω12,Ω13,的转置; 分别是,G12,G14,G24,G410,的转置;分别是M1,M2,M3,M4,M5的转置;N1,N2,N3,N4,N5分别是的转置;和分别是第一个、第二个、第三个、第四个和第五个相关比例系数,0表示的是矩阵块的元素均为0。They are D k , K k , E t,k , C t,k , ΔA k , H k , ΔB k , ΔA k , The transpose of E k , K k , C k , R 3k , is the first matrix of the defined left and right intervals, is the second matrix of the defined left and right intervals, is the third matrix of the defined left and right intervals, 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; It represents the sum of ι=1 to k+1. represents the known channel attenuation matrix at the kth time, diag{·} represents the diagonal matrix, 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; Θ 1k T , They are Θ 1k , 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; is the first weight matrix at the kth moment; is the second weight matrix at the kth moment; is the third weight matrix at the kth moment; is the transpose of R 3k at the kth moment; is the first real matrix of known appropriate dimension of the first component at time k, is the second real matrix of known appropriate dimension of the second component at time k; is the state estimate of the nonlinear activation function at the kth moment; is the first metric matrix of known appropriate dimension of the first component at time k; is the second metric matrix of known appropriate dimension of the second component at time k; is the third metric matrix of known appropriate dimension of the third component at time k; 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. is the estimated neuron state at the kth moment, is the semi-positive definite matrix at the kth moment; is the semi-positive definite matrix at the kth moment; is the semi-positive definite matrix at the kdth moment; 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, and are all known real-valued weight matrices, is an unknown matrix and satisfies yes The transpose of , γ is a given positive scalar; For the given semi-positive definite matrix one; They are Ω 12 , Ω 13 , The transpose of They are G12 , G14 , G24 , G410 , The transpose of 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 The transpose of and 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:
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, C 2k = [-0.2-0.1sin(3k)] T ,
被控输出调节矩阵为:The controlled output adjustment matrix is:
Hk=[-0.01-0.01sin(2k)]H k =[-0.01-0.01sin(2k)]
状态权重矩阵为:The state weight matrix is:
非线性函数的权重矩阵和调节参数为:The weight matrix and adjustment parameters of the nonlinear function are:
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:
式中,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,半正定矩阵一号上界矩阵{Ωk}1≤k≤N=diag{0.2,0.2}和协方差初始状态通道参数Ct,s=0.38和Et,s=0.12,分别取传感器到中继信道和中继到估计器信道的噪声协方差和 Perturbation attenuation level γ = 0.7, semi-positive definite matrix No. The upper bound matrix {Ω k } 1≤k≤N = diag{0.2,0.2} and the covariance Initial state 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. and
利用递推线性矩阵不等式,求解线性矩阵不等式(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:
情形二(Case II):Case II:
状态估计器效果: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)
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)
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)
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 |
-
2022
- 2022-12-06 CN CN202211559637.5A patent/CN116227324B/en active Active
Patent Citations (4)
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)
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 |