Abstract:In the present work, an attempted was made to develop a numerical algorithm by the use of new orthogonal hybrid functions formed from hybrid of piecewise constant orthogonal sample-and-hold functions and piecewise linear orthogonal triangular functions to obtain the numerical solution of multi-order fractional differential equations. The construction of numerical algorithm involves a formula, consisting of generalized one-shot operational matrices, expressing explicitly the Riemann-Liouville fractional order integral of the new orthogonal hybrid functions in terms of new orthogonal hybrid functions themselves. A set of test problems comprising of linear and nonlinear multi-order fractional differential equations with constant and variable coefficients was considered to demonstrate the accuracy and computational efficiency of our algorithm and to compare our results with those acquired by some other well-known methods used for solving multi-order fractional differential equations in the literature. The comparison highlighted that our algorithm exhibits superior performance to those methods.

What problem does this paper attempt to address?

### Problems the paper attempts to solve The paper aims to solve the numerical solution problem of multi - order fractional differential equations (FDEs). Specifically, the author has developed a new numerical algorithm, using a new type of hybrid orthogonal functions (HFs) composed of piecewise - constant orthogonal sample - and - hold functions (SHFs) and piecewise - linear orthogonal triangular functions (TFs) to obtain the numerical solutions of multi - order fractional differential equations. ### Background and motivation In fractional calculus, solving differential equations containing non - integer - order derivatives or integrals has always been a research hotspot. Traditional analytical methods such as Fourier transform and Laplace transform can only be used to solve linear fractional differential equations, and it is difficult to find analytical solutions for nonlinear or complex linear equations. Therefore, it is necessary to develop efficient and accurate numerical methods. ### Proposed method The author proposes the following steps to solve the problem: 1. **Construct new hybrid orthogonal functions**: By combining SHFs and TFs, new hybrid orthogonal functions HFs are constructed. 2. **Derive the generalized single - operation matrix**: Using the generalized single - operation matrix, the form of Riemann - Liouville fractional - order integral under the new orthogonal hybrid functions is clearly expressed. 3. **Convert to algebraic equations**: Convert multi - order fractional differential equations into algebraic equations, thereby simplifying the solution process. 4. **Testing and verification**: Through a series of test problems of linear and nonlinear multi - order fractional differential equations, the accuracy and computational efficiency of the proposed algorithm are verified. ### Results and conclusions By comparing with existing methods, the paper shows that its proposed numerical algorithm based on new hybrid orthogonal functions has higher precision and computational efficiency in solving multi - order fractional differential equations. Specifically: - It shows superior performance in different types of multi - order fractional differential equations. - It is applicable to linear and nonlinear equations, as well as equations with constant and variable coefficients. - Compared with other commonly used methods (such as Adomian decomposition method, Podlubny numerical method, etc.), it can provide more accurate results. In conclusion, the paper proposes a novel and effective numerical algorithm, providing new ideas and tools for solving multi - order fractional differential equations.

Novel hybrid function operational matrices of fractional integration: An application for solving multi-order fractional differential equations

Piecewise linear approximate solution of fractional order non-stiff and stiff differential-algebraic equations by orthogonal hybrid functions

New orthogonal hybrid function based numerical method to solve system of fractional order differential equations

Novel derivative operational matrix in Caputo sense with applications

A New Operational Matrices-Based Spectral Method For Multi-Order Fractional Problems

Numerical Technique Based on Operational Matrices of Fractional Integration Using ψ ‐Shifted Chebyshev Polynomials

A numerical technique for solving multi-dimensional fractional optimal control problems using fractional wavelet method

Novel operational matrices-based method for solving fractional-order delay differential equations via shifted Gegenbauer polynomials

Numerical solution of nonlinear partial quadratic integro‐differential equations of fractional order via hybrid of block‐pulse and parabolic functions

Analytical and Numerical Methods for Solving Second-Order Two-Dimensional Symmetric Sequential Fractional Integro-Differential Equations

A Robust Scheme Based on Novel‐operational Matrices for Some Classes of Time‐fractional Nonlinear Problems Arising in Mechanics and Mathematical Physics

An approximate method for solving fractional optimal control problems by hybrid functions

Numerical study of generalized modified Caputo fractional differential equations

Treatment of fractional multi-order/multi-term differential equations: utilizing fractional shifted Lucas polynomials

A Novel Method for Linear Systems of Fractional Ordinary Differential Equations with Applications to Time-Fractional PDEs

Comparison of fractional-order generalized wavelets and orthonormal wavelets methods for solving multi-dimensional fractional optimal control problems

The numerical solution of fractional integral equations via orthogonal polynomials in fractional powers

Enhanced shifted Jacobi operational matrices of integrals: spectral algorithm for solving some types of ordinary and fractional differential equations

Linearized Novel Operational Matrices-Based Scheme for Classes of Nonlinear Time-Space Fractional Unsteady Problems in 2D

Solution of System of Linear Fractional Differential Equations with Modified derivative of Jumarie Type

Comparison of Two Hybrid Functions for Numerical Solution of Nonlinear Mixed Partial Integro-Differential Equations