In this paper, the optimal approximation algorithm is proposed to simplify non-linear functions and/or discrete data as piecewise polynomials by using the constrained least squares. In time-sensitive applications or in embedded systems with limited resources, the runtime of the approximate function is as crucial as its accuracy. The proposed algorithm search to find the Optimal Piecewise Polynomial (OPP) with minimum computational cost while ensuring the error below a specified threshold. This was accomplished by using smooth piecewise polynomials with optimal order and number of intervals. The computational cost only depends on polynomial complexity, i.e., the order and the number of intervals at runtime function call. For optimal approximation, computational costs for all the possible combinations of piecewise polynomials were calculated and tabulated as ascending order for the specific target CPU off-line. Each combination was optimized through constrained least squares and random selection method for given sample points. Afterward, whether the approximation error is below the predetermined value is examined. When the error is permissible, the combination is selected as the optimal approximation or the next combination was examined. To verify the performance, several representative functions were examined and analyzed.

A complete MAPLE procedure is designed to implement effectively an algorithm for approximating the trigonometric functions. The algorithm gives a piecewise polynomial approximation on an arbitrary interval, presenting a special partition that we can get its parts, subintervals with ending points of finite rational numbers, together with corresponding approximate polynomials. The procedure takes a sequence of pairs of interval-polynomial as its output that we can easily explore in some useful ways. Examples on calculating approximate values of the sine function with arbitrary accuracy for both of rational and irrational arguments as well as drawing the graph of the piecewise approximate functions will be presented. Moreover, from the approximate integration of integrands of the form $x^m\sin x$ on $[a,b]$, another MAPLE procedure is proposed to find the desired polynomial estimates in norm for the best $L^2$-approximation of the sine function in the vector space $\mathcal{P}_{\ell}$ of polynomials of degree at most $\ell$, a subspace of $L^2(a,b)$.

In this study, we introduce the dynamics of a Hepatitis B virus (HBV) model with the class of asymptomatic carriers and conduct a comprehensive analysis to explore its theoretical aspects and examine the crossover effect within the HBV model. To investigate the crossover behavior of the operators, we divide the study interval into two subintervals. In the first interval, the classical derivative is employed to study the qualitative properties of the proposed system, while in the second interval, we utilize the ABC fractional differential operator. Consequently, the study is initiated using the piecewise Atangana-Baleanu derivative framework for the systems. The HBV model is then analyzed to determine the existence, Hyers-Ulam (HU) stability, and disease-free equilibrium point of the model. Moreover, we showcase the application of an Adams-type predictor-corrector (PC) technique for Atangana-Baleanu derivatives and an extended Adams-Bashforth-Moulton (ABM) method for Caputo derivatives through numerical results. Subsequently, we employ computational methods to numerically solve the models and visually present the obtained outcomes using different fractional-order values. This network is designed to provide more precise information for disease modeling, considering that communities often interact with one another, and the rate of disease spread is influenced by this factor.

The continuing public health issue known as COVID-19 (the 2019 Novel Corona virus infection) has a global emphasis. Despite (or perhaps because of) the fact that there are significant gaps in our understanding of COVID-19 epidemiology, transmission dynamics, research methods, and management breakout poses a new kind of global hazard. The good news is that there is currently enough knowledge about the epidemic process to allow for the creation of mathematical forecasting models. We modify a conventional SEIR epidemic model to the unique dynamic compartments and epidemic features of COVID 19 as it spreads in a population with a diverse age structure. Although many US states and other nations around the world followed lockdown and reopening processes, we perform some analysis on using some techniques of the epidemic course. A new perspective of fractional calculus known as piecewise derivatives of fractional order is used to study the proposed model. Sufficient conditions are established to show the existence theory. In addition, a numerical scheme based on Newton’s polynomials is established to simulate the approximate solutions of the proposed model by using various fractional orders. Some real data results are also shown with comparison of the numerical results.

