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

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

Axioms, Volume 12, Issue 4 (April 2023) – 90 articles

Cover Story (view full-size image): Piecewise differential systems have been extensively studied in recent decades, mainly due to their applications. Of particular interest is the generalization of Hilbert's problem 16 to piecewise polynomial systems of degree n. Here, we are interested in determining the maximum number of algebraic limit cycles of the piecewise linear differential systems separated by one straight line. First, we classified all the linear differential systems with a polynomial first integral. Using this information, we have proved that if in one of the two pieces of the piecewise linear differential systems there is a Hamiltonian system of degree one, the piecewise differential system has at most one limit cycle. 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:
23 pages, 382 KiB  
Article
Fractional Step Scheme to Approximate a Non-Linear Second-Order Reaction–Diffusion Problem with Inhomogeneous Dynamic Boundary Conditions
by Constantin Fetecău and Costică Moroşanu
Axioms 2023, 12(4), 406; https://doi.org/10.3390/axioms12040406 - 21 Apr 2023
Cited by 3 | Viewed by 1782
Abstract
Two main topics are addressed in the present paper, first, a rigorous qualitative study of a second-order reaction–diffusion problem with non-linear diffusion and cubic-type reactions, as well as inhomogeneous dynamic boundary conditions. Under certain assumptions about the input data: [...] Read more.
Two main topics are addressed in the present paper, first, a rigorous qualitative study of a second-order reaction–diffusion problem with non-linear diffusion and cubic-type reactions, as well as inhomogeneous dynamic boundary conditions. Under certain assumptions about the input data: gd(t,x), gfr(t,x), U0(x) and ζ0(x), we prove the well-posedness (the existence, a priori estimates, regularity and uniqueness) of a solution in the space Wp1,2(Q)×Wp1,2(Σ). Here, we extend previous results, enabling new mathematical models to be more suitable to describe the complexity of a wide class of different physical phenomena of life sciences, including moving interface problems, material sciences, digital image processing, automatic vehicle detection and tracking, the spread of an epidemic infection, semantic image segmentation including U-Net neural networks, etc. The second goal is to develop an iterative splitting scheme, corresponding to the non-linear second-order reaction–diffusion problem. Results relating to the convergence of the approximation scheme and error estimation are also established. On the basis of the proposed numerical scheme, we formulate the algorithm alg-frac_sec-ord_dbc, which represents a delicate challenge for our future works. The benefit of such a method could simplify the process of numerical computation. Full article
(This article belongs to the Special Issue Recent Advances in Fractional Calculus)
14 pages, 331 KiB  
Article
The Cauchy–Optimal Stability Results for Cauchy–Jensen Additive Mappings in the Fuzzy Banach Space and the Unital Fuzzy Banach Space
by Zahra Eidinejad, Reza Saadati and Hari M. Srivastava
Axioms 2023, 12(4), 405; https://doi.org/10.3390/axioms12040405 - 21 Apr 2023
Viewed by 1160
Abstract
In this article, we apply a new class of fuzzy control functions to approximate a Cauchy additive mapping in fuzzy Banach space (FBS). Further, considering the unital FBS (UFBS), we will investigate the isomorphisms defined in this space. By introducing several specific functions [...] Read more.
In this article, we apply a new class of fuzzy control functions to approximate a Cauchy additive mapping in fuzzy Banach space (FBS). Further, considering the unital FBS (UFBS), we will investigate the isomorphisms defined in this space. By introducing several specific functions and choosing the optimal control function from among these functions, we evaluate the Cauchy–Optimal stability (C–O-stability) for all defined mappings. Full article
(This article belongs to the Special Issue Mathematical Analysis and Applications IV)
19 pages, 1057 KiB  
Article
Multi-Objective Optimization of the Robustness of Complex Networks Based on the Mixture of Weighted Surrogates
by Junfeng Nie, Zhuoran Yu and Junli Li
Axioms 2023, 12(4), 404; https://doi.org/10.3390/axioms12040404 - 21 Apr 2023
Viewed by 1845
Abstract
Network robustness is of paramount importance. Although great progress has been achieved in robustness optimization using single measures, such networks may still be vulnerable to many attack scenarios. Consequently, multi-objective network robustness optimization has recently garnered greater attention. A complex network structure plays [...] Read more.
Network robustness is of paramount importance. Although great progress has been achieved in robustness optimization using single measures, such networks may still be vulnerable to many attack scenarios. Consequently, multi-objective network robustness optimization has recently garnered greater attention. A complex network structure plays an important role in both node-based and link-based attacks. In this paper, since multi-objective robustness optimization comes with a high computational cost, a surrogate model is adopted instead of network controllability robustness in the optimization process, and the Dempster–Shafer theory is used for selecting and mixing the surrogate models. The method has been validated on four types of synthetic networks, and the results show that the two selected surrogate models can effectively assist the multi-objective evolutionary algorithm in finding network structures with improved controllability robustness. The adaptive updating of surrogate models during the optimization process leads to better results than the selection of two surrogate models, albeit at the cost of longer processing times. Furthermore, the method demonstrated in this paper achieved better performance than existing methods, resulting in a marked increase in computational efficiency. Full article
(This article belongs to the Special Issue Complex Networks, Evolutionary Computation and Machine Learning)
Show Figures

Figure 1

Figure 1
<p>The structure of the surrogate ensemble for estimating network robustness.</p>
Full article ">Figure 2
<p>D-S theory weighted evaluation of surrogate model steps.</p>
Full article ">Figure 3
<p>Non-dominated solutions of different methods for four types of complex networks.</p>
Full article ">Figure 4
<p>HV values of different methods under four types of complex networks.</p>
Full article ">Figure 5
<p><math display="inline"><semantics> <mrow> <mi>MOEA</mi> <mo>_</mo> <mi>Two</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>MOEA</mi> <mo>_</mo> <mi>Two</mi> <mo>_</mo> <mi>Adapt</mi> </mrow> </semantics></math> under four types of complex network.</p>
Full article ">Figure 6
<p>HV values of <math display="inline"><semantics> <mrow> <mi>MOEA</mi> <mo>_</mo> <mi>Two</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>MOEA</mi> <mo>_</mo> <mi>Two</mi> <mo>_</mo> <mi>Adapt</mi> </mrow> </semantics></math> in four types of network.</p>
Full article ">
17 pages, 900 KiB  
Article
Aperiodic Sampled-Data Control for Anti-Synchronization of Chaotic Nonlinear Systems Subject to Input Saturation
by Meixuan Li and Yingjie Fan
Axioms 2023, 12(4), 403; https://doi.org/10.3390/axioms12040403 - 21 Apr 2023
Cited by 2 | Viewed by 1549
Abstract
This paper studies the aperiodic sampled-data (SD) control anti-synchronization issue of chaotic nonlinear systems under the effects of input saturation. At first, to describe the simultaneous existence of the aperiodic SD pattern and the input saturation, a nonlinear closed-loop system model is established. [...] Read more.
This paper studies the aperiodic sampled-data (SD) control anti-synchronization issue of chaotic nonlinear systems under the effects of input saturation. At first, to describe the simultaneous existence of the aperiodic SD pattern and the input saturation, a nonlinear closed-loop system model is established. Then, to make the anti-synchronization analysis, a relaxed sampling-interval-dependent Lyapunov functional (RSIDLF) is constructed for the resulting closed-loop system. Thereinto, the positive definiteness requirement of the RSIDLF is abandoned. Due to the indefiniteness of RSIDLF, the discrete-time Lyapunov method (DTLM) then is used to guarantee the local stability of the trivial solutions of the modeled nonlinear system. Furthermore, two convex optimization schemes are proposed to expand the allowable initial area (AIA) and maximize the upper bound of the sampling period (UBSP). Finally, two examples of nonlinear systems are provided to illustrate the superiority of the RSIDLF method over the previous methods in expanding the AIA and enlarging the UBSP. Full article
(This article belongs to the Topic Advances in Nonlinear Dynamics: Methods and Applications)
Show Figures

Figure 1

Figure 1
<p>Chaotic behavior of <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> on (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>α</mi> <mn>1</mn> </msub> <mo>−</mo> <msub> <mi>α</mi> <mn>2</mn> </msub> </mrow> </semantics></math> and (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>α</mi> <mn>2</mn> </msub> <mo>−</mo> <msub> <mi>α</mi> <mn>3</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>The allowable initial area (AIAs) based on Theorem 1 and Corollary 2 (method in [<a href="#B22-axioms-12-00403" class="html-bibr">22</a>,<a href="#B35-axioms-12-00403" class="html-bibr">35</a>]).</p>
Full article ">Figure 3
<p>Schematic of the saturated aperiodic sampled-data (SD) controller <math display="inline"><semantics> <mrow> <mi>sat</mi> <mo>(</mo> <mi mathvariant="normal">u</mi> <mo>(</mo> <mi mathvariant="normal">t</mi> <mo>)</mo> <mo>)</mo> </mrow> </semantics></math> of Example I.</p>
Full article ">Figure 4
<p>Schematic of the master system <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> and the slave system <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> under the saturated aperiodic SD controller of Example I.</p>
Full article ">Figure 5
<p>Schematic of the error system <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> under the saturated aperiodic SD controller of Example I.</p>
Full article ">Figure 6
<p>Schematic of the saturated aperiodic SD controller <math display="inline"><semantics> <mrow> <mi>sat</mi> <mo>(</mo> <mi mathvariant="normal">u</mi> <mo>(</mo> <mi mathvariant="normal">t</mi> <mo>)</mo> <mo>)</mo> </mrow> </semantics></math> of Example II.</p>
Full article ">Figure 7
<p>Schematic of the master system <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> and the slave system <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> under the saturated aperiodic SD controller of Example II.</p>
Full article ">Figure 8
<p>Schematic of the error system <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> under the saturated aperiodic SD controller of Example II.</p>
Full article ">
22 pages, 3141 KiB  
Article
Default Prediction with Industry-Specific Default Heterogeneity Indicators Based on the Forward Intensity Model
by Zhengfang Ni, Minghui Jiang and Wentao Zhan
Axioms 2023, 12(4), 402; https://doi.org/10.3390/axioms12040402 - 21 Apr 2023
Viewed by 1612
Abstract
When predicting the defaults of a large number of samples in a region, this will be affected by industry default heterogeneity. To build a credit risk model that is more suitable for Chinese-listed firms, which have highly industry-specific default heterogeneity, we extend the [...] Read more.
When predicting the defaults of a large number of samples in a region, this will be affected by industry default heterogeneity. To build a credit risk model that is more suitable for Chinese-listed firms, which have highly industry-specific default heterogeneity, we extend the forward intensity model to predict the defaults of Chinese-listed firms with information about the default heterogeneity of industries. Compared with the original model, we combine the Bayes approach with the forward intensity model to generate time-varying industry-specific default heterogeneity indicators. Our model can capture co-movements of different industries that cannot be observed based on the original forward intensity model so that the model can flexibly adjust the firm’s PD according to the industry. In addition, we also consider the impact of default heterogeneity in other industries by studying the influence of the level and trends of other industries’ default heterogeneity on a firm’s credit risk. Finally, we compute PDs for 4476 firms from January 2001 to December 2019 for 36 prediction horizons. The extended model improves the prediction accuracy ratios both for the in-sample and out-of-sample firm’s PDs for all 36 horizons. Almost all the accuracy ratios of the prediction horizons’ PDs are increased by more than 6%. In addition, our model also reduces the gap between the aggregated PDs and the realized number of defaults. Our industry-specific default heterogeneity indicator is helpful to improve the model’s performance, especially for predicting defaults in a large portfolio, which is of significance for credit risk management in China and other regions. Full article
(This article belongs to the Special Issue Big Data Analytics and Mathematical Methods in Digital Economy)
Show Figures

Figure 1

Figure 1
<p>The number of all industry defaults per year. Source: NUS-CRI database.</p>
Full article ">Figure 2
<p>The frequency of all industry defaults per year. Source: NUS-CRI database.</p>
Full article ">Figure 3
<p>This figure shows a comparison of the realized number of defaults, original aggregated accumulated PDs, and new aggregated accumulated PDs for a 12-month prediction horizon in the sample.</p>
Full article ">Figure 4
<p>This figure shows a comparison of the realized number of defaults, original aggregated accumulated PDs, and new aggregated accumulated PDs for a 12-month prediction horizon out of sample.</p>
Full article ">Figure 5
<p>This figure shows the in-sample cumulative accuracy profiles of new PDs from January 2000 to December 2019 for different prediction horizons.</p>
Full article ">Figure 6
<p>This figure shows the out-of-sample cumulative accuracy profiles of new PDs from January 2000 to December 2019 for different prediction horizons.</p>
Full article ">Figure 7
<p>This figure compares the in-sample cumulative accuracy profiles of new PDs and original PDs for horizons of 1–36 months.</p>
Full article ">Figure 8
<p>This figure compares the out-of-sample cumulative accuracy profiles of new PDs and original PDs for horizons of 1–36 months.</p>
Full article ">
14 pages, 313 KiB  
Article
On the Fixed Circle Problem on Metric Spaces and Related Results
by Nabil Mlaiki, Nihal Özgür, Nihal Taş and Dania Santina
Axioms 2023, 12(4), 401; https://doi.org/10.3390/axioms12040401 - 20 Apr 2023
Cited by 3 | Viewed by 1652
Abstract
The fixed-circle issue is a geometric technique that is connected to the study of geometric characteristics of certain points, and that are fixed by the self-mapping of either the metric space or of the generalized space. The fixed-disc problem is a natural resultant [...] Read more.
The fixed-circle issue is a geometric technique that is connected to the study of geometric characteristics of certain points, and that are fixed by the self-mapping of either the metric space or of the generalized space. The fixed-disc problem is a natural resultant that arises as a direct outcome of this problem. In this study, our goal is to examine new classes of self-mappings that meet a new particular sort of contraction in a metric space. The common geometrical characteristic of the set of fixed points of any element of these classes is that a circle or even a disc, that is either termed the fixed circle or even the fixed disc of the appropriate self-map, is included within that set. In order to accomplish this, we establish two new classifications of contraction mapping: Fc-contractive mapping and Fc-expanding mapping. In the investigation of neural networks, activation functions with either fixed circles (or even fixed discs) are observed frequently. This demonstrates how successful our results with the fixed-circle (respectively, the fixed-disc) model were. Full article
(This article belongs to the Special Issue Fixed Point Theory and Its Related Topics IV)
Show Figures

Figure 1

Figure 1
<p>The graph of the circle <math display="inline"><semantics> <msub> <mi>C</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> </semantics></math>.</p>
Full article ">
13 pages, 1294 KiB  
Article
Two Novel Computational Techniques for Solving Nonlinear Time-Fractional Lax’s Korteweg-de Vries Equation
by Nidhish Kumar Mishra, Mashael M. AlBaidani, Adnan Khan and Abdul Hamid Ganie
Axioms 2023, 12(4), 400; https://doi.org/10.3390/axioms12040400 - 20 Apr 2023
Cited by 13 | Viewed by 1411
Abstract
This article investigates the seventh-order Lax’s Korteweg–de Vries equation using the Yang transform decomposition method (YTDM) and the homotopy perturbation transform method (HPTM). The physical phenomena that emerge in physics, engineering and chemistry are mathematically expressed by this equation. For instance, the KdV [...] Read more.
This article investigates the seventh-order Lax’s Korteweg–de Vries equation using the Yang transform decomposition method (YTDM) and the homotopy perturbation transform method (HPTM). The physical phenomena that emerge in physics, engineering and chemistry are mathematically expressed by this equation. For instance, the KdV equation was constructed to represent a wide range of physical processes involving the evolution and interaction of nonlinear waves. In the Caputo sense, the fractional derivative is considered. We employed the Yang transform, the Adomian decomposition method and the homotopy perturbation method to obtain the solution to the time-fractional Lax’s Korteweg–de Vries problem. We examined and compared a particular example with the actual result to verify the approaches. By utilizing these methods, we can construct recurrence relations that represent the solution to the problem that is being proposed, and we are then able to present graphical representations that enable us to visually examine all of the results in the proposed case for different fractional order values. Furthermore, the results of the current approach exhibit a good correlation with the precise solution to the problem being studied. Furthermore, the present study offers an example of error analysis. The numerical outcomes obtained by applying the provided approaches demonstrate that the techniques are easy to use and have superior computational performance. Full article
(This article belongs to the Special Issue Fractional Calculus and the Applied Analysis)
Show Figures

Figure 1

Figure 1
<p>Graphical depiction of our techniques and the accurate solution.</p>
Full article ">Figure 2
<p>Graphical depiction of our techniques solution at <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>=</mo> <mn>0.8</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mn>0.6</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Graphical depiction of our techniques solution for various values of <math display="inline"><semantics> <mi>σ</mi> </semantics></math>.</p>
Full article ">Figure 4
<p>Graphical depiction our techniques solution in terms of error.</p>
Full article ">
19 pages, 928 KiB  
Article
Properties of Convex Fuzzy-Number-Valued Functions on Harmonic Convex Set in the Second Sense and Related Inequalities via Up and Down Fuzzy Relation
by Muhammad Bilal Khan, Željko Stević, Abdulwadoud A. Maash, Muhammad Aslam Noor and Mohamed S. Soliman
Axioms 2023, 12(4), 399; https://doi.org/10.3390/axioms12040399 - 20 Apr 2023
Cited by 2 | Viewed by 1187
Abstract
In this paper, we provide different variants of the Hermite–Hadamard (HH) inequality using the concept of a new class of convex mappings, which is referred to as up and down harmonically s-convex fuzzy-number-valued functions (UDH [...] Read more.
In this paper, we provide different variants of the Hermite–Hadamard (HH) inequality using the concept of a new class of convex mappings, which is referred to as up and down harmonically s-convex fuzzy-number-valued functions (UDH s-convex FNVM) in the second sense based on the up and down fuzzy inclusion relation. The findings are confirmed with certain numerical calculations that take a few appropriate examples into account. The results deal with various integrals of the 2ρσρ+σ type and are innovative in the setting of up and down harmonically s-convex fuzzy-number-valued functions. Moreover, we acquire classical and new exceptional cases that can be seen as applications of our main outcomes. In our opinion, this will make a significant contribution to encouraging more research. Full article
20 pages, 8061 KiB  
Article
Reconstructing Loads in Nanoplates from Dynamic Data
by Alexandre Kawano and Antonino Morassi
Axioms 2023, 12(4), 398; https://doi.org/10.3390/axioms12040398 - 20 Apr 2023
Viewed by 1198
Abstract
It was recently proved that the knowledge of the transverse displacement of a nanoplate in an open subset of its mid-plane, measured for any interval of time, allows for the unique determination of the spatial components [...] Read more.
It was recently proved that the knowledge of the transverse displacement of a nanoplate in an open subset of its mid-plane, measured for any interval of time, allows for the unique determination of the spatial components {fm(x,y)}m=1M of the transverse load m=1Mgm(t)fm(x,y), where M1 and {gm(t)}m=1M is a known set of linearly independent functions of the time variable. The nanoplate mechanical model is built within the strain gradient linear elasticity theory, according to the Kirchhoff–Love kinematic assumptions. In this paper, we derive a reconstruction algorithm for the above inverse source problem, and we implement a numerical procedure based on a finite element spatial discretization to approximate the loads {fm(x,y)}m=1M. The computations are developed for a uniform rectangular nanoplate clamped at the boundary. The sensitivity of the results with respect to the main parameters that influence the identification is analyzed in detail. The adoption of a regularization scheme based on the singular value decomposition turns out to be decisive for the accuracy and stability of the reconstruction. Full article
Show Figures

Figure 1

Figure 1
<p>First six eigenfunctions obtained by the FEM with a <math display="inline"><semantics> <mrow> <mn>6</mn> <mo>×</mo> <mn>5</mn> </mrow> </semantics></math> mesh.</p>
Full article ">Figure 2
<p>Comparison of the numerical time histories of the displacement function at the point <math display="inline"><semantics> <mrow> <mo stretchy="false">(</mo> <mfrac> <msub> <mi>L</mi> <mi>X</mi> </msub> <mn>2</mn> </mfrac> <mo>,</mo> <mfrac> <msub> <mi>L</mi> <mi>Y</mi> </msub> <mn>3</mn> </mfrac> <mo stretchy="false">)</mo> </mrow> </semantics></math>. (<b>a</b>) Comparison between the results for <math display="inline"><semantics> <mrow> <mn>4</mn> <mo>×</mo> <mn>3</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mn>6</mn> <mo>×</mo> <mn>5</mn> </mrow> </semantics></math> meshes. (<b>b</b>) Comparison between the results for <math display="inline"><semantics> <mrow> <mn>4</mn> <mo>×</mo> <mn>3</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mn>6</mn> <mo>×</mo> <mn>5</mn> </mrow> </semantics></math> meshes with the exact solution.</p>
Full article ">Figure 3
<p>Spatial loads. (<b>a</b>) Load <math display="inline"><semantics> <msub> <mi>f</mi> <mn>1</mn> </msub> </semantics></math>. (<b>b</b>) Load <math display="inline"><semantics> <msub> <mi>f</mi> <mn>2</mn> </msub> </semantics></math>.</p>
Full article ">Figure 4
<p>Target loads. (<b>a</b>) Load <math display="inline"><semantics> <msub> <mi>f</mi> <mn>1</mn> </msub> </semantics></math>. (<b>b</b>) Load <math display="inline"><semantics> <msub> <mi>f</mi> <mn>2</mn> </msub> </semantics></math>.</p>
Full article ">Figure 5
<p>Comparison between the exact solution (solid line) and the approximate solution obtained by the FEM. Displacement history at the point <math display="inline"><semantics> <mrow> <mo stretchy="false">(</mo> <mfrac> <msub> <mi>L</mi> <mi>X</mi> </msub> <mn>2</mn> </mfrac> <mo>,</mo> <mfrac> <msub> <mi>L</mi> <mi>Y</mi> </msub> <mn>3</mn> </mfrac> <mo stretchy="false">)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Comparison of the reconstruction results when we vary the size of the observation set <math display="inline"><semantics> <mover accent="true"> <mo>Ω</mo> <mo stretchy="false">˜</mo> </mover> </semantics></math>. (<b>a</b>) <math display="inline"><semantics> <mrow> <mover accent="true"> <mo>Ω</mo> <mo stretchy="false">˜</mo> </mover> <mo>=</mo> <mrow> <mo>]</mo> <msub> <mi>x</mi> <mi>c</mi> </msub> <mo>−</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>4</mn> <mo>,</mo> <msub> <mi>x</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>4</mn> <mrow> <mo>[</mo> <mo>×</mo> <mo>]</mo> </mrow> <msub> <mi>y</mi> <mi>c</mi> </msub> <mo>−</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>4</mn> <mo>,</mo> <msub> <mi>y</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>4</mn> <mo>[</mo> </mrow> </mrow> </semantics></math>. (<b>b</b>) <math display="inline"><semantics> <mrow> <mover accent="true"> <mo>Ω</mo> <mo stretchy="false">˜</mo> </mover> <mo>=</mo> <mrow> <mo>]</mo> <msub> <mi>x</mi> <mi>c</mi> </msub> <mo>−</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>8</mn> <mo>,</mo> <msub> <mi>x</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>8</mn> <mrow> <mo>[</mo> <mo>×</mo> <mo>]</mo> </mrow> <msub> <mi>y</mi> <mi>c</mi> </msub> <mo>−</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>8</mn> <mo>,</mo> <msub> <mi>y</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>8</mn> <mo>[</mo> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>Comparison of the reconstruction results when we vary the position of the observation set <math display="inline"><semantics> <mover accent="true"> <mo>Ω</mo> <mo stretchy="false">˜</mo> </mover> </semantics></math>. (<b>a</b>) <math display="inline"><semantics> <mrow> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>3</mn> <mo>,</mo> <msub> <mi>L</mi> <mi>Y</mi> </msub> <mo>/</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </semantics></math>. (<b>b</b>) <math display="inline"><semantics> <mrow> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mi>X</mi> </msub> <mo>/</mo> <mn>4</mn> <mo>,</mo> <msub> <mi>L</mi> <mi>Y</mi> </msub> <mo>/</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 8
<p>Comparison of the results when we vary the position of the observation set <math display="inline"><semantics> <mover accent="true"> <mo>Ω</mo> <mo stretchy="false">˜</mo> </mover> </semantics></math>. (<b>a</b>) Attainable target load with two peaks. (<b>b</b>) Recovered load.</p>
Full article ">Figure 9
<p>Recovery of a discontinuous load. (<b>a</b>) Target load. (<b>b</b>) Attainable target load. (<b>c</b>) Recovered load (error level <math display="inline"><semantics> <mrow> <msub> <mi>ϵ</mi> <mi>e</mi> </msub> <mo>=</mo> <mn>0.05</mn> </mrow> </semantics></math>).</p>
Full article ">Figure 10
<p>Recovery of a discontinuous load. (<b>a</b>) Sections of the recovered load on the line <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <msub> <mi>L</mi> <mi>Y</mi> </msub> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math>. Error level <math display="inline"><semantics> <mrow> <msub> <mi>ϵ</mi> <mi>e</mi> </msub> <mo>=</mo> <mn>0.05</mn> </mrow> </semantics></math>. (<b>b</b>) Sections of the recovered load on the line <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <msub> <mi>L</mi> <mi>Y</mi> </msub> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math>. Error level <math display="inline"><semantics> <mrow> <msub> <mi>ϵ</mi> <mi>e</mi> </msub> <mo>=</mo> <mn>0.02</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 11
<p>Application of the SVD-based regularization strategy. (<b>a</b>) Recovery of the load corresponding to <a href="#axioms-12-00398-f007" class="html-fig">Figure 7</a>a when the regularization strategy is applied with <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mrow> <mi>T</mi> <mi>r</mi> <mi>u</mi> <mi>n</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>20</mn> </mrow> </semantics></math>. (<b>b</b>) Recovery of the load corresponding to <a href="#axioms-12-00398-f007" class="html-fig">Figure 7</a>b when the regularization strategy is applied with <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mrow> <mi>T</mi> <mi>r</mi> <mi>u</mi> <mi>n</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>20</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 12
<p>Recovery of loads when there is no regularization. (<b>a</b>) Recovery of <math display="inline"><semantics> <msub> <mi>f</mi> <mn>1</mn> </msub> </semantics></math> when there is no regularization. (<b>b</b>) Recovery of <math display="inline"><semantics> <msub> <mi>f</mi> <mn>2</mn> </msub> </semantics></math> when there is no regularization.</p>
Full article ">Figure 13
<p>Recovery of loads when the SVD-based regularization is applied, with <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mrow> <mi>T</mi> <mi>r</mi> <mi>u</mi> <mi>n</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>12</mn> </mrow> </semantics></math>. (<b>a</b>) Recovery of <math display="inline"><semantics> <msub> <mi>f</mi> <mn>1</mn> </msub> </semantics></math> with SVD-based regularization. <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mrow> <mi>T</mi> <mi>r</mi> <mi>u</mi> <mi>n</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>12</mn> </mrow> </semantics></math>. (<b>b</b>) Recovery of <math display="inline"><semantics> <msub> <mi>f</mi> <mn>2</mn> </msub> </semantics></math> with SVD-based regularization. <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mrow> <mi>T</mi> <mi>r</mi> <mi>u</mi> <mi>n</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>12</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 14
<p>Influence of the thickness. (<b>a</b>) Recovered load. Plate thickness <math display="inline"><semantics> <msub> <mi>h</mi> <mn>0</mn> </msub> </semantics></math> = 5 <math display="inline"><semantics> <msub> <mi>l</mi> <mn>0</mn> </msub> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.091</mn> </mrow> </semantics></math>. (<b>b</b>) Recovered load. Plate thickness <math display="inline"><semantics> <mrow> <msub> <mi>h</mi> <mn>0</mn> </msub> <mo>=</mo> </mrow> </semantics></math>7 <math display="inline"><semantics> <msub> <mi>l</mi> <mn>0</mn> </msub> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.11</mn> </mrow> </semantics></math>. (<b>c</b>) Recovered load. Plate thickness <math display="inline"><semantics> <mrow> <msub> <mi>h</mi> <mn>0</mn> </msub> <mo>=</mo> </mrow> </semantics></math>2 <math display="inline"><semantics> <msub> <mi>l</mi> <mn>0</mn> </msub> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.05</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 15
<p>Influence of the density. (<b>a</b>) Recovered load. Mass per unit area <math display="inline"><semantics> <mrow> <msub> <mi>ρ</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>4.444</mn> <mo>×</mo> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> <mi>kg</mi> <mo>/</mo> <msup> <mi mathvariant="normal">m</mi> <mn>2</mn> </msup> </mrow> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.09</mn> </mrow> </semantics></math>. (<b>b</b>) Recovered load. Mass per unit area 0.2 <math display="inline"><semantics> <msub> <mi>ρ</mi> <mn>0</mn> </msub> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.013</mn> </mrow> </semantics></math>. (<b>c</b>) Recovered load. Mass per unit area 0.5 <math display="inline"><semantics> <msub> <mi>ρ</mi> <mn>0</mn> </msub> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.04</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 16
<p>Influence of the aspect ratio. (<b>a</b>) Recovered load. Aspect ratio <math display="inline"><semantics> <mrow> <mfrac> <msub> <mi>L</mi> <mi>Y</mi> </msub> <msub> <mi>L</mi> <mi>X</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <mn>3</mn> <mn>4</mn> </mfrac> </mrow> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.09</mn> </mrow> </semantics></math>. (<b>b</b>) Recovered load. Aspect ratio <math display="inline"><semantics> <mrow> <mfrac> <msub> <mi>L</mi> <mi>Y</mi> </msub> <msub> <mi>L</mi> <mi>X</mi> </msub> </mfrac> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.04</mn> </mrow> </semantics></math>. (<b>c</b>) Recovered load. Aspect ratio <math display="inline"><semantics> <mrow> <mfrac> <msub> <mi>L</mi> <mi>Y</mi> </msub> <msub> <mi>L</mi> <mi>X</mi> </msub> </mfrac> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>. Relative error <math display="inline"><semantics> <mrow> <msub> <mi>Error</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>2.06</mn> </mrow> </semantics></math>.</p>
Full article ">
16 pages, 365 KiB  
Article
On the Existence of Solutions to a Boundary Value Problem via New Weakly Contractive Operator
by Rhoda Chiroma, Trad Alotaibi, Mohammed Shehu Shagari, Awad A. Bakery, OM Kalthum S. K. Mohammed and Arafa O. Mustafa
Axioms 2023, 12(4), 397; https://doi.org/10.3390/axioms12040397 - 20 Apr 2023
Cited by 1 | Viewed by 1242
Abstract
In this paper, the notion of generalized quasi-weakly contractive operators in metric-like spaces is introduced, and new conditions for the existence of fixed points for such mappings are investigated. A non-trivial example which highlights the novelty of our principal idea is constructed. It [...] Read more.
In this paper, the notion of generalized quasi-weakly contractive operators in metric-like spaces is introduced, and new conditions for the existence of fixed points for such mappings are investigated. A non-trivial example which highlights the novelty of our principal idea is constructed. It is observed comparatively that the proposed concepts herein subsume some important results in the corresponding literature. As an application, one of our obtained findings is utilized to setup novel criteria for the existence of solutions to two-point boundary value problems of a second order differential equation. To attract new researchers in the directions examined in this article, a significant number of corollaries are pointed out and discussed. Full article
(This article belongs to the Special Issue Fixed Point Theory and Its Related Topics IV)
Show Figures

Figure 1

Figure 1
<p>Illustration of contractive inequality (<a href="#FD7-axioms-12-00397" class="html-disp-formula">7</a>) using Example 2.</p>
Full article ">Figure 2
<p>Illustration of contractive inequality (<a href="#FD24-axioms-12-00397" class="html-disp-formula">24</a>) using Example 3.</p>
Full article ">
12 pages, 289 KiB  
Article
Fixed-Point Theorems for Nonlinear Contraction in Fuzzy-Controlled Bipolar Metric Spaces
by Gunaseelan Mani, Arul Joseph Gnanaprakasam, Santosh Kumar, Ozgur Ege and Manuel De la Sen
Axioms 2023, 12(4), 396; https://doi.org/10.3390/axioms12040396 - 19 Apr 2023
Cited by 4 | Viewed by 1144
Abstract
In this paper, we introduce the concept of fuzzy-controlled bipolar metric space and prove some fixed-point theorems in this space. Our results generalize and expand some of the literature’s well-known results. We also provide some applications of our main results to integral equations. [...] Read more.
In this paper, we introduce the concept of fuzzy-controlled bipolar metric space and prove some fixed-point theorems in this space. Our results generalize and expand some of the literature’s well-known results. We also provide some applications of our main results to integral equations. Full article
(This article belongs to the Special Issue Nonlinear Functional Analysis in Natural Sciences)
19 pages, 1079 KiB  
Article
Multi-Objective ABC-NM Algorithm for Multi-Dimensional Combinatorial Optimization Problem
by Muniyan Rajeswari, Rajakumar Ramalingam, Shakila Basheer, Keerthi Samhitha Babu, Mamoon Rashid and Ramar Saranya
Axioms 2023, 12(4), 395; https://doi.org/10.3390/axioms12040395 - 19 Apr 2023
Cited by 2 | Viewed by 1326
Abstract
This article addresses the problem of converting a single-objective combinatorial problem into a multi-objective one using the Pareto front approach. Although existing algorithms can identify the optimal solution in a multi-objective space, they fail to satisfy constraints while achieving optimal performance. To address [...] Read more.
This article addresses the problem of converting a single-objective combinatorial problem into a multi-objective one using the Pareto front approach. Although existing algorithms can identify the optimal solution in a multi-objective space, they fail to satisfy constraints while achieving optimal performance. To address this issue, we propose a multi-objective artificial bee colony optimization algorithm with a classical multi-objective theme called fitness sharing. This approach helps the convergence of the Pareto solution set towards a single optimal solution that satisfies multiple objectives. This article introduces multi-objective optimization with an example of a non-dominated sequencing technique and fitness sharing approach. The experimentation is carried out in MATLAB 2018a. In addition, we applied the proposed algorithm to two different real-time datasets, namely the knapsack problem and the nurse scheduling problem (NSP). The outcome of the proposed MBABC-NM algorithm is evaluated using standard performance indicators such as average distance, number of reference solutions (NRS), overall count of attained solutions (TNS), and overall non-dominated generation volume (ONGV). The results show that it outperforms other algorithms. Full article
Show Figures

Figure 1

Figure 1
<p>Workflow of MBABC-NM.</p>
Full article ">Figure 2
<p>Performance of MBABC-NM w.r.t NRS.</p>
Full article ">Figure 3
<p>Performance of MBABC-NM with respect to TNS.</p>
Full article ">Figure 4
<p>Performance of MBABC-NM with respect to average distance.</p>
Full article ">
21 pages, 313 KiB  
Article
Nonuniform Dichotomy with Growth Rates of Skew-Evolution Cocycles in Banach Spaces
by Ariana Găină, Mihail Megan and Rovana Boruga (Toma)
Axioms 2023, 12(4), 394; https://doi.org/10.3390/axioms12040394 - 18 Apr 2023
Cited by 1 | Viewed by 1181
Abstract
This paper presents integral charaterizations for nonuniform dichotomy with growth rates and their correspondents for the particular cases of nonuniform exponential dichotomy and nonuniform polynomial dichotomy of skew-evolution cocycles in Banach spaces. The connections between these three concepts are presented. Full article
18 pages, 1382 KiB  
Article
A Quasiconformal-Based Geometric Model for Craniofacial Analysis and Its Application
by Ming-Hei Wong, Meixi Li, King-Man Tam, Hoi-Man Yuen, Chun-Ting Au, Kate Ching-Ching Chan, Albert Martin Li and Lok-Ming Lui
Axioms 2023, 12(4), 393; https://doi.org/10.3390/axioms12040393 - 18 Apr 2023
Viewed by 1967
Abstract
We address the problem of craniofacial morphometric analysis using geometric models, which has important clinical applications for the diagnosis of syndromes associated with craniofacial dysmorphologies. In this work, a novel geometric model is proposed to analyze craniofacial structures based on local curvature information [...] Read more.
We address the problem of craniofacial morphometric analysis using geometric models, which has important clinical applications for the diagnosis of syndromes associated with craniofacial dysmorphologies. In this work, a novel geometric model is proposed to analyze craniofacial structures based on local curvature information and Teichmüller mappings. A key feature of the proposed model is that its pipeline starts with few two-dimensional images of the human face captured at different angles, from which the three-dimensional craniofacial structure can be reconstructed. The 3D surface reconstruction from 2D images is based on a modified 3D morphable model (3DMM) framework. Geometric quantities around important feature landmarks according to different clinical applications can then be computed on each three-dimensional craniofacial structure. Together with the Teichmüller mapping, the landmark-based Teichmüller curvature distances (LTCDs) for every classes can be computed, which are further used for three-class classification. A composite score model is used and the parameter optimization is carried out to further improve the classification accuracy. Our proposed model is applied to study the craniofacial structures of children with and without the obstructive sleep apnoea (OSA). Sixty subjects, with accessible multi-angle photography and polysomnography (PSG) data, are divided into three classes based on the severity of OSA. Using our proposed model, our proposed model achieves a high 90% accuracy, which outperforms other existing models. This demonstrates the effectiveness of our proposed geometric model for craniofacial analysis. Full article
Show Figures

Figure 1

Figure 1
<p>Demonstration of the 3D facial model reconstructed from 2D photos.</p>
Full article ">Figure 2
<p>Demonstration the facial landmark points (red dots) superimposed on a sample 3D face model.</p>
Full article ">Figure 3
<p>Automatic landmark selection by the Euclidean Distance Method.</p>
Full article ">Figure 4
<p>Landmarks having high classification power chosen by the bagging algorithm.</p>
Full article ">Figure 5
<p>Reasonable expectation of the distribution of the template feature vector <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mn>0</mn> </msub> <mo>,</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>A new subject <span class="html-italic">K</span> that is equidistant from <math display="inline"><semantics> <msub> <mi>T</mi> <mn>0</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>T</mi> <mn>1</mn> </msub> </semantics></math> is needed to be classified.</p>
Full article ">Figure 7
<p>Take <math display="inline"><semantics> <msub> <mi>d</mi> <mrow> <mi>k</mi> <mn>2</mn> </mrow> </msub> </semantics></math> into the consideration of classification.</p>
Full article ">Figure 8
<p>Classification accuracy against number of feature used for the binary-class model.</p>
Full article ">Figure 9
<p>Classification accuracy against number of feature used for the three-class model.</p>
Full article ">
26 pages, 1010 KiB  
Article
Theoretical Validation of New Two-Dimensional One-Variable-Power Copulas
by Christophe Chesneau
Axioms 2023, 12(4), 392; https://doi.org/10.3390/axioms12040392 - 18 Apr 2023
Cited by 2 | Viewed by 1384
Abstract
One of the most effective ways to illustrate the relationship between two quantitative variables is to describe the corresponding two-dimensional copula. This approach is acknowledged as practical, nonredundant, and computationally manageable in the context of data analysis. Modern data, however, contain a wide [...] Read more.
One of the most effective ways to illustrate the relationship between two quantitative variables is to describe the corresponding two-dimensional copula. This approach is acknowledged as practical, nonredundant, and computationally manageable in the context of data analysis. Modern data, however, contain a wide variety of dependent structures, and the copulas now in use may not provide the best model for all of them. As a result, researchers seek to innovate by building novel copulas with appealing properties that are also based on original methodologies. The foundations are theoretical; for a copula to be validated, it must meet specific requirements, which frequently dictate the constraints that must be placed on the relevant parameters. In this article, we make a contribution to the understudied field of one-variable-power copulas. We first identify the specific assumptions that, in theory, validate copulas of such nature. Some other general copulas and inequalities are discussed. Our general results are illustrated with numerous examples depending on two or three parameters. We also prove that strong connections exist between our assumptions and well-established distributions. To highlight the importance of our findings, we emphasize a particular two-parameter, one-variable-power copula that unifies the definition of some other copulas. We reveal its versatile shapes, related functions, various symmetry, Archimedean nature, geometric invariance, copula ordering, quadrant dependence, tail dependence, correlations, and distribution generation. Numerical tables and graphics are produced to support some of these properties. The estimation of the parameters based on data is discussed. As a complementary contribution, two new, intriguing one-variable-power copulas beyond the considered general form are finally presented and studied. Full article
(This article belongs to the Special Issue Statistical Modeling of Modern Multivariate Data)
Show Figures

Figure 1

Figure 1
<p>Plots of the TE copula for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math> (<b>left</b>) and <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>20</mn> </mrow> </semantics></math> (<b>right</b>).</p>
Full article ">Figure 2
<p>Plots of the IP copula for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> (<b>left</b>) and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math> (<b>right</b>).</p>
Full article ">Figure 3
<p>Plots of the IP copula for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> (<b>left</b>) and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math> (<b>right</b>).</p>
Full article ">Figure 4
<p>Plots of the IP copula density for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> (<b>left</b>) and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math> (<b>right</b>).</p>
Full article ">Figure 5
<p>Plots of the IP copula density for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> (<b>left</b>) and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math> (<b>right</b>).</p>
Full article ">Figure 6
<p>Scatter plot of the data (<b>left</b>) and plot of the estimated IP copula density (<b>right</b>).</p>
Full article ">Figure 7
<p>Scatter plot of the data (<b>left</b>) and plot of the estimated IP copula density (<b>right</b>).</p>
Full article ">Figure 8
<p>Plots of the Int1 copula (<b>left</b>) and the Int1 copula density (<b>right</b>).</p>
Full article ">Figure 9
<p>Plots of the Int2 copula (<b>left</b>) and the Int2 copula density (<b>right</b>) for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mi>γ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 10
<p>Plots of the Int2 copula (<b>left</b>) and the Int2 copula density (<b>right</b>) for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>2</mn> <mo>/</mo> <mn>3</mn> </mrow> </semantics></math>.</p>
Full article ">
14 pages, 749 KiB  
Article
Dynamics of a Double-Impulsive Control Model of Integrated Pest Management Using Perturbation Methods and Floquet Theory
by Fahad Al Basir, Jahangir Chowdhury and Delfim F. M. Torres
Axioms 2023, 12(4), 391; https://doi.org/10.3390/axioms12040391 - 18 Apr 2023
Cited by 1 | Viewed by 1589
Abstract
We formulate an integrated pest management model to control natural pests of the crop through the periodic application of biopesticide and chemical pesticides. In a theoretical analysis of the system pest eradication, a periodic solution is found and established. All the system variables [...] Read more.
We formulate an integrated pest management model to control natural pests of the crop through the periodic application of biopesticide and chemical pesticides. In a theoretical analysis of the system pest eradication, a periodic solution is found and established. All the system variables are proved to be bounded. Our main goal is then to ensure that pesticides are optimized, in terms of pesticide concentration and pesticide application frequency, and that the optimum combination of pesticides is found to provide the most benefit to the crop. By using Floquet theory and the small amplitude perturbation method, we prove that the pest eradication periodic solution is locally and globally stable. The acquired results establish a threshold time limit for the impulsive release of various controls as well as some valid theoretical conclusions for effective pest management. Furthermore, after a numerical comparison, we conclude that integrated pest management is more effective than single biological or chemical controls. Finally, we illustrate the analytical results through numerical simulations. Full article
(This article belongs to the Special Issue Mathematical Methods in the Applied Sciences)
Show Figures

Figure 1

Figure 1
<p>Impact of biopesticide application in impulsive mode on system (<a href="#FD1-axioms-12-00391" class="html-disp-formula">1</a>). Evolution of (<b>a</b>) crop; (<b>b</b>) uninfected pest; (<b>c</b>) infected pest; (<b>d</b>) biopesticide. The set of parameters are <math display="inline"><semantics> <mrow> <mi>r</mi> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>0.2</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>0.003</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>0.8</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.6</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>c</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.8</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>0.15</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.2</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>=</mo> <mn>0.05</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>κ</mi> <mo>=</mo> <mn>100</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>s</mi> <mo>=</mo> <mn>0.3</mn> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>λ</mi> <mo>=</mo> <mn>0.35</mn> </mrow> </semantics></math>. Here, the time interval is <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math> days and the rates of biopesticide release are <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> (black line), <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math> (red line), and <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>12</mn> </mrow> </semantics></math> (blue line).</p>
Full article ">Figure 2
<p>Impact of biopesticide on system (<a href="#FD1-axioms-12-00391" class="html-disp-formula">1</a>) for different impulse intervals and rates. Evolution of (<b>a</b>) crop; (<b>b</b>) uninfected pest; (<b>c</b>) infected pest; (<b>d</b>) biopesticide. Red line indicates <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>, green line indicates <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>12</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>, and blue line indicates <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>12</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Impact of both biopesticide and chemical pesticide on system (<a href="#FD1-axioms-12-00391" class="html-disp-formula">1</a>) with the same impulse interval, <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>τ</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>. Evolution of (<b>a</b>) crop; (<b>b</b>) uninfected pest; (<b>c</b>) infected pest; (<b>d</b>) biopesticide; (<b>e</b>) chemical pesticide. The rates of impulses are <math display="inline"><semantics> <mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0.15</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math> for black dotted color; <math display="inline"><semantics> <mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math> for red dashed line; and <math display="inline"><semantics> <mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0.05</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math> for blue solid line.</p>
Full article ">Figure 4
<p>Impact of both biopesticide and chemical pesticide on system (<a href="#FD1-axioms-12-00391" class="html-disp-formula">1</a>) with same impulse interval, <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>3</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>τ</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, and where the rates of impulses are <math display="inline"><semantics> <mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0.15</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math>. Evolution of (<b>a</b>) crop; (<b>b</b>) uninfected pest; (<b>c</b>) infected pest; (<b>d</b>) pesticides (biopesticide in blue, chemical pesticide in red).</p>
Full article ">
14 pages, 1643 KiB  
Article
Stability Switching in Lotka-Volterra and Ricker-Type Predator-Prey Systems with Arbitrary Step Size
by Shamika Kekulthotuwage Don, Kevin Burrage, Kate J. Helmstedt and Pamela M. Burrage
Axioms 2023, 12(4), 390; https://doi.org/10.3390/axioms12040390 - 17 Apr 2023
Cited by 1 | Viewed by 1707
Abstract
Dynamical properties of numerically approximated discrete systems may become inconsistent with those of the corresponding continuous-time system. We present a qualitative analysis of the dynamical properties of two-species Lotka-Volterra and Ricker-type predator-prey systems under discrete and continuous settings. By creating an arbitrary time [...] Read more.
Dynamical properties of numerically approximated discrete systems may become inconsistent with those of the corresponding continuous-time system. We present a qualitative analysis of the dynamical properties of two-species Lotka-Volterra and Ricker-type predator-prey systems under discrete and continuous settings. By creating an arbitrary time discretisation, we obtain stability conditions that preserve the characteristics of continuous-time models and their numerically approximated systems. Here, we show that even small changes to some of the model parameters may alter the system dynamics unless an appropriate time discretisation is chosen to return similar dynamical behaviour to what is observed in the corresponding continuous-time system. We also found similar dynamical properties of the Ricker-type predator-prey systems under certain conditions. Our results demonstrate the need for preliminary analysis to identify which dynamical properties of approximated discretised systems agree or disagree with the corresponding continuous-time systems. Full article
(This article belongs to the Special Issue Mathematical Modeling with Differential Equations)
Show Figures

Figure 1

Figure 1
<p>Demonstration of the different dynamics that can arise with the discrete Ricker-type system (1) when the parameter <math display="inline"><semantics> <mi>α</mi> </semantics></math> is varied. This system is solved with <math display="inline"><semantics> <mrow> <mi>K</mi> <mo>=</mo> <mn>2500</mn> <mo>,</mo> <mi>γ</mi> <mo>=</mo> <mn>0.01</mn> <mo>,</mo> <mi>c</mi> <mo>=</mo> <mn>0.2</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>r</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>. For three different <math display="inline"><semantics> <mi>α</mi> </semantics></math> values <math display="inline"><semantics> <mrow> <mn>0.05</mn> <mo>,</mo> <mn>0.048</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mn>0.04</mn> </mrow> </semantics></math>, the predator-prey populations diverge, converge very slowly and converge, respectively. Here, system (1) is derived for a unit step size, which is similar to system (2) when <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>Different predator-prey dynamics of discrete Ricker-type and Lotka-Volterra models with slightly varying parameter <span class="html-italic">c</span> as <math display="inline"><semantics> <mrow> <mi>c</mi> <mo>+</mo> <mi>ζ</mi> </mrow> </semantics></math> by <math display="inline"><semantics> <mrow> <mi>ζ</mi> <mo>=</mo> <mo>{</mo> <mn>0</mn> <mo>,</mo> <mn>0.01</mn> <mo>,</mo> <mn>0.05</mn> <mo>,</mo> <mn>0.1</mn> <mo>,</mo> <mn>0.2</mn> <mo>}</mo> </mrow> </semantics></math> values where <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mi>K</mi> <mo>=</mo> <mn>2500</mn> <mo>,</mo> <mi>γ</mi> <mo>=</mo> <mn>0.01</mn> <mo>,</mo> <mi>α</mi> <mo>=</mo> <mn>0.05</mn> <mo>,</mo> <mi>c</mi> <mo>=</mo> <mn>0.2</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>r</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Stability regions of the discrete Ricker-type model and discrete Lotka-Volterra model as a function of <span class="html-italic">h</span> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mi>α</mi> <mi>γ</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>K</mi> <mo>=</mo> <mn>2500</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>r</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>. The fixed-point convergence region is bounded by <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mfrac> <mi>c</mi> <mi>K</mi> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mi>h</mi> </mfrac> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mfrac> <mi>c</mi> <mi>K</mi> </mfrac> </mrow> </semantics></math>, and boundary changes are marked in red and black lines for different <math display="inline"><semantics> <mrow> <mi>c</mi> <mo>=</mo> <mo>{</mo> <mn>0.1</mn> <mo>,</mo> <mn>0.2</mn> <mo>,</mo> <mn>0.3</mn> <mo>}</mo> </mrow> </semantics></math>. The stability regions are coloured for <math display="inline"><semantics> <mrow> <mi>c</mi> <mo>=</mo> <mn>0.2</mn> </mrow> </semantics></math>, as represented in solid lines, and are represented as dashed lines for <math display="inline"><semantics> <mrow> <mi>c</mi> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>c</mi> <mo>=</mo> <mn>0.3</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>The fixed-point convergence region as a function of <math display="inline"><semantics> <mi>β</mi> </semantics></math> and <span class="html-italic">c</span> for the discrete Ricker-type model and discrete Lotka-Volterra model where <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mi>α</mi> <mi>γ</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>K</mi> <mo>=</mo> <mn>2500</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>r</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>. The fixed-point convergence region is bounded by <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mfrac> <mi>c</mi> <mi>K</mi> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mi>h</mi> </mfrac> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mfrac> <mi>c</mi> <mi>K</mi> </mfrac> <mo>,</mo> <mo>∀</mo> <mi>h</mi> <mo>&gt;</mo> <mn>0</mn> </mrow> </semantics></math>. For <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>, the stability regions are coloured, and the upper and lower boundaries of the fixed-point region are plotted for <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mfrac> <mi>c</mi> <mi>K</mi> </mfrac> <mo>+</mo> <mn>1</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mfrac> <mi>c</mi> <mi>K</mi> </mfrac> </mrow> </semantics></math>, as displayed in solid red and black lines, respectively. The upper boundary of fixed-point region moves upward with decreasing step size, marked as red dashed lines. Note that the lower boundary of the fixed-point convergence region is valid for any <span class="html-italic">h</span>.</p>
Full article ">Figure 5
<p>Special behaviour of predator-prey populations for Ricker-type and Lotka-Volterra discrete models if <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mi>θ</mi> </mfrac> <mo>=</mo> <mfrac> <mn>20</mn> <mn>21</mn> </mfrac> </mrow> </semantics></math>, where <math display="inline"><semantics> <mrow> <mi>K</mi> <mo>=</mo> <mn>2500</mn> <mo>,</mo> <mi>γ</mi> <mo>=</mo> <mn>0.01</mn> <mo>,</mo> <mi>α</mi> <mo>=</mo> <mn>0.05</mn> <mo>,</mo> <mi>c</mi> <mo>=</mo> <mn>0.2</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>r</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>. Predator-prey populations seem to converge to a fixed point at the beginning; however, after a long time, the populations oscillate around the fixed point. Note that this exceptional case occurs only at the upper bound of the fixed-point convergence region.</p>
Full article ">
25 pages, 539 KiB  
Article
Efficiency of Orthogonal Matching Pursuit for Group Sparse Recovery
by Chunfang Shao, Xiujie Wei, Peixin Ye and Shuo Xing
Axioms 2023, 12(4), 389; https://doi.org/10.3390/axioms12040389 - 17 Apr 2023
Cited by 2 | Viewed by 1511
Abstract
We propose the Group Orthogonal Matching Pursuit (GOMP) algorithm to recover group sparse signals from noisy measurements. Under the group restricted isometry property (GRIP), we prove the instance optimality of the GOMP algorithm for any decomposable approximation norm. Meanwhile, we show the robustness [...] Read more.
We propose the Group Orthogonal Matching Pursuit (GOMP) algorithm to recover group sparse signals from noisy measurements. Under the group restricted isometry property (GRIP), we prove the instance optimality of the GOMP algorithm for any decomposable approximation norm. Meanwhile, we show the robustness of the GOMP under the measurement error. Compared with the P-norm minimization approach, the GOMP is easier to implement, and the assumption of γ-decomposability is not required. The simulation results show that the GOMP is very efficient for group sparse signal recovery and significantly outperforms Basis Pursuit in both scalability and solution quality. Full article
(This article belongs to the Special Issue Applied Mathematical Modeling and Optimization)
Show Figures

Figure 1

Figure 1
<p>Examples of group 50-sparse signals in different dimensions: (<b>a</b>) An example of the group 50-sparse signal in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>128</mn> </mrow> </semantics></math>. (<b>b</b>) An example of the group 50-sparse signal in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>256</mn> </mrow> </semantics></math>. (<b>c</b>) An example of the group 50-sparse signal in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>512</mn> </mrow> </semantics></math>. (<b>d</b>) An example of the group 50-sparse signal in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>1024</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>The recovery of an input signal via GOMP in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>512</mn> </mrow> </semantics></math> with group sparsity level <math display="inline"><semantics> <mrow> <mi>K</mi> <mo>=</mo> <mn>50</mn> </mrow> </semantics></math> and number of measurements <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>200</mn> </mrow> </semantics></math> under different noises: (<b>a</b>) The recovery of an input signal via GOMP under the noise <math display="inline"><semantics> <mrow> <mi>e</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>. (<b>b</b>) The recovery of an input signal via GOMP under a Gaussian noise <span class="html-italic">e</span> from <math display="inline"><semantics> <mrow> <mi mathvariant="script">N</mi> <mo>(</mo> <mn>0</mn> <mo>,</mo> <msup> <mrow> <mn>0.1</mn> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>The recovery of an input signal via BP in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>512</mn> </mrow> </semantics></math> with group sparsity level <math display="inline"><semantics> <mrow> <mi>K</mi> <mo>=</mo> <mn>50</mn> </mrow> </semantics></math> and number of measurements <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>200</mn> </mrow> </semantics></math> under different noises: (<b>a</b>) The recovery of an input signal via BP under the noise <math display="inline"><semantics> <mrow> <mi>e</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>. (<b>b</b>) The recovery of an input signal via BP under a Gaussian noise <span class="html-italic">e</span> from <math display="inline"><semantics> <mrow> <mi mathvariant="script">N</mi> <mo>(</mo> <mn>0</mn> <mo>,</mo> <msup> <mrow> <mn>0.1</mn> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>The percentage of the support of 100 input signals correctly recovered as a function of number <span class="html-italic">N</span> of Gaussian measurements for different group sparsity levels <span class="html-italic">K</span> in dimension <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>=</mo> <mn>1024</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>MSE and running time of the GOMP in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>1024</mn> </mrow> </semantics></math> with number of measurements <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>500</mn> </mrow> </semantics></math> under the noise <math display="inline"><semantics> <mrow> <mi>e</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>MSE and running time of the GOMP in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>1024</mn> </mrow> </semantics></math> with number of measurements <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>500</mn> </mrow> </semantics></math> under a Gaussian noise <span class="html-italic">e</span> from <math display="inline"><semantics> <mrow> <mi mathvariant="script">N</mi> <mo>(</mo> <mn>0</mn> <mo>,</mo> <msup> <mrow> <mn>0.1</mn> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>MSE and running time of BP in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>1024</mn> </mrow> </semantics></math> with number of measurements <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>500</mn> </mrow> </semantics></math> under the noise <math display="inline"><semantics> <mrow> <mi>e</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 8
<p>MSE and running time of BP in dimension <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>1024</mn> </mrow> </semantics></math> with number of measurements <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>500</mn> </mrow> </semantics></math> under a Gaussian noise <span class="html-italic">e</span> from <math display="inline"><semantics> <mrow> <mi mathvariant="script">N</mi> <mo>(</mo> <mn>0</mn> <mo>,</mo> <msup> <mrow> <mn>0.1</mn> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">
12 pages, 431 KiB  
Article
A Numerical Approach of Handling Fractional Stochastic Differential Equations
by Iqbal M. Batiha, Ahmad A. Abubaker, Iqbal H. Jebril, Suha B. Al-Shaikh and Khaled Matarneh
Axioms 2023, 12(4), 388; https://doi.org/10.3390/axioms12040388 - 17 Apr 2023
Cited by 9 | Viewed by 2201
Abstract
This work proposes a new numerical approach for dealing with fractional stochastic differential equations. In particular, a novel three-point fractional formula for approximating the Riemann–Liouville integrator is established, and then it is applied to generate approximate solutions for fractional stochastic differential equations. Such [...] Read more.
This work proposes a new numerical approach for dealing with fractional stochastic differential equations. In particular, a novel three-point fractional formula for approximating the Riemann–Liouville integrator is established, and then it is applied to generate approximate solutions for fractional stochastic differential equations. Such a formula is derived with the use of the generalized Taylor theorem coupled with a recent definition of the definite fractional integral. Our approach is compared with the approximate solution generated by the Euler–Maruyama method and the exact solution for the purpose of verifying our findings. Full article
(This article belongs to the Special Issue Stochastic Processes: Theory and Application)
Show Figures

Figure 1

Figure 1
<p>The Wiener process.</p>
Full article ">Figure 2
<p>The numerical solution (<a href="#FD29-axioms-12-00388" class="html-disp-formula">29</a>) of the FSDE (<a href="#FD28-axioms-12-00388" class="html-disp-formula">28</a>) according to different fractional-order values.</p>
Full article ">Figure 3
<p>Comparison between the numerical solution (<a href="#FD29-axioms-12-00388" class="html-disp-formula">29</a>), Euler–Maruyama’s solution, and exact solution of problem (<a href="#FD28-axioms-12-00388" class="html-disp-formula">28</a>).</p>
Full article ">Figure 4
<p>Absolute error generated by the numerical solution (<a href="#FD29-axioms-12-00388" class="html-disp-formula">29</a>) and Euler–Maruyama’s solution.</p>
Full article ">Figure 5
<p>The numerical solution (<a href="#FD31-axioms-12-00388" class="html-disp-formula">31</a>) of the FSDE (<a href="#FD30-axioms-12-00388" class="html-disp-formula">30</a>) according to different fractional-order values.</p>
Full article ">Figure 6
<p>Comparison between the numerical solution (<a href="#FD31-axioms-12-00388" class="html-disp-formula">31</a>), Euler–Maruyama’s solution and exact solution of problem (<a href="#FD30-axioms-12-00388" class="html-disp-formula">30</a>).</p>
Full article ">Figure 7
<p>Absolute error generated by the numerical solution (<a href="#FD31-axioms-12-00388" class="html-disp-formula">31</a>) and Euler–Maruyama’s solution.</p>
Full article ">Figure 8
<p>The numerical solution (<a href="#FD33-axioms-12-00388" class="html-disp-formula">33</a>) of the FSDE (<a href="#FD32-axioms-12-00388" class="html-disp-formula">32</a>) according to different fractional-order values.</p>
Full article ">Figure 9
<p>Comparison between the numerical solution (<a href="#FD33-axioms-12-00388" class="html-disp-formula">33</a>), Euler–Maruyama’s solution and exact solution of problem (<a href="#FD32-axioms-12-00388" class="html-disp-formula">32</a>).</p>
Full article ">Figure 10
<p>Absolute error generated by the numerical solution (<a href="#FD33-axioms-12-00388" class="html-disp-formula">33</a>) and Euler–Maruyama’s solution.</p>
Full article ">
18 pages, 2678 KiB  
Article
Application of the Optimal Homotopy Asymptotic Approach for Solving Two-Point Fuzzy Ordinary Differential Equations of Fractional Order Arising in Physics
by Ali Fareed Jameel, Dulfikar Jawad Hashim, Nidal Anakira, Osama Ababneh, Ahmad Qazza, Abedel-Karrem Alomari and Khamis S. Al Kalbani
Axioms 2023, 12(4), 387; https://doi.org/10.3390/axioms12040387 - 17 Apr 2023
Cited by 2 | Viewed by 1491
Abstract
This work focuses on solving and analyzing two-point fuzzy boundary value problems in the form of fractional ordinary differential equations (FFOBVPs) using a new version of the approximation analytical approach. FFOBVPs are useful in describing complex scientific phenomena that include heritable characteristics and [...] Read more.
This work focuses on solving and analyzing two-point fuzzy boundary value problems in the form of fractional ordinary differential equations (FFOBVPs) using a new version of the approximation analytical approach. FFOBVPs are useful in describing complex scientific phenomena that include heritable characteristics and uncertainty, and obtaining exact or close analytical solutions for these equations can be challenging, especially in the case of nonlinear problems. To address these difficulties, the optimal homotopy asymptotic method (OHAM) was studied and extended in a new form to solve FFOBVPs. The OHAM is known for its ability to solve both linear and nonlinear fractional models and provides a straightforward methodology that uses multiple convergence control parameters to optimally manage the convergence of approximate series solutions. The new form of the OHAM presented in this work incorporates the concepts of fuzzy sets theory and some fractional calculus principles to include fuzzy analysis in the method. The steps of fuzzification and defuzzification are used to transform the fuzzy problem into a crisp problem that can be solved using the OHAM. The method is demonstrated by solving and analyzing linear and nonlinear FFOBVPs at different values of fractional derivatives. The results obtained using the new form of the fuzzy OHAM are analyzed and compared to those found in the literature to demonstrate the method’s efficiency and high accuracy in the fuzzy domain. Overall, this work presents a feasible and efficient approach for solving FFOBVPs using a new form of the OHAM with fuzzy analysis. Full article
(This article belongs to the Special Issue Mathematical Methods in the Applied Sciences)
Show Figures

Figure 1

Figure 1
<p>The three-dimensional approximate solution of Equation (31) given by the third-order FF-OHAM over all <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>∈</mo> <mfenced open="[" close="]" separators="|"> <mrow> <mn>0,0.5</mn> </mrow> </mfenced> <mo>,</mo> </mrow> </semantics></math> and for all <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>∈</mo> <mfenced open="[" close="]" separators="|"> <mrow> <mn>0,1</mn> </mrow> </mfenced> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>The three-dimensional approximate solution for Equation (31) given using the fifth-order FF-OHAM over all <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>∈</mo> <mfenced open="[" close="]" separators="|"> <mrow> <mn>0,0.5</mn> </mrow> </mfenced> <mo>,</mo> </mrow> </semantics></math> and for all <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>∈</mo> <mo>[</mo> <mn>0,1</mn> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Comparison of the lower approximate solution for Equation (31) using the fifth-order FF-OHAM and the fifth-order SCM for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∀</mo> <mi>x</mi> <mo>∈</mo> <mo>[</mo> <mn>0,1</mn> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Comparison of the upper approximate solution for Equation (31) using the fifth-order FF-OHAM and the fifth-order SCM for <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∀</mo> <mi>x</mi> <mo>∈</mo> <mo>[</mo> <mn>0,1</mn> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>The three-dimensional approximate solution for Equation (47) given using the tenth-order FF-OHAM over all <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>∈</mo> <mfenced open="[" close="]" separators="|"> <mrow> <mn>0,0.1</mn> </mrow> </mfenced> </mrow> </semantics></math> at <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">ω</mi> <mo>=</mo> <mn>1.9</mn> </mrow> </semantics></math> and for <math display="inline"><semantics> <mrow> <mi>η</mi> <mo>=</mo> <mn>0.6</mn> <mo>,</mo> </mrow> </semantics></math> and for all <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>∈</mo> <mo>[</mo> <mn>0,1</mn> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">
24 pages, 15344 KiB  
Article
Developing a Deep Learning-Based Defect Detection System for Ski Goggles Lenses
by Dinh-Thuan Dang and Jing-Wein Wang
Axioms 2023, 12(4), 386; https://doi.org/10.3390/axioms12040386 - 17 Apr 2023
Cited by 7 | Viewed by 2178
Abstract
Ski goggles help protect the eyes and enhance eyesight. The most important part of ski goggles is their lenses. The quality of the lenses has leaped with technological advances, but there are still defects on their surface during manufacturing. This study develops a [...] Read more.
Ski goggles help protect the eyes and enhance eyesight. The most important part of ski goggles is their lenses. The quality of the lenses has leaped with technological advances, but there are still defects on their surface during manufacturing. This study develops a deep learning-based defect detection system for ski goggles lenses. The first step is to design the image acquisition model that combines cameras and light sources. This step aims to capture clear and high-resolution images on the entire surface of the lenses. Next, defect categories are identified, including scratches, watermarks, spotlight, stains, dust-line, and dust-spot. They are labeled to create the ski goggles lenses defect dataset. Finally, the defects are automatically detected by fine-tuning the mobile-friendly object detection model. The mentioned defect detection model is the MobileNetV3 backbone used in a feature pyramid network (FPN) along with the Faster-RCNN detector. The fine-tuning includes: replacing the default ResNet50 backbone with a combination of MobileNetV3 and FPN; adjusting the hyper-parameter of the region proposal network (RPN) to suit the tiny defects; and reducing the number of the output channel in FPN to increase computational performance. Our experiments demonstrate the effectiveness of defect detection; additionally, the inference speed is fast. The defect detection accuracy achieves a mean average precision (mAP) of 55%. The work automatically integrates all steps, from capturing images to defect detection. Furthermore, the lens defect dataset is publicly available to the research community on GitHub. The repository address can be found in the Data Availability Statement section. Full article
(This article belongs to the Special Issue Various Deep Learning Algorithms in Computational Intelligence)
Show Figures

Figure 1

Figure 1
<p>Some lens samples of ski goggles in different sizes, curvatures, and colors.</p>
Full article ">Figure 2
<p>Some standard sizes of ski goggles lenses. The width ranges from 92 to 205 mm, and the height ranges from 60 to 126 mm.</p>
Full article ">Figure 3
<p>Some main components in the optical inspection system.</p>
Full article ">Figure 4
<p>The flowchart of the proposed method. There are two modules for processing images. The first one is the image acquisition to capture the raw image, extract regions of interest, and data labeling. The second module is defect detection, which involves training data and inferencing defects. This module combines Faster R-CNN, MobileNetV3, and FPN to create the customized end-to-end model. It is compatible with data of the ski goggles lenses.</p>
Full article ">Figure 5
<p>The ski goggles’ lens sample. It is wide and curved; thus, we mark its surface to be easily controlled by cameras.</p>
Full article ">Figure 6
<p>The most that a custom light source meets the curvature of the ski goggles lens. Five dot matrix LED modules are connected by an angle of 125°.</p>
Full article ">Figure 7
<p>Design diagram of the image acquisition system. Five cameras are placed at the top. The custom light source (yellow) is placed at the bottom. The ski goggles’ lens sample (red) is positioned in the middle and held on two sides.</p>
Full article ">Figure 8
<p>The actual model of image acquisition. The bottom is a custom light source that shines through the ski goggles lens surface to the cameras. The computer controls the five cameras through the acquisition card. The developed program will access the card’s interface to capture images simultaneously.</p>
Full article ">Figure 9
<p>Input and Output of image acquisition. Input is some lens samples of ski goggles, and output is raw images from five cameras.</p>
Full article ">Figure 10
<p>Crop regions of interest from raw images. Five cameras capture five images in the first row. To facilitate data labeling, extracting parts of good images is necessary. The <b>bottom</b> row is five regions of interest.</p>
Full article ">Figure 11
<p>Some defect samples on the surface of ski goggles lenses. Column (<b>a</b>): scratch defects, (<b>b</b>): watermark, (<b>c</b>): spotlight, (<b>d</b>): stain, (<b>e</b>): dust-line, (<b>f</b>): dust-spot.</p>
Full article ">Figure 12
<p>The GUI of LabelMe: the image annotation tool used to label defects on the surface of ski goggles lenses.</p>
Full article ">Figure 13
<p>Synthetic image generated from the flip technique.</p>
Full article ">Figure 14
<p>This is an instant of the inverted-residual block in MobileNetV3 architecture. First is a convolutional expand layer that widens channels from 24 to 72. Second is the convolutional depth-wise layer for better efficiency than traditional convolution. Its input and output channels are equal to 72, and the striking attribute of convolution halves the resolution. Next is the squeeze and excitation module to improve the power of features in the network. The final convolutional projection layer presents features in the lower dimension space, from 72 to 40.</p>
Full article ">Figure 15
<p>Backbone: the feature pyramid network and MobileNetV3 backbone together. The input is the image of size <math display="inline"><semantics> <mrow> <mi>H</mi> <mo>,</mo> <mi>W</mi> </mrow> </semantics></math>. Firstly, The MobileNetV3 extracts the image to many multi-scale feature maps <math display="inline"><semantics> <mrow> <mo>{</mo> <msub> <mrow> <mi>C</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mrow> <mi>C</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>,</mo> <msub> <mrow> <mi>C</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msub> <mo>,</mo> <msub> <mrow> <mi>C</mi> </mrow> <mrow> <mn>4</mn> </mrow> </msub> <mo>,</mo> <msub> <mrow> <mi>C</mi> </mrow> <mrow> <mn>5</mn> </mrow> </msub> <mo>}</mo> </mrow> </semantics></math>. Secondly, the multi-scale output of MobileNetV3 is the input for FPN. The final result is the feature maps at multiple levels <math display="inline"><semantics> <mrow> <mo>{</mo> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>,</mo> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msub> <mo>,</mo> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mn>4</mn> </mrow> </msub> <mo>,</mo> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mn>5</mn> </mrow> </msub> <mo>}</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 16
<p>Illustrate how to create the box anchors in the region proposal network. Left is the image containing some red ground-truth boxes. RPN generates the reference boxes called “anchors” to map to ground-truth boxes. The multi-scale anchors are generated on the right image at various positions. They are the rectangular boxes marked with white, black, yellow, green, and blue colors.</p>
Full article ">Figure 17
<p>Box regression transforms the proposal anchor set to ground-truth label set.</p>
Full article ">Figure 18
<p>Illustration of transformation <math display="inline"><semantics> <mrow> <mi>δ</mi> </mrow> </semantics></math> from the anchor <math display="inline"><semantics> <mrow> <mi>a</mi> </mrow> </semantics></math> to the ground-truth box <math display="inline"><semantics> <mrow> <mi>b</mi> </mrow> </semantics></math>. The formula is in Equation (12).</p>
Full article ">Figure 19
<p>Selected examples of defect detection results. Defects are marked by the red rectangular boxes and the label above is the defect category. (<b>a</b>–<b>c</b>) display a variety of defect types, including dust-spot and dust-line. (<b>d</b>) showcases the spotlight defect type, while (<b>e</b>) highlights the stain defect type. (<b>f</b>–<b>h</b>) and h feature the scratch defect type.</p>
Full article ">
17 pages, 3610 KiB  
Article
A New Perspective for Multivariate Time Series Decision Making through a Nested Computational Approach Using Type-2 Fuzzy Integration
by Martha Ramirez and Patricia Melin
Axioms 2023, 12(4), 385; https://doi.org/10.3390/axioms12040385 - 17 Apr 2023
Cited by 3 | Viewed by 1540
Abstract
The integration of key indicators from the results of the analysis of time series represents a constant challenge within organizations; this could be mainly due to the need to establish the belonging of each indicator within a process, geographic region or category. This [...] Read more.
The integration of key indicators from the results of the analysis of time series represents a constant challenge within organizations; this could be mainly due to the need to establish the belonging of each indicator within a process, geographic region or category. This paper thus illustrates how both primary and secondary indicators are relevant for decision making, and why they need to be integrated by making new final fuzzy indicators. Thus, our proposal consists of a type-2 fuzzy integration of multivariate time series, such as OECD country risk classification, inflation, population and gross national income (GNI) by using multiple type-1 fuzzy inference systems to perform time series classification tasks. Our contribution consists of the proposal to integrate multiple nested type-1 fuzzy inference systems using a type-2 fuzzy integration. Simulation results show the advantages of using the proposed method for the fuzzy classification of multiple time series. This is done in order so the user can have tools that allow them to understand the environment and generate comparative analyses of multiple information sources, and finally use it during the process prior to decision making considering the main advantage of modeling the inherent uncertainty. Full article
Show Figures

Figure 1

Figure 1
<p>OECD member countries: total annual population.</p>
Full article ">Figure 2
<p>OECD members: average annual GNI.</p>
Full article ">Figure 3
<p>OECD members: inflation annual percentage.</p>
Full article ">Figure 4
<p>Illustration of the proposed method.</p>
Full article ">Figure 5
<p>First Type-1 fuzzy system population and GNI variables.</p>
Full article ">Figure 6
<p>Second Type-1 fuzzy system GNI variables.</p>
Full article ">Figure 7
<p>Third Type-1 fuzzy system inflation and OECD country risk variables.</p>
Full article ">Figure 8
<p>Type-2 fuzzy criteria integration.</p>
Full article ">
22 pages, 1406 KiB  
Article
Barrier Options and Greeks: Modeling with Neural Networks
by Nneka Umeorah, Phillip Mashele, Onyecherelam Agbaeze and Jules Clement Mba
Axioms 2023, 12(4), 384; https://doi.org/10.3390/axioms12040384 - 17 Apr 2023
Cited by 3 | Viewed by 2867
Abstract
This paper proposes a non-parametric technique of option valuation and hedging. Here, we replicate the extended Black–Scholes pricing model for the exotic barrier options and their corresponding Greeks using the fully connected feed-forward neural network. Our methodology involves some benchmarking experiments, which result [...] Read more.
This paper proposes a non-parametric technique of option valuation and hedging. Here, we replicate the extended Black–Scholes pricing model for the exotic barrier options and their corresponding Greeks using the fully connected feed-forward neural network. Our methodology involves some benchmarking experiments, which result in an optimal neural network hyperparameter that effectively prices the barrier options and facilitates their option Greeks extraction. We compare the results from the optimal NN model to those produced by other machine learning models, such as the random forest and the polynomial regression; the output highlights the accuracy and the efficiency of our proposed methodology in this option pricing problem. The results equally show that the artificial neural network can effectively and accurately learn the extended Black–Scholes model from a given simulated dataset, and this concept can similarly be applied in the valuation of complex financial derivatives without analytical solutions. Full article
(This article belongs to the Special Issue Various Deep Learning Algorithms in Computational Intelligence)
Show Figures

Figure 1

Figure 1
<p>Visualization plot.</p>
Full article ">Figure 2
<p>Train/test MAE values for Models A1, A2 and A3; (<b>a</b>) MAE—Model A1; (<b>b</b>) MAE—Model A2; (<b>c</b>) MAE—Model A3.</p>
Full article ">Figure 3
<p>Train/Test MAE values for Models B1, B2 and B3; (<b>a</b>) MAE—Model B1; (<b>b</b>) MAE—Model B2; (<b>c</b>) MAE—Model B3.</p>
Full article ">Figure 4
<p>Train/Test LOSS values for Models A1, A2 and A3; (<b>a</b>) LOSS—Model A1; (<b>b</b>) LOSS—Model A2; (<b>c</b>) MAE—Model A3.</p>
Full article ">Figure 5
<p>Train/Test LOSS values for Models B1, B2 and B3; (<b>a</b>) LOSS—Model B1; (<b>b</b>) LOSS—Model B2; (<b>c</b>) LOSS—Model B3.</p>
Full article ">Figure 6
<p>Option values visualization for Models A1, A2, A3; (<b>a</b>) Model A1: Regression plot; (<b>b</b>) Model A1: Histogram plot; (<b>c</b>) Model A1: Box plot; (<b>d</b>) Model A2: Regression plot; (<b>e</b>) Model A2: Histogram plot; (<b>f</b>) Model A2: Box plot; (<b>g</b>) Model A3: Regression plot; (<b>h</b>) Model A3: Histogram plot; (<b>i</b>) Model A3: Box plot.</p>
Full article ">Figure 7
<p>Option values visualization for Models B1, B2, B3; (<b>a</b>) Model B1: Regression plot; (<b>b</b>) Model B1: Histogram plot; (<b>c</b>) Model B1: Box plot; (<b>d</b>) Model B2: Regression plot; (<b>e</b>) Model B2: Histogram plot; (<b>f</b>) Model B2: Box plot; (<b>g</b>) Model B3: Regression plot; (<b>h</b>) Model B3: Histogram plot; (<b>i</b>) Model B3: Box plot.</p>
Full article ">Figure 8
<p>Ranking of models for optimality.</p>
Full article ">Figure 9
<p>Option values and Greeks; (<b>a</b>) DO option value; (<b>b</b>) DO delta; (<b>c</b>) DO gamma; (<b>d</b>) DO vega.</p>
Full article ">
27 pages, 613 KiB  
Article
A Compound Class of Inverse-Power Muth and Power Series Distributions
by Leonardo Barrios-Blanco, Diego I. Gallardo, Héctor J. Gómez and Marcelo Bourguignon
Axioms 2023, 12(4), 383; https://doi.org/10.3390/axioms12040383 - 16 Apr 2023
Cited by 1 | Viewed by 1631
Abstract
This paper introduces the inverse-power Muth power series model, which is a composition of the inverse-power Muth and the class of power series distributions. The use of the Bell distribution in this context is emphasized for the first time in the literature. Probability [...] Read more.
This paper introduces the inverse-power Muth power series model, which is a composition of the inverse-power Muth and the class of power series distributions. The use of the Bell distribution in this context is emphasized for the first time in the literature. Probability density, survival and hazard functions are studied, as well as their moments. Using the stochastic representation of the model, the maximum-likelihood estimators are implemented by the use of the expectation-maximization algorithm, while standard errors are calculated using Oakes’ method. Monte Carlo simulation studies are conducted to show the performance of the maximum-likelihood estimators in finite samples. Two applications to real datasets are shown, where our proposal is compared with some models based on power series compositions. Full article
(This article belongs to the Special Issue Probability, Statistics and Estimation)
Show Figures

Figure 1

Figure 1
<p>PDF and hazard function for the IPM-P distribution with <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>PDF and hazard function for the IPM-B distribution with <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>EM algorithm scheme for the IPM-PS distribution.</p>
Full article ">Figure 4
<p>PDF of the IPM-B, WP, GP, and W models on the IMC dataset (<b>left panel</b>) and the empirical CDF compared with the estimated CDF for the same models (<b>right panel</b>).</p>
Full article ">Figure 5
<p>QQ plot for the QR for the IPM-B, BBXII, WP, GP, IPM, and W distributions for the beta-carotene data.</p>
Full article ">Figure 5 Cont.
<p>QQ plot for the QR for the IPM-B, BBXII, WP, GP, IPM, and W distributions for the beta-carotene data.</p>
Full article ">Figure 6
<p>PDF of IPM-P, WG, GP, GG, IPM, and W models for failure time dataset (<b>left panel</b>) and the empirical CDF versus the estimated CDF for the same models (<b>right panel</b>).</p>
Full article ">Figure 7
<p>QQ plot of the QR for IPM-P, BBXII, WG, GG, IPM, and W distributions for the failure time dataset.</p>
Full article ">Figure 7 Cont.
<p>QQ plot of the QR for IPM-P, BBXII, WG, GG, IPM, and W distributions for the failure time dataset.</p>
Full article ">Figure A1
<p>PDF and hazard function for the IPM-G distribution, with <math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure A2
<p>PDF and hazard function for the IPM-L distribution, with <math display="inline"><semantics> <mrow> <mi>θ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">
30 pages, 457 KiB  
Article
Infinite Series and Logarithmic Integrals Associated to Differentiation with Respect to Parameters of the Whittaker Wκ,μ(x) Function II
by Alexander Apelblat and Juan Luis González-Santander
Axioms 2023, 12(4), 382; https://doi.org/10.3390/axioms12040382 - 16 Apr 2023
Viewed by 1308
Abstract
In the first part of this investigation, we considered the parameter differentiation of the Whittaker function Mκ,μx. In this second part, first derivatives with respect to the parameters of the Whittaker function Wκ,μx are [...] Read more.
In the first part of this investigation, we considered the parameter differentiation of the Whittaker function Mκ,μx. In this second part, first derivatives with respect to the parameters of the Whittaker function Wκ,μx are calculated. Using the confluent hypergeometric function, these derivatives can be expressed as infinite sums of quotients of the digamma and gamma functions. Furthermore, it is possible to obtain these parameter derivatives in terms of infinite integrals, with integrands containing elementary functions (products of algebraic, exponential, and logarithmic functions), from the integral representation of Wκ,μx. These infinite sums and integrals can be expressed in closed form for particular values of the parameters. Finally, an integral representation of the integral Whittaker function wiκ,μx and its derivative with respect to κ, as well as some reduction formulas for the integral Whittaker functions Wiκ,μx and wiκ,μx, are calculated. Full article
29 pages, 433 KiB  
Article
Infinite Series and Logarithmic Integrals Associated to Differentiation with Respect to Parameters of the Whittaker Mκ,μ(x) Function I
by Alexander Apelblat and Juan Luis González-Santander
Axioms 2023, 12(4), 381; https://doi.org/10.3390/axioms12040381 - 16 Apr 2023
Viewed by 1398
Abstract
In this paper, first derivatives of the Whittaker function Mκ,μx are calculated with respect to the parameters. Using the confluent hypergeometric function, these derivarives can be expressed as infinite sums of quotients of the digamma and gamma functions. Moreover, [...] Read more.
In this paper, first derivatives of the Whittaker function Mκ,μx are calculated with respect to the parameters. Using the confluent hypergeometric function, these derivarives can be expressed as infinite sums of quotients of the digamma and gamma functions. Moreover, from the integral representation of Mκ,μx it is possible to obtain these parameter derivatives in terms of finite and infinite integrals with integrands containing elementary functions (products of algebraic, exponential, and logarithmic functions). These infinite sums and integrals can be expressed in closed form for particular values of the parameters. For this purpose, we have obtained the parameter derivative of the incomplete gamma function in closed form. As an application, reduction formulas for parameter derivatives of the confluent hypergeometric function are derived, along with finite and infinite integrals containing products of algebraic, exponential, logarithmic, and Bessel functions. Finally, reduction formulas for the Whittaker functions Mκ,μx and integral Whittaker functions Miκ,μx and miκ,μx are calculated. Full article
16 pages, 326 KiB  
Article
Pascu-Rønning Type Meromorphic Functions Based on Sălăgean-Erdély–Kober Operator
by Sheza M. El-Deeb, Gangadharan Murugusundaramoorthy, Kaliappan Vijaya and Alhanouf Alburaikan
Axioms 2023, 12(4), 380; https://doi.org/10.3390/axioms12040380 - 16 Apr 2023
Cited by 1 | Viewed by 1264
Abstract
In the present investigation, we introduce a new class of meromorphic functions defined in the punctured unit disk Δ*:={ϑC:0<|ϑ|<1} by making use of the Erdély–Kober operator [...] Read more.
In the present investigation, we introduce a new class of meromorphic functions defined in the punctured unit disk Δ*:={ϑC:0<|ϑ|<1} by making use of the Erdély–Kober operator Iς,ϱτ,κ which unifies well-known classes of the meromorphic uniformly convex function with positive coefficients. Coefficient inequalities, growth and distortion inequalities, in addition to closure properties are acquired. We also set up a few outcomes concerning convolution and the partial sums of meromorphic functions in this new class. We additionally state some new subclasses and its characteristic houses through specializing the parameters that are new and no longer studied in association with the Erdély–Kober operator thus far. Full article
(This article belongs to the Special Issue Mathematical Analysis and Applications IV)
29 pages, 1913 KiB  
Article
A New Extension of the Kumaraswamy Exponential Model with Modeling of Food Chain Data
by Eman A. Eldessouky, Osama H. Mahmoud Hassan, Mohammed Elgarhy, Eid A. A. Hassan, Ibrahim Elbatal and Ehab M. Almetwally
Axioms 2023, 12(4), 379; https://doi.org/10.3390/axioms12040379 - 16 Apr 2023
Cited by 6 | Viewed by 2605
Abstract
Statistical models are useful in explaining and forecasting real-world occurrences. Various extended distributions have been widely employed for modeling data in a variety of fields throughout the last few decades. In this article we introduce a new extension of the Kumaraswamy exponential (KE) [...] Read more.
Statistical models are useful in explaining and forecasting real-world occurrences. Various extended distributions have been widely employed for modeling data in a variety of fields throughout the last few decades. In this article we introduce a new extension of the Kumaraswamy exponential (KE) model called the Kavya–Manoharan KE (KMKE) distribution. Some statistical and computational features of the KMKE distribution including the quantile (QUA) function, moments (MOms), incomplete MOms (INMOms), conditional MOms (COMOms) and MOm generating functions are computed. Classical maximum likelihood and Bayesian estimation approaches are employed to estimate the parameters of the KMKE model. The simulation experiment examines the accuracy of the model parameters by employing Bayesian and maximum likelihood estimation methods. We utilize two real datasets related to food chain data in this work to demonstrate the importance and flexibility of the proposed model. The new KMKE proposed distribution is very flexible, more so than numerous well-known distributions. Full article
(This article belongs to the Special Issue Mathematical Modelling in Sustainable Global Supply Chain Management)
Show Figures

Figure 1

Figure 1
<p>A detailed graphic representation of the article.</p>
Full article ">Figure 2
<p>Plots of pdf and hrf for the KMKE model.</p>
Full article ">Figure 3
<p>Plots of the pdf for the KMKE model in 3D.</p>
Full article ">Figure 3 Cont.
<p>Plots of the pdf for the KMKE model in 3D.</p>
Full article ">Figure 4
<p>Plots of the hrf for the KMKE model in 3D.</p>
Full article ">Figure 5
<p>Plots of mean, var, SK, KU, CV and ID in 3D for the KMKE model.</p>
Full article ">Figure 6
<p>Heatmaps of MSE values for parameters of the KMWE distribution with different sample cases: <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> <mi>β</mi> <mo>=</mo> <mn>0.4</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>Heatmaps of MSE values for parameters of the KMWE distribution with different sample cases: <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> <mi>β</mi> <mo>=</mo> <mn>1.5</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 8
<p>Heatmaps of MSE values for parameters of KMWE distribution with different sample cases: <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> <mi>β</mi> <mo>=</mo> <mn>1.5</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>Profile MLE of the KMWE model for food chain data.</p>
Full article ">Figure 10
<p>MLE of cdf, and pdf with empirical and histogram, QQ and PP of the KMWE model for food chain data.</p>
Full article ">Figure 11
<p>MCMC plot and convergence line for parameters of the KMWE model for food chain data.</p>
Full article ">Figure 12
<p>Histogram plot with normal curve for parameters of the KMWE model for food chain data.</p>
Full article ">Figure 13
<p>Profile MLE of the KMWE model for food and drink wholesaling data.</p>
Full article ">Figure 14
<p>MLE of cdf, and pdf with empirical and histogram, QQ and PP of the KMWE model for food and drink wholesaling data.</p>
Full article ">Figure 15
<p>MCMC plot and convergence line for parameters of the KMWE model for food and drink wholesaling data.</p>
Full article ">Figure 16
<p>Histogram plot with normal curve for parameters of the KMWE model for food and drink wholesaling data.</p>
Full article ">
13 pages, 424 KiB  
Article
Reinsurance Policy under Interest Force and Bankruptcy Prohibition
by Yangmin Zhong and Huaping Huang
Axioms 2023, 12(4), 378; https://doi.org/10.3390/axioms12040378 - 16 Apr 2023
Viewed by 1566
Abstract
In this paper, we solve an optimal reinsurance problem in the mathematical finance area. We assume that the surplus process of the insurance company follows a controlled diffusion process and the constant interest rate is involved in the financial model. During the whole [...] Read more.
In this paper, we solve an optimal reinsurance problem in the mathematical finance area. We assume that the surplus process of the insurance company follows a controlled diffusion process and the constant interest rate is involved in the financial model. During the whole optimization period, the company has a choice to buy reinsurance contract and decide the reinsurance retention level. Meanwhile, the bankruptcy at the terminal time is not allowed. The aim of the optimization problem is to minimize the distance between the terminal wealth and a given goal by controlling the reinsurance proportion. Using the stochastic control theory, we derive the Hamilton-Jacobi-Bellman equation for the optimization problem. Via adopting the technique of changing variable as well as the dual transformation, an explicit solution of the value function and the optimal policy are shown. Finally, several numerical examples are shown, from which we find several main factors that affect the optimal reinsurance policy. Full article
Show Figures

Figure 1

Figure 1
<p>The optimal value function <span class="html-italic">s</span> with respect to <span class="html-italic">z</span> at time <math display="inline"><semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>The optimal reinsurance policy with respect to <span class="html-italic">z</span> at time <math display="inline"><semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>The optimal reinsurance policy with respect to <span class="html-italic">z</span> at time <math display="inline"><semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>The optimal reinsurance policy with respect to <span class="html-italic">z</span> under different interest rates <math display="inline"><semantics> <mrow> <mi>i</mi> <mo>=</mo> <mn>0.05</mn> <mo>,</mo> <mn>0.1</mn> <mo>,</mo> <mn>0.15</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>The optimal reinsurance policy with respect to <span class="html-italic">z</span> under different volatility rates <math display="inline"><semantics> <mrow> <mi>n</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mn>1.5</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>The optimal reinsurance policy with respect to <span class="html-italic">z</span> with different reinsurance safety loading <math display="inline"><semantics> <mrow> <mi>ϱ</mi> <mo>=</mo> <mn>0.4</mn> <mo>,</mo> <mn>0.5</mn> <mo>,</mo> <mn>0.6</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>The optimal reinsurance policy under different expected losses in unit time <math display="inline"><semantics> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>1.5</mn> <mo>,</mo> <mn>2</mn> </mrow> </semantics></math>.</p>
Full article ">
12 pages, 323 KiB  
Article
A Note on Finite Coarse Shape Groups
by Ivan Jelić and Nikola Koceić-Bilan
Axioms 2023, 12(4), 377; https://doi.org/10.3390/axioms12040377 - 14 Apr 2023
Cited by 1 | Viewed by 1118
Abstract
In this paper, we investigate properties concerning some recently introduced finite coarse shape invariants—the k-th finite coarse shape group of a pointed topological space and the k-th relative finite coarse shape group of a pointed topological pair. We define the notion [...] Read more.
In this paper, we investigate properties concerning some recently introduced finite coarse shape invariants—the k-th finite coarse shape group of a pointed topological space and the k-th relative finite coarse shape group of a pointed topological pair. We define the notion of finite coarse shape group sequence of a pointed topological pair X,X0,x0 as an analogue of homotopy and (coarse) shape group sequences and show that for any pointed topological pair, the corresponding finite coarse shape group sequence is a chain. On the other hand, we construct an example of a pointed pair of metric continua whose finite coarse shape group sequence fails to be exact. Finally, using the aforementioned pair of metric continua together with a pointed dyadic solenoid, we show that finite coarse-shape groups, in general, differ from both shape and coarse-shape groups. Full article
(This article belongs to the Special Issue Advances in General Topology and Its Application)
Previous Issue
Next Issue
Back to TopTop