Research Article Open Access
Investigation of Stabilities and Instabilities at Tokamak Plasma Behavior and Machine Learning with Big Data
Danilo Rastovic*
Danilo Rastovic, Control Systems Group, Cetingradska 23, 10110 Zagreb, Croatia.
*Corresponding author: Danilo Rastovic, Control Systems Group, Cetingradska 23, 10110 Zagreb, Croatia; E-mail: @
Received: September 29, 2017; Accepted: October 27, 2017; Published: November 14, 2017
Citation: Danilo Rastovic (2017) Investigation of Stabilities and Instabilities at Tokamak Plasma Behavior and Machine Learning with Big Data. Int J Mol Theor Phy 1(1): 1-5.
Abstract Top
We investigate the problem of stability and instability at tokamak plasma behavior. Generally, Jaynes maximum entropy method and Bayesian decision can be applied for recognizing the shape of the plasma.In the case of the power law behavior and the instabilities of plasma we introduce a new method. The maximization of mathematical expectations for events and fuzzy entropy is used for applications of fuzzy Bayesian neural networks for optimization and simulation without assumption on recurrence. In this case, it is possible to consider also the non-Gibbsian probability distribution functions with the power law case. The new calibration method for the non-equilibrium systems has been given.

Keywords: Stability; Instability; Tokamak; Bayesian learning;
Introduction
Whether driven by magnetic or electrostatic turbulence, anomalous transport causes a large increase in the energy losses, and poses severe constraint on the amount of auxiliary power required to heat up the plasma to the high temperatures required for thermonuclear ignition.

The possible way to simulate the turbulence is given by a system of equations. We can consider all possible types of motions: regular motion, Gaussian diffusion and Levy walk diffusion. At non-equilibrium systems the entropy can decrease with the step of the time.

We introduce a new kind of entropy that is always monotonic and increasing as it is the case with the classical entropy for the equilibrium systems.

In many cases actual plasma are far from thermal equilibrium, being inhomogeneous or having non-Maxwellian velocity distributions, features which significantly complicate the analysis of experimental data.

Understanding and predicting turbulent transport is a key issue on the way toward commercially, viable fusion reactors. Indeed, turbulence controls the confinement control of any magnetically confined plasmas.

Externally applied, non-axisymmetric magnetic fields form the basis of several relatively simple and direct methods to control Magnetohydrodynamic (MHD) instabilities in a tokamak, and most present and planned tokamaks now include a set of non-axisymmetric control coils for application of fields with low toroidal mode numbers. Saturated magnetic islands are manipulated directly with non-axisymmetric fields [1].
The Properties of Turbulence
Dynamical system can reach ‘spontaneously’ a steady state exibiting some anologies with a thermodynamic system at a critical point. The system is driven to a stationary state where avalanches are distributed according to a power law. This phenomenon is called Self-Organized-Criticality (SOC) [2].

The development of fault diagnosis and sensor validation technique must face: the non-linearity of the process, the uncertain knowledge about the phenomena, often expressed in a linguistic fashion. The central problems are: plasma position, current and shape control.

The variational principle states that the entropy is maximal at equilibrium. There is no ‘natural ‘ Gibbs distribution nor free energy in SOC system. It is not possible to analyze the out-of-equilibrium stationary state via a local equilibrium hypothesis where one decomposes the system into mesoscopic cells , locally at equilibrium.

Step- by-step new regimes of plasma operation have been identified, wherely turbulence can be externally controlled, which led to better and better confinement. Turbulence as a whole is not an equilibrium phenomenon. In the framework of mesoscopic non-equilibrium thermodynamics has been obtained a generalized Fokker-Planck equation, incorporating memory effects through time-dependent coefficients.

Anomalous transport has been found where turbulent processes of various kind are responsable for enhanced transport. Fluctuations are in many cases driven by nonlinear instability mechanisms, not by the linear instability mechanisms. Turbulent states do not satisfy the equipartition law. Therefore, even if plasmas are in thermodynamical equilibrium, the conventional statistical theory may not be sufficient. One formal way of the statistical theory is to extend the definition of entropy. We consider the fuzzy entropy.

For the most complicated cases of long interactions of the tokamak plasma, the approach with long memory effects can be applied via the theory of fuzzy neural Bayesian networks. The goal is, with localized control, to overcome the problems of non-equilibrium in plasma. In the case of stochastic turbulence, we define the fuzzy information entropy by S= i=0 j=0 p ij Δ ln p ij Δ       (1) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaiabg2 da9iabgkHiTmaaqahabaWaaabCaeaacaWGWbWaa0baaSqaaiaadMga caWGQbaabaGaeyiLdqeaaaqaaiaadQgacqGH9aqpcaaIWaaabaGaey OhIukaniabggHiLdaaleaacaWGPbGaeyypa0JaaGimaaqaaiabg6Hi LcqdcqGHris5aOGaciiBaiaac6gacaWGWbWaa0baaSqaaiaadMgaca WGQbaabaGaeyiLdqeaaOGaaeiiaiaabccacaqGGaGaaeiiaiaabcca caqGGaGaaeikaiaabgdacaqGPaaaaa@5619@ Where S is the spatiotemporal fuzzy Shannon Ent,

p ij Δ MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadchadaqhaa WcbaGaamyAaiaadQgaaeaacqGHuoaraaaaaa@3A0B@ is the fuzzy component that has a membership function.

As final criteria we take the maximal amount of mathematical expectations over the tree of events. Non-extensive statistical mechanics appears as a powerfull way to describe complex systems. With S(i) we shall denote for each time step i, the fuzzy spatial entropy, S(i)= j=0 p ij Δ ln p ij Δ      (2) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaiaacI cacaWGPbGaaiykaiabg2da9iabgkHiTmaaqahabaGaamiCamaaDaaa leaacaWGPbGaamOAaaqaaiabgs5aebaaaeaacaWGQbGaeyypa0JaaG imaaqaaiabg6HiLcqdcqGHris5aOGaciiBaiaac6gacaWGWbWaa0ba aSqaaiaadMgacaWGQbaabaGaeyiLdqeaaOGaaeiiaiaabccacaqGGa GaaeiiaiaabccacaqGOaGaaeOmaiaabMcaaaa@515D@ Turbulent transport in tokamak is very challenging scientific problem. Our goal is to develop quantitative predictions of tokamak turbulence and to use them in design of future fusion reactors.

We can consider the control system of the linearized Vlasov-Poisson-Fokker-Planck equations f t = v x f( E + v × B + M ) v f MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaalaaabaGaey OaIyRaamOzaaqaaiabgkGi2kaadshaaaGaeyypa0JaeyOeI0YaaCbi aeaacaWG2baaleqabaGaeyOKH4kaaOGaeyyXICTaey4bIe9aaSbaaS qaaiaadIhaaeqaaOGaamOzaiabgkHiTiaacIcadaWfGaqaaiaadwea aSqabeaacqGHsgIRaaGccqGHRaWkdaWfGaqaaiaadAhaaSqabeaacq GHsgIRaaGccqGHxdaTdaWfGaqaaiaadkeaaSqabeaacqGHsgIRaaGc cqGHRaWkdaWfGaqaaiaad2eaaSqabeaacqGHsgIRaaGccaGGPaGaey yXICTaey4bIe9aaSbaaSqaaiaadAhaaeqaaOGaamOzaaaa@5DEC@ df dt =Af+Bu,u=Ff+w,w= F 1 f+ w 1      (3) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGKbGaamOzaaqaaiaadsgacaWG0baaaiabg2da9iaadgeacaWGMbGa ey4kaSIaamOqaiaadwhacaGGSaGaamyDaiabg2da9iaadAeacaWGMb Gaey4kaSIaam4DaiaacYcacaWG3bGaeyypa0JaamOramaaBaaaleaa caaIXaaabeaakiaadAgacqGHRaWkcaWG3bWaaSbaaSqaaiaaigdaae qaaOGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGOaGaae4maiaa bMcaaaa@52BD@ Where A is Af=(v x )f,BFf= v (( E σ βv)f)  and MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyqaiaadA gacqGH9aqpcqGHsislcaGGOaGaamODaiabgEGirpaaBaaaleaacaWG 4baabeaakiaacMcacaWGMbGaaiilaiaadkeacaWGgbGaamOzaiabg2 da9iabgkHiTiabgEGirpaaBaaaleaacaWG2baabeaakiaacIcacaGG OaGaamyramaaBaaaleaacqaHdpWCaeqaaOGaeyOeI0IaeqOSdiMaam ODaiaacMcacaWGMbGaaiykaiaabccacaqGGaGaaeyyaiaab6gacaqG Kbaaaa@5553@ B F 1 f=σ v v f        (4) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOqaiaadA eadaWgaaWcbaGaaGymaaqabaGccaWGMbGaeyypa0Jaeq4WdmNaeyyX ICTaey4bIe9aaSbaaSqaaiaadAhaaeqaaOGaeyyXICTaey4bIe9aaS baaSqaaiaadAhaaeqaaOGaeyyXICTaamOzaiaabccacaqGGaGaaeii aiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeikaiaabsdacaqGPa aaaa@5089@ E(t, x) is the force field acting on particles. The strictly positive parameters β MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqOSdigaaa@3796@ and σ MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdmhaaa@37B8@ model a certain type of interaction between particles. The control of plasma transport, laminar and turbulent has been investigated.

We have the optimal control of tokamak resistive wall mode in the presence of noise. The effect of the finitely conducting wall is to reduce the ideal mode to a slow growing resistive wall mode, which needs to be feedback stabilized for high performance. Instead of “state” feedback one may consider use of multiple independent simple single mode (N) controllers acting on spatially filtered several N modes. However, these are often coupled at tokamak.

We can introduce feature regularization during feature selection for value function approximation. The smoothness prior is effective in the feature selection setting and present closed form smoothness, regularizers for the Fourier base [3].

We have general formalism of optimal control. Optimal feedback design that is based on deterministic model, is minimization of total fluctuation energy of the instabilities, as well as minimization of control power J( t f )= 1 2 0 t f [ f N T (t)Q f N (t)+ u T (t)Ru(t) ] dt[  4  ]      (5) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadQeacaGGOa GaamiDamaaBaaaleaacaWGMbaabeaakiaacMcacqGH9aqpdaWcaaqa aiaaigdaaeaacaaIYaaaamaapehabaWaamWaaeaacaWGMbWaa0baaS qaaiaad6eaaeaacaWGubaaaOGaaiikaiaadshacaGGPaGaeyyXICTa amyuaiabgwSixlaadAgadaWgaaWcbaGaamOtaaqabaGccaGGOaGaam iDaiaacMcacqGHRaWkcaWG1bWaaWbaaSqabeaacaWGubaaaOGaaiik aiaadshacaGGPaGaeyyXICTaamOuaiabgwSixlaadwhacaGGOaGaam iDaiaacMcaaiaawUfacaGLDbaaaSqaaiaaicdaaeaacaWG0bWaaSba aWqaaiaadAgaaeqaaaqdcqGHRiI8aOGaamizaiaadshadaWadaqaaa baaaaaaaaapeGaaeiiaiaaisdacaqGGaaapaGaay5waiaaw2faaiaa bccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqG1aGaae ykaaaa@6C51@ For Lagrangian L( f N , f N ' )= f N T Q f N (t)+ u T (t)Ru(t)+ λ T (t)[ A f N (t)+Bu(t) f N ' (t) ]      (6) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKaaGiaadYeaca GGOaGaamOzaOWaaSbaaKqaGeaacaWGobaabeaajaaicaGGSaGaamOz aOWaa0baaKqaGeaacaWGobaabaGaai4jaaaajaaicaGGPaGaeyypa0 JaamOzaOWaa0baaKqaGeaacaWGobaabaGaamivaaaajaaicqGHflY1 caWGrbGaeyyXICTaamOzaOWaaSbaaKqaGeaacaWGobaabeaajaaica GGOaGaamiDaiaacMcacqGHRaWkcaWG1bGcdaahaaqcbasabeaacaWG ubaaaKaaGiaacIcacaWG0bGaaiykaiabgwSixlaadkfacqGHflY1ca WG1bGaaiikaiaadshacaGGPaGaey4kaSIaeq4UdWMcdaahaaqcbasa beaacaWGubaaaKaaGiaacIcacaWG0bGaaiykaOWaamWaaKaaGeaaca WGbbGaamOzaOWaaSbaaKqaGeaacaWGobaabeaajaaicaGGOaGaamiD aiaacMcacqGHRaWkcaWGcbGaamyDaiaacIcacaWG0bGaaiykaiabgk HiTiaadAgakmaaDaaajeaibaGaamOtaaqaaiaacEcaaaqcaaIaaiik aiaadshacaGGPaaacaGLBbGaayzxaaGccaqGGaGaaeiiaiaabccaca qGGaGaaeiiaiaabccacaqGOaGaaeOnaiaabMcaaaa@78DE@ It holds that L f N d dt ( L f N ' )=0      (7) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaalaaabaGaey OaIyRaamitaaqaaiabgkGi2kaadAgadaWgaaWcbaGaamOtaaqabaaa aOGaeyOeI0YaaSaaaeaacaWGKbaabaGaamizaiaadshaaaGaaiikam aalaaabaGaeyOaIyRaamitaaqaaiabgkGi2kaadAgadaqhaaWcbaGa amOtaaqaaiaacEcaaaaaaOGaaiykaiabg2da9iaaicdacaqGGaGaae iiaiaabccacaqGGaGaaeiiaiaabccacaqGOaGaae4naiaabMcaaaa@4E57@ with optimal control u(t)= R 1 B T λ(t)     (8) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadwhacaGGOa GaamiDaiaacMcacqGH9aqpcqGHsislcaWGsbWaaWbaaSqabeaacqGH sislcaaIXaaaaOGaeyyXICTaamOqamaaCaaaleqabaGaamivaaaaki abgwSixlabeU7aSjaacIcacaWG0bGaaiykaiaabccacaqGGaGaaeii aiaabccacaqGGaGaaeikaiaabIdacaqGPaaaaa@4D4C@ which is easily obtained from Hamiltonian H. The Hamiltonian represents the total energy of the system, which is the sum of kinetic and potential energy. According to the Pontrjagin’s principle. H u =0 MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq GHciITcaWGibaabaGaeyOaIyRaamyDaaaacqGH9aqpcaaIWaaaaa@3C58@

When a magnetically confined plasma is heated strongly and a threshold heating power level is exited, it may make spontaneously transition from a low-confinement (or L-mode) state to a high confinement (or H-mode) state. In the H-mode, the energy confinement time is significantly enhanced i.e., typically by a factor of 2 or more. H-mode profiles have a characteristic edge pedestal. Improved confinement regime, or H-mode,first discovered by Fritz Wagner on ASDEX tokamak, is primarily due to the formation of an insulating region at the plasma edge, where step pressure gradients can form. A feature of H-mode operation is a series of explosive plasma events, called Edge Localized Modes or ELMs. There is the present understanding of these ELMs in terms of ideal magneto-hydrodynamic instabilities.

The Magneto-Hydrodynamic (MHD) stability of burning plasma is a key issue for operations of ITER. One key of MHD instability that can affected by alpha particles is the sawtooth oscilation. The prevention of core impurity accumulation is to apply Neutral Beam Injection (NBI) heating in the opposite direction to the plasma current.

At the ASDEX tokamak experiment, a new feedback control loop is under construction with the aim of stabilizing MHD instabilities, such as neoclassical tearing modes and sawteeth. At the sawtooth behavior in the comiter – NBI regime is attributed to a subtle balance of the competing stabilization from the toroidal rotation and destabilization from the presence of energetic ions [5]. Broadening the internal transport barriers can improve MHD stability.

Stochasticity of magnetic field lines is paramount to understanding particle transport and mixing. In the paper is studied magnetic field separation in turbulent plasmas, which feature so-called super-diffusion, where field lines separate faster than diffusively [6].
The Stability of Equilibrium
The starting point of theoretical treatment of plasma confinement in tokamak is the equilibrium configuration. Any magnetic field that breaks the exact axisymmetric of the equilibrium tokamak field is a non-integrable perturbation [7]. The KAM (Kolmogorov-Arnold-Moser) theory predicts that, for those irrational surfaces with safety factors sufficiently far from rational m/n, the topology is preserved and the surfaces are only slightly deformed from the unperturbed tori (KAM surfaces). Similar effects can be obtained by the theory of Landau damping.

The equilibrium code solves the Grad-Shafranov equation : Δ * ψ= μ 0 R j ϕ     (9) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiabgs5aenaaCa aaleqabaGaaiOkaaaakiabeI8a5jabg2da9iabgkHiTiabeY7aTnaa BaaaleaacaaIWaaabeaakiaadkfacqGHflY1caWGQbWaaSbaaSqaai abew9aMbqabaGccaqGGaGaaeiiaiaabccacaqGGaGaaeikaiaabMda caqGPaaaaa@4905@ j ϕ =R dp dψ G μ 0 R dG dψ      (10) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadQgadaWgaa WcbaGaeqy1dygabeaakiabg2da9iabgkHiTiaadkfadaWcaaqaaiaa dsgacaWGWbaabaGaamizaiabeI8a5baacqGHsisldaWcaaqaaiaadE eaaeaacqaH8oqBdaWgaaWcbaGaaGimaaqabaGccaWGsbaaamaalaaa baGaamizaiaadEeaaeaacaWGKbGaeqiYdKhaaiaabccacaqGGaGaae iiaiaabccacaqGGaGaaeikaiaabgdacaqGWaGaaeykaaaa@4FB0@ where ψ MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiYdKhaaa@37C3@ is the poloidal flux function, j ϕ MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOAamaaBa aaleaacqaHvpGzaeqaaaaa@38D8@ is the toroidal current density, p is the pressure, and G=R B ϕ MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadEeacqGH9a qpcaWGsbGaamOqamaaBaaaleaacqaHvpGzaeqaaaaa@3B09@ is the toroidal field function.

The ‘standard model’ of classical plasma physics is the Vlasov-Poisson-Fokker-Planck equation, here written with periodic boundary conditions and in dimensionless unit f t +v x f+M| f | v f= logΛ 2πΛ Q (f.f) L      (11) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaalaaabaGaey OaIyRaamOzaaqaaiabgkGi2kaadshaaaGaey4kaSIaamODaiabgEGi rpaaBaaaleaacaWG4baabeaakiaadAgacqGHRaWkcaWGnbWaaqWaae aacaWGMbaacaGLhWUaayjcSdGaey4bIe9aaSbaaSqaaiaadAhaaeqa aOGaamOzaiabg2da9maalaaabaGaciiBaiaac+gacaGGNbGaeu4MdW eabaGaaGOmaiabec8aWjabfU5ambaacaWGrbWaaSraaSqaaiaadYea aeqaaOGaaiikaiaadAgacaGGUaGaamOzaiaacMcacaqGGaGaaeiiai aabccacaqGGaGaaeiiaiaabIcacaqGXaGaaeymaiaabMcaaaa@5E3F@ where f = f(t, x, v) is the electron distribution function.

The nonlinear Landau damping for general interactions has been obtained by C.Mouhot and C.Villani [8].

The unique solution of the nonlinear Vlasov equation satisfies ρ(t,) ρ C r ( T d ) Cδ e 2π λ ' | t |     (12) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiablwIiqjabeg 8aYjaacIcacaWG0bGaaiilaiabgwSixlaacMcacqGHsislcqaHbpGC daWgaaWcbaGaeyOhIukabeaakiablwIiqnaaBaaaleaacaWGdbWaaW baaWqabeaacaWGYbaaaSGaaiikaiaadsfadaahaaadbeqaaiaadsga aaWccaGGPaaabeaakiabgsMiJkaadoeacqGHflY1cqaH0oazcqGHfl Y1caWGLbWaaWbaaSqabeaacqGHsislcaaIYaGaeqiWdaNaeq4UdW2a aWbaaWqabeaacaGGNaaaaSWaaqWaaeaacaWG0baacaGLhWUaayjcSd aaaOGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqGXaGaaeOmaiaa bMcaaaa@6213@ where ρ(t,x)= f(t,x,v)dv , MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdiNaai ikaiaadshacaGGSaGaamiEaiaacMcacqGH9aqpdaWdbaqaaiaadAga caGGOaGaamiDaiaacYcacaWG4bGaaiilaiaadAhacaGGPaGaamizai aadAhaaSqabeqaniabgUIiYdGccaGGSaaaaa@47E8@ ρ = f i (x,v)dvdx    (13) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aaS baaSqaaiabg6HiLcqabaGccqGH9aqpdaWdcaqaaiaadAgadaWgaaWc baGaamyAaaqabaaabeqab0Gaey4kIiVaey4kIipakiaacIcacaWG4b GaaiilaiaadAhacaGGPaGaamizaiaadAhacaWGKbGaamiEaiaabcca caqGGaGaaeiiaiaabccacaqGOaGaaeymaiaabodacaqGPaaaaa@4D45@ Again, here we have some kind of exponential behavior.

Let us consider fuzzy entropy for controlled fuzzy stochastic process S= i=0 j=0 p ij Δ ln p ij Δ     (14) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaiabg2 da9iabgkHiTmaaqahabaWaaabCaeaacaWGWbWaa0baaSqaaiaadMga caWGQbaabaGaeyiLdqeaaaqaaiaadQgacqGH9aqpcaaIWaaabaGaey OhIukaniabggHiLdaaleaacaWGPbGaeyypa0JaaGimaaqaaiabg6Hi LcqdcqGHris5aOGaciiBaiaac6gacaWGWbWaa0baaSqaaiaadMgaca WGQbaabaGaeyiLdqeaaOGaaeiiaiaabccacaqGGaGaaeiiaiaabIca caqGXaGaaeinaiaabMcaaaa@558A@ Generally, even for non-equilibrium systems, the maximum of fuzzy entropy should be obtained without techniques of Jaynes’ constraints, with Bayesian learning rules for infinite fuzzy logic controller

IF p( x i ) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiCaiaacI cacaWG4bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaaaa@3A64@ is E 1(i) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaaIXaGaaiikaiaadMgacaGGPaaabeaaaaa@39ED@ AND p( x i+1 / x i ) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiCaiaacI cacaWG4bWaaSbaaSqaaiaadMgacqGHRaWkcaaIXaaabeaakiaac+ca caWG4bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaaaa@3ED5@ is E 2(i) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaaIYaGaaiikaiaadMgacaGGPaaabeaaaaa@39EE@ THEN a i MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGPbaabeaaaaa@37F5@ is U (i) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyvamaaBa aaleaacaGGOaGaamyAaiaacMcaaeqaaaaa@3942@ (i=0, 1, 2,…..) (15) [9].

The goal is, with localized control, overcome the problem of non-equilibrium in plasma. The aim is converting the anomalous diffusion to the exponentially stabilizable diffusion, when it is possible [10].

We also show the passage from the non-Markovian to the Markovian behaviour in the normal diffusion regime using memory effects, appropriate sensors and actuators. Many current models for long-memory are non-Markovian processes. Some processes can be Markovian.

The Monte Carlo algorithm we are concerned with can be understood as a numerical method of solving the drift kinetic equation . Fuzzy probability theory, a theory in the process of development, is just standard probability theory with the generalized concept of random variables

More investigations are needed to establish the nature of instabilities and driving mechanism involved. For a detailed analysis of turbulence, adequate analysis tools are required. In the core of tokamak plasmas, the turbulence is believed to govern transport processes.

It is shown that for some cases, the control system of the linearized transport equation with a controllable unstable part is exponentially stabilizable [11]. As generalization, the asymptotic stability of a certain class of integrated semigroups is discussed by means of Lyapunov functionals [12]. In this case, we obtain the exponentially bounded behavior in the sense that

S(t)x K 1 t N exp(ωt) | x | X MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaauWaaeaaca WGtbGaaiikaiaadshacaGGPaGaamiEaaGaayzcSlaawQa7amaaBaaa leaaaeqaaOGaeyizImQaam4samaaCaaaleqabaGaeyOeI0IaaGymaa aakiaadshadaahaaWcbeqaaiaad6eaaaGcciGGLbGaaiiEaiaaccha caGGOaGaeqyYdCNaamiDaiaacMcadaabdaqaaiaadIhaaiaawEa7ca GLiWoadaWgaaWcbaGaamiwaaqabaaaaa@5002@ for t0 MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiDaiabgw MiZkaaicdaaaa@396E@ and xX MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiEaiabgI GiolaadIfaaaa@3953@ (where K denotes n factorial) (16)

The possible interconnection between of the equilibrium and of the anomalous transport has been discussed.

A nice example of construction of fuzzy probabilistic distribution and the fuzzy entropy has been given in [13].

There is a method for solving the Fokker-Planck equation in tokamak plasmas [14]. The code provides a set of powerful tools for studying radio frequency wave heating and current drive in tokamak plasmas.
Fuzzy Bayesian Computation
Approximate Bayesian Computation is a family of computational techniques which offer an almost automated solution in situations where evaluation of the posterior likelihood is computationally prohibitive or whenever suitable likelihoods are not available [15]. We consider the fuzzy Bayesian computation. The basic idea is that using a representative enough summary statistic coupled with a small enough tolerance level ε should produce a good approximation of the posterior distribution.

Power laws are well known from the work by Pareto for the distribution of income. For example (short- time) particle travel behavior can be well described by a decreasing power law of distance x, that is, f(x)¨ x d MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaiaacI cacaWG4bGaaiykaiabg2Hi1kaacIkacaWG4bWaaWbaaSqabeaacqGH sislcaWGKbaaaaaa@3EE2@ with positive decay parameter d. An important characteristic of this power law is its slow convergence to zero, which in application enables long-range actions in addition to principal short-range actions. They describe the distribution of earthquake magnitudes and many other natural phenomena. An alternative spatial interaction function of an isotropic Gaussian kernel f(x)=exp( x 2 2 σ 2 ) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaiaacI cacaWG4bGaaiykaiabg2da9iGacwgacaGG4bGaaiiCaiaacIcacqGH sisldaWcaaqaaiaadIhadaahaaWcbeqaaiaaikdaaaaakeaacaaIYa Gaeq4Wdm3aaWbaaSqabeaacaaIYaaaaaaakiaacMcaaaa@44CF@ with scale parameter σ allows for occasional long-range actions : a power law.The shifted version known as Pareto type II has density kernel f(x)= (x+σ) d (1+ x σ ) d      (17a)  or MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaiaacI cacaWG4bGaaiykaiabg2da9iaacIcacaWG4bGaey4kaSIaeq4WdmNa aiykamaaCaaaleqabaGaeyOeI0Iaamizaaaakiabg2Hi1kaacIcaca aIXaGaey4kaSYaaSaaaeaacaWG4baabaGaeq4WdmhaaiaacMcadaah aaWcbeqaaiabgkHiTiaadsgaaaGccaqGGaGaaeiiaiaabccacaqGGa GaaeiiaiaabIcacaqGXaGaae4naiaabggacaqGPaGaaeiiaiaabcca caqGVbGaaeOCaaaa@549C@ f(x)= ( x 2 + σ 2 ) d (1+ ( x σ ) 2 ) d      (17b) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rq0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadAgacaGGOa GaamiEaiaacMcacqGH9aqpcaGGOaGaamiEamaaCaaaleqabaGaaGOm aaaakiabgUcaRiabeo8aZnaaCaaaleqabaGaaGOmaaaakiaacMcada ahaaWcbeqaaiabgkHiTiaadsgaaaGccqGHDisTcaGGOaGaaGymaiab gUcaRiaacIcadaWcaaqaaiaadIhaaeaacqaHdpWCaaGaaiykamaaCa aaleqabaGaaGOmaaaakiaacMcadaahaaWcbeqaaiabgkHiTiaadsga aaGccaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqGXaGaae 4naiaabkgacaqGPaaaaa@5552@ Extreme Learning Machine (ELM) is a new approach for learning the parameters of the hidden layers of a multilayer neural network (as a multilayer perceptron or the radial basis function neural network) . Its main advantage is the lower computational cost. The research on Bayesian methods for neural models has become very intense currently. These methods introduce a probability distribution on the network parameters and the committed errors. The Bayesian ELM has the advantages of both ELM and Bayesian models [16].

Bayesian linear regression is optimizing the weights of the output layer. Any Bayesian modeling is carried out in two steps:

1. Inference of the posterior distribution
2. Calculation of the output distribution

The Bayesian approach involves the use of some parameters (hyper-parameters) that allow regularization. The use of confidence intervals increases the reliability of a model’s output. A priori knowledge can be introduced in the models by means of error distribution and parameter distributions that must be defined when applying Bayes’ theorem. This knowledge can improve the performance of the model.

Normal distributions emerge naturally when formulating a central limit theorem, as solution of linear Fokker-Planck (FP) equations. Similarly, q-Gaussian distributions are related to generalized central limit theorem, are solutions of a kind of nonlinear FP equations and maximize non-extensive Tsallis entropy and a generalized second order moment constraints to obtain oscillatory behavior. There is the equivalent formulation of such phenomena inside of the theory of fuzzy Jaynes entropy.

Fuzzy controlled Markovian process can be obtained from general non-Markovian stochastic process via artificial fuzzy neural networks in some cases. Many current models for long-memory are non-Markovian. The turbulent flux constitutes a random process with significant deviations from Gaussian statistics. Detailed information concerning the turbulent fluxes can be obtained by conditioned averaging.

We understood entropy as a measure of uncertainity of Markov process. The “uncertainity of uncertainity” is unavoidable in non-equilibrium systems. It is the case which is obtained by fuzzy entropy. For controlled fuzzy Markov chain we use the following infinite fuzzy logic controller with Bayesian learning rules

IF P Δ ( x i ) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuamaaCa aaleqabaGaeyiLdqeaaOGaaiikaiaadIhadaWgaaWcbaGaamyAaaqa baGccaGGPaaaaa@3BE2@ c is E 1(i) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaaIXaGaaiikaiaadMgacaGGPaaabeaaaaa@39ED@ AND P Δ ( x i+1 / x i ) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuamaaCa aaleqabaGaeyiLdqeaaOGaaiikaiaadIhadaWgaaWcbaGaamyAaiab gUcaRiaaigdaaeqaaOGaai4laiaadIhadaWgaaWcbaGaamyAaaqaba GccaGGPaaaaa@4053@ is E 2(i) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaaIYaGaaiikaiaadMgacaGGPaaabeaaaaa@39EE@ THEN a i Δ MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaDa aaleaacaWGPbaabaGaeyiLdqeaaaaa@395D@ is U (i) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyvamaaBa aaleaacaGGOaGaamyAaiaacMcaaeqaaaaa@3942@ (i=0, 1, 2,…) (18)

As “defuzzification “ on G i = E 1(i) E 2(i) U (i) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4ramaaBa aaleaacaWGPbaabeaakiabg2da9iaadweadaWgaaWcbaGaaGymaiaa cIcacaWGPbGaaiykaaqabaGccqGHPiYXcaWGfbWaaSbaaSqaaiaaik dacaGGOaGaamyAaiaacMcaaeqaaOGaeyykICSaamyvamaaBaaaleaa caGGOaGaamyAaiaacMcaaeqaaaaa@4779@ we can take the following reason: the actions δu(k)=( a 0 , a 1 , a 2 ,...) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiTdqMaam yDaiaacIcacaWGRbGaaiykaiabg2da9iaacIcacaWGHbWaaSbaaSqa aiaaicdaaeqaaOGaaiilaiaadggadaWgaaWcbaGaaGymaaqabaGcca GGSaGaamyyamaaBaaaleaacaaIYaaabeaakiaacYcacaGGUaGaaiOl aiaac6cacaGGPaaaaa@46E7@ (19)

are given on such a way that the maximum of information fuzzy entropy S= i=0 j=0 p ij Δ ln p ij Δ MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaiabg2 da9iabgkHiTmaaqahabaWaaabCaeaacaWGWbWaa0baaSqaaiaadMga caWGQbaabaGaeyiLdqeaaaqaaiaadQgacqGH9aqpcaaIWaaabaGaey OhIukaniabggHiLdaaleaacaWGPbGaeyypa0JaaGimaaqaaiabg6Hi LcqdcqGHris5aOGaciiBaiaac6gacaWGWbWaa0baaSqaaiaadMgaca WGQbaabaGaeyiLdqeaaaaa@5032@ is obtained, where P i Δ = j=0 N i p ij Δ MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuamaaDa aaleaacaWGPbaabaGaeyiLdqeaaOGaeyypa0ZaaabCaeaacaWGWbWa a0baaSqaaiaadMgacaWGQbaabaGaeyiLdqeaaaqaaiaadQgacqGH9a qpcaaIWaaabaGaamOtamaaBaaameaacaWGPbaabeaaa0GaeyyeIuoa aaa@4596@ is possible uncertain clustering , for global optimization in time domain [17].

For actions a i MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGPbaabeaaaaa@37F5@ we can take fuzzy neural network to memorize the past (history) at given moment i.

For times longer than the relaxation times, the correlation function for anomalous diffusion becomes a power law for broad-band noise. Anomalous transport to be induced by the plasma turbulence caused various instabilities. In the fusion-related research using toroidal plasmas, the suppression of turbulence has been main problem to make confinement time longer.

Two major challenges remain in clustering methods. One is efficient sampling. The other major challenge is posterior inference. We can use Bayesian inference.

Let us consider an ill-posed inverse problem where a physical quantity needs to be recovered from an indirect measurement modeled by m=Af+ε MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBaiabg2 da9iaadgeacaWGMbGaey4kaSIaeqyTdugaaa@3C27@ , where ε MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdugaaa@379C@ is white noise taking values in R k MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOuamaaCa aaleqabaGaam4Aaaaaaaa@37E9@ .

In Bayesian inversion, the quantities m n , f (n) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBamaaBa aaleaacaWGUbaabeaakiaacYcacaWGMbWaaSbaaSqaaiaacIcacaWG UbGaaiykaaqabaaaaa@3C23@ and ε MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdugaaa@379C@ , ε MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdugaaa@379C@ where is white noise, taking values in R k MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOuamaaCa aaleqabaGaam4Aaaaaaaa@37E9@ . These quantities are considering as random variables. We use the fuzzy Bayes formula post Δ ( f (n) / m n )= Δ ( f n ) Δ (( m n )/ f (n) ) Δ ( m n )       (20) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaebmaeaaca GGOaGaamOzamaaBaaaleaacaGGOaGaamOBaiaacMcaaeqaaaqaaiaa dchacaWGVbGaam4CaiaadshaaeaacqGHuoara0Gaey4dIunakiaac+ cacaWGTbWaaSbaaSqaaiaad6gaaeqaaOGaaiykaiabg2da9maalaaa baGaey4dIu9aaWbaaSqabeaacqGHuoaraaGccaGGOaGaamOzamaaBa aaleaacaWGUbaabeaakiaacMcacqGHpis1daahaaWcbeqaaiabgs5a ebaakiaacIcacaGGOaGaamyBamaaBaaaleaacaWGUbaabeaakiaacM cacaGGVaGaamOzamaaBaaaleaacaGGOaGaamOBaiaacMcaaeqaaOGa aiykaaqaaiabg+GivpaaCaaaleqabaGaeyiLdqeaaOGaaiikaiaad2 gadaWgaaWcbaGaamOBaaqabaGccaGGPaaaaiaabccacaqGGaGaaeii aiaabccacaqGGaGaaeiiaiaabIcacaqGYaGaaeimaiaabMcaaaa@65B5@ Given a realization m MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaCbiaeaaca WGTbaaleqabaGaeS4SLyfaaaaa@38EF@ , after appropriate defuzzification, the maximum a posteriori estimate f (n) MAP MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaDa aaleaacaGGOaGaamOBaiaacMcaaeaacaWGnbGaamyqaiaadcfaaaaa aa@3BC6@ defined by post Δ ( f (n) MAP / m )=max[ post Δ ( f (n) / m): f (n) T n (x) ]      (21) MathType@MTEF@5@5@+= feaagGart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaebmaKaaGe aacaGGOaGaamOzaOWaa0baaKqaGeaacaGGOaGaamOBaiaacMcaaeaa caWGnbGaamyqaiaadcfaaaaabaGaamiCaiaad+gacaWGZbGaamiDaa qaaiabgs5aebqcdaIaey4dIunajaaicaGGVaGcdaWfGaqcaasaaiaa d2gaaKqaGeqabaGaey4jIKnaaKaaGiaacMcacqGH9aqpciGGTbGaai yyaiaacIhakmaadmaajaaibaGcdaqeWaqcaasaaiaacIcacaWGMbGc daWgaaqcbasaaiaacIcacaWGUbGaaiykaaqabaqcaaIaai4laOWaaC biaKaaGeaacaWGTbGaaiykaiaacQdacaWGMbGcdaWgaaqcbasaaiaa cIcacaWGUbGaaiykaaqabaqcaaIaeyicI4SaamivaOWaaSbaaKqaGe aacaWGUbaabeaajaaicaGGOaGaamiEaiaacMcaaKqaGeqabaGaey4j IKnaaaqaaiaadchacaWGVbGaam4CaiaadshaaeaacqGHuoaraKWaGi abg+GivdaajaaicaGLBbGaayzxaaGccaqGGaGaaeiiaiaabccacaqG GaGaaeiiaiaabccacaqGOaGaaeOmaiaabgdacaqGPaaaaa@70C7@ The Bayesian inversion generalizes deterministic regularization and yields information about the structure of ill posed inverse problem.

The plasma evolution solver involves more physics of uncertainity than the equilibrium solver [18].

The aim was to obtain the maximum of information fuzzy entropy by fuzzy Bayesian rules, that generalizes the Jaynes method for equilibriums to the case of non-equilibriums [19].

In a tokamak reactor, several magnetic fields have to be applied to confine the plasma. These are combined toroidal and poloidal magnetic fields. They are necessary to both: achieve the desired plasma configuration and to control the plasma shape and position [20].

With the methods of artificial neural networks with back-propagation it is possible recognize the solutions of Grad-Shafranov equations [21].
Conclusion
We have described how to determine the optimal a posteriori probability distribution and the appropriate actions. Examples of particularly popular estimation approaches include regularized variants of the Expectation-Maximization algorithm. We have explored the effects of specifying white noise (no memory), autoregressive (short memory) and long-memory correlation structures for error processes. The prior level allows for Bayesian inference by providing prior distributions for all unknown parameters in the data-and process- levels [22].

The Bayesian ,modeling framework is closely related to the stochastic filtering methods, for example classical Kalman filtering.

The Bayesian approach allows for all parameters to be estimated while avoiding the known estimation difficulties inherent to filtering. When estimating the spline coefficients, no information on levels or trends during the observation period was exchanged if it was possible [23]. In the case of changes in plasma evolution it has been shown that with Bayesian learning and maximum of fuzzy entropy can be reached all probabilistic distributions, including non-Gibbsians distributions [24, 25]. As reference papers concerning the diagnostic of high-temperature plasma conditions in tokamak reactors under the effect of plasma-wall interaction see [26-29].
ReferencesTop
  1. EJ Strait. Magnetic Control of Magnetohydrodynamic Instabilities in Tokamaks.  Phys. Plasmas. 2014;22(2). Doi: 10.1063/1.4902126
  2. B Cessac, Ph Blanchard, T Kruger and JL Meunier. Self-Organized Criticality and Thermodynamic Formalism. Jour. of Stat.Phys.2004;115(5-6):1283-1326.
  3. DS Wookey and GD Konidaris. Regularized feature selection in reinforcement learning.  Mach. Learn. 2015;100(2-3):655-676.
  4. AK Sen, M Nagashima and RW Longman. Optimal control of tokamak resistive wall modes in the presence of noise. Phys.Plasmas.2003;10(11).
  5. D Rastovic. Fuzzy Scaling and Stability of Tokamaks.  J.Fusion Energy. 2009;28(1):101-106.
  6. A Beresnyak. Asymmetric diffusion of magnetic field lines in turbulence. Phys.Scr. 2014;89(11).
  7. PR Garabedian. A unified theory of tokamaks and stellarators.  Comm.Pure Appl.Math. 1994;47(3):281-292.
  8. C Mouhot and C Villani. Landau damping,  J.Math.Phcys. 2010;51(1).
  9. D Rastovic. Targeting and Synchronization of Tokamak with Recurrent Artificial Neural Networks. J.Neural Comp.&Appl. 2012;21(5):1065-1069.
  10. D Rastovic. On Stochastic Control of Tokamak and Artificial Intelligence.  J.Fusion Energy. 2007;26(4):337-342.
  11. D. Rastovic. Feedback Stabilization of Some Classes of Nonlinear Transport Systems. Rend.Circ.. Mat. Palermo. 2002;51(2):325-332.
  12. D Rastovic. A Note on Stability Properties of Integrated Semigroups, Acta Math. Inf. Univ. Ostraviensis. 1995;3(1):61-65.
  13. A Kaufmann. Introduction to the Theory of Fuzzy Subsets, Academic Press. Inc. New York. 1975.
  14. Zheng Pingwei, Gong Xueyu, Yu Jun and Yu Jun. Fully Implicit Iterative Solving Method for  the Fokker-Planck Equation in Tokamak Plasmas, Plasma Sci.Technol. 2014;16(11):1000-1009.
  15. G Biau, F Cerou and A Guyader. New insights into Approximate Bayesian Computation, Annal.Inst.Henri Poincare, Prob.-Stat. 2015;51(1):376-403.
  16. E Soria-Olivas, Gómez-Sanchis J, Martín JD, Vila-Francés J, Martínez M and Magdalena JR, et al. BELM: Bayesian Extreme Learning Machine. IEEE Trans.Neural Networks. 2011;22(3):505-510. doi: 10.1109/TNN.2010.2103956
  17. D Rastovic. From Non-Markovian Processes to Stochastic Real Time Control for Tokamak Plasma Turbulence Via Artificial Intelligence Techniques.  J.Fusion Energy. 2015;34(2):207-215.
  18. FC Morabito. Equilibrium parameters recovery for experimental data in ASDEX upgrade elongated plasmas. Neural Networks, Proc.IEEE Int.Conf. 1995.
  19. D Rastovic. Infinite fuzzy logic controller and maximum entropy principle. InterStat, Virginia. 2003.
  20. G Ambrosino M Ariola, GD Tommasi and A pironti. Plasma Vertical Stabilization in the ITER Tokamak via Constrained Static Output Feedback. IEEEE Trans. Control Systems Technology. 2011;19(2):376-389.
  21. JB Lister and H Schnurrenberger. Fast non-linear extraction of plasma equilibrium parameters using a neural network mapping. Nucl.Fusion. 1991;31(7):1291-1300.
  22. Dong Li, Yi Liu, J. Svensson, YQ Liu, XM Song and Rui Mao, et al., Bayesian soft x-ray tomography and MHD mode analysis on HL-2A. Nucl.Fusion. 2016;56(3):8.
  23. D Rastovic. Approach with different entropies to the problem of plasma confinement. Journal of Physics& Astronomy. 2016;4(1):1-5.
  24. D Rastovic. Physics, electronics and entropies for plasma confinement. Journal of Physics & Astronomy. 2013;2(3):75-80.
  25. D Rastovic. Tokamak design as one sustainable system. Neural Network World. 2011;6(11):493-504.
  26. K Slabkowska, M Polasik, E Szymanska, J Starosta, L Syrocki and J Rzadkiewicz, et al. Modeling of the L and M  x-ray  line structures  for tungsten in high-temperature tokamak plasmas, Physica Scripta. 2014;T161.
  27. K Slabkowska, J Rzadkiewicz, L Syrocki, E Szymanska, A Shumack and M polasik, et al. On the interpretation of high-resolution x-ray spectra from JET with an ITER-like wall.  J. Phy B: Atom. Mol. Opt. 2015;48(14).
  28. K Slabkowska et al., The individual M x-ray line contributions originating from Cu- and Co-like tungsten for various plasma temperatures. Nucl. Instr. Meth. Phys. Res. B. 2017;408:265-270.
  29. Ł Syrocki, K Słabkowska, E Węder, J tarosta-Sztuczka and M Polasik. Modeling of the soft N, M and L x-ray lines from tungsten relevant to plasma parameters in the WEST tokamak.  Nucl. Instr. Meth. Phys. Res. B. 2017;408:257-264 (2017).
 
Listing : ICMJE   

Creative Commons License Open Access by Symbiosis is licensed under a Creative Commons Attribution 3.0 Unported License