There are many uncertainties and risks in residential electricity consumption during the economic development. Knowledge of the relationship between residential electricity consumption and its key determinant—income—are important to the sustainable development of electric power industry. Using panel data from 30 provinces for the 1995-2012 period, this study investigates how residential electricity consumption changes as incomes increase in China. Previous studies typically used linear or quadratic double-logarithmic models imposing ex ante restrictions on the indistinct relationship between residential electricity consumption and income. Contrary to those models, we employed a reduced piecewise linear model that is self-adaptive and highly flexible and circumvents the problem of “prior restrictions.” Robust tests of different segment specifications and regression methods are performed to ensure the conservatism of the research. The results provide strong evidence that the income elasticity was approximately one, and it remained stable throughout the estimation period. The income threshold at which residential electricity consumption automatically remains stable or slows has not been reached. To ensure the sustainable development of the electric power industry, introducing higher energy efficiency standards for electrical appliances and improving income levels are vital. And government should emphasize electricity conservation in industrial sector rather than in residential sector.

In the present study, piecewise integrated composite (PIC) bumper beam for passenger cars was proposed and design optimisation process for composite bumper beam against IIHS test was carried out with the help of machine learning. Several elements in IIHS bumper FE model have been assigned to be references, in order to collect training data which, allow the machine learning model to study the method of predicting loading types of each finite element. 2-D and 3-D implementations were provided by machine learning models, which determined stacking sequences of each finite element in PIC bumper beam. It was found that the PIC bumper beam, which was designed by machine learning model has direct impact on reducing the possibility of failure as well as increasing bending strength effectively than conventional composite bumper beam. Moreover, 3-D implementation produced better results compared with 2-D implementation since it was preferable to choose loading type information which was achieved from surroundings when the target elements were located either at corner or junction of planes instead of using information came from the same plane of target.

Only cross-sectional studies have demonstrated muscle deficits in children with spastic cerebral palsy (SCP). The impact of functional impairments on altered muscle growth remains unclear. This prospective longitudinal study modelled the morphological muscle growth in 87 children with SCP (age range 6 months to 11 years, Gross Motor Function Classification System [GMFCS] level I/II/III=47/22/18). Repeated ultrasound assessments were performed during 2-year (y) follow-up with an interval of minimal 6 months. Three-dimensional freehand ultrasound was applied to assess medial gastrocnemius muscle volume (MV), mid-belly cross-sectional area (CSA) and muscle belly length (ML). Non-linear mixed models compared trajectories of (normalized) muscle growth between GMFCS-I versus GMFCS-II&III. MV and CSA growth trajectories showed piecewise model with 2 breakpoints, with the highest growth before 2y and negative growth rates after 6-9y. Before 2y, children with GMFCS-II&III already showed lower growth rates compared to GMFCS-I. From 2-9y, the growth rates did not differ between GMFCS-levels. After 9y, a more pronounced reduction in normalized CSA was observed in GMFCS-II&III. Different trajectories in ML growth were shown between the GMFCS-level subgroups. These longitudinal trajectories could be used to monitor the SCP muscle pathology during childhood and to optimize treatment planning and goals aiming to stimulate muscle growth

Estimating the vehicle crashworthiness parameters experimentally is expensive and time consuming. For these reasons different modelling approaches are utilized to predict the vehicle behaviour and reduce the need for full-scale crash testing. The earlier numerical methods used for vehicle crashworthiness analysis were based on the use of lumped parameters models (LPM), a combination of masses and nonlinear springs interconnected in various configurations. Nowadays, the explicit nonlinear finite element analysis (FEA) is probably the most widely recognized modelling technique. Although informative, finite element models (FEM) of vehicle crash are expensive both in terms of man-hours put into assembling the model and related computational costs. A simpler analytical tool for early analysis of vehicle crashworthiness could greatly assist the modelling and save time. In this paper a simple piecewise LPM composed of a mass-spring-damper system, is used to estimate the vehicle crashworthiness parameters, focusing on the dynamic crush and the acceleration severity index (ASI). The model is first calibrated against a full-scale crash test and a FEM, post-processed with the LS-DYNA software, at an impact velocity of 56 km/h. The genetic algorithm is used to calibrate the model by estimating the piecewise lumped parameters (stiffness and damping of the front structure of the vehicle). After calibration, the LPM is applied to a range of velocities (40, 48, 64 and 72 km/h). The predictions for crashworthiness parameters from the LPM were compared with the predictions from the FEA and the results are much similar. It is shown that the LPM can assist in crash analysis, since LPM has some predictive capabilities and requires less computation time in comparison with the explicit nonlinear FEA.

In this work, an innovative numerical approach for polylinear approximation (polylinearization) of non-self-intersecting compact sensor characteristics (transfer functions) specified either pointwise or analytically is introduced. The goal is to optimally partition the sensor characteristic, i.e. to select the vertices of the approximating polyline (approximant) along with their positions on, the sensor characteristics so that the distance (i.e. the separation) between the approximant and the characteristic is rendered below a certain problem-specific tolerance. To achieve this goal, two alternative non-linear optimization problems are solved, whose essential difference is in the adopted quantitative measure of the separation between the transfer function and the approximant. In the first problem, which relates to absolutely integrable sensor characteristics (their energy is not necessarily finite, but they can be represented in terms of convergent Fourier series), the polylinearization is constructed by numerical minimization of the L^1--metric (a distance-based separation measure), concerning the number of polyline vertices and their locations. In the second problem which covers the quadratically integrable sensor characteristics (whose energy is finite, but they do not necessarily admit a representation in terms of convergent Fourier series), the polylinearization is constructed by minimizing numerically the L^2-metric (area-, or energy-based separation measure) for the same set of optimization variables –the locations, and the number of polyline vertices.

This paper focus on the modeling and time-domain analyzing method of piecewise-smooth fractional-order circuit systems. Firstly, an analytical framework is constructed based on a generalized piecewise-smooth two-port impedance network model, which takes the topological characteristics of Dickson type switched-capacitor converters (SCCs) as examples, and considers the fractional-order characteristics of class-2 X7R ceramic capacitors, which are widely used in SCCs. Secondly, the time-domain analyzing method for such a piecewise-smooth fractional-order circuit system is developed, numerical solutions are obtained, and the influence of the fractional-order capacitor on the power conversion efficiency of the system are revealed. Finally, in order to validate the analysis, both numerical simulation and experiments are provided in the work, the results of which confirm that, the proposed fractional-order analytical framework can capture the power loss mechanism of piecewise-smooth fractional-order circuit systems more effectively than the traditional integer-order based approach, thus benefiting the design of circuit systems.

The concepts of fractional calculus and the related techniques have been gaining momentum in circuit system fields, and an increasing body of evidence suggests that fractional-order characteristics are widely distributed in electronic components. In this paper, we consider the fractional characteristics of non-solid aluminum electrolytic capacitors and focus on the nonlinear dynamic characteristics of peak-current-mode controlled (PCMC) dc-dc boost converters with such components. In the work, we establish a piecewise smooth non-commensurate fractional-order model to describe the converter, in which the fractional-order equivalent impedance model of capacitors is adopted and the detailed control loop is considered. In order to analyze the converter, we use the fractional Adams-Bashforth-Moulton typed method (F-ABM) and stroboscopic map technique to obtain the time-domain solutions and the bifurcation diagram. The results of bifurcation diagrams indicate that, with the change of different parameters, the converter enters chaotic state through period doubling bifurcation. In order to provide empirical evidence for the nonlinear dynamic analysis, we consider both circuit-level simulation and experiments in the work, the results of which further confirm the results of theoretical analysis.

Since the conventional split-merge algorithm is sensitive to the object scale variance and splitting starting point, a piecewise split-merge polygon approximation method is proposed to extract the object contour features. Specifically, the contour corner is used as the starting point for the contour piecewise approximation to reduce the sensitivity of the contour segment on the starting point; then, the split-merge algorithm is used to implement the polygon approximation for each contour segments. Both the distance ratio and the arc length ratio instead of the distance error are used as the iterative stop condition to improve the robustness to the object scale variance. Both the angle and length as two features describe the shape of the contour polygon, and affect each other along the contour order relationship. Since they have a strong coupling relationship. To improve the description correction of the contour, these two features are combined to construct a Coupled Hidden Markov Model to detect the object by calculating the probability of the contour feature. The proposed algorithm is validated on ETHZ Shape Classes and INRIA Horses standard datasets. Compared with other contour-based object detection algorithms, the proposed algorithm reduces the complexity of contour description, improves the robustness of contour features to scale variance, and has a higher object detection rate.

This paper proposes training improvement methods of artificial neural networks (ANN) trajectory predictors. First, a dynamic power system time-series trajectory is split into several different segments to simplify the original ANN training problem. Moreover, the time-derivative of the trajectory is included to obtain an augmented loss function. Compared to previous studies which mainly focused on increasing the prediction accuracy, the aim of these novel techniques is to reduce the computational burden where the ANN output performance is still acceptable. The effectiveness of the developed methods is validated based on the WSCC three-machine nine-bus and IEEE 39-bus system models. The mean absolute error (MAE) and trajectory prediction results are analyzed, in which the numbers of neurons, hidden layers, and training epochs are constrained during the ANN training process. Rotor-angle difference between generators and the system frequency are investigated as the dynamic trajectories of the power system models. It is revealed that when the ANN architecture and epochs are constrained, the MAE results can be reduced, and the ANN training results can better reflect the original trajectory using the improvement approaches.

Refinery planning optimization is a challenging problem as regards handling the nonconvex bilinearity mainly due to pooling operations in processes such as crude oil distillation and product blending. This work investigates the performance of several representative piecewise-linear (or piecewise-affine) relaxation schemes (referred to as McCormick, bm, nf5, nf6t, and de (which is a new approach proposed based on eigenvector decomposition) that mainly give rise to mixed-integer optimization programs to convexify a bilinear term using predetermined univariate partitioning for instances of uniform and non-uniform partition sizes. Computational results show that applying these schemes give improved relaxation tightness than only applying convex and concave envelopes as estimators. Uniform partition sizes typically perform better in terms of relaxation solution quality and convergence behavior. It is also seen that there is a limit on the number of partitions that contributes to relaxation tightness, which does not necessarily correspond to a larger number of partitions, while a direct relation between relaxation size and tightness does not always hold for non-uniform partition sizes.

This study reports the development of a numerical method to simulate two-phase flows of Newtonian fluids that are incompressible, immiscible, and isothermal. The interface in the simulation was located and reconstructed using the geometric volume of fluid (VOF) method. The implementation of the Piecewise-Linear Interface Calculation (PLIC) scheme of the VOF method was performed to solve the three-dimensional (3D) interface transport during the dynamics of two-phase flows. In this method, the interface is approximated by a line segment in each interfacial cell. The balance of forces at the interface was accounted for using the continuum interfacial force (CSF) model. To solve the Navier-Stokes equations, meshless finite difference schemes from the HiG-Flow computational fluid dynamics software were employed. The 3D PLIC-VOF HiG-Flow algorithm was used to simulate several benchmark two-phase flows for the purpose of validating the numerical implementation. First, the performance of the PLIC implementation was evaluated by conducting two standard advection numerical tests: the 3D shearing flow test and the 3D deforming field test. Good agreement is obtained for the 3D interface shape using both the 3D PLIC-VOF HiG-Flow algorithm and those found on the scientific literature. In addition, the absolute error of the volume tracking advection calculation obtained by our 3D PLIC-VOF HiG-Flow algorithm was found to be smaller than the one found in the scientific literature for both the 3D shearing and 3D deforming flow tests. Lastly, the 3D bubble rising problem was simulated for different fluid densities and viscosity’s ratios. Again, good agreement is obtained for the 3D interface shape using both the newly implemented algorithm and the OpenFOAM, DROPS and NaSt3D software. The bubble’s rise velocity and evolution of the bubble’s center of mass is also computed with the 3D PLIC-VOF HiG-Flow algorithm and found to be in agreement with these software.

In this paper we address the approximation of the coupling problem for the wave equation and Maxwell's equations of electromagnetism in the time domain in terms of the electric field, by means of a nodal linear finite-element discretization in space, combined with a classical'explicit finite-difference scheme for the time discretization. Our study applies to the particular case where the dielectric permittivity has a constant value outside a sub-domain, whose closure does not intersect the boundary of the domain where the problem is defined. Inside this sub-domain Maxwell's equations hold. Outside this sub-domain the wave equation holds, which may correspond to the Maxwell's equations with a constant permittivity under certain conditions. We consider as a model the case of first order absorbing boundary conditions. Optimal error estimates that hold in natural norms under reasonable assumptions are given, among which lies a typical CFL condition for hyperbolic equations.

In real-life scenario, there are many mathematical tools to handle the incomplete and imprecise data. One of them is the fuzzy approach. This article aims to contribute to the literature of fuzzy multi-objective dynamic programming issues involving the fuzzy objective functions. The piecewise quadratic fuzzy numbers characterize these fuzzy parameters. Some basic notions in the problem under the α-pareto optimal solution concept is redefined and analyzed to study the stability of the problem. Furthermore, a technique, named as enhancing decomposition approach, is presented for achieving a subset for the parametric space that contains the sameα-pareto optimal solution. For a better understanding and comprehension of the suggested concept, a numerical example is provided.