摘要
利用指数加权在线核序列极限学习机(exponential weighted online sequential extreme learning machine with kernel, EW-KOSELM)辨识算法,开展了针对混沌动力学系统的动态重构研究. EW-KOSELM算法将核递归最小二乘(kernel recursive least squares, KRLS)算法直接延伸至在线ELM (extreme learning machine)框架中,通过引入遗忘因子削弱了旧数据的影响,并基于"固定预算(fixed-budget, FB)"内存技术,应对在线核学习算法所固有的规模不断增长的计算困难.将所提辨识算法应用于Duffing-Ueda振子的混沌动力学系统数值仿真实例中,对基于FB-EW-KOSELM的辨识模型与原系统的动态性能进行了定性与定量的分析校验,定性校验准则是基于对比辨识模型与原系统吸引子(轨迹嵌入)、庞加莱映射、分岔图、极限环完成的,定量校验准则包括对比辨识模型与原系统的李雅普诺夫指数与关联维.进一步将其分别应用于来自测量蔡氏电路产生双涡卷吸引子与螺旋吸引子的实测数据实验及某一实际混沌电路所产生的时间序列中,对于具有低信噪比的实测电压或电流数据还需进行了小波降噪预处理.通过分析辨识模型重构吸引子,实验结果表明,FB-EW-KOSELM算法具有良好的动态重构性能,能精确地再生出展示混沌动态行为的过程非线性模型,且具有与原混沌系统非常接近的动态不变性指标.
For the dynamic reconstruction of the chaotic dynamical system, a method of identifying an exponential weighted online sequential extreme learning machine with kernel(EW-KOSELM) is proposed. The kernel recursive least square(KRLS) algorithm is directly extended to an online sequential ELM framework, and weakens the effect of old data by introducing a forgetting factor. Meanwhile, the proposed algorithm can deal with the ever-increasing computational difficulties inherent in online kernel learning algorithms based on the‘fixed-budget' memory technique. The employed EW-KOSELM identification method is firstly applied to the numerical example of Duffing-Ueda oscillator for chaotic dynamical system based on simulated data, the qualitative and quantitative analysis for various validation tests of the dynamical properties of the original system as well as the identification model are carried out. A set of qualitative validation criteria is implemented by comparing chaotic attractors i.e. embedding trajectories, computing the corresponding Poincare mapping,plotting the bifurcation diagram, and plotting the steady-state trajectory i.e. the limit cycle between the original system and the identification model. Simultaneously, the quantitative validation criterion which includes computing the largest positive Lyapunov exponent and the correlation dimension of the chaotic attractors is also calculated to measure the closeness i.e. the approximation error between the original system and the identification model. The employed method is further applied to a practical implementation example of Chua's circuit based on the experimental data which are generated by sampling and recording the measured voltage across a capacitor, the inductor current from the double-scroll attractor, the measured voltage across a capacitor from the Chua's spiral attractor and an experimental time series from a chaotic circuit. The digital filtering technique is then used as a preprocessing approach, on the basis of wavelet denoising the measured data with lower signal-noise ratio(SNR) which can produce the double-scroll attractor or the spiral attractor,the reconstruction attractor of the identification model is compared with the reconstruction attractor from the experimental data for original system. The above experimental results confirm that the EW-FB-KOSELM identification method has a better performance of dynamic reconstruction, which can produce an accurate nonlinear model of process exhibiting chaotic dynamics. The identification model is dynamically equivalent or system approximation to the original system.
引文
[1]Chen G R,Dong X N 1998 From Chaos to Order:Methodologies,Perspectives and Applications(Singapore:World Scientific)pp21-134
[2]Strogatz S H 2015 Nonlinear Dynamics and Chaos:with Applications to Physics,Biology,Chemistry,and Engineering2nd(Boulder,CO:Westview press)pp244-411
[3]Aguirre L A,Letellier C 2009 Math.Probl.Eng.2009 1
[4]Ahmad T A,Sundarapandian V 2015 Chaos Modeling and Control Systems Design(Berlin:Springer-Verlag)pp59-72
[5]Feng J C 2008 Reconstruction of Chaotic Signals with Applications to Chaos-Based Communications(Singapore:World Scientific)pp27-164
[6]Mattera D,Haykin S 1999 Advances in Kernel Methods:Support vector learning(Cambridge:MIT Press)p211
[7]Ishii S,Sato M A 2001 Neural Networks 14 1239
[8]Zhang Z,Wang T,Liu X 2014 Neurocomputing 131 368
[9]Li R G,Zhang H L,Fan W H,Wang Y 2015 Acta Phys.Sin.64 200506(in Chinese)[李瑞国,张宏立,范文慧,王雅2015物理学报64 200506]
[10]Wang X Y,Han M 2015 Acta Phys.Sin.64 070504(in Chinese)[王新迎,韩敏2015物理学报64 070504]
[11]Wen S,Zeng Z,Huang T,Chen Y 2013 Phys.Lett.A 3772016
[12]Miranian A,Abdollahzade M 2013 IEEE Trans.Neural Netw.Learn.Syst.24 207
[13]Shaw P K,Saha D,Ghosh S,Janaki M S,Lyengar A S 2015Chaos Soliton.Fract.78 285
[14]Billings S A,Coca D 1999 Int.J.Bifurcat.Chaos 9 1263
[15]Sanchez L,Infante S 2013 Chil.J.Stat.4 35
[16]Aguirre L A,Teixeira B O S,Torres L A B 2005 Phys.Rev.E 72 026226
[17]Li J,Dong H Y 2008 Acta Phys.Sin.57 4756(in Chinese)[李军,董海鹰2008物理学报57 4756]
[18]Huang G B,Zhu Q Y,Siew C K 2006 Neurocomputing 70 489
[19]Liang N Y,Huang G B,Saratchandran P,Sundararajan N2006 IEEE Trans.Neural Netw.17 1411
[20]Scardapane S,Comminiello D,Scarpiniti M,Uncini A 2015IEEE Trans.Neural Netw.Learn.Syst.26 2214
[21]Vaerenbergh S V,Santamaria I,Liu W,Principe J C 2010The 35th IEEE International Conference on Acoustics Speech&Signal Processing Dallas,Texas,USA,March 14-19,2010p1882
[22]Liu W,Principe J C,Haykin S 2010 Kernel Adaptive Filtering:A Comprehensive Introduction pp108-110
[23]Engel Y,Mannor S,Meir R 2004 IEEE Trans.Signal Process.52 2275
[24]se Kruif B J,de Vries T J A 2003 IEEE Trans.Neural Netw.14 696
[25]Ueda Y 1985 Int.J.Nonlin.Mech.20 481
[26]Yang Z,Gao Y,Gao Y,Zhang J 2009 Chin.Phys.Lett.26060506
[27]Odavic J,Mali P,Tekic J,Pantic M,Pavkov-Hrvojevic M2017 Commun.Nonlinear Sci.Numer.Simul.47 100
[28]Suykens J A K,Vandewalle J 1995 IEEE Trans.Circ.Syst.Fund.Theor.Appl.42 499
[29]Cannas B,Cincotti S,Marchesi M,Pilo F 2001 Chaos Soliton.Fract.12 2109
[30]Aguirre L A,Rodrigues G G,Mendes E M 1997 Int.J.Bifurcat.Chaos 7 1411
[31]Kennedy M P 1992 Frequenz 46 66
[32]Timmer J,Rust H,Horbelt W,Voss H U 2000 Phys.Lett.A274 123