Nothing Special   »   [go: up one dir, main page]

Next Issue
Volume 23, May
Previous Issue
Volume 23, March
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
entropy-logo

Journal Browser

Journal Browser

Entropy, Volume 23, Issue 4 (April 2021) – 112 articles

Cover Story (view full-size image): T cells recognize pathogen-derived structures (antigens) through randomly recombined surface T-cell receptors (TCRs). T-cell development comprises extensive proliferation and sequential steps of quality control based on the affinity of TCRs for self-antigens, interrogated by successive interaction with antigen-presenting cells, leading to substantial cell death and, thus, generating a pool of T cells with a broad, but not self-reactive, TCR repertoire. Here, we review experimental and mathematical strategies to infer the dynamic properties of T-cell development in the thymus across multiple scales: cell cycle, population dynamics and their regulations, and how physiological T-cell development emerges from cellular interactions. View this paper
  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Section
Select all
Export citation of selected articles as:
14 pages, 2106 KiB  
Article
Subgraphs of Interest Social Networks for Diffusion Dynamics Prediction
by Valentina Y. Guleva, Polina O. Andreeva and Danila A. Vaganov
Entropy 2021, 23(4), 492; https://doi.org/10.3390/e23040492 - 20 Apr 2021
Viewed by 2361
Abstract
Finding the building blocks of real-world networks contributes to the understanding of their formation process and related dynamical processes, which is related to prediction and control tasks. We explore different types of social networks, demonstrating high structural variability, and aim to extract and [...] Read more.
Finding the building blocks of real-world networks contributes to the understanding of their formation process and related dynamical processes, which is related to prediction and control tasks. We explore different types of social networks, demonstrating high structural variability, and aim to extract and see their minimal building blocks, which are able to reproduce supergraph structural and dynamical properties, so as to be appropriate for diffusion prediction for the whole graph on the base of its small subgraph. For this purpose, we determine topological and functional formal criteria and explore sampling techniques. Using the method that provides the best correspondence to both criteria, we explore the building blocks of interest networks. The best sampling method allows one to extract subgraphs of optimal 30 nodes, which reproduce path lengths, clustering, and degree particularities of an initial graph. The extracted subgraphs are different for the considered interest networks, and provide interesting material for the global dynamics exploration on the mesoscale base. Full article
Show Figures

Figure 1

Figure 1
<p>The framework of the appropriate subgraphs search. Sampling methods are used for subgraph extraction, then motif extraction techniques are used for topological verification: Motif distributions for sub- and supergraphs are compared; the regression model on the base of sample motif distribution is used for transient time prediction for the SI diffusion model. The best subgraphs and corresponding methods for topological and functional criteria are evaluated.</p>
Full article ">Figure 2
<p>Mean squared error between motif distribution of extracted subgraphs and corresponding supergraphs for different sampling techniques.</p>
Full article ">Figure 3
<p>Coefficient of determination for regression model over subgraph <span class="html-italic">M</span> motif distribution, predicting <math display="inline"><semantics> <mrow> <mi>τ</mi> <mo>(</mo> <mi>G</mi> <mo>,</mo> <mi mathvariant="script">D</mi> <mo>)</mo> </mrow> </semantics></math> for different sampling techniques.</p>
Full article ">Figure 4
<p>Resulting subgraphs on the corresponding supergraphas and their single versions for the selected topics: (<b>a</b>) Counter strike, (<b>b</b>) Free donut—i like it!, (<b>c</b>) feminism, (<b>d</b>) introvert, (<b>e</b>) pizza, and (<b>f</b>) stop smoking. Node size correspond to its degree, blue nodes highlight subgraphs extracted.</p>
Full article ">
17 pages, 1401 KiB  
Article
A New Two-Stage Algorithm for Solving Optimization Problems
by Sajjad Amiri Doumari, Hadi Givi, Mohammad Dehghani, Zeinab Montazeri, Victor Leiva and Josep M. Guerrero
Entropy 2021, 23(4), 491; https://doi.org/10.3390/e23040491 - 20 Apr 2021
Cited by 34 | Viewed by 3563
Abstract
Optimization seeks to find inputs for an objective function that result in a maximum or minimum. Optimization methods are divided into exact and approximate (algorithms). Several optimization algorithms imitate natural phenomena, laws of physics, and behavior of living organisms. Optimization based on algorithms [...] Read more.
Optimization seeks to find inputs for an objective function that result in a maximum or minimum. Optimization methods are divided into exact and approximate (algorithms). Several optimization algorithms imitate natural phenomena, laws of physics, and behavior of living organisms. Optimization based on algorithms is the challenge that underlies machine learning, from logistic regression to training neural networks for artificial intelligence. In this paper, a new algorithm called two-stage optimization (TSO) is proposed. The TSO algorithm updates population members in two steps at each iteration. For this purpose, a group of good population members is selected and then two members of this group are randomly used to update the position of each of them. This update is based on the first selected good member at the first stage, and on the second selected good member at the second stage. We describe the stages of the TSO algorithm and model them mathematically. Performance of the TSO algorithm is evaluated for twenty-three standard objective functions. In order to compare the optimization results of the TSO algorithm, eight other competing algorithms are considered, including genetic, gravitational search, grey wolf, marine predators, particle swarm, teaching-learning-based, tunicate swarm, and whale approaches. The numerical results show that the new algorithm is superior and more competitive in solving optimization problems when compared with other algorithms. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart of the TSO algorithm.</p>
Full article ">Figure 2
<p>Plots of the objective function average with y-axis in logarithm scale for the indicated algorithm and function.</p>
Full article ">Figure 2 Cont.
<p>Plots of the objective function average with y-axis in logarithm scale for the indicated algorithm and function.</p>
Full article ">Figure 2 Cont.
<p>Plots of the objective function average with y-axis in logarithm scale for the indicated algorithm and function.</p>
Full article ">
34 pages, 1413 KiB  
Article
Understanding the Variability in Graph Data Sets through Statistical Modeling on the Stiefel Manifold
by Clément Mantoux, Baptiste Couvy-Duchesne, Federica Cacciamani, Stéphane Epelbaum, Stanley Durrleman and Stéphanie Allassonnière
Entropy 2021, 23(4), 490; https://doi.org/10.3390/e23040490 - 20 Apr 2021
Cited by 2 | Viewed by 3380
Abstract
Network analysis provides a rich framework to model complex phenomena, such as human brain connectivity. It has proven efficient to understand their natural properties and design predictive models. In this paper, we study the variability within groups of networks, i.e., the structure of [...] Read more.
Network analysis provides a rich framework to model complex phenomena, such as human brain connectivity. It has proven efficient to understand their natural properties and design predictive models. In this paper, we study the variability within groups of networks, i.e., the structure of connection similarities and differences across a set of networks. We propose a statistical framework to model these variations based on manifold-valued latent factors. Each network adjacency matrix is decomposed as a weighted sum of matrix patterns with rank one. Each pattern is described as a random perturbation of a dictionary element. As a hierarchical statistical model, it enables the analysis of heterogeneous populations of adjacency matrices using mixtures. Our framework can also be used to infer the weight of missing edges. We estimate the parameters of the model using an Expectation-Maximization-based algorithm. Experimenting on synthetic data, we show that the algorithm is able to accurately estimate the latent structure in both low and high dimensions. We apply our model on a large data set of functional brain connectivity matrices from the UK Biobank. Our results suggest that the proposed model accurately describes the complex variability in the data set with a small number of degrees of freedom. Full article
(This article belongs to the Special Issue Approximate Bayesian Inference)
Show Figures

Figure 1

Figure 1
<p>One thousand samples of three von Mises–Fisher distributions on <math display="inline"><semantics> <msub> <mi mathvariant="script">V</mi> <mrow> <mn>3</mn> <mo>,</mo> <mn>2</mn> </mrow> </msub> </semantics></math>. The mode of the distribution is represented by two red arrows along the <span class="html-italic">x</span> and <span class="html-italic">z</span> axes, and the two vectors in each matrix by two blue points. The concentration parameters are set to <math display="inline"><semantics> <mrow> <mrow> <mo>|</mo> </mrow> <msub> <mi>f</mi> <mi>z</mi> </msub> <mrow> <mo>|</mo> <mo>=</mo> <mn>10</mn> </mrow> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mrow> <mo>|</mo> </mrow> <msub> <mi>f</mi> <mi>x</mi> </msub> <mrow> <mo>|</mo> <mo> </mo> <mo>∈</mo> <mo> </mo> <mrow> <mo>[</mo> <mn>10</mn> <mo>,</mo> <mn>100</mn> <mo>,</mo> <mn>500</mn> <mo>]</mo> </mrow> </mrow> </mrow> </semantics></math> (from <b>left</b> to <b>right</b>). Samples are drawn with an adaptive Metropolis–Hastings sampler using the transition kernel described in <a href="#sec4-entropy-23-00490" class="html-sec">Section 4</a>. A stronger concentration of the <span class="html-italic">x</span> vector impacts the spread of the <span class="html-italic">z</span> vector.</p>
Full article ">Figure 2
<p>Graphical model for a data set of adjacency matrices <math display="inline"><semantics> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>…</mo> <mo>,</mo> <msub> <mi>A</mi> <mi>N</mi> </msub> </mrow> </semantics></math>. The variables <math display="inline"><semantics> <mi>π</mi> </semantics></math> and <math display="inline"><semantics> <msup> <mi>z</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> </semantics></math> can be added to get a mixture model.</p>
Full article ">Figure 3
<p>True latent variables <math display="inline"><semantics> <msup> <mi>X</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> </semantics></math> and their posterior MCMC mean estimation. The red arrows represent the true <math display="inline"><semantics> <mrow> <msub> <mi>π</mi> <mi>V</mi> </msub> <mrow> <mo>(</mo> <mi>F</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> parameter and its estimate <math display="inline"><semantics> <mrow> <msub> <mi>π</mi> <mi>V</mi> </msub> <mrow> <mo>(</mo> <mover accent="true"> <mi>F</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> </mrow> </semantics></math>. (<b>a</b>) The true mode and samples. (<b>b</b>) Mode and samples estimates when <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <mi>ε</mi> </msub> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math>. (<b>c</b>) Mode and samples estimates when <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <mi>ε</mi> </msub> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>. The columns are rearranged using Algorithm 2 to ease visualization. The latent variables are accurately estimated when the noise is small. A stronger noise causes the estimated latent variables to spread over the Stiefel manifold.</p>
Full article ">Figure 4
<p>Convergence of the concentration parameters <math display="inline"><semantics> <mrow> <mo>(</mo> <mo>|</mo> <msub> <mi>f</mi> <mn>1</mn> </msub> <mo>|</mo> <mo>,</mo> <mo>⋯</mo> <mo>,</mo> <mo>|</mo> <msub> <mi>f</mi> <mi>p</mi> </msub> <mo>|</mo> <mo>)</mo> </mrow> </semantics></math> (<b>left</b>) and the mean eigenvalues <math display="inline"><semantics> <mi>μ</mi> </semantics></math> (<b>right</b>) over the SAEM iterations. The red lines represent the values of the parameters along the iterations. The black dotted lines represent the true values, which are grouped in batches to ease visualization. The convergence is fastest for the large eigenvalues and concentration parameters. At the start of the algorithm, the biggest changes in the parameters come from the greedy permutation performed every 5 iterations. As explained in the text, the concentration parameters are underestimated. However, they keep the right order of magnitude, which allows interpreting the output of the algorithm in practice.</p>
Full article ">Figure 5
<p>Von Mises-Stiefel distribution parameter <span class="html-italic">F</span> and its estimation <math display="inline"><semantics> <mover accent="true"> <mi>F</mi> <mo>^</mo> </mover> </semantics></math>. (<b>Top row</b>): the two parameters and their difference. (<b>Bottom row</b>): mode of the true distribution (given by <math display="inline"><semantics> <mrow> <msub> <mi>π</mi> <mi>V</mi> </msub> <mrow> <mo>(</mo> <mi>F</mi> <mo>)</mo> </mrow> </mrow> </semantics></math>), mode of the estimated distribution <math display="inline"><semantics> <mrow> <msub> <mi>π</mi> <mi>V</mi> </msub> <mrow> <mo>(</mo> <mover accent="true"> <mi>F</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> </mrow> </semantics></math> and their difference. The images show each matrix as an array of coefficients, with pixel color corresponding to coefficient amplitude. Since the matrix columns are orthonormal, the projection just consists of normalizing the columns. The columns are sorted by decreasing the concentration parameter. The normalized columns of <span class="html-italic">F</span> corresponding to the smallest concentration parameters are estimated with less precision.</p>
Full article ">Figure 6
<p>Relative RMSE of parameters <span class="html-italic">F</span> and <math display="inline"><semantics> <mi>μ</mi> </semantics></math> after 100 MCMC-SAEM iterations depending on the number of MCMC steps per SAEM iteration. Results are averaged over 10 experiments to reduce the variance. The shaded areas indicate the extremal values across the repeated experiments. When using the greedy permutation, the rRMSE decreases rapidly when the number of MCMC steps increases before stabilizing. On the other hand, without the permutation step, the performance stays poor for any number of MCMC steps per maximization, as the parameters cannot be estimated correctly. In this experiment only, the latent variables are initialized at random to highlight the result.</p>
Full article ">Figure 7
<p>Result for missing link inference using the posterior distribution. (<b>a</b>) Ground truth input matrix <span class="html-italic">A</span>. (<b>b</b>) Posterior mean of the masked coefficients. (<b>c</b>) MAP estimator. (<b>d</b>) Mean of model samples for comparison. The area of masked edges is highlighted by a black square. Above each matrix is the rRMSE with the ground truth. Both the posterior mean and the MAP give a reasonable estimation for the missing weights, significantly better than the empirical mean of all adjacency matrices, which is the base reference for missing data imputation. The images show each matrix as an array of coefficients, with pixel color corresponding to coefficient amplitude.</p>
Full article ">Figure 8
<p>True latent variables <math display="inline"><semantics> <msup> <mi>X</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> </semantics></math> and their posterior mean estimation for the clustering problem. (<b>Top row</b>): the plots (<b>a</b>–<b>c</b>) represent the true vMF modes (in red), as well as the true <math display="inline"><semantics> <msup> <mi>X</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> </semantics></math> samples (in green) in their true class. (<b>Bottom row</b>): the plots (<b>d</b>–<b>f</b>) represent the three estimated vMF central modes (in red) and the estimated <math display="inline"><semantics> <msup> <mi>X</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </msup> </semantics></math> in their estimated class (in blue). The cluster centers are well recovered, as well as the concentration parameters. In particular, the two first clusters, which mainly differ by their concentration parameters, are correctly separated.</p>
Full article ">Figure 9
<p>Functional connectivity matrices (<math display="inline"><semantics> <mrow> <mn>21</mn> <mo>×</mo> <mn>21</mn> </mrow> </semantics></math> ) of 25 UK Biobank subjects. The connectivity structure changes a lot depending on the subject, with various patterns expressing with different weights. The matrices in the data set have no diagonal coefficients; hence, the diagonals are shown as zero.</p>
Full article ">Figure 10
<p>Normalized rank-one connectivity patterns. The matrix <span class="html-italic">i</span> represents <math display="inline"><semantics> <mrow> <mi>sign</mi> <mrow> <mo>(</mo> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <msubsup> <mi>f</mi> <mi>i</mi> <mo>⊤</mo> </msubsup> <mo>/</mo> <msup> <mfenced separators="" open="&#x2225;" close="&#x2225;"> <msub> <mi>f</mi> <mi>i</mi> </msub> </mfenced> <mn>2</mn> </msup> </mrow> </semantics></math>. The caption above each pattern gives the related concentration parameter and mean eigenvalue. The diagonal coefficients are set to zero, as they do not correspond to values in the data set. The images show each matrix as an array of coefficients, with pixel color corresponding to coefficient amplitude.</p>
Full article ">Figure 11
<p>(<b>a</b>) UK Biobank connectivity matrices for 5 subjects. (<b>b</b>) Corresponding posterior mean value of <math display="inline"><semantics> <mrow> <mi>λ</mi> <mo>·</mo> <mi>X</mi> </mrow> </semantics></math> estimated by the MCMC-SAEM. (<b>c</b>) Projection of the true connectivity matrices onto the subspace of the first five PCA components. The posterior mean matrix achieves a better rRMSE than PCA by capturing the main patterns of each individual matrix. As in <a href="#entropy-23-00490-f010" class="html-fig">Figure 10</a>, the diagonal cofficients are set to zero.</p>
Full article ">Figure 12
<p>From left to right: (<b>a</b>) True connectivity matrix <span class="html-italic">A</span>. (<b>b</b>) MAP estimator for the masked coefficients framed in a black square. (<b>c</b>) Linear model prediction for the masked coefficients. (<b>d</b>) Rank 5 truncation of the matrix <span class="html-italic">A</span> with masked coefficients set to zero. (<b>e</b>) Mean of all data set matrices. Above each matrix is the rRMSE with the ground truth.</p>
Full article ">Figure A1
<p>Normalized connectivity patterns when <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, computed as in <a href="#entropy-23-00490-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure A2
<p>Normalized connectivity patterns when <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>, computed as in <a href="#entropy-23-00490-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure A3
<p>(<b>a</b>) UK Biobank connectivity matrices for 5 subjects. (<b>b</b>) M10 posterior mean value of <math display="inline"><semantics> <mrow> <mi>λ</mi> <mo>·</mo> <mi>X</mi> </mrow> </semantics></math>. (<b>c</b>) M5 posterior mean value of <math display="inline"><semantics> <mrow> <mi>λ</mi> <mo>·</mo> <mi>X</mi> </mrow> </semantics></math>. (<b>d</b>) M2 posterior mean value of <math display="inline"><semantics> <mrow> <mi>λ</mi> <mo>·</mo> <mi>X</mi> </mrow> </semantics></math>. The rRMSE coherently increases as <span class="html-italic">p</span> decreases.</p>
Full article ">Figure A4
<p>Frontal, sagittal, and transverse cuts of the brain for the UK Biobank fMRI brain regions analyzed in this paper. As explained in <a href="#sec5dot2-entropy-23-00490" class="html-sec">Section 5.2</a>, region 1 comprises part of the Default Mode Network of the brain, which characterizes its activity at rest. Region 3, which is anti-correlated to region 1, is related to sensory functions. Regions 2, 4, 8, 9, and 19 are involved in the visual functions. Regions 10, 11, 12 correspond to motor control. Region 17 is involved in memory and spatial navigation. The L/R letters distinguish the left and right hemispheres. The black axes on each view give the three-dimensional position of the cut. The color strength corresponds to the truncated ICA weight.</p>
Full article ">
20 pages, 6542 KiB  
Article
Numerical Simulation Study on Flow Laws and Heat Transfer of Gas Hydrate in the Spiral Flow Pipeline with Long Twisted Band
by Yongchao Rao, Lijun Li, Shuli Wang, Shuhua Zhao and Shidong Zhou
Entropy 2021, 23(4), 489; https://doi.org/10.3390/e23040489 - 20 Apr 2021
Cited by 5 | Viewed by 2112
Abstract
The natural gas hydrate plugging problems in the mixed pipeline are becoming more and more serious. The hydrate plugging has gradually become an important problem to ensure the safety of pipeline operation. The deposition and heat transfer characteristics of natural gas hydrate particles [...] Read more.
The natural gas hydrate plugging problems in the mixed pipeline are becoming more and more serious. The hydrate plugging has gradually become an important problem to ensure the safety of pipeline operation. The deposition and heat transfer characteristics of natural gas hydrate particles in the spiral flow pipeline have been studied. The DPM model (discrete phase model) was used to simulate the motion of solid particles, which was used to simulate the complex spiral flow characteristics of hydrate in the pipeline with a long twisted band. The deposition and heat transfer characteristics of gas hydrate particles in the spiral flow pipeline were studied. The velocity distribution, pressure drop distribution, heat transfer characteristics, and particle settling characteristics in the pipeline were investigated. The numerical results showed that compared with the straight flow without a long twisted band, two obvious eddies are formed in the flow field with a long twisted band, and the velocities are maximum at the center of the vortices. Along the direction of the pipeline, the two vortices move toward the pipe wall from near the twisted band, which can effectively carry the hydrate particles deposited on the wall. With the same Reynolds number, the twisted rate was greater, the spiral strength was weaker, the tangential velocity was smaller, and the pressure drop was smaller. Therefore, the pressure loss can be reduced as much as possible with effect of the spiral flow. In a straight light flow, the Nusselt number is in a parabolic shape with the opening downwards. At the center of the pipe, the Nusselt number gradually decreased toward the pipe wall at the maximum, and at the near wall, the attenuation gradient of the Nu number was large. For spiral flow, the curve presented by the Nusselt number was a trough at the center of the pipe and a peak at 1/2 of the pipe diameter. With the reduction of twist rate, the Nusselt number becomes larger. Therefore, the spiral flow can make the temperature distribution more even and prevent the large temperature difference, resulting in the mass formation of hydrate particles in the pipeline wall. Spiral flow has a good carrying effect. Under the same condition, the spiral flow carried hydrate particles at a distance about 3–4 times farther than that of the straight flow. Full article
Show Figures

Figure 1

Figure 1
<p>The physical model.</p>
Full article ">Figure 2
<p>Twisted band model.</p>
Full article ">Figure 3
<p>Twisted rate schematic diagram of twisted band.</p>
Full article ">Figure 4
<p>Calculation grid for long twisted band pipe section.</p>
Full article ">Figure 5
<p>The velocity distribution on pipeline section at Z = 5D with different mesh sizes.</p>
Full article ">Figure 6
<p>The effect of Reynolds number on pressure drop.</p>
Full article ">Figure 7
<p>The effect of Reynolds number on Nusselt number.</p>
Full article ">Figure 8
<p>Velocity distribution clouds at different sections under Re = 15,000 and different twisted rates.</p>
Full article ">Figure 9
<p>Velocity and vector distribution nephogram in the long twisted band Y = 6.2.</p>
Full article ">Figure 10
<p>Relationship between Reynolds number and pressure drop.</p>
Full article ">Figure 11
<p>Relationship between twisted rate and pressure drop.</p>
Full article ">Figure 12
<p>Temperature nephogram of each section in long twisted band pipeline with different twisted rate at Re = 20,000.</p>
Full article ">Figure 12 Cont.
<p>Temperature nephogram of each section in long twisted band pipeline with different twisted rate at Re = 20,000.</p>
Full article ">Figure 13
<p>Variation curve of Nusselt number with Re in a long twisted band pipeline.</p>
Full article ">Figure 14
<p>Nusselt number distribution curves of different sections under different twisted rates at Re = 20,000.</p>
Full article ">Figure 15
<p>Nusselt number distribution curve on Z = 5D section at Re = 20,000.</p>
Full article ">Figure 16
<p>Distribution of particle concentration at different sections of the pipeline.</p>
Full article ">Figure 16 Cont.
<p>Distribution of particle concentration at different sections of the pipeline.</p>
Full article ">Figure 17
<p>Volume fraction distribution curve of particles along pipeline at Y = 6.2.</p>
Full article ">Figure 18
<p>Volume fraction distribution curve of particles along pipeline at Y = 6.2, α<sub>0</sub> = 8%.</p>
Full article ">Figure 19
<p>Volume fraction distribution curve of different twisted band particles along pipeline at Re = 4000, α<sub>0</sub> = 8%.</p>
Full article ">
15 pages, 946 KiB  
Article
Study of Dependence of Kinetic Freezeout Temperature on the Production Cross-Section of Particles in Various Centrality Intervals in Au–Au and Pb–Pb Collisions at High Energies
by Muhammad Waqas and Guang-Xiong Peng
Entropy 2021, 23(4), 488; https://doi.org/10.3390/e23040488 - 20 Apr 2021
Cited by 8 | Viewed by 2551
Abstract
Transverse momentum spectra of π+, p, Λ, Ξ or Ξ¯+, Ω or Ω¯+ and deuteron (d) in different centrality intervals in nucleus–nucleus collisions at the center of mass energy are analyzed by [...] Read more.
Transverse momentum spectra of π+, p, Λ, Ξ or Ξ¯+, Ω or Ω¯+ and deuteron (d) in different centrality intervals in nucleus–nucleus collisions at the center of mass energy are analyzed by the blast wave model with Boltzmann Gibbs statistics. We extracted the kinetic freezeout temperature, transverse flow velocity and kinetic freezeout volume from the transverse momentum spectra of the particles. It is observed that the non-strange and strange (multi-strange) particles freezeout separately due to different reaction cross-sections. While the freezeout volume and transverse flow velocity are mass dependent, they decrease with the resting mass of the particles. The present work reveals the scenario of a double kinetic freezeout in nucleus–nucleus collisions. Furthermore, the kinetic freezeout temperature and freezeout volume are larger in central collisions than peripheral collisions. However, the transverse flow velocity remains almost unchanged from central to peripheral collisions. Full article
Show Figures

Figure 1

Figure 1
<p>Transverse momentum spectra of <math display="inline"><semantics> <msup> <mi>π</mi> <mo>+</mo> </msup> </semantics></math>, <span class="html-italic">p</span>, <math display="inline"><semantics> <mo>Λ</mo> </semantics></math>, <math display="inline"><semantics> <msup> <mover accent="true"> <mo>Ξ</mo> <mo>¯</mo> </mover> <mo>+</mo> </msup> </semantics></math> and <math display="inline"><semantics> <msup> <mover accent="true"> <mo>Ω</mo> <mo>¯</mo> </mover> <mo>+</mo> </msup> </semantics></math> rapidity at <math display="inline"><semantics> <mrow> <mo>|</mo> <mi>y</mi> <mo>|</mo> <mo>&lt;</mo> <mn>0.1</mn> </mrow> </semantics></math>, and deuteron <math display="inline"><semantics> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> </semantics></math> at rapidity <math display="inline"><semantics> <mrow> <mo>|</mo> <mi>y</mi> <mo>|</mo> <mo>&lt;</mo> <mn>0.3</mn> </mrow> </semantics></math>, produced in different centrality intervals in Au–Au collisions at 62.4 GeV. Different symbols represent the <math display="inline"><semantics> <msub> <mi>p</mi> <mi>T</mi> </msub> </semantics></math> spectra of different particles measured by the STAR collaboration [<a href="#B21-entropy-23-00488" class="html-bibr">21</a>,<a href="#B55-entropy-23-00488" class="html-bibr">55</a>,<a href="#B56-entropy-23-00488" class="html-bibr">56</a>] and the curves are our fitted results with the blast wave model with Boltzmann Gibbs statistics (BGBW). The corresponding results of the data/fit are presented in each panel.</p>
Full article ">Figure 2
<p>Transverse momentum spectra of <math display="inline"><semantics> <msup> <mi>π</mi> <mo>+</mo> </msup> </semantics></math>, <span class="html-italic">p</span>, <math display="inline"><semantics> <mo>Λ</mo> </semantics></math>, <math display="inline"><semantics> <mo>Ξ</mo> </semantics></math>, <math display="inline"><semantics> <mo>Ω</mo> </semantics></math> and deuteron <math display="inline"><semantics> <mrow> <mo>(</mo> <mi>d</mi> <mo>)</mo> </mrow> </semantics></math> produced in different centrality intervals in Pb–Pb collisions at 2.76 TeV at rapidity <math display="inline"><semantics> <mrow> <mo>|</mo> <mi>y</mi> <mo>|</mo> <mo>&lt;</mo> <mn>0.5</mn> </mrow> </semantics></math>. Different symbols represent the <math display="inline"><semantics> <msub> <mi>p</mi> <mi>T</mi> </msub> </semantics></math> spectra of different particles measured by the ALICE collaboration [<a href="#B57-entropy-23-00488" class="html-bibr">57</a>,<a href="#B58-entropy-23-00488" class="html-bibr">58</a>,<a href="#B59-entropy-23-00488" class="html-bibr">59</a>] and the curves are our fitted results with the BGBW model. The corresponding results of the data/fit are presented in each panel.</p>
Full article ">Figure 3
<p>Dependence of <math display="inline"><semantics> <msub> <mi>T</mi> <mn>0</mn> </msub> </semantics></math> on the centrality class (<math display="inline"><semantics> <mrow> <mi>C</mi> <mo>%</mo> </mrow> </semantics></math>) and resting mass (<math display="inline"><semantics> <msub> <mi>m</mi> <mn>0</mn> </msub> </semantics></math>) of the particle.</p>
Full article ">Figure 4
<p>Dependence of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>T</mi> </msub> </semantics></math> on the centrality class (<math display="inline"><semantics> <mrow> <mi>C</mi> <mo>%</mo> </mrow> </semantics></math>) and resting mass (<math display="inline"><semantics> <msub> <mi>m</mi> <mn>0</mn> </msub> </semantics></math>) of the particle.</p>
Full article ">Figure 5
<p>Dependence of <span class="html-italic">V</span> on the centrality class (<math display="inline"><semantics> <mrow> <mi>C</mi> <mo>%</mo> </mrow> </semantics></math>) and resting mass (<math display="inline"><semantics> <msub> <mi>m</mi> <mn>0</mn> </msub> </semantics></math>) of the particle.</p>
Full article ">
18 pages, 9820 KiB  
Article
Research on Dynamic Evolution Model and Method of Communication Network Based on Real War Game
by Tongliang Lu, Kai Chen, Yan Zhang and Qiling Deng
Entropy 2021, 23(4), 487; https://doi.org/10.3390/e23040487 - 20 Apr 2021
Cited by 7 | Viewed by 2578
Abstract
Based on the data in real combat games, the combat System-of-Systems is usually composed of a large number of armed equipment platforms (or systems) and a reasonable communication network to connect mutually independent weapons and equipment platforms to achieve tasks such as information [...] Read more.
Based on the data in real combat games, the combat System-of-Systems is usually composed of a large number of armed equipment platforms (or systems) and a reasonable communication network to connect mutually independent weapons and equipment platforms to achieve tasks such as information collection, sharing, and collaborative processing. However, the generation algorithm of the combat system in the existing research is too simple and not suitable for reality. To overcome this problem, this paper proposes a communication network generation algorithm by adopting the joint distribution strategy of power law distribution and Poisson distribution to model the communication network. The simulation method is used to study the operation under continuous attack on communication nodes. The comprehensive experimental results of the dynamic evolution of the combat network in the battle scene verify the rationality and effectiveness of the communication network construction. Full article
Show Figures

Figure 1

Figure 1
<p>The diagram of relational hierarchy of the network structure of the combat SoS.</p>
Full article ">Figure 2
<p>The generation of the combat network.</p>
Full article ">Figure 3
<p>Hierarchical diagram of combat network topology.</p>
Full article ">Figure 4
<p>C2 subnetwork degree distribution. (<b>a</b>) The network topology. (<b>b</b>). Histogram of the node digrees of C2 subnetwork. (<b>c</b>) Log distribution of the node digrees of C2 subnetwork.</p>
Full article ">Figure 5
<p>Sensors subnetwork degree distribution. (<b>a</b>) The network topology. (<b>b</b>). Histogram of the node digrees of sensors subnetwork. (<b>c</b>) Log distribution of the node digrees of sensors subnetwork.</p>
Full article ">Figure 6
<p>Degree distribution of the communication subnetwork. (<b>a</b>) The network topology. (<b>b</b>). Histogram of the node digrees of communication subnetwork. (<b>c</b>) Log distribution of the node digrees of communication subnetwork.</p>
Full article ">Figure 7
<p>Attack subnetwork degree distribution. (<b>a</b>) The network topology. (<b>b</b>). Histogram of the node digrees of attack subnetwork. (<b>c</b>) Log distribution of the node digrees of attack subnetwork.</p>
Full article ">Figure 8
<p>Controlled experimental random network dynamic attack evolution diagram of combat system in continuous time series.</p>
Full article ">Figure 9
<p>Changes in the number of nodes and edges of differently distributed combat networks under <span class="html-italic">α</span> = 0.2, <span class="html-italic">β</span> = 0.5 attack.</p>
Full article ">Figure 10
<p>Changes in the robustness of differently distributed combat networks under <span class="html-italic">α</span> = 0.2, <span class="html-italic">β</span> = 0.5 attack.</p>
Full article ">Figure 10 Cont.
<p>Changes in the robustness of differently distributed combat networks under <span class="html-italic">α</span> = 0.2, <span class="html-italic">β</span> = 0.5 attack.</p>
Full article ">Figure 11
<p>Error analysis of robustness of differently distributed combat networks under <span class="html-italic">α</span> = 0.2, <span class="html-italic">β</span> = 0.5 attack.</p>
Full article ">Figure 11 Cont.
<p>Error analysis of robustness of differently distributed combat networks under <span class="html-italic">α</span> = 0.2, <span class="html-italic">β</span> = 0.5 attack.</p>
Full article ">Figure 12
<p>Changes in various indicators of the same distributed combat network under <span class="html-italic">β</span> = 0.5.</p>
Full article ">Figure 12 Cont.
<p>Changes in various indicators of the same distributed combat network under <span class="html-italic">β</span> = 0.5.</p>
Full article ">Figure 13
<p>Changes in the indicators of the same distributed combat network under <span class="html-italic">α</span> = 0.2.</p>
Full article ">Figure 13 Cont.
<p>Changes in the indicators of the same distributed combat network under <span class="html-italic">α</span> = 0.2.</p>
Full article ">Figure A1
<p>The relationship between the sea and air joint strike command organization.</p>
Full article ">
18 pages, 9188 KiB  
Article
Porcelain Insulator Crack Location and Surface States Pattern Recognition Based on Hyperspectral Technology
by Yiming Zhao, Jing Yan, Yanxin Wang, Qianzhen Jing and Tingliang Liu
Entropy 2021, 23(4), 486; https://doi.org/10.3390/e23040486 - 20 Apr 2021
Cited by 8 | Viewed by 3327
Abstract
A porcelain insulator is an important part to ensure that the insulation requirements of power equipment can be met. Under the influence of their structure, porcelain insulators are prone to mechanical damage and cracks, which will reduce their insulation performance. After a long-term [...] Read more.
A porcelain insulator is an important part to ensure that the insulation requirements of power equipment can be met. Under the influence of their structure, porcelain insulators are prone to mechanical damage and cracks, which will reduce their insulation performance. After a long-term operation, crack expansion will eventually lead to breakdown and safety hazards. Therefore, it is of great significance to detect insulator cracks to ensure the safe and reliable operation of a power grid. However, most traditional methods of insulator crack detection involve offline detection or contact measurement, which is not conducive to the online monitoring of equipment. Hyperspectral imaging technology is a noncontact detection technology containing three-dimensional (3D) spatial spectral information, whereby the data provide more information and the measuring method has a higher safety than electric detection methods. Therefore, a model of positioning and state classification of porcelain insulators based on hyperspectral technology is proposed. In this model, image data were used to extract edges to locate cracks, and spectral information was used to classify the surface states of porcelain insulators with EfficientNet. Lastly, crack extraction was realized, and the recognition accuracy of cracks and normal states was 96.9%. Through an analysis of the results, it is proven that the crack detection method of a porcelain insulator based on hyperspectral technology is an effective non-contact online monitoring approach, which has broad application prospects in the era of the Internet of Things with the rapid development of electric power. Full article
(This article belongs to the Special Issue Reliability of Modern Electro-Mechanical Systems)
Show Figures

Figure 1

Figure 1
<p>Crack location diagram: (<b>a</b>) structure of hyperspectral camera; (<b>b</b>) schematic diagram of hyperspectral imaging principle.</p>
Full article ">Figure 2
<p>Comparison of filtering methods: (<b>a</b>) image after adding Gaussian noise; (<b>b</b>) image after mean filtering; (<b>c</b>) image after box filtering; (<b>d</b>) image after Gaussian filtering.</p>
Full article ">Figure 3
<p>Crack location diagram: (<b>a</b>) image after Gaussian filtering; (<b>b</b>) gradient amplitude map; (<b>c</b>) image after Non-Maximum Suppression(NMS) processing; (<b>d</b>) image after double threshold detection and edge connection.</p>
Full article ">Figure 4
<p>The principle of crack detection based on the hyperspectral graph.</p>
Full article ">Figure 5
<p>EfficientNet-B0 network architecture.</p>
Full article ">Figure 6
<p>The overall framework of the insulator crack positioning and identification model.</p>
Full article ">Figure 7
<p>Hyperspectral detection platform: (<b>a</b>) diagram of the hyperspectral detection platform; (<b>b</b>) hyperspectral detection experimental platform.</p>
Full article ">Figure 8
<p>Hyperspectral three-dimensional (3D) data map of insulator.</p>
Full article ">Figure 9
<p>Black-and-white correction and multiplicative scatter correction (MSC): (<b>a</b>) normal sample spectrum after black-and-white correction; (<b>b</b>) normal sample spectrum after MSC; (<b>c</b>) spectrum of crack samples after black-and-white correction; (<b>d</b>) crack sample spectrum after correction.</p>
Full article ">Figure 10
<p>Schematic of pure pixel exponential end-member extraction algorithm.</p>
Full article ">Figure 11
<p>Spectral dimension reduction results.</p>
Full article ">Figure 12
<p>Crack localization of porcelain insulator sample 1: (<b>a</b>) image after Gaussian filtering; (<b>b</b>) gradient amplitude map; (<b>c</b>) image after NMS processing; (<b>d</b>) image after double threshold detection and edge connection.</p>
Full article ">Figure 13
<p>Crack localization of porcelain insulator sample 2: (<b>a</b>) image after Gaussian filtering; (<b>b</b>) gradient amplitude map; (<b>c</b>) image after NMS processing; (<b>d</b>) image after double threshold detection and edge connection.</p>
Full article ">Figure 14
<p>Spectral sample collection: (<b>a</b>) crack region spectrum of sample 1; (<b>b</b>) crack region spectrum of sample 2; (<b>c</b>) normal region spectrum of sample 3; (<b>d</b>) normal region spectrum of sample 4.</p>
Full article ">Figure 15
<p>The schematic diagram of 5-fold cross-validation.</p>
Full article ">Figure 16
<p>Porcelain insulator state identification accuracy results. SVM, support vector machine; DT, decision tree; RF, random forest; BPNN, back propagation neural networks.</p>
Full article ">
23 pages, 602 KiB  
Article
Knowledge Discovery for Higher Education Student Retention Based on Data Mining: Machine Learning Algorithms and Case Study in Chile
by Carlos A. Palacios, José A. Reyes-Suárez, Lorena A. Bearzotti, Víctor Leiva and Carolina Marchant
Entropy 2021, 23(4), 485; https://doi.org/10.3390/e23040485 - 20 Apr 2021
Cited by 73 | Viewed by 8789
Abstract
Data mining is employed to extract useful information and to detect patterns from often large data sets, closely related to knowledge discovery in databases and data science. In this investigation, we formulate models based on machine learning algorithms to extract relevant information predicting [...] Read more.
Data mining is employed to extract useful information and to detect patterns from often large data sets, closely related to knowledge discovery in databases and data science. In this investigation, we formulate models based on machine learning algorithms to extract relevant information predicting student retention at various levels, using higher education data and specifying the relevant variables involved in the modeling. Then, we utilize this information to help the process of knowledge discovery. We predict student retention at each of three levels during their first, second, and third years of study, obtaining models with an accuracy that exceeds 80% in all scenarios. These models allow us to adequately predict the level when dropout occurs. Among the machine learning algorithms used in this work are: decision trees, k-nearest neighbors, logistic regression, naive Bayes, random forest, and support vector machines, of which the random forest technique performs the best. We detect that secondary educational score and the community poverty index are important predictive variables, which have not been previously reported in educational studies of this type. The dropout assessment at various levels reported here is valid for higher education institutions around the world with similar conditions to the Chilean case, where dropout rates affect the efficiency of such institutions. Having the ability to predict dropout based on student’s data enables these institutions to take preventative measures, avoiding the dropouts. In the case study, balancing the majority and minority classes improves the performance of the algorithms. Full article
Show Figures

Figure 1

Figure 1
<p>Scheme of the four different models proposed to predict student retention/dropout, where PSU indicates the university selection test.</p>
Full article ">Figure 2
<p>Scheme of the KDD methodology.</p>
Full article ">Figure 3
<p>Confusion matrix and performance metrics.</p>
Full article ">Figure 4
<p>Dashboard of enrollment in the UCM.</p>
Full article ">
33 pages, 12047 KiB  
Article
A Volatility Estimator of Stock Market Indices Based on the Intrinsic Entropy Model
by Claudiu Vințe, Marcel Ausloos and Titus Felix Furtună
Entropy 2021, 23(4), 484; https://doi.org/10.3390/e23040484 - 19 Apr 2021
Cited by 8 | Viewed by 5848
Abstract
Grasping the historical volatility of stock market indices and accurately estimating are two of the major focuses of those involved in the financial securities industry and derivative instruments pricing. This paper presents the results of employing the intrinsic entropy model as a substitute [...] Read more.
Grasping the historical volatility of stock market indices and accurately estimating are two of the major focuses of those involved in the financial securities industry and derivative instruments pricing. This paper presents the results of employing the intrinsic entropy model as a substitute for estimating the volatility of stock market indices. Diverging from the widely used volatility models that take into account only the elements related to the traded prices, namely the open, high, low, and close prices of a trading day (OHLC), the intrinsic entropy model takes into account the traded volumes during the considered time frame as well. We adjust the intraday intrinsic entropy model that we introduced earlier for exchange-traded securities in order to connect daily OHLC prices with the ratio of the corresponding daily volume to the overall volume traded in the considered period. The intrinsic entropy model conceptualizes this ratio as entropic probability or market credence assigned to the corresponding price level. The intrinsic entropy is computed using historical daily data for traded market indices (S&P 500, Dow 30, NYSE Composite, NASDAQ Composite, Nikkei 225, and Hang Seng Index). We compare the results produced by the intrinsic entropy model with the volatility estimates obtained for the same data sets using widely employed industry volatility estimators. The intrinsic entropy model proves to consistently deliver reliable estimates for various time frames while showing peculiarly high values for the coefficient of variation, with the estimates falling in a significantly lower interval range compared with those provided by the other advanced volatility estimators. Full article
(This article belongs to the Special Issue Information Theory and Economic Network)
Show Figures

Figure 1

Figure 1
<p>S&amp;P 500 volatility for a 20 day time interval produced by the intrinsic entropy-based estimator.</p>
Full article ">Figure 2
<p>S&amp;P 500 volatility for a 20 day time interval produced by the Yang–Zhang estimator.</p>
Full article ">Figure 3
<p>S&amp;P 500 volatility for a 60 day time interval produced by the intrinsic entropy-based estimator.</p>
Full article ">Figure 4
<p>S&amp;P 500 volatility for a 60 day time interval produced by the Yang–Zhang estimator.</p>
Full article ">Figure 5
<p>Mean, Var, and CV of S&amp;P 500 volatility estimates for 5 day time windows.</p>
Full article ">Figure 6
<p>Mean, Var, and CV of S&amp;P 500 volatility estimates for 10 day time windows.</p>
Full article ">Figure 7
<p>Mean, Var, and CV of S&amp;P 500 volatility estimates for 15 day time windows.</p>
Full article ">Figure 8
<p>Mean, Var, and CV of S&amp;P 500 volatility estimates for 20 day time windows.</p>
Full article ">Figure 9
<p>Mean, Var, and CV of S&amp;P 500 volatility estimates for 30 day time windows.</p>
Full article ">Figure 10
<p>Mean, Var, and CV of S&amp;P 500 volatility estimates for 60 day time windows.</p>
Full article ">Figure 11
<p>The volatility estimators’ efficiency for the S&amp;P 500 index over a 30 day time window.</p>
Full article ">Figure 12
<p>The volatility estimators’ efficiency for the NYSE Composite index over a 30 day time window.</p>
Full article ">Figure 13
<p>Yang–Zhang and intrinsic entropy-based estimates for the S&amp;P 500 stock market index over a 260 day time window.</p>
Full article ">Figure 14
<p>Yang–Zhang and intrinsic entropy-based estimates for the NYSE Composite stock market index over a 260 day time window.</p>
Full article ">Figure A1
<p>Volatility estimators’ efficiency for the Dow 30 index over a 30 day time window.</p>
Full article ">Figure A2
<p>Volatility estimators’ efficiency for the NASDAQ Composite index over a 30 day time window.</p>
Full article ">Figure A3
<p>Volatility estimators’ efficiency for the Russell 2000 index over a 30 day time window.</p>
Full article ">Figure A4
<p>Volatility estimators’ efficiency for the Nikkei 225 index over a 30 day time window.</p>
Full article ">Figure A5
<p>Volatility estimators’ efficiency for the Hang Seng Index over a 30 day time window.</p>
Full article ">Figure A6
<p>Yang–Zhang estimates and the intrinsic entropy-based estimates for the Dow 30 stock market index over a 260 day time window.</p>
Full article ">Figure A7
<p>Yang–Zhang estimates and the intrinsic entropy-based estimates for the Russel 2000 stock market index over a 260 day time window.</p>
Full article ">Figure A8
<p>Yang–Zhang estimates and the intrinsic entropy-based estimates for the NASDAQ Composite stock market index over a 260 day time window.</p>
Full article ">Figure A9
<p>Yang–Zhang estimates and the intrinsic entropy-based estimates for the Nikkei 225 stock market index over a 260 day time window.</p>
Full article ">Figure A10
<p>Yang–Zhang estimates and the intrinsic entropy-based estimates for the Hang Seng stock market index over a 260 day time window.</p>
Full article ">
22 pages, 985 KiB  
Article
Resultant Information Descriptors, Equilibrium States and Ensemble Entropy
by Roman F. Nalewajski
Entropy 2021, 23(4), 483; https://doi.org/10.3390/e23040483 - 19 Apr 2021
Cited by 5 | Viewed by 2020
Abstract
In this article, sources of information in electronic states are reexamined and a need for the resultant measures of the entropy/information content, combining contributions due to probability and phase/current densities, is emphasized. Probability distribution reflects the wavefunction modulus and generates classical contributions to [...] Read more.
In this article, sources of information in electronic states are reexamined and a need for the resultant measures of the entropy/information content, combining contributions due to probability and phase/current densities, is emphasized. Probability distribution reflects the wavefunction modulus and generates classical contributions to Shannon’s global entropy and Fisher’s gradient information. The phase component of molecular states similarly determines their nonclassical supplements, due to probability “convection”. The local-energy concept is used to examine the phase equalization in the equilibrium, phase-transformed states. Continuity relations for the wavefunction modulus and phase components are reexamined, the convectional character of the local source of the resultant gradient information is stressed, and latent probability currents in the equilibrium (stationary) quantum states are related to the horizontal (“thermodynamic”) phase. The equivalence of the energy and resultant gradient information (kinetic energy) descriptors of chemical processes is stressed. In the grand-ensemble description, the reactivity criteria are defined by the populational derivatives of the system average electronic energy. Their entropic analogs, given by the associated derivatives of the overall gradient information, are shown to provide an equivalent set of reactivity indices for describing the charge transfer phenomena. Full article
(This article belongs to the Special Issue Entropic and Complexity Measures in Atomic and Molecular Systems)
Show Figures

Figure 1

Figure 1
<p>Classical (probability) and quantum (wavefunction) information schemes in molecular QM. The quantum mapping {<b><span class="html-italic">r</span></b> → <span class="html-italic">ψ</span>(<b><span class="html-italic">r</span></b>)} implies both the classical {<b><span class="html-italic">r</span></b>→<span class="html-italic">p</span>(<b><span class="html-italic">r</span></b>)} and nonclassical attributions {<b><span class="html-italic">r</span></b>→[<span class="html-italic">φ</span>(<b><span class="html-italic">r</span></b>), <b><span class="html-italic">j</span></b>(<b><span class="html-italic">r</span></b>) or <b><span class="html-italic">V</span></b>(<b><span class="html-italic">r</span></b>)]}.</p>
Full article ">Figure 2
<p>Local “vertical” (<span class="html-italic">v</span>) and “horizontal” (<span class="html-italic">h</span>) directions.</p>
Full article ">Figure 3
<p>Schematic diagrams of atomic and molecular vortices of “horizontal” flows of electronic probability density in atomic fragments of diatomic promolecule M<sup>0</sup>, the polarized system M<sup>+</sup>, and in molecule M.</p>
Full article ">
3 pages, 184 KiB  
Editorial
Information Theory in Molecular Evolution: From Models to Structures and Dynamics
by Faruck Morcos
Entropy 2021, 23(4), 482; https://doi.org/10.3390/e23040482 - 19 Apr 2021
Cited by 1 | Viewed by 2132
Abstract
Historically, information theory has been closely interconnected with evolutionary theory [...] Full article
18 pages, 1965 KiB  
Article
Performance of Portfolios Based on the Expected Utility-Entropy Fund Rating Approach
by Daniel Chiew, Judy Qiu, Sirimon Treepongkaruna, Jiping Yang and Chenxiao Shi
Entropy 2021, 23(4), 481; https://doi.org/10.3390/e23040481 - 18 Apr 2021
Viewed by 2222
Abstract
Yang and Qiu proposed and reframed an expected utility–entropy (EU-E) based decision model. Later on, a similar numerical representation for a risky choice was axiomatically developed by Luce et al. under the condition of segregation. Recently, we established a fund rating approach based [...] Read more.
Yang and Qiu proposed and reframed an expected utility–entropy (EU-E) based decision model. Later on, a similar numerical representation for a risky choice was axiomatically developed by Luce et al. under the condition of segregation. Recently, we established a fund rating approach based on the EU-E decision model and Morningstar ratings. In this paper, we apply the approach to US mutual funds and construct portfolios using the best rating funds. Furthermore, we evaluate the performance of the fund ratings based on the EU-E decision model against Morningstar ratings by examining the performance of the three models in portfolio selection. The conclusions show that portfolios constructed using the ratings based on the EU-E models with moderate tradeoff coefficients perform better than those constructed using Morningstar. The conclusion is robust to different rebalancing intervals. Full article
(This article belongs to the Special Issue Entropy Method for Decision Making)
Show Figures

Figure 1

Figure 1
<p>S&amp;P 500 index for 2002-08–2015-07. The date is shown on the horizontal axis and the level of the S&amp;P 500 index, in terms of basis points, on the vertical axis. GFC and EDC denote the global financial crisis and the European debt crisis periods, respectively.</p>
Full article ">Figure 2
<p>Number of rating upgrades by year. Note: This figure illustrates the number of rating upgrades from lower to higher grades (e.g., from 1 to 2 stars, 2 to 3 stars, etc.) using the ratings based on Morningstar and EU-E (λ = 0), EU-E (λ = 0.25), EU-E (λ = 0.50), EU-E (λ = 0.75), and EU-E (λ = 1).</p>
Full article ">Figure 3
<p>Summary of portfolio rebalancing periods. Each portfolio is rebalanced, and the performance of each portfolio is recorded at each interval, as illustrated by the interval number. <span class="html-italic">N</span> refers to the total number of portfolios constructed using each of the ratings based on Morningstar and the EU-E model, where λ takes a value of 0.25 and 0.75.</p>
Full article ">Figure 4
<p>Average abnormal returns of portfolios through time. The abnormal return of the rating measure relative to the benchmark as a % is illustrated on the vertical axis and the date on the horizontal axis.</p>
Full article ">Figure 5
<p>Summary of statistically significant abnormal returns over 12-month rebalancing intervals. This figure shows the 12-month interval on the horizontal axis and the number of statistically significant outperformances as a percentage on the vertical axis. Statistical significance is tested at the 10% level.</p>
Full article ">Figure 6
<p>Cumulative average abnormal returns during the sample period. The vertical axis shows the cumulative average abnormal return of the portfolios constructed using the rating models relative to the benchmark as a percentage.</p>
Full article ">
13 pages, 3015 KiB  
Article
Novel Features for Binary Time Series Based on Branch Length Similarity Entropy
by Sang-Hee Lee and Cheol-Min Park
Entropy 2021, 23(4), 480; https://doi.org/10.3390/e23040480 - 18 Apr 2021
Cited by 2 | Viewed by 2782
Abstract
Branch length similarity (BLS) entropy is defined in a network consisting of a single node and branches. In this study, we mapped the binary time-series signal to the circumference of the time circle so that the BLS entropy can be calculated for the [...] Read more.
Branch length similarity (BLS) entropy is defined in a network consisting of a single node and branches. In this study, we mapped the binary time-series signal to the circumference of the time circle so that the BLS entropy can be calculated for the binary time-series. We obtained the BLS entropy values for “1” signals on the time circle. The set of values are the BLS entropy profile. We selected the local maximum (minimum) point, slope, and inflection point of the entropy profile as the characteristic features of the binary time-series and investigated and explored their significance. The local maximum (minimum) point indicates the time at which the rate of change in the signal density becomes zero. The slope and inflection points correspond to the degree of change in the signal density and the time at which the signal density changes occur, respectively. Moreover, we show that the characteristic features can be widely used in binary time-series analysis by characterizing the movement trajectory of Caenorhabditis elegans. We also mention the problems that need to be explored mathematically in relation to the features and propose candidates for additional features based on the BLS entropy profile. Full article
Show Figures

Figure 1

Figure 1
<p>Definition of BLS entropy in a simple network composed of one node and several branches, and an example of BLS entropy change according to the change in the branch length.</p>
Full article ">Figure 2
<p>Time circle generated by mapping signals of random binary time-series with length <span class="html-italic">L</span> = 400 on the circumference of a circle and BLS entropy profile obtained from the “1” signals on the time circle.</p>
Full article ">Figure 3
<p>(<b>A</b>) Heterogeneous landscape according to the control variable, <span class="html-italic">H</span>, and its binary images through the binarization process, (<b>B</b>) binary time-series for the image, and (<b>C</b>) BLS entropy profiles corresponding to the time series. (<b>D</b>) Signal pattern with the enlarged blue square (upper figure) and red square (lower figure) on the BLS entropy profile.</p>
Full article ">Figure 4
<p>(<b>A</b>) Two binary time-series, <span class="html-italic">Q</span><sub>12</sub>(<span class="html-italic">t</span>) and <span class="html-italic">Q</span><sub>123</sub>(<span class="html-italic">t</span>), created by combining <span class="html-italic">Q</span><sub>1</sub>(<span class="html-italic">t</span>), <span class="html-italic">Q</span><sub>2</sub>(<span class="html-italic">t</span>), and <span class="html-italic">Q</span><sub>3</sub>(<span class="html-italic">t</span>). <span class="html-italic">Q</span><sub>12</sub>(<span class="html-italic">t</span>) was created by the sequential combination of <span class="html-italic">Q</span><sub>1</sub>(<span class="html-italic">t</span>) and <span class="html-italic">Q</span><sub>2</sub>(<span class="html-italic">t</span>), and <span class="html-italic">Q</span><sub>123</sub>(<span class="html-italic">t</span>) was generated by the combination of <span class="html-italic">Q</span><sub>1</sub>(<span class="html-italic">t</span>), <span class="html-italic">Q</span><sub>2</sub>(<span class="html-italic">t</span>), and <span class="html-italic">Q</span><sub>3</sub>(<span class="html-italic">t</span>). (<b>B</b>) BLS entropy profiles corresponding to <span class="html-italic">Q</span><sub>12</sub>(<span class="html-italic">t</span>) and <span class="html-italic">Q</span><sub>123</sub>(<span class="html-italic">t</span>), and (<b>C</b>) binary time-series, <span class="html-italic">Q</span>(<span class="html-italic">t</span>), in which the signal density linearly decreases and then increases, its BLS entropy profile, and the derivative function of the entropy profile (<math display="inline"><semantics> <mrow> <mi>S</mi> <mrow> <mo>(</mo> <mover accent="true"> <mi>t</mi> <mo>˜</mo> </mover> <mo>)</mo> </mrow> </mrow> </semantics></math>).</p>
Full article ">Figure 5
<p>Heterogeneous binary time-series with <span class="html-italic">H</span> = 0.3, <span class="html-italic">Q</span>(<span class="html-italic">t</span>), and its local maximum (minimum) points (marked by a black triangle) and inflection points (marked by a red triangle) on the BLS entropy profile for the time series.</p>
Full article ">Figure 6
<p>Crawling trajectories of <span class="html-italic">C. elegans</span> for the control, benzene-treated (0.5 ppm), and formaldehyde-treated (0.5 ppm) groups, and BLS entropy profiles for the trajectories. The triangles represent the locations of the local maximum points on the entropy profile.</p>
Full article ">Figure 7
<p>Binary time-series consisting of heterogeneously distributed “1” signals, <span class="html-italic">Q</span><sub>1</sub>(<span class="html-italic">t</span>) (<span class="html-italic">H</span> = 0.1), <span class="html-italic">Q</span><sub>2</sub>(<span class="html-italic">t</span>) (<span class="html-italic">H</span> = 0.2), <span class="html-italic">Q</span><sub>3</sub>(<span class="html-italic">t</span>) (<span class="html-italic">H</span> = 0.3), and their BLS entropy profiles. Here, <span class="html-italic">H</span> indicates the heterogeneity of the binary time-series. ED and <span class="html-italic">ρ</span> represent the similarity based on the ED and the similarity based on the BLS entropy profile, respectively.</p>
Full article ">Figure A1
<p>(<b>A</b>) Illustration showing the relationship between the BLS entropy value and the signal density for a network consisting of two line segments, <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <mi>p</mi> <mi>q</mi> </mrow> <mo stretchy="true">¯</mo> </mover> </mrow> </semantics></math> (of length <span class="html-italic">L</span><sub>1</sub>) and <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mo stretchy="true">¯</mo> </mover> </mrow> </semantics></math> (of length <span class="html-italic">L</span><sub>2</sub>) on the time circle, (<b>B</b>) Change in the BLS entropy value for the change in the <span class="html-italic">L</span><sub>2</sub> value, with a fixed <span class="html-italic">L</span><sub>1</sub>. The inlet graph shows the BLS entropy values when the point <span class="html-italic">p</span> has <span class="html-italic">n</span> line segments, and the length (<span class="html-italic">x</span>) of one of them changes.</p>
Full article ">
20 pages, 879 KiB  
Article
Entanglement and Non-Locality in Quantum Protocols with Identical Particles
by Fabio Benatti, Roberto Floreanini and Ugo Marzolino
Entropy 2021, 23(4), 479; https://doi.org/10.3390/e23040479 - 18 Apr 2021
Cited by 5 | Viewed by 3318
Abstract
We study the role of entanglement and non-locality in quantum protocols that make use of systems of identical particles. Unlike in the case of distinguishable particles, the notions of entanglement and non-locality for systems whose constituents cannot be distinguished and singly addressed are [...] Read more.
We study the role of entanglement and non-locality in quantum protocols that make use of systems of identical particles. Unlike in the case of distinguishable particles, the notions of entanglement and non-locality for systems whose constituents cannot be distinguished and singly addressed are still debated. We clarify why the only approach that avoids incongruities and paradoxes is the one based on the second quantization formalism, whereby it is the entanglement of the modes that can be populated by the particles that really matters and not the particles themselves. Indeed, by means of a metrological and of a teleportation protocol, we show that inconsistencies arise in formulations that force entanglement and non-locality to be properties of the identical particles rather than of the modes they can occupy. The reason resides in the fact that orthogonal modes can always be addressed while identical particles cannot. Full article
(This article belongs to the Special Issue Quantum Information and Quantum Optics)
23 pages, 650 KiB  
Article
Excitation Functions of Tsallis-Like Parameters in High-Energy Nucleus–Nucleus Collisions
by Li-Li Li, Fu-Hu Liu and Khusniddin K. Olimov
Entropy 2021, 23(4), 478; https://doi.org/10.3390/e23040478 - 18 Apr 2021
Cited by 21 | Viewed by 2663
Abstract
The transverse momentum spectra of charged pions, kaons, and protons produced at mid-rapidity in central nucleus–nucleus (AA) collisions at high energies are analyzed by considering particles to be created from two participant partons, which are assumed to be contributors from the collision system. [...] Read more.
The transverse momentum spectra of charged pions, kaons, and protons produced at mid-rapidity in central nucleus–nucleus (AA) collisions at high energies are analyzed by considering particles to be created from two participant partons, which are assumed to be contributors from the collision system. Each participant (contributor) parton is assumed to contribute to the transverse momentum by a Tsallis-like function. The contributions of the two participant partons are regarded as the two components of transverse momentum of the identified particle. The experimental data measured in high-energy AA collisions by international collaborations are studied. The excitation functions of kinetic freeze-out temperature and transverse flow velocity are extracted. The two parameters increase quickly from ≈3 to ≈10 GeV (exactly from 2.7 to 7.7 GeV) and then slowly at above 10 GeV with the increase of collision energy. In particular, there is a plateau from near 10 GeV to 200 GeV in the excitation function of kinetic freeze-out temperature. Full article
Show Figures

Figure 1

Figure 1
<p>Transverse mass spectra of charged pions, kaons, and protons produced in 0–5% Au-Au collisions at <math display="inline"><semantics> <mrow> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> <mo>=</mo> </mrow> </semantics></math> (<b>a</b>) 2.7, (<b>b</b>) 3.32, (<b>c</b>) 3.84, (<b>d</b>) 4.3, and (<b>e</b>) 5.03 GeV, and in 0–5% Pb-Pb collisions at <math display="inline"><semantics> <mrow> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> <mo>=</mo> </mrow> </semantics></math> (<b>f</b>) 6.3 GeV. In panel (<b>f</b>), the factor <math display="inline"><semantics> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> <mi>π</mi> </mrow> </semantics></math> does not appear, which causes different normalization from other panels. The symbols represent the experimental data at mid-<span class="html-italic">y</span> measured by the E866, E895, and E802 Collaboration at the AGS [<a href="#B31-entropy-23-00478" class="html-bibr">31</a>,<a href="#B32-entropy-23-00478" class="html-bibr">32</a>,<a href="#B33-entropy-23-00478" class="html-bibr">33</a>,<a href="#B34-entropy-23-00478" class="html-bibr">34</a>,<a href="#B35-entropy-23-00478" class="html-bibr">35</a>] and by the NA49 Collaboration at the SPS [<a href="#B36-entropy-23-00478" class="html-bibr">36</a>,<a href="#B37-entropy-23-00478" class="html-bibr">37</a>]. The solid and dashed curves are our results, fitted by using Equation (<a href="#FD10-entropy-23-00478" class="html-disp-formula">10</a>) due to Equations (7) and (9), with <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>μ</mi> <mi>B</mi> </msub> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math>, respectively.</p>
Full article ">Figure 2
<p>Transverse momentum spectra of charged pions, kaons, and protons produced in 0–5% Au-Au collisions at <math display="inline"><semantics> <mrow> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> <mo>=</mo> </mrow> </semantics></math> (<b>a</b>) 7.7, (<b>b</b>) 11.5, (<b>c</b>) 14.5, (<b>d</b>) 19.6, (<b>e</b>) 27, and (<b>f</b>) 39 GeV. In panel (<b>c</b>), the factor <math display="inline"><semantics> <mrow> <mn>1</mn> <mo>/</mo> <msub> <mi>N</mi> <mrow> <mi>E</mi> <mi>V</mi> </mrow> </msub> </mrow> </semantics></math>, i.e., the number of events is included on the vertical axis, which can be omitted. The symbols represent the experimental data at mid-<span class="html-italic">y</span> measured by the STAR Collaboration at the RHIC [<a href="#B38-entropy-23-00478" class="html-bibr">38</a>,<a href="#B39-entropy-23-00478" class="html-bibr">39</a>,<a href="#B40-entropy-23-00478" class="html-bibr">40</a>]. The solid and dashed curves are our results, fitted by using Equation (<a href="#FD10-entropy-23-00478" class="html-disp-formula">10</a>) due to Equations (7) and (9), with <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>μ</mi> <mi>B</mi> </msub> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math>, respectively.</p>
Full article ">Figure 3
<p>Transverse momentum spectra of charged pions, kaons, and protons produced in 0–5% Au-Au collisions at <math display="inline"><semantics> <mrow> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> <mo>=</mo> </mrow> </semantics></math> (<b>a</b>) 62.4, (<b>b</b>) 130, and (<b>c</b>) 200 GeV; in 0–5% Pb-Pb collisions at <math display="inline"><semantics> <mrow> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> <mo>=</mo> </mrow> </semantics></math> (<b>d</b>) 2.76 and (<b>e</b>) 5.02 TeV; and in 0–5% Xe-Xe collisions at <math display="inline"><semantics> <mrow> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> <mo>=</mo> </mrow> </semantics></math> (<b>f</b>) 5.44 TeV. In panels (<b>c</b>,<b>d</b>,<b>f</b>), the factor <math display="inline"><semantics> <mrow> <mn>1</mn> <mo>/</mo> <msub> <mi>N</mi> <mrow> <mi>E</mi> <mi>V</mi> </mrow> </msub> </mrow> </semantics></math> is included on the vertical axis, which can be omitted. In panels (<b>e</b>,<b>f</b>), the item <math display="inline"><semantics> <msup> <mrow> <mo>(</mo> <mn>2</mn> <mi>π</mi> <msub> <mi>p</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </semantics></math> is not included on the vertical axis, which results in different calculation for vertical values from other panels in the normalization. The symbols represent the experimental data at mid-<span class="html-italic">y</span> measured by the STAR Collaboration at the RHIC [<a href="#B38-entropy-23-00478" class="html-bibr">38</a>,<a href="#B39-entropy-23-00478" class="html-bibr">39</a>,<a href="#B40-entropy-23-00478" class="html-bibr">40</a>] and by the ALICE Collaboration at the LHC [<a href="#B41-entropy-23-00478" class="html-bibr">41</a>,<a href="#B42-entropy-23-00478" class="html-bibr">42</a>,<a href="#B43-entropy-23-00478" class="html-bibr">43</a>]. The solid and dashed curves are our results, fitted by using Equation (<a href="#FD10-entropy-23-00478" class="html-disp-formula">10</a>) due to Equations (7) and (9), with <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>μ</mi> <mi>B</mi> </msub> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math>, respectively.</p>
Full article ">Figure 4
<p>Dependences of (<b>a</b>) effective temperature <span class="html-italic">T</span>, (<b>b</b>) entropy index <span class="html-italic">q</span>, and (<b>c</b>) revised index <math display="inline"><semantics> <msub> <mi>a</mi> <mn>0</mn> </msub> </semantics></math> on energy <math display="inline"><semantics> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> </semantics></math>, where the closed and open symbols are cited from <a href="#entropy-23-00478-t001" class="html-table">Table 1</a> and <a href="#entropy-23-00478-t002" class="html-table">Table 2</a> which are obtained from the fittings with <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> (solid curves) and <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>μ</mi> <mi>B</mi> </msub> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math> (dashed curves) in <a href="#entropy-23-00478-f001" class="html-fig">Figure 1</a>, <a href="#entropy-23-00478-f002" class="html-fig">Figure 2</a> and <a href="#entropy-23-00478-f003" class="html-fig">Figure 3</a>, respectively. The triangles, circles, squares, and pentagrams represent the results for charged pions, kaons, protons, and the average by weighting different yields, respectively.</p>
Full article ">Figure 5
<p>Dependences of <span class="html-italic">T</span> on <math display="inline"><semantics> <msub> <mi>m</mi> <mn>0</mn> </msub> </semantics></math>. Different symbols represent the results from identified particles produced in central AA collisions at different energies shown in panels (<b>a</b>)–(<b>f</b>). The lines are the results fitted by the least square method, where the intercepts are regarded as <math display="inline"><semantics> <msub> <mi>T</mi> <mn>0</mn> </msub> </semantics></math>. The closed and open symbols (the solid and dashed curves) correspond to the results for <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>μ</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>μ</mi> <mi>B</mi> </msub> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math> respectively.</p>
Full article ">Figure 6
<p>Different symbols represent the results from identified particles produced in central AA collisions at different energies shown in panels (<b>a</b>)–(<b>f</b>). Same as for <a href="#entropy-23-00478-f005" class="html-fig">Figure 5</a>, but showing the dependences of <math display="inline"><semantics> <mrow> <mo>〈</mo> <msub> <mi>p</mi> <mi>T</mi> </msub> <mo>〉</mo> </mrow> </semantics></math> on <math display="inline"><semantics> <mover> <mi>m</mi> <mo>¯</mo> </mover> </semantics></math>. The lines are the results fitted by the least square method, where the slopes are regarded as <math display="inline"><semantics> <msub> <mi>β</mi> <mi>T</mi> </msub> </semantics></math>.</p>
Full article ">Figure 7
<p>Dependences of (<b>a</b>) <math display="inline"><semantics> <msub> <mi>T</mi> <mn>0</mn> </msub> </semantics></math> on <math display="inline"><semantics> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <msub> <mi>β</mi> <mi>T</mi> </msub> </semantics></math> on <math display="inline"><semantics> <msqrt> <msub> <mi>s</mi> <mrow> <mi>N</mi> <mi>N</mi> </mrow> </msub> </msqrt> </semantics></math>, and (<b>c</b>) <math display="inline"><semantics> <msub> <mi>T</mi> <mn>0</mn> </msub> </semantics></math> on <math display="inline"><semantics> <msub> <mi>β</mi> <mi>T</mi> </msub> </semantics></math>. The parameter values are obtained from <a href="#entropy-23-00478-t003" class="html-table">Table 3</a> and <a href="#entropy-23-00478-t004" class="html-table">Table 4</a>, which are from the linear fittings in <a href="#entropy-23-00478-f005" class="html-fig">Figure 5</a> and <a href="#entropy-23-00478-f006" class="html-fig">Figure 6</a>.</p>
Full article ">
14 pages, 1292 KiB  
Article
Hybrid Basketball Game Outcome Prediction Model by Integrating Data Mining Methods for the National Basketball Association
by Wei-Jen Chen, Mao-Jhen Jhou, Tian-Shyug Lee and Chi-Jie Lu
Entropy 2021, 23(4), 477; https://doi.org/10.3390/e23040477 - 17 Apr 2021
Cited by 20 | Viewed by 6734
Abstract
The sports market has grown rapidly over the last several decades. Sports outcomes prediction is an attractive sports analytic challenge as it provides useful information for operations in the sports market. In this study, a hybrid basketball game outcomes prediction scheme is developed [...] Read more.
The sports market has grown rapidly over the last several decades. Sports outcomes prediction is an attractive sports analytic challenge as it provides useful information for operations in the sports market. In this study, a hybrid basketball game outcomes prediction scheme is developed for predicting the final score of the National Basketball Association (NBA) games by integrating five data mining techniques, including extreme learning machine, multivariate adaptive regression splines, k-nearest neighbors, eXtreme gradient boosting (XGBoost), and stochastic gradient boosting. Designed features are generated by merging different game-lags information from fundamental basketball statistics and used in the proposed scheme. This study collected data from all the games of the NBA 2018–2019 seasons. There are 30 teams in the NBA and each team play 82 games per season. A total of 2460 NBA game data points were collected. Empirical results illustrated that the proposed hybrid basketball game prediction scheme achieves high prediction performance and identifies suitable game-lag information and relevant game features (statistics). Our findings suggested that a two-stage XGBoost model using four pieces of game-lags information achieves the best prediction performance among all competing models. The six designed features, including averaged defensive rebounds, averaged two-point field goal percentage, averaged free throw percentage, averaged offensive rebounds, averaged assists, and averaged three-point field goal attempts, from four game-lags have a greater effect on the prediction of final scores of NBA games than other game-lags. The findings of this study provide relevant insights and guidance for other team or individual sports outcomes prediction research. Full article
(This article belongs to the Special Issue Complex and Fractional Dynamics II)
Show Figures

Figure 1

Figure 1
<p>Flowchart of the proposed basketball game score prediction scheme.</p>
Full article ">Figure 2
<p>Example of the designed features for variable <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mrow> <mi>i</mi> <mo>,</mo> <mo> </mo> <mi>t</mi> </mrow> </msub> </mrow> </semantics></math> in different game-lags.</p>
Full article ">Figure 3
<p>Evaluation results of the selection of different numbers of important features for modeling the two-stage models: (<b>a</b>) T-ELM, (<b>b</b>) T-MARS, (<b>c</b>) T-XGBoost, (<b>d</b>) T-SGB, (<b>e</b>) T-KNN.</p>
Full article ">
16 pages, 7568 KiB  
Article
Refined Composite Multi-Scale Reverse Weighted Permutation Entropy and Its Applications in Ship-Radiated Noise
by Yuxing Li, Bo Geng and Shangbin Jiao
Entropy 2021, 23(4), 476; https://doi.org/10.3390/e23040476 - 17 Apr 2021
Cited by 10 | Viewed by 2494
Abstract
Ship-radiated noise is one of the important signal types under the complex ocean background, which can well reflect physical properties of ships. As one of the valid measures to characterize the complexity of ship-radiated noise, permutation entropy (PE) has the advantages of high [...] Read more.
Ship-radiated noise is one of the important signal types under the complex ocean background, which can well reflect physical properties of ships. As one of the valid measures to characterize the complexity of ship-radiated noise, permutation entropy (PE) has the advantages of high efficiency and simple calculation. However, PE has the problems of missing amplitude information and single scale. To address the two drawbacks, refined composite multi-scale reverse weighted PE (RCMRWPE), as a novel measurement technology of describing the signal complexity, is put forward based on refined composite multi-scale processing (RCMP) and reverse weighted PE (RWPE). RCMP is an improved method of coarse-graining, which not only solves the problem of single scale, but also improves the stability of traditional coarse-graining; RWPE has been proposed more recently, and has better inter-class separability and robustness performance to noise than PE, weighted PE (WPE), and reverse PE (RPE). Additionally, a feature extraction scheme of ship-radiated noise is proposed based on RCMRWPE, furthermore, RCMRWPE is combined with discriminant analysis classifier (DAC) to form a new classification method. After that, a large number of comparative experiments of feature extraction schemes and classification methods with two artificial random signals and six ship-radiated noise are carried out, which show that the proposed feature extraction scheme has better performance in distinguishing ability and stability than the other three similar feature extraction schemes based on multi-scale PE (MPE), multi-scale WPE (MWPE), and multi-scale RPE (MRPE), and the proposed classification method also has the highest recognition rate. Full article
(This article belongs to the Special Issue Multiscale Entropy Approaches and Their Applications II)
Show Figures

Figure 1

Figure 1
<p>The origin of RCMRWPE.</p>
Full article ">Figure 2
<p>A pattern in PE and the corresponding three possible patterns in RWPW.</p>
Full article ">Figure 3
<p>The functions <math display="inline"><semantics> <mrow> <msub> <mi>H</mi> <mrow> <mi>P</mi> <mi>E</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>H</mi> <mrow> <mi>R</mi> <mi>W</mi> <mi>P</mi> <mi>E</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>The functions <math display="inline"><semantics> <mrow> <mn>1</mn> <mo>−</mo> <msub> <mi>H</mi> <mrow> <mi>P</mi> <mi>E</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>H</mi> <mrow> <mi>R</mi> <mi>W</mi> <mi>P</mi> <mi>E</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>.</mo> </mrow> </semantics></math></p>
Full article ">Figure 5
<p>The coarse-graining for MPE and RCMP for RCMRWPE.</p>
Full article ">Figure 6
<p>The flow chart of the proposed classification method.</p>
Full article ">Figure 7
<p>The mean and STD entropy curves of WGN and <math display="inline"><semantics> <mrow> <mn>1</mn> <mo>/</mo> <mi>f</mi> </mrow> </semantics></math> noise.</p>
Full article ">Figure 8
<p>The normalized time-domain waveforms and probability density estimation function for six ship-radiated noise.</p>
Full article ">Figure 8 Cont.
<p>The normalized time-domain waveforms and probability density estimation function for six ship-radiated noise.</p>
Full article ">Figure 9
<p>The mean and STD entropy curves of different ship-radiated noise.</p>
Full article ">Figure 10
<p>The recognition rate of four feature extraction schemes based on DAC.</p>
Full article ">Figure 11
<p>The confusion matrix under the combination of features with the highest recognition rate for each classification method.</p>
Full article ">Figure 11 Cont.
<p>The confusion matrix under the combination of features with the highest recognition rate for each classification method.</p>
Full article ">
23 pages, 1239 KiB  
Article
Extended Lattice Boltzmann Model
by Mohammad Hossein Saadat, Benedikt Dorschner and Ilya Karlin
Entropy 2021, 23(4), 475; https://doi.org/10.3390/e23040475 - 17 Apr 2021
Cited by 10 | Viewed by 2837
Abstract
Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation [...] Read more.
Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation that restores Galilean invariance and the isotropy of the stress tensor by introducing an extended equilibrium. This modification extends lattice Boltzmann models to simulations with higher values of the flow velocity and can be used at temperatures that are higher than the lattice reference temperature, which enhances computational efficiency by decreasing the number of required time steps. Furthermore, the extended model also remains valid for stretched lattices, which are useful when flow gradients are predominant in one direction. The model is validated by simulations of two- and three-dimensional benchmark problems, including the double shear layer flow, the decay of homogeneous isotropic turbulence, the laminar boundary layer over a flat plate and the turbulent channel flow. Full article
(This article belongs to the Section Statistical Physics)
Show Figures

Figure 1

Figure 1
<p>Numerical measurement of viscosity for axis-aligned setup at temperature <math display="inline"><semantics> <mrow> <mi>T</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math> for different velocities. The exact solution corresponds to <math display="inline"><semantics> <mrow> <msub> <mi>ν</mi> <mrow> <mi>n</mi> <mi>u</mi> <mi>m</mi> </mrow> </msub> <mo>/</mo> <mi>ν</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>Numerical measurement of viscosity for axis-aligned setup at Mach number <math display="inline"><semantics> <mrow> <mi>M</mi> <mi>a</mi> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math> for different temperatures. The exact solution corresponds to <math display="inline"><semantics> <mrow> <msub> <mi>ν</mi> <mrow> <mi>n</mi> <mi>u</mi> <mi>m</mi> </mrow> </msub> <mo>/</mo> <mi>ν</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Numerical measurement of viscosity for rotated setup at temperature <math display="inline"><semantics> <mrow> <mi>T</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math> for different velocities and stretching ratios. The exact solution corresponds to <math display="inline"><semantics> <mrow> <msub> <mi>ν</mi> <mrow> <mi>n</mi> <mi>u</mi> <mi>m</mi> </mrow> </msub> <mo>/</mo> <mi>ν</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Velocity magnitude in lattice units for the decaying homogeneous isotropic turbulence at <math display="inline"><semantics> <mrow> <msub> <mi>Ma</mi> <mi mathvariant="normal">t</mi> </msub> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>Re</mi> <mi mathvariant="sans-serif">Λ</mi> </msub> <mo>=</mo> <mn>72</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msup> <mi>t</mi> <mo>*</mo> </msup> <mo>=</mo> <mn>1.0</mn> </mrow> </semantics></math> with temperature <math display="inline"><semantics> <mrow> <mi>T</mi> <mo>=</mo> <mn>0.55</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Time evolution of the turbulent kinetic energy for decaying isotropic turbulence at <math display="inline"><semantics> <mrow> <msub> <mi>Ma</mi> <mi mathvariant="normal">t</mi> </msub> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>Re</mi> <mi mathvariant="sans-serif">Λ</mi> </msub> <mo>=</mo> <mn>72</mn> </mrow> </semantics></math>. Lines: present model; symbol: DNS [<a href="#B36-entropy-23-00475" class="html-bibr">36</a>].</p>
Full article ">Figure 6
<p>Time evolution of the Taylor microscale Reynolds number for decaying isotropic turbulence at <math display="inline"><semantics> <mrow> <msub> <mi>Ma</mi> <mi mathvariant="normal">t</mi> </msub> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>Re</mi> <mi mathvariant="sans-serif">Λ</mi> </msub> <mo>=</mo> <mn>72</mn> </mrow> </semantics></math>. Lines: present model; symbol: direct numerical simulations (DNS) [<a href="#B36-entropy-23-00475" class="html-bibr">36</a>].</p>
Full article ">Figure 7
<p>Vorticity field for double shear layer flow at <math display="inline"><semantics> <mrow> <msup> <mi>t</mi> <mo>*</mo> </msup> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> with regular lattice (<b>left</b>) and stretched lattice (<b>right</b>). Vorticity magnitude is normalized by its maximum value.</p>
Full article ">Figure 8
<p>Evolution of kinetic energy (<b>left</b>) and enstrophy (<b>right</b>) for double shear layer flow at <math display="inline"><semantics> <mrow> <mi>R</mi> <mi>e</mi> <mo>=</mo> <msup> <mn>10</mn> <mn>4</mn> </msup> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>Comparison of the velocity profile at <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <msub> <mi>L</mi> <mi>x</mi> </msub> </mrow> </semantics></math> for flow over a flat plate at different stretching ratios. Lines: present model; symbols: Blasius solution.</p>
Full article ">Figure 10
<p>Comparison of the skin friction coefficient for flow over a flat plate at different stretching ratio. Lines: present model; symbols: analytical solution.</p>
Full article ">Figure 11
<p>Snapshot of the velocity magnitude in lattice units for turbulent channel flow at <math display="inline"><semantics> <mrow> <msub> <mi>Re</mi> <mi>τ</mi> </msub> <mo>=</mo> <mn>180</mn> </mrow> </semantics></math> with <math display="inline"><semantics> <mrow> <msub> <mi>λ</mi> <mi>x</mi> </msub> <mo>=</mo> <mn>1.4</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 12
<p>Comparison of the mean velocity profile in a turbulent channel flow at <math display="inline"><semantics> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>τ</mi> </msub> <mo>=</mo> <mn>180</mn> </mrow> </semantics></math> with <math display="inline"><semantics> <mrow> <msub> <mi>λ</mi> <mi>x</mi> </msub> <mo>=</mo> <mn>1.4</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 13
<p>Comparison of the rms of the velocity fluctuations in a turbulent channel flow at <math display="inline"><semantics> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>τ</mi> </msub> <mo>=</mo> <mn>180</mn> </mrow> </semantics></math> with <math display="inline"><semantics> <mrow> <msub> <mi>λ</mi> <mi>x</mi> </msub> <mo>=</mo> <mn>1.4</mn> </mrow> </semantics></math>. Symbols: present model; lines: DNS [<a href="#B40-entropy-23-00475" class="html-bibr">40</a>].</p>
Full article ">Figure A1
<p>Spectral dissipation of acoustic modes for different models. Red symbols: LBGK; black symbols: extended LBM (<a href="#FD33-entropy-23-00475" class="html-disp-formula">33</a>); blue symbols: LC LBM [<a href="#B9-entropy-23-00475" class="html-bibr">9</a>]; dashed line: Navier–Stokes. The velocity and temperature are set to <math display="inline"><semantics> <mrow> <msub> <mi>u</mi> <mi>x</mi> </msub> <mo>=</mo> <mn>0.3</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>T</mi> <mo>=</mo> <msub> <mi>T</mi> <mi>L</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure A2
<p>Comparison of density profile for shock tube problem at density ratio <math display="inline"><semantics> <mrow> <msub> <mi>ρ</mi> <mi>l</mi> </msub> <mo>/</mo> <msub> <mi>ρ</mi> <mi>r</mi> </msub> <mo>=</mo> <mn>3</mn> </mrow> </semantics></math>, after 500 iterations. Solid line: LBGK; dashed line: extended LBM (<a href="#FD33-entropy-23-00475" class="html-disp-formula">33</a>); symbols: LC LBM [<a href="#B9-entropy-23-00475" class="html-bibr">9</a>].</p>
Full article ">Figure A3
<p>Mach number profile, <math display="inline"><semantics> <mrow> <mi>Ma</mi> <mo>=</mo> <mi>u</mi> <mo>/</mo> <msqrt> <msub> <mi>T</mi> <mi>L</mi> </msub> </msqrt> </mrow> </semantics></math>, for the shock tube problem at density ratio <math display="inline"><semantics> <mrow> <msub> <mi>ρ</mi> <mi>l</mi> </msub> <mo>/</mo> <msub> <mi>ρ</mi> <mi>r</mi> </msub> <mo>=</mo> <mn>3</mn> </mrow> </semantics></math>, after 500 iterations. Solid line: LBGK; dashed line: extended LBM (<a href="#FD33-entropy-23-00475" class="html-disp-formula">33</a>); symbols: LC LBM [<a href="#B9-entropy-23-00475" class="html-bibr">9</a>].</p>
Full article ">
17 pages, 3153 KiB  
Article
Low-Frequency Seismic Noise Properties in the Japanese Islands
by Alexey Lyubushin
Entropy 2021, 23(4), 474; https://doi.org/10.3390/e23040474 - 16 Apr 2021
Cited by 12 | Viewed by 4912
Abstract
The records of seismic noise in Japan for the period of 1997–2020, which includes the Tohoku seismic catastrophe on 11 March 2011, are considered. The following properties of noise are analyzed: The wavelet-based Donoho–Johnston index, the singularity spectrum support width, and the entropy [...] Read more.
The records of seismic noise in Japan for the period of 1997–2020, which includes the Tohoku seismic catastrophe on 11 March 2011, are considered. The following properties of noise are analyzed: The wavelet-based Donoho–Johnston index, the singularity spectrum support width, and the entropy of the wavelet coefficients. The question of whether precursors of strong earthquakes can be formulated on their basis is investigated. Attention is paid to the time interval after the Tohoku mega-earthquake to the trends in the mean properties of low-frequency seismic noise, which reflect the constant simplification of the statistical structure of seismic vibrations. Estimates of two-dimensional probability densities of extreme values are presented, which highlight the places in which extreme values of seismic noise properties are most often realized. The estimates of the probability densities of extreme values coincide with each other and have a maximum in the region: 30° N  Lat  34° N, 136° E  Lon 140° E. The main conclusions of the conducted studies are that the preparation of a strong earthquake is accompanied by a simplification of the structure of seismic noise. It is shown that bursts of coherence between the time series of the day length and the noise properties within annual time window precede bursts of released seismic energy. The value of the lag in the release of seismic energy relative to bursts of coherence is about 1.5 years, which can be used to declare a time interval of high seismic hazard after reaching the peak of coherence. Full article
(This article belongs to the Special Issue Complex Systems Time Series Analysis and Modeling for Geoscience)
Show Figures

Figure 1

Figure 1
<p>Positions of 78 seismic stations in Japan. The blue dashed line shows the position of the Nankai Deep Trench.</p>
Full article ">Figure 2
<p>Daily number of working stations.</p>
Full article ">Figure 3
<p>Averaged maps of seismic noise parameters <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>α</mi> </mrow> </semantics></math> (<b>a1</b>–<b>a3</b>), <math display="inline"><semantics> <mrow> <mi>E</mi> <mi>n</mi> </mrow> </semantics></math> (<b>b1</b>–<b>b3</b>) and <math display="inline"><semantics> <mi>γ</mi> </semantics></math> (<b>c1</b>–<b>c3</b>), calculated for three time intervals (1 January 1997–25 September 2003, 26 September 2003–10 October 2011, and 14 March 2011–31 March 2021). Black stars indicate hypocenters of two strong earthquakes: 25 September 2003, M = 8.3 and 11 March 2011, M = 9.1. The spatial distribution of seismic noise properties is shown only in the vicinity of the Japanese Islands in the union of circles with a radius of 250 km, built around each seismic station.</p>
Full article ">Figure 4
<p>Left panel presents graphs of daily median values of seismic noise parameters <math display="inline"><semantics> <mi>γ</mi> </semantics></math>, <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>α</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>E</mi> <mi>n</mi> </mrow> </semantics></math>, green lines are graphs of running average of the length 57 days. Right panel presents corresponding results of deep smoothing of median daily values by Gaussian kernel with bandwidth 182 days, red lines present linear trends of smoothed values after 11 March 2011. Horizontal blue lines present mean values of seismic noise statistics for three time intervals: 1 January 1997–25 September 2003, 26 September 2003–3 October 2011, and 14 March 2011–31 March 2021.</p>
Full article ">Figure 5
<p>Maps of probability densities of extreme values, minimums for <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>α</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mi>γ</mi> </semantics></math> and maximums <math display="inline"><semantics> <mrow> <mi>E</mi> <mi>n</mi> </mrow> </semantics></math>, estimated for three time intervals (1 January 1997–25 September 2003, 26 September 2003–10 October 2011, and 14 March 2011–31 March 2021). Black stars indicate hypocenters of two strong earthquakes: 25 September 2003, M = 8.3 and 11 March 2011, M = 9.1. The distribution of probability densities of extreme values of seismic noise properties is shown only in the vicinity of the Japanese Islands in the union of circles with a radius of 250 km, built around each seismic station.</p>
Full article ">Figure 6
<p>(<b>a</b>) Plot of length of day (LOD); (<b>b</b>) graph of first principal component of three seismic noise parameters calculated within the moving time window of 365 days; green line is a graph of running average of the length of 57 days; (<b>c</b>) time-frequency map of quadratic coherence between LOD and principal component (<b>b</b>) in a moving time window of the length of 365 days with a mutual shift of 3 days; (<b>d</b>) graph of the maximum values with respect to the frequencies of the squared coherence between LOD and the first principal component; (<b>e</b>) plot of the decimal logarithm of seismic energy released in the vicinity of the Japan Islands; (<b>f</b>) correlation function between the values of the logarithm of the released seismic energy and the maximums of coherence between the day length and the first principal component. Negative values of time shifts on graph (<b>f</b>) correspond to the lag in the release of seismic energy relative to bursts of coherence between LOD and seismic noise first principal component. Graphs (<b>d</b>) and (<b>e</b>) are plotted in dependence of the right end of time window of 365 days with an offset of 3 days.</p>
Full article ">
24 pages, 8988 KiB  
Article
Estimation of Feeding Composition of Industrial Process Based on Data Reconciliation
by Yusi Luan, Mengxuan Jiang, Zhenxiang Feng and Bei Sun
Entropy 2021, 23(4), 473; https://doi.org/10.3390/e23040473 - 16 Apr 2021
Cited by 5 | Viewed by 2131
Abstract
For an industrial process, the estimation of feeding composition is important for analyzing production status and making control decisions. However, random errors or even gross ones inevitably contaminate the actual measurements. Feeding composition is conventionally obtained via discrete and low-rate artificial testing. To [...] Read more.
For an industrial process, the estimation of feeding composition is important for analyzing production status and making control decisions. However, random errors or even gross ones inevitably contaminate the actual measurements. Feeding composition is conventionally obtained via discrete and low-rate artificial testing. To address these problems, a feeding composition estimation approach based on data reconciliation procedure is developed. To improve the variable accuracy, a novel robust M-estimator is first proposed. Then, an iterative robust hierarchical data reconciliation and estimation strategy is applied to estimate the feeding composition. The feasibility and effectiveness of the estimation approach are verified on a fluidized bed roaster. The proposed M-estimator showed better overall performance. Full article
(This article belongs to the Special Issue Complex Dynamic System Modelling, Identification and Control)
Show Figures

Figure 1

Figure 1
<p>Characteristicfunctions of three estimators at distinct relative efficiency levels.</p>
Full article ">Figure 1 Cont.
<p>Characteristicfunctions of three estimators at distinct relative efficiency levels.</p>
Full article ">Figure 2
<p>Objective functions for Xie, Welsch, Cauchy, Fair and Proposed estimator.</p>
Full article ">Figure 3
<p>Influence functions for Xie, Welsch, Cauchy, Fair and Proposed estimator.</p>
Full article ">Figure 4
<p>Diagram of the measurement network.</p>
Full article ">Figure 5
<p>Statistical diagram of gross errors detection process for distinct methods.</p>
Full article ">Figure 5 Cont.
<p>Statistical diagram of gross errors detection process for distinct methods.</p>
Full article ">Figure 6
<p>Results of <span class="html-italic">OP</span> and <span class="html-italic">AVTI</span> with different methods.</p>
Full article ">Figure 7
<p>Structure diagram of fluidized bed roaster.</p>
Full article ">Figure 8
<p>The framework of iterative robust hierarchical data reconciliation and composition estimation.</p>
Full article ">Figure 9
<p>Reconciled results for the first feed material quantity.</p>
Full article ">Figure 10
<p>Reconciled results for the second feed material quantity.</p>
Full article ">Figure 11
<p>Reconciled results for the blast velocity.</p>
Full article ">Figure 12
<p>Reconciled results for the temperature of blast.</p>
Full article ">Figure 13
<p>Reconciled results for the temperature of gas.</p>
Full article ">Figure 14
<p>Reconciled results for the temperature of calcine.</p>
Full article ">Figure 15
<p>The standard deviation of reconciled results and measurements.</p>
Full article ">Figure 16
<p>The percentage of ZnS before data reconciliation.</p>
Full article ">Figure 17
<p>The percentage of ZnS after data reconciliation.</p>
Full article ">
11 pages, 2171 KiB  
Article
Empirical Mode Decomposition-Derived Entropy Features Are Beneficial to Distinguish Elderly People with a Falling History on a Force Plate Signal
by Li-Wei Chou, Kang-Ming Chang, Yi-Chun Wei and Mei-Kuei Lu
Entropy 2021, 23(4), 472; https://doi.org/10.3390/e23040472 - 16 Apr 2021
Cited by 8 | Viewed by 2524
Abstract
Fall risk prediction is an important issue for the elderly. A center of pressure signal, derived from a force plate, is useful for the estimation of body calibration. However, it is still difficult to distinguish elderly people’s fall history by using a force [...] Read more.
Fall risk prediction is an important issue for the elderly. A center of pressure signal, derived from a force plate, is useful for the estimation of body calibration. However, it is still difficult to distinguish elderly people’s fall history by using a force plate signal. In this study, older adults with and without a history of falls were recruited to stand still for 60 s on a force plate. Forces in the x, y and z directions (Fx, Fy, and Fz) and center of pressure in the anteroposterior (COPx) and mediolateral directions (COPy) were derived. There were 49 subjects in the non-fall group, with an average age of 71.67 (standard derivation: 6.56). There were also 27 subjects in the fall group, with an average age of 70.66 (standard derivation: 6.38). Five signal series—forces in x, y, z (Fx, Fy, Fz), COPX, and COPy directions—were used. These five signals were further decomposed with empirical mode decomposition (EMD) with seven intrinsic mode functions. Time domain features (mean, standard derivation and coefficient of variations) and entropy features (approximate entropy and sample entropy) of the original signals and EMD-derived signals were extracted. Results showed that features extracted from the raw COP data did not differ significantly between the fall and non-fall groups. There were 10 features extracted using EMD, with significant differences observed among fall and non-fall groups. These included four features from COPx and two features from COPy, Fx and Fz. Full article
(This article belongs to the Special Issue Entropy in Biomedical Applications)
Show Figures

Figure 1

Figure 1
<p>Academic papers related to center of pressure. The x-axis is publication years; the y-axis is the number of publications. Data are retrieved from PubMed.</p>
Full article ">Figure 2
<p>Center of pressure (COP) data and related intrinsic mode function (IMF) decompositions.</p>
Full article ">Figure 2 Cont.
<p>Center of pressure (COP) data and related intrinsic mode function (IMF) decompositions.</p>
Full article ">Figure 2 Cont.
<p>Center of pressure (COP) data and related intrinsic mode function (IMF) decompositions.</p>
Full article ">Figure 3
<p>Experiment flowchart.</p>
Full article ">
18 pages, 786 KiB  
Article
Effect of Inter-System Coupling on Heat Transport in a Microscopic Collision Model
by Feng Tian, Jian Zou, Lei Li, Hai Li and Bin Shao
Entropy 2021, 23(4), 471; https://doi.org/10.3390/e23040471 - 16 Apr 2021
Cited by 5 | Viewed by 2278
Abstract
In this paper we consider a bipartite system composed of two subsystems each coupled to its own thermal environment. Based on a collision model, we mainly study whether the approximation (i.e., the inter-system coupling is ignored when modeling the system–environment interaction) is valid [...] Read more.
In this paper we consider a bipartite system composed of two subsystems each coupled to its own thermal environment. Based on a collision model, we mainly study whether the approximation (i.e., the inter-system coupling is ignored when modeling the system–environment interaction) is valid or not. We also address the problem of heat transport unitedly for both excitation-conserving system–environment interactions and non-excitation-conserving system–environment interactions. For the former interaction, as the inter-system interaction strength increases, at first this approximation gets worse as expected, but then counter-intuitively gets better even for a stronger inter-system coupling. For the latter interaction with asymmetry, this approximation gets progressively worse. In this case we realize a perfect thermal rectification, and we cannot find an apparent rectification effect for the former interaction. Finally and more importantly, our results show that whether this approximation is valid or not is closely related to the quantum correlations between the subsystems, i.e., the weaker the quantum correlations, the more justified the approximation and vice versa. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic sketch of a bipartite system <span class="html-italic">S</span> made up of two interacting subsystems connected to two independent subenvironments. In the <span class="html-italic">n</span>th round of the dynamics, after a free evolution of the whole system, <math display="inline"><semantics> <msub> <mi>S</mi> <mn>1</mn> </msub> </semantics></math> interacts with <math display="inline"><semantics> <msubsup> <mi>E</mi> <mi>n</mi> <mn>1</mn> </msubsup> </semantics></math> and next <math display="inline"><semantics> <msub> <mi>S</mi> <mn>2</mn> </msub> </semantics></math> interacts with <math display="inline"><semantics> <msubsup> <mi>E</mi> <mi>n</mi> <mn>2</mn> </msubsup> </semantics></math>. In the <math display="inline"><semantics> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </semantics></math>th round of the dynamics, also after a free evolution of the whole system, <math display="inline"><semantics> <msub> <mi>S</mi> <mn>1</mn> </msub> </semantics></math> interacts with <math display="inline"><semantics> <msubsup> <mi>E</mi> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>1</mn> </msubsup> </semantics></math> and next <math display="inline"><semantics> <msub> <mi>S</mi> <mn>2</mn> </msub> </semantics></math> interacts with <math display="inline"><semantics> <msubsup> <mi>E</mi> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> </semantics></math>. The system then moves to the <math display="inline"><semantics> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>2</mn> <mo>)</mo> </mrow> </semantics></math>th round and this process is repeated over and over.</p>
Full article ">Figure 2
<p>Steady heat currents <math display="inline"><semantics> <msub> <mi>J</mi> <mi>h</mi> </msub> </semantics></math> as a function of <math display="inline"><semantics> <mi>δ</mi> </semantics></math> for different <math display="inline"><semantics> <mi>γ</mi> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.2</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (blue), <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.5</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (red), and <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.8</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (green).The solid and dashed lines correspond to cases without and with the decoupling approximation, respectively. The inset is the magnified <math display="inline"><semantics> <msub> <mi>J</mi> <mi>h</mi> </msub> </semantics></math> for <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.2</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>. In both cases the system is initialized in <math display="inline"><semantics> <mrow> <mo>|</mo> <mn>11</mn> <mo>〉</mo> </mrow> </semantics></math> and each ancilla is initialized in its thermal state. The plots are obtained for <math display="inline"><semantics> <mrow> <msub> <mi>ω</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>ω</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>ω</mi> <mn>0</mn> </msub> <mi>τ</mi> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>5</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Trace distance <math display="inline"><semantics> <msub> <mi mathvariant="script">D</mi> <mi mathvariant="script">T</mi> </msub> </semantics></math> between the density matrices obtained from cases with and without the decoupling approximation against <math display="inline"><semantics> <mi>δ</mi> </semantics></math> at steady state. The blue line, the red line and the green line correspond to <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.2</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.5</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.8</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, respectively. The other parameters are the same as those in <a href="#entropy-23-00471-f002" class="html-fig">Figure 2</a>. Inset shows <math display="inline"><semantics> <msup> <mi>δ</mi> <mo>′</mo> </msup> </semantics></math> for the maximum of <math display="inline"><semantics> <msub> <mi mathvariant="script">D</mi> <mi mathvariant="script">T</mi> </msub> </semantics></math> as a function of <math display="inline"><semantics> <mi>γ</mi> </semantics></math>.</p>
Full article ">Figure 4
<p>Quantum discord <math display="inline"><semantics> <mi mathvariant="script">D</mi> </semantics></math> of the qubit system as a function of <math display="inline"><semantics> <mi>δ</mi> </semantics></math>, without approximation. The blue line and red line correspond to <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.2</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.5</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, respectively. Inset a is the magnified <math display="inline"><semantics> <mi mathvariant="script">D</mi> </semantics></math> for <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.2</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>. Inset b corresponds to <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.8</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>. The remaining parameters are the same as those in <a href="#entropy-23-00471-f002" class="html-fig">Figure 2</a>.</p>
Full article ">Figure 5
<p><math display="inline"><semantics> <msub> <mi>J</mi> <mi>h</mi> </msub> </semantics></math> as a function of <math display="inline"><semantics> <mi>δ</mi> </semantics></math> with various temperature difference without the decoupling approximation. We use solid line for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and different <math display="inline"><semantics> <msub> <mi>T</mi> <mn>2</mn> </msub> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.1</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (green), <math display="inline"><semantics> <mrow> <mn>4</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (red), and <math display="inline"><semantics> <mrow> <mn>8</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (blue); dashed line for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and different <math display="inline"><semantics> <msub> <mi>T</mi> <mn>1</mn> </msub> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>0.1</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (green), <math display="inline"><semantics> <mrow> <mn>4</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (red), and <math display="inline"><semantics> <mrow> <mn>8</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (blue). The two qubits are off-resonant with <math display="inline"><semantics> <mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>ω</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>ω</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>. Here we set <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.3</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, and the other parameters are the same as those in <a href="#entropy-23-00471-f002" class="html-fig">Figure 2</a>. Insets a, b, and c show the corresponding rectification factor <span class="html-italic">R</span> for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and different <math display="inline"><semantics> <msub> <mi>T</mi> <mn>2</mn> </msub> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>8</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mn>4</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mn>0.1</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, respectively.</p>
Full article ">Figure 6
<p><math display="inline"><semantics> <msub> <mi>J</mi> <mi>h</mi> </msub> </semantics></math> as a function of <math display="inline"><semantics> <mi>δ</mi> </semantics></math> with various temperature difference without the decoupling approximation. The black dashed line corresponds to <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>. We use blue line for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>5</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>; red and green lines for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>5</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mn>0.05</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, respectively. The other parameters are the same as those in <a href="#entropy-23-00471-f005" class="html-fig">Figure 5</a>.</p>
Full article ">Figure 7
<p>Trace distance <math display="inline"><semantics> <msub> <mi mathvariant="script">D</mi> <mi mathvariant="script">T</mi> </msub> </semantics></math> between the density matrices obtained from cases with and without the decoupling approximation against <math display="inline"><semantics> <mi>δ</mi> </semantics></math> for steady state, and the corresponding plot is obtained for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and different <math display="inline"><semantics> <msub> <mi>T</mi> <mn>2</mn> </msub> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.1</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (green), <math display="inline"><semantics> <mrow> <mn>4</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (red), and <math display="inline"><semantics> <mrow> <mn>8</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (blue).</p>
Full article ">Figure 8
<p>Quantum discord <math display="inline"><semantics> <mi mathvariant="script">D</mi> </semantics></math> of the qubit system as a function of <math display="inline"><semantics> <mi>δ</mi> </semantics></math>, without the decoupling approximation. The plots are obtained for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and different <math display="inline"><semantics> <msub> <mi>T</mi> <mn>2</mn> </msub> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.1</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (green), <math display="inline"><semantics> <mrow> <mn>4</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (red), and <math display="inline"><semantics> <mrow> <mn>8</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> (blue). The remaining parameters are the same as those in <a href="#entropy-23-00471-f005" class="html-fig">Figure 5</a>.</p>
Full article ">Figure 9
<p><math display="inline"><semantics> <msub> <mi>J</mi> <mi>h</mi> </msub> </semantics></math> as a function of inter-system coupling strength <math display="inline"><semantics> <mi>δ</mi> </semantics></math>. The solid and dashed lines correspond to without and with the decoupling approximation, respectively. We use red line for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>. Blue line is for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> in the case without the approximation. Here we set <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.5</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>. The other parameters are the same as those in <a href="#entropy-23-00471-f002" class="html-fig">Figure 2</a>. Inset shows corresponding rectification factor in the case without the decoupling approximation.</p>
Full article ">Figure 10
<p>Quantum discord <math display="inline"><semantics> <mi mathvariant="script">D</mi> </semantics></math> of the qubit system as a function of <math display="inline"><semantics> <mi>δ</mi> </semantics></math>, without the decoupling approximation. The plot is obtained for <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10</mn> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>ω</mi> <mn>0</mn> </msub> </mrow> </semantics></math>. The remaining parameters are the same as those in <a href="#entropy-23-00471-f009" class="html-fig">Figure 9</a>.</p>
Full article ">
12 pages, 1478 KiB  
Article
Applying the Horizontal Visibility Graph Method to Study Irreversibility of Electromagnetic Turbulence in Non-Thermal Plasmas
by Belén Acosta-Tripailao, Denisse Pastén and Pablo S. Moya
Entropy 2021, 23(4), 470; https://doi.org/10.3390/e23040470 - 16 Apr 2021
Cited by 17 | Viewed by 3120
Abstract
One of the fundamental open questions in plasma physics is the role of non-thermal particles distributions in poorly collisional plasma environments, a system that is commonly found throughout the Universe, e.g., the solar wind and the Earth’s magnetosphere correspond to natural plasma physics [...] Read more.
One of the fundamental open questions in plasma physics is the role of non-thermal particles distributions in poorly collisional plasma environments, a system that is commonly found throughout the Universe, e.g., the solar wind and the Earth’s magnetosphere correspond to natural plasma physics laboratories in which turbulent phenomena can be studied. Our study perspective is born from the method of Horizontal Visibility Graph (HVG) that has been developed in the last years to analyze time series avoiding the tedium and the high computational cost that other methods offer. Here, we build a complex network based on directed HVG technique applied to magnetic field fluctuations time series obtained from Particle In Cell (PIC) simulations of a magnetized collisionless plasma to distinguish the degree distributions and calculate the Kullback–Leibler Divergence (KLD) as a measure of relative entropy of data sets produced by processes that are not in equilibrium. First, we analyze the connectivity probability distribution for the undirected version of HVG finding how the Kappa distribution for low values of κ tends to be an uncorrelated time series, while the Maxwell–Boltzmann distribution shows a correlated stochastic processes behavior. Subsequently, we investigate the degree of temporary irreversibility of magnetic fluctuations that are self-generated by the plasma, comparing the case of a thermal plasma (described by a Maxwell–Botzmann velocity distribution function) with non-thermal Kappa distributions. We have shown that the KLD associated to the HVG is able to distinguish the level of reversibility that is associated to the thermal equilibrium in the plasma, because the dissipative degree of the system increases as the value of κ parameter decreases and the distribution function departs from the Maxwell–Boltzmann equilibrium. Full article
Show Figures

Figure 1

Figure 1
<p>Construction of Horizontal Visibility Graph. Top, a time series where the degree <math display="inline"><semantics> <msub> <mi>k</mi> <mi>in</mi> </msub> </semantics></math> for in-going links and <math display="inline"><semantics> <msub> <mi>k</mi> <mi>out</mi> </msub> </semantics></math> for out-going links of each of the <math display="inline"><semantics> <mrow> <mi>n</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> nodes are detailed. Bottom, probability distribution <span class="html-italic">P</span> in relation to degree <span class="html-italic">k</span>, where <math display="inline"><semantics> <msub> <mi>n</mi> <mi>in</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>n</mi> <mi>out</mi> </msub> </semantics></math> correspond to the frequency of appearance of the degrees <math display="inline"><semantics> <msub> <mi>k</mi> <mi>in</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>k</mi> <mi>out</mi> </msub> </semantics></math>, respectively, defining <math display="inline"><semantics> <msub> <mi>P</mi> <mi>in</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>P</mi> <mi>out</mi> </msub> </semantics></math>.</p>
Full article ">Figure 2
<p>(<b>Left</b>) Average magnetic field energy density fluctuations <math display="inline"><semantics> <msup> <mrow> <mo>(</mo> <mi>δ</mi> <mi>B</mi> <mo>/</mo> <msub> <mi>B</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </semantics></math> as a function of time obtained from Particle In Cell (PIC) simulations for Maxwell–Boltzmann (where MB represents <math display="inline"><semantics> <mrow> <msub> <mi>κ</mi> <mi>e</mi> </msub> <mo>→</mo> <mo>∞</mo> </mrow> </semantics></math>) and Kappa distributions considering different values of the <math display="inline"><semantics> <msub> <mi>κ</mi> <mi>e</mi> </msub> </semantics></math> parameter. (<b>Right</b>) Detrended average magnetic field energy density magnitude.</p>
Full article ">Figure 3
<p>Semi-log plot of the degree distributions of HVG associated to Kappa and Maxwell–Boltzmann distribution. There is an exponential behavior <math display="inline"><semantics> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>∼</mo> <mo form="prefix">exp</mo> <mfenced separators="" open="(" close=")"> <mo>−</mo> <mi>γ</mi> <mi>k</mi> </mfenced> </mrow> </semantics></math> and the <math display="inline"><semantics> <mi>γ</mi> </semantics></math> value is shown for each distribution. The left panel corresponds to the results for the magnetic field of the trend data from <a href="#entropy-23-00470-f002" class="html-fig">Figure 2</a> (left), while the right panel for the detrended data from <a href="#entropy-23-00470-f002" class="html-fig">Figure 2</a> (right).</p>
Full article ">Figure 4
<p>KL-Divergence (<span class="html-italic">D</span>) of magnetic field for different Kappa distributions. (<b>Left</b>) Horizontal Visibility Graph (HVG) method applied on the original data. (<b>Right</b>) HVG on the detrended data. The technique used to determine whether the data represent a reversible process consists of applying the HVG algorithm to randomly disordered copies of the data, obtaining the standard deviation <math display="inline"><semantics> <mi>σ</mi> </semantics></math> around the average divergence computed using the disordered data (black dot and vertical lines).</p>
Full article ">Figure 5
<p>Temporal evolution of the KL-divergence considering a moving window that covers 8000 data overlapping every 1000 data on the magnetic time series. (<b>Left</b>) HVG method applied on the original data and (<b>Right</b>) on the detrended data.</p>
Full article ">
11 pages, 922 KiB  
Article
Comparison between Highly Complex Location Models and GAMLSS
by Thiago G. Ramires, Luiz R. Nakamura, Ana J. Righetto, Renan J. Carvalho, Lucas A. Vieira and Carlos A. B. Pereira
Entropy 2021, 23(4), 469; https://doi.org/10.3390/e23040469 - 16 Apr 2021
Cited by 6 | Viewed by 2982
Abstract
This paper presents a discussion regarding regression models, especially those belonging to the location class. Our main motivation is that, with simple distributions having simple interpretations, in some cases, one gets better results than the ones obtained with overly complex distributions. For instance, [...] Read more.
This paper presents a discussion regarding regression models, especially those belonging to the location class. Our main motivation is that, with simple distributions having simple interpretations, in some cases, one gets better results than the ones obtained with overly complex distributions. For instance, with the reverse Gumbel (RG) distribution, it is possible to explain response variables by making use of the generalized additive models for location, scale, and shape (GAMLSS) framework, which allows the fitting of several parameters (characteristics) of the probabilistic distributions, like mean, mode, variance, and others. Three real data applications are used to compare several location models against the RG under the GAMLSS framework. The intention is to show that the use of a simple distribution (e.g., RG) based on a more sophisticated regression structure may be preferable than using a more complex location model. Full article
Show Figures

Figure 1

Figure 1
<p>Densities of <span class="html-italic">y</span> for each voltage level, disregarding censored observations.</p>
Full article ">Figure 2
<p>For the voltage data: (<b>a</b>) the estimated and empirical survival function from the generalized additive models for location, scale, and shape (GAMLSS) model based on the reverse Gumbel (RG) distribution considering smoothing functions and (<b>b</b>) the worm plot (WP).</p>
Full article ">Figure 3
<p>Densities of <span class="html-italic">y</span> for each temperature level.</p>
Full article ">Figure 4
<p>For class-H data: (<b>a</b>) the estimated and empirical survival function from RG and (<b>b</b>) the WP.</p>
Full article ">Figure 5
<p>Dispersion plot and boxplots for heart data as a function of the explanatory variables.</p>
Full article ">Figure 6
<p>Worm plot of the fitted fully parametrc GAMLSS model based on the RG distribution.</p>
Full article ">
10 pages, 876 KiB  
Article
The Quality of Statistical Reporting and Data Presentation in Predatory Dental Journals Was Lower Than in Non-Predatory Journals
by Pentti Nieminen and Sergio E. Uribe
Entropy 2021, 23(4), 468; https://doi.org/10.3390/e23040468 - 16 Apr 2021
Cited by 9 | Viewed by 4130
Abstract
Proper peer review and quality of published articles are often regarded as signs of reliable scientific journals. The aim of this study was to compare whether the quality of statistical reporting and data presentation differs among articles published in ‘predatory dental journals’ and [...] Read more.
Proper peer review and quality of published articles are often regarded as signs of reliable scientific journals. The aim of this study was to compare whether the quality of statistical reporting and data presentation differs among articles published in ‘predatory dental journals’ and in other dental journals. We evaluated 50 articles published in ‘predatory open access (OA) journals’ and 100 clinical trials published in legitimate dental journals between 2019 and 2020. The quality of statistical reporting and data presentation of each paper was assessed on a scale from 0 (poor) to 10 (high). The mean (SD) quality score of the statistical reporting and data presentation was 2.5 (1.4) for the predatory OA journals, 4.8 (1.8) for the legitimate OA journals, and 5.6 (1.8) for the more visible dental journals. The mean values differed significantly (p < 0.001). The quality of statistical reporting of clinical studies published in predatory journals was found to be lower than in open access and highly cited journals. This difference in quality is a wake-up call to consume study results critically. Poor statistical reporting indicates wider general lower quality in publications where the authors and journals are less likely to be critiqued by peer review. Full article
(This article belongs to the Special Issue Statistical Methods for Medicine and Health Sciences)
Show Figures

Figure 1

Figure 1
<p>Scatter plot of statistical reporting and data presentation (SRDP) score with CONSORT score in 50 clinical trial articles published in dental open access journals.</p>
Full article ">Figure 2
<p>Distribution of statistical reporting and data presentation quality score (SRDP) by journal group.</p>
Full article ">
21 pages, 1153 KiB  
Article
Toward a Comparison of Classical and New Privacy Mechanism
by Daniel Heredia-Ductram, Miguel Nunez-del-Prado and Hugo Alatrista-Salas
Entropy 2021, 23(4), 467; https://doi.org/10.3390/e23040467 - 15 Apr 2021
Cited by 1 | Viewed by 2698
Abstract
In the last decades, the development of interconnectivity, pervasive systems, citizen sensors, and Big Data technologies allowed us to gather many data from different sources worldwide. This phenomenon has raised privacy concerns around the globe, compelling states to enforce data protection laws. In [...] Read more.
In the last decades, the development of interconnectivity, pervasive systems, citizen sensors, and Big Data technologies allowed us to gather many data from different sources worldwide. This phenomenon has raised privacy concerns around the globe, compelling states to enforce data protection laws. In parallel, privacy-enhancing techniques have emerged to meet regulation requirements allowing companies and researchers to exploit individual data in a privacy-aware way. Thus, data curators need to find the most suitable algorithms to meet a required trade-off between utility and privacy. This crucial task could take a lot of time since there is a lack of benchmarks on privacy techniques. To fill this gap, we compare classical approaches of privacy techniques like Statistical Disclosure Control and Differential Privacy techniques to more recent techniques such as Generative Adversarial Networks and Machine Learning Copies using an entire commercial database in the current effort. The obtained results allow us to show the evolution of privacy techniques and depict new uses of the privacy-aware Machine Learning techniques. Full article
(This article belongs to the Special Issue Machine Learning Ecosystems: Opportunities and Threats)
Show Figures

Figure 1

Figure 1
<p>Spatial representation of dataset.</p>
Full article ">Figure 2
<p>Feature correlation heatmap for the raw dataset.</p>
Full article ">Figure 3
<p>Evolution of the (<b>a</b>) Information Loss and the (<b>b</b>) Disclosure Risk with Noise addition filter.</p>
Full article ">Figure 4
<p>Evolution of the (<b>a</b>) Information Loss and the (<b>b</b>) Disclosure Risk with the Microagreggation filter.</p>
Full article ">Figure 5
<p>Evolution of the (<b>a</b>) Information Loss and the (<b>b</b>) Disclosure Risk with the Rank Swapping filter.</p>
Full article ">Figure 6
<p>Evolution of the (<b>a</b>) Information Loss and the (<b>b</b>) Disclosure Risk with the Laplacian mechanism.</p>
Full article ">Figure 7
<p>Evolution of the (<b>a</b>) Information Loss and the (<b>b</b>) Disclosure Risk with the Exponential mechanism.</p>
Full article ">Figure 8
<p>Evolution of (<b>a</b>) Information Loss and (<b>b</b>) Disclosure Risk measures for GAN architectures.</p>
Full article ">Figure 9
<p>Evolution of Information Loss and Disclosure Risk measures for Knowledge Distillation.</p>
Full article ">
40 pages, 4144 KiB  
Article
“Exact” and Approximate Methods for Bayesian Inference: Stochastic Volatility Case Study
by Yuliya Shapovalova
Entropy 2021, 23(4), 466; https://doi.org/10.3390/e23040466 - 15 Apr 2021
Cited by 4 | Viewed by 4342
Abstract
We conduct a case study in which we empirically illustrate the performance of different classes of Bayesian inference methods to estimate stochastic volatility models. In particular, we consider how different particle filtering methods affect the variance of the estimated likelihood. We review and [...] Read more.
We conduct a case study in which we empirically illustrate the performance of different classes of Bayesian inference methods to estimate stochastic volatility models. In particular, we consider how different particle filtering methods affect the variance of the estimated likelihood. We review and compare particle Markov Chain Monte Carlo (MCMC), RMHMC, fixed-form variational Bayes, and integrated nested Laplace approximation to estimate the posterior distribution of the parameters. Additionally, we conduct the review from the point of view of whether these methods are (1) easily adaptable to different model specifications; (2) adaptable to higher dimensions of the model in a straightforward way; (3) feasible in the multivariate case. We show that when using the stochastic volatility model for methods comparison, various data-generating processes have to be considered to make a fair assessment of the methods. Finally, we present a challenging specification of the multivariate stochastic volatility model, which is rarely used to illustrate the methods but constitutes an important practical application. Full article
(This article belongs to the Special Issue Approximate Bayesian Inference)
Show Figures

Figure 1

Figure 1
<p>Graphical representation of stochastic volatility model. Observations <math display="inline"><semantics> <msub> <mi>y</mi> <mi>t</mi> </msub> </semantics></math> represented by shaded edges depend at each time point on the state of the latent volatility process <math display="inline"><semantics> <msub> <mi>h</mi> <mi>t</mi> </msub> </semantics></math>.</p>
Full article ">Figure 2
<p>Variance of the estimated likelihood in different points of the parameter space for <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>S</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a>.</p>
Full article ">Figure 3
<p>Variance of the estimated likelihood in different points of the parameter space for <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>S</mi> <mo>=</mo> <mn>3</mn> </mrow> </semantics></math> from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a>.</p>
Full article ">Figure 4
<p>Variance of the estimated likelihood in different points of the parameter space for <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>S</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math> from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a>.</p>
Full article ">Figure 5
<p>Results of the sampling from the posterior distribution with PMH and RMHMC for <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>S</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a>. The first column corresponds to the trace plots, the middle column to histograms obtained with the samples from the posterior distribution, and the last column corresponds to autocorrelation function for the samples.</p>
Full article ">Figure 6
<p>Illustration of the Fixed-form Variational Bayes in comparison to RMHMC and INLA. Subfigures illustrate the posterior distributions estimated with different methods for the different data-generating processes. (<b>a</b>–<b>c</b>) correspond to Experiment 1 from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a> and <a href="#entropy-23-00466-t006" class="html-table">Table 6</a>, (<b>d</b>–<b>f</b>) correspond to Experiment 2, (<b>g</b>–<b>i</b>) correspond to Experiment 3, (<b>j</b>–<b>l</b>) correspond to experiment 4, and (<b>m</b>–<b>o</b>) correspond to Experiment 5. Red vertical lines indicate true parameter values.</p>
Full article ">Figure 7
<p>Comparison of PMH (pink), VB (blue), INLA (green), and RMHMC (yellow) on the weekly log-returns of the Australian dollar against the US dollar. Subfigures illustrate the posterior distributions for different parameters of the model obtained with different methods. (<b>a</b>) Corresponds to the posterior distribution for the parameter <math display="inline"><semantics> <mi>β</mi> </semantics></math>. (<b>b</b>) Corresponds to the posterior distribution of the parameter <math display="inline"><semantics> <mi>ϕ</mi> </semantics></math>. (<b>c</b>) Corresponds to the posterior distribution of the parameter <math display="inline"><semantics> <msub> <mi>σ</mi> <mi>η</mi> </msub> </semantics></math>.</p>
Full article ">Figure 8
<p>Comparison of PMH (pink), VB (blue), INLA (green), and RMHMC (yellow) on the daily log-returns of DAX index. Subfigures illustrate the posterior distributions for different parameters of the model obtained with different methods. (<b>a</b>,<b>b</b>) correspond to the posterior distribution of the parameter <math display="inline"><semantics> <mi>β</mi> </semantics></math>. (<b>c</b>) corresponds to the posterior distribution of the parameter <math display="inline"><semantics> <mi>ϕ</mi> </semantics></math>. (<b>d</b>) corresponds to the posterior distribution of the parameter <math display="inline"><semantics> <msub> <mi>σ</mi> <mi>η</mi> </msub> </semantics></math>.</p>
Full article ">Figure A1
<p>Results of the sampling from the posterior distribution with PMH and RMHMC for <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>S</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math> from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a>. The first column corresponds to the trace plots, the middle column to histograms obtained with the samples from the posterior distribution, and the last column to autocorrelation function for the samples.</p>
Full article ">Figure A2
<p>Results of the sampling from the posterior distribution with PMH and RMHMC for <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>S</mi> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math> from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a>. The first column corresponds to the trace plots, the middle column to histograms obtained with the samples from the posterior distribution, and the last column to autocorrelation function for the samples.</p>
Full article ">Figure A3
<p>Results of the sampling from the posterior distribution with PMH and RMHMC for <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>S</mi> <mo>=</mo> <mn>9</mn> </mrow> </semantics></math> from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a>. The first column corresponds to the trace plots, the middle column to histograms obtained with the samples from the posterior distribution, and the last column to autocorrelation function for the samples.</p>
Full article ">Figure A4
<p>Results of the sampling from the posterior distribution with PMH and RMHMC for <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>S</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> from <a href="#entropy-23-00466-t002" class="html-table">Table 2</a>. The first column corresponds to the trace plots, the middle column to histograms obtained with the samples from the posterior distribution, and the last column to autocorrelation function for the samples.</p>
Full article ">
9 pages, 283 KiB  
Article
On the Locally Polynomial Complexity of the Projection-Gradient Method for Solving Piecewise Quadratic Optimisation Problems
by Agnieszka Prusińska, Krzysztof Szkatuła and Alexey Tret’yakov
Entropy 2021, 23(4), 465; https://doi.org/10.3390/e23040465 - 15 Apr 2021
Viewed by 1718
Abstract
This paper proposes a method for solving optimisation problems involving piecewise quadratic functions. The method provides a solution in a finite number of iterations, and the computational complexity of the proposed method is locally polynomial of the problem dimension, i.e., if the initial [...] Read more.
This paper proposes a method for solving optimisation problems involving piecewise quadratic functions. The method provides a solution in a finite number of iterations, and the computational complexity of the proposed method is locally polynomial of the problem dimension, i.e., if the initial point belongs to the sufficiently small neighbourhood of the solution set. Proposed method could be applied for solving large systems of linear inequalities. Full article
(This article belongs to the Section Complexity)
28 pages, 1106 KiB  
Article
On a Variational Definition for the Jensen-Shannon Symmetrization of Distances Based on the Information Radius
by Frank Nielsen
Entropy 2021, 23(4), 464; https://doi.org/10.3390/e23040464 - 14 Apr 2021
Cited by 18 | Viewed by 9206
Abstract
We generalize the Jensen-Shannon divergence and the Jensen-Shannon diversity index by considering a variational definition with respect to a generic mean, thereby extending the notion of Sibson’s information radius. The variational definition applies to any arbitrary distance and yields a new way to [...] Read more.
We generalize the Jensen-Shannon divergence and the Jensen-Shannon diversity index by considering a variational definition with respect to a generic mean, thereby extending the notion of Sibson’s information radius. The variational definition applies to any arbitrary distance and yields a new way to define a Jensen-Shannon symmetrization of distances. When the variational optimization is further constrained to belong to prescribed families of probability measures, we get relative Jensen-Shannon divergences and their equivalent Jensen-Shannon symmetrizations of distances that generalize the concept of information projections. Finally, we touch upon applications of these variational Jensen-Shannon divergences and diversity indices to clustering and quantization tasks of probability measures, including statistical mixtures. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Illustrating several cases of the relative Jensen-Shannon divergence based on whether <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>∈</mo> <mi mathvariant="script">R</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>q</mi> <mo>∈</mo> <mi mathvariant="script">R</mi> </mrow> </semantics></math> or not.</p>
Full article ">Figure 2
<p>Three equivalent expressions of the ordinary (skewed) Jensen-Shannon divergence which yield three different generalizations.</p>
Full article ">
27 pages, 7136 KiB  
Article
Adaptive Diagnosis for Fault Tolerant Data Fusion Based on α-Rényi Divergence Strategy for Vehicle Localization
by Khoder Makkawi, Nourdine Ait-Tmazirte, Maan El Badaoui El Najjar and Nazih Moubayed
Entropy 2021, 23(4), 463; https://doi.org/10.3390/e23040463 - 14 Apr 2021
Cited by 9 | Viewed by 2363
Abstract
When applying a diagnostic technique to complex systems, whose dynamics, constraints, and environment evolve over time, being able to re-evaluate the residuals that are capable of detecting defaults and proposing the most appropriate ones can quickly prove to make sense. For this purpose, [...] Read more.
When applying a diagnostic technique to complex systems, whose dynamics, constraints, and environment evolve over time, being able to re-evaluate the residuals that are capable of detecting defaults and proposing the most appropriate ones can quickly prove to make sense. For this purpose, the concept of adaptive diagnosis is introduced. In this work, the contributions of information theory are investigated in order to propose a Fault-Tolerant multi-sensor data fusion framework. This work is part of studies proposing an architecture combining a stochastic filter for state estimation with a diagnostic layer with the aim of proposing a safe and accurate state estimation from potentially inconsistent or erroneous sensors measurements. From the design of the residuals, using α-Rényi Divergence (α-RD), to the optimization of the decision threshold, through the establishment of a function that is dedicated to the choice of α at each moment, we detail each step of the proposed automated decision-support framework. We also dwell on: (1) the consequences of the degree of freedom provided by this α parameter and on (2) the application-dictated policy to design the α tuning function playing on the overall performance of the system (detection rate, false alarms, and missed detection rates). Finally, we present a real application case on which this framework has been tested. The problem of multi-sensor localization, integrating sensors whose operating range is variable according to the environment crossed, is a case study to illustrate the contributions of such an approach and show the performance. Full article
Show Figures

Figure 1

Figure 1
<p>Detailed diagram for the proposed approach.</p>
Full article ">Figure 2
<p>Diagram to create two probability distributions.</p>
Full article ">Figure 3
<p>Extracting estimated means and variances from real values.</p>
Full article ">Figure 4
<p>Trajectory C3 reference.</p>
Full article ">Figure 5
<p>Positioning estimation without Fault Detection and Isolation (FDI) algorithm vs reference in three/two-dimensional (3D/2D) views.</p>
Full article ">Figure 6
<p>Elevation of each satellites during the whole trajectory.</p>
Full article ">Figure 7
<p><math display="inline"><semantics> <mi>α</mi> </semantics></math>-Rényi Residuals without FDI.</p>
Full article ">Figure 8
<p>Partials <math display="inline"><semantics> <mrow> <mi>R</mi> <msub> <mi>D</mi> <mrow> <mn>0.5</mn> </mrow> </msub> </mrow> </semantics></math> observation for identification.</p>
Full article ">Figure 9
<p><math display="inline"><semantics> <mrow> <mi>R</mi> <msub> <mi>D</mi> <mrow> <mn>0.5</mn> </mrow> </msub> </mrow> </semantics></math> divergence without/with FDI with adaptive threshold.</p>
Full article ">Figure 10
<p>P0 during the whole trajectory.</p>
Full article ">Figure 11
<p>Partials Rényi observation for identification.</p>
Full article ">Figure 12
<p>Weighted Alpha balance.</p>
Full article ">Figure 13
<p>Estimated and real means &amp; variances for faulty and non-faulty distributions.</p>
Full article ">Figure 14
<p>Rényi divergence without/with FDI with adaptive threshold.</p>
Full article ">Figure 15
<p>P0 during the whole trajectory.</p>
Full article ">Figure 16
<p>index of isolated satellites using fixed and balanced <math display="inline"><semantics> <mi>α</mi> </semantics></math>.</p>
Full article ">Figure 17
<p>The difference in decisions between fixed <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>α</mi> </semantics></math> balance.</p>
Full article ">Figure 18
<p>The effect of the decisions between fixed <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>α</mi> </semantics></math> balance on the position.</p>
Full article ">Figure 19
<p>The difference in decisions between fixed <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>α</mi> </semantics></math> balance.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